Mohamed K Salah

Mohammad Alqudah

Christian David

Petrophysical and Acoustic Assessment of Carbonate Rocks, Zahle Area, Central Lebanon

Keywords: Zahle Area, Central Lebanon, Transport Properties, Elastic Properties, Shallow Marine Carbonates. major challenges

Transport and acoustic measurements have been conducted on fifteen rock samples collected from central Lebanon to characterize these rock properties and to assess the impact of lithology, pore types, pore sizes, and textural parameters on the storage capacity and elastic characteristics. The spontaneous imbibition of selected samples was examined by detecting the moving capillary front across a sample. The coefficient of capillarity derived from the spontaneous imbibition is positively correlated with permeability. Seismic velocities in the studied rocks vary widely and are positively correlated with rock density. Matrix-supported rocks with a larger proportion of micropores have lower seismic wave velocities, whereas the grain-supported carbonates exhibit low to moderate porosity, and higher acoustic velocities. The elastic properties of the studied carbonates indicate a granular, rather than a crystalline, texture of the majority of the studied rocks. The widely-variable elastic behavior of basically carbonate rocks implies that seismic reflection profiles in carbonate sequences may contain seismic reflections that do not result from non-carbonate intercalations but which result mainly from variable porosity, pore types, and shapes in the carbonate rocks themselves. This study provides important relationships between the textural, petrophysical, and elastic properties of carbonate rocks which are important for the assessment and understanding of the seismic reflection sections in subsurface carbonate reservoirs.

Introduction

Carbonate rocks are important reservoirs hosting 50-60% of the hydrocarbon reserves around the world (e.g., [START_REF] Burchette | Carbonate rocks and petroleum reservoirs: a geological perspective fromthe industry[END_REF]. They are characterized by complex fabric and pore geometry which together represent rocks having a higher proportion of macropores; and thus more permeable and much more significant fluid flow and hydrocarbon recovery. In the giant oil fields of the Middle East, microporous carbonate reservoirs are not uncommon and represent a significant component of reservoir rocks in the region [START_REF] Nurmi | Carbonates, the inside story[END_REF][START_REF] Ehrenberg | An overview of reservoir quality in producing Cretaceous strata of the Middle East[END_REF]). These rocks were not adequately studied because, compared to granular limestones, they are less permeable, ineffectively contributing to fluid flow, and hinder the extraction of their hydrocarbon content [START_REF] Ehrenberg | An overview of reservoir quality in producing Cretaceous strata of the Middle East[END_REF][START_REF] Vanorio | Laboratory measurements of the acoustic and transport properties of carbonate rocks and their link with the amount of microcrystalline matrix[END_REF]. These rocks, although porous with sometimes porosities between 8% and 30%, are characterized generally by low (rarely moderate) permeability [START_REF] Kaldi | Diagenetic microporosity (chalky porosity), Middle Devonian Kee Scarp reef complex, Norman Wells, Northwest Territories, Canada[END_REF][START_REF] Moshier | Development of microporosity in amicritic limestone reservoir, Lower Cretaceous, Middle East[END_REF]; Deville de [START_REF] De Periere | Morphometry of micrite particles in cretaceous microporous limestones of the Middle East: Influence on reservoir properties[END_REF]. The intercrystalline micropores in these carbonate rocks occur usually in a lithified matrix of microcrystalline calcite which is known as micritie [START_REF] Folk | A review of grain-size parameters[END_REF]).

