Weak discontinuity in porous media: an enriched EFG method for fully coupled layered porous media

In this paper, a series of multimaterial benchmark problems in saturated and partially saturated two‐phase and three‐phase deforming porous media are addressed. To solve the process of fluid flow in partially saturated porous media, a fully coupled three‐phase formulation is developed on the basis of available experimental relations for updating saturation and permeabilities during the analysis. The well‐known element free Galerkin mesh‐free method is adopted. The partition of unity property of MLS shape functions allows for the field variables to be extrinsically enriched by appropriate functions that introduce existing discontinuities in the solution field. Enrichment of the main unknowns including solid displacement, water phase pressure, and gas phase pressure are accounted for, and a suitable enrichment strategy for different discontinuity types are discussed. In the case of weak discontinuity, the enrichment technique previously used by Krongauz and Belytschko [Int. J. Numer. Meth. Engng., 1998; 41:1215–1233] is selected. As these functions possess discontinuity in their first derivatives, they can be used for modeling material interfaces, generating only minor oscillations in derivative fields (strain and pressure gradients for multiphase porous media), as opposed to unenriched and constrained mesh‐free methods. Different problems of multimaterial poro‐elasticity including fully saturated, partially saturated one, and two‐phase flows under the assumption of fully coupled extended formulation of Biot are examined. As a further development, problems involved with both material interface and impermeable discontinuities, where no fluid exchange is permitted across the discontinuity, are considered and numerically discussed. Copyright © 2014 John Wiley & Sons, Ltd.     

Post-2012 CDM multi-criteria analysis of industries in six Asian countries: Iranian case study

The prospects of the Clean Development Mechanism (CDM) and for carbon income, up to and beyond 2012, in the industrial sectors of Iran and five other Asian countries are investigated. The attractiveness and suitability of each host country, the status of their industrial sectors (based on four post-2012 scenarios), and the post-2012 potential of the CDM (or similar carbon projects) in these sectors are all examined. A multi-criteria analysis of Iran, Saudi Arabia, the UAE, Qatar, China, and India, based on seven sets of criteria (institutional, regulatory, economic, political, social, CDM experience, and energy production/consumption), is conducted, and the post-2012 potential carbon incomes of each country – based on CO₂e emissions of industrial processes – are calculated. Finally, the Iranian industrial sector and the impact of deregulation of energy prices are examined. The post-2012 potential savings in the Iranian industrial sector are calculated based on energy savings, carbon income, and environmental savings. The results indicate that there is strong demand for investment and new technology in this sector to combat several-fold energy price increases. Moreover, high-priced carbon credits could play a meaningful role in post-2012 energy policies in this sector.

Policy relevance

This research is the first study to quantify the carbon market potentials in the industrial sectors of the selected Organization of the Petroleum Exporting Countries (OPEC) members. The Kyoto Protocol is considered by most OPEC countries to be a mixed bag of threats and opportunities and they have shown ambivalence towards it, mainly due to the threat a reduction of fossil fuel consumption poses to their economies. On the other hand, energy efficiency is a desirable goal for their industrial sectors. Iran, as an OPEC member country with vast energy resources, has mostly ignored the CDM during the first commitment period of the Kyoto Protocol and has performed poorly on CDM implementation. However, the current deregulation of energy prices in Iran, with profound cuts in energy subsidies, would definitely alter the perspective of its industrial decision makers on the post-2012 carbon potentials.     

A hydrochemical and dye- tracing investigation in the Posht-e-Naz Karstic Aquifer, Alburz Mountain, Northern Iran

A hydrochemical and tracer study (uranine injection) was conducted in Jurassic limestone of the Posht-e-Naz area in the Alborz belt to evaluate hydraulic relations between a large diameter (about 100 m) sinkhole and springs and aquifer parameters. A main goal of the project was to find out the source of turbidity of the Emarate drinking water supply spring (SP4) in rainy season. Springs discharge were measured and hydrochemical investigation were carried out. The uranine tracer was injected and eight springs, three wells and the Neka River were selected and totally 989 samples in 107 days were collected. Hydrochemical data demonstrated a relative connection between sinkhole (Sh1) and spring (SP4). The results of the tracing by sampling water indicated only a hydraulic connection between Sange- Nou spring (SP8) and injection point, while the charcoal packets analysis revealed tracer exits from spring numbers SP1, SP3, SP4, SP5, SP8, in wells W1 and W2, and in the Neka River. The concentration — time curves of charcoal packets for qualitative analysis and exit tracer for quantitative analysis is also assessed.     

Control of Cambrian evaporites on fracturing in fault-related anticlines in the Zagros fold-and-thrust belt

Orientation and distribution of fractures in the Oligocene–Early Miocene Asmari Formation (a major reservoir rock of the Zagros petroleum system) were investigated in two anticlines of the Zagros fold-and-thrust belt. The Sim and Kuh-e-Asmari anticlines developed in the areas of the Zagros characterized by the occurrence and absence of Cambrian evaporites at the bottom of the stratigraphic pile, respectively. The aim was to outline major differences in terms of fracture spacing and saturation. Organic matter maturity and clay minerals-based geothermometers suggest that the depth of deformation for the top of the Asmari Formation in the Kuh-e-Asmari anticline was in the range of 1.5–2.7 km assuming a geothermal gradient of 22.5 °C/km. The Asmari Formation in the Sim anticline probably experienced a slightly deeper sedimentary burial (maximum 3 km) with a geothermal gradient of 20 °C/km. The spacing of fractures is generally 2–3 times larger (i.e., strain accommodated by fracturing is smaller) in the Sim anticline than in the Kuh-e-Asmari anticline. This is consistent with regional geological studies, analogue, and numerical models that suggest that thrust faults geometry and related folds are markedly different in the absence or presence of a weak decòllement (evaporites). The larger spacing in the Sim anticline is also consistent with higher temperature predicted for the Asmari Formation in this area. By contrast, the orientation of fractures with respect to the fold axes is the same in both anticlines. The fracture systems are rather immature in both anticlines. The amount and density of fractures in the twofolds are controlled by regional (occurrence/absence of salt and probably different burial), rather than local features (fold geometry).     

