Factors affecting fluid flow in strike–slip fault systems: coupled deformation and fluid flow modelling with application to the western Mount Isa Inlier, Australia

Deformation and focused fluid flow within a mineralized system are critical in the genesis of hydrothermal ore deposits. Dilation and integrated fluid flux due to coupled deformation and fluid flow in simple strike–slip fault geometries were examined using finite difference analysis in three dimensions. A series of generic fault bend and fault jog geometries consistent with those seen in the western Mount Isa Inlier were modelled in order to understand how fault geometry parameters influence the dilation and integrated fluid flux. Fault dip, fault width, bend/jog angle, and length were varied, and a cross-cutting fault and contrasting rock types were included. The results demonstrate that low fault dips, the presence of contrasts in rock type, and wide faults produce highest dilation and integrated fluid flux values. Increasing fault bend lengths and angles increases dilation and integrated fluid flux, but increasing fault jog length or angle has the opposite effect. There is minimal difference between the outputs from the releasing and restraining fault bend and jog geometries. Model characteristics producing greater fluid flows and/or gradients can be used in a predictive capacity in order to focus exploration on regions with more favorable fault geometries, provided that the mineralized rocks had Mohr–Coulomb rheologies similar to the ones used in the models.     

Timescales for reaching steady-state fluid flow in systems perturbed by formation of highly permeable faults

We consider the case of an isothermal, fluid‐saturated, homogeneous rock layer with transverse fluid flow driven by an imposed constant fluid pressure gradient. A rupture in the centre of the rock layer generates a highly permeable fault and results in a change of the initially homogeneous permeability distribution. This leads to a perturbation of the fluid flow field and its gradual transition to a new steady‐state corresponding to the new permeability distribution. An examination of this transitional process permits us to obtain an analytical estimation of the transition stage duration. The application of the results obtained to km‐scale faults in crystalline rock bodies leads to the conclusion that the evolution of the fluid velocity field is rather rapid compared with geological timescales.     

Structural controls on hydrothermal flow in a segmented rift system, Taupo Volcanic Zone, New Zealand

Vigorous hydrothermal convection transfers 10 times the average continental heat flow through the central Taupo Volcanic Zone (TVZ), a region of active extension (approximately 8 mm year^?1) and productive rhyolitic volcanism. Over 20 high‐temperature (>250°C) geothermal fields occur within Quaternary pyroclastic basins, with convective circulation to depths of 7–8 km presumably extending through basement rocks. Parallel‐striking normal faults, fractures and dikes dissect the convective regime, interacting with fluids to either enhance or restrict flow according to the relative permeability of structure and host rock. In the basement, high bulk permeability is maintained by focussed flow through faults and associated fractures well oriented for reactivation in the prevailing stress field. In contrast, distributed flow through fault‐bounded compartments prevails within Quaternary basins, masking any signal of deeper structural control. Exceptions occur where more competent rocks are exposed at the surface. As in narrow magmatic rifts elsewhere, the extensional fabric is partitioned into discrete rift segments linked along strike by accommodation zones. Eighty per cent of TVZ geothermal fields correlate spatially with rift architecture, with 60% located in accommodation zones. We suggest that segmented rift fabrics generate bulk permeability anisotropy that is to some extent predictable, with rift segments characterized by enhanced axial flow, and accommodation zones characterized by locally enhanced vertical permeability that is tectonically maintained. This provides a plausible explanation for the common occurrence of geothermal fields within accommodation zones and their notable absence within densely faulted rift segments. Maintenance of structural permeability in zones of active hydrothermal precipitation necessarily requires repeated brittle failure. Geothermal plumes therefore exploit tectonically maintained permeability within accommodation zones, with rift segments functioning mostly as drawdown regions. The influence of rift architecture on flow paths has important implications for geothermal extraction and epithermal mineral exploration within the TVZ and other structurally segmented hydrothermal systems, both active and extinct.     

Contrasting paleofluid systems in the continental basement: a fluid inclusion and stable isotope study of hydrothermal vein mineralization,Schwarzwald district,Germany

