Claire Bossennec
email: claire.bossennec@univ-lorraine.fr

Yves Géraud
email: yves.geraud@univ-lorraine.fr

Isabelle Moretti
email: isabelle.moretti@engie.com

Luca Mattioni
email: luca.mattioni@engie.com

Didier Stemmelen
email: didier.stemmelen@univ-lorraine.fr

Bossennec Claire

Mattioni Luca

Stemmelen Didier

Georessources

Pore network properties of sandstones in a fault damage zone

Keywords: Damage zone, normal fault, pore network, fracture pattern, Buntsandstein sandstones

niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

The understanding of fluid-flow characteristics and their controlling factors has a large spectrum of applications. These include natural risk evaluation, an understanding of mechanics of faulting, resource prospection and production (ore, hydrocarbons, geothermal resources). Much recent research has focused on understanding the characteristics and parameters (geometrical, structural, chemical) of pore, fracture and fault networks, along with their impact on flow properties [START_REF] Caine | Fault zone architecture and permeability structure[END_REF][START_REF] Caine | Internal structure, fault rocks, and inferences regarding deformation, fluid flow, and mineralization in the seismogenic Stillwater normal fault, Dixie Valley, Nevada[END_REF][START_REF] Wibberley | Recent advances in the understanding of fault zone internal structure: A review[END_REF][START_REF] Faulkner | A review of recent developments concerning the structure, mechanics and fluid flow properties of fault zones[END_REF][START_REF] Loveless | Fault architecture and deformation processes within poorly lithi fi ed rift sediments , Central Greece[END_REF][START_REF] Bense | Fault zone hydrogeology[END_REF].

A fault zone is defined as a lithological heterogeneous and structurally anisotropic discontinuity in the rock. This is caused by deformation and displacement of blocks of rock, delimiting the footwall and the hanging wall [START_REF] Caine | Fault zone architecture and permeability structure[END_REF]. In a fault zone, deformation of the rock results in different size objects from pluri-metric fractures or minor faults, to micrometric cracks. The term "fracture" or "fault" will be reserved here for metric to pluri-centimetric discontinued planes. In contrast, discontinuities of less than a centimeter will be referred as "cracks".

The overall global structure of a fault zone includes a fault core surrounded by damage zone domains [START_REF] Chester | Implications for Mechanical Properties of Brittle Faults from Observations of the Punchbowl Fault Zone, California[END_REF][START_REF] Antonellini | Effect of faulting on fluid flow in porous sandstones:petrophysical properties[END_REF][START_REF] Caine | Fault zone architecture and permeability structure[END_REF][START_REF] Caine | Internal structure, fault rocks, and inferences regarding deformation, fluid flow, and mineralization in the seismogenic Stillwater normal fault, Dixie Valley, Nevada[END_REF][START_REF] Heynekamp | Controls on Fault-Zone Architecture in Poorly Lithified Sediments , Rio Grande Rift , New Mexico : Implications for Fault-Zone Permeability and Fluid Flow[END_REF][START_REF] Jourde | Computing permeability of fault zones in eolian sandstone from outcrop measurements[END_REF][START_REF] Fossen | Fault interaction in porous sandstone and implications for reservoir management; examples from southern Utah[END_REF][START_REF] Faulkner | A review of recent developments concerning the structure, mechanics and fluid flow properties of fault zones[END_REF][START_REF] Bense | Fault zone hydrogeology[END_REF]. The damage zone has characteristics that vary according to the rock type: 1) For brittle consolidated rocks, in the model of [START_REF] Caine | Fault zone architecture and permeability structure[END_REF] the damage zone is defined as an increased fracture density and contains minor faults compared to the protolith. 2) For porous sandstones, models of the damage zone from [START_REF] Heynekamp | Controls on Fault-Zone Architecture in Poorly Lithified Sediments , Rio Grande Rift , New Mexico : Implications for Fault-Zone Permeability and Fluid Flow[END_REF] modified by [START_REF] Loveless | Fault architecture and deformation processes within poorly lithi fi ed rift sediments , Central Greece[END_REF], include minor faults and deformation bands [START_REF] Antonellini | Effect of faulting on fluid flow in porous sandstones:petrophysical properties[END_REF][START_REF] Flodin | Representation of Fault Zone Permeability in Reservoir Flow Models[END_REF][START_REF] Labaume | Diagenesis-dependence of cataclastic thrust fault zone sealing in sandstones. Example from the Bolivian Sub-Andean Zone[END_REF][START_REF] Jourde | Computing permeability of fault zones in eolian sandstone from outcrop measurements[END_REF][START_REF] Odling | Permeability scaling properties of fault damage zones in siliclastic rocks[END_REF][START_REF] Fossen | Deformation bands and their influence on fluid flow[END_REF][START_REF] Fossen | Deformation bands formed during soft-sediment deformation: Observations from SE Utah[END_REF], but fractures are mostly absent.

In porous rocks, fluid flow properties are strongly influenced by the pore network geometry and its evolution during diagenesis. Significant modifications occur if these are affected by fault zones. Diagenetic processes are all post-depositional ones (e.g. physical, chemical and biological reactions to reach a physico-chemical equilibrium between grains and pore fluids). These lead to modification of sediments [START_REF] Worden | Sandstone Diagenesis: The Evolution of Sand to Stone, Sandstone Diagenesis: Recent and Ancient[END_REF]. The processes are ubiquitous in sandstones. They are active during burial and the possible uplift cycle of the basin; including compaction, temperature increase and lithification. They occur in both un-deformed protolith (interstitial pore fluids) and deformation zones (fractures, deformation bands, or at a larger scale on a fault core), and depend on the chemistry of the fluids, mineralogy and morphology of grains and temperature. The diagenetic processes that affect sandstones are detailed in [START_REF] Worden | Sandstone Diagenesis: The Evolution of Sand to Stone, Sandstone Diagenesis: Recent and Ancient[END_REF], and are merely summarized here as a reminder: 1) Authigenesis through quartz and feldspar precipitation in the pore space. Sources for silica can be assorted, such as quartz intergranular dissolution, grain chemical compaction, feldspar dissolution, smectite or kaolinite illitization [START_REF] Worden | Quartz Cementation in Oil Field Sandstones: A Review of the Key Controversies[END_REF][START_REF] Worden | Sandstone Diagenesis: The Evolution of Sand to Stone, Sandstone Diagenesis: Recent and Ancient[END_REF][START_REF] Tournier | Relationship between Deep Diagenetic Quartz Cementation and Sedimentary Facies in a Late Ordovician Glacial Environment (Sbaa Basin, Algeria)[END_REF]. 2) Mechanical compaction. 3) Dissolutions from chemical compaction and fluid leaching. 4) Neomorphysm, recrystallization and replacement of clays and feldspars: sheets of illite partially fill and compartmentalize the void spaces, especially those of micro-and meso-pores. Illite is either detrital clay mineral, or is also formed during advanced diagenesis (burial and fluid-rock interaction). The latter is issued from the transformation of smectite into illite, through intermediate mixed-layer illite/smectite; and/or is a product of K-feldspar albitization with the presence of kaolinite [START_REF] Gaupp | Diagenesis and fluid evolution of deeply buried Permian (Rotliegende) gas reservoirs, northwest Germany[END_REF][START_REF] Mckinley | Smectite in Sandstones: A Review of the Controls on Occurrence and Behaviour During Diagenesis[END_REF][START_REF] Worden | Clay Mineral Cements in Sandstones[END_REF][START_REF] Fischer | Products and timing of diagenetic processes in Upper Rotliegend sandstones from Bebertal (North German Basin, Parchim Formation, Flechtingen High, Germany)[END_REF].

These processes interact with one another, because of their potential co-eval formation [START_REF] Clauer | Thermal history of Triassic sandstones from the Vosges Mountains-Rhine Graben rifting area, NE France, based on K-Ar illite dating[END_REF][START_REF] Worden | Clay Mineral Cements in Sandstones[END_REF].They have an impact on the pore network and reservoir quality of sandstones, though their quantification is rarely determined [START_REF] Bjorlykke | Quartz cementation in sedimentary basins[END_REF][START_REF] Bjorkum | How Important is Pressure in causing Dissolution of Quartz in Sandstones?[END_REF][START_REF] Worden | Sandstone Diagenesis: The Evolution of Sand to Stone, Sandstone Diagenesis: Recent and Ancient[END_REF][START_REF] Eichhubl | Structural and diagenetic control of fluid migration and cementation along the Moab fault , Utah[END_REF][START_REF] Henares | Diagenetic constraints on heterogeneous reservoir quality assessment : A Triassic outcrop analog of meandering fl uvial reservoirs[END_REF].

