Numerical modeling of the origin of calcite mineralization in the Refugio‐Carneros fault,Santa Barbara Basin,California

Many faults in active and exhumed hydrocarbon‐generating basins are characterized by thick deposits of carbonate fault cement of limited vertical and horizontal extent. Based on fluid inclusion and stable isotope characteristics, these deposits have been attributed to upward flow of formation water and hydrocarbons. The present study sought to test this hypothesis by using numerical reactive transport modeling to investigate the origin of calcite cements in the Refugio‐Carneros fault located on the northern flank of the Santa Barbara Basin of southern California. Previous research has shown this calcite to have low δ¹³C values of about ?40 to ?30‰PDB, suggesting that methane‐rich fluids ascended the fault and contributed carbon for the mineralization. Fluid inclusion homogenization temperatures of 80–125°C in the calcite indicate that the fluids also transported significant quantities of heat. Fluid inclusion salinities ranging from fresh water to seawater values and the proximity of the Refugio‐Carneros fault to a zone of groundwater recharge in the Santa Ynez Mountains suggest that calcite precipitation in the fault may have been induced by the oxidation of methane‐rich basinal fluids by infiltrating meteoric fluids descending steeply dipping sedimentary layers on the northern basin flank. This oxidation could have occurred via at least two different mixing scenarios. In the first, overpressures in the central part of the basin may have driven methane‐rich formation waters derived from the Monterey Formation northward toward the basin flanks where they mixed with meteoric water descending from the Santa Ynez Mountains and diverted upward through the Refugio‐Carneros fault. In the second scenario, methane‐rich fluids sourced from deeper Paleogene sediments would have been driven upward by overpressures generated in the fault zones because of deformation, pressure solution, and flow, and released during fault rupture, ultimately mixing with meteoric water at shallow depth. The models in the present study were designed to test this second scenario, and show that in order for the observed fluid inclusion temperatures to be reached within 200 m of the surface, moderate overpressures and high permeabilities were required in the fault zone. Sudden release of overpressure may have been triggered by earthquakes and led to transient pulses of accelerated fluid flow and heat transport along faults, most likely on the order of tens to hundreds of years in duration. While the models also showed that methane‐rich fluids ascending the Refugio‐Carneros fault could be oxidized by meteoric water traversing the Vaqueros Sandstone to form calcite, they raised doubts about whether the length of time and the number of fault pulses needed for mineralization by the fault overpressuring mechanism were too high given existing geologic constraints.     

An exhumed palaeo-hydrocarbon migration fairway in a faulted carrier system, Entrada Sandstone of SE Utah, USA

The Moab Anticline, east‐central Utah, is an exhumed hydrocarbon palaeo‐reservoir which was supplied by hydrocarbons that migrated from the Moab Fault up‐dip towards the crest of the structure beneath the regional seal of the Tidwell mudstone. Iron oxide reduction in porous, high permeability aeolian sandstones records the secondary migration of hydrocarbons, filling of traps against small sealing faults and spill pathways through the Middle Jurassic Entrada Sandstone. Hydrocarbons entered the Entrada Sandstone carrier system from bends and other leak points on the Moab Fault producing discrete zones of reduction that extend for up to 400 m from these leak points. They then migrated in focused stringers, 2–5 m in height, to produce accumulations on the crest of the anticline. Normal faults on the anticline were transient permeability barriers to hydrocarbon migration producing a series of small compartmentalized accumulations. Exsolution of CO₂ as local fault seals were breached resulted in calcite cementation on the up‐dip side of faults. Field observations on the distribution of iron oxide reduction and calcite cements within the anticline indicate that the advancing reduction fronts were affected neither by individual slip bands in damage zones around faults nor by small faults with sand: sand juxtapositions. Faults with larger throws produced either sand: mudstone juxtapositions or sand: sand contacts and fault zones with shale smears. Shale‐smeared fault zones provided seals to the reducing fluid which filled the structural traps to spill points.     

Capillary seal capacity of faults under hydrodynamic conditions