An integrated fluid inclusion and stable isotope study was carried out on hydrothermal veins (Sb‐bearing quartz veins, metal‐bearing fluorite–barite–quartz veins) from the Schwarzwald district, Germany. A total number of 106 Variscan (quartz veins related to Variscan orogenic processes) and post‐Variscan deposits were studied by microthermometry, Raman spectroscopy, and stable isotope analysis. The fluid inclusions in Variscan quartz veins are of the H₂O–NaCl–(KCl) type, have low salinities (0–10 wt.% eqv. NaCl) and high T_h values (150–350°C). Oxygen isotope data for quartz range from +2.8‰ to +12.2‰ and calculated δ¹⁸O_H2O values of the fluid are between ?12.5‰ and +4.4‰. The δD values of water extracted from fluid inclusions vary between ?49‰ and +4‰. The geological framework, fluid inclusion and stable isotope characteristics of the Variscan veins suggest an origin from regional metamorphic devolatilization processes. By contrast, the fluid inclusions in post‐Variscan fluorite, calcite, barite, quartz, and sphalerite belong to the H₂O–NaCl–CaCl₂ type, have high salinities (22–25 wt.% eqv. NaCl) and lower T_h values of 90–200°C. A low‐salinity fluid (0–15 wt.% eqv. NaCl) was observed in late‐stage fluorite, calcite, and quartz, which was trapped at similar temperatures. The δ¹⁸O values of quartz range between +11.1‰ and +20.9‰, which translates into calculated δ¹⁸O_H2O values between ?11.0‰ and +4.4‰. This range is consistent with δ¹⁸O_H2O values of fluid inclusion water extracted from fluorite (?11.6‰ to +1.1‰). The δD values of directly measured fluid inclusion water range between ?29‰ and ?1‰, ?26‰ and ?15‰, and ?63‰ and +9‰ for fluorite, quartz, and calcite, respectively. Calculations using the fluid inclusion and isotope data point to formation of the fluorite–barite–quartz veins under near‐hydrostatic conditions. The δ¹⁸O_H2O and δD data, particularly the observed wide range in δD, indicate that the mineralization formed through large‐scale mixing of a basement‐derived saline NaCl–CaCl₂ brine with meteoric water. Our comprehensive study provides evidence for two fundamentally different fluid systems in the crystalline basement. The Variscan fluid regime is dominated by fluids generated through metamorphic devolatilization and fluid expulsion driven by compressional nappe tectonics. The onset of post‐Variscan extensional tectonics resulted in replacement of the orogenic fluid regime by fluids which have distinct compositional characteristics and are related to a change in the principal fluid sources and the general fluid flow patterns. This younger system shows remarkably persistent geochemical and isotopic features over a prolonged period of more than 100 Ma.     

Transient fluid flow in and around a fault

We present the results of simple numerical experiments in which we study the evolution with time of fluid flow around and within a permeable fault embedded in a less permeable porous medium. Fluid movement is driven by an imposed vertical pressure gradient. The results show that fluid flow is controlled by two timescales: τ_f = Sl²/κ_F and τ_F = Sl²/κ_M, where S is the specific storage of the porous material, l the length of the fault, and κ_M and κ_F are the hydraulic conductivities of the porous material and the fault, respectively. Fluid flow and the associated fluid pressure field evolve through three temporal stages: an early phase [t < τ_f] during which the initial fluid pressure gradient within the fault is relaxed; a second transient stage [τ_f < t < τ_F] when fluid is rapidly expelled at one end of the fault and extracted from the surrounding rocks at the other end leading to a reduction in the pressure gradient in the intact rock; a third phase [t < τ_F] characterized by a steady‐state flow. From the numerical experiments we derive an expression for the steady‐state maximum fluid velocity in the fault and the values of the two timescales, τ_f and τ_F. A comparison indicates excellent agreement of our results with existing asymptotic solutions. For km‐scale faults, the model results suggest that steady‐state is unlikely to be reached over geological timescales. Thus, the current use of parameters such as the focusing ratio defined under the assumption of steady‐state conditions should be reconsidered.     

