Numerical estimation of carbonate rock properties using multiscale images

Characterizing the pore space of rock samples using three‐dimensional (3D) X‐ray computed tomography images is a crucial step in digital rock physics. Indeed, the quality of the pore network extracted has a high impact on the prediction of rock properties such as porosity, permeability and elastic moduli. In carbonate rocks, it is usually very difficult to find a single image resolution which fully captures the sample pore network because of the heterogeneities existing at different scales. Hence, to overcome this limitation a multiscale analysis of the pore space may be needed. In this paper, we present a method to estimate porosity and elastic properties of clean carbonate (without clay content) samples from 3D X‐ray microtomography images at multiple resolutions. We perform a three‐phase segmentation to separate grains, pores and unresolved porous phase using 19 μm resolution images of each core plug. Then, we use images with higher resolution (between 0.3 and 2 μm) of microplugs extracted from the core plug samples. These subsets of images are assumed to be representative of the unresolved phase. We estimate the porosity and elastic properties of each sample by extrapolating the microplug properties to the whole unresolved phase. In addition, we compute the absolute permeability using the lattice Boltzmann method on the microplug images due to the low resolution of the core plug images. In order to validate the results of the numerical simulations, we compare our results with available laboratory measurements at the core plug scale. Porosity average simulations for the eight samples agree within 13%. Permeability numerical predictions provide realistic values in the range of experimental data but with a higher relative error. Finally, elastic moduli show the highest disagreements, with simulation error values exceeding 150% for three samples.     

Seismic design factors for RC special moment resisting frames in Dubai, UAE

This study investigates the seismic design factors for three reinforced concrete (RC) framed buildings with 4, 16 and 32-stories in Dubai, UAE utilizing nonlinear analysis. The buildings are designed according to the response spectrum procedure defined in the 2009 International Building Code (IBC’09). Two ensembles of ground motion records with 10% and 2% probability of exceedance in 50 years (10/50 and 2/50, respectively) are used. The nonlinear dynamic responses to the earthquake records are computed using IDARC-2D. Key seismic design parameters are evaluated; namely, response modification factor (R), deflection amplification factor (Cd), system overstrength factor (Ωo), and response modification factor for ductility (Rd) in addition to inelastic interstory drift. The evaluated seismic design factors are found to significantly depend on the considered ground motion (10/50 versus 2/50). Consequently, resolution to the controversy of Dubai seismicity is urged. The seismic design factors for the 2/50 records show an increase over their counterparts for the 10/50 records in the range of 200%-400%, except for the Ωo factor, which shows a mere 30% increase. Based on the observed trends, period-dependent R and Cd factors are recommended if consistent collapse probability (or collapse prevention performance) in moment frames with varying heights is to be expected.     

Geographic analysis of manufacturing industry in Saudi Arabia

The study focuses on the analysis of manufacturing industry in Saudi Arabia by examining development, structure and spatial distribution. The analysis has shown that Saudi manufacturing is basically of recent development and is highly capital- and energy-intensive, particularly the hydrocarbon-based groups. Perhaps the most conspicuous weakness of Saudi industrialization is the almost total dependence on expatriate labour. Foreign capital and expertise share in the building and development of the industrial sector by providing technical, managerial, and marketing skills. Saudi industry is heavily subsidized by the government, and the public sector plays a commanding role in industrial development. Large-scale basic industries are owned and run by the government, while light manufacturing of the import-substitution type is dominated by the private sector. One of the main characteristics of manufacturing industry is that it is scantily diversified, modestly concentrated and remarkably linked geographically but not functionally or structurally. Although the authorities are exerting tremendous efforts to achieve a regional balance in industrial development, Saudi manufacturing industry is significantly polarized, urban-oriented and highly associated with the major axes of population scatter. Like all developing countries beginning to industrialize, Saudi manufacturing industry faces certain problems which have to be resolved so as to lay the foundations for a sound and prosperous industrial future. The paper goes on to indicate problems and to suggest solutions.     

A Genetic Programming–Based Model for Colloid Retention in Fractures

Understanding the behavior of colloids in groundwater is critical as some are pathogenic while others may facilitate or inhibit the transport of dissolved contaminants. Colloid behavior in saturated fractured aquifers is governed by the physical and chemical properties of the groundwater-particle-fracture system. The interaction between these properties is nonlinear, and there is a need for a mathematical model describing the relationship between them to advance the mechanistic understanding of colloid transport in fractures and facilitate modeling in fractured environments. This paper coupled genetic programming and linear regression within a multigene genetic programming framework to develop a robust mathematical model describing the relationship between colloid retention in fractures and the physical and chemical parameters that describe the system. The data employed for model development and validation were collected from a series of 75 laboratory-scale colloid tracer experiments conducted under a range of conditions in three laboratory-induced discrete dolomite fractures and their epoxy replicas. The model sufficiently reproduced the observed data with coefficients of determination (R²) of 0.92 and 0.80 for model development and validation, respectively. A cross-validation demonstrated the model generality to 86% of the observed data. A variance-based global sensitivity analysis confirmed that attachment is the primary retention mechanism in the systems employed in this work. The model developed in this study provides a tool describing colloid retention in factures, which furthers the understanding of groundwater-particle-fracture system conditions contributing to the retention of colloids and can aid in the design of groundwater remediation strategies and development of groundwater management plans.     

