Reprint of: Maturity modelling integrated with apatite fission-track dating: Implications for the thermal history of the Mid-Polish Trough (Poland)

The Mid-Polish Trough (MPT) is situated in the easternmost part of the Central European Basin System (CEBS) and stretches NW–SE across the Polish Basin. It was characterised by pronounced subsidence and thick sediment accumulation between the Permian and the Late Cretaceous. Late Cretaceous–early Paleogene basin inversion led to the formation of the Mid-Polish Swell (MPS). The study area is located within the Pomeranian segment of the MPT/MPS (NW Poland) and experienced up to 7 km Permian-Mesozoic subsidence. PetroMod 1-D modelling was performed on several well-sections in order to study Permian to recent burial-uplift evolution. The modelling was calibrated with new vitrinite reflectance (VR_r) data and allowed to constrain the magnitude of uplift and related erosion as well as provided a first overview of the temperature history. The base of the studied Permian–Mesozoic successions attained maximum burial depths of 4800–5400 m before the onset of the inversion, less than in the axial trough area. The thickness of pre- and most probably also syn-inversion Upper Cretaceous deposits is estimated as 300 m. Erosion associated with inversion processes removed between 900 and 1400 m of the Mesozoic sediments, i.e. 1000–1500 m less than in the most inverted central part of the trough. VR_r data suggest constant Permian–Mesozoic heat flows corresponding to present-day values (40–45 mW/m²). Apatite fission-track (AFT) ages modelled with the PetroMod module PetroTracks show a good fit with AFT ages directly measured on well samples, and further support the assumption of steady heat flow in the range 40–45 mW/m². Palaeotemperatures appear to have decreased towards the East European Craton margin, which is compatible with the present day distribution of heat flow. Thermal history modelling shows a relatively simple Permian–Mesozoic heat flow pattern in the Pomeranian segment of the MPT. Such a scenario implies that the present-day heat flow distribution has not changed essentially since Mesozoic times.     

Evidence of a large upper-Cretaceous depocentre across the Continent-Ocean boundary of the Congo-Angola basin. Implications for palaeo-drainage and potential ultra-deep source rocks

The analysis of 2D deep-seismic-reflection profiles across the slope and abyssal plain of the Angola oceanic basin reveals the existence of a significant and formerly unknown depocentre beneath the giant Cenozoic Congo deep-sea fan, between 7000 m and 9000 m depth, deposited directly onto the Aptian oceanic crust. The unit, which is up to 2.5 km thick and extends for more than 200 km basinwards of the Continent-Ocean boundary, is probably aged Albian–Turonian. Its radial fan-shaped depocentre is centred on the present-day Congo River outlet and contains at least 0.2 Mio km³ of sediments. These observations and the results from flexural modelling indicate that (1) the location of the Congo River's outlet has remained fairly stable since the Late Cretaceous, and (2) the basal unit was indeed sourced by a palaeo-Congo River probably located nearby the present-day one. Thus, the Atlantic sedimentary system related to the exoreism of the Congo River is much older than previously thought. Thermal modelling indicates that the maturation history of this upper-Cretaceous deposits is highly influenced by the interaction between the initial high heat flow of the young oceanic crust and further increase in sediment supply due to the progradation of the overlying Tertiary deep-sea fan during the Miocene. Hence, despite low present-day heat-flow values, should the unit have source rock potential, its basal section may be currently generating hydrocarbons.     

Source rock characteristics and hydrocarbon generation modelling of Upper Cretaceous Mukalla Formation in the Jiza-Qamar Basin,Eastern Yemen