Many fault bound traps are underfilled despite the top seal capacity being secure. The hydrocarbon sealing performance of faults themselves can be compromised either by mechanical or capillary process. Capillary process can be important either due to juxtaposition or to fine‐grained clay or cataclastic material within the fault zone itself. There is debate about how important each of these mechanisms is over geological timescales of hydrocarbon trapping. Recent work has provided insights into fine‐tuning capillary‐related fault seal calibration methodologies. Over the last 15 years, vigorous scientific debate with multiple published laboratory experiments and modelling studies has led some researchers and industry technologists to theorise that for water‐wet conventional hydrocarbon reservoirs, the relative water permeability in the reservoir (towards the top of the hydrocarbon column) may become very small, but in practice never reach zero. While not advocating for either side in this debate, the importance of accounting for hydrodynamic conditions regardless of the capillary sealing mechanism is demonstrated. Additionally, it is noted that nonzero relative water permeability has implications on how a seal's capillary threshold pressure for the nonwetting hydrocarbon phase is estimated from field data. In the particular case where there are pressure differences between unproduced hydrocarbon reservoirs on either side of a fault, then the hydrocarbon saturation must be discontinuous across the fault. For hydrocarbon leakage to occur across the entire thickness of the fault zone, the hydrocarbon pressure must exceed the threshold pressure on the side of the fault zone with the highest formation water hydraulic head. This approach to estimating across‐fault pressure difference will result in an improved calibration data set used for predrill estimation of capillary fault seal capacity.     

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.     

The Tune Viking Ship Reconsidered

The Tune Viking ship has been a riddle for more than 150 years, since being found within a burial in the Oslo fjord area in 1867. It was long thought that the ship's freeboard was too low for it to have crossed the North Sea. Advances in documentation methods and a detailed study of the preserved parts of the ship have provided new data, and this article outlines a new proposal for how the ship looked when it was built in the early 10th century AD. The Tune ship is reinterpreted as a seagoing vessel, in no way inferior to the Oseberg or Gokstad Viking ships.     

Seismic evidence for fluid migration pathways from an overpressured system in the South China Sea

Overpressured systems and intense, anomalously hot fluid expulsion in the Yinggehai Basin of the South China Sea offer an opportunity to understand the history of fluid flow and the process of hydrocarbon accumulation in overpressured environments. Fluid migration pathways from overpressured compartments in the basin are largely controlled by the distribution of faults and fractures. Episodic opening of these faults are related to the dynamics of an overpressured system and tectonic movements during basin evolution. At the crests of diapiric structures, fluid expulsion is seismically imaged as chimney‐ or plume‐like features, low to middle seismic amplitudes, and intermittently chaotic and blank reflecting seismic facies. These fluid pathways are controlled by vertical faults, which commonly penetrate overpressured and overlying normally pressured zones. Fluid expulsion is also observed near the main faults, such as the No. 1 Fault at the north‐eastern margin of the basin. Investigation by sidescan sonar on onshore and offshore Hainan Island indicates that there are more than 100 gas seepages adjacent to the No. 1 Fault. Migration pathways in the diapiric structures are controlled by three types of fault and fracture. Penetrative faults formed by dextral strike‐slip movement of the Red River faults commonly occur in the centre of the diapirs, and may have been a triggering factor for the diapirism, and controlled their distribution. Hydrofractures occur in certain mud‐rich layers and may have been generated by hydraulic fracturing. Radial normal faults occur at the top of diapirs and were formed by the intrusive process. These fluid migration pathways played an important role in regional hydrocarbon accumulation.     

Effects of episodic fluid flow on hydrocarbon migration in the Newport‐Inglewood Fault Zone,Southern California