Episodic fluid flow within continental rift basins: some insights from field data and mathematical models of the Rhinegraben

Vitrinite reflectance data from a petroleum exploration well in the northern Upper Rhinegraben show an unusual vertical maturity trend. Above and below a 500 m thick marl layer the vitrinite reflectance levels are consistent with modern, conductive, geothermal gradients. Between about 1000 and 1500 m depth, however, vitrinite reflectance levels are significantly elevated (about 0.6%Ro). This anomaly cannot be explained with one‐dimensional conductive or conductive–convective heat transfer models, and thermal effects of sedimentation or igneous intrusion seem implausible for this geological setting. The thermal anomaly that formed this maturation anomaly must have been hydrothermal in origin, two‐dimensional in nature, and persisted long enough to elevate the vitrinite reflectance values within this marl unit, yet it must have dissipated before the thermal perturbation would have altered the organic matter below and above the unit. In this study, we propose that the vitrinite reflectance anomalies were caused by a transient thermal inversion induced by episodic, lateral flow of hot (130–160°C) groundwater along conductive fractures and bedding planes. Heat flow constraints suggest that fluids must have moved rapidly up a vertical feeder fault from a depth of at least 3.6 km before migrating laterally. To test this hypothesis, we present a suite of simple, idealized mathematical models of groundwater flow, heat transfer, thermal degradation of kerogen and vitrinite systematics to explore the episodic flow that could have produced the observed thermal anomaly. In these simulations, a single, horizontal aquifer is sandwiched between two less permeable units: the total dimensions of the vertical section model are 4 km thick by 10 km long. The top of the aquifer coincides with the position of the observed thermal maturity anomaly in the Rhinegraben. Boundary conditions along the left edge of this aquifer were varied through time to allow for the migration of hot fluids out into the basin. Inflow temperature, horizontal velocity, duration and frequency of flow and thickness of the aquifer were varied. We found that a thermal maturity anomaly could only be produced by a rather restrictive set of hydrothermal conditions. It was possible to produce the observed vitrinite reflectance anomaly by a single hydrothermal flow event of 130°C fluid migrating laterally into the aquifer at a rate of 1 m a^?1 for about 10 000 years. The anomaly is spatially confined to near the left edge of the basin, near the feeder fault. If the flow event lasted longer than 100 000 years, then the maturation anomaly disappeared as the lower confining unit approached steady‐state thermal conditions. It is possible that such an event occurred about 5 million years ago in response to increases in fault permeability associated with far field Alpine tectonism.     

Numerical modelling of 3D fluid flow and oxygen isotope exchange in fractured media: spatial distribution of isotope patterns

An understanding of fluid flow, mass transport and isotopic exchange in fractured rock is required to understand the origin of several geological processes including hydrothermal mineral deposits. The numerical model HydroGeoSphere simulates 3D advection, molecular diffusion, mechanical dispersion and isotopic exchange in a discretely fractured porous media, and can be used to better understand the processes of mass transport and isotopic exchange in fractured rocks. Study of ¹⁸O isopleth patterns for different types of fractures and fracture networks with a range of structural complexity and hydraulic properties shows that fracture properties and geometry control mass transport and isotopic exchange. The hydraulic properties, as well as the density, spacing, and connectivity of fractures determine the isotopic patterns. Asymmetries in the geometry of oxygen isotope patterns could be used to determine the direction of hydrothermal fluid flow.     

Fracture mineralization and fluid flow evolution: an example from Ordovician–Devonian carbonates,southwestern Ontario,Canada