In addition to the value of limestones as common reservoir rocks allover the world, the carbonate rocks of Lebanon are good karst aquifers and are highly-significant engineering and construction materials (e.g., [START_REF] Nader | Petroleum prospects of Lebanon: some remarks from sedimentological and diagenetic studies of Jurassic carbonates[END_REF]. In this study, we are concerned with the analysis of 15 rock samples collected from the Paleogene strata exposed in Zahle, Central Lebanon to investigate their petrography, nannofossil content, as well as other thin-section derived parameters. In addition, we measure the basic fluid transport and elastic properties including porosity, permeability, density, electrical conductivity, capillarity, and the acoustic wave velocities; then interpret the observed measurements utilizing the results of the petrographic and scanning electron microscope (SEM) investigations. The study of the spontaneous imbibition of selected porous samples provides valuable information on the fluid flow across reservoir rocks. The present study provides also important clues on the factors affecting the propagation of elastic waves and the physical characteristics of microporous limestones, hence upgrading our current understanding on the interconnections between the petrophysical and acoustic properties as well as the microstructural features of carbonate rocks.

Geologic Setting and Stratigraphy

During the Paleogene, the Arabian Plate was located between the Tropic of Cancer and the Equator [START_REF] Beydoun | Arabian Plate Hydrocarbon, Geology and Potential -A Plate Tectonic Approach[END_REF][START_REF] Alsharhan | Tertiary of the Arabian Gulf: Sedimentology and Hydrocarbon Potential[END_REF][START_REF] Sharland | Arabian Plate sequence stratigraphy[END_REF]. During this time, the opening of both the Atlantic and the Indian Oceans has resumed, resulting in a continuous northward movement of Arabia toward Eurasia and rapidly closing the Neo-Tethys Ocean; an event that inherited a deep foreland basin near the plate margin [START_REF] Sharland | Arabian Plate sequence stratigraphy[END_REF]. Thus, the early Paleogene times of Lebanon are dominated by deep carbonate platform of a monotonous chalky, fine-grained character, which is partially rich in organic matter.

The lower and middle Eocene rocks in central Lebanon overly the thick chalky Paleocene successions, and are represented mainly by two facies; one is a chalky limestone, often with chert concretions which is common in basin deposits, while the other is a nummulitic limestone which is found on the basin margins. The chalky limestone facies is geographically limited due to the influence of the Syrian Arc event on the region [START_REF] Alqudah | Calcareous nannofossil, nummulite and ostracod assemblages from Paleocene to Miocene successions in the Bekaa Valley (Lebanon) and its paleogeographic implications[END_REF]. This rock sequence has a total thickness of 400m in Zahle (Fig. 1) where it is well-exposed [START_REF] Walley | The lithostratigraphy of Lebanon: a review[END_REF]. The encountered facies reflect clearly the influence of tectonics on the area during the Paleogene [START_REF] Walley | The lithostratigraphy of Lebanon: a review[END_REF][START_REF] Alqudah | Calcareous nannofossil, nummulite and ostracod assemblages from Paleocene to Miocene successions in the Bekaa Valley (Lebanon) and its paleogeographic implications[END_REF].

Tectonic events, mainly the collision between Arabia and Eurasia, affect the sedimentation pattern in the region causing a regional uplift, which resulted in the formation of the paleolake of Zahle [START_REF] Beydoun | Evolution and development of the Levant (Dead Sea Rift) Transform System: a historicalchronological review of a structural controversy[END_REF][START_REF] Sharland | Arabian Plate sequence stratigraphy[END_REF]. Neogene rocks of the Bekaa valley include deposits of non-marine origin containing nonmarine microfossils such as charophytes, gastropods and ostracods [START_REF] Sanjuan | Charophyte flora from the Miocene of Zahle (Beeka Valley, Lebanon): Biostratigraphic, palaeoenvironmental and palaeobiogeographical implications[END_REF]. The preservation and diversity of these microfossils imply that the paleolake of Zahle was well-established and characterized by a permanent, shallow oligotrophic freshwater.

Materials and Methods

The studied samples were collected from two closely-spaced exposures in Zahle; the first of which is the early and middle Eocene sequence where the chalky limestones enriched in planktonic foraminifera and calcareous nannofossils (Fig. 2, and Fig. 3A-C) prograde to the middle Eocene shallow marine platform rocks which are represented by a nummulitic bank (Fig. 3D). The second location is about one km away from the first one; the sediments of which belong mainly to the Miocene (Fig. 1). The exposed rocks in this location are composed essentially of mixtures of marl to marly limestones with some intercalations of ash and lignites (Fig. 3E andF).

Fifteen rock samples were collected from the Paleogene and Neogene deposits in Zahle. The number of samples that were obtained from the chalky limestones, the nummulitic limestones, and from the Neogene deposits, are ten, two, and three, respectively. Core samples were extracted from fresh outcrops using a portable, gasoline-powered, core-driller (Model D261-C). Thin sections were prepared from all collected samples, then examined under the polarized microscope to identify the standard microfacies, differentiate between the primary and secondary porosity, investigate the different post-depositional diagenetic processes, and to relate these features to the prevailing paleoenvironmental conditions. The procedures described in [START_REF] Salah | Effects of porosity and composition on seismic wave velocities and elastic moduli of Lower Cretaceous rocks, central Lebanon[END_REF], and references therein, were followed her to characterize porosity, its origin, grain types, cement and matrix nature, sorting, and any observed diagenetic process. We also examined the calcareous nannofossil content of the rock samples, which provides additional constraints on the paleoecology. The obtained results are listed in Table 1.

Field-drilled plugs were subsequently slivered in the lab to optimum lengths of 2 -5 cm long, cleaned with distilled water, and then dried at a temperature of 80 °C for 24 hours [START_REF] Mavko | The Rock Physics Handbook[END_REF]. Porosity (φ) was determined by the triple weight method using water as the saturating liquid. This method gives the effective porosity of a rock sample since it is only the connected pore spaces that are occupied by the saturating fluid.

The bulk density (ρb, g/cm 3 ) is calculated from the dry mass of the core sample divided by its bulk volume obtained during the porosity measurement, whereas the grain density (ρg, g/cm 3 ) is calculated from the porosity and bulk density using the relation: ρg = ρb/(1-φ).

Permeability is measured by applying Darcy's law using nitrogen gas and converted automatically to the liquid permeability (KL) as proposed by [START_REF] Klinkenberg | The permeability of porous media to liquids and gases[END_REF]. The electrical conductivity was measured twice; first after samples were saturated with tap water (water conductivity; σw = 0.66 mS/m) and then with saline water (σw = 58 mS/m). The formation resistivity factor (F) , at the two different brine concentrations, was calculated from the electrical conductivity of the rock sample (σr) and the water conductivity as: F = σw/σr. The obtained petrophysical and electrical conductivity data are listed in Table 2.

Many methods are available to study the capillary imbibition and to monitor the moving capillary front in rocks such as the standard weight increase measurement, X-ray imaging, and ultrasonic monitoring [START_REF] David | Detection of moving capillary front in porous rocks using X-ray and ultrasonic methods[END_REF].

In this study, we used the routine mass increase method due to the water intake with time to determine the coefficient of capillarity and examine its relationship with permeability and the rock micro-texture. This was accomplished through continuously weighing a sample while monitoring the height of the water front during imbibition [START_REF] Siebold | Capillary rise for thermodynamic characterization of solid particle surface[END_REF]. The resulting data are presented by plotting the water absorption as a function of the radical of time. In this way, a straight line would be obtained if the porous rock is replaced by one pipe of constant radius r. Because Earth materials are often inhomogeneous, this linear relationship is never discerned for normal rocks. However, the capillarity (C) parameter can be estimated from the average slope of the plot between the mass of water versus the square root of time plot: C = 1/πR 2 .dm/d(t) 0.5 , with R being the sample radius, m is the water mass intake, and t is time.

Primary and secondary wave velocities (Vp, and Vs) were first measured for the dried plugs at ambient conditions using Panametrics Pulser-Receiver (Model 5058PR) and an Agilent DSO-X-2014A Digital Storage Oscilloscope (100 MHz). Two transducers are attached firmly to the plug ends and the travel time of the acoustic pulse across the core is read on a screen with an accuracy of 0.01 µs, and then used, along with the core length, to calculate the corresponding seismic velocity. Three measurements of velocity at perpendicular directions have been taken for each core and the results averaged in order to account for any possible anisotropy and to have a more representative velocity for each rock sample (e.g., [START_REF] Soete | Acoustic properties in travertines and their relation to porosity and pore types[END_REF]. Acoustic velocities are strongly affected by fluids due to scattering, rock-fluid interaction, and local grain-scale flow effects [START_REF] Mavko | Estimating grain-scale fluid effects on velocity dispersion in rocks[END_REF]. Fluid saturation causes also velocity dispersion where the measured velocity increases with the measurement frequency, whereas the seismic velocity in dry rocks are insensitive to the frequency [START_REF] Mavko | The rock physics handbook[END_REF][START_REF] Khazanehdari | Variation in dynamic shear elastic modulus of sandstone upon fluid saturation and substitution[END_REF]. Therefore, both Vp and Vs have been re-measured for the watersaturated samples. The obtained bulk density and acoustic wave velocities were utilized to compute the acoustic impedance and the Poisson's (ν) ratio of the studied rocks. The ν ratio is the ratio between transverse contraction over longitudinal stretching of a stressed sample and controls the speed and reflection of seismic waves and, therefore, is substantial for anticipating the quality of rocks (e.g., [START_REF] Rafavich | The relationship between acoustic properties and the petrographic character of carbonate rocks[END_REF][START_REF] Soete | Acoustic properties in travertines and their relation to porosity and pore types[END_REF]. It is calculated from Vp and Vs as follows (e.g., [START_REF] Mavko | The rock physics handbook[END_REF]:

ν = (Vp 2 -2Vs 2 )/2(Vp 2 -Vs 2 )
In seismology, the ν ratio is used as an indicator of the physical properties and the fluid content in rocks, as well as the seismogenic behavior of the Earth's crust [START_REF] Zhao | Tomography of the source area of the 1995 Kobe earthquake: evidence for fluids at the hypocenter?[END_REF][START_REF] Zhao | The 1995 Kobe earthquake: seismic image of the source zone and its implications for the rupture nucleation[END_REF][START_REF] Kearey | An Introduction to Geophysical Exploration[END_REF][START_REF] Salah | 3-D seismic structure of Kii Peninsula in southwest Japan: evidence for slab dehydration in the forearc[END_REF][START_REF] Salah | Seismic velocity and Poisson´s ratio tomography of the crust beneath East Anatolia[END_REF]. The elastic moduli and the ν ratio relate stresses applied to rocks to resulting strains, and hence are good measures of the rock stiffness and mechanical properties [START_REF] Siegesmund | Physical and mechanical properties of rocks[END_REF][START_REF] Yu | Effects of porosity on seismic velocities, elastic moduli and Poisson's ratios of solid materials and rocks[END_REF]. The ν ratio and Young's modulus, in particular, are used to assess the mechanical behavior of rocks used in construction purposes [START_REF] Gerçek | Poisson's ratio values for rocks[END_REF][START_REF] Hudson | Faults in crystalline rock and the estimation of their mechanical properties at the Olkiluoto site, western Finland[END_REF]. The combinations of either the shear and bulk moduli or the ν ratio and the Young's modulus are enough to fully characterize the mechanical properties of Earth's materials. The measured and calculated petrophysical and elastic parameters are listed in Table 2.

Results and Discussion

Petrography

Based on the petrography of the collected rocks, the Eocene limestones of Zahle area can be divided into two facies; fine-grained limestones (samples E1-E10) varying from wack-to packstones, and coarse-grained limestones (samples E-11 and E12). The collected Miocene rocks (samples M1-M3) are mainly packstones with minor wackstone intercalations (Table 1; Fig. 2).

Wackstone

Samples E1, E9, and E10 are dominantly matrix-supported with about 10% bioclasts of planktonic and benthic foraminifera which are generally well-preserved (Fig. 4A andB; Fig. 5 D-I), in addition to some shark bone remains (Fig. 5G-H). Non-bioclasts in these samples are composed mainly of calcite; matrix is almost calcitic although it contains some amount of iron oxides. Calcium carbonate sparite either partially or completely fills the intragranular pores (Fig. 4A andB). In the same time, secondary calcium carbonate precipitations fill the spaces in the fractures (Fig. 4C).

Wack-/packstone

This microfacies is detected in samples E4 -E7 where the percentage of bioclasts increased slightly. Both well-and poorly-preserved planktonic and benthic foraminifera are detected; reflecting shallower marine settings during the deposition of these rocks (Fig. 4L, M, N, P, Q, R, andS). Non-bioclasts are represented mainly by calcium carbonate crystals of irregular shapes. All bioclasts and non-bioclasts are floating in a groundmass of micrite composed mainly of calcium carbonate as well as other components like iron oxides (Fig. 4L andM).

Porosity resulting from the intraparticle pores inside the foraminifera are rare in this case since they have been filled with sparite (Fig. 4L). Intergranular pores are also recognized (Fig. 4M). Although the secondary precipitation of sparite diminished the primary porosity in some cases; the pore space volume has increased in others due to the high dissolution of calcium carbonate producing sometimes very large pores (Fig. 4O).

Therefore, porosity varies from moderate to high according to the degree of dissolution.

Packstone

This microfacies was detected in the majority of samples and can be divided into three, slightly different, sub-facies: the first includes samples E2, E3 and E8; the second comprises samples E11 and E12; whereas the third one includes samples M1, M2, and M3.

Samples of the first subfacies contain both bio-and non-bioclasts with a good percentage of matrix and micro-sparite in their texture. Bioclasts contain benthic and planktic foraminifera, gastropods and bone fragments (Fig. 4D-K, and Fig. 5A-C). Non-bioclasts are represented by well-sorted micro-calcium carbonate crystals (Fig. 5 A-C). These micro-crystals resulted probably from secondary dissolution and re-precipitation of calcium carbonate as sparite crystals. Porosity is represented by intragranular spaces (Fig. 5C), with clear dissolution effects in the samples (Fig. 5A andB). The SEM imaging helped to check some processes that might have affected porosity. In this regard, two features have been observed which comprise the free formation of calcium carbonate crystals inside the foraminifera chambers as a secondary precipitation (Fig. 6-1), and the development of fenestral porosity (Fig. 6-5). These two processes have contrasting effects on porosity; the first diminishes porosity, while the second would enhance it.

The second subfacies is found in samples E11 and E12, in which the nummulite sp. is the major bioclastic component. The nummulites mark the margin between the continent and ocean setting pointing to a very shallow marine deposition. In such setting, the bioclasts, i.e nummulite sp., are dominant, while micrite is barely found (Fig. 5K, L, andM). Porosity is very high, and is represented by intergranular primary pores and fenestral secondary pores, which are filled partially with cementation followed by later dissolution (Fig. 5J andL).

The third type of packstone subfacies is detected in the Miocene lacustrine sediments; specifically in samples M1, M2, and M3, where a significant amount of gastropods and ostracods is encountered (Fig. 5N, andO). These bioclasts are surrounded by fine lime particles partially filling the interior of the bioclasts (Fig. 5P, andQ). The sediments are barely affected by cementation and other diagenetic processes where most of the bioclasts are not filled with precipitates of calcium carbonate although some of the grains suffered dissolution (Fig. 5Q). Some samples, however, show an evidence of fracturing which enhanced the porosity initially; later on, a secondary calcium carbonate precipitated in the fracture space (Fig. 5R). In the same way, a secondary calcium carbonate precipitated also inside the gastropod chambers; thus reducing the porosity of the rock sample (Fig. 5S).

Calcareous nannofossil assemblages and porosity

Calcareous nannofossils are more common in the pelagic settings rather than in the neritic zone. Due to this fact, calcareous nannofossil assemblages are common in the majority of the Eocene samples from E1 to E10, whereas samples E11 and E12 and those of the Miocene (M1 to M3) are barren. Marker calcareous nonnofossil species recognized in the samples such as T. serotinus and D. barbadiensis assign an early Eocene age for these sediments (Fig. 6-3, and8).

On a micro-scale, both intergranular and intragranular pores were observed (Figure 6-2, 6, 9 and 10), due to the good preservation of microfossils in the deep pelagic settings. Nannofossils, which are defined traditionally by their rims and central area, are commonly empty or have a sieve, cross, or bar shapes. In all aforementioned cases, micro-porosity still exists if the assemblages are quiet well-preserved or where there is no clay intrusion in the central area 7,9,and 11). In few cases, like in Discoaster spp., the central area is closed meaning that no microporosity can be afforded (Fig. 678, 10, and 12).

Transport and elastic properties

Measured values of φ, ρ b, ρg, the Klinkenberg permeability, F1 and F2, Vp, Vs, the ν ratio, and the calculated elastic moduli are all listed in Table 2. The measured porosity varies widely between from 0.01 to 0.40, with a mean value of 0.23. The lowest-porosity samples (E11 and E12) are grain-supported packstones (nummulitic) enriched in benthic foraminifera with around 20% micro-calcite crystals. Although the measured Klinkenberg permeability of the studied samples varies over three orders of magnitude (from 0.005 to 9.013 mD), it is generally low with an average of 0.94 mD. The two low-porosity samples have also the lowest permeability (0.005 and 0.007 mD, respectively). The low permeability originates mainly from the large proportion of micropores and the existence of micritic and micro-calcite crystals, which do not contribute effectively to the flow of fluids. The highest porosity was observed in the matrix-supported wackstone facies (Tables 1, and2).

The bulk and grain densities vary between 1.62 and 2.68, and 2.62 and 2.72, with averages of 2.07 and 2.69 g/cm 3 , respectively. The unified lithology of the studied rocks is reflected in the very narrow range of the grain density data, whereas the complex diagenetic processes affecting these samples induce wide variations in their porosity and accordingly in their bulk densities. The electrical conductivity data vary over more than three orders of magnitude for a given measurement cycle (Table 2). The calculated values of the formation resistivity factor (F1 and F2) are generally small because of the high conductivity of most rock samples induced by their generally moderate/high porosity.

Acoustic wave velocities in the dry samples vary widely from 2281 to 6449 m/s and from 1438 to 3292 m/s for Vp and Vs, respectively (Table 2). The average Vp and Vs are 3878 and 2213 m/s; whereas the average Vp/Vs and ν ratios are 1.71, and 0.23, respectively (Table 2). The ν ratio varies, however, between 0.16 and 0.33; whereas the former varies from 1.56 to 1.98 (Table 2). The acoustic impedance (AI) of a medium through which the seismic waves are propagating is the product of Vp and ρb. It controls the reflectivity of rocks and is an essential parameter in the synthesis of seismic sections and their correlation with the observed seismic profiles (e.g., [START_REF] Soete | Acoustic properties in travertines and their relation to porosity and pore types[END_REF]. It is affected by many parameters such as lithology, porosity, fluids, and pore geometry [START_REF] Zhang | Improving porosity-velocity relations using carbonate pore types[END_REF]. The Vp and ρb data of the dry and water-saturated samples listed in Table 2, result in AI values that vary, respectively, from 3.7 to 17.25*10 6 kg/m 2 .s, and from 3.38 to 17.3*10 6 kg/m 2 .s.

Discussion

Porosity versus other petrophysical properties

Porosity and bulk density are two closely-related parameters which are inversely correlated to each other (e.g., [START_REF] Ali | Depositional constraints and diagenetic pathways controlling petrophysics of Middle Miocene shallow-water carbonate reservoirs (Leitha limestones), Central Paratethys, Austria-Hungary[END_REF]. Their robust negative correlation (Fig. 7A) is a good indicator of the reliability of the measurements. On the other hand, porosity is directly proportional to permeability (the poro-perm relationship) with a coefficient of correlation (R 2 ) of 0.78 (Fig. 7B). The observed slight scatter indicates that, in addition to porosity, there are other factors affecting permeability in carbonate rocks including the grain size, pore volume, pore throat size, and pore geometry. Few samples have different permeabilities at a given porosity such as samples E2 and M1, which have different permeabilities of 0.02 and 0.572 mD, but similar porosities of 0.16 and 0.17, respectively (Table 2, Fig. 7B). Another fact is that although most studied samples are porous, with sometimes porosity of up to 0.40, permeability is generally very low (few fractions of a millidarcy); perhaps because of the small size and the poor connectivity of many pores, as well as the dominance of foraminifera with its cul-de-sac (pores closed at one end) porosity. Carbonate rocks have primarily complex pore fabrics and are prone to diagenetic alterations such as dissolution, recrystallization, cementation, and dolomitization, which additionally complicate the poro-perm relationship [START_REF] Nurmi | Carbonates, the inside story[END_REF][START_REF] Lapponi | From outcrop to 3D modelling: a case study of a dolomitized carbonate reservoir, Zagros Mountains, Iran[END_REF][START_REF] Dewit | Reservoir characteristics of fault-controlled hydrothermal dolomite bodies: Ramales Platform case study[END_REF].

Porosity is strongly correlated with the formation resistivity factor (F1 and F2 data) by power relationships and high correlation coefficients (Fig. 7C, andD). The greater improvement in the electrical conductivity of brine-saturated carbonate rocks (the second run of measurements) is expected if they have negligible clay content or any other electricity-conducting minerals [START_REF] Casteleyn | An integrated study of the petrophysical properties of carbonate rocks from the "Oolithe Blanche" formation in the Paris Basin[END_REF][START_REF] Regnet | Acoustic and reservoir properties of microporous carbonate rocks: Implication of micrite particle size and morphology[END_REF]. The obtained average values of the tortuosity (a) factor, in the equation: F = a.φ -m [START_REF] Archie | Electrical resistivity log as an aid in determining some reservoir characteristics[END_REF][START_REF] Guéguen | Introduction to the Physics of Rocks[END_REF] at the two measurement runs are 0.56, 1.83, while those of the cementation (m) exponent are 1.72 and 1.89, respectively. [START_REF] Regnet | Acoustic and reservoir properties of microporous carbonate rocks: Implication of micrite particle size and morphology[END_REF] found, in his log-log plot of porosity versus the formation resistivity factor, that the data are not clustered by texture because transmitting properties are greatly affected by the pore diversity even within the same setting. They, however, found a clear clustering by micrite microtextures revealing a strong relationship between micrite properties and micrite morphology. In our study, it is difficult to show such correlation because of the small number of samples.

Permeability and other rock properties

Porosity is inversely correlated with bulk density but directly correlated with permeability, respectively (Fig. 7A, andB). Therefore, an inverse relationship is expected between bulk density and permeability (Fig. 8A). Some scatter in the points is seen, and the correlation coefficient is fair (R 2 = 0.72). The difficulty to measure permeability in the field and the relative easiness of measuring the electrical properties have encouraged many researchers to study the relationship between permeability and the formation resistivity factor. The formation resistivity factor is inversely related to permeability usually with a power relationship (e.g., [START_REF] Ravalec | Transport properties and microstructural characteristics of a thermally cracked mylonite[END_REF][START_REF] Regnet | Acoustic and reservoir properties of microporous carbonate rocks: Implication of micrite particle size and morphology[END_REF]. According to [START_REF] Walsh | The effect of pressure on porosity and the transport properties of rock[END_REF], the two parameters are related as: K α F -p , where 1 ≤ p ≤ 3. The values obtained in our study are 0.75 and 0.76, which are very low, indicating the poor correlation between permeability and formation resistivity factor and that estimating permeability from electrical measurements is subject to large possible errors [START_REF] Regnet | Acoustic and reservoir properties of microporous carbonate rocks: Implication of micrite particle size and morphology[END_REF]. The improvement in the coefficient of correlation from 0.58 at the first run of brine concentration to 0.73 at the second run (Fig. 8B, andC) points to the effect of factors other than lithology, such as the type and geometry of the pore spaces as well as the fluid nature which dominate the electrical properties of the investigated samples at the higher brine concentration.

Permeability is inversely correlated with the acoustic wave velocities but with fair to moderate correlations (Fig. 8D, andE). Because permeability is sensitive to many other parameters rather than to the total volume of pores only, the relationships between permeability and acoustic wave velocities exhibit higher scatter and lower correlations compared to the φ-V relationship. The measured Vp and Vs for the dry and water-saturated samples are, however, strongly related to each other (Fig. 8F) because they are basically influenced equally by most rock parameters such as lithology, porosity, pore types and shapes although they may respond differently to the fluid type and degree of saturation. Both relationships have higher correlations indicating the consistency of the measurements. For a given dry rock sample, the Vs is 57% of the corresponding Vp value, but decreased to only 48% of the corresponding Vp after saturation (Table 2; Fig. 8F). This is induced mainly by the stronger reduction of shear wave velocities under saturation conditions [START_REF] Anselmetti | Acoustic properties of Neogene carbonates and siliciclastics from the subsurface of the Florida Keys: implications for seismic reflectivity[END_REF][START_REF] Kearey | An Introduction to Geophysical Exploration[END_REF][START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF].

Capillary imbibition

It is not easy to achieve complete saturation when conducting capillary imbibition tests because the driving mechanism of the capillary front in substances with pores of variable sizes induces air-trapping. Thus, estimating the exact amount of saturation (by water) in the middle phases is not that easy. In contrast, it is much easier to obtain the overall saturation when the front of water reaches the topmost end of the investigated core (David et al. 2015). As stated before, we measured the mass increase (due to water intake) with time for eight samples, then plot the mass intake versus radical time data to calculate the coefficient of capillarity C. The imbibition curves of five plugs are illustrated in Fig. 9A, revealing the satisfactory fit between the mass intake and the square root of time. The obtained coefficient of capillarity C, the strength of the fit judged from the value of the coefficient of correlation, porosity and permeability of the tested samples are listed in Table 3. The final saturation at the end of the experiment is also reported. In their study of the X-ray imaging of water motion during capillary imbibition in porous media, [START_REF] David | X-ray imaging of water motion during capillary imbibition in porous rocks. Part 1: methods, microstructures and mechanical deformation[END_REF] predicted a positive relationship between capillarity coefficient (C) and permeability (K) because both of them are controlled by the motion of fluids in the connected pore network of the material. The obtained relationship is illustrated in Figure 9B where a direct trend is observed between the two parameters which are related to each other as: K = 165.97C 1.6 (R 2 =0.82). The gradually rising water front is monitored during the experiment as shown for sample 9 in Figure 9C-F. Recently, [START_REF] Benavente | Predicting water permeability in sedimentary rocks from capillary imbibition and pore structures[END_REF] have arrived to an analogous relationship, which uses both the capillary coefficient and the rock microstructural parameters to calculate permeability as:

K = C 2 .ηr/4φρ 2 γcosθ
This equation includes information on water flow rate (the coefficient of capillarity C) water properties (viscosity η, and density ρ), wetting condition (interfacial tension γ, and the angle of contact θ) and pore structure (porosity φ, and radius 'r' of the pore). If these parameters are accurately known, a better correlation between permeability and the capillarity coefficient could be obtained. The scatter in our plot (Fig. 9B), and the moderate correlation coefficient are indeed induced by the variability of these parameters especially porosity and the pore radii. This possibility of predicting important rock parameters from simple experiments justifies the new interest in introducing the inexpensive capillary imbibition tests in the lab [START_REF] David | Detection of moving capillary front in porous rocks using X-ray and ultrasonic methods[END_REF].

Acoustic velocities and porosity

The relationship between Vp and the Vp/Vs ratio is shown in Fig. 10A. The obtained Vp/Vs and the ν ratios (Table 2) are in general lower than the commonly observed values for limestones (ν =0.3;Schlumberger, 2013) and may result from the existence of crystalline, rather than granular, textures in the studied carbonate rocks [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF]. This inference is supported by the petrographic study which indicates the dominance of crystalline textures in the studied samples (Figures 4, and5). Upon saturation, a great decrease in the shear wave velocity was observed, where the average Vp and Vs are 3723 and 1801 m/s (Fig. 8F); and averages Vp/Vs and ν ratios of 2.10 and 0.35, respectively. In this case, the Vs value is only 48% of the Vp value of the water-saturated sample. The Vp/Vs ratios of the water-saturated samples vary widely between 1.70 to 2.57 (Fig. 10B) compared to only from 1.56 to 1.98 for the dry samples. The higher Vp/Vs and ν ratios of the water-saturated rocks are induced mainly by the greater drop in the shear wave velocity. Although the Vp of the water-saturated samples are generally higher than those of the dry rocks, some previous results revealed a mixed behavior where by the Vp may increase or decrease with saturation (e.g., [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF]. Such behavior was also observed for our samples (Table 2).

Acoustic velocities are strongly dependent on bulk densities [START_REF] Kearey | An Introduction to Geophysical Exploration[END_REF] where velocities generally display linear increase with increasing density and decreasing porosity (e.g., [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF][START_REF] Ali | Depositional constraints and diagenetic pathways controlling petrophysics of Middle Miocene shallow-water carbonate reservoirs (Leitha limestones), Central Paratethys, Austria-Hungary[END_REF]; a fact that is clearly supported by the strong positive correlation between ρ b and both the Vp and Vs of the presently-investigated samples (Fig. 10C, andD). On the other hand, the AI, calculated for the dry and watersaturated rocks, is inversely related to porosity with strong correlations (Fig. 10E, andF).

Total porosity plays the major role in controlling the seismic wave velocity in rocks of similar mineralogical composition where both Vp and Vs decrease with increasing porosity (e.g., [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF][START_REF] Regnet | Acoustic and reservoir properties of microporous carbonate rocks: Implication of micrite particle size and morphology[END_REF]. Similar to the usually-observed concave-upward trends between porosity and acoustic velocity [START_REF] Neto | Evaluation of carbonate pore system under texture control for prediction of microporosity aspect ratio and shear wave velocity[END_REF], our seismic velocity data display inverse polynomial trends with porosity (Fig. 11A, andB). Velocity variation at a given porosity is higher for the water-saturated samples than for the dry rocks. The good fit of the V-φ relationship and the low/moderate dispersion of velocity data at a given porosity have also been observed in previous studies (e.g., [START_REF] Baechle | Effects of microporosity on sonic velocity in carbonate rocks[END_REF][START_REF] Fournier | Elastic properties of microporous cemented grainstones[END_REF], although the exact pattern may differ for a larger number of samples. Few samples, especially the water-saturated, display different seismic velocities at a given porosity. For example, samples E2 and M1 have similar porosities of 0.162 and 0.169 but differ largely in their Vp which are 3903 and 4654 m/s, respectively (Table 2, Fig. 11A) even though both are grain-supported packstones (Table 1). This behavior has also been observed by other researchers and arises possibly due to different pore types in the carbonate rock [START_REF] Eberli | Factors controlling elastic properties in carbonate sediments and rocks[END_REF]. Pores occurring in a rigid rock framework induce higher velocities than interparticle pores, for example. Tightly-cemented marine carbonate grainstones with significant amounts of micropores [START_REF] Fournier | Elastic properties of microporous cemented grainstones[END_REF][START_REF] Fournier | Pore space evolution and elastic properties of platform carbonates (Urgonian limestone, BarremianeAptian, SE France)[END_REF] show lower seismic wave velocities at a given porosity compared to travertines [START_REF] Soete | Acoustic properties in travertines and their relation to porosity and pore types[END_REF], which are deposited in thin water films under greatly variable sub-aerial conditions . Previous studies revealed that not only porosity, but also pore types and shapes, affect the seismic velocity in carbonate rocks and any increase in the proportion of micropores reduces the velocity at a given porosity. On the other hand, stiff (such as mouldic) pores would induce a more rigid framework and thus higher seismic velocity at a given porosity [START_REF] Anselmetti | Acoustic properties of Neogene carbonates and siliciclastics from the subsurface of the Florida Keys: implications for seismic reflectivity[END_REF][START_REF] Eberli | Factors controlling elastic properties in carbonate sediments and rocks[END_REF][START_REF] Weger | Quantification of pore structure and its effect on sonic velocity and permeability in carbonates[END_REF][START_REF] Xu | Modeling elastic properties in carbonate rocks[END_REF].

For the dry samples, the calculated ν ratio is inversely related to porosity with a high coefficient of correlation (Fig. 11C). [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF] noticed that the highest porosity rocks have Poisson's ratios lower than 0.28, which is consistent with our observations. [START_REF] Neto | Evaluation of carbonate pore system under texture control for prediction of microporosity aspect ratio and shear wave velocity[END_REF] found also an inverse relationship between porosity and the Vp/Vs ratio of a suite of carbonate rocks from Brazil. Samples showing this behavior are characterized by high-porosity recrystallized textures with a frame that is supporting the shear wave propagation. The ν ratio of the water-saturated samples increases significantly (Table 2, and Fig. 11C) because of the large decrease in Vs with saturation especially at higher porosities. A similar behavior has also been observed for the Neogene carbonates and siliciclastics from the subsurface of the Florida Keys by [START_REF] Anselmetti | Acoustic properties of Neogene carbonates and siliciclastics from the subsurface of the Florida Keys: implications for seismic reflectivity[END_REF]; hence the higher Vs sensitivity to the induced shear weakening.

Porosity is inversely related to the elastic moduli of the dry and water-saturated samples (Fig. 11D, E, andF). They are related to each other by polynomial trends that are clearly concave upward. In addition, there is a clear decrease in the values of Young's and shear moduli of the water-saturated samples (Fig. 11D, andF), although the Gassmann's theory assumes that the shear modulus remains essentially the same during fluid substitution. Results obtained in this study as well as other previous studies (e.g., [START_REF] Adam | Gassmann's fluid substitution and shear modulus variability in carbonates at laboratory seismic and ultrasonic frequencies[END_REF] indicate that this does not apply for limestones. The rocks are usually influenced by rigidity changes (up or down); hence restricting the application of Gassmann's theory to predict the seismic wave speeds in carbonate rocks. The shear modulus decreases with saturation in granular carbonate rocks, but increases in crystalline and cemented carbonates [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF]. Therefore, the elastic properties of a carbonate rock alone may enable a quick estimation of its granular or crystalline texture and can help in predicting the rock nature from seismic data acquired across boreholes. On the other hand, values of the bulk moduli have slightly increased upon saturation (Fig. 11E).

In addition to porosity, acoustic velocities in carbonate rocks are controlled also both by the depositional lithology and the subsequent diagenetic alterations [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF]. In this way, primary fabrics may be quickly overprinted by diagenetic changes such as recrystallization and cementation. Most of the matrixsupported wackstones have higher porosities and lower seismic wave velocities. On the other hand, samples E5, E11, E12, M1, and M3 have higher-than-average seismic velocities and are either matrix-supported wackstone/packstone or grain supported packstones. [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF] used the texture class assignment proposed by [START_REF] Kenter | Acoustic behavior of sedimentary rocks: Geologic properties versus Poisson's ratios[END_REF], whether granular or crystalline, to relate particular ranges of diagnostic acoustic velocities for each texture. Granular textures refer to depositional grains where grain-to-grain contacts support the rock framework. Rocks with granular textures have little matrix or cement and result typically in low Vp, relatively low Vs, and higher ν ratios under saturation conditions. Rocks with crystalline textures have been recrystallized or cemented after deposition and have altered their original grain-dominated fabric to a crystaldominated one. These rocks have generally low to high Vp, high Vs, resulting in low ν ratios. Highly-cemented, sparite-forming and intensive secondary crystal-fill samples of low porosity and higher seismic wave velocity such as E5, E11, E12, M1, and M3 belong to this category (Tables 1 and2).

Mineralogy and the ν ratio in carbonate rocks

The mineral composition has a larger impact on velocity variations in siliciclastic sediments, which exhibit wide variations in mineralogy from quartz to clay, than in carbonates [START_REF] Christensen | Seismic properties and the origin of reflectivity from a classic Paleozoic sedimentary sequence, valley and ridge province, southern Appalachians[END_REF][START_REF] Braaksma | Case history controls on acoustic properties of Upper Jurassic siliciclastic rocks (Boulonnais, northern France)[END_REF]. Consistent with earlier findings of [START_REF] Anselmetti | Controls on sonic velocity in carbonates[END_REF]; [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF] found no significant effect of the dolomite fraction on the measured seismic velocities of carbonates which indicate that the observed variations in seismic velocities are induced mainly by the rock texture rather than variations in mineralogy. Based on the results of [START_REF] Anselmetti | Controls on sonic velocity in carbonates[END_REF], 2001[START_REF] Kenter | Acoustic properties of "Young" carbonate rocks, Odp Leg 166 and Boreholes Clino and Unda, Western Great Bahama Bank[END_REF] porosity is the main control on the seismic velocity of carbonate rocks, although pore type and the percentage of micropores play also a role and this is why some rocks have different velocities at a given porosity. However, the details of the coupling between pore types and seismic velocity are not straightforward and very difficult to quantify, in contrast to those of the φ-V relationship because pore type characterization is a qualitative process and because many pores are combinations of different types and their description remains subject to interpretation (Lønøy 2006;[START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF]. Increasing quartz content reduces the acoustic velocities of the carbonate rocks because it resists dissolution, hindering cementation and the associated increase of the seismic wave velocities. On the other hand, acoustic wave velocities increase slightly with increasing dolomite content. [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF] evaluated the specific effects of rock texture on a crossplot of Vp and the ν ratio.

Samples which plot in the lower left of the plot (with Vp < 4500 m/s and ν ratio < 0.28) are classified as granular texture, whereas samples located on the remaining space of the plot (Vp > 4500 m/s and ν ratio > 0.28 or Vp < 4500 m/s and ν ratio > 0.28) have crystalline texture. In this way, we can distinguish the samples of granular texture (E1, E2, E3, E4, E6-E10, and M2) from those of crystalline textures (E5, E11, E12, M1, and M3) on the plot of Vp versus the ν ratio (Fig. 12A). Moreover, the crossplot of Vp versus ν ratio for the saturated samples displays a V-shaped cluster in which the granular textures occupy the upper part of the left wing, whereas the crystalline textures occupy the lower left and the right wing of the crossplot (Fig. 12B). This means that rocks of granular textures show large increase of ν ratio upon saturation with a mean value of 0.15, whereas rocks of the crystalline textures show only slight ν ratio increase [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF]. The transition between the granular and crystalline textures occurs at a Vp of 4500 m/s and a ν ratio of 0.28. Early diagenesis and compaction of rocks having granular textures increase their Vp and decrease their ν ratio (e.g., [START_REF] Wilkens | The ratio Vp/Vs as a discriminant of composition for siliceous limestones[END_REF]. With further maturation and cementation, the Vp will increase accompanied by a slight increase of ν ratio. Toward zero porosity (high cementation and high Vp), the ν ratio of carbonate rocks is approximated by the ν ratio of pure calcite and dolomite at 0.31 and 0.29, respectively [START_REF] Simmons | Single crystal elastic constants and calculated aggregate properties. A handbook[END_REF]. Therefore, samples E11 and E12 (grain-supported packstone) have very low porosity (1.1%), high Vp (above 6300 m/s) and high ν ratios of 0.32. The processes of dolomitization, recrystallization, dissolution, and cementation have differential impacts on the Vp and Vs, and, accordingly, on the ν ratio as explained in detail by [START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF]. Any diagenetic process may act on a rock whose Vp and ν ratio have a certain position on the plot but each process has its specific trend on the path [START_REF] Anselmetti | Subsurface geology of a prograding carbonate platform margin, Great Bahama Bank, Results of the Bahamas drilling project[END_REF].

In order to compare the presently-obtained trend between porosity and velocity with previous results, we plot our data superimposed on the trends obtained by the Raymer-Hunt-Gardner (RHG) model [START_REF] Raymer | An improved sonic transit time-to-porosity transform[END_REF], and the Wyllie's time average (WTA) equation [START_REF] Wyllie | Elastic wave velocities in heterogeneous and porous media[END_REF][START_REF] Wyllie | An experimental investigation of factors affecting elastic wave velocities in porous media[END_REF] in Figure 12C. In addition, the density-velocity trends of [START_REF] Gardner | Formation velocity and density-the diagnostic basics for stratigraphic traps[END_REF] (GRD) and [START_REF] Anselmetti | Controls on sonic velocity in carbonates[END_REF], (ANEB) equations are also superimposed on the crossplot of bulk density versus velocity (Fig. 12D). At a given porosity, velocities obtained in our study are lower than those of the RHG model but close to the results of the WTA equation. On the other hand, the obtained velocities are higher than velocities of the GRD model but closer the results of ANEB at a given bulk density. In general, depositional conditions, composition, the sequence of diagenetic alterations, and the subsequent deformation are all factors leading to the observed differences (Anselmetti andEberli 1993, 2001;[START_REF] Eberli | Factors controlling elastic properties in carbonate sediments and rocks[END_REF][START_REF] Verwer | Acoustic properties of carbonates: Effects of rock texture and implications for fluid substitution[END_REF]).

Conclusions

Transport and elastic properties (porosity, density, permeability, electrical conductivity, and acoustic velocities) coupled with thin sections, calcareous nannofossil and SEM investigations were conducted on fifteen carbonate samples collected from the Paleogene strata exposed at Zahle, central Lebanon. The electrical conductivity was in two successive runs of different brine concentrations, whereas the acoustic wave velocities were measured for the dry and water-saturated samples. We also examined the self-imbibition of selected samples by monitoring the moving capillary front across a sample. From the measured properties, we further calculated other petrophysical and elastic parameters such as the formation resistivity factor, the cementation exponent, the Vp/Vs and ν ratios, the acoustic impedance, as well as the rocks' elastic moduli. This integrated analysis provides a deep insight into the different parameters controlling the petrophysical and elastic behavior of the investigated rocks.

Porosity is negatively and positively correlated with bulk density and permeability, respectively. Imbibition experiments showed good linear fits between the mass increase and the square root of time. The resulting capillary coefficients exhibit good correlation with permeability. Although porosity has the major control on both the transport (permeability and electrical conductivity) and elastic properties of the studied rocks, the amount of microporosity plays a role as well. Acoustic velocities decrease with the increase in porosity because the ultrasonic wave propagation is impeded in porous media. On average, the shear wave velocity is 57% and 48% of the primary wave velocity of a sample under dry and water-saturation condition, respectively, revealing the great drop in the values of the shear modulus of rocks with saturation and the consequent increase of the Vp/Vs and ν ratios. Some samples in this study display different seismic wave velocities at equal porosities indicating that in addition to porosity, pore types, sizes and shapes as well as the rock diagenetic history control also the acoustic wave velocities. Higher velocities may result not only from lower porosity but also from the development of stiff and cohesive micrite microstructures. In addition, velocity increase may also result from increase of grain size and growing contact cementation. Grain size, packing, matrix, cement, and the crystal shapes are therefore additional factors controlling the propagation velocity of seismic waves in carbonate rocks. Acoustic impedance data of both the dry and water-saturated samples vary widely because of the wide variations in the dry ρb and the Vp. Therefore, seismic reflection sections in carbonate sequences may contain reflections produced from different pore types and shapes rather than the existence of intercalations of other non-carbonate minerals. The present study revealed also that seismic measurements such as the changes in the shear modulus, and hence the Vp/Vs and ν ratios with saturation are useful tools for discriminating between the different textures of carbonate rocks. The ability to extract the rock texture from seismic measurements allows the prediction of porosity types and its evolution from sonic logs acquired through similar sediments.

Results obtained in this study give a deeper insight into the understanding and interpretation of the various petrophysical properties of microporous rocks which are valuable in seismic-based exploration and evaluation of marine platform carbonate rocks. However, more detailed petrographic and pore fabric analyses, coupled with petrophysical and acoustic measurements on a large number of dry and water-saturated core samples are required for deeper understanding of the nature of these rocks and for extending the obtained results to the interpretation of the seismic profiles acquired above the corresponding subsurface occurrences of the studied formations.

Figure Captions

Fig. 1 .

 1 Fig. 1. Geological map of Lebanon in the Eastern Mediterranean region (A), with an inset (B) representing the details of the studied rock exposures in Zahle area which are shown by the red rectangle in (A).

Fig. 2 .

 2 Fig. 2. Columnar section of the Paleogene succession exposed in Zahle area.

Fig. 3 .

 3 Fig. 3. Field photographs representing the exposed strata in the studied area: A, B, and C belong to the early Eocene rocks, D represents middle Eocene outcrops, whereas E and F represent the Neogene Zahle Paleolake deposits.

Fig. 4 .

 4 Fig. 4. A plate showing the different types of microfacies encountered in the studied samples of the early Eocene with their porosity content.

Fig. 5 .

 5 Fig. 5. A plate showing the different types of microfacies encountered in the studied samples of the late Eocene and early Miocene with their porosity content.

Fig. 6 .

 6 Fig. 6. A plate showing the different calcareous nannofossil assemblages found in the studied Paleocene and Eocene strata.

Fig. 7 .

 7 Fig. 7. Crossplots of porosity versus bulk density (A), permeability (B), and the formation resistivity factor (F1 and F2) calculated from electrical conductivity measurements at two different brine concentrations (C and D).

Fig. 8 .

 8 Fig. 8. Crossplots of bulk density versus permeability (A); formation resistivity factors versus permeability (B and C); permeability versus acoustic wave velocities (D and E) and between Vp and Vs of the dry (solid circles) and water-saturated (colored squares) samples (F).

Fig. 9 .

 9 Fig. 9. (A) Water mass intake evolution of five carbonate samples as a function of the square root of time; and (B) The relationship between the coefficient of capillarity and permeability of the studied eight samples. Height evolution of the water front across sample E9 at different time intervals (C-F). The experiment started at 09:19:31 am and ended at 13:31:51 (i.e., lasted for 4 hours, 12 minutes, and 20 seconds).

Fig. 10 .

 10 Fig. 10. Crossplots of Vp versus the Vp/Vs ratio of the dry (A) and water-saturated (B) samples; ρb versus Vp and Vs (C and D); and φ versus the acoustic impedance (AI) of the dry (E) and water-saturated (F) samples.

Fig. 11 .

 11 Fig. 11. Crossplots of φ versus Vp of the dry (solid circles) and water-saturated (colored squares) samples (A); φ versus Vs of the dry (solid circles) and water-saturated (colored squares) samples (B); φ versus ν ratio of the dry (solid circles) and water-saturated (colored squares) samples (C); φ versus Young's modulus of the dry (solid circles) and water-saturated (colored squares) samples (D); φ versus bulk modulus of the dry (solid circles) and water-saturated (colored squares) samples (E); and φ versus shear modulus of the dry (solid circles) and watersaturated (colored squares) samples (F).

Fig. 12 .

 12 Fig. 12. Crossplots of Vp versus the ν ratio of the dry (A), and water-saturated (B) samples. Plots of φ versus Vp (C) and ρb versus Vp (D) are shown with the models of WTA, RHG, GRD, and ANEB superimposed (see text for details).

  Figure 6

Table 1 .

 1 Summary of the main results of the petrographic investigation of the studied rocks.

Acknowledgments

This research has been partially financed through a URB (University Research Board) Grant from the American University of Beirut (Award No. 103009; Project No. 22759).