A laboratory approach for the calibration of seismic data in the western part of the Swiss Molasse Basin: the case history of well Humilly-2 (France) in the Geneva area

A collection of 81 plugs were obtained from the Humilly-2 borehole (France), that reached the Permo-Carboniferous sediments at a depth of 3051 m. Experimental measurements of physical parameters and mineralogical analysis were performed to explore the links between sedimentary facies and seismic characteristics and provide a key tool in the interpretation of seismic field data in terms of geological formations. The plugs, cylinders of 22.5 mm in diameter and ∼\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}30 mm in length were collected parallel and perpendicular to the bedding in order to explore their anisotropy. Ultrasound wave propagation was measured at increasing confining pressure conditions up to 260 MPa, a pressure where all micro-fractures are considered closed. The derivatives of velocities with pressure were established, allowing the simulation of lithological transitions at in-situ conditions. At room conditions, measured grain densities (kg/m3) range from 2630 to 2948 and velocities vary from 4339 to 6771 and 2460 to 3975 m/s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text {m/s}}$$\end{document} for P- and S-waves propagation modes, respectively. The largest seismic-reflections coefficients were calculated for the interface between the evaporitic facies of the Keuper (Lettenkohle) and the underlying Muschelkalk carbonates (Rc=\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{\text {c}}=$$\end{document} 0.3). The effective porosity has the range of 0.23% to 16.65%, while the maximum fluid permeability m2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( \text {m}^{\text {2}}\right)$$\end{document} is 9.1 ×10-16\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times \text {10}^{-{\text {16}}}$$\end{document}. A positive correlation between porosity and ultrasound velocity has been observed for P- and S-waves. The link between velocities and modal content of quartz, dolomite, calcite, and micas has been explored. This paper presents a unique set of seismic parameters potentially useful for the calibration of seismic data in the Geneva Molasse Basin.


Introduction
Seismic reflection imaging is a fundamental tool for evaluating the geometry of the deep reflectors within the Earth's crust (Brown 1986;Claerbout 1985;Schilt et al. 1979). Nevertheless, seismic velocity distribution at depth cannot be resolved only with the aid of seismic reflection methods and other independent geophysical methods need to be applied. Furthermore, if the formations that compose the geological section under scrutiny have a similar mineralogical composition and/ or if bedding and layering are not pronounced, imaging the geological sequence by using the reflection seismic method can be quite problematic. One possibility to help understanding the origin of seismic reflectors is to use the laboratory experiments at ultrasonic frequency and under high-confining pressure, which simulate insitu subsurface conditions. This method has been successfully applied in numerous studies and in different environments (e.g. Almqvist et al. 2013; Barberini et al. 2007;Birch 1960;Wenning et al. 2016;Zappone et al. 2000). In the present work we apply this approach for the geological interpretation of the seismic data recorded in the Southwestern part of the Swiss Molasse Basin.
The Geneva Basin (GB), in the westernmost part of Swiss Molasse Basin, has been the object of intense seismic reflection exploration for hydrocarbon (Gorin et al. 1993). Between 1958 and, about 300 km of deep reflection seismic have been acquired, and a 16 km long, composite seismic transect was presented (Gorin et al. 1993;Moscariello 2016). Following the unsuccessful hydrocarbon exploration campaigns, interest in the GB has been renewed by geothermal investigations. A geothermal exploration program in Canton Geneva aims at producing geothermal energy from hot water in Mesozoic sediments (Moscariello 2019). In order to do that, the project has to evaluate the potential reservoirs through the integration of available geological and geophysical data. It focuses on the area stretching from the folded Jura Mountains in the north to the Mount Salève ridge in the south (Moscariello 2016).
More recently a transnational French-Swiss project within the European project GeoMol (2014-2021) developed new studies over an area of 2800 km 2 in the Savoy-Geneva region, and used 985 km of available 2D seismic dataset (800 km over French territory, and 185 km over the Swiss territory) together with logs of 43 available deep boreholes loggings (including bottom hole temperature; BHT) to reconstruct a 3D geological model and temperature distribution model for the sedimentary basin. Available logs from the wells Humilly-2 and Thônex-1 were used to calibrate stratigraphic boundaries with seismic horizons in the GB (Fig. 1). Thônex-1 is a deviated, geothermal exploration well drilled in 1993. It reached a total depth of 2530 m in marly Jurassic limestones (Jenny et al. 1995). The vertical Humilly-2 (HU2) borehole is the deepest well drilled in the GB during a 1960s hydrocarbon exploration campaign in nearby France. Among all boreholes penetrating the entire Mesozoic succession in the area, HU2 is the only one where several stratigraphic intervals were cored down to a total depth of 3051 m in the top of the Permo-Carboniferous ( Fig. 1; Marti 1969). HU2 (46 • 6′ 53″ N and 6 • 1′ 30″ E) is located on the composite seismic profile shown in Fig. 1 (see also Sommaruga et al. 2012). Well logs, including gamma ray and sonic, were run over the entire section penetrated in HU2 (Fig. 2) and used to tie with the seismic data.
Six regional reflectors are commonly interpreted in the GB (Gorin et al. 1993;Clerc et al. 2015;Rusillon 2018). From top to bottom, they are: top Urgonian (base Tertiary), top argillaceous Upper Jurassic, top Middle Jurassic, top Lower Jurassic, top evaporitic Triassic and top carbonate Triassic (Fig. 2). The basal Triassic unit of the Buntsandstein is too thin to be seismically interpreted. HU2 reached total depth in Permo-Carboniferous sandstones characterized by a transparent seismic facies. The latter sequence is often underlain by an interval of high-amplitude negative reflections interpreted as coalbearing Permo-Carboniferous sediments (Gorin et al. 1993). Where this reflective sequence comes close to the Triassic, it highlights the angular unconformity between Triassic and Permo-Carboniferous sediments (Signer and Gorin 1995).
This paper summarizes the results of elastic wave velocity measurements combined with both thin section observations in optical microscopy and quantitative micro-textural analysis using the automated QEMSCAN technique, carried out on core samples from HU2. The purpose of this paper is to quantify the effect of porosity, mineral composition and micro-textural characteristics, such as banding and layering of minerals, of sedimentary rocks on seismic properties. The final aims of this paper are twofold: (1) to provide a potential contribution to the interpretation of seismic reflectivity zones in a 3D seismic survey to be acquired for geothermal exploration in the GB; (2) to investigate the possible use of these measurements in the calibration of borehole logs when cores are missing.