Petrography, geochemistry (stable and radiogenic isotopes), and fluid inclusion microthermometry of matrix dolomite, fracture‐filling calcite, and saddle dolomite in Ordovician to Devonian carbonates from southwestern Ontario, Canada, provide useful insights into fluid flow evolution during diagenesis. The calculated δ¹⁸O_fluid, ΣREE, and REE_SN patterns of matrix and saddle dolomite suggest diverse fluids were involved in dolomitization and/or recrystallization of dolomite. The ⁸⁷Sr/⁸⁶Sr ratios of dolomite of each succession vary from values in the range of coeval seawater to values more radiogenic than corresponding seawater, which indicate diagenetic fluids were influenced by significant water/rock interaction. High salinities (22.4–26.3 wt. % NaCl + CaCl₂) of Silurian and Ordovician dolomite–hosted fluid inclusions indicate involvement of saline waters from dissolution of Silurian evaporites. High fluid inclusion homogenization temperatures (>100°C) in all samples from Devonian to Ordovician show temperatures higher than maximum burial (60–90°C) of their host strata and suggest involvement of hydrothermal fluids in precipitation and/or recrystallization of dolomite. A thermal anomaly over the mid‐continent rift during Devonian to Mississippian time likely was the source of excess heat in the basin. Thermal buoyancy resulting from this anomaly was the driving force for migration of hydrothermal fluids through regional aquifers from the center of the Michigan Basin toward its margin. The decreasing trend of homogenization temperatures from the basin center toward its margin further supports the interpreted migration of hydrothermal fluids from the basin center toward its margin. Hydrocarbon‐bearing fluid inclusions in late‐stage Devonian to Ordovician calcite cements with high homogenization temperatures (>80°C) and their ¹³C‐depleted values (approaching ?32‰ PDB) indicate the close relationship between hydrothermal fluids and hydrocarbon migration.     

Deep geothermal groundwater flow in the Seferihisar–Balçova area,Turkey: results from transient numerical simulations of coupled fluid flow and heat transport processes

The Seferihisar–Balçova Geothermal system (SBG) is characterized by complex temperature and hydrochemical anomalies. Previous geophysical and hydrochemical investigations suggest that hydrothermal convection in the faulted areas of the SBG and recharge flow from the Horst may be responsible for the observed patterns. A numerical model of coupled fluid flow and heat transport processes has been built in order to study the possible fluid dynamics of deep geothermal groundwater flow in the SBG. The results support the hypothesis derived from interpreted data. The simulated scenarios provide a better understanding of the geophysical conditions under which the different fluid dynamics develop. When recharge processes are weak, the convective patterns in the faults can expand to surrounding reservoir units or below the seafloor. These fault‐induced drag forces can cause natural seawater intrusion. In the Melange of the Seferihisar Horst, the regional flow is modified by buoyant‐driven flow focused in the series of vertical faults. As a result, the main groundwater divide can shift. Sealing caprocks prevent fault‐induced cells from being overwhelmed by vigorous regional flow. In this case, over‐pressured, blind geothermal reservoirs form below the caprocks. Transient results showed that the front of rising hot waters in faults is unstable: the tip of the hydrothermal plumes can split and lead to periodical temperature oscillations. This phenomenon known as Taylor–Saffman fingering has been described in mid‐ocean ridge hydrothermal systems. Our findings suggest that this type of thermal pulsing can also develop in active, faulted geothermal systems. To some extent, the role of an impervious fault core on the flow patterns has been investigated. Although it is not possible to reproduce basin‐scale transport processes, this first attempt to model deep groundwater geothermal flow in the SBG qualitatively supported the interpreted data and described the different fluid dynamics of the basin. Geofluids (2010) 10 , 388–405     

Evolution of porosity, permeability and fluid pressure in dilatant faults post-failure: implications for fluid flow and mineralization

Mineralization of brittle fault zones is associated with sudden dilation, and the corresponding changes in porosity, permeability and fluid pressure, that occur during fault slip events. The resulting fluid pressure gradients cause fluid to flow into and along the fault until it is sealed. The volume of fluid that can pass through the deforming region depends on the degree of dilation, the porosity and permeability of the fault and wall rocks, and the rate of fault sealing. A numerical model representing a steep fault cutting through a horizontal seal is used to investigate patterns of fluid flow following a dilatant fault slip event. The model is initialized with porosity, permeability and fluid pressure representing the static mechanical state of the system immediately after such an event. Fault sealing is represented by a specified evolution of porosity, coupled to changes in permeability and fluid pressure, with the rate of porosity reduction being constrained by independent estimates of the rate of fault sealing by pressure solution. The general pattern of fluid flow predicted by the model is of initial flow into the fault from all directions, followed by upward flow driven by overpressure beneath the seal. The integrated fluid flux through the fault after a single failure event is insufficient to account for observed mineralization in faults; mineralization would require multiple fault slip events. Downward flow is predicted if the wall rocks below the seal are less permeable than those above. This phenomenon could at least partially explain the occurrence of uranium deposits in reactivated basement faults that cross an unconformity between relatively impermeable basement and overlying sedimentary rocks.     

