Correction to: Multilevel Graph Partitioning for Three-Dimensional Discrete Fracture Network Flow Simulations

Mathematical Geosciences - The publication of this article unfortunately contained a mistake. The assignment to the affiliations of author Satish Karra was not correct     

Investigation into the use of 2D elastic waveform inversion from look‐ahead walk‐away VSP surveys

Pore pressure in sediments beneath salt in the Gulf of Mexico varies widely creating a potentially dangerous and difficult drilling challenge. Estimating elastic parameters of sediments beneath salt is key to the prediction of pore pressure and reducing the drilling risk in exiting the base of the salt. In this paper we investigate the ability of 2D waveform inversion to recover the elastic parameters in the sedimentary layer beneath the salt from a walk‐away VSP (vertical seismic profile) carried out with the receivers in the salt, with the objective of estimating pore pressure at the base of the salt (to be estimated using traditional methods). We propose an effective method for performing the inversion and apply this method to a blind test of a large and realistic synthetic dataset. To facilitate the design of a VSP survey suitable to this type of inversion we also present an analysis into the effects of the receiver array receiver spacing. It is shown that the resulting velocity estimates are sufficiently accurate to predict the pore pressure within the limits required for drilling.     

Mechanical decoupling and thermal structure at the East Pacific Rise axis 9°N: Constraints from axial magma chamber geometry and seafloor structures

We study the relationships between the seafloor structures and the axial magma chamber geometry in the 9°N overlapping spreading center (OSC) area on the fast spreading East Pacific Rise (EPR). Our observations are based on a new high resolution bathymetric map of the 9°N OSC area derived from picks of the seafloor arrival on 3D seismic data, and on previously published data that constrain the presence and distribution of melt below the 9°N OSC. Differences in the orientation of structures between the seafloor and the magma chamber indicate a sharp change in principal stress directions with depth, suggesting that the brittle crust above the melt sill is decoupled from the melt sill itself and the ductile crust underlying it. The stress-field within the brittle upper crust results from a local interaction of the two overlapping spreading centers, whereas the stress-field in the crust below the melt sill corresponds to the regional stress-field imposed by plate separation. Given this mechanical structure of the crust, the melt sill shape and location appear to be controlled by the following factors: the location of the deep melt source below the melt sill, the ambient stress-field at the depth of the melt sill, and the stress-field in the brittle upper crust above the melt sill, which thermally shapes the roof of the melt sill through repeated eruptions.     

Structural imaging of Mesozoic sediments of Kachchh,India, and their hydrocarbon prospects

Travel-time modelling of first and second arrivals has been used to produce P-wave velocity models of the shallow sub-surface structure of Kachchh, India, using new wide-angle reflection/refraction profiles. Results obtained from the 2-D tomographic modelling have also been corroborated with magnetotelluric and borehole data within the proximity of these profiles. Based on multiple geophysical data, a composite sub-surface model is presented in this article. The prominent features of the model include the presence of more than 3 km thick Mesozoic sediments residing above the Precambrian basement. Its correlation with available lithological knowledge shows that the Mesozoic rocks have been deposited from early Triassic to late Cretaceous until the eruption of the Deccan Flood basalts at the Cretaceous-Tertiary boundary. The new results suggest that the Mesozoic rocks of Kachchh can be sub-divided in two major parts which correspond to late and early Mesozoic deposits respectively. The sub-surface models successfully demarcate the regional variations in the basaltic layer across the region and show maximum thickness of 1.2 km. The modelling results also exhibit that the Precambrian basement in this region varies between 4 and 6 km depth. The new data provide a much clearer picture than was hitherto available of the entire sedimentary succession in the basins that underlie the Kachchh region. They show a significant variability in thickness of the various sedimentary layers, and in particular illustrate a marked topographic irregularity at the base of the sedimentary succession.     

Estimates of sediment and nutrient loads in 10 major catchments draining to the Great Barrier Reef during 2006-2009

