Thermo-Physical and Mechanical Properties of Al Hashimiyya Basaltic Rocks, Jordan ()
1. Introduction
Over the last four decades, several investigations on geothermal energy utilizations have been done in Jordan. Most of such studies were accompanied by Natural Resources Authority (NRA) incorporated with different public and private companies and institutions. They deduced that Jordan has wide spread low enthalpy geothermal resources. In Jordan, there are 108 hot springs discharging annually about 25 million cubic meters of thermal water into the Dead Se from lower Cretaceous Sandstone [1] . In addition, various thermal water wells aligned NE Jordan along NW-SE faults. Recorded water temperature ranges from 42˚C to 50˚C. In this region, thermal wells discharge water from the upper Cretaceous limestone and lower Cretaceous sandstone. Basaltic rocks, the scope of this study, represent the upper unit of aquifer system in NE Jordan. In some areas, geothermal gradient is considered to be high and reaches 50˚C/km. In central Jordan, the estimated maximum reservoir temperature predicted by geochemical studies is around 115˚C [2] .
Most of the published work gave an attention to the Dead Sea Rift area where more than 50 hot springs were evaluated and considered [3] [4] . Thermal water in north and middle Jordan has been used directly as therapeutic water, e.g. Zara, Zarqa-Ma’in, Afra and North Shunah hot springs of an annual use about 1540 GWh [5] [6] . Other utilizations including bathing and irrigation had earlier been used. It is expected in the near future that shallow geothermal water will have greenhouses heating for flowers planting utilizations. In addition, fish farming is another future application to provide the local market with fresh fish [2] . In addition, [7] proposed a future geothermal utilization for air conditioning and heating of the Queen Alia Airport. On the other hand, [8] deduced that geothermal energy in Jordan could be utilized for electricity generation. Various locations in Jordan are suitable for shallow geothermal storage due to the recorded underground temperature within the upper 100 m in subsurface. [9] investigated geothermal cooling potential of basaltic reservoir in NE Jordan. The previous studies concluded that the origin of thermal waters in Jordan is of meteoric type, while the source of heat is the deep circulation of cold and hot mixed water. In NE Jordan, geothermal potential for different utilizations is expected. Geothermal resources in this area did not receive detailed investigations. Future exploration in NE Jordan needs further evaluations and assessment of this resource.
Correlation between reservoir rocks thermal conductivity and their physical properties had previously been modeled. Numerous models have been developed to determine thermal conductivity based on mineral constituents, density, permeability, porosity, uniaxial compressive strength and P-wave velocity (e.g. [10] - [16] ). Nevertheless, some investigations include oceanic basalt thermo-physical properties [17] have been done, such results could not be applied for continental flood basalts in Jordan. Physical properties are highly affecting the heat efficiency in any geothermal reservoir evaluation. Permeability is one of the most important properties which influence thermal conductivity. [18] stated that thermal conductivity and permeability correlation is only feasible where both properties are measured for the same sample; consequently, anisotropic factors should be considered. Comprehensive investigations were reported by [19] for thermal conductivity and permeability interrelation. He specified different factors controlling thermal conductivity and permeability interrelation such like; mineral geometry, crystal size, vesicularity includes vesicles internal geometry and microstructures. [19] investigated mineral crystal size and microstructure for their effect on thermal conductivity and permeability correlation.
[14] presented the data which expected to improve the statistical confidence on geophysical and thermo-hydro-mechanical numeric models data. In addition, he examined the ability to predict rock properties at greater depths of the Taupo Volcanic Zone in New Zealand. Mechanical properties correlation with thermo-physical properties were modeled, too. [20] found that the ratio of secondary minerals to primary minerals (SEC/PR) of the studied ultramafic rocks have a good correlations with their physicochemical and mechanical properties. They suggested that alteration has a negative effect on the engineering performance of the ultramafic rocks. [21] determined thermal conductivity and thermal diffusivity with empirical relations from physical properties of limestone, dolomite, and siliciclastic rock samples. He found that the Grüneisen parameter and melting temperature decrease nonlinearly but the Debye temperature decreases linearly with increasing porosity. In addition, [22] studied the petrophysical properties (density, P- and S-wave velocity, porosity) and ultrasonic P- and S-wave attenuation mechanism for the saturated Deccan basalt cores. They found an interesting relationship between P-wave velocity and porosity.
This study aimed at examines the rock physical properties influence on thermal conductivity. On the other hand, a relationship between examined thermo-physical and mechanical properties will be correlated. This will develop a better understanding of basalts’ thermo-physical and mechanical properties interrelations. The properties involved in this work are permeability, porosity, density, compressional wave velocity, point load, indirect tensile strength and unconfined compressive strength parallel with thermal conductivity. The results obtained in this work are expected to develop data confidence on thermo-physical models for the studied rocks. Therefore, it will improve the ability to predict basaltic rocks properties from limited input data at greater depths.
2. Geologic Settings
Successive basaltic flows were erupted from Cenozoic; Miocene to Pliocene time in the magmatic activity within Arabian plate [23] . This basaltic flood which covers an area of about 12,000 km2 is known as Harrat Al Shaam [24] . The Basaltic rocks in Al Hashimiyya are of Abed Olivine Basalt of Pliocene age [23] . The studied flow is considered as the western extinction of Harrat eruption (Figure 1).
[25] divided Jordanian basalts into six major flows (named B1-B6) and one eruption of tephra (given a name of B’t). B1-B3 basaltic flows are not visible in Jordan, but are known from borehole data [26] . [27] proposed a new classification based on K-Ar dating. They divided volcanic rocks into three major episodes: Oligocene to early Miocene (26 - 22 Ma), middle to late Miocene (13-8 Ma), and late Miocene to Pleistocene (7 Ma to < 0.1 Ma) Al Hashimiyya basalts belong to youngest eruption phase with an age of ca. 2 Ma [27] . Basalts of ca. 400 m thickness of successive flows are found in the NE part of the study area, while less than 50 m is found in the southern parts of the study area [26] .
The studied basalts show slight vertical and lateral variations in morphology. Basaltic flows in Al Hashimiyya are composed of a series of gently dipping to sub
![]()
Figure 1. Geologic map of the study area (after [23] ). Lower right; location map of the study area, Jordanian Harrat is modified after [24] . Lower left; lower 10 m of the studied basalt flow (H4 to H6 sub-flows).
horizontal lava flows. This basalt flow come from Al-Khalidiyya town, locates at south-eastern of Al-Mafraq city, then it extends to fill Wadi Ad Dhulayl through Manshiatt Al-Ulayyan to the west. It crossing Al-Zarqa-Al-Mafraq highway after that it’s partitioning into two flows, the northern one ends in Ar Rusaifa ash Shamaliya crossed Al Hijaz Railway. The second flow to the south is passing Al-Kherba as Samra ending south west at Al Hashimiyyia city down (Figure 1). At Wadi Al-Sukhneh the eruption divided into two parts, the highest one shows a step-like structure in Wadi Wad’a. While the lower extends along Wadi Az Zarqa up to city of Jerash (Sarrut village) [28] . Al Hashimiyya studied flow is measured to be around 20 m thickness. The basalts in the study area are lithologically well differentiated, and structurally developed. It is easy to identify each basalt flow from its special lithology and primary magmatic structures. The study area basalt cover about 30 km2 northern Zarqa city. The studied flow locates at latitude of 32˚08'50"N and longitude of 36˚05'15"E. Al Hashimiyya flow subdivided into six sub-flows dependence on lithological and structural aspects. The vertical extinction of each sub-flow eruption is indicated by about 30 cm of vesicular basalt. Al Hashimiyya basalts eruptions form successive sheet-like flows with stepped cliffs.
The basalts erupted in this area parallel to Wadi Sirhan fault system extending NW-SE (Figure 1). Al Hashimiyya flow was selected after taking different determining criteria into concern; 1) location as it’s the south western utmost end of basaltic flow in the region, 2) the outcrop accessibility, 3) structural aspects of presenting moderate (upper sub-flows) to low (lower sub-flows) joints and fractures, 4) basalt alteration is uncommon under microscope as in the other flows in the eastern area [16] .
The total studied thickness is about 20 m, the bottom sub-flow elevation is 512 m a.s.l. The upper two sub-flows are characterized with humpy ridges structure and hexagonal cooling joints. The blocky structure (of Aa lava type) is exhibited by the middle and the lower sub-flows. In general, Al Hashimiyya flow is characterized by tectonic fractures parallel with the preferred orientation of main fault system of NW-SE direction.
3. Methodology
Total 20 m thickness of basaltic sub- flows was investigated. Two main complementary sections were taken, along the flows with 10 to 14 samples from each sub-flow. The samples were prepared into predefined dimensions or milled for the thermo-physical and mechanical characterization.
Optical Scanning Method was used to measure thermal conductivity [19] . Thermal conductivity was conducted under standard laboratory conditions on oven dried samples. Optical Scanning Method is based on scanning of sample surface (curved or flat) with a straight and constantly run mobile heat source of 4˚C temperature. This heat source is associated with a temperature sensor. Such device is one of the quick thermal conductivity measuring apparatus. However, different sample types can be measured considering thermal conductivity anisotropy. The measuring conception is as follows; a 1 mm diameter mobile heat source is heating with scanning the examined surface. The sample is few centimeters above the mobile heat source. The accompanied (with fixed distance a part) temperature sensor is moving parallel with the heat source in a constant speed. Sample thermal conductivity is calculated as the arithmetic mean of narrow successive measurements on the sample surface. Before and after heat measurements are controlling the local scanned measurements. On the scan line, two known thermal conductivity standards are located prior and after the sample. The standards have almost the same thermal conductivity (λR) as the sample. The comparison between the studied sample and the generated temperature differences of the standards will determine the sample thermal conductivity λ (Equation (1)) with an accuracy of 3%. The measured electrical potential for both samples UR and U is dependent on the ratio of temperature rise for standard and sample
and
, as described below [19] :
(1)
The heat capacity was calculated out of the thermal conductivity, temperature diffusivity and bulk density. Permeability measurements were carried out using a pressure air driven gas-permeameter. The gas pressure permeameter allows an infinite measurements, this can be allocated to distinct layers of a rock sample. The air is pressed into the sample surface through a vertical detached ram. Compressed air is driven through an outlet to initiate a specific flow-through. This will lead to initiate a volumetric flow rate. Together with the dynamic viscosity and adaptor radius, the permeability is calculated as a result of volumetric flow rate (Equations (2) & (3)) [29] , the error is limited to a maximum deviation of 5%.
(2)
where
is the volumetric flow rate,
is the atmospheric pressure,
is the injection pressure and
is the mass flow rate. The permeability (k) is calculated then as follow
(3)
where
is the dynamic viscosity, r is the adapter radius and
is a geometric factor.
Porosity was measured using pycnometer. This device measures the porosity considering a known envelope volume, bulk density and the particle density of an examined cylindrical sample. In general, with an accuracy of 0.03%, the device computes both densities by measuring the displacement generated after applying a pressure on the sample including pore volume. This is measured in a cylindrical cell, filled with powder, ends with a mobile ram. The concept of measuring is to calculate the difference between two readings; the first is the measured distance that powder needs to reach a predefined pressure and second, is the measured distance needed to reach the previous pressure for both powder and sample together. The pores volume is calculated by integrating the sample oven dry weight with the particle density [30] .
Unconfined compression strength (
) is the most important strength tests among rock physical properties. The maximum applied uniaxial stress that a rock sample can stand without failure is its uniaxial compressive strength; it is calculated using the following (Equation (4)):
(4)
where:
= The maximum applied load, KN.
A = Area of the specimen cross section, cm2.
For ultrasonic wave velocity measurement, a sample is aligned between two transducers. As the device turned on, wave pulses are transmitted from one transducer to the other through the sample. Travel time of the transmission is given in microseconds, and thus the velocity can be calculated considering the length of the specimen.
One of the basic rock tests which indicates the failure strength is tensile strength test. Before conducting the tensile strength test, basalt specimens were prepared according the (ISRM) standards. The prepared sample was placed between two platens so that its diameter is perpendicular to the platens plane. The load was applied at a rate is of 0.22 (KN/sec), until failure was reached. Consequently, cracks starting at this region and propagate along the plane parallel to the direction of load application until the basalt sample split into two parts. The indirect tensile compressive strength (
) can be represented mathematically using the following relationship (Equation (5)):
(5)
where:
P = Maximum load applied, Kg.
D = Core diameter, cm.
l = Core thickness, cm.
Point load test is that test examined rock’s response on tension stress, thus it belongs to the tensile strength tests family. This test has three variations; diametric, axial or irregular lumps, the choice of which depends on the available specimen geometry [31] . In this test basalt samples were loaded horizontally between two thin ends of the device those are compressed perpendicular to the specimen until fracture occurred. The sample then is split ruptured into four pieces. The load strength index (Is) of the examined basalt was calculated as follows (Equation (6)):
(6)
where:
P = failure load, KN.
D2 = Core diameter, cm.
4. Results and Discussion
Correlation between thermal conductivity, permeability, porosity, thermal diffusivity, density, specific heat capacity, unconfined compressive strength, compressional wave velocity, indirect tensile strength and point load are illustrated in Figure 2. Basalts sub-flows are given names H1 to H6. Vertical axis in this figure is not to scale. Table 1 gives the ranges and mean values (with standard deviation) of the measured properties. In addition, a statistical Box-Whisker plot of thermo-physical and mechanical properties are shown in Figure 3.
From the results illustrated above it can be concluded that basalts thermal conductivity correlates with several physical properties; thermal diffusivity, specific heat capacity, unconfined compressive strength and point load. In contrast, it’s inversely proportioned with permeability, porosity and indirect tensile strength. Density and compressional wave velocity does not show obvious correlation with thermal conductivity. Thermo-physical properties of Al Hashimiyya
![]()
Table 1. Mean and ranges; minimum and maximum values of the examined basalts properties.
![]()
Figure 2. Correlations of thermo-physical and mechanical properties of Al Hashimiyya basalts.
![]()
Figure 3. Box-Whisker plots of thermo-physical and mechanical properties of the studied basalts.
basalts are in line with eastern Jordanian Harrat basalt properties reported by [9] . Thermal conductivity is proportionally correlates with some mineral constituents of basalts. A predicting model which link basalt thermal conductivity with its mineral constituents in this area has been approached [16] .
Unconfined compression strength is the most important index strength test as well as the most important property needed for rock classification. It’s found that average value of unconfined strength at Al-Hashimiyya basalts of 104.9 MPa is match very well with the value of 105.8 MPa reported by [32] in Northern Jordan basalt. Both values are in the range of continental basalt unconfined compressive strength values reported earlier. On the other hand, most rocks exhibit a weak tension strength value compared with compressive one, this weakness is very important in structure and excavation constructions. Therefore, indirect tensile strength is found to be one fifth of the unconfined compressive strength of Al Hashimiyya basalts. In addition, compressional wave velocity is in the range of basalts reported by [33] , which is of (0.3940 - 0.5150) km/s. Compressional wave velocity may be attributed to the values of porosity, density of micro fractures, mineral composition and dry density. [34] obtained highly correlated models (R2 1/4 0.717 - 0.890) between P-wave velocity and density, porosity, uniaxial compressive strength, Brazilian tensile strength, modulus of elasticity and Poisson’s ratio for basaltic rocks in other region. [32] reported that the ultrasonic wave velocity in basalt is 0.5399 km/s, it’s slightly lower than those of Al-Hashimiyya basalt which averages 0.5548 km/s. However, both readings are within the continental basalt range of ultrasonic wave velocity which is (0.35 - 0.55) km/s [33] . Point load strength averages 6.08 Mpa. Relatively, low point load strength values are due to the stress impact distribution over the objected area to the device platens during the test. This test is conducted on un-prepared ends cores, where the diameter is the most important parameter in determining the point load strength.
A number of factors may affect the results; rock weathering degree, micro fractures and porosity. In addition, other factors, such as moisture content, sample surface flatness, loading rate, specimen dimensions and surface shape, might have also great impact on the results. In addition, for aggregates’ mechanical and physical properties are affected by texture, size, shape, chemical, and mineralogical compositions ( [35] [36] [37] ).
Thermo-physical and mechanical rock properties were determined at oven dried samples with an ambient temperature and pressure. This will enable a direct comparison of all examined rock samples. The actual conditions on depth do not reflected by the laboratory conditions. Deep conditions include; higher pressure, higher temperature, water saturation, steam with other gases contents and surrounding hydrochemistry. In compare, higher temperatures at depths, parallel with increasing pressure, may change the values of porosity, permeability, unconfined compressive strength, indirect tensile strength and point load strength from the laboratory measurements. This effect is caused by the development of micro fractures at higher depths. Rock thermal expansion, parallel with pressure effect, will form micro and macro cracking. This will develop a second porosity, increasing the total porosity value. In addition, permeability is expected to be increased, too. In contrast, rock strength may be reduced, as well ( [38] [39] [40] [41] [42] ). The effect of such cracks on density, heat and compressional wave velocity is less noticeable in saturated rocks. This is probably caused by contact resistance reduction by the fluids [14] . [43] and [44] stated that if the porosity is increased by a factor of 3 to 4 then the strength will be decreased by a factor of 8. If this pores and cracks getting saturated with water then the thermal conductivity and compressional wave propagation are expected to be increased as water conduction is higher than air.
Basalt rocks, as a crystalline igneous rock, have tightly interlocked crystals which affect heat and compressional waves. This can be seen obviously comparing with other rock types [14] . However, the increasing in grains contact resistance between mineral constituents may decrease heat and ultrasonic wave propagation by up to 50%. Consequently, heat and compressional wave conduction is expected to duplicate in saturated porous reservoir rocks ( [41] [42] [45] [46] [47] ).
However, studies on thermo-physical rock properties evaluations under high pressure, temperature and fluid saturation are very rare, due to the technical effort needed. According to [48] and [49] , the existing apparatus are typically capable of reaching sufficient pressure conditions, but the maximum temperature is often limited by 200˚C, and only selected rock properties such as permeability, acoustic velocity and rock strength can be measured.
In addition to temperature, pressure and water saturation, the scale dependence limitations affect results, too. Laboratory measurements will considerably differs from in situ reservoir measurements due to the small size and relative homogeneity of the examined small samples. In-situ measurements cover a much larger volume of heterogeneous rocks [45] .
Furthermore, geothermal systems may modify reservoir rocks properties due to the continuous physical and chemical condition changes. Hydrothermal alteration, as an example, can either increase thermal conductivity, wave velocity, permeability and rock strength by precipitates secondary minerals in pores and fractures, or, decrease them by dissolving the primary mineral constituents leaving cavities filled with new different secondary minerals. In depth, it’s expected for the basalts porosity and permeability to be increased by micro and macro fractures. However, hydrothermal alteration may seal those cracks and reduce their effect on thermo-physical and mechanical properties. Porosity and fluid saturation will control thermo-physical properties of basalts, too. This is because of relatively high porosity exhibited by the examined basalts flows ( [13] [43] [44] [50] [51] ). [52] discussed the saturation effects on thermal properties. They found that thermal conductivity and thermal diffusion coefficient of rocks increases under water saturated conditions compared to dry conditions, but the specific heat capacity decreases. In addition, [53] indicated that the rock lithology is the predominant control on the physical and mechanical properties of the geotechnical units. They suggested that this correlation is due to mineral precipitation within fractures and pores in the brecciated lava margins. They concluded that permeability of the unaltered breccia facilitated efficient hydrothermal fluid circulation and mineral precipitation.
Therefore, thermo-physical and mechanical rock properties changes with depth are expected to be highly variable. Hence, predicting physical and mechanical properties of reservoir rocks at depth become very complicated and un applicable, only estimations can be done.
5. Conclusions
This work examined thermo-physical and mechanical rock properties evaluation of basaltic rocks in Al Hashimiyya area. Samples covered about 20 m thickness of six successive sub-flows. This study provides the most comprehensive thermo-physical and mechanical properties data suite on the Jordanian Harrat basaltic rock to date. Previous work investigated some of the properties separately. Thus, such developed correlations between thermal conductivity and other rock properties provide a unique and new reference data set for less studied thermo-physical basalt properties in Al Hashimiyya. The results outlined in this work are expected to develop a new approach for predicting rock properties into a pre-defined thermo-physical model. Such models may be expanded when the input data are limited, due to depth conditions, based on different patterns of correlations given in this work.
Study outcomes reveal that thermo-physical and mechanical properties of Al Hashimiyya basalts are of regional ranges of continental basalts. Al Hashimiyya basalts thermal conductivity averages at 1.55 W・m−1・K−1. Porosity ranges between 12.76% and 7.5%. The basalts density in the study area averages 2.68 kg・m−3. The indirect tensile strength averages 5.16 Mpa and is one fifth the unconfined compressive strength which ranges from 168.7 to 41.2 Mpa.
The results indicated that thermal conductivity is proportionally correlated with thermal diffusivity, specific heat, unconfined compressive strength and point load. While other properties like permeability, porosity and indirect tensile strength show an inverse proportion with thermal conductivity. Density and compressional wave velocity does not exhibit a direct correlation with thermal conductivity. Different factors, which affect the results in laboratory and depth, had been discussed. Permeability and porosity are expected to be the most controlling factors on thermal conductivity results. They also affect other mechanical properties of basalts. In addition, temperature, pressure, fluid saturation and hydrothermal solutions are expected to change the results values in depth.
In depth, in addition to the microcracks, the relative heterogeneity in pores shape, volume and distribution of basaltic sub-flows would absolutely affect the rock properties. Hydrothermal solutions will definitely change pores and fractures effects in two different ways: cavities improvement or existing pores and cracks filling. Therefore, reservoir conditions will limit the ability to predict thermophysical and mechanical properties in depth.