Fault permeability may vary through time due to tectonic deformations, transients in pore pressure and effective stress, and mineralization associated with water‐rock reactions. Time‐varying permeability will affect subsurface fluid migration rates and patterns of petroleum accumulation in densely faulted sedimentary basins such as those associated with the borderland basins of Southern California. This study explores the petroleum fluid dynamics of this migration. As a multiphase flow and petroleum migration case study on the role of faults, computational models for both episodic and continuous hydrocarbon migration are constructed to investigate large‐scale fluid flow and petroleum accumulation along a northern section of the Newport‐Inglewood fault zone in the Los Angeles basin, Southern California. The numerical code solves the governing equations for oil, water, and heat transport in heterogeneous and anisotropic geologic cross sections but neglects flow in the third dimension for practical applications. Our numerical results suggest that fault permeability and fluid pressure fluctuations are crucial factors for distributing hydrocarbon accumulations associated with fault zones, and they also play important roles in controlling the geologic timing for reservoir filling. Episodic flow appears to enhance hydrocarbon accumulation more strongly by enabling stepwise build‐up in oil saturation in adjacent sedimentary formations due to temporally high pore pressure and high permeability caused by periodic fault rupture. Under assumptions that fault permeability fluctuate within the range of 1–1000 millidarcys (10^?15–10^?12 m²) and fault pressures fluctuate within 10–80% of overpressure ratio, the estimated oil volume in the Inglewood oil field (approximately 450 million barrels oil equivalent) can be accumulated in about 24 000 years, assuming a seismically induced fluid flow event occurs every 2000 years. This episodic petroleum migration model could be more geologically important than a continuous‐flow model, when considering the observed patterns of hydrocarbons and seismically active tectonic setting of the Los Angeles basin.     

Simulation of the impact of faults on CO2 injection into sandstone reservoirs

Numerical simulations of multiphase CO₂ behavior within faulted sandstone reservoirs examine the impact of fractures and faults on CO₂ migration in potential subsurface injection systems. In southeastern Utah, some natural CO₂ reservoirs are breached and CO₂‐charged water flows to the surface along permeable damage zones adjacent to faults; in other sites, faulted sandstones form barriers to flow and large CO₂‐filled reservoirs result. These end‐members serve as the guides for our modeling, both at sites where nature offers ‘successful’ storage and at sites where leakage has occurred. We consider two end‐member fault types: low‐permeability faults dominated by deformation‐band networks and high‐permeability faults dominated by fracture networks in damage zones adjacent to clay‐rich gouge. Equivalent permeability (k) values for the fault zones can range from <10^?14 m² for deformation‐band‐dominated faults to >10^?12 m² for fracture‐dominated faults regardless of the permeability of unfaulted sandstone. Water–CO₂ fluid‐flow simulations model the injection of CO₂ into high‐k sandstone (5 × 10^?13 m²) with low‐k (5 × 10^?17 m²) or high‐k (5 × 10^?12 m²) fault zones that correspond to deformation‐band‐ or fracture‐dominated faults, respectively. After 500 days, CO₂ rises to produce an inverted cone of free and dissolved CO₂ that spreads laterally away from the injection well. Free CO₂ fills no more than 41% of the pore space behind the advancing CO₂ front, where dissolved CO₂ is at or near geochemical saturation. The low‐k fault zone exerts the greatest impact on the shape of the advancing CO₂ front and restricts the bulk of the dissolved and free CO₂ to the region upstream of the fault barrier. In the high‐k aquifer, the high‐k fault zone exerts a small influence on the shape of the advancing CO₂ front. We also model stacked reservoir seal pairs, and the fracture‐dominated fault acts as a vertical bypass, allowing upward movement of CO₂ into overlying strata. High‐permeability fault zones are important pathways for CO₂ to bypass unfaulted sandstone, which leads to reduce sequestration efficiency. Aquifer compartmentalization by low‐permeability fault barriers leads to improved storativity because the barriers restrict lateral CO₂ migration and maximize the volume and pressure of CO₂ that might be emplaced in each fault‐bound compartment. As much as a 3.5‐MPa pressure increase may develop in the injected reservoir in this model domain, which under certain conditions may lead to pressures close to the fracture pressure of the top seal.     

Migration of Mn-rich fluids through normal faults and fine-grained terrigenous sediments during early development of the Neogene Vallès-Penedès half-graben (NE Spain)