Transfer properties within the fault zone depend on the properties of the fault core and of the damage zones. In the latter, the permeability parallel to the fault can be enhanced using the fracture network (Fossen et al., 2007;[START_REF] Bense | Fault zone hydrogeology[END_REF].

Hydraulic behavior of the damage zone depends on distinct factors, which are scalerelated:

1) Concerning the fracture network scale; namely its geometry, length and aperture distribution of fractures, orientation, density, connectivity, and infilling of micro-structures. The drain effect of the fracture network in the damage zone has to be nuanced by taking into account the local heterogeneities enhanced by diagenetic processes specific to fault zones. Faults and fractures can be cemented (partly or completely), and lead to a compartmentalization of the reservoir [START_REF] Jolley | Reservoir compartmentalization: an introduction[END_REF].

2) At the scale of the pore network, faulting results in changes of grain shape and size (physical and chemical compaction, cataclasis) [START_REF] Schueller | Spatial distribution of deformation bands in damage zones of extensional faults in porous sandstones: Statistical analysis of field data[END_REF][START_REF] Farrell | Anisotropic pore fabrics in faulted porous sandstones[END_REF], and increased cementation (a higher surface of quartz is available to cement nucleation). In sandstones, where intergranular spaces represent the large majority of the total porosity, cataclasis and compaction, combined with cementation, can change the amount of porosity, and its pore size distribution. This directly impacts on the permeability of the media.

The fault models presented above are not applicable to all rock types, and do not take into consideration the multi-scale impact of both deformation and diagenetic processes. For example, the model for porous sandstones [START_REF] Loveless | Fault architecture and deformation processes within poorly lithi fi ed rift sediments , Central Greece[END_REF] indicates deformation bands in the damage zone. These are not observed in our case study. We propose here another type of fault damage zone model for lithified porous sandstones in a normal faulting context, for the faulted sandstones of Cleebourg (GPS coordinates: 32U, 418833 E, 5429554 N) (Fig. 1).

In this paper, we integrate the characterization of sedimentary and diagenetic heterogeneities, and the spatial organization of the structural features. Thus we can estimate their control on hydraulic and geo-mechanical properties. Our results suggest that transfer properties of the fault zone depend on macrostructure organization (fracture corridors, micro-faults). This organization controls the microstructure of rock (sedimentary and diagenetic characteristics), and associated grain scale pore network properties. Here, we show that the distance to the fault-core, and the fracture network characterization at the fault scale, are not sufficient to describe and understand transfer and geo-mechanical property distributions. The characterization of the evolution of the pore network properties requires investigation at both intermediate scale of a single fracture and that of grain scale (variable diagenesis of the pore network). We also demonstrate that this multiple embedded boxes approach is useful to understand the distribution pattern of hydraulic and geo-mechanical properties within the damage zone, and its implications in terms of fault zone behavior.

After defining the geological context of the case study, the results obtained by applied methodology will be presented. These will be specifically on aspects of petrography, porous network properties (combined X-Ray, MIP and NMR quantification), crack pattern, P-waves, and Young's Modulus properties. These results are discussed in terms of multi-parameter integration. Thus the effects of properties and damage-zone compartmentalization, of sedimentary and diagenetic heritage, of fractures and modification of the porous network through cracks, and fluid/rock interactions on fluid flow can be deciphered.

Geological settings

2.1.1.Regional geology

The Upper Rhine Graben (URG), with a 300 km length and 40-50 km width in a NNE-SSW to N-S trend (Fig. 1), is a typical example of syn-orogenic Cenozoic continental foreland rifting [START_REF] Ziegler | European Cenozoic rift system[END_REF][START_REF] Merle | The formation of the West European rift: A new model as exemplified by the Massif Central area[END_REF][START_REF] Bourgeois | Separation of rifting and lithospheric folding signatures in the NW-Alpine foreland[END_REF][START_REF] Edel | The Variscan tectonic inheritance of the Upper Rhine Graben: evidence of reactivations in the Lias, Late Eocene-Oligocene up to the recent[END_REF][START_REF] Ford | Kinematic and geometric analysis of fault-related folds in a rift setting: The Dannemarie basin, Upper Rhine Graben, France[END_REF][START_REF] Ziegler | Cenozoic uplift of Variscan Massifs in the Alpine foreland: Timing and controlling mechanisms[END_REF]. The URG is limited by the Vosges and Black Forest massifs, respectively constituting the western and eastern shoulders. To the north, the URG is bounded by the Rhenish Massif, and to the south by the Jura Massif, a folded alpine structure. This basin is segmented into several subbasins, which are delimited toward the URG axis by oblique faults. The stratigraphic column of the URG area is presented in Figure 2, with associated insight into the regional tectonic calendar.

The group studied is the Buntsandstein Gp., also known locally as "Grès Vosgien Fm", which dates from the Early Triassic (Olenekian) period. Due to its high porosity, this group constitutes often a good reservoir for both geothermal and hydrocarbon resources [START_REF] Haffen | Caractéristiques géothermiques du réservoir gréseux du Buntsandstein d'Alsace[END_REF][START_REF] Haffen | Geothermal, structural and petrophysical characteristics of Buntsandstein sandstone reservoir in the Upper Rhine Graben[END_REF][START_REF] Böcker | An overview on source rocks and the petroleum system of the central Upper Rhine Graben[END_REF]. These Early Triassic sandstones were deposited in fluvial and playa-lake systems, at the southern border of the Western Germany sag basin [START_REF] Bourquin | Lower Triassic sequence stratigraphy of the western part of the Germanic Basin (west of Black Forest): Fluvial system evolution through time and space[END_REF][START_REF] Bourquin | Braided rivers within an arid alluvial plain (example from the lower triassic, western german basin): Recognition criteria and expression of stratigraphic cycles[END_REF][START_REF] Soyk | Diagenesis and reservoir quality of the Lower and Middle Buntsandstein (Lower Triassic)[END_REF].

The Buntsandstein Gp. forms the bottom of the sedimentary pile topping the basement. These sandstones have undergone a complex burial history and associated diagenetic events [START_REF] Clauer | Thermal history of Triassic sandstones from the Vosges Mountains-Rhine Graben rifting area, NE France, based on K-Ar illite dating[END_REF][START_REF] Dresmann | The south-western Black Forest and the Upper Rhine Graben main border fault: Thermal history and hydrothermal fluid flow[END_REF][START_REF] Blaise | Reconstructing fluid-flow events in Lower-Triassic sandstones of the eastern Paris Basin by elemental tracing and isotopic dating of nanometric illite crystals[END_REF]. The area of the future URG was in subsidence from the Triassic to Jurassic eras [START_REF] Ziegler | Collision related intra-plate compression deformations in Western and Central Europe[END_REF][START_REF] Ziegler | European Cenozoic rift system[END_REF]. This phase was followed by an uplift-erosion phase from the Cretaceous to Early Eocene periods. At the Late Eocene period, a first phase of left-lateral shears started, along the Paleozoic inherited fault system [START_REF] Ziegler | Post-Variscan evolution of the lithosphere in the Rhine Graben area: constraints from subsidence modelling[END_REF][START_REF] Edel | The Variscan tectonic inheritance of the Upper Rhine Graben: evidence of reactivations in the Lias, Late Eocene-Oligocene up to the recent[END_REF].

Several sub-basins developed and nucleated, with a NNE-SSW orientation. These were sub-parallel to σ1, the major horizontal stress. During Oligocene, the major rift phase took place. During Early Oligocene, there was an activation (or even reactivation) of regional fault systems, in a normal displacement [START_REF] Schumacher | Upper Rhine Graben: Role of preexisting structures during rift evolution[END_REF].

