An Automatic Well Planner for Complex Well Trajectories

Mathematical Geosciences - A data-driven automatic well planner procedure is implemented to develop complex well trajectories by efficiently adapting to near-well reservoir properties and geometry....     

Development and validation of a finite element morphological model for shallow water basins

Recognising the importance of understanding sediment dynamics to evaluate the status of a coastal lagoon environment, this work has been focused on the investigation of the hydrodynamic and sediment transport processes occurring in such basins. In order to describe the lagoon system, a modelling approach combining hydrodynamics, waves and sediment dynamics has been developed. The framework of the numerical model consists of a finite element hydrodynamic model, a third generation finite element spectral wave model and a sediment transport and morphodynamic model for both cohesive and non-cohesive sediments. The model adopts the finite element technique for spatial integration, which has the advantage to describe more accurately complicated bathymetry and irregular boundaries for shallow water areas. The developed model has been applied to test cases and to a very shallow tidal lagoon, the Venice Lagoon, Italy. Numerical results show good agreement with water level, waves and turbidity measurements collected in several monitoring stations inside the Lagoon of Venice. Such a model represents an indispensable tool in analysing coastal problems and assessing morphological impacts of human interference.     

Fluid flow and pore pressure development throughout the evolution of a trough mouth fan,western Barents Sea

Using a combination of geophysical and geotechnical data from Storfjorden Trough Mouth Fan off southern Svalbard, we investigate the hydrogeology of the continental margin and how this is affected by Quaternary glacial advances and retreats over the continental shelf. The geotechnical results show that plumites, deposited during the deglaciation, have high porosities, permeabilities and compressibilities with respect to glacigenic debris flows and tills. These results together with margin stratigraphic models obtained from seismic reflection data were used as input for numerical finite element models to understand focusing of interstitial fluids on glaciated continental margins. The modelled evolution of the Storfjorden TMF shows that tills formed on the shelf following the onset of glacial sedimentation (ca. 1.5 Ma) acted as aquitards and therefore played a significant role in decreasing the vertical fluid flow towards the sea floor and diverting it towards the slope. The model shows that high overpressure ratios (up to λ ca. 0.6) developed below the shelf edge and on the middle slope. A more detailed model for the last 220 kyrs accounting for ice loading during glacial maxima shows that the formation of these aquitards on the shelf focused fluid flow towards the most permeable plumite sediments on the slope. The less permeable glacigenic debris flows that were deposited during glacial maxima on the slope hinder fluid evacuation from plumites allowing high overpressure ratios (up to λ ca. 0.7) to develop in the shallowest plumite layers. These high overpressures likely persist to the Present and are a critical precondition for submarine slope failure.     

Unique achondrite Northwest Africa 11042: Exploring the melting and breakup of the L chondrite parent body

Northwest Africa (NWA) 11042 is a heavily shocked achondrite with medium‐grained cumulate textures. Its olivine and pyroxene compositions, oxygen isotopic composition, and chromium isotopic composition are consistent with L chondrites. Sm‐Nd dating of its primary phases shows a crystallization age of 4100 ± 160 Ma. Ar‐Ar dating of its shocked mineral maskelynite reveals an age of 484.0 ± 1.5 Ma. This age coincides roughly with the breakup event of the L chondrite parent body evident in the shock ages of many L chondrites and the terrestrial record of fossil L chondritic chromite. NWA 11042 shows large depletions in siderophile elements (<0.01×CI) suggestive of a complex igneous history involving extraction of a Fe‐Ni‐S liquid on the L chondrite parent body. Due to its relatively young crystallization age, the heat source for such an igneous process is most likely impact. Because its mineralogy, petrology, and O isotopes are similar to the ungrouped achondrite NWA 4284 (this work), the two meteorites are likely paired and derived from the same parent body.     

Experimental impact cratering: A summary of the major results of the MEMIN research unit

This paper reviews major findings of the Multidisciplinary Experimental and Modeling Impact Crater Research Network (MEMIN). MEMIN is a consortium, funded from 2009 till 2017 by the German Research Foundation, and is aimed at investigating impact cratering processes by experimental and modeling approaches. The vision of this network has been to comprehensively quantify impact processes by conducting a strictly controlled experimental campaign at the laboratory scale, together with a multidisciplinary analytical approach. Central to MEMIN has been the use of powerful two-stage light-gas accelerators capable of producing impact craters in the decimeter size range in solid rocks that allowed detailed spatial analyses of petrophysical, structural, and geochemical changes in target rocks and ejecta. In addition, explosive setups, membrane-driven diamond anvil cells, as well as laser irradiation and split Hopkinson pressure bar technologies have been used to study the response of minerals and rocks to shock and dynamic loading as well as high-temperature conditions. We used Seeberger sandstone, Taunus quartzite, Carrara marble, and Weibern tuff as major target rock types. In concert with the experiments we conducted mesoscale numerical simulations of shock wave propagation in heterogeneous rocks resolving the complex response of grains and pores to compressive, shear, and tensile loading and macroscale modeling of crater formation and fracturing. Major results comprise (1) projectile–target interaction, (2) various aspects of shock metamorphism with special focus on low shock pressures and effects of target porosity and water saturation, (3) crater morphologies and cratering efficiencies in various nonporous and porous lithologies, (4) in situ target damage, (5) ejecta dynamics, and (6) geophysical survey of experimental craters.     