A study of the behaviour of brittle rocks subjected to confined stress based on the Mohr–Coulomb failure criterion

In this study the relationship between brittle rocks’ behaviour and uniaxial compression stress is obtained based on the Mohr-Coulomb failure criterion and the behavioural characteristics of brittle rocks under uniaxial compression stress are investigated. According to the laboratory uniaxial compression tests on granite rocks, and also available results in the literature, required parameters in the study are obtained. It is indicated that the behaviour of brittle rocks which is affected by uniaxial compression is function of cohesion. In the other words, cohesion changes as the uniaxial compression increases. Also, in yield point where its stress is equivalent to the uniaxial compressive strength, the cohesion maximises. By suggesting a new viewpoint to the Mohr’s circle, normal stresses on failure plane, shear stresses tangential to the failure plane and hydrostatic stresses are investigated. Results show that, normal stress on the failure plane in yield point of the behavioural curve equals zero and shear stress tangential to the failure plane is maximum which is equal to the maximum rocks’ natural strength. Also, in this point the strength is equi-pressure, therefore the stress is of hydrostatic type.     

Review of Laboratory Scale Models of Karst Aquifers: Approaches,Similitude, and Requirements

This review focuses on investigations of groundwater flow and solute transport in karst aquifers through laboratory scale models (LSMs). In particular, LSMs have been used to generate new data under different hydraulic and contaminant transport conditions, testing of new approaches for site characterization, and providing new insights into flow and transport processes through complex karst aquifers. Due to the increasing need for LSMs to investigate a wide range of issues, associated with flow and solute migration karst aquifers this review attempts to classify, and introduce a framework for constructing a karst aquifer physical model that is more representative of field conditions. The LSMs are categorized into four groups: sand box, rock block, pipe/fracture network, and pipe-matrix coupling. These groups are compared and their advantages and disadvantages highlighted. The capabilities of such models have been extensively improved by new developments in experimental methods and measurement devices. Newer technologies such as 3D printing, computed tomography scanning, X-rays, nuclear magnetic resonance, novel geophysical techniques, and use of nanomaterials allow for greater flexibilities in conducting experiments. In order for LSMs to be representative of karst aquifers, a few requirements are introduced: (1) the ability to simulate heterogeneous distributions of karst hydraulic parameters, (2) establish Darcian and non-Darcian flow regimes and exchange between the matrix and conduits, (3) placement of adequate sampling points and intervals, and (4) achieving some degree of geometric, kinematic, and dynamic similitude to represent field conditions.     

Risk assessment of workers exposed to noise pollution in a textile plant

In order to assess the risks associated with worker’s hearing loss, due to exposure to noise pollution, at the Boroujerd Textile Factory a cross sectional study was conducted. 60 workers from the Spinning and Weaving workshops and official staff were randomly selected as case and control groups and their hearing were tested by audiometry. The audiometric results revealed that amongst workers being exposed to the noise pollution in the Spinning and Weaving Workshops there is a distinct increase in the number of cases of hearing loss at high frequencies. There also appears to be a slight bias towards hearing loss in the left ear, in preference to the right in addition, those who have worked for longer than 16 years are subjected to hearing loss even at low frequencies conversely, the office workers (the control group) seem to suffer little or no hearing loss at low frequencies at all, but after approximately 10 years of service there may be some hearing loss, but only at higher frequencies, which could be attributed to presbycusis. Also, the t-test (statistical hypothesis test) results verified the significant difference among both groups at high frequencies. The questionnaire results indicated that workers feel symptoms such as headache, no sense centralization, excitement, nervousness, vertigo. Hence, controlling procedures seem to be essential to protect workers from noise disorders.     

Comparing different statistical models for assessing Fe-contaminated soils based on VNIR/SWIR spectral data

Two statistical models including partial least squares regression (PLSR) and principal component regression were comparatively utilized to determine the predictive accuracy of visible–near-infrared and short-wave infrared reflectance spectroscopy in quantifying the Fe concentration in contaminated soils. Two scenarios were applied to select the best model: Scenario I included all wavelengths (400–2450 nm) and Scenario II encompassed characteristic bands of Fe. Pre-processing techniques used to select the best model included: first and second derivatives (FD and SD), multiplicative scatter correction (MSC) and standard normal variate. The abilities of the predictive models were evaluated by splitting soil samples into two random groups (80 and 20%). The first group (80%) was used to evaluate calibration and validation sets by employing the cross‐validation method, and the second group (20%) was applied to test the models. The coefficient of determination (R ²), root mean square error and residual prediction deviation were calculated to evaluate the models. Applying Scenario I indicated that the PLSR model with SD pre-processing was a more accurate technique for predicting the Fe concentration, whereas in the Scenario II, the PLSR model with MSC pre-processing had a better performance. Comparing Scenarios I and II indicated that the more reliable models for predicting the soil Fe content could be constructed by the PLSR model with the SD pre-processing techniques and all wavelengths. The modeling results produced by the PLSR model with the SD pre-processing could be used to detect, map and monitor Fe-contaminated soils by proximal and remote sensing in the mining areas.