The minimum principal stress σ3 is estimated to be horizontal and oriented WNW-ESE [START_REF] Schumacher | Upper Rhine Graben: Role of preexisting structures during rift evolution[END_REF][START_REF] Derer | The northern Upper Rhine Graben: basin geometry and early syn-rift tectono-sedimentary evolution[END_REF], with subsidence in the URG Basin, and uplift and erosion on the Graben shoulders [START_REF] Timar-Geng | Low-temperature thermochronology of the flanks of the southern Upper Rhine Graben[END_REF][START_REF] Dresmann | The south-western Black Forest and the Upper Rhine Graben main border fault: Thermal history and hydrothermal fluid flow[END_REF][START_REF] Brockamp | Mesozoic hydrothermal impact on Rotliegende and Bunter immature sandstones of the High Rhine trough and its adjacent eastern area (southern Black Forest, Germany)[END_REF]. During Late Oligocene, the regional stress field evolved towards a σ3 oriented NE-SW, accompanied by sinistral shearing displacement along the major faults oriented NS. From Miocene to current times, the global stress field returned to its initial configuration, accompanied by left-lateral shear, and the development of small pull-apart basins. It maintained the uplift of shoulders of the Vosges -Black Forest Massif areas [START_REF] Schumacher | Upper Rhine Graben: Role of preexisting structures during rift evolution[END_REF][START_REF] Ziegler | Post-Variscan evolution of the lithosphere in the Rhine Graben area: constraints from subsidence modelling[END_REF][START_REF] Timar-Geng | Low-temperature thermochronology of the flanks of the southern Upper Rhine Graben[END_REF]. These listed tectonic events initiated and/or reactivated fault zones and fracture networks, and involved superficial and deep fluid flows [START_REF] Lampe | Episodic fluid flow within continental rift basins: Some insights from field data and mathematical models of the Rhine Graben[END_REF][START_REF] Stober | Hydraulic and hydrochemical properties of deep sedimentary reservoirs of the Upper Rhine Graben[END_REF].

'Insert Figure 1 here' 'Insert Figure 2 here'

2.1.2.Outcrop localization and structural context

On the western border of the URG, the Saverne Fracture Field (SFF) is delimited by the Vosgian Fault and Rhenish Fault [START_REF] Eisbacher | Karlsruhe und seine Region -Nordschwarzwald, Kraichgau, Neckartal, südlicher Odenwald, Oberrhein-Graben, Pfälzerwald und westliche Schwäbische Alb[END_REF]. These two major faults delimit normal conjugate to left-lateral shear fault systems, forming a set of small horsts and grabens, tilted eastward. This configuration allows the outcropping of deeply buried reservoir analogs, from Permian up to Jurassic deposits.

On the N-E extremity of the SFF, the Lower Triassic formations outcrop on the western part of the Hochwald Horst (Fig. 1.B, C). The outcrop studied is the quarry of Cleebourg that is exploited for ornamental stones. In terms of stratigraphy, the "Grès Vosgien Fm" outcrop of Cleebourg belongs to the Lower Buntsandstein Formation [START_REF] Ménillet | Notice explicative, Carte géol[END_REF][START_REF] Bourquin | Lower Triassic sequence stratigraphy of the western part of the Germanic Basin (west of Black Forest): Fluvial system evolution through time and space[END_REF] (Fig. 2). These sandstones are described as fine to medium grained sandstones, with a relatively good sorting. Selected strata for sampling show a fully pervasive bleaching, described as "Haardtrand Fazies" [START_REF] Soyk | Diagenesis and reservoir quality of the Lower and Middle Buntsandstein (Lower Triassic)[END_REF], and belong to playa-lake sand sheet deposits.

The northern face of the quarry consists of a large fault zone, well described by [START_REF] Bauer | Architecture, fracture system, mechanical properties and permeability structure of a fault zone in Lower Triassic sandstone, Upper Rhine Graben[END_REF], with a 14 meters thick inner fault core, surrounded by a 5 meters thick outer fault core. This outer fault core zone is surrounded by a 40 meters thick damage zone, from which the external limit is not outcropping. The inner core zone is constituted of multiple slip surfaces, deformation bands and fractures, separating different units with different lithologies (clays, silts, and sands) and in various consolidation states. The outer core zone is characterized by high densities of fractures separating less deformed host-rock lenses. The border between the fault core and the damage zone is clearly defined by a slip surface. Despite the high porosity of the sandstone, the deformation in the damage zone is mostly brittle, with fractures distributed in corridors, with higher fracture intensity, and some associated minor faults. Some fractures present slight to strong brown iron hydroxide staining.

Sampling and analytical procedure

Damage zone features and sampling

The word matrix as used here describes the porous media in the rock from which fracture porosity is excluded (rock matrix) [START_REF] Warren | The Behavior of Naturally Fractured Reservoirs[END_REF]. To decipher the impact of fractures on matrix properties (in terms of diagenetic heterogeneities and hydraulic behavior), the area selected for sampling was restricted to two continuous sedimentary layers in the damage zone of the fault zone. The core zone (comprising the inner and outer fault core) was not considered, due to weak outcropping conditions and poor lithification.

In the damage zone, 22 samples were collected as oriented plugs and blocks. Aiming to link diagenetic process variations with hydraulic and mechanical rock property changes, we performed the following measurements and observations: petrography, petro-physical, and elastic rock properties determination. Uncertainty remained on initial properties of the protolith in the area, as the outer border of the damage zone was not outcropping.

Petrography for diagenetic sequence determination

Standard 30 µm-thick polished thin sections were prepared for petro-structural analysis. Special focus has been directed towards quartz cementation quantification, through cathode-luminescence observations performed on a cold cathode cathodyn OPEA with a 10-12kV, intensity 6, and 250-300 µa. Cathodoluminescence allows generations of quartz cements, grain deformation, and recrystallization to be determined easily. At the pore scale, the deformation of grains was described in terms of microstructure. Combined secondary electron microscopy (SEM) analysis of micro textural habits and X-Ray Diffraction results were used to decipher clay types and proportions. Attention was paid to petrography, to determine spatial heterogeneities in terms of diagenetic processes, both in matrices and fractures. Mineralogy is expressed in percentage of total surface area on thin section, and determined by computer assisted point-counting (minimum 300 per picture, 4 pictures per sample).

Fracture network characterization methods:

Fracture network characterization at the outcrop scale is based on the results collected on cross-section H1 by [START_REF] Bauer | Architecture, fracture system, mechanical properties and permeability structure of a fault zone in Lower Triassic sandstone, Upper Rhine Graben[END_REF] and the methods described therein.

Crack networks at the thin-section scale

Crack intensities were characterized in two ways, adapted from [START_REF] Nelson | 1 -Evaluating Fractured Reservoirs: Introduction[END_REF]:

• Number of cracks per meter (i.e. P10). This was quantified by counting, with a constant scale of observation (magnification = 200), the number of cracks crossing a 2 mm long scanline on a thin-section picture with the same orientation as the outcrop scanline.

• Total length of cracks per surface unit (i.e. P21). This was measured by selecting 4 spots of 1 mm diameter for each thin section, all the cracks being visible at a constant scale of observation (magnification = 200).

Porosity and permeability

3.5.1.Porosity

Porosity was determined by different methods:

Porosity using water and gas saturation was determined using a helium pycnometer and triple mass methods. Water and Helium porosity gave a value of connected porosity.

Mercury Injection Porosimetry (MIP) [START_REF] Ritter | Pressure Porosimeter and Determination of Complete Macropore-Size Distributions. Pressure Porosimeter and Determination of Complete Macropore-Size Distributions[END_REF] was performed to compute the distribution of porosity versus threshold size. The experimental process was divided into three steps: (1) a first drainage with varying pressure from 0.003 MPa to 220 MPa, allowing the measure of pore throat radii distribution to the smallest pore throats of 0.001 µm; (2) imbibition; (3) second drainage. This protocol was applied to 3 to 4 cm 3 small cores from all the samples, using Micromeritics AutoPore IV, with pore throat radii measure accuracy from 0.003 to 500 µm. Equilibrium time settings were 20 s for low pressure cell (0.003 MPa to 0.2 MPa), and 60 s for high pressure cell (0.2 to 220 MPa).

3.5.2.Permeability

Permeability measurements were carried out using two methods: 1) On the outcrop directly with a transportable air-permeameter Tiny Perm (II) 2) in the laboratory, on sampled plugs using a Nitrogen permeameter, with a Hassler cell, at a confining pressure of 1.5 MPa and stabilized atmospheric foot pressure, establishing a steady state gas flow. At least 8 steps of head pressure were measured per sample, and the Klinkenberg correlation was applied to correct measurements from gas-slippage [START_REF] Klinkenberg | The permeability of porous media to liquids and gases[END_REF]. Permeabilities were estimated in millidarcy (mD).

Elastic rock properties and deduced mechanical properties

To estimate mechanical rock properties, propagation velocities of P-and S-waves were measured, on 16 dry samples, with 250 kHz electrodes installed on the PunditLab Proceeq ® tool. Dynamic Young's Modulus, Poisson's ratio and the fissuration ratio were calculated, to estimate variations of rigidity and compressibility of the media, and to estimate the effective amount of fractures in the media. To determine Poisson's ratio (ν), Eq. ( 1), and Young's Modulus (E), Eq.( 2), formulae are the following (Homand etal., 2000):

𝜈 = ! " !" #$ #% # " $!" #$ #% # " (1) 𝐸 = 𝜌 %&²($)*)($!,*) $!* (2)
These dynamic moduli differ from static moduli by their method of establishment. Static moduli are defined according to the results of a mechanical experiment as compression/shear tests, whereas dynamic moduli are based on acoustic properties.