The Miocene siliciclastic sediments infilling the Vallès‐Penedès half‐graben are affected by two sets of structures developed during the extensional tectonics that created the basin. The first set, represented by extension fractures infilled with mud and sands, is attributed to seismically induced liquefaction. The second set, represented by normal faults, corresponds to a high‐permeability horsetail extensional fracture mesh developed near the surface in the hanging walls of normal faults. The incremental character of the vein‐fills indicates episodic changes in the tectonic stress state and fault zone permeability. Two episodes of fluid migration are recorded. The first episode occurred prior to consolidation and lithification when shallow burial conditions allowed oxidizing meteoric waters to flow horizontally through the more porous and permeable sandy layers. Development of clastic dikes allowed local upward flow and dewatering of the sandy beds. Liquefaction and expulsion of fluids were probably driven by seismic shaking. During the first episode of fluid migration there was no cementation of the sandstone or within the fractures, probably because little fluid was mobilized by the predominantly compaction‐driven flow regime. The second episode of fluid migration occurred synchronously with normal fault development, during which time the faults acted as fluid conduits. Fluids enriched in manganese, probably leached from local manganese oxyhydroxides soon after sedimentation, moved laterally and produced cementation in the sandstone layers, eventually arriving at the more porous and permeable fault pathways that connected compartments of different porosities and permeabilities. Carbonate probably precipitated in fractures saturated with meteoric water near the ground surface at a transitional redox potential. Once the faults became occluded by calcite cement, shortly after fault development, they became barriers to both vertical and horizontal fluid flow.     

Thermal convection in faulted extensional sedimentary basins: theoretical results from finite-element modeling

Thermal convection has the potential to be a significant and widespread mechanism of fluid flow, mass transport, and heat transport in rift and other extensional basins. Based on numerical simulation results, large‐scale convection can occur on the scale of the basin thickness, depending on the Rayleigh number for the basin. Our analysis indicates that for syn‐rift and early post‐rift settings with a basin thickness of 5 km, thermal convection can occur for basal heat flows ranging from 80 to 150 mW m^?2, when the vertical hydraulic conductivity is on the order of 1.5 m year^?1 and lower. The convection cells have characteristic wavelengths and flow patterns depending on the thermal and hydraulic boundary conditions. Steeply dipping extensional faults can provide pathways for vertical fluid flow across large thicknesses of basin sediments and can modify the dynamics of thermal convection. The presence of faults perturbs the thermal convective flow pattern and can constrain the size and locations of convection cells. Depending on the spacing of the faults and the hydraulic properties of the faults and basin sediments, the convection cells can be spatially organized to align with adjacent faults. A fault‐bounded cell occurs when one convection cell is constrained to occupy a fault block so that the up‐flow zone converges into one fault zone and the down‐flow zone is centred on the adjacent fault. A fault‐bounded cell pair occurs when two convection cells occupy a fault block with the up‐flow zone located between the faults and the down‐flow zones centred on the adjacent faults or with the reverse pattern of flow. Fault‐bounded cells and cell pairs can be referred to collectively as fault‐bounded convective flow. The flow paths in fault‐bounded convective flow can be lengthened significantly with respect to those of convection cells unperturbed by the presence of faults. The cell pattern and sense of circulation depend on the fault spacing, sediment and fault permeabilities, lithologic heterogeneity, and the basal heat flow. The presence of fault zones also extends the range of conditions for which thermal convection can occur to basin settings with Rayleigh numbers below the critical value for large‐scale convection to occur in a basin without faults. The widespread potential for the occurrence of thermal convection suggests that it may play a role in controlling geological processes in rift basins including the acquisition and deposition of metals by basin fluids, the distribution of diagenetic processes, the temperature field and heat flow, petroleum generation and migration, and the geochemical evolution of basin fluids. Fault‐bounded cells and cell pairs can focus mass and heat transport from longer flow paths into fault zones, and their discharge zones are a particularly favourable setting for the formation of sediment‐hosted ore deposits near the sea floor.     

U/Th-rich bitumen in Archean granites and Palaeoproterozoic metasediments, Rum Jungle Mineral Field, Australia: implications for mineralizing fluids