Regional geology
The Molasse Basin (MB) is located at the northern front of the Alps. It stretches from France (Savoy) through Switzerland up to southern Germany (Bavarian foreland) with a regional WSW-ENE trend and connects in Austria with the Vienna and Pannonian Basins (Fig. 1). The MB is a Cenozoic foreland basin which formed as a response to flexural subsidence during the collision of Eurasia and Adriatic-Africa plates (Pfiffner 1986;Véron 2005). Its substratum comprises Mesozoic (locally Permo-Carboniferous) sediments and crystalline metamorphic (See figure on next page.) Fig. 1 Geological map and structural framework of the western Swiss Molasse Basin (modified after Clerc et al. 2015;Makhloufi et al. 2018). The seismic profile AB (located on the map) crosses the Geneva Basin and intersects the Humilly-2 well location; it shows the main stratigraphic units. The geological units have been produced by time-to-depth conversion using the available seismic transect and the seismic velocity derived from the sonic log of Humilly-2 borehole (after Gorin et al. 1993;Moscariello 2016) basement rocks. The MB is divided into two geomorphological entities; the deformed and undeformed sub-units. The undeformed "platform" Molasse comprises 90% of the basin, while the folded and deformed "subalpine" Molasse is located along the southernmost flank of the basin (Véron 2005).
The Swiss Molasse Basin (SMB) occupies the western and central part of the MB, stretching over approximately 300 km between Savoy in the SW (30 km wide) and Lake Constance in the NE, where it reaches a width of 80 km. The Savoy-Geneva area is located at the southwesternmost end of the SMB. It is the result of a geological evolution initiated in the Palaeozoic (Schegg et al. 1999, and references therein), forming a narrow (max. 40 km in width) elongated ( ∼110 km in length) NNE-SSW depression. This depression is limited to the northwest by the Jura fold belt, to the southeast by the thrusted anticlinal ridge of the Salève and to the west by the Vuache anticline. The Salève ridge is intersected by several NW-SE trending sinistral wrench faults, which acted as transfer faults during thrust tectonics and influenced the distribution of the Tertiary deposition within SMB (Fig. 1).