Within the theory of elasticity, both methods should have given identical results.

Though for rock media, which are not perfectly elastic, differences were observed between static and dynamic moduli [START_REF] Homand | Manuel de mécanique des roches -Tome 1 Fondements[END_REF]. These differences were explained by a strain value effect. The strain value (a dimensionless value, which is a function of displacement divided by the wavelength) of both methods differed from 10 - 8 to 10 -6 for dynamic to a minimum of 10 -4 for static methods [START_REF] Zinszner | A geoscientist's guide to petrophysics[END_REF].

The P-wave propagation is perturbed by crack porosity. The measure of the P-wave propagation velocity was used to estimate the rate of open discontinuities in the sample. Fissuration rate (FR) determination was based on the continuity indicator (IC), defined as the ratio of measured P-wave propagation velocity Vp on Vp*. Vp* is the Pwave propagation velocity calculated from the mineralogical content, as the arithmetic weighted mean of Vi, propagation velocity in each mineral phase, and the coefficient of weight Ci, as a proportion of the mineral phase i in the whole rock, Eq.( 3), [START_REF] Homand | Manuel de mécanique des roches -Tome 1 Fondements[END_REF]. For non-cracked rocks, IC, written ICp, was defined by the experimental relation, Eq.( 4), with Φ the porosity being determined by thin section analysis. FR was defined in Eq.( 5), as the relative difference between the ICp (for non-cracked rocks) (Eq.( 4)) and IC calculated from Eq.(3).

𝐼𝐶(%) = 100 %& %& * , with 𝑉 & * = ∑ 𝑐 . 𝑉 &,.
.

(3)

𝐼𝐶𝑝 = 100(1 -1.4𝜙) (4) 𝐹𝑅 = 1 - 01 01&
(5)

Nuclear Magnetic Resonance (NMR) relaxometry

Nuclear Magnetic Resonance (NMR) methods are based on the measurement of the magnetization decay of water molecule hydrogen nuclei and were helpful to estimate the real pore size distribution [START_REF] Dunn | Nuclear Magnetic Resonance: Petrophysical and Logging Applications[END_REF][START_REF] Zinszner | A geoscientist's guide to petrophysics[END_REF][START_REF] Yao | Petrophysical characterization of coals by low-field nuclear magnetic resonance[END_REF][START_REF] Vincent | NMR relaxation of neritic carbonates: An integrated petrophysical and petrographical approach[END_REF][START_REF] Faÿ-Gomord | New insight into the microtexture of chalks from NMR analysis[END_REF].

In this study, T2, e.g. the transverse relaxation time along a plan perpendicular to B0, was used, and measured with a Carr-Purcell-Meiboom-Gill sequence (CPMG), at regular time intervals 2τ (or TE) of 150 µs. The diffusion processes were supposed to be fast, and the inter-echo times (TE) used in experiments were very low. Thus, the term of diffusivity can be neglected in Eq. ( 6) of T2 concerning V/S, (respectively the compartment of porous media volume and surface):

$ 2 " = $ 2 "'()* + 𝜌 3 4 % + (2567) " 8 $, (6)
where T2 bulk is the bulk relaxation time of free water (2.6 s at 25°C), and ρs the specific surface relaxivity (constant for each sample, determined by calibration with X-ray tomography data, see section 3.9). TE is the inter-echo time, γ the gyromagnetic ratio, G the average local magnetic field gradient, and D the auto-diffusion coefficient of the fluid.

Considering a model of spherical particles, pore radii were approximated from equation ( 7). Combining equations ( 6) and ( 7), the void radius was expressed as a proxy of relaxation time, as written in equation ( 8). (8)

X-ray micro-tomography

The porosity quantification by crossed use of X-ray tomography and MIP is strongly dependent on the scanning resolution [START_REF] Peng | Using X-ray computed tomography in pore structure characterization for a Berea sandstone: Resolution effect[END_REF]. X-ray tomography imaging was thus performed on a Phoenix Nanotom S ®, at the highest resolution, i.e. with a magnification of 250. This allowed a resolution of 1 voxel = 1 µm 3 . Five representative samples of the different facies, with small chips of 3*5*3 mm, were analyzed to calibrate the meso-pore network characteristics (connectivity, percentage of porosity, tortuosity, and equivalent pore radius).

The 3D image of each sample was treated with VG Studio software (Volume Graphics), to extract the porosity and estimate volume, surface, and sphericity. Equivalent radii could be determined from the volume and sphericity for each connected volume, assuming a pore to be a deformed sphere. These connected volumes are added from the largest to the smallest, to build an equivalent of an X-ray cumulated porosity curve.

Combination of pore access size (MIP) and pore volume quantifications (NMR and X-ray tomography)

A major challenge was the building of a porous network model that integrated the throat distribution obtained with MIP, and the void size determined from NMR measurements and X-ray tomography images. Previous studies presented computed data on pore networks of X-ray, tomography, SEM [START_REF] Bera | Understanding the micro structure of Berea Sandstone by the simultaneous use of micro-computed tomography (micro-CT) and focused ion beam-scanning electron microscopy (FIB-SEM)[END_REF], and energy dispersive X-ray spectroscopy (EDS) [START_REF] De Boever | Datafusion of high resolution X-ray CT, SEM and EDS for 3D and pseudo-3D chemical and structural characterization of sandstone[END_REF]. The methodology we proposed combined MIP, NMR, SEM and X-Ray tomography data on pore fabric, and is illustrated in Figure 3.

NMR T2 distribution was converted into an equivalent pore radius distribution and was superposed onto X-ray pore size distributions that fitted large pore sized MIP cumulative porosity curves (>20µm). The largest pores observed by X-ray tomography corresponded to the NMR and MIP first peak. Thus, using this calibration, quantification of smaller pore size by NMR data was then possible. Such pores had not been directly detected by X-ray tomography. A comparison of NMR and MIP results allowed the quantification of several small pore families and pore network parameters (e.g. pore throat size, void size, tortuosity, and surface). Quality control with SEM imagery confirmed the correct size for the void dimension determined by NMR and MIP.

'Insert Figure 3 here'

Results

The total rock mineralogy, diagenetic steps, and pore network morphology were described through thin section petrography, MIP and SEM analysis. Detailed mineralogical compositions were determined by point counting on thin sections (Table 1). First, a description of the sedimentary and diagenetic facies is presented. The structural aspect of the petro-facies will then be described, along with associated pore network, crack patterns and petro-physical properties. As detailed by [START_REF] Bauer | Architecture, fracture system, mechanical properties and permeability structure of a fault zone in Lower Triassic sandstone, Upper Rhine Graben[END_REF], the damage zone outcropping at the north wall of the quarry is affected by a fracture network, in which three corridors are identified (Fig. 4). Fractures are mostly fault synthetic (dip ESE; strike-direction: NNE-SSW), with a secondary minor orthogonal fracture set (strike N-S and E-W). Corridors are defined here when fracture intensity is higher than 4, for two successive one-meter intervals.