Fluids associated with hydrothermal dolomitization in St. George Group,western Newfoundland,Canada

Dolomite reservoirs are increasingly recognized as an important petroleum exploration target, although the application of a hydrothermal dolomite exploration model to these reservoirs remains controversial. The St. George Group of western Newfoundland consists of a sequence of dolomitised carbonates, with significant porosity development (up to 30%) and petroleum accumulations. Fluid inclusion microthermometry and bulk fluid leach analyses indicated that fluids responsible for matrix dolomitization (associated with intercrystalline porosity) and later saddle dolomitization are CaCl₂ ± MgCl₂ rich, high salinity (up to 26 eq. wt% NaCl) brines. Integration of fluid inclusion data with thermal maturation histories from the St. George Group show that these dolomites formed at temperatures higher than the ambient rock temperature, and are therefore hydrothermal in origin. Bulk leach analyses show that dolomitization is associated with influxes of postevaporitic brines (±Cl enriched magmatic fluids) late in the diagenetic history of these carbonates. This dolomitization is possibly Devonian in age, during a period of significant magmatic activity, extensional tectonics and development of hypersaline basins. Petrographic and geochemical similarities between Paleozoic hosted hydrothermal dolomitization in western Newfoundland, eastern Canada and the northeastern United States are consistent with a regional‐scale hydrothermal dolomitization event late in the diagenetic history of these carbonates. Geofluids (2010) 10 , 422–437     

Quartz recrystallization and fluid flow during contact metamorphism: a cathodoluminescence study

Cathodoluminescence (CL) images of quartz grains in the Appin Quartzite from the aureole of the Ballachulish Igneous Complex (Scotland) reveal a textural complexity that we interpret in the light of published models of the evolution of the contact aureole. Five distinct generations of quartz can be discriminated in CL. The oldest of these is a dark luminescing mottled quartz (Type 1 quartz) that occurs in the centres of pre‐existing grains, in samples collected from 210 m to 0.1 m from the contact. Dark mottled quartz is interpreted to be unrecrystallized material and has a regional metamorphic CL spectral signature. The onset of contact metamorphism resulted in grain growth visible in CL as a series of fine‐scale alternating bands of bright and dark luminescing material (Type 2 quartz), which we attribute to infiltration of repeated pulses of small amounts of H₂O along grain boundaries. Close to the intrusion, a subgrain‐scale network of intragranular, bright luminescing features could have resulted from either intragranular microcrack‐controlled infiltration of H₂O at high temperatures or intergranular cracking followed by grain growth (Type 3 quartz). Broad bands of bright material on grain boundaries in samples that are inferred to have undergone partial melting are interpreted as quartz crystallized from the melt phase (Type 4 quartz). The final stage in the textural development is marked by a series of aligned fractures, detected in CL by nonluminescing material (Type 5 quartz) and corresponding closely with trails of fluid inclusions. These fractures are interpreted as the pathways for late‐stage, low‐temperature, retrogressive fluids.     

Evolution of deformation and fault‐related fluid flow within an ancient methane seep system (Eocene,Varna, Bulgaria)