Uranium/thorium (U/Th)‐rich bitumen has been discovered within both Palaeoproterozoic black pelites and the Archean granitic basement of the Rum Jungle Mineral Field, Northern Territory, Australia. Granite‐hosted bitumen occurs as small (up to 400‐µm diameter) discrete individual nodules, which exhibit many morphological similarities to those observed in Phanerozoic siliciclastic rocks. Thorium, the dominant radioelement, occurs primarily as a hydrated Th–Y–Si–P phase. Uranium‐rich inclusions are rare, and correspond to a hydrated U–Th–Y–Si–P phase, identified as coffinite–thorogummite. Metasediment‐hosted bitumen is more variable in morphology, occurring as massive (<2 cm in width) veins that cross‐cut all foliations, as discrete individual nodules or as elongate seams (up to 500 µm in length), interpreted to represent a series of coalesced individual nodules. In all examples, uranium, the dominant radioelement in the metasediment‐hosted bitumen, is present as Th‐poor uraninite, with variable Y₂O₃ contents (up to 3.21 wt.%). Raman investigation of all types of bitumen indicates that it is a poorly organized carbonaceous matter, which has not been subjected to metamorphism. Consequently, a post‐metamorphic timing for hydrocarbon emplacement can be inferred and a magmatic origin can be precluded. Potential source rocks for the bitumen are black shales of the Whites Formation (up to 8 wt.% total organic carbon (TOC)) and the Koolpin Formation (approximately 13 wt.% TOC). Post‐metamorphic sericitization of rocks within the Whites Formation is accompanied by a near‐complete removal of organic matter. Alteration was possibly the catalyst for hydrocarbon generation. The Th–Y–Si–P phase within the granite‐hosted nodules is interpreted to be the result of the alteration of antecedent monazite. During this alteration, U, LREE and P were fractionated and removed, while Th, Y and Si remained immobile, and recombined to form a hydrated Th–Y–Si phase. This pervasive alteration within the basement U/Th‐rich granites is proposed as a genetic model for the formation of uranium deposits in the Rum Jungle Mineral Field and possibly unconformity associated uranium deposits on a global scale.     

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     

Downward hydrocarbon migration predicted from numerical modeling of fluid overpressure in the Paleozoic Anticosti Basin,eastern Canada

The Anticosti Basin is a large Paleozoic basin in eastern Canada where potential source and reservoir rocks have been identified but no economic hydrocarbon reservoirs have been found. Potential source rocks of the Upper Ordovician Macasty Formation overlie carbonates of the Middle Ordovician Mingan Formation, which are underlain by dolostones of the Lower Ordovician Romaine Formation. These carbonates have been subjected to dissolution and dolomitization and are potential hydrocarbon reservoirs. Numerical simulations of fluid‐overpressure development related to sediment compaction and hydrocarbon generation were carried out to investigate whether hydrocarbons generated in the Macasty Formation could migrate downward into the underlying Mingan and Romaine formations. The modeling results indicate that, in the central part of the basin, maximum fluid overpressures developed above the Macasty Formation due to rapid sedimentation. This overpressured core dissipated gradually with time, but the overpressure pattern (i.e. maximum overpressure above source rock) was maintained during the generation of oil and gas. The downward impelling force associated with fluid‐overpressure gradients in the central part of the basin was stronger than the buoyancy force for oil, whereas the buoyancy force for gas and for oil generated in the later stage of the basin is stronger than the overpressure‐related force. Based on these results, it is proposed that oil generated from the Macasty Formation in the central part of the basin first moved downward into the Mingan and Romaine formations, and then migrated laterally up‐dip toward the basin margin, whereas gas throughout the basin and oil generated in the northern part of the basin generally moved upward. Consequently, gas reservoirs are predicted to occur in the upper part of the basin, whereas oil reservoirs are more likely to be found in the strata below the source rocks. Geofluids (2010) 10 , 334–350     

Potential of palaeofluid analysis for understanding oil charge history

