Reconstructing fluvial bar surfaces from compound cross‐strata and the interpretation of bar accretion direction in large river deposits

The interpretation of fluvial styles from the rock record is based for a significant part on the identification of different types of fluvial bars, characterized by the geometric relationship between structures indicative of palaeocurrent and surfaces interpreted as indicative of bar form and bar accretion direction. These surfaces of bar accretion are the boundaries of flood‐related bar increment elements, which are typically less abundant in outcrops than what would be desirable, particularly in large river deposits in which each flood mobilizes large volumes of sediment, causing flood‐increment boundary surfaces to be widely spaced. Cross‐strata set boundaries, on the other hand, are abundant and indirectly reflect the process of unit bar accretion, inclined due to the combined effect of the unit bar surface inclination and the individual bedform climbing angle, in turn controlled by changes in flow structure caused by local bar‐scale morphology. This work presents a new method to deduce the geometry of unit bar surfaces from measured pairs of cross‐strata and cross‐strata set boundaries. The method can be used in the absence of abundant flood‐increment bounding surfaces; the study of real cases shows that, for both downstream and laterally accreting bars, the reconstructed planes are very similar to measured bar increment surfaces.     

Assessment of small-diameter shallow wells for managed aquifer recharge at a site in southern Styria,Austria

An approach to establish the recharge component of managed aquifer recharge (MAR) has recently been proposed that uses small-diameter shallow wells installed using relatively inexpensive drilling methods such as direct push. As part of further development of that approach, a generalized procedure is presented for a technical and economic assessment of the approach’s potential in comparison to other systems. Following this procedure, the use of small-diameter wells was evaluated both experimentally and numerically for a site located in southern Styria, Austria. MAR is currently done at the site using a horizontal pipe infiltration system, and system expansion has been proposed with a target rate of 12 l/s using small-diameter wells as one possible option. A short-duration single-well field recharge experiment (recharge rate 1.3–3.5 l/s) was performed (recharge by gravity only). Numerical modeling of the injection test was used to estimate hydraulic conductivity (K). Quasi-steady-state, single-well recharge simulations for different locations, as well as a long-term transient simulation, were performed using the K value calibrated from the field injection test. Results indicate that a recharge capacity of 4.1 l/s was achievable with a maximum head rise of 0.2 m at the injection well. Finally, simulations were performed for three different well fields (4, 6 and 8 wells, respectively) designed to infiltrate a target rate of 12 l/s. The experimental and numerical assessments, supported by a cost analysis of the small-diameter wells, indicate that the small-diameter wells are a viable, cost-effective recharge approach at this and other similar sites.     

Indirect unstructured hex-dominant mesh generation using tetrahedra recombination

Corner-point gridding is widely used in reservoir and basin modeling but generally yields approximations in the representation of geological interfaces. This paper introduces an indirect method to generate a hex-dominant mesh conformal to 3D geological surfaces and well paths suitable for finite-element and control-volume finite-element simulations. By indirect, we mean that the method first generates an unstructured tetrahedral mesh whose tetrahedra are then merged into primitives (hexahedra, prisms, and pyramids). More specifically, we focus on determining the optimal set of primitives that can be recombined from a given tetrahedral mesh. First, we detect in the tetrahedral mesh all the feasible volumetric primitives using a pattern-matching algorithm (Meshkat and Talmor Int. J. Numer. Meth. Eng. 49(1-2), 17–30 2000) that we re-visit and extend with configurations that account for degenerated tetrahedra (slivers). Then, we observe that selecting the optimal set of primitives among the feasible ones can be formalized as a maximum weighted independent set problem (Bomze et al. 1999), known to be \(\mathcal {N}\mathcal {P}\)-Complete. We propose several heuristic optimizations to find a reasonable set of primitives in a practical time. All the tetrahedra of each selected primitive are then merged to build the final unstructured hex-dominant mesh. This method is demonstrated on 3D geological models including a faulted and folded model and a discrete fracture network.