Thermal and seismic hints for chimney type cross-stratal fluid flow in onshore basins.

When modelling onshore sedimentary basins, modellers generally assume that semi-permeable layers (aquitards) greatly restrict vertical flow between aquifers. Aquitards are therefore considered as confining media and vertical flow is assumed to take place mainly within localised permeable faults, if any. In the offshore context, however, interpretation of seismic data frequently provides evidence of fluid flow between sedimentary layers via structurally disrupted formations (pervasive fractures) recognised as zones of reduced seismic amplitude and generically called "chimneys". Here we show that chimneys are also present onshore, and that they crosscut confining layers. In the Anglo-Paris Basin, seismic data suggest 1 to 2 km wide zones of disrupted seismic signal spatially correlated to a hitherto unexplained major temperature anomaly of 20 °C. When included in geothermal models using a five-order increase in permeabilities with respect to confining layers, we find that fluid flows vertically through aquifers and confining layers, thereby explaining this major temperature anomaly. Despite the importance of their hydrodynamic and thermal impacts, chimneys - less obvious than faults - have been overlooked as fluid flow paths in many onshore sedimentary basins exploited for their resources. This indicates a clear need for better understanding of pervasive flow paths, especially as the resources and properties of basins (i.e. conventional and unconventional hydrocarbons, geothermal potential, CO2 storage, nuclear waste repository, drinking water, etc.) are increasingly being harnessed.

Generation of 3D Spatially Variable Anisotropy for Groundwater Flow Simulations.

Sedimentary units generally present anisotropy in their hydraulic properties, with higher hydraulic conductivity along bedding planes, rather than perpendicular to them. This common property leads to a modeling challenge if the sedimentary structure is folded. In this paper, we show that the gradient of the geological potential used by implicit geological modeling techniques can be used to compute full hydraulic conductivity tensors varying in space according to the geological orientation. For that purpose, the gradient of the potential, a vector normal to the bedding, is used to construct a rotation matrix that allows the estimation of the 3D hydraulic conductivity tensor in a single matrix operation. A synthetic 2D cross section example is used to illustrate the method and show that flow simulations performed in such a folded environment are highly influenced by this rotating anisotropy. When using the proposed method, the streamlines follow very closely the folded formation. This is not the case with an isotropic model.