Swiss Molasse Basin evolution
The present-day configuration of the SMB and its substratum underwent four major tectonic events (Schegg et al. 1999, and references therein): (a) Variscan orogeny and the subsequent formation of a series of narrow, deep, SWS-ENE and WNW-ESE trending Permo-Carboniferous grabens, filled by clastic sediments (Thury et al. 1994;(b) extension related to the sag-basin development and thermal subsidence during the Mesozoic (Loup 1992; (c) closure of Late Cretaceous Alpine flysch basins and development of the North Alpine Foreland Basin as a flexural depression during the late Cretaceous-early Miocene (Lihou and Allen 1996; detachment of Mesozoic and Cenozoic sediments on Triassic evaporites, deformation and uplift/inversion during the Miocene to Pliocene (Schegg and Leu 1996).
The tectonic elements responsible for the Geneva basin structuration are controlled by four major NW-SE trending wrench-fault systems. Their continuations across the Geneva Basin were revealed on the 2D seismic data and link the subalpine units with the Jura Mountains across the Bornes Plateau, the Salève ridge and the Geneva Basin (Clerc et al. 2015, and references therein). Detailed information of the basin evolution, and about the depositional environments are described in Brentini (2018).

Stratigraphy of the GB
The stratigraphy of GB consists of a 3-5 km thick sequence of Mesozoic and Cenozoic sedimentary rocks overlying Permo-Carboniferous sediments and Variscan crystalline basement. The Mesozoic and Cenozoic sequence dips gently (1°-3 • ) southeastwards towards the Alpine foreland (Fig. 2). A detailed study of the stratigraphic framework at basin-scale is presented in Brentini (2018).
The basement is mainly composed of biotite-rich gneiss, green schist and porphyritic granite resulting from the Palaeozoic Variscan orogeny (Matte 2001). Permo-Carboniferous siliciclastic sediments are locally present above this basement. They represent the infill of several SW-NE narrow Palaeozoic (mostly Hercynian) extensional basins developed in the crystalline basement (Wilson et al. 2004;McCann et al. 2006). They present a large variation in lateral thickness, depending on the location and depth of the troughs. The lithology consists of tight brownish conglomerates and arkosic sandstones, interlayered with silty, and organic-rich clays deposited mainly in a lacustrine environment in association with fluvial deposits resulting from the erosion of the crystalline basement. The Mesozoic succession represent the thickest part of the sediments. Generally, in the SMB the main lithologies (limestones and marls, except for the Triassic) are correlated to the Helvetic domain, which corresponds to the distal part of the European margin at the time of deposition. However, the GB shows proximal facies, which are associated with the Jura domain. The top of the basement is characterized by an angular unconformity on which Mesozoic sediments were deposited (Signer and Gorin 1995;Sommaruga et al. 2017).
The Triassic is characterized by a marine transgression leading to the deposition of thick evaporitic series in an epicontinental sea environment (Disler 1914). The Buntsandstein (Lower Triassic) is mainly composed of sandstone. Conglomeratic and silty layers are intercalated in this formation (c.a 15 m in thickness in HU2). It is overlain by the Muchelkalk (Middle Triassic), composed of marly limestones, anhydrites and dolomites and the Keuper (Upper Triassic) formed by a thick evaporite deposits (gypsum/anhydrite and salt). Given the ductility of salt, the Keuper constitute the basal décollement for the foreland thrust system (Diepolder et al. 2015) and generates Fig. 2 Stratigraphic column of well Humilly-2 with location of the potential reservoirs investigated. Location of investigated samples HU taken from cores and position of the six interpreted regional seismic reflectors. Cores were extracted in four depth intervals (I, II, III, and V in column "Core"). The depth interval IV is partially represented by samples from literature while the depth interval V corresponds to 3 sub-intervals: Permo-Carboniferous (VII), Lower Triassic (VI), and Middle Triassic (V) (see text for explanations). Anis. Anisian, Hauteriv. Hauterivian, A. Aleian, Batho. Bathonian, Muschelk. Muschelkalk (See figure on next page.) thrusts and duplexes, which produce large lateral variations in thickness (240 to 350 m; Brentini 2018).
The transgressive trend is still present during the Lower Jurassic leading to the deposition of limestones overlain by marls (Signer and Gorin 1995). The Middle Jurassic succession is characterized by a first regressive trend and a shift towards shallower conditions during the Bajocian and Bathonian (Strasser 1994). A second regressive trend affects the Upper Jurassic which exhibits the development of a shallow carbonate platform, extending towards the northwest. During the Kimmeridgian, existing structural highs saw the growth of patch-reefs (Meyer 2000;Makhloufi et al. 2018) while the sealing of inter-reef depressions by tidal deposits occurred during the Tithonian (Strasser 1994).
The Early Cretaceous saw shallow and warm-water conditions prevailing (Debelmas et al. 1961;Charollais et al. 2013) with low amplitude sea-level changes inducing several episodes of emersion and drowing characterized by the deposition of bioclastic limestones, bioturbated limestones and organic-rich marls. The Lower Cretaceous deposits are not recorded in the GB due to their emersion and later erosion during the early Cenozoic (Eocene-Oligocene). This event led to large-scale karstification of the Urgonian limestones. These karsts were later filled with Eocene lateritic sediments (Becker et al. 2013) and reworked Aptian-Albian Sediments. The Early Cretaceous limestones and Eocene lateritic sediments are overlain by alternations of continental sandstones and marlstones composing the Lower Freshwater Molasse. Both the Upper Marine and Upper Freshwater Molasse units are not preserved in the GB resulting from either the uplift of the Jura during the Miocene-Pliocene and/ or by the glacial advances of the Pleistocene (Signer and Gorin 1995;Schegg and Leu 1996;Charollais et al. 2007).

Methods
Twenty-three samples from the five aquifers have been recovered from HU2 borehole; multiple cylindrical plugs with 25.4 mm diameter and 29.9-52.4 mm length were cut from the samples. Fig. 2 shows the position of the samples within the stratigraphic succession. Plugs were cut parallel and perpendicular to bedding in order to study the anisotropy of the elastic properties in each rock type. Because at Humilly the layering is sub-horizontal and the borehole is not deviated, the direction perpendicular to bedding (Z) coincides with the drilling axis. One or two mutually perpendicular plugs (X and Y directions) were cut in the bedding plane (Fig. 3). In addition, two plugs were cored at the depth of 107.35-108.5 m from an anhydrite rich interval of the BEX borehole, Bohrung-3 (Canton Vaud). A lathe machine was used to guarantee planar and parallel surfaces at the ends of the cylindrical plugs, with ±5 µ m precision. The plugs were dried in an oven at 70 • C for 24 h to remove free water in pores and fractures, and stored in a desiccator, before performing further measurements.

Compositional and mineralogical characterization
Thin sections for optical microscopy were prepared from the heads of the cylinders. In order to distinguish calcite from dolomite in carbonate samples, thin sections were stained using a proportional mixture of Alizarin Red S and Potassium Ferricyanide, using the adapted methodology of Dickson (1965). All thin sections were examined and photographed under an optical microscope using a camera Nikon DS-Fi2 on a Nikon Eclipse E600 Polarizing Light Microscope (POL) equipped with a NIS Element D software package.
Optical microscopy analysis was complemented by scanning electron microscopy analysis performed on thin sections at the University of Geneva (QEMSCAN QUANTA 650F), that provided high-resolution mineralogical maps and automated quantitative mineralogical compositions identification. Mineral identification was based on the combination of back-scattered electron values, energy-dispersive X-ray spectra, and X-ray count rates. Scanning was performed at 15 kV with a resolution of 7.5 µ m on a 1 × 1 cm area. All these tools were integrated to assess the mineralogical composition, rock type, and vertical facies distribution. Orientation of the plugs relative to bedding and to borehole core. The arrows show the propagation direction of the seismic waves. Therefore V p (X) , V p (Y) and V p (Z) correspond to the velocities of compressional waves along X, Y and Z directions, respectively. For shear waves, we used the following convention: the two letters identify the polarization plane while the direction of propagation is indicated by the second letter. For example, V s (XZ) is the velocity of a shear wave polarized in the XZ plane and traveling in the Z direction

Density, porosity and permeability measurements
Bulk density was calculated by weighting the plugs after dehydration in oven for 24 hours, with a highly precise balance ( ±5.0 × 10 −7 kg tolerance), and determining the bulk volume of the plugs, using a caliper ( ±2 × 10 −5 m accuracy) to measure their length and diameter. The grain density was measured by using a gas displacement He-pycnometer apparatus (Accupyc II 1340, Micromeritics), that measures the matrix volume of the plugs at room conditions. The standard deviation of each measurement was lower than 5%. The effective porosity has been calculated as the difference between grain and bulk density. Fluid permeability measurements were performed using an oil-based hydrostatic pressure medium apparatus (permeameter) developed "in house" based on the transient step method (Brace et al. 1968). Using an automated permeameter-porosimeter (AP-608, Coretest Systems, INC.) at Geneva University, the fluid permeability (Y-plugs) measurements are performed based on the unsteady pressure drop method and Boyle's law (Jones 1972).

Seismic velocity measurements
Measurements of ultrasound waves propagation were performed at the Rock Physics and Mechanics Laboratory-ETH Zurich, following the pulse transmission method (Birch 1960(Birch , 1961 to measure separately the velocities of compressional-and shear-wave as a function of confining pressure (Fig. 4). The pulse transmission method, which is described in detail by Christensen (1985), is based on measuring the travel time of a mechanical impulse through a rock plug, placed between two piezoelectric transducers.
The experiments were performed under two successive cycles of pressurization-depressurization using an oil-medium hydrostatic pressure vessel. The confining pressure was generated using an air-driven fluid pump in conjunction with a compressed air system and was measured within ±2 MPa precision at the highest pressure. The pressure was incrementally increased up to 260 MPa and subsequently decreased with steps of ∼ 20 MPa. The experiments were performed at room temperature. The plug was placed between two lead zirconate titanate piezoceramic transducers, which were mounted into metal-heads, containing buffer rods, in order to optimize the transmission of mechanical energy towards the plug. The heads and the plug were encapsulated with a double polyolefin heat shrink tube (with an undrained pore space) to avoid contamination of the plug with the oil and to hold the plug in the correct alignment throughout the experiment. Separate Pand S-wave emitters were employed.
We used transducer with 1 MHz resonance frequency, therefore compressional and shear wavelengths were an order of magnitude larger than the grain size observed in the thin sections. The impedance used was 1 M and range was ± 500 mV. The waveforms were collected with a time step of 8 × 10 −8 s, and sampling rate of 25 MHz, and the electric noise was minimized by averaging the tracks. All acoustic waveforms during both pressurization and depressurization cycles were transferred to the computer for further processing without loss of resolution. Cables, transducers, and interfaces in the electronic system introduce a delay in the observed traveling time (i.e., t observed = t rock + t system ), therefore a calibration was performed before velocities measurements and repeated after the data collection phase.
The calibration was performed by measuring the travel time of the ultrasonic impulses (P and S) through steel cylinders of 24 mm, 29 mm, 39 mm and 44 mm in length, Fig. 4 Schematic setup of the seismic apparatus used to measure ultrasonic velocities up to 260 MPa confining pressure and at room temperature. The core assembly is accommodated inside a hydrostatic vessel, where the oil is the confining medium plotting the travel times against the cylinders lengths, interpolating linearly the results, and calculating the time value corresponding to zero length. Measurements were conducted at confining pressures of 50 MPa, 100 MPa, and 200 MPa. The traveling time through the system ( t system ) was obtained by averaging the results obtained for the three different pressures.

Permo-Carboniferous sediments (samples HU27 and HU29)
The Permo-Carboniferous is characterized by an alternation of sandstones and shales layers. The sandstone mineralogy exhibits kaolinite sometime affected by illitization (Fig. 5). The large crystal morphology observed are typical of burial diagenesis under high temperatures. Coal chips as well as palynological fossils are commonly found in these sandstone layers.

Lower Triassic (Buntsandstein siliciclastic unit; samples HU25 and HU26)
The Buntsandstein unit is a thin Lower Triassic sandstone ( Figure 2) showing many heterogeneities at the microscopic scale (Table 1). It is characterized by green, subrounded, coarse quartz grain. Associated with quartz, the mineralogical composition displays low amounts of K-feldspar (7%) and illite (4%; Fig. 6e). Rarely, quartz grains can present micro-fractures that are filled by diagenetic minerals (Fig. 6f ).

Middle Triassic (Muschelkalk and Lettenkhole units; samples HU17-HU24)
The Muschelkalk and Lettenkhole series present alternations of grey dolomitized grainstone with anhydrites nodules and dark, laminated, anhydrite layers. The base of the Muschelkalk is marked by a decrease in clastic sediment input and the presence of glauconite and evaporites. Two microfacies can be observed: (1) a dolomitized grainstone/packstone with peloids, ooids and anhydrite nodules displaying fluid inclusions and (2) peloidal laminated anhydritic layers (HU23, Fig. 6c, d). The dolomitic intervals are also characterized by frequent microcracks and sub-vertical fractures filled with anhydrite.

Upper Triassic (Keuper Formation; Samples BX107, CL126, and GDHlt)
The anhydrite-clay-salt interval (Keuper formation; 2548 to 2896 m depth) was originally not cored at Humilly. Nevertheless, we thought that such horizon could give interesting seismic results and we wanted to test it in laboratory. Therefore, in addition to the twenty-eight samples from HU2 borehole, we obtained an anhydrite sample from Bex-3 borehole in the Canton Vaud (depth of 107.35-108.5 m; Ultrahelvetic Zone, Bex-Laubhorn-Nappe, Switzerland). The horizontal plug shows vein networks filled by anhydrite (0.3-0.7 mm), cleavage surface and brecciation. For the clay component we considered literature data referring to sample CL126 (Freund 1992); it is a claystone from Salzwedel 2/64 borehole,

Middle Jurassic (Dogger limestone unit; samples HU07-HU12)
The Dogger unit is about 315 m thick in HU2 but the Liassic-Dogger transition is arbitrarily set at the clear distinction between the Upper Liassic unit of bioclastic Mudstone-Wackstone and the first presence of micarich limestones (Fig. 7c). The middle and upper part of  Sinnemurian-very fine-grained homogeneous micrite (> 91% calcite). A few of shell fragment (benthonic forams) which are preserved as a mold (MF) and filled with mostly diagenetic calcite, are embedded within very fine-grained homogeneous micrite. c sample HU21; Triassic Muschelkalk shows the microcrystalline micrite-dominated dolostone (fine-crystalline dolostone (Dol). d sample HU24; Triassic Muschelkalk presents the abundance of fibrous anhydrite (Anh) as the evaporite minerals in dolomitic (Dol) carbonate-support. e sample HU25; Triassic Buntsandstein. The kaolinization (Kln) of K-feldspar (Kfs) is occurring to form the clay cement for pore-filling and grain-coating on Quartz (Qz) minerals. Micro-fractures in quartz grains are not so common and filled by later diagenetic minerals. f sample HU29; Permo-Carboniferous unit. It shows slight laminations of fine-quartz (Qz) grains associated with kaolinite and biotite (Bt). The preservation of detrital feldspar and clay minerals (background) significantly contribute in the reduction of the porosity (0.87 ± 0.16%) the Dogger are mainly composed of bioclastic grainstone with dominant heterozoan assemblages, and sparse detritic quartz grains (maximum 1% of total allochem amount, Fig. 7d, e).

Upper Jurassic (Kimmeridgian Reef Complex unit; samples HU01-HU06)
The Upper Jurassic is more than 1000 m thick in HU2. The lithology is heterogeneous, ranging from marls to pure carbonates. This unit is often affected by various types of dolomitization from simple slight dolomitization in the form of fine dolomite crystals to fabric-destructive, sucrosic, dolomitization. The base of the Upper Jurassic displays a dark grey, bioclastic, marly limestone with few darker marl intervals. The Upper Jurassic limestone displays a white, massive, bioclastic limestone, with a 15 cm thick dolomitic interval at the base. The successive microfacies observed include at the base microbial boundstone, alternating with bioclastic grainstone to rudstone, whose bioclasts show large micrite envelopes that can sometimes prevent their identification (Fig. 7a, b) but are mainly composed of echinoderms, lamellibranchs, gastropods and foraminifers. Microbial constructions are then progressively replaced by coral framestone, still intercalated with bioclastic grainstone to rudstone. Moldic dissolution features, intra/inter-particle pores, and fractures are systematically filled with calcitic cement.

Modal composition, density, porosity and permeability
The results of the mineralogical mapping, using QEM-SCAN analysis are shown in Table 1 [the data available in Hefny et al. (2020)]. The mineral composition of the anhydrite-clay-halite section derived from literature (Freund 1992;Müller and Briegel 1977;Yan et al. 2016) are also reported in Table 1. Grain density, effective porosity, and fluid permeability data are given in Table 2. The full dataset is made available by Hefny et al. (2020). The Upper Jurassic section is characterized by monomineralic limestones while the Middle and Lower Jurassic limestones show a more complex mineralogy, with quartz in variable percentages up to 9%, and minor dolomite, illite and biotite. The Triassic dolomitic section of Muschelkalk/Lettenkhole present generally abundant Gypsum/Anhydrite, up to 54% in sample HU20, and minor quartz (around 1%). The siliciclastic units in the deepest part of the borehole are represented either by quartz rich sandstones (Buntsandstein), or fine grained to argillaceous quartz rich sediments (Permo-Carboniferous unit).
In the Upper Jurassic limestone Reef Complex unit, the grain density average is 2720 ± 19 kg/m 3 , reflecting the homogeneous mineralogical composition of the samples. Effective porosity ranges from 10.15 ± 0.4 to 16.65 ± 0.15% while permeability ranges from 3.9 × 10 −16 to 9.08 × 10 −16 m 2 , representing the highest values we could measure in our sequence. No relevant variation is observed between permeability measurements parallel or perpendicular to bedding. Porosity and permeability increase from top to the bottom of the Upper Jurassic limestone section.
In the Dogger limestone unit, the average grain density is 2714 ± 19 kg/m 3 , therefore similar to the limestones of the Upper Jurassic. The porosity is significantly lower than in the Upper Jurassic limestone Reef Complex, varying from 1.53 ± 0.15 to 2.72 ± 0.16%. The maximum permeability measured 3.6 × 10 −16 m 2 , was found in the X direction of sample HU07, in the upper part of the sequence, but most of the samples in the sequence show permeability values lower than the minimum value detectable with our instruments.
The Liassic shales at 2242.7 m depth has an average grain density of 2657 kg/m 3 , therefore slightly lower than the limestone in the same sequence, but equal porosity (2.1 ± 0.2%) and permeability below detectability. The Liassic limestones are also characterized by a low porosity; the maximum value is 2.42 ± 0.4%, observed along the bedding plane of sample HU15 (argillaceous limestone). Permeability is also very low: 10 −18 m 2 along the bedding plane and below the detection capability of our instrument perpendicular to bedding. Density (kg/m 3 ) values are typical of limestone and vary from 2678 ± 2 to 2723 ± 3. The mineral composition of the anhydrite sample (BX126) is 75.4% anhydrite, 15% magnesite, 5% quartz and 2% gypsum (Müller and Briegel 1977). Data reported in literature Gassmann et al. (1979) for density and porosity are respectively 2910 ± 10 kg/m 3 and 0.58 ± 0.2%.
The Muschelkalk-Lettenkohle anhydrite-dolomites show the lowest effective porosity of the whole samples set, with values ranging from 0.28 ± 0.11 to 1.65 ± 0.09%. The only value of the permeability we could determine was 4.0 ± 0.7 × 10 −19 m 2 corresponding to plug HU24X at a depth of 2967 m. The presence of abundant dolomite and anhydrite determines a relatively high grain density compared to the other units in the sedimentary sequence, with average values of 2887 ± 34 kg/m 3 .
The density values of Buntsandstein sandstones unit are lower than the overlying Muschelkalk formation, while porosity and permeability are higher. The average density is 2655 ± 18 kg/m 3 . Porosity ranges from 2.09 ± 0.27 to 4.50 ± 0.21%. The permeability values range from 0.4 ± 0.1 × 10 −18 to 54.0 × 10 −18 m 2 , with values perpendicular to banding one order of magnitude lower than the measured values in the bedding plane.

Seismic velocities
The seismic velocities were calculated using Eq. (1), where L is the plug's length, t rock is the traveling time of either P-or S-wave through the rock plug, and δV p,s is the uncertainty of V p or V p propagation through the experiments, respectively. No correction for length variation during pressurization was necessary, samples length was measured before and after pressurization to make sure that no permanent deformation was induced. The experimental traveling time of either P-or S-wave through the rock plug under loading/unloading cycles (i.e. up to 260 MPa confining pressure measurements) are available in Hefny et al. (2020). The first arrival times of P-and S-wave were manually picked using a MATLAB code developed for that purpose (Grab et al. 2015), with ±2 × 10 −9 second precision. Figure 8 is an example of a typical experimental curve of compressional velocity variation during a two pressurization/depressurization cycles. The seismic velocities-confining pressure curves show a typical non-linear relationship at lower confining pressure, typically below 150-180 MPa, (in the example of Fig. 8 it is 175 MPa), attributed to crack closure. The linear portion, at higher pressure, reflects the intrinsic seismic properties in the selected direction of the rock under consideration (Birch 1960(Birch , 1961. The linear regression of the data measured between 150 and 250 MPa was used to calculate velocity reference values at room pressure for P-and S-waves respectively. The intercept of the linear regression is referred as V p (0) and V s (0) , i.e. V p,s at zero pressure, and (1) V p,s = L t rock ± δV p,s , Table 2 Summary of the grain density; ρ grain (kg/m 3 ), effective porosity; φ eff (%), and permeability; k (m 2 ) for the studied plugs Standard deviations inside parentheses are expressed in units of the last digit. The permeability measurement performed at Geneva University are in italic NA not available, BL below limit (lower limit = 9.87 × 10 −20 m 2 ), NM Not Measured According to Burke and Fountain (1990), only measurements taken during depressurisation are reproducible within error limits. Therefore, we used only the depressurization cycles to calculate V 0 . We estimated the error ( δV p,s ) for each velocity by using Equ. (2) (Taylor 1997). (2) t observed − t system 2 δt system 2 ∂V p,s /∂L , ∂V p,s /∂t observed , and ∂V p,s /∂t system are the partial derivative of the seismic velocities with respect to the plug's length (L), the observed travel time ( t observed ), and the travel time through the system ( t system ), respectively. δL is the standard deviation of the plug's length, measured using a digital caliper with 2 × 10 −5 m accuracy. δt observed is the uncertainty in the picking of the first arrival on the waveform, and it is the main source of error in the determination of seismic velocity. We estimated a typical uncertainty in the picking both for P-and S-wave in the order of magnitude of δt observed ≃ ±2 × 10 −9 s. From the calibration of seismic rig (Sect. 3.3), we calculated a δt system = ±1.4794 × 10 −8 s, and ±1.9341 × 10 −8 for V p and V s , respectively. The uncertainty has nearly the same size as the symbol used in the Figures, at most ∼ 27 m/s, and therefore not visible in the diagrams. Velocity data are presented in Table 3. Average V p0 velocities (m s −1 ) range from 4897 to 5897 for limestone with the highest values in the middle Jurassic samples, from 5979 to 6648 for dolomite, and from 5303 to 5646 in sandstones. V s0 velocities (m s −1 ) range from 2759 to 3267 for limestone, from 3293 to 3737 for dolomite, and from 3336 to 3591 in sandstones. The only shale in our collection has a maximum value of 4775 and 2901 for V p0 and V s0 , respectively. The anhydrite rich sample from BEX has an average velocity of 4813 and 2859 for V p0 and V s0 , respectively .

Seismic anisotropy and Vp/Vs ratio
Velocities measured along mutually perpendicular directions, parallel and perpendicular to banding, were used for calculating the seismic anisotropy AV p,s (%) of the rocks as in Eq. (3).
(3) AV p,s = V max p,s − V min p,s V mean p,s × 100 ± δAV p,s , Fig. 8 Example of two cycles of pressurization-depressurization of compressional-waves (sample HU25X). The non-linear velocity/ pressure correlation at low pressure (up to ca. 175 MPa) is due to the closure of cracks/pores. The linear velocity/pressure correlation is due to intrinsic compressibility elastic moduli of the mineral skeleton of the rock. The velocity measurements during the depressurization run (filled symbols) show larger values than those obtained during the pressurization runs (empty symbols). This effect (hysteresis) is interpreted as due to the only partial re-opening of fractures and microcracks that close during pressurization. The intrinsic velocity V p (0) was calculated by the linear regression (red line) intersection of the high-pressure part of the second depressurization 2nd run back to the zero-confining pressure. These coefficients of the linear regression are corresponding to the pressure derivative of velocities. The fitting curve was calculated on the 2nd complete run using the Wepfer and Christensen (1991) equation  Table 3 presents the anisotropies calculated on V p0 Vp0 and V s0 Vs0, and the respective errors. The seismic velocities measured along the horizontal X-direction plugs are mostly higher than those measured along the vertical Z-direction (perpendicular to the bedding planes). Nonetheless , in some cases, the vertical seismic velocities are sporadically higher than the horizontal one. In general, the samples show very weak anisotropy. The only exception is anhydrite, in fact the maximum anisotropy was found in BX107 with 19.68% and 6.59% for the compressional-and shear-wave seismic anisotropy, respectively. The dolomite-anhydrite samples of the Muschelkalk section also show relatively high P-wave anisotropy (4.77 to 5.85%). Some limestones such as HU01 and HU07 show a relatively high P-wave anisotropy (5.14 and 5.39% respectively), while in general limestone have a lower average anisotropy for P-and S-wave. The Lettenkohle dolomite sample shows the lowest value of 0.26% and 0.51% for the compressional-and shear-wave anisotropy, respectively. Like seismic velocities, anisotropy tends towards linearity for confining pressures above 150-180 MPa. Figure 9 shows one example of anisotropy calculated from the fitting curves of the sample HU24. The anisotropy at low pressure is generally much higher than the highpressure zone, indicating that the cracks give a contribution to the anisotropy of the solid skeleton of the rocks.
In addition to anisotropy, we calculated the parameters ǫ , and γ as defined by Thomsen (1986) for a Vertical Transverse Isotropy (VTI) symmetry case. Thomsen parameters are often used in oil and gas industry literature to describe low anisotropy formations in a synthetic and simplified way. P-and S-wave moduli are represented in the general form c xx = ρV 2 , where c is the rock stiffness. The P-wave anisotropy parameter ( ǫ ) measures the difference of P-wave velocity between vertical and horizontal directions; Eq. (4). Similarly, the S-wave anisotropy parameter ( γ ) measures the difference in S-wave propagation in the vertical and horizontal directions; Eq. (5). The P-wave moduli for the minimum velocity direction perpendicular to the foliation  Wepfer and Christensen (1991). The plot is split into two areas, on the basis of the linearity of the curves. The low-pressure zone, where the curves are non-linear, represents the pressure interval where the seismic behavior is dominating by crack and pore closures (e.g., Birch 1960). It shows an abrupt decrease due to the closure of the aligned microcracks which reinforce the anisotropy induced by the lattice-preferred orientation (LPO) of anisotropic minerals (sample: HU24). Above 150 MPa, the curves become linear; the elastic behavior of the rocks is determined by the elastic moduli of the rock-forming minerals (Z) and maximum direction parallel to banding (X) plugs are represented by c 33 and c 11 , respectively. In order to calculate the c 11 component, when we used the maximum velocity, which in most of the cases is measured in the bedding plane. The bedding parallel and perpendicular S-wave moduli are represented by c 66 and c 44 , respectively.
The calculated Thomsen parameters are reported in Table 3 and plotted in Fig. 10b. Figure 10a shows the relation between anisotropy of compressional and shear waves. While sandstones tend to be more anisotropic in the shear component (with the exception sample HU25), the anhydrite (BX-107) and dolomites show a higher anisotropy of P waves (with exception of sample HU24). Limestones are scattered in the plot. The V p /V s ratio is reported in Table 3 and has been calculated from the V 0 of compressional and shear waves. Dolomites and limestones show the highest ratios and sandstones the lowest, with the lowest values for the Buntsandstein formation.

In-situ conditions
In order to calculate seismic parameters at in-situ conditions we calculated lithostatic and geothermal gradients for the Humilly borehole. Lithostatic pressure, as well as the geothermal gradient, were calculated for the depth interval of the borehole. Equation (6) was used to determine the derivative of the lithostatic pressure (MPa), and to calculate the in-situ pressure for each plug at depth The maximum lithostatic pressure calculated at the bottom of the borehole is 87 MPa. The geothermal gradient in the Geneva Basin, following Chelle-Michou et al. (2017), has an average value of 27.5 • C/km, therefore the calculated temperature at the bottom of the borehole is around 84 • C. Geothermal and lithostatic gradients are needed to calculate the seismic velocities at in situ conditions. Derivatives for pressure, measured from our experimental data were applied to calculated velocities at borehole conditions. Derivatives for temperature, available in the literature for limestone (Christensen 1985) and sandstone (Kern and Wenk 1990), are about 0.001 km s −1• C −1 , therefore we considered that the effect of temperature is negligible in the range of temperatures of our borehole, and no correction for temperature was applied.
Velocity data measured during pressurization up to 250 MPa and depressurization cycles, were interpolated (solid line in Fig. 8) using the four-parameter equation of Wepfer and Christensen (1991) (Eq. 7). V p,s is the seismic velocity (m/s), P confining pressure (MPa), and A (m/s), a, B (m/s), and b (1/MPa) all are the four adjustable parameters.
The equation was then used to calculate velocities at in situ conditions. The Wepfer and Christensen fitting  parameters are presented in Table 4, together with the P-and S-wave velocities calculated for the lithostatic pressure conditions corresponding to the depth of the samples. Table NA presents also the seismic impedance ( − Z ), i.e. the product of velocity and density, calculated for each sample at in situ conditions. Dolomite from the Muschelkalk formation shows the highest values of impedance, and the lowest in the limestones of the Upper Jurassic series. The impedance was then used to calculate the reflection coefficients as explained in the following section.

Elastic moduli calculation
The relationship between the seismic velocities compressional-and shear-wave was used to calculate the Poisson ratio (Eq. 8), where ver means that we used velocities in the Z-direction while for hor we used velocity values either X-or Y-directions. As for anisotropy, the dynamic Poisson ratio shows considerable variations as a function of the confining pressure, especially in the low confining pressures range (< 150 MPa) until all microcracks have been completely closed. Within the mineral compressibility region, the Poisson ratio stabilizes and have a gentle slope (Fig. 9).
The Poisson ratio, together with other elastic moduli that calculated along the horizontal and vertical directions for all the samples are as illustrated Fig. 11. All samples tend to concentrate along the 1:1 line, as expected in isotropic rocks. The only exception is sample BX-107 reflecting the anisotropic character of this sample. In Fig. 11

Influence of porosity on velocities
The observation on texture and mineral content at the plug-scale provides useful information for better understanding the seismic characteristics of the studied rock types. In Fig. 12 the average seismic velocities calculated at 40 MPa confining pressure are plotted against effective porosity (a) and average permeability (b). We chose 40 MPa as a reference pressure for those diagrams in order to consider the effect of partially open fractures at an average in situ pressure condition. The correlation between porosity and permeability is positive (Fig. 12b), as expected by a general common behavior. The only sample that significantly deviates from the regression line is HU07, a limestone from the middle Jurassic section, that shows a low permeability/ porosity ratio. It has to be noted that the permeability of the sample is strongly anisotropic, with k values measured within the bedding plane three orders of magnitude higher than the one measured perpendicular to the bedding. The effect of the extremely low permeability in the vertical direction reduces the average permeability and is responsible for the deviation of the sample from the general trend.
Regarding the link between velocities and porosity, even if a negative trend is recognizable, the scatter of the data is quite large, meaning that for rocks of similar porosity we can have very different P and S velocities. Looking at closer detail we notice that for limestones, the monomineralic Upper Jurassic samples show a linear negative correlation between porosity and both P-and S-waves velocities. The same is observed in the argillaceous limestones of Middle and Lower Jurassic age, but with a steeper negative correlation, indicating that the decrease of velocity with increasing porosity is enhanced by other factors, probably associated with the effects of mineralogical composition. On the contrary in the Muschelkalk Formation, the porosity/ The in-situ confining pressure; (MPa) was derived using Equation (6) in text. A, B, a, and b are the parameters of the curve-fit of Wepfer and Christensen (1991) [Eq. (7) in text]. The acoustic impedance at borehole conditions (Pa s/m) for compressional-( − Z p ) and shear-waves ( − Z s ) was calculated by the multiplication of seismic velocities and grain densities (Table 2) NA Not available a Lithostatic pressure for the sample BEX was considered as for a depth of 2630 m velocity trend is positive indicating that porosity has a little influence on seismic properties, and other factors are more relevant. For sandstones, again the correlation is quite complex: the velocities nicely follow the general decreasing trend, for compressional-wave, but not for shear-wave. In summary, the seismic velocities and elastic properties depend not only on porosity as it is evident within the large scatter-type of data for the velocity-porosity relationships. The constituent minerals have a major contribution for those variations, together with other factors that are only partially addressed in this work such as cementation, pore geometry, and grain size.

Influence of mineral contents on velocities
The modal composition obtained from QEMSCAN analysis was used to interpret the influence of the main mineral components on the seismic properties of the rocks. The average of the seismic velocities recalculated at room pressure conditions ( V p0 and V s0 ) are plotted against the main rock forming minerals of the studied samples and shown in Fig. 13. We choose to plot V p0 and V s0 against modal components fractions, because they represent the intrinsic seismic properties i.e. the elastic properties related to the mineral skeleton and not to the void space (pores, fractures).
The Permo-Carboniferous-Lower Triassic sandstones are primarily composed of quartz (39-86%). A variable amount of clay minerals (3-19% of illite and kaolinite), is present either along the grain boundaries, in the occluded pores (sample HU25, Fig. 6e) or as multiple thinly laminated layering of fine-grained bands (HU29, Fig. 6f ). Even small amounts of clay can significantly change the elastic properties of the sandstone (Castagna et al. 1985;Han et al. 1986;Klimentos 1991). In fact, the average V p of quartz (Table 5) Wepfer and Christensen (1991). A clear negative trend confirms the general observation that the seismic velocity decreases with increasing effective porosity. The significant dispersion in the velocity-porosity relationship, particularly for the shear velocity, can be related to other controlling factors such as the mineralogy, grain size, grain shape and pore geometry. Error bars are within symbols a b c d e Fig. 13 Average seismic velocity of compressional and shear-waves are plotted against the content of various rock-forming minerals. It shows that the compressional velocities decrease with increasing quartz and clay content, while for shear waves the correlation is opposite for quartz (decrease) than for clay (increase). Despite a certain dispersion of the data, that we interpret as due to factors such as alterations and diagenesis, we observe a positive correlation between velocities (P-and S-waves) and dolomite content. Less evident is the correlation of velocities with calcite content and with evaporitic minerals. Error bars are within symbols and references therein). Similarly, averaged V s are much lower in clay minerals than in quartz. Therefore, the presence of clay in sandstone has the effect to lower both velocities, with greater effect on shear component. The increase in V p and decrease in V s within Bundsandstein samples (HU25 and HU27) is attributed to the decreasing percentage of quartz, as shown in Fig. 13a. Siderite (not shown in Fig. 13) is another mineral phase present in the Permo-Carboniferous sandstones (samples HU27 and HU29) with concentrations similar to those in clay minerals. The elastic moduli of siderite are much higher than the moduli of clay and quartz, therefore the presence of siderite has the effect to increase the elastic moduli, and consequently V p and V s of the whole rock. Siderite and clay in the sandstones have then competing effects.
The effect of clay in reducing both V p and V s is also visible in the shale (HU13, 30% illite) and in the claystone (CL126; 87% illite + kaolinite); Fig. 13b. The anisotropy of the only shale in our collection is relatively low ( AV p = 4.16%, and AV s = 5.80%) compared with data in literature (e.g. in Valcke et al. 2006), probably due to the scarce iso-orientation of the grains and to the competing effect of quartz and calcite; Fig. 6a.
There is not a clear influence of calcite on the velocity of the limestones in our collection (Fig. 13c). The same is observed for the dolomite (Fig. 13d), even if a slight tendency towards increased velocities is observed especially within the sandstones. The positive correlation with velocity is more marked for anhydrites in the  Simmons (1965) Muschelkalk formation (Fig. 13e): the higher the anhydrite content, the higher V p and V s . The influence of mineral content for both V p and V s is not evident, V p anisotropy tends to be higher in carbonates (high calcite content) than in sandstones (high quartz content), but it clearly does not depend on whether the rock is clastic or carbonates but seems to rely on the nature of the heterogeneities and microtexture within the rock (i.e., clay banding, laminations, microcracks, etc.). The influence of mineral content on V p /V s ratio is summarized in Fig. 14a, where shear-vs. compressional-wave of our samples are plotted together with the seismic velocities values for single crystals of quartz, calcite, and dolomite, i.e. the main components of our rocks.
The sandstones fall either on the line of constant V p /V s ratio of quartz, or in an intermediate region between the lines of quartz and clay minerals, nicely following the progressive content of clay minerals from top to bottom of the section. The Muschelkalk samples fall in the narrow field between lines of constant V p /V s anhydrite and dolomite. The Dogger limestones fall together with the Upper Jurassic limestone, all below the line V p /V s for calcite. This result is interpreted for Upper Jurassic limestone as the effect of porosity, lowering more the shear-than the compressionalwave (Toksöz et al. 1976), and for Dogger limestones as the effect of quartz content.
A summary of single crystal elastic properties of the main minerals, averaged along all propagation directions, is reported in Table 5. as a reference for the considerations described below.

Seismic reflectivity
The seismic impedance ( − Z ) of an elastic medium controls the seismic reflection process because seismic energy can be reflected at interfaces only where there is a change in impedance across the interface. It is defined as the product of velocity and density. The Reflection Coefficient ( R c ) is the difference in seismic impedance over the sum of seismic impedance of two rock types; Eq. (9). R c will be then a positive number when V 2 > V 1 , and a negative number when V 2 < V 1 , which means equal or inverse polarity between incident and reflected waves, respectively. The larger the contrast in seismic impedance, the larger the amount of incident energy that is reflected (and the smaller the amount that is transmitted). The seismic impedance at the interface of different rock-types, calculated from intrinsic velocity as well as at its corresponding borehole conditions, is provided in Table 6. We tested R c of the interfaces that occur repeatedly in the sequence crosscut by the Humilly borehole. We could not retrieve samples from each interface (above and below interfaces), as big part of the sequence was not sampled, but we tried to sample each litho-type and combine them in order to simulate most of the interfaces that occur at Humilly. That is why we also inserted data from (9) Fig. 14 A cross-plot of average compressional-versus average shear-wave velocities with isolines for Vp/Vs ratio of different common minerals. The compressional and shear wave velocities for such minerals were after (1) Carmichael (1990), (2) Humbert and Plicque (1972) literature for clays and evaporates because they were not sampled in Humilly. For the seismic wave to be reflected at normal incidence angle, the interface between two different rock types must generate a significant impedance contrasts and have a minimum width equal to onequarter of seismic wavelength (Christensen 1989;Nanda 2016;Warner 1990). According to Sheriff and Geldart (1995), an R c = 0.1 represents a strong reflector, whereas good reflections originate from rock interfaces with R c ≥ 0.04. In Figure 15, the compressional-and shear-wave velocities of the samples at borehole conditions have been plotted against the grain density. In order to investigate the effect of the seismic anisotropy on the acoustic impedance, and therefore on reflection coefficient, the velocities measured along the three structural directions have been used. Isolines of acoustic impedance (Pa s/m 3 ) have been drawn on the plots and the reflection coefficient arising from adjacent impedance isolines is indicated.
Rock-types plotting at a distance greater than two isolines are considered as a potential reflector. Figure 15 shows that the dolomite-rich Muschelkalk plots in the fast V p,s and higher density corner (upper right) referring to high acoustic impedance compared to the halite sample of Keuper which plot in the lower left corner of the figure. The stratigraphic interface between the Muschelkalk and the overlying Keuper Halite, therefore, gives a strong reflection.
The measurements along the Z-direction were used to simulate the normal incidence of seismic waves encountering the sub-horizontal layering during a near-vertical incidence reflection survey. The strongest reflecting interfaces were found within the evaporitic facies of the Keuper (Lettenkohle) between the overlying Liassic carbonates and underlying Muschelkalk

Table 6 Seismic Impedance (Pa s/m) for compressional-waves propagating perpendicular to the bedding planes of the sedimentary formations (plugs in Z-direction)
Seismic Impedance is calculated as a function of intrinsic seismic velocities, − Z 0 , as well as at the Bottomhole Pressure, − Z BHP . The in-situ seismic velocity was derived following Wepfer and Christensen (1991) relationship (Eq. 7) in text and using the parameters A, B, a, and b provided in Table 4 . 16 Comparison between grain density, seismic velocities and impedance measured in the laboratory with corresponding physical parameters from borehole logging. Compressional velocities along horizontal (red square) and vertical (blue circle) plugs have been calculated at in-situ condition applying the velocity/pressure derivatives using empirical formula of Wepfer and Christensen (1991). The vertical acoustic impedance has been used to calculate the reflection coefficient to simulate the normal incidence of seismic waves. Error bars are within symbols. Samples of Core IV are after: a Müller and Briegel (1977), b Freund (1992), c Yan et al. (2016) (See figure on next page.) unit ( R c = 0.27) and between the Dogger and the Upper Jurassic carbonates ( R c = 0.09). Reflective interfaces could exist between Keuper and Sinemurian carbonates. Figure 16 shows a comparison between density, seismic velocities, and acoustic impedance measured in laboratory and the logging data. For grain density, the results of the laboratory measurements are significantly higher only in the Upper Jurassic carbonates, reflecting the role of the high porosity in that section. Velocities are only slightly higher for sample HU01. Slightly higher values have been measured also in the limestones of the lower Jurassic section both for density and velocities, a difference that could also arise from the used resonance frequency (1 MHz transducers in Laboratory, and 10 kHz for acoustic log).

Conclusions
We have presented new experimental data on seismic properties using a collection of cylindrical samples obtained from the Humilly-2 borehole (France), which reached a depth of 3051 m in the top of Permo-Carboniferous sediments. Experimental measurements of physical rock properties and mineralogical analysis were performed to explore the links between sedimentary facies and seismic characteristics and provide a tool in the interpretation of seismic field data in terms of geological formations.
Among the studied formations, grain densities (kg m −3 ) measured at ambient conditions show a wide range from 2200 to 2935. Seismic velocities (m s −1 ) at the atmospheric pressure are in the range 3955 to 6648 and 2426 to 3737 for P-and S-wave propagation modes, respectively. The seismic anisotropy seems to be quite low ( AV p,s = 0.3-6%) in the limestone of Jurassic age, Muschelkalk, Sandstone (Triassic and Permo-Carboniferous), and Liassic Shale, while anhydrite exhibits an unusually high anisotropy value ( AV p = 20%, AV s = 7%). The variations in seismic characteristics are due to mineral content, porosity and micro-textural features.
In the Upper Jurassic limestone, monomineralic and highly porous limestones, the seismic behaviors are linked to porosity, that varies in volume and size over several orders of magnitude. These rocks are almost isotropic, in a good agreement with the homogeneity of the samples (absence of bedding, compositional layering, and shape preferred orientation of grains) observed both at hand-specimen scale and in thin section. The various generations of calcite, either primary or diagenetic, presumably have no optical preferred orientation. Limestones with the same characteristics as the Upper Jurassic section might generate good reflectors if in contact with dolomitic/anhydrite rocks such as the Muschelkalk formation. Also, the stratigraphic contact between Upper Jurassic limestone and underlying Dogger would give rise to a fairly high reflection.
In the Dogger limestones, the compressional-and shear-wave velocities are on average slightly higher than in Upper Jurassic limestone carbonates, as an effect of the low porosity. The mineral composition of these carbonates is more complex than in the Upper Jurassic limestones: clay minerals, dolomite and quartz are present and might be the reason of a V p /V s ratio slightly lower than for calcite single crystals. Their impedance might give rise to a reflective contact with the dolomitic Muschelkalk or with a shale such as the Liassic shale unit.
The only shale in our collection, in the Liassic shale unit, shows one of the lowermost values for compressional-wave velocity, in agreement with the relatively high content of clay minerals. The low V p /V s ratio can be influenced by the presence of quartz. The anisotropy, due to iso-orientation of phyllosilicate, is probably lowered by contrasting effects of quartz grains with no optical preferred orientation. The shale might produce good reflectors if in contact with sandstones, or dolomitic carbonates, and also with some limestones. The Liassic shale unit would give rise to a good reflection at the stratigraphic contact with the underlying Liassic carbonates.
The highest values of the seismic velocity and density were measured in the Muschelkalk; therefore, those rocks might give rise to good reflectors if in contact with almost all the other lithotypes. That Muschelkalk anhydrite/dolomite interface produces high-energy seismic reflections has been known since decades. What's new in our work is that we've data from real rock samples that justify these observations from seismic logs. The high velocities might be related to the high content of dolomite, a mineral with relatively high velocity compared to other common components such as calcite or quartz. The V p /V s ratio is also similar to the one reported for dolomite single crystal. The Muschelkalk anisotropy is relatively high (> 5%), due to the small-scale compositional layering, visible at hand specimen and thin section scale.
The sandstones of the Lower Triassic and Permo-Carboniferous section present a quite wide range of V p and V s values, that can be correlated with the volume of dolomite and siderite in the rocks, contributing to increase the velocities. The V p /V s ratio is also different from the pure quartz value in samples with more complex modal composition. The highest anisotropy values have been reached on the only sample with high content of clay minerals, distributed in thin layers.
The comparison of lab data and log data show a relatively good agreement, with exception for density in those rocks with high porosity. The lab measurements, performed at dry conditions, slightly underestimate the log values for bulk density and overestimate it for grain density, an effect that can be clearly explained by the presence of filled pores in the in-situ measurements.
The reflection coefficients calculated for the stratigraphic sequence in Humilly show possible good reflectors at the top of the Muschelkalk and at the top of the Keuper (with the caveat that we used samples from different boreholes or literature due to absence of Keuper samples in Humilly borehole). A good reflector can also exist at the contact between the Buntsandstein and Muschelkalk units.