The Great Barrier Reef (GBR) catchment area has been monitored simultaneously for sediment and nutrient exports from 10 priority catchments discharging into the GBR lagoon between 2006 and 2009. This allows GBR catchment-wide exports to be estimated and spatially compared within a discrete time-frame. Elevated levels of sediment and nutrient exports were recorded in all monitored catchments as compared to pre-European estimates, but vary around previous estimates of mean annual loads. During the period of monitoring, the Burdekin and Fitzroy catchments contributed the highest sediment and nutrient exports, however when loads were normalised for area, these catchments produced the lowest unit yields. In contrast, the highest yields were produced in the wetter and proportionately more intensively cultivated Johnstone, O'Connell, and Pioneer catchments particularly for dissolved nitrogens. This assessment offers the necessary scientific foundation for future monitoring, assessment, and management of sediment and nutrient loads entering the GBR.     

Early Carboniferous HP metamorphism in the Hida Gaien Belt,Japan: Implications for the Palaeozoic tectonic history of proto‐Japan

We report two new eclogite localities (at Kanayamadani and Shinadani) in the high‐P (HP) metamorphic rocks of the Omi area in the western most region of Niigata Prefecture, Japan, which form part of the Hida Gaien Belt, and determine metamorphic conditions and pressure–temperature (P–T) paths. The metamorphic evolution of the eclogites is characterized by a tight hairpin‐shaped P–T path from prograde epidote–blueschist facies to peak eclogite facies and then retrograde blueschist facies. The prograde metamorphic stage is characterized by various amphibole (winchite, barroisite, glaucophane) inclusions in garnet, whereas the peak eclogite facies assemblage is characterized by omphacite, garnet, phengite and rutile. Peak P–T conditions of the eclogites were estimated to be ~600°C and up to 2.0 GPa by conventional cation‐exchange thermobarometry, Ti‐in‐zircon thermometry and quartz inclusion Raman barometry respectively. However, the Raman spectra of carbonaceous material thermometry of metapelites associated with the eclogites gave lower peak temperatures, possibly due to metamorphism at different conditions before being brought together during exhumation. The blueschist facies overprint following the peak of metamorphism is recognized by the abundance of glaucophane in the matrix. Zircon grains in blueschist facies metasedimentary samples from two localities adjacent to the eclogites have distinct oscillatory‐zoned cores and overgrowth rims. Laser ablation inductively coupled plasma mass spectrometry U–Pb ages of the detrital cores yield a wide range between 3,200 and 400 Ma, with a peak at 600–400 Ma. In the early Palaeozoic, proto‐Japan was located along the continental margin of the South China craton, providing the source of the older population of detrital zircon grains (3,200–600 Ma) deposited in the trench‐fill sediments. In addition, subduction‐related magmatism c. 500–400 Ma is recorded in the crust below proto‐Japan, which might have been the source for the younger detrital zircon grains. The peak metamorphic age was constrained by SHRIMP dating of the overgrowth rims, yielding Tournaisian ages of 347 ± 4 Ma, suggesting subduction in the early Carboniferous. Our results provide clear constraints on the initiation of subduction, accretion and the development of an arc‐trench system along the active continental margin of the South China craton and help to unravel the Palaeozoic tectonic history of proto‐Japan.     

Sub‐sample time shift and horizontal displacement measurements using phase‐correlation method in time‐lapse seismic

Hydrocarbon production and fluid injection affect the level of subsurface stress and physical properties of the subsurface, and can cause reservoir‐related issues, such as compaction and subsidence. Monitoring of oil and gas reservoirs is therefore crucial. Time‐lapse seismic is used to monitor reservoirs and provide evidence of saturation and pressure changes within the reservoir. However, relative to background velocities and reflector depths, the time‐lapse changes in velocity and geomechanical properties are typically small between consecutive surveys. These changes can be measured by using apparent displacement between migrated images obtained from recorded data of multiple time‐lapse surveys. Apparent displacement measurements by using the classical cross‐correlation method are poorly resolved. Here, we propose the use of a phase‐correlation method, which has been developed in satellite imaging for sub‐pixel registration of the images, to overcome the limitations of cross‐correlation. Phase correlation provides both vertical and horizontal displacements with a much better resolution. After testing the method on synthetic data, we apply it to a real dataset from the Norne oil field and show that the phase‐correlation method can indeed provide better resolution.