Faults are often important in fuelling methane seep systems; however, little is known on how different components in fault zones control subsurface fluid circulation paths and how they evolve through time. This study provides insight into fault‐related fluid flow systems that operated in the shallow subsurface of an ancient methane seep system. The Pobiti Kamani area (NE Bulgaria) encloses a well‐exposed, fault‐related seep system in unconsolidated Lower Eocene sandy deposits of the Dikilitash Formation. The Beloslav quarry and Beloslav N faults displace the Dikilitash Formation and are typified by broad, up to 80 m wide, preferentially lithified hanging wall damage zones, crosscut by deformation bands and deformation band zones, smaller slip planes and fault‐related joints. The formation of a shallow plumbing system and chimney‐like concretions in the Dikilitash Formation was followed by at least two phases of fault‐related methane fluid migration. Widespread fluid circulation through the Dikilitash sands caused massive cementation of the entire damage zones in the fault hanging walls. During this phase, paths of ascending methane fluids were locally obstructed by decimetre‐thick, continuous deformation band zones that developed in the partly lithified sands upon the onset of deformation. Once the entire damage zone was pervasively cemented, deformation proceeded through the formation of slip planes and joints. This created a new network of more localized conduits in close vicinity to the main fault plane and around through‐going slip planes. ¹³C‐depleted crustiform calcite cements in several joints record the last phase of focused methane fluid ascent. Their formation predated Neogene uplift and later meteoric water infiltration along the joint network. This illustrates how fault‐related fluid pathways evolved, over time, from ‘plumes’ in unconsolidated sediments above damage zones, leading to chimney fields, over widespread fluid paths, deflected by early deformation structures, to localized paths along fracture networks near the main fault.     

Inferring intermediate‐scale fluid flow in a heterogeneous metasedimentary multilayer sequence during progressive deformation: evidence from the Monts d’Arrée slate belt (Brittany,France)

Quartz veins in the early Variscan Monts d’Arrée slate belt (Central Armorican Terrane, Western France), have been used to determine fluid‐flow characteristics. A combination of a detailed structural analysis, fluid inclusion microthermometry and stable isotope analyses provides insights in the scale of fluid flow and the water–rock interactions. This research suggests that fluids were expelled during progressive deformation and underwent an evolution in fluid chemistry because of changing redox conditions. Seven quartz‐vein generations were identified in the metasedimentary multilayer sequence of the Upper Silurian to Lower Devonian Plougastel Formation, and placed within the time frame of the deformation history. Fluid inclusion data of primary inclusions in syn‐ to post‐tectonic vein generations indicate a gradual increase in methane content of the aqueous–gaseous H₂O–CO₂–NaCl–CH₄–N₂ fluid during similar P–T conditions (350–400°C and 2–3.5 kbar). The heterogeneous centimetre‐ to metre‐scale multilayer sequence of quartzites and phyllites has a range of oxygen‐isotope values (8.0–14.1‰ Vienna Standard Mean Ocean Water), which is comparable with the range in the crosscutting quartz veins (10.5–14.7‰ V‐SMOW). Significant differences between oxygen‐isotope values of veins and adjacent host rock (Δ = ?2.8‰ to +4.9‰ V‐SMOW) suggest an absence of host‐rock buffering on a centimetre scale, but based on the similar range of isotope values in the Plougastel Formation, an intraformational buffering and an intermediate‐scale fluid‐flow system could be inferred. The abundance of veins, their well‐distributed and isolated occurrence, and their direct relationship with the progressive deformation suggests that the intermediate‐scale fluid‐flow system primarily occurred in a dynamically generated network of temporarily open fractures.     

Computing geochemical mass transfer and water/rock ratios in submarine hydrothermal systems: implications for estimating the vigour of convection

A method of calculating chemical water/rock ratios is presented that enables the estimation of fluid velocities in open, flow‐through hydrologic systems. The approach is based on relating the gain/loss of a chemical species per kilogram of solid phase to the loss/gain of that species in the fluid phase, integrated across a specified length of the flowpath. After examining the underlying approximations of the approach using a one‐dimensional model of seawater moving through a basalt under nonisothermal conditions, the method is applied to representative zones within a two‐dimensional hydrothermal convective system. The method requires that regions within the flow system can be identified in which the direction of flow is steady for an extended period of time. Estimates of fluid velocity are spatial and temporal averages for the length of the flowpath used in the calculation. The location within the flow system and the nature of the alteration reactions determine which species can provide reliable values of the chemical water/rock ratio and useful estimates of fluid velocities. Over the length of the flowpath considered, the calculation of water/rock ratios works best when a species is controlled by a single reaction. Accurate estimates are favoured if the concentration profile of a species along the flowpath increases or decreases monotonically. If the length of the flowpath extends over more than one reaction zone, then erroneous estimates of the water/rock ratio and fluid velocity are more likely. Model calculations suggest that the quartz/silica system should provide reliable estimates for fluid velocity under a wide range of temperature and flow conditions, in particular in those regions of a system at or near quartz equilibrium, so that the aqueous silica concentration is buffered by quartz and correlated with the temperature distribution.     