Fluid inclusion data, particularly the distribution of hydrocarbon fluid inclusions and their chemistry, can provide insights into oil charge in a petroleum-prospective region. Examples from the UK Atlantic margin show how we can understand thermal regime, timing and chemistry of oil charge. Data from the UK Atlantic margin based on fluid inclusion temperature profiles shows anomalously high temperatures which are highest at the top of the Triassic–Eocene sequence. This is interpreted as a product of hot fluid flow, probably reflecting hydrothermal activity related to intrusion of sills at depth. The preservation of high temperatures also implies rapid migration from depth through fracture systems. Ar–Ar analysis of oil-bearing K-feldspar cements, and petrographic studies of oil inclusion distribution help delimit timing and migration pathways for the hot fluid charge and later fluid migration events. Coupled with compositional data for oils measured destructively (organic geochemistry) or non-destructively (fluorescence), these approaches allow the development of oil charge histories based on real data rather than theoretical modelling.     

Faults and fault properties in hydrocarbon flow models

The petroleum industry uses subsurface flow models for two principal purposes: to model the flow of hydrocarbons into traps over geological time, and to simulate the production of hydrocarbon from reservoirs over periods of decades or less. Faults, which are three-dimensional volumes, are approximated in both modelling applications as planar membranes onto which predictions of the most important fault-related flow properties are mapped. Faults in porous clastic reservoirs are generally baffles or barriers to flow and the relevant flow properties are therefore very different to those which are important in conductive fracture flow systems. A critical review and discussion is offered on the work-flows used to predict and model capillary threshold pressure for exploration fault seal analysis and fault transmissibility multipliers for production simulation, and of the data from which the predictions derive. New flow simulation models confirm that failure of intra-reservoir sealing faults can occur during a reservoir depressurization via a water-drive mechanism, but contrary to anecdotal reports, published examples of production-induced seal failure are elusive. Ignoring the three-dimensional structure of fault zones can sometimes have a significant influence on production-related flow, and a series of models illustrating flow associated with relay zones are discussed.     

Communication between overpressured compartments

Overpressure build‐up in compartments, and communication between overpressured compartments across faults are studied with simple analytical and numerical models. It is shown that the excess pressure in a (vertical) one‐dimensional, one‐compartment model can be written as the sum of the excess pressure generated in the seal above the compartment, and a second part, which is due to the expulsion of fluid from the compartment and the rocks below. The one‐compartment model is generalized to a two‐compartment model, which accounts for the fluid communication between the compartments through a fault zone. The volume rates of flow through the seals and the fault zone are shown to be the weighted mean of the volume rates of the one‐dimensional, one‐compartment model. The normalized weights are given by dimensionless numbers, called fault–seal numbers, which control the communication between the compartments. A fault–seal number much less than unity implies that the fault is a stronger barrier for the fluid flow than the seal. A fault–seal number larger than unity implies the opposite: that the seal is a stronger barrier than the fault. The conditions for isolated compartments and other regimes are identified in terms of the fault–seal numbers. It is discussed how the compartment fault–seal numbers can be computed when the permeability is given in the fault zone. The results given by the analytical compartment models are demonstrated and validated with two‐dimensional numerical (finite element) simulations.     

Quantifying secondary migration efficiencies

Exploration success relies on properly risking the hydrocarbon system relevant for each prospect. Accurate risking of secondary migration efficiencies has been difficult due to lack of simple procedures that relate rock properties such as permeability and entry pressures to migration velocities, oil stringer heights and saturations. In order to achieve improved estimates of charge probabilities, equations for the secondary migration process are formulated based upon the Darcy flow and buoyancy conditions. An analytical solution of the formulated equations is shown, making it possible to construct charts for efficiently assessing the column height of secondary migration hydrocarbon stringers. The average oil (hydrocarbon) saturation of the migrating stringer can be computed, making it easy to compute the permeability related, secondary migration losses. Inputs to the chart are hydrocarbon flow‐rates and flow‐path width, hydrocarbon viscosity and density, carrier bed dip, permeability and entry pressures. Outputs are stringer heights, hydrocarbon saturation, relative permeability, migration velocities and migration losses. A procedure for including the new equations into existing basin scale fluid flow simulators is outlined and a Java applet for calculating the properties is described. The Java applet is useful for sensitivity studies, and can also be used to test results from basin simulators with the new migration efficiency equations. The analytical solution suggests that many published methods for calculating hydrocarbon migration in fluid flow simulators will over‐estimate hydrocarbon saturations and therefore losses. Calculated migration velocities will also be too low.     