Due to the presence of protolith lenses on the scan line, some intervals in the corridors can have lower fracture densities. The corridor ( 1 The first facies Fa1 (Fig. 5, Table 1), is characterized by a medium to fine sub-mature sandstone, with average grain size ranging from 220 to 350 µm, and porosity from 11.43 to 13.4%. The framework is mainly composed of quartz (58.5 to 67 %) and Kfeldspar (8.2 to 18.1 %) sub-angular to sub-rounded grains with well-developed clay coating and rare thin syntaxial cements (0.6 to 2.6%) on both quartz and feldspar grains (average cement thinner than 20 µm). Illite contents vary between 7.6 and 7.9%.

Illite is present as a coating of grains and as pore infill, and has two morpho-types mixed with varying iron hydroxide contents (varying intensities on SEM, and semi quantitative EDS from respectively 4.3 and 16.7 Fe2O3 wt.%). Kaolinite was observed in very low amounts (less than 3%) on feldspar rims and inside cracks in the feldspar grains and cements.

The second facies, Fa2 (Fig. 6, Table 1), is a medium to coarse sandstone, with an average grain size ranging from 350 to 490 µm. Fa2 has thinner illite coatings, and a larger amount of quartz and feldspar diagenetic syntaxic cements (from 5.3 to 9.9 %, with an average thickness for cement between 30 and 60 µm). As for Fa1, the framework is made of sub-rounded grains of quartz (62.3 to 71.4%) and K-feldspar (9.7 to 16.1%). The illite content remains between 1.9 and 4.9%. The same illite morphotypes as described for Fa1 are present in Fa2. The first two morphotypes of illite are relatively similar, consisting mostly of tangled ribbons associated with iron and titanium hydroxides (SEM -XRD) (Fig. 5). A third illite morphotype present in fractured samples has a fibrous habitus.

'Insert Figure 5' 'Insert Figure 6'

4.2.2.Pore network and hydraulic properties

The distributions of pore throats (MIP) and pore size (NMR, tomography) for all facies are both bi-modal, suggesting at least two families of pore voids (P1) and (P2) and two families of throats (S1) and (S2) with the following characteristics (Fig. 7, Pore network properties of fractured Fa1 samples are similar to undeformed facies in terms of bi-modal distribution of pore and pore throat diameters, with the following features (Fig. 7): here the f index indicates fractured facies.

(S1f) has a throat diameter of 0.69 to 1.87 µm; (P1f) a pore diameter of 5 to 12 µm.

This family represents 25 to 77 % of the total porosity. (S2f) has a throat diameter of 9.85 to 12.57 µm; (P2f) a pore diameter of 38 to 50 µm. This family represents 23 to 75 % of the total porosity. Total porosities of fractured and undeformed Fa1 facies are similar, at around 12 %. The permeability for fractured Fa1 sandstone varies from 42 mD up to 235 mD (average 142 mD), and is higher than for low-fractured Fa1 sandstones (with permeability ranging from 28 to 60 mD, average 40 mD).

Similar to Fa1, undeformed Fa2 sample distributions of pore throat and void sizes are bi-modal (Fig. 8, Table 2). However, the ranges of families (P1) and (P2) differ slightly for Fa2 facies, with the following values: (S1) has a throat diameter of 0.45 to 0.66 µm;

(P1) a pore size of 5.9 to 8.

4.2.3.Geo-mechanical properties

Mechanical properties were calculated for each sample from measurements of Vp and Vs (for propagation velocities, see Table 3). The profiles are represented in Figure 9, and values are listed in Table 3. For Fa1 samples, dynamic Young's Moduli E values range from 6.28 to 8.62 GPa, and for Fa2 samples, from 5.36 to 9.29 GPa. With respect to their location within the damage zone, E range from 5.36 GPa in fracture corridors up to 9.29 GPa in low fractured blocks (Fig. 9) and globally decrease toward the fault core. Poisson's ratios are relatively constant (at around 0.3), with no clear relationship between micro-fracture intensity or macro fracture intensity, nor with distance from the fault core. Both facies have fissuration ratio values of between 29 and 47 %. For Fa1the fissuration ratios are on average higher than for Fa2 (respectively av. 42 % and av.

35%)).

'Insert Figure 9' 'Insert Table 3'

Discussion

The impact on rock properties of the presence of fracture corridors and differential diagenetic evolution in the damage zone can be described by a three embedded scales model: 1) at the grain scale, micro-fractures modify pore network properties; 2) at the scale of one fracture and its own volume of influence, dynamic mechanical properties of rocks are affected by micro-fracturing; 3) at the scale of damage zone, overprinting of each corridor occurs. Furthermore, deformation has an impact on rock properties at the scale of the entire damage zone, leading to heterogeneous hydraulic behavior.

Diagenetic para-sequences

Through observations of thin sections (using optical microscopy and SEM), a diagenetic sequence has been built up (Fig. 10) that is consistent with diagenetic sequences previously reported in the literature [START_REF] Haffen | Caractéristiques géothermiques du réservoir gréseux du Buntsandstein d'Alsace[END_REF][START_REF] Blaise | Reconstruction of low burial (< 100 °C) in sedimentary basins: A comparison of geothermometer sensitivity in the intracontinental Paris Basin[END_REF][START_REF] Soyk | Diagenesis and reservoir quality of the Lower and Middle Buntsandstein (Lower Triassic)[END_REF].

During the early stages of burial, a clay-ferruginous grain coating is established, (phase (I)). The nature of the original clay mineral remains unknown. However, according to depositional environment and other mineral phases observed, it might be a kaolinite coating as described by [START_REF] Surma | Microstructures d' un grès affecté par une faille normale: Anisotropie de connectivité et de perméabilité[END_REF], or a smectite/illite coating, which was then fully transformed into Illite ( 1), as proposed in other studies [START_REF] Clauer | Thermal history of Triassic sandstones from the Vosges Mountains-Rhine Graben rifting area, NE France, based on K-Ar illite dating[END_REF][START_REF] Blaise | Reconstructing fluid-flow events in Lower-Triassic sandstones of the eastern Paris Basin by elemental tracing and isotopic dating of nanometric illite crystals[END_REF].

Phase (II) is a cementation phase, occurring on feldspar grains with the formation of feldspar syntaxial cement, followed by quartz cementation (quartz overgrowth around the grains surrounded by thin illite coatings) (Fig. 5,Fig. 6).

Phase (III) consists of the illitization of several mineral phases, and is divided into several steps for each mineralogical contemporaneous reaction with the same reactants and products, and the same temperature required (above 80°C). It is probably preceded by a kaolinitization phase of feldspar grains and cements.

Morphotypes of illite ( 2) and (3) are relatively similar, consisting mostly of tangled ribbons associated with iron and titanium hydroxides (SEM -XRD) (Fig. 5). All previously described phases are present in all samples, irrespective of facies and fracture density. It is thus possible to group phases I, II, and III, as prior-to-fracture diagenetic events.

Phase (IV) is only described for fractured facies, where several later diagenetic processes occurred. For samples from low fractured blocks, post-fracture illite precipitations are observed in pores. This morphotype of illite (4a) is fibrous with a completely different morphology to illites ( 1), ( 2), and (3) and is not combined with any hydroxides.

For highly fractured samples, feldspar and quartz cements are partly dissolved, forming their current border morphology. Cements present combined micro-fractures and dissolved zones. In these cavities a later precipitation of small sheets of illite crystals, a morphotype of illite (4b), occurs. In several samples from highly-fractured domains, generation (4a) of illite has also been observed. Unfortunately, no geometrical arrangement allowed the discrimination of which of generations (4a) or (4b) occurred first. A latter pore-infilling kaolinite is observed at a very low quantity (<3 %) in 4 samples containing a macro-fracture.

'Insert Figure 10'

Microstructural and petro-physical property relationships

To better decipher the interdependency between fracture distribution visible at a macroscopic scale, pore network features, microstructural characteristics, and permeability, a Principal Component Analysis (PCA) is performed. This concerns two datasets: Fa1 and Fa2 samples (Fig. 11). The main conclusion drawn is that the control of pore network properties ensured by the structure of the damage zone is not directly linked to the fault core distance. However, at the macroscopic scale, strain in the damage zone induces the development of corridors with decreasing fracture intensity, lowering the thickness of a corridor with its distance from the fault core. The control scale of petrophysical properties depends on the organization within the corridors, including a dependency toward the distance from the closest fracture. Fracture and crack intensities increase toward the fault core, but the distance from the fault core is not a controlling factor for microstructures or matrix permeability. Effectively, PCA

shows that there is a poor correlation between the distance to the fault core and pore dimension (S1: Micro-pore throat size, S2: Meso-pore throat size, P1: Micro-pore void size, P2: Meso-pore void size).

The petrophysical parameters can be divided into two categories, with specific microstructural controls operating on each:

The fissuration ratio, permeability, dynamic Young's Modulus, % of illite all have a property relationship with linked behavior for both facies. Thus, their evolution is only controlled by the distance from the closest fracture. The fissuration ratio is positively correlated with an increase in permeability. This suggests that an increased number of open discontinuities toward the fracture plan facilitates the connectivity of the crack network, and an improved permeability (Fig. 12). Furthermore, the decreasing percentage of illite simplifies the percolation of fluids through the intergranular voids, which also improves the permeability.

'Insert figure 11' 'Insert figure 12'

For both facies, porosity, cementation rate, throat diameters (S1 and S2), and micropore void diameters (P1) show different behavior depending on the facies considered.

These properties are dependent on the distance from fracture and facies, and they influence the type of pore network and permeability.

For Fa1, the pore network dimensions increase (P1, S1, P2, and S2) when permeability increases (Fig. 13). Pore voids and throats are affected by cataclastic deformation. The increase of crack intensity favors the grain crushing and pore collapse observed in SEM imaging (Figs. 5E and5F). This reduces pore network dimensions and global porosity. For this facies, permeability is mostly controlled by the connectivity of cracks, as illustrated by the positive correlation between fissuration ratio and permeability.