The Upper Cretaceous Mukalla coals and other organic-rich sediments which are widely exposed in the Jiza-Qamar Basin and believed to be a major source rocks, were analysed using organic geochemistry and petrology. The total organic carbon (TOC) contents of the Mukalla source rocks range from 0.72 to 79.90% with an average TOC value of 21.50%. The coals and coaly shale sediments are relatively higher in organic richness, consistent with source rocks generative potential. The samples analysed have vitrinite reflectance in the range of 0.84–1.10 %Ro and pyrolysis T_max in the range of 432–454 °C indicate that the Mukalla source rocks contain mature to late mature organic matter. Good oil-generating potential is anticipated from the coals and coaly shale sediments with high hydrogen indices (250–449 mg HC/g TOC). This is supported by their significant amounts of oil-liptinite macerals are present in these coals and coaly shale sediments and Py-GC (S₂) pyrograms with n-alkane/alkene doublets extending beyond nC₃₀. The shales are dominated by Type III kerogen (HI < 200 mg HC/g TOC), and are thus considered to be gas-prone.One-dimensional basin modelling was performed to analysis the hydrocarbon generation and expulsion history of the Mukalla source rocks in the Jiza-Qamar Basin based on the reconstruction of the burial/thermal maturity histories in order to improve our understanding of the of hydrocarbon generation potential of the Mukalla source rocks. Calibration of the model with measured vitrinite reflectance (Ro) and borehole temperature data indicates that the present-day heat flow in the Jiza-Qamar Basin varies from 45.0 mW/m² to 70.0 mW/m² and the paleo-heat flow increased from 80 Ma to 25 Ma, reached a peak heat-flow values of approximately 70.0 mW/m² at 25 Ma and then decreased exponentially from 25 Ma to present-day. The peak paleo-heat flow is explained by the Gulf of Aden and Red Sea Tertiary rifting during Oligocene-Middle Miocene, which has a considerable influence on the thermal maturity of the Mukalla source rocks. The source rocks of the Mukalla Formation are presently in a stage of oil and condensate generation with maturity from 0.50% to 1.10% Ro. Oil generation (0.5% Ro) in the Mukalla source rocks began from about 61 Ma to 54 Ma and the peak hydrocarbon generation (1.0% Ro) occurred approximately from 25 Ma to 20 Ma. The modelled hydrocarbon expulsion evolution suggested that the timing of hydrocarbon expulsion from the Mukalla source rocks began from 15 Ma to present-day.     

Combining 2D-basin and structural modelling to constrain heat transport along the Muroto Transect,Nankai Trough,Japan

