Modeling three‐dimensional hydraulic fracture propagation using virtual multidimensional internal bonds

Propagation of fractures, especially those emanating from wellbores and closed natural fractures, often involves Mode I and Mode II, and at times Mode III, posing significant challenges to its numerical simulation. When an embedded inclined fracture is subjected to compression, the fracture edge is constrained by the surrounding materials so that its true propagation pattern cannot be simulated by 2D models. In this article, a virtual multidimensional internal bond (VMIB) model is presented to simulate three‐dimensional (3D) fracture propagation. The VMIB model bridges the processes of macro fracture and micro bond rupture. The macro 3D constitutive relation in VMIB is derived from the 1D bond in the micro scale and is implemented in a 3D finite element method. To represent the contact and friction between fracture surfaces, a 3D element partition method is employed. The model is applied to simulate fracture propagation and coalescence in typical laboratory experiments and is used to analyze the propagation of an embedded fracture. Simulation results for single and multiple fractures illustrate 3D features of the tensile and compressive fracture propagation, especially the propagation of a Mode III fracture. The results match well with the experimental observation, suggesting that the presented method can capture the main features of 3D fracture propagation and coalescence. Moreover, by developing an algorithm for applying pressure on the fracture surfaces, propagation of a natural fracture is also simulated. The result illustrates an interesting and important phenomenon of Mode III fracture propagation, namely the fracture front segmentation. Copyright © 2012 John Wiley & Sons, Ltd.     

Neuro-fuzzy models employing wavelet analysis for suspended sediment concentration prediction in rivers

Abstract

The study of sediment load is important for its implications to the environment and water resources engineering. Four models were considered in the study of suspended sediment concentration prediction: artificial neural networks (ANNs), neuro-fuzzy model (NF), conjunction of wavelet analysis and neuro-fuzzy (WNF) model, and the conventional sediment rating curve (SRC) method. Using data from a US Geological Survey gauging station, the suspended sediment concentration predicted by the WNF model was in satisfactory agreement with the measured data. Also the proposed WNF model generated reasonable predictions for the extreme values. The cumulative suspended sediment load estimated by this model was much higher than that predicted by the other models, and is close to the observed data. However, in the current modelling, the ANN, NF and SRC models underestimated sediment load. The WNF model was successful in reproducing the hysteresis phenomenon, but the SRC method was not able to model this behaviour. In general, the results showed that the NF model performed better than the ANN and SRC models.

Citation Mirbagheri, S. A., Nourani, V., Rajaee, T. & Alikhani, A. (2010) Neuro-fuzzy models employing wavelet analysis for suspended sediment concentration prediction in rivers. Hydrol. Sci. J. 55(7), 1175–1189.     

Isotopic analysis to quantify the role of the Indian monsoon on water resources of selected river basins in the Himalayas

Western disturbances (WDs) and Indian summer monsoon (ISM) led precipitation play a central role in the Himalayan water budget. Estimating their contributions to water resource is although a challenging but essential for hydrologic understanding and effective water resource management. In this study, we used stable water isotope data of precipitation and surface waters to estimate the contribution of ISM and WDs to the water resources in three mountainous river basins - Indus, Bhagirathi and Teesta river basins of western, central and Eastern Himalayas. The study reveals distinct seasonality in isotope characteristics of precipitation and surface waters in each river basin is due to changes in moisture source, hydrometeorology and relief. Despite steady spatial variance in the slope and intercept of regression lines from the Teesta to Indus and the Bhagirathi river basins, the slope and intercept are close to the global meteoric water line and reported local meteoric water line of other regions in the Himalayas and the Tibetan Plateau. The two-component end-member mixing method using d-excess as tracer were used to estimate the contribution from ISM and WD led precipitation to surface water in aforementioned river basins. The results suggest that the influence of the ISM on the water resources is high (>72% to annual river flow) in Teesta river basin (eastern Himalayas), while as the WDs led precipitation is dominantly contributing (>70% average annual river flow) to the surface waters in the Indus river basin (western Himalayas). The contribution of ISM and WD led precipitation in Bhagirathi river basin is 60% and 40%, respectively. The findings demonstrate that the unusual changes in the ISM and WD moisture dynamics have the potential to affect the economy and food security of the region, which is dependent on the availability of water resources. The obtained results are of assistance to policy makers/mangers to make use of the information for better understanding hydrologic response amid unusual behaviour of the dual monsoon system over the region.     