For Fa2, permeability increases toward the fracture, and is associated with an increase in throat diameter (S1 and S2), and of micro-pore voids (P1). Different strain mechanisms explain the pore network evolution with the distance to the fracture. The P2 pore size decreases towards the fracture and this can be explained by pore crushing (as for Fa1). Here this process is counterbalanced by the widening of pore throats upon dissolution of cements around grains, affecting S1, S2 and P1

dimensions. This is illustrated by SEM pictures (Figs. 6E and6F), and the decrease in the cementation rate close to a fracture. In fractured Fa2 samples, K-feldspar grains and cements are partly transformed into illite. K-feldspar and quartz cement rims have a developed system of etches and concavities, interpreted as chemical corrosion.

Permeability of Fa2 is controlled by connectivity (fissuration ratio and permeability correlate positively), but also by throat size widening.

From a microstructural point of view, the deformation mechanisms are different for each facies (Fig. 12): For Fa1, the deformation tends to be expressed as a cataclastic flow, with strong grain crushing, lowered Young's Modulus (E), and reduction of pore show more strain than the undeformed Fa1c and Fa2c facies. The slightly higher ratio between Fa2 curves than Fa1 ones also suggests that the strain will be more developed for Fa2 than for Fa1.

'Insert figure 13'

Hydraulic and geo-mechanical properties at the damage zone scale

To understand the hydraulic and geo-mechanical properties at the damage zone scale, an intermediate scale must be introduced. This is the single fracture and its surrounding matrix scale, a link between the grain scale and the damage zone scale.

At this intermediate scale of the single fracture, the evolution of petro-physical and mechanical properties varies with the distance from the fracture (Fig. 13). It is possible to define a damage zone for each fracture plane, as deformation processes impact the grain scale features within a limited distance around the fracture (Fig. 14). For Fa1, with cataclastic deformation, the average limit of the fracture damage zone can be settled for this case study at around 0.3 m. Above this distance from the fracture plane, dynamic Young's Moduli are similar to areas far from the fracture plane. Pore crushing and clogging don't affect pore dimensions significantly. For Fa2, where fracturing is combined with dissolution, this limit of the fracture damage zone is far higher (above 0.8 m). Pore dimensions and dynamic Young's Moduli values are not completely stabilized, because dissolution affects a wider volume than grain crushing. For the two facies, this limit of the fracture damage zone varies because mechanisms of deformation and secondary diagenetic sequences produce two different porous networks. These have different controlling factors, one on permeability, and the other on geo-mechanical properties.

The changes in dynamic Young's Moduli with the distance from the fracture are also representative of geo-mechanical behavior. At the damage zone scale, it will be expressed by a preferential reactivation of fracture corridors when stress is applied.

The media tends to have the properties of fractured facies.

This suggests dependency between the thickness of the fracture damage zone and the deformation type and diagenetic processes affecting the matrix around a single fracture. Strain processes are not only open cracks and reduced mechanical properties, as suggested by [START_REF] Bauer | Architecture, fracture system, mechanical properties and permeability structure of a fault zone in Lower Triassic sandstone, Upper Rhine Graben[END_REF], but also cement dissolution and weakening of contacts between grains. Thus, diagenetic processes control microstructures of rocks and also affect Young's Modulus and P-wave propagation properties, as has already been observed in other studies [START_REF] Major | The effect of CO2-related diagenesis on geomechanical failure parameters; fracture testing of CO2-altered reservoir and seal rocks from a natural analog at Crystal Geyser[END_REF][START_REF] Faÿ-Gomord | New insight into the microtexture of chalks from NMR analysis[END_REF][START_REF] Korneva | Deformation mechanisms and petrophysical properties of chert and limestone fault rocks within slope-to-basin succession (Gargano Promontory, Southern Italy)[END_REF].

'Insert figure 14'

These interpretations can be integrated into a larger fault zone scale model (Fig. 15.A). Considering the fracture network corridor pattern, it is possible to combine this with the notion of fracture damage zones. The impact of the fractured corridors on compartmentalization of the damage zone results in permeability and geo-mechanical anisotropy. In this case study, three types of behavior can be differentiated: structural position, sedimentary inheritance and fluid/rock interaction:

1) In slightly fractured areas, fractures modify the pore network properties at the fracture scale, with permeability increase, and weakening of geo-mechanical resistance in a limited area around the fracture plan (Fig. 15B).

The fracture intensity decreases far from the fault core, causing an increase of the spacing between fractures. Thus a decrease in the coalescence of the fracture damage zones is evidenced, which directly affects the permeability.

Two types of behavior are thus identified for fracture corridors, depending on their proximity to the fault core.

2) Far from the fault core, the medium is more likely to present a behavior similar to that observed for samples close to fracture planes, as damage zones of single fractures interfere and enhance hydraulic properties (Fig. 15C). These freshly opened corridors have the following features: a highly permeable fracture corridor, with a high fissuration ratio. A medium fracture intensity, suggesting mixed fracturing and cataclastic flow deformation. In the case of Fa2 facies, cement dissolution increases micropore throat and void sizes, along with meso-pore throat size. Pore collapse mainly affects meso-pore voids as a side effect of cataclasic flow. These two mechanisms of deformation on sandstones had already been recognized in previous studies [START_REF] Ballas | Shear-enhanced compaction bands formed at shallow burial conditions; implications for fluid flow (Provence, France)[END_REF][START_REF] Saillet | Permeability and flow impact of faults and deformation bands in high-porosity sand reservoirs: Southeast Basin, France, analog[END_REF][START_REF] Soliva | A model of strain localization in porous sandstone as a function of tectonic setting, burial and material properties; new insight from Provence (southern France)[END_REF][START_REF] Farrell | Anisotropy of permeability in faulted porous sandstones[END_REF][START_REF] Farrell | Anisotropic pore fabrics in faulted porous sandstones[END_REF].

3) Close to the fault core, fracture intensity and connectivity increase [START_REF] Bauer | Architecture, fracture system, mechanical properties and permeability structure of a fault zone in Lower Triassic sandstone, Upper Rhine Graben[END_REF]. In these corridors, relationships of permeability versus fracture intensity at the fault scale remain valid, with a bulk permeability of the media increased in high fracture intensity corridors. Fracture network in these corridors is due to proximity to the fault core, and the geo-mechanical properties of the fracture walls. These features lead to a high to medium permeability in fracture corridors, with high fracture intensity (structural inheritance). During reactivation, deformation becomes more cataclastic. Grain crushing is increased, with pore collapse affecting meso-pore voids and throats, and illite clogging is a side effect of fracturing. This explains the slight decrease in permeability, despite the increase of fracture intensity for corridors that touch the fault core.

'Insert figure 15'

Conclusion

Understanding the hydraulic behavior of fractured sandstones in a normal fault damage zone requires the precise description of the fracture network at different scales. The fault-model proposed here is one with a strong impact of interactions between fracture network and pore network features on the overall permeability estimation. To access the pore network properties, an intermediate scale of observation was used to estimate the interaction volume between a single fracture and its surrounding rock-volume. For both facies studied, this threshold varies because mechanisms of deformation and secondary diagenetic sequences produce two distinct porous networks. This is consistent with a previous study that quantified the pore fabric of deformed sandstones [START_REF] Farrell | Anisotropic pore fabrics in faulted porous sandstones[END_REF]. Different behavior of the corridors has been defined, due to differential rock properties (initial facies, fracture intensity, and diagenetic processes). All these factors control reservoir compartmentalization and increase the complexity of fluid flow pathways in the damage zone. In this case study, the conduit behavior of the damage zone in sandstones is favored within fracture corridors [START_REF] Eichhubl | Structural and diagenetic control of fluid migration and cementation along the Moab fault , Utah[END_REF][START_REF] Williams | Diagenetic controls on the evolution of fault-zone architecture and permeability structure : Implications for episodicity of fault-zone fluid transport in extensional basins[END_REF], with a dynamic architecture also influenced by a sedimentary heritage and pre-deformation diagenesis (consistent with conclusions from [START_REF] Farrell | Anisotropic pore fabrics in faulted porous sandstones[END_REF] for the pore fabric aspects, and with those of [START_REF] Knipe | Faulting , fault sealing and fluid flow in hydrocarbon reservoirs : an introduction[END_REF]; [START_REF] Odling | Permeability scaling properties of fault damage zones in siliclastic rocks[END_REF]; [START_REF] Pei | A review of fault sealing behaviour and its evaluation in siliciclastic rocks[END_REF]; [START_REF] Peacock | A broader classi fi cation of damage zones[END_REF] at the damage zone scale). This approach and the issuing conceptual model of the fault damage zone permeability pattern can be used in a double porosity model, for these types of reservoir or analogues, involving a strong knowledge of structural and pore network features. Graben, showing major variscan and cenozoic faults [START_REF] Schumacher | Upper Rhine Graben: Role of preexisting structures during rift evolution[END_REF]. BL:

Figures

Badenweiler -Lenzkirch fault system; LB: Lalaye-Lubine-Baden-Baden fault system; SHB: South Hunsrück Taunus border fault system.

B) Zoom of the area of study, in the northern part of the Saverne facture field (Saverne FF). White dots are main towns in the area, and the grey star shows the site on an N150 fault of the outcrop studied, connected to the Rhenish fault, on the western border of the Upper Rhine Graben [START_REF] Bauer | Architecture, fracture system, mechanical properties and permeability structure of a fault zone in Lower Triassic sandstone, Upper Rhine Graben[END_REF]. The red line represents the traverse line of the simplified cross-section presented on C. C) Simplified cross-section of the western rift shoulder in the Cleebourg area [START_REF] Ménillet | Notice explicative, Carte géol[END_REF]. [START_REF] Schumacher | Upper Rhine Graben: Role of preexisting structures during rift evolution[END_REF][START_REF] Blaise | Reconstructing fluid-flow events in Lower-Triassic sandstones of the eastern Paris Basin by elemental tracing and isotopic dating of nanometric illite crystals[END_REF], and regionally known hydrothermal events (he) associated with major fault zones [START_REF] Cathelineau | Dating multistage paleofluid percolations: A K-Ar and O -18 /O -16 study of fracture illites from altered Hercynian plutonites at the basement/cover interface (Poitou High, France)[END_REF][START_REF] Cathelineau | A major Late Jurassic fluid event at the basin/basement unconformity in western France: 40 Ar/ 39 Ar and K-Ar dating, fluid chemistry, and related geodynamic context[END_REF][START_REF] Clauer | Thermal history of Triassic sandstones from the Vosges Mountains-Rhine Graben rifting area, NE France, based on K-Ar illite dating[END_REF][START_REF] Blaise | Reconstructing fluid-flow events in Lower-Triassic sandstones of the eastern Paris Basin by elemental tracing and isotopic dating of nanometric illite crystals[END_REF]. The paleo-stress orientations are represented in black-σ1, shaded -σ2, or white -σ3. For each facies, the increase of the fissuration ratio is correlated with the increase of permeability. On the global dataset, this log-relationship is not clear, because it is influenced by intrinsic properties of each facies.

B) Permeability versus Young's Modulus. The two facies seem to have similar Young's Moduli for high to very high permeability samples. The lower Young's Modulus is, the more permeable the sample.

Figure 13: Evolution of porous network dimensions with increasing distance from the macro-fracture, with best line data fits (BF) for each property and each facies.

A) Throat size variations with distance from the fracture (S2 (black) -large pore throats; S1 (red) -small pore throats); for Fa1 (triangles), a weak decrease of throat size with vicinity from the fracture is seen, while for Fa2 (squares), there is an increase in both ranges of throats with proximity to the fracture plane: enlargement through dissolution occurs.

B) Void size variations with distance from the fracture. (P2 (black) -large pore voids; P1 (red) -small pore voids); for Fa1, the pore size is decreased (crushing and pore collapse, clogging by enhanced illite). For Fa2, meso-pore size is decreased by same processes as for Fa1, but an increase in micropore size occurs due to enhanced dissolution.

C) Permeability evolution with distance from the fracture.

.

4. 1 .

 1 Fracture network in the damage zone 4.1.1.Fracture network at outcrop scale (Fracture P10)

  ) starts in direct contact with a slip surface delimiting the outer fault core. This corridor starts at the distance of 19 m, and finishes at the distance of 21 m from the inner fault core datum. The deformation is the most developed in this corridor, with up to 12 fractures/m, and an average of 10 fractures/m. This fracture corridor contains two minor antithetic faults, and one minor synthetic fault, relative to the major fault. Displacements of the minor faults ranges from 2 to 6 cm. The second corridor, located between distance 24 m and distance 30 m contains intervals with a fracture density of between 2 and 11 fractures/m, and average density of 5.85 fractures/m. This corridor also includes a minor fault, synthetic to the major fault core. The third corridor, located between distance 37 m and distance 44 m has a fracture density from 1 to 8 fractures/m, and an average density of 4.28 fractures/m. Connectivity between fractures increases toward the fault core and is enhanced in the corridors (with more than 50% of connected fractures). The fracture corridors have an average lateral thickness of 5 meters, and the distance between corridors tends to decrease toward the fault core.'Insert Figure4here'4.1.2.Crack intensity as number of cracks (/m) (P10)P10 Crack intensity values range from 3.9x10 5 to 65.7x10 5 cracks/m (Fig.4.C). The crack intensity is higher for samples from fractured corridors than in samples from low fracture intensity blocks, where intensities are less than 10x10 5 cracks/m. Three highcrack intensity corridors are identified on both crack and fracture intensity profiles (Fig.4C).4.1.3.Total crack length (m/m²) (P21)Total crack length intensity (P21) values range from 1345 to 18530 m/m ² . Three highfracture intensity corridors are recognized, with P21 above 2000 m/m ² .Crack intensity distributions (P21 and P10), represented in Figure4C, show the same trend, with slightly higher peaks for P21, in samples where cracks are well-developed and propagated through longer distances, with intergranular cracks. The envelopes for crack and macro-fracture P10 intensities are similar.4.2. Damage zone characterization:4.2.1.PetrographySamples can be divided in two lithofacies (Fa1 and Fa2) with respect to their mineralogical and pore network features. Both petro-facies are present in the two strata sampled, equally in low fractured areas and in fracture corridors.

  network dimensions. The Fa2 facies shows a mixed brittle and cataclastic flow behavior (low (E) near fracture planes), with fractures and cracks remaining open. For the same fissuration ratio, Fa2 samples have higher permeabilities than Fa1, suggesting a preservation of open discontinuities. This type of deformation allows fluid flow and dissolution of cements, leading to the modification of the pore fabric. For the same stress, Fa1f and Fa2f fractured facies have a lower Young's Modulus, and will

Figure 1 :

 1 Figure 1: Regional geological localization

Figure 2 :

 2 Figure 2: Chronostratigraphic chart of the Upper Rhine Graben area (modified after

Figure 3 :

 3 Figure 3: Methodology applied to compare pore network parameter measurements. A. NMR T2 raw data distribution. B. MIP throat size distribution, which provides the distribution of pore throat sizes. C. NMR T2 distribution is converted into a distribution of void size radius compared to porosity. D. Combined cumulated porosities measured by MIP (black curve) and X-ray tomography (black dots).E. Combined data from graphs B, C and D. The largest void is determined by Xray tomography and is used to verify that the conversion from T2 to the void radius is calibrated. Each void radius (VR) peak can be associated to its corresponding throat radius (TR). Thus the features of each pore family are quantified.F. Representation of parameters measured. Lines are delimited by squares, indicating the void size radius measured by NMR, and grey lines are delimited by white triangles showing the throat voids, which are measured by MIP.

Figure 4 :

 4 Figure 4: Cleebourg's outcrop characterization (see location on Figs. 1.B and 1.C). A) Outcrop panorama, with sampling location and fault zone structure shown, along with damage zone, outer and inner fault core. The last contains the datum point 0. B) Drawing of the panorama, with sample location, facies distribution and corridor positions. C) Fracture and crack density profiles along the damage zone, from distance 19 to 50 m from the inner core zone. Fracture corridors are shown in pink. Triangles and squares are the locations of samples along the profile.Fracture intensities and connectivities are shown as grey, light green and dark green bars. Data are from[START_REF] Bauer | Architecture, fracture system, mechanical properties and permeability structure of a fault zone in Lower Triassic sandstone, Upper Rhine Graben[END_REF]. Crack intensities (P21 and P10) are from thin section observations.

Figure 5 :

 5 Figure 5: Fa1 facies thin sections A) Low crack density Fa1 facies (Plane polarized light): Quartz cement bound by automorph face (white arrow), Feldspar cement rims (black arrow). B) Low crack density Fa1 facies (cold cathode cathodoluminescence image): Quartz cement bound by automorph face (white arrow), Feldspar cement rims (black arrow). C) High crack density Fa1 facies (Plane polarized light):Quartz cement (plain white arrow), Illite clogging pores, compaction contact between grains (dashed white arrow) D) High crack density Fa1 facies (cold cathode cathodoluminescence image): Quartz cement rims are partly dissolved. E) High crack density Fa2 (Back-scattered electron microscopy imaging): Quartz grain and generation Ilt (2) of very iron oxide rich illite, Iron oxide rich mesh illite (Ilt (3)), Ribbon and plate later generations of illite (Ilt 4b) and (Ilt 4a). Black rectangle represents the detailed zone F. F) Detail back-scattered electron microscopy imaging: Cracks (red) affect the detrital grain and generations 2 and 3 of illite Generations 4 of illite are not easily identified.