High‐pressure phase transitions of α‐quartz under nonhydrostatic dynamic conditions: A reconnaissance study at PETRA III

Hypervelocity collisions of solid bodies occur frequently in the solar system and affect rocks by shock waves and dynamic loading. A range of shock metamorphic effects and high‐pressure polymorphs in rock‐forming minerals are known from meteorites and terrestrial impact craters. Here, we investigate the formation of high‐pressure polymorphs of α‐quartz under dynamic and nonhydrostatic conditions and compare these disequilibrium states with those predicted by phase diagrams derived from static experiments under equilibrium conditions. We create highly dynamic conditions utilizing a mDAC and study the phase transformations in α‐quartz in situ by synchrotron powder X‐ray diffraction. Phase transitions of α‐quartz are studied at pressures up to 66.1 and different loading rates. At compression rates between 0.14 and 1.96 GPa s^?1, experiments reveal that α‐quartz is amorphized and partially converted to stishovite between 20.7 GPa and 28.0 GPa. Therefore, coesite is not formed as would be expected from equilibrium conditions. With the increasing compression rate, a slight increase in the transition pressure occurs. The experiments show that dynamic compression causes an instantaneous formation of structures consisting only of SiO₆ octahedra rather than the rearrangement of the SiO₄ tetrahedra to form a coesite. Although shock compression rates are orders of magnitude faster, a similar mechanism could operate in impact events.     

Combining shock barometry with numerical modeling: Insights into complex crater formation—The example of the Siljan impact structure (Sweden)

Siljan, central Sweden, is the largest known impact structure in Europe. It was formed at about 380 Ma, in the late Devonian period. The structure has been heavily eroded to a level originally located underneath the crater floor, and to date, important questions about the original size and morphology of Siljan remain unanswered. Here we present the results of a shock barometry study of quartz‐bearing surface and drill core samples combined with numerical modeling using iSALE. The investigated 13 bedrock granitoid samples show that the recorded shock pressure decreases with increasing depth from 15 to 20 GPa near the (present) surface, to 10–15 GPa at 600 m depth. A best‐fit model that is consistent with observational constraints relating to the present size of the structure, the location of the downfaulted sediments, and the observed surface and vertical shock barometry profiles is presented. The best‐fit model results in a final crater (rim‐to‐rim) diameter of ~65 km. According to our simulations, the original Siljan impact structure would have been a peak‐ring crater. Siljan was formed in a mixed target of Paleozoic sedimentary rocks overlaying crystalline basement. Our modeling suggests that, at the time of impact, the sedimentary sequence was approximately 3 km thick. Since then, there has been around 4 km of erosion of the structure.     

Preprocessing of gravity gradients at the GOCE high-level processing facility

One of the products derived from the gravity field and steady-state ocean circulation explorer (GOCE) observations are the gravity gradients. These gravity gradients are provided in the gradiometer reference frame (GRF) and are calibrated in-flight using satellite shaking and star sensor data. To use these gravity gradients for application in Earth scienes and gravity field analysis, additional preprocessing needs to be done, including corrections for temporal gravity field signals to isolate the static gravity field part, screening for outliers, calibration by comparison with existing external gravity field information and error assessment. The temporal gravity gradient corrections consist of tidal and nontidal corrections. These are all generally below the gravity gradient error level, which is predicted to show a 1/f behaviour for low frequencies. In the outlier detection, the 1/f error is compensated for by subtracting a local median from the data, while the data error is assessed using the median absolute deviation. The local median acts as a high-pass filter and it is robust as is the median absolute deviation. Three different methods have been implemented for the calibration of the gravity gradients. All three methods use a high-pass filter to compensate for the 1/f gravity gradient error. The baseline method uses state-of-the-art global gravity field models and the most accurate results are obtained if star sensor misalignments are estimated along with the calibration parameters. A second calibration method uses GOCE GPS data to estimate a low-degree gravity field model as well as gravity gradient scale factors. Both methods allow to estimate gravity gradient scale factors down to the 10⁻³ level. The third calibration method uses high accurate terrestrial gravity data in selected regions to validate the gravity gradient scale factors, focussing on the measurement band. Gravity gradient scale factors may be estimated down to the 10⁻² level with this method.