Geothermometry and geobarometry of overpressured environments in Qiongdongnan Basin,South China Sea*

We demonstrate the use of PVT fluid inclusion modelling in the calculation of palaeofluid formation pressures, using samples from the YC21‐1‐1 and YC21‐1‐4 wells in the YC21‐1 structural closure, Qiongdongnan Basin, South China Sea. Homogenisation temperatures and gas/liquid ratios were measured in aqueous fluid inclusions, and associated light hydrocarbon/CO₂‐bearing inclusions, and their compositions were determined using a crushing technique. The vtflinc software was used to construct P–T phase diagrams that enabled derivation of the minimum trapping pressure for each order of fluid inclusion. Through the projection of average homogenisation temperatures (155, 185.5 and 204.5°C) for three orders of fluid inclusion on the thermal‐burial history diagram of the Oligocene Yacheng and Lingshui formations, their trapping times were constrained at 4.3, 2.1 and 1.8 Ma, respectively. The formation pressure coefficient, the ratio of fluid pressure/hydrostatic pressure established by PVT modelling coupled with DST data, demonstrates that one and a half cycles of pressure increase–discharge developed in the Yacheng and Lingshui formations for about 4.3 Ma. In comparison, the residual formation pressure determined by 2D numerical modelling in the centre of LeDong depression shows two and a half pressure increase–discharge cycles for about 28 Ma. The two different methods suggest that a high fluid potential in the Oligocene reservoir of the YC21‐1 structure developed at two critical stages for regional oil and natural gas migration and accumulation (5.8 and 2.0 Ma, respectively). Natural gas exploration in this area is therefore not advisable.     

What do aqueous geothermometers really tell us?

The application of chemical geothermometry to shallow groundwaters or spring discharge assumes that there is minimal mixing or re-equilibration of water as it travels from depth to the surface. In this study, we examine the potential for mixing and re-equilibration by examining heat and fluid flow along crustal-scale faults in tectonic geothermal systems. Numerical modeling results indicate that maximum in situ temperatures could be under-predicted by up to 30% due to mixing of fluids that enter the fault at different depths. This, coupled with the depression of isotherms by downward groundwater flow in the hanging wall, could cause underestimates of maximum circulation depth of greater than 80% in extreme cases. Kinetics does not favor re-equilibration in the shallower portions of faults due to low temperatures and higher fluid velocities. However, in areas of deeper circulation or higher heat flow such reactions are possible.     

Impacts of Pleistocene glacial loading on abnormal pore‐water pressure in the eastern Michigan Basin

The hydromechanical effects of Pleistocene glacial loading on the Michigan Basin are assessed using numerical analysis based on coupled stress and pore‐water pressure. The two‐dimensional model domain included the Basin cross section and extended 10 km into the Precambrian. In the analysis, we considered the effects of the number of glacial loading cycles, the presence and connectedness of a deep Cambrian aquifer, the direction of glacial advance, the effect of a wet versus dry glacier/soil interface, topographic effects, density‐driven flow effects, and lithosphere flexure on the development of anomalous pressures. The modeling results were compared with data collected from highly instrumented wells completed in the eastern margin of the Basin. The present‐day results define regions of significant underpressure in the upper Ordovician and lower Silurian formations characterized by very low hydraulic conductivity and regions of overpressure where hydraulic conductivity is higher. To achieve the degree of underpressure observed in the instrumented wells using the model, a specific loading cycle applied over 100 000 years was repeated four times. As the number of loading cycles increased, the Paleozoic formations reached a state where the underpressures remain constant. Our results also illustrate the difference in poroelastic modeling between the application of mechanical loads on the land surface and the application of an equivalent hydraulic head, where the latter developed overpressures rather than the observed underpressures. The modeling also shows that the overpressures observed in the Cambrian formations are most likely to be the result of density‐driven flow. Finally, the simulations show that the effects of lithosphere flexure in the hydromechanical models results in the development of lateral stresses that generate overpressures rather than underpressures in the southern half of the domain. As there are no suitable observation points, these results remain unconfirmed, and further study is warranted.