Figure 6 :

 6 Figure 6: Fa2 facies thin sections A) Low crack density Fa2 facies (Plane polarized light): Quartz cement rims automorph faces (white arrow), Feldspar cement rims (black arrow). B) Low crack density Fa2 facies (cold cathode cathodoluminescence image): Quartz cement rims automorph faces (white arrow), Feldspar cement rims (black arrow). C) High crack density Fa2 facies (Plane polarized light): Cement rims are partly dissolved (quartz -broken white arrow; K-felspar -plain red arrow). D) High crack density Fa2 facies (cold cathode cathodoluminescence image): Cement rims are partly dissolved (quartz -broken white arrow; K-felspar -plain red arrow). E) High crack density Fa2 (Back-scattered electron microscopy imaging), with dissolved Feldspar grains and cement rims (red arrows).F) Back-scattered electron imaging: Illite generations, pre-(broken white arrow) and postfracture (solid white arrow), with various morphologies and iron hydroxide contents (variable density of greyscale). The dissolved Feldspar grains and cement rims (red arrows) are also present.

Figure 7 :

 7 Figure 7: Fa1 porous network properties: incremental NMR, X-ray tomography

Figure 8 :

 8 Figure 8: Fa2 properties: See captions of figure 7 for other explanations.

Figure 9 :

 9 Figure 9: Dynamic and static moduli for the Cleebourg's damage zone from P and S wave measurements (this study), compared with Bauer et al., (2015). The fracture corridor locations are shown in pink.

Figure 10 :

 10 Figure 10: Diagenetic sequence for Cleebourg's samples

Figure 12 :

 12 Figure 12: Permeability versus mechanical moduli A) Permeability versus fissuration ratio relationship.

Figure 14 :

 14 Figure 14: Schematic figure showing fracture damage, with areas where deformation

Figure 15 :

 15 Figure 15: Fault damage zone schematic profile.

Table 2

 2

	),

pore size of 8 to 20 µm. This family represents 24 to 28 % of the total porosity; (S2) has a throat diameter of 16.15 to 22.59 µm; (P2) a pore size of 55-105 µm. For Fa1, this family represents 72 to 76 % of total porosity.

Table 1 :

 1 Mineralogical composition of Cleebourg's samples, expressed in percentage of total surface area on thin section, and determined by computer assisted point-counting (minimum 300 per picture, 4 pictures per sample).

			Sample	Porosity (%) Illite (%) Iron oxides (%)	Kaolinite (%)	K-feldspar (%) Quartz (%)
	Fa1	undeformed	K466 K300 K500	13.2 12.8 12.3	7.6 7.9 7.8	3.2 2.1 1.5	1 2 1	11.9 8.2 17.8	63.1 67.0 59.7
			K430	12.3	6.3	3.1	0.1	15.1	63.2
	Fa1 fractured	K444 M386 M195	11.4 13.1 12.0	6.1 6.3 7.1	2.9 3 1.2	2 0.5 1	16.1 10.1 18.1	61.5 67.0 60.6
			M423	11.8	6.9	3.1	0.5	14.9	62.8
	Fa2	undeformed	M225 M352 M340	17.0 14.5 14.8	4.1 3.8 4.9	2.1 1.1 2.9	0.5 2 2	12.1 16.1 8.2	64.2 62.5 67.2
			K390	14.7	1.95	1.1	0.5	10.4	71.4
	Fa2 fractured	K393 M475 M250	14.3 18.4 18.7	2.2 1.9 2.9	1.8 0.9 1.9	1 3 2	9.7 12.1 10.1	71.0 63.7 64.4
			M295	16.6	4.8	1.8	0.5	14	62.3

Table 2 :

 2 Pore network properties measured on Cleebourg's samples

			Sample	Distance to fracture (m)	Porosity MIP (%)	Permeability (mD)	Average mesopore diameter (µm)	Average meso-throat diameter (MIP) (µm)	Average micropore diameter (µm)	Average micropore throat diameter min (µm)	Average micropore throat diameter max (µm)	Average micro-throat diameter (MIP) (µm)	Ratio of micro-porosity on total porosity
	Fa1 undeformed	K300 K466 K500	0.3 0.41 1	12.8 13.18 13.4	33 60 28	55 90 105	(P2)	16.15 22.26 22.59	(S2)	8 15 20	(P1)	0.42 0.86 0.69	3.06 5.024 4.2	1.14 2.39 1.65	(S1)	0.28 0.236 0.252
	Fa1 fractured	M195 M386 M423 K430 K444	0.2 0.05 0.25 0.08 0.08	12 13.06 11.83 12.25 11.43	42 235 58.5 157 220	45 40 50 45 38	(P2f)	9.85 12.57 11.59 11.62 9.93	(S2f)	12 10 12 12 5	(P1f)	0.06 0.029 0.09 0.33 0.011	6.68 3.05 3.06 3.93 4	1.45 1.54 1.13 1.87 0.69	(S1f)	0.774 0.44 0.25 0.38 0.572
	Fa2	undeformed	M225 M340 M352	0.8 0.8 0.85	16.98 14.8 14.5	164 100 145	120 135 120	(P2)	25.5 18.8 25.2	(S2)	8.2 7.1 5.9	(P1)	0.2 0.022 0.03	1.46 3.08 1.23	0.54 0.66 0.45	(S1)	0.16 0.38 0.12
	Fa2 fractured	M250 M295 K390 K393 M475	0.25 0.7 0.15 0.15 0.2	18.7 16.6 14.66 14.29 18.44	340 139 554 512 845	85 130 60 65 45	(P2f)	29.5 22.3 34.5 30.6 29.5	(S2f)	11.2 11.2 17.1 15 15	(P1f)	0.05 0.001 0.42 0.42 0.02	2.15 2.38 6.511 6.5 5.12	1.8 1.23 3.08 3.1 3.17	(S1f)	0.25 0.44 0.15 0.12 0.18

Table 3 :

 3 Acoustic properties measured and calculated on Cleebourg's samples

			Sample Vp (m/s)	Vp* (m/s)	Vs (m/s) IC (unitless)	ICp (unitless)	Poisson's ratio (unitless)	Fissuration Ratio (unitless)	Young's Modulus (Dynamic) (GPa)
	Fa1	undeformed	K300 K466 K500	2654 2684 2741	5328 5299 5242	1350 1460 1420	49.81 50.65 52.29	82.08 81.55 81.24	0.325 0.29 0.317	0.39 0.38 0.36	8.02 8.56 8.62
	Fa1 fractured	M195 M386 M423 K430 K444	2541 2265 2612 2438 2382	5318 5301 5362 5383	1210 1240 1358 1240 1240	47.78 42.73 48.71 45.74 44.25	83.20 81.72 83.44 82.85 84.00	0.353 0.286 0.315 0.326 0.314	0.43 0.48 0.42 0.45 0.47	7.08 6.29 7.78 6.71 6.69
	Fa2	undeformed	M225 M340 M352	2589 2900 2845	5032 5202 5160	1440 1450 1475	51.45 55.75 55.14	76.23 79.28 79.70	0.276 0.333 0.316	0.33 0.30 0.31	7.86 9.25 9.30
	Fa2 fractured	M250 M295 K390 K393 M475	2241 2684 2367 2440 2150	4934 5048 5153 5192 4929	1280 1360 1240 1250 1410	45.42 53.17 45.93 47.00 43.62	73.82 76.76 79.48 79.99 74.18	0.258 0.327 0.311 0.322 0.123	0.38 0.31 0.42 0.41 0.41	5.89 8.06 6.05 6.61 5.37

Acknowledgements

This work was supported by a PhD grant from ENGIE and the University of Lorraine.

We would like to thank Rauscher Company for access to the quarry and sampling authorizations. All the GeoRessources and LEMTA laboratories and ENGIE collaborators are also thanked for numerous advice, settings supplies and help with measurements. Special thanks are due to Dr. Sebastien Leclerc from LEMTA, Dr. Marc Diraison from GeoRessources, and Dr. Andrei Lecomte and Christophe Morlot, from SCMEM (Université de Lorraine), for SEM and X-ray tomography data treatment respectively. We are also grateful to the editor Dr. Toru Takeshita, the reviewers Dr.

Fabrizio Agosta, Dr. Nathalie Farrell, and the third anonymous reviewer for their time and suggestions for improvement of the manuscript.