Dynamic vertical compliance of strip foundations in layered soils

In this paper, results of a detailed investigation on the dynamic response of rigid strip foundations in viscoelastic soils under vertical excitation are presented. An advanced boundary element algorithm developed by incorporating isoparametric quadratic elements and a sophisticated self-adapting numerical integration scheme has been used for this investigation. Foundations supported on three types of soil profiles, namely, homogeneous half-space, stratum-over-half-space and stratum-over-bedrock are considered. The influence of material properties like Poisson's ratio and material damping as well as the influence of geometrical properties such as depth of embedment and layer thickness are studied. The effect of the type of contact at the soil-foundation interface is also investigated.     

Stiffness identification of a tall building during construction period using ambient tests

Ambient response measurements were made on an eighteen-storey building at three different stages of construction to detect any changes in the frequencies, mode shapes and stiffness with construction. The first nine frequencies and corresponding mode shapes for each stage of construction are found. A comparison is made among these mode shapes and frequencies and with the mode shapes and frequencies of an analytical model incorporating beams, columns, shear walls, panels and diagonal elements. The added effects, on frequencies and mode shapes, of non-structural elements such as stairs, elevators, claddings and partition walls are studied. Using Improved Statistical Structural Identification, an attempt is made to study the stiffening effect of non-structural elements by updating the stiffness matrix of the building.     

Areal extent of the 1988 flood in Bangladesh: How much did the satellite imagery show?

This paper presents the methods and results of visual interpretation of satellite imagery for estimating areal extent of the 1988 flood in Bangladesh. The main limitation of this method was that cloud-free images were not available for the peak flood period in the first week of September. Relatively cloud-free images on three selected dates in the last three weeks of September showed that areas of inundation on those dates ranged from 44000 to 60000 km² ± 5% (31 to 42% of Bangladesh). The Government of Bangladesh official estimates of the areal extent of flooding were tentative in nature, since these were based on a choropleth map of relative proportions (percentages) of flooding in different parts of the country. The official estimate of the maximum extent of flooding was 82000 km² (57% of Bangladesh), i.e. 15 to 26% larger than the area shown on satellite imagery. The actual extent of flooding was certainly larger than the estimates from satellite imagery but, perhaps, smaller than the official estimates, since a reduction of 15 to 26% of flood area in one to three weeks seemed unlikely. An alternative method of mapping flood-affected areas by using newspaper-interpreted data was attempted, but the method had limited values because of reporting bias.     

Effects of search parameters on kriged reserve estimates

Summary Reliable ore reserve estimates for deposits with highly skewed grade distributions are difficult tasks to perform. Although some recent geostatistical techniques are available to handle problems with these estimations, ordinary kriging or conventional interpolation methods are still widely used to estimate the ore reserves for such deposits. The estimation results can be very sensitive to the search parameters used during the interpolation of grades with these methods.This paper compares the ore reserve estimates from ordinary kriging using several cases in which certain search parameters are varied. The comparisons are extended to different mineralizations to show the changing effects of these parameters.     

Protected areas in Saudi Arabia: Sustainable use of natural resources

The present paper reviews the conservation movement in Saudi Arabia as measured against the established protected areas, as well as the basic philosophy regarding natural resource management. The degree of representation of the biophysical diversity of Saudi Arabia in the established protected areas is discussed by using three areas as case studies: Harrat Al-Harrah (lava field), Urug Bani Mu'arid (Cuesta and Sand), and Raydah Escarpment (High Mountain).     

Cadmium content, accumulation and toxicity in dog whelks collected around the Welsh coastline

The concentration of cadmium in sexually mature Nucella lapillus was investigated around the coast of Wales. Higher concentrations of cadmium were found in whelks located in the Bristol Channel compared with Cardigan Bay Nucella. Accumulation and toxicity of cadmium were investigated at 50–600 μg Cd l⁻¹ over a period of 21 days. At the lower concentrations of 50–100 μg Cd l⁻¹ all Nucella survived the 21 days of experimentation. The rate of cadmium uptake after exposure to cadmium was inconsistent and not related to the degree of site contamination. At all sites the concentration of cadmium in Nucella, after treatment with 50–600 μg Cd l⁻¹ for 21 days, was significantly different from controls, but after 7 days it varied with the site and cadmium concentration. Only samples collected at Fishguard, Dale, and Swansea survived 21 days at 600 μg Cd l⁻¹. Tissue cadmium concentrations did not significantly effect the whole body concentration of zinc.