Seismic constraints on the effects of gas hydrate on sediment physical properties and fluid flow: a review

The formation of gas hydrates in marine sediments changes their physical properties and hence influences fluid flow. Here, we review seismic indicators of gas hydrates and relate these indicators to gas hydrate formation and fluid migration. Analyses of seismic data from sediments containing gas and gas hydrates in a variety of locations have shown that the characteristic bottom‐simulating reflector (BSR), which commonly marks the hydrate phase boundary is caused mainly by the presence of gas beneath the gas hydrate stability zone (GHSZ). The amplitude of the BSR is also dependent on the hydrate concentration and on the porosity of the sediment. The presence of gas hydrate alters the elastic properties of sediments, particularly if it cements sediment grains. However, multifrequency studies in various geological provinces show that any loss of reflectivity or blanking observed within the GHSZ is dependent on both the nature of the sediments and concentration of hydrate present. Gas beneath the BSR may cause amplitude anomalies and may result in bright spots and enhanced reflections. The presence of gas beneath the BSR is the primary cause of observed amplitude versus offset (AVO) anomalies, but the amplitude of these anomalies is also dependent on the amount of cementation brought by the gas hydrates within the GHSZ. Fluid migration appears to play an important role in the formation and dissociation of gas hydrates in both active and passive margin settings. Fluid migration in accretionary prisms influences hydrate accumulation and may therefore control the spatial distribution of BSRs. Fluid migration may influence also the type of hydrate formed by bringing thermogenic gas containing higher order hydrocarbons to the GHSZ from below. Fluid advection may cause local dissociation of gas hydrates by bringing heat from below, thus shifting the gas hydrate phase boundary. Fluid flow within the GHSZ is limited by the formation of hydrate in the pore space, which reduces the permeability of the sediment. Features such as pockmarks, acoustic masking and acoustic turbidity are indirect indicators of fluid flow and identification of these features in seismic sections within and beneath the GHSZ may also suggest the formation of gas hydrate.     

Unravelling the petroleum system by enhancing fluid migration paths in seismic data using a neural network based pattern recognition technique

Understanding the hydrocarbon migration system in the sub‐surface is a key aspect of oil and gas exploration. It is well known that conventional 3D seismic data contains information about hydrocarbon accumulations. Less known is the fact that 3D seismic data also contains information about hydrocarbon migration paths in the form of vertical noise trails. A method has been developed to highlight vertical noise trails in seismic data semi‐automatically, using assemblies of directive multi‐trace seismic attributes and neural network technology. The results of this detection method yield valuable information about the origin of hydrocarbons, about migration paths from source to prospect and about leakage or spillage from these prospects to shallow gas pockets or to the sea bed. Besides, the results reveal the sealing quality of faults, provide information on overpressure and whether prospects are charged or not. All these aspects are useful information for basin modelling studies and for an increased understanding of the petroleum system.     

Combining finite element and finite volume methods for efficient multiphase flow simulations in highly heterogeneous and structurally complex geologic media

The permeability of the Earth's crust commonly varies over many orders of magnitude. Flow velocity can range over several orders of magnitude in structures of interest that vary in scale from centimeters to kilometers. To accurately and efficiently model multiphase flow in geologic media, we introduce a fully conservative node‐centered finite volume method coupled with a Galerkin finite element method on an unstructured triangular grid with a complementary finite volume subgrid. The effectiveness of this approach is demonstrated by comparison with traditional solution methods and by multiphase flow simulations for heterogeneous permeability fields including complex geometries that produce transport parameters and lengths scales varying over four orders of magnitude.