The Nankai Trough located southeast of Shikoku Island, Japan, exhibits a zone of exceptionally high heat flow. In the central part of the Nankai Trough the fossil spreading centre of the Shikoku Basin is subducted beneath the southwest Japan arc. We have modelled the temperature and maturation history along the Muroto Transect reaching from the tip of the thrust zone out into nearly undeformed Quaternary and Tertiary sediments seawards of Nankai Trough. We used two balanced cross-sections defining the sections before and after overthrusting as input for 2D-basin modelling. We can show that rapid burial and overthrusting during the Quaternary in combination with a heat flow history following the cooling curve of a 15 Ma old oceanic plate is not sufficient to explain the measured maturity of organic material in the sediments. Several heat flow scenarios derived from theoretical concepts [Yamano, M., Kinoshita, M., Goto, S., Matsubayashi, O., 2003. Extremely high heat flow anomaly in the middle part of the Nankai Trough. Physics and Chemistry of the Earth, Parts A/B/C 28, 487–497.] and previous modelling approaches [e.g. Brown, K.M., Saffer, D.M., Bekins, B.A., 2001. Smectite diagenesis, pore water freshening, and fluid flow at the toe of the Nankai wedge. Earth and Planetary Science Letters 194, 97–109; Spinelli, G.A., Underwood, M.B., 2005. Modeling thermal history of subducting crust in Nankai Trough: constraints from in situ sediment temperature and diagenetic reaction progress. Geophysical Research Letters 32(L09301): doi:10.1029/2005GL022793; Steurer, J., Underwood, M.B., 2003. Clay mineralogy of mudstones from the Nankai Trough reference sites 1173 and 1177 and frontal accretionary prism site 1174. In: H. Mikada et al. (Eds.), pp. 1–37. Available from: <http://www-odp.tamu.edu/publications/190196SR/VOLUME/CHAPTERS/211.PDF>] were tested. The best match between observed maturity levels, temperature and heat flow measurements is reached for a heat flow history which initially assumes the cooling of a 15 Ma old oceanic lithosphere but is reheated to 170–180 mW/m² during the phase of rapid burial in the Quaternary. This can be achieved either by assuming the onset of hydrothermal circulation in the cooling crust or by reheating caused by off-axis volcanism at about 6 Ma [Yamano, M., Kinoshita, M., Goto, S., Matsubayashi, O., 2003. Extremely high heat flow anomaly in the middle part of the Nankai Trough. Physics and Chemistry of the Earth, Parts A/B/C 28, 487–497.].     

pT-effects of Pleistocene glacial periods on permafrost,gas hydrate stability zones and reservoir of the Mittelplate oil field,northern Germany

During the past two million years low surface temperatures as well as episodically advancing ice sheets from Scandinavia acted on the subsurface pT-regime of northern Germany. Their likely effects on the petroleum system of Schleswig-Holstein were investigated. For the entire Quaternary mean annual ground temperature (MAGT) was reconstructed at a resolution of 1000 years by calibrating oxygen isotope records from ODP-site 659 to the climate of northern Germany of the past 120 kyr. The resulting MAGT trend served as input to an ice sheet model and a permafrost model along a 2D section crossing the petroleum bearing south-western part of Schleswig-Holstein. Here advances and retreats of the Scandinavian ice sheet during Saalian and Elsterian glaciation Stages were reconstructed. Maximum ice thicknesses of up to 1700 m and up to 20 periods of regional permafrost in northern Germany were reconstructed for the past 1.25 million years. Based on a basal heat flow of 50 mW/m² permafrost thicknesses exceeded 100 m during most of these periods, temporarily extending down to depths of more than 300 m. Favourable surface temperatures and long durations of cold periods provided favourable conditions for onshore gas hydrate stability zones at Mittelplate. Implementing these glacial dynamics into 2D basin modelling (PetroMod, IES, Aachen, Germany) of the Mittelplate oil field reveals five phases of gas hydrate stability at depths down to 750 m. The latest of these events occurred during the Weichselian about 20 kyr ago. The effect of the ice sheets on pore pressure in the subsurface strongly depends on the hydraulic boundary conditions at the ice base (e.g. frozen vs. temperate ice sheet base). Excess pore pressure in the reservoir of more than 10 MPa during ice overriding is possible and probable. The calculated temperature effect of the Pleistocene cooling on the Mittelplate reservoir is in the range of 3–7 °C. Even today temperature in the reservoir is still lowered by about 4 °C in comparison to pre-Pleistocene times. Despite the fact that a significant influence of glacial effects on petroleum generation can be ruled out at Mittelplate, we state that pT-effects in reservoirs related to glacial processes in formerly glaciated areas have been underestimated in the past.     

Latest Caledonian to Present tectonomorphological development of southern Norway

The regional resultant stress field of the northeastern North Atlantic has shifted significantly throughout the Phanerozoic. In Fennoscandian parts of the Caledonian orogen, mountain building, which was characterized by NW-SE contraction (reference to present North), was followed by a collapse with transport both parallel and transverse to the mountain chain. The Late Palaeozoic – Mesozoic saw several stages of E-W to NW-SE extension, varying in time and position. Local episodes of inversion are traceable in some cases, particularly in connection with deep-seated and long-lived zones of weakness. The Cenozoic has to a larger degree been affected by compression, including folding and basin inversion. Again some of the more pronounced effects of local inversion are related to pre-existing fault systems. Neogene uplift of the western mountainous area in Scandinavia can be unravelled by potential field study, AFT data and reflection seismic sections. Assuming that the region is close to isostatic equilibrium, the uplifted areas must be supported at depth by substantial volumes of low-density material within the crust or the mantle, close to the crust/mantle interface or close to the lithosphere/asthenosphere interfaces.     

Basin modeling of the Late Miocene Zeit source rock in the Sudanese portion of Red Sea Basin: Implication for hydrocarbon generation and expulsion history

The Late Miocene Zeit Formation is exposed in the Red Sea Basin of Sudan and represents an important oil-source rock. In this study, five (5) exploratory wells along Red Sea Basin of Sudan are used to model the petroleum generation and expulsion history of the Zeit Formation. Burial/thermal models illustrate that the Red Sea is an extensional rift basin and initially developed during the Late Eocene to Oligocene. Heat flow models show that the present-day heat flow values in the area are between 60 and 109 mW/m². The variation in values of the heat flow can be linked to the raise in the geothermal gradient from margins of the basin towards offshore basin. The offshore basin is an axial area with thick burial depth, which is the principal heat flow source.The paleo-heat flow values of the basin are approximately from 95 to 260 mW/m², increased from Oligocene to Early Pliocene and then decreased exponentially prior to Late Pliocene. This high paleo-heat flow had a considerable effect on the source rock maturation and cooking of the organic matter. The maturity history models indicate that the Zeit Formation source rock passed the late oil-window and converted the oil generated to gas during the Late Miocene.The basin models also indicate that the petroleum was expelled from the Zeit source rock during the Late Miocene (>7 Ma) and it continues to present-day, with transformation ratio of more than 50%. Therefore, the Zeit Formation acts as an effective source rock where significant amounts of petroleum are expected to be generated in the Red Sea Basin.     

Variations of methane induced pyrite formation in the accretionary wedge sediments offshore southwestern Taiwan

The accretionary wedge of offshore southwestern Taiwan contains abundant deposits of gas hydrate beneath the sea floor. High concentrations of methane in pore waters are observed at several locations with little data concerning historical methane venting available. To understand temporal variation of methane venting in sediments over geologic time, a 23-m-long Calypso piston core (MD05-2911) was collected on the flank of the Yung-An Ridge. Pore water sulfate, dissolved sulfide, dissolved iron, methane, sedimentary pyrite, acid volatile sulfide, reactive iron, organic carbon and nitrogen as well as carbonate δ¹³C were analyzed.Three zones with markedly different pyrite concentration were found at the study site. Unit I sediments (>20 mbsf) were characterized with a high amount of pyrite (251–380 μmol/g) and a δ¹³C-depleted carbonate, Unit II sediments (15–20 mbsf) with a low pyrite (15–43 μmol/g) and a high content of iron oxide mineral and Unit III sediments (<10 mbsf) by a present-day sulfate–methane interface (SMI) at 5 m with a high amount of pyrite (84–221 μmol/g) and a high concentration of dissolved sulfide.The oscillation records of pyrite concentrations are controlled by temporal variations of methane flux. With an abundant supply of methane to Unit I and III, anaerobic methane oxidation and associated sulfate reduction favor diagenetic conditions conducive for significant pyrite formation. No AOM signal was found in Unit II, characterized by typical organically-limited normal marine sediments with little pyrite formation. The AOM induced pyrite formation near the SMI generates a marked pyrite signature, rendering such formation of pyrite as a useful proxy in identifying methane flux oscillation in a methane flux fluctuate environment.     

A ten-year data set for fetch- and depth-limited wave growth

This paper presents the key results from a ten-year data set for Lake IJssel and Lake Sloten in The Netherlands, containing information on wind, storm surges and waves, supplemented with SWAN 40.51 wave model results. The wind speeds U₁₀, effective fetches x and water depths d for the data set ranged from 0–24 m s^− 1, 0.8–25 km and 1.2–6 m respectively. For locations with non-sloping bottoms, the range in non-dimensional fetch x? ( = gxU₁₀^− 2) was about 25–80,000, while the range in dimensionless depth d? ( = g d U₁₀^− 2) was about 0.03–1.7. Land–water wind speed differences were much smaller than the roughness differences would suggest. Part of this seems due to thermal stability effects, which even play a role during near-gale force winds. For storm surges, a spectral response analysis showed that Lake IJssel has several resonant peaks at time scales of order 1 h. As for the waves, wave steepnesses and dimensionless wave heights H? ( = gH_m0U₁₀^− 2) agreed reasonably well with parametric growth curves, although there is no single curve to which the present data fit best for all cases. For strongly depth-limited waves, the extreme values of d? (0.03) and H_m0 / d (0.44) at the 1.7 m deep Lake Sloten were very close to the extremes found in Lake George, Australia. For the 5 m deep Lake IJssel, values of H_m0 / d were higher than the depth-limited asymptotes of parametric wave growth curves. The wave model test cases of this study demonstrated that SWAN underestimates H_m0 for depth-limited waves and that spectral details (enhanced peak, secondary humps) were not well reproduced from H_m0 / d = 0.2–0.3 on. SWAN also underestimated the quick wave response (within 0.3–1 h) to sudden wind increases. For the remaining cases, the new [Van der Westhuysen, A.J., Zijlema, M., and Battjes, J.A., 2007. Nonlinear saturation-based whitecapping dissipation in SWAN for deep and shallow water, Coast. Eng., 54, 151–170] SWAN physics yielded better results than the standard physics of Komen, G.J., Hasselmann, S., Hasselmann, K., 1984. On the existence of a fully developed wind-sea spectrum. J. Phys. Oceanogr. 14, 1271–1285, except for persistent overestimations that were found for short fetches. The present data set contains many interesting cases for detailed model validation and for further studies into the evolution of wind waves in shallow lakes.     

Particle dynamics of Be,Pb and the implications of sedimentation of heavy metals in the Wenjiao/Wenchang and Wanquan River estuaries,Hainan, China

Radionuclides (i.e., ⁷Be and ²¹⁰Pb) can be used to trace particle and sediment dynamics and to quantify coastal oceanic processes with time scales ranging from a few days to a hundred years. Here, we study the settling dynamics of suspended particles and the implication by sedimentary heavy metals in the Wenjiao/Wenchang River and Wanquan River estuaries through the measurement of the particulate ⁷Be and ²¹⁰Pb nuclides. Activity in the particulate phase had a range of 2.1–54.5 and 4.6–67.9 Bq kg⁻¹ for ⁷Be and excess ²¹⁰Pb (²¹⁰Pb_xs), respectively, in the Wenjiao/Wenchang River estuary. In the Wanquan River estuary, activity is in the range of 1.2–43.5 Bq kg⁻¹ for ⁷Be and 6.2–194.5 Bq kg⁻¹ for ²¹⁰Pb_xs. At the same time, activity in the dissolved phase had a range of 0.46–1.26 and 0.30–1.17 Bq m⁻³ for ⁷Be and ²¹⁰Pb, respectively, in the Wenjiao/Wenchang River estuary; ranges of 0.10–2.31 and 0.09–1.87 Bq m⁻³ for ⁷Be and ²¹⁰Pb, respectively, were observed in the Wanquan River estuary. The distribution coefficients (K_d) for the two nuclides decreased within increased in suspended particle matters (SPM) concentration and/or salinity in Wanquan River estuary. The residence times of particulate ⁷Be and ²¹⁰Pb_xs had ranges of 0.4–1.6 and 1.65–5.15 days, respectively, in the Wenjiao/Wenchang River estuary; and ranges of 0.02–3.2 and 0.61–4.44 days, respectively in the Wanquan River estuary. All residence times for the two nuclides increased in the seaward direction. In the Wenjiao/Wenchang River estuary, we found that 11.8–21.0% of Cu, 3.0–9.0% of Zn and 43.2–69.9% for Cd is removed from the water column and deposited into the estuary, and 24.2–34.8% for Cu, 7.2–23.8% for Zn, and 70.0–82.5% for Cd in the Wanquan River estuary, respectively.     

Sediment mounds and other sedimentary features related to hydrate occurrences in a columnar seismic blanking zone of the Ulleung Basin, East Sea, Korea

A mound related to a cold vent in a columnar seismic blanking zone (CSBZ) was formed around site UBGH1-10 in the central Ulleung Basin (2077 m water depth), East Sea, Korea. The mound is 300–400 m wide and 2–3 m high according to multi-beam bathymetry, 2–7 kHz sub-bottom profiler data, and multi-channel reflection seismic data. Seafloor topography and characteristics were investigated using a remotely operated vehicle (ROV) around site UBGH1-10, which is located near the northern part of the mound. The origin of the mound was investigated through lithology, mineralogy, hydrate occurrence, and sedimentary features using dive cores, piston cores, and a deep-drilling core. The CSBZ extends to ∼265 ms two-way traveltime (TWT) below the seafloor within a mass-transport deposit (MTD) unit. Gas hydrate was entirely contained 6–141 m below the seafloor (mbsf) within hemipelagic deposits intercalated with a fine-grained turbidite (HTD) unit, characteristically associated with high resistivity values at site UBGH1-10. The hydrate is commonly characterized by veins, nodules, and massive types, and is found within muddy sediments as a fracture-filling type. Methane has been produced by microbial reduction of CO₂, as indicated by C₁/C₂₊, δ¹³CCH4, and δD₄CH analyses. The bowl-shaped hydrate cap revealed at 20–45 ms TWT below the seafloor has very high resistivity and high salinity, suggesting rapid and recent gas hydrate formation. The origin of the sediment mound is interpreted as a topographic high formed by the expansion associated with the formation of the gas hydrate cap above the CSBZ. The lower sedimentation rate of the mound sediments may be due to local enhancement of bottom currents by topographic effects. In addition, no evidence of gas bubbles, chemosynthetic communities, or bacterial mats was observed in the mound, suggesting an inactive cold vent.     

New heat flow measurements in the Ulleung Basin, East Sea (Sea of Japan): relationship to local BSR depth, and implications for regional heat flow distribution

In July 2007, new marine heat flow data were collected at ten sites (HF01–10) in the central and southwestern sectors of the Ulleung Basin (East Sea or Sea of Japan) as part of regional gas hydrate research. In addition, cores were collected at five of these sites for laboratory analysis. The results show that the geothermal gradient ranged from 103–137 mK/m, and the in-situ thermal conductivity from 0.82–0.95 W/m·K. Laboratory measurements of thermal conductivity were found to deviate by as much as 40% from the in-situ measurements, despite the precautions taken to preserve the cores. Based on the in-situ conductivity, the heat flow was found to increase with water depth toward the center of the basin, ranging from 84–130 mW/m². Using a simple model, we estimated the heat flow from the depths of the BSR, and compared this with the observed heat flow. In our study area, the two sets of values were quite consistent, the observed heat flows being slightly higher than the BSR-derived ones. The evaluation of regional pre-1994 data revealed that the heat flow varied widely from 51–157 mW/m² in and around the basin. Due to a large scatter in these older data, a clear relationship between heat flow and water depth was not evident, in contrast to what would be expected for a rifted sedimentary basin. This raises the question as to whether the pre-1994 data represent the true background heat flow from the underlying basin crust since the basin opening, and/or whether they contain large measurement errors. In fact, evidence in support of the latter explanation exists. BSRs are generally found in the deep parts of the basin, and vary by only ±15 m in depth below the seafloor. From the average BSR depth, we inferred the background heat flow using a simple model, which in the case of the Ulleung Basin is approximately 120 and 80 mW/m² for 2.5 and 1 km below sea level, respectively.     

Charging of the Penglai 9-1 oil field,Bohai Bay basin,China: Functions of the delta on accumulating petroleum

The Penglai 9-1 (PL9-1) oil field, which contains China's third largest offshore oil accumulation (in-place reserves greater than 2.28 × 10⁸ ton or 1.49 × 10⁹ bbl), was found in shallow reservoirs (700–1700 m, 2297–5577 ft) within the most active fault zone in east China. The PL9-1 field contains two oil-bearing series, the granite intrusions in Mesozoic (Mz) and both the sandstone reservoirs in Neogene Guantao (Ng) and Neogene Minghuazhen (Nm) Formation. The origins of the PL9-1 field, both in terms of source rock intervals and generative kitchens, were determined by analyzing biomarker distributions for 61 source rock samples and 33 oil samples. The Mesozoic granite intrusions, which hold more than 80% of the oil reserves in the field, were charged in the west by oil generated from the third member (Es₃) of the Shahejie Formation in the Bodong depression. The Neogene reservoirs of the PL9-1 field were charged in the west by oil generated from the third member (Es₃) of the Shahejie Formation in the Bodong depression and in the south by oil generated from the first member (Es₁) of the Shahejie Formation in the Miaoxibei depression. Interactive contact between the large fan delta and the mature source rocks residing in the Es₃ Formation of the Bodong depression resulted in a high expulsion efficiency from the source rocks and rapid oil accumulation in the PL9-1 field, which probably explains how can this large oil field accumulate and preserve within the largest and most active fault zone in east China.     

Effect of sedimentary heterogeneities on hydrocarbon accumulations in the Permian Shanxi Formation,Ordos Basin,China: Insight from an integrated stratigraphic forward and petroleum system modelling

Sedimentary heterogeneities are ubiquitous in nature and occur over a range of scales from core, reservoir to basin scales. They may thus exert significant influences on hydrocarbon generation, migration and accumulation. The sedimentary heterogeneities of the Permian Shanxi Formation in the Ordos Basin, China were modelled using Sedsim, a stratigraphic forward modelling program. The simulation results were then used to construct a 3D petroleum system model using PetroMod. The effects of sedimentary heterogeneities on hydrocarbon accumulations were evaluated by comparing the integrated Sedsim-PetroMod model with the classic 3D basin model. The Sedsim simulation shows that considerable sedimentary heterogeneities are present within the Shanxi Formation, as a result of the interplay of the initial topography, tectonic subsidence, base level change and sediment inputs. A variety of lithologies were developed both laterally and vertically within the Shanxi Formation at kilometre and metre scales, respectively, with mudstones mainly developed in the depositional centre, while sandstones developed in the southern and northern margin areas. A typical source-ward retrogradation is well developed within the Lower Shanxi Formation.A base-case classic 3D basin model was constructed to quantify the Permian petroleum system in the Ordos Basin. The geological and thermal models were calibrated using Vr and borehole temperature data. The source rocks of the Upper Paleozoic became mature (R_o > 0.5%) and high mature (R_o > 1.2%) in the late Triassic and late Jurassic, respectively, in the central and southern areas. During the Early Cretaceous, a tectonically induced geothermal event occurred in the southern Ordos Basin. This caused the source rocks to reach over maturity (R_o > 2.0%) quite rapidly in the early Late Cretaceous in the central and southern areas. All the source rock transformation ratios (TR) at present are greater than 70% in the P₁ coal and P₁ mudstone layers with TR values approaching 100% in the central and southern areas. The transformation ratios of the P₁ limestone are close to 100% over the entire interval.In the base-case model, a large amount of hydrocarbons appear to have been expelled and migrated into the Shanxi Formation, but only a minor amount was accumulated to form reservoirs. In the model, the Shanxi Formation sandstone layer was set to be homogeneous vertically and there was no regional seal rocks present at the top of the Shanxi Formation. Therefore hydrocarbons could not be trapped effectively with only minor accumulations in some local structural highs where hydrocarbons are trapped both at the top and in the up-dip direction by the adjacent mudstone facies. In contrast, the integrated Sedsim-PetroMod model takes into account of the internal lithological and sedimentary facies heterogeneities within the Shanxi Formation, forming complex contiguous sandstone-mudstone stacking patterns. Hydrocarbons were found to have accumulated in multiple intervals of lithological traps within the Shanxi Formation. The results indicate that lithological distinctions, controlled by sedimentary heterogeneities in three dimensions can provide effective sealing in both the top and up-dip directions for hydrocarbon accumulations, with gas being mainly accumulated near the depocentre where lithological traps usually formed due to frequent oscillations of the lake level.     

Heat flow derived from BSR and its implications for gas hydrate stability zone in Shenhu Area of northern South China Sea

Bottom simulating reflectors (BSRs), known as the base of gas hydrate stability zone, have been recognized and mapped using good quality three-dimensional (3D) pre-stack migration seismic data in Shenhu Area of northern South China Sea. Additionally, seismic attribute technique has been applied to better constrain on the distribution of gas hydrate. The results demonstrate that gas hydrate is characterized by “blank” zone (low amplitude) in instantaneous amplitude attribute. The thickness of gas hydrate stability zone inferred from BSR ranges from 125 to 355 m with an average of 240 m at sea water depth from 950 to 1,600 m in this new gas hydrate province. The volume of gas in-place bound in hydrate is estimated from 1.7 × 10⁹ to 4.8 × 10⁹ m³, with the most likely value of around 3.3 × 10⁹ m³, using Monte Carlo simulation. Furthermore, geothermal gradient and heat flow are derived from the depths of BSRs using a conductive heat transfer model. The geothermal gradient varies from 35 to 95°C km⁻¹ with an average of 54°C km⁻¹. Corresponding heat flow values range from 43 to 105 mW m⁻² with an average of 64 mW m⁻². By comparison with geological characteristics, we suggest that the distribution of gas hydrate and heat flow are largely associated with gas chimneys and faults, which are extensively distributed in Shenhu Area, providing easy pathways for fluids migrating into the gas hydrate stability zone for the formation of gas hydrate. This study can place useful constraints for modeling gas hydrate stability zone from measured heat flow data and understanding the mechanism of gas hydrate formation in Shenhu Area.     

The thermohaline structure and evolution of the deep waters in the Canada Basin,Arctic Ocean

Below the sill depth (at about 2400 m) of the Alpha-Mendeleyev ridge complex, the waters of the Canada Basin (CB) of the Arctic Ocean are isolated, with a ¹⁴C isolation age of about 500 yr. The potential temperature θ decreases with depth to a minimum θ_m≈−0.524°C near 2400 m, increases with depth through an approximately 300 m thick transition layer to θ_h≈−0.514°C, and then remains uniform from about 2700 m to the bottom at 3200–4000 m. The salinity increases monotonically with depth through the deep θ_m and transition layer from about 34.952 to about 34.956 and then remains uniform in the bottom layer. A striking staircase structure, suggestive of double-diffusive convection, is observed within the transition layer. The staircase structure is observed for about 1000 km across the basin and has been persistent for more than a decade. It is characterized by 2–3 mixed layers (10–60 m thick) separated by 2–16 m thick interfaces. Standard formulae, based on temperature and salinity jumps, suggest a double-diffusive heat flux through the staircase of about 40 mW m⁻², consistent with the measured geothermal heat flux of 40–60 mW m⁻². This is to be expected for a scenario with no deep-water renewal at present as we also show that changes in the bottom layer are too small to account for more than a small fraction of the geothermal heat flux. On the other hand, the observed interfaces between mixed layers in the staircase are too thick to support the required double-diffusive heat flux, either by molecular conduction or by turbulent mixing, as there is no evidence of sufficiently vigorous overturns within the interfaces. It therefore seems, that while the staircase structure may be maintained by a very weak heat flux, most of the geothermal heat flux is escaping through regions of the basin near lateral boundaries, where the staircase structure is not observed. The vertical eddy diffusivity required in these near-boundary regions is O(10⁻³) m² s⁻¹. This implies Thorpe scales of order 10 m. We observe what may be Thorpe scales of this magnitude in boundary-region potential temperature profiles, but cannot tell if they are compensated by salinity. The weak stratification of the transition layer means that the large vertical mixing rate implies a local dissipation rate of only O(10⁻¹⁰) W kg⁻¹, which is not ruled out by plausible energy budgets. In addition, we discuss an alternative scenario of slow, continuous renewal of the CB deep water. In this scenario, we find that some of the geothermal heat flux is required to heat the new water and vertical fluxes through the transition layer are reduced.     

Pressure cells and pressure seals in the UK Central Graben

The Central Graben of the North Sea is characterised by high levels of overpressure (up to 40 MPa overpressure at 4500 m depth). We present pressure data for Cenozoic and Mesozoic reservoirs. Palaeocene sandstones control pressures in Tertiary mudstones and Cretaceous Chalk by acting as a regional ‘drain’. We divide the Jurassic into 18 pressure cells. The rift structure of the Graben controls the magnitude of pressure in each cell. Lateral hydraulic communication exists over 10 km distance between deeply-buried terraces (> 5000 m depth) and shallow structural highs (< 4500 m depth). Lateral communication increases pressure in the structurally-elevated sandstones to the minimum stress. This dynamic process produces zones of vertical fluid flow on the Forties-Montrose High, termed Leak Points. Vertical flow at Leak Points produces a 20 MWm⁻² heat flow anomaly and controls hydrocarbon retention. Leak Points are water-wet, while deep terraces in hydraulic communication with Leak Points are condensate-bearing. The Kimmeridge Clay Fm. forms the pressure seal in deep terraces.     

Simultaneous solution for water,heat and salt balances in a Mediterranean coastal lagoon (Mar Menor,Spain)

A modelling approach is proposed to evaluate the environmental dynamics of coastal lagoons. The water, heat and salt balances are addressed simultaneously, providing a better estimation of evaporation and water exchanges. Compared to traditional approaches, the model presented accounts for the effects of water salinity, heat storage and net energy advected in the water body. The model was applied daily to the Mar Menor coastal lagoon (SE Spain) from 2003 through 2006. Water exchanges with the Mediterranean Sea were estimated based on the monthly trend of the lagoon salinity and were correlated with monthly averages of wind speed. The mean daily water exchange with the sea was 1.77 hm³ d⁻¹. This exchange accounted for only 1% of the heat losses in the lagoon heat balance, and it is the most important flow in the water balance. The mean annual evaporation flux amounted to 101.3 W m⁻² (3.55 mm d⁻¹), while the sensible heat flux amounted to 19.7 W m⁻², leading to an annual Bowen ratio on the order of 0.19. To validate the model, daily water temperatures were predicted based on the daily heat balance of the water body and were compared with remote sensing data from water surface standard products.     

Temporal changes of fault seal and early charge of the Maui Gas-condensate field,Taranaki Basin,New Zealand

Fault seal due to juxtaposition or the generation of low-permeability fault rock has the potential to change through time with displacement accumulation. Temporal variations in cross-fault flow of hydrocarbons have been assessed for the Cape Egmont Fault (CEF), Taranaki Basin New Zealand, using displacement backstripping, juxtaposition and Shale Gouge Ratio (SGR) analysis. The timing of hydrocarbon migration and charge of the giant Maui Gas-condensate Field across the CEF have been assessed using seismic reflection lines (2D & 3D), coherency cubes, VShale curves from the Maui-2 well and PetroMod modelling. Displacement–backstripping analysis suggests that between the Late Miocene and early Pleistocene (5.5 and 2.1 Ma) sandstone reservoir units of the Maui Field (Mangahewa, Kaimiro and Farewell Formations) and underlying source rocks (Rakopi Formation) were partly juxtaposed across the CEF with low SGRs (< 0.2) present in the fault zone. Following 2.1 Ma SGRs increased to 0.2–0.55 adjacent to the Eocene–Palaeocene reservoir succession which was not in juxtaposed contact with source rocks. PetroMod modelling using these SGR values and juxtaposition relationships supports cross-fault flow prior to 2.1 Ma with later charge across the fault being less likely. Gas chimneys and the gas–water contact in the Eocene reservoir proximal to the fault suggest that despite limited cross-fault flow, upward leakage of hydrocarbons from the reservoir occurred after 2.1 Ma, possibly associated with active fault movement or fracturing related to faulting, and may account for the loss of an early oil phase.