Gravity map of Kalabsha area, northwest of Aswan Lake, and its structural significance

A detailed gravity survey was carried out in one of the seismo-active areas at the northwestern region of the High Dam Lake (Kalabsha area) to study its subsurface structure. In order to understand the seismicity of the area, the establishment of a geodynamic model from geological and geodetic data is of great importance. After a series of adjustments and corrections of the measured gravity data, free-air and Bouguer anomaly maps were constructed for the Kalabsha area, and several interpretation techniques were applied to analyse these anomalies. The results of the analysis indicate that the Kalabsha area is affected by several faults trending mainly E-W and N-S. The active area located west of Gebel Marawa is bounded by a set of faults striking NE-SW, N-S and E-W. The throws of these faults range from 160 to 370 m.The minimum depth to the basement complex is about 200 m and its maximum depth is around 600 m. The thickness of the sedimentary column (Nubia sandstone) in the Kalabsha area decreases due west and increases toward the southern and eastern parts of the area. The results explain the tectonic framework of the area well.     

Influence of Degree of Saturation in the Electric Resistivity–Hydraulic Conductivity Relationship

The relationship between aquifer hydraulic conductivity and aquifer resistivity, either measured on the ground surface by vertical electrical sounding (VES) or from resistivity logs, or measured in core samples have been published for different types of aquifers in different locations. Generally, these relationships are empirical and semi-empirical, and confined in few locations. This relation has a positive correlation in some studies and negative in others. So far, there is no potentially physical law controlling this relation, which is not completely understood. Electric current follows the path of least resistance, as does water. Within and around pores, the model of conduction of electricity is ionic and thus the resistivity of the medium is controlled more by porosity and water conductivity than by the resistivity of the rock matrix. Thus, at the pore level, the electrical path is similar to the hydraulic path and the resistivity should reflect hydraulic conductivity. We tried in this paper to study the effect of degree of groundwater saturation in the relation between hydraulic conductivity and bulk resistivity via a simple numerical analysis of Archie’s second law and a simplified Kozeny-Carmen equation. The study reached three characteristic non-linear relations between hydraulic conductivity and resistivity depending on the degree of saturation. These relations are: (1) An inverse power relation in fully saturated aquifers and when porosity equals water saturation, (2) An inverse polynomial relation in unsaturated aquifers, when water saturation is higher than 50%, higher than porosity, and (3) A direct polynomial relation in poorly saturated aquifers, when water saturation is lower than 50%, lower than porosity. These results are supported by some field scale relationships.     

Shear waves in elastic medium with void pores welded between vertically inhomogeneous and anisotropic magnetoelastic semi-infinite media

The paper intends to study the propagation of horizontally polarized shear waves in an elastic medium with void pores constrained between a vertically inhomogeneous and an anisotropic magnetoelastic semi-infinite media. Elasto-dynamical equations of elastic medium with void pores and magnetoelastic solid have been employed to investigate the shear wave propagation in the proposed three-layered earth model. Method of separation of variables has been incorporated to deduce the dispersion relation. All possible special cases have been envisaged and they fairly comply with the corresponding results for classical cases. The role of inhomogeneity parameter, thickness of layer, angle with which the wave crosses the magnetic field and anisotropic magnetoelastic coupling parameter for three different materials has been elucidated and represented by graphs using MATHEMATICA.     

Regularized Bayesian estimation for GEV-B-splines model

Large observed datasets are not stationary and/or depend on covariates, especially, in the case of extreme hydrometeorological variables. This causes the difficulty in estimation, using classical hydrological frequency analysis. A number of non-stationary models have been developed using linear or quadratic polynomial functions or B-splines functions to estimate the relationship between parameters and covariates. In this article, we propose regularised generalized extreme value model with B-splines (GEV-B-splines models) in a Bayesian framework to estimate quantiles. Regularisation is based on penalty and aims to favour parsimonious model especially in the case of large dimension space. Penalties are introduced in a Bayesian framework and the corresponding priors are detailed. Five penalties are considered and the corresponding priors are developed for comparison purpose as: Least absolute shrinkage and selection (Lasso and Ridge) and smoothing clipped absolute deviations (SCAD) methods (SCAD1, SCAD2 and SCAD3). Markov chain Monte Carlo (MCMC) algorithms have been developed for each model to estimate quantiles and their posterior distributions. Those approaches are tested and illustrated using simulated data with different sample sizes. A first simulation was made on polynomial B-splines functions in order to choose the most efficient model in terms of relative mean biais (RMB) and the relative mean-error (RME) criteria. A second simulation was performed with the SCAD1 penalty for sinusoidal dependence to illustrate the flexibility of the proposed approach. Results show clearly that the regularized approaches leads to a significant reduction of the bias and the mean square error, especially for small sample sizes (n < 100). A case study has been considered to model annual peak flows at Fort-Kent catchment with the total annual precipitations as covariates. The conditional quantile curves were given for the regularized and the maximum likelihood methods.