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.     

Influence of surface and vegetation characteristics on C-band radar measurements for soil moisture content

Soil moisture estimation using microwave remote sensing faces challenges of the segregation of influences mainly from roughness and vegetation. Under static surface conditions, it was found that Radarsat C-band SAR shows reasonably good correlation and sensitivity with changing soil moisture. Dynamic surface and vegetation conditions are supposed to result in a substantial reduction in radar sensitivity to soil moisture. A C-band scatterometer system (5.2 GHz) with a multi-polarization and multi-angular configuration was used 12 times to sense the soil moisture over a tall vegetated grass field. A score of vegetation and soil parameters were recorded on every occasion of the experiment. Three radar backscattering models Viz., Integral Equation Model (IEM), an empirical model and a volume scattering model, have been used to predict the backscattering phenomena. The volume scattering model, using the Distorted Born Approximation, is found to predict the backscattering phenomena reasonably well. But the surface scattering models are expectedly found to be inadequate for the purpose. The temporal variation of soil moisture does show good empirical relationship with the observed radar backscattering. But as the vegetation biomass increases, the radar shows higher sensitivity to the vegetation parameters compared to surface characteristics. A sensitivity analysis of the volume scattering model for all the parameters also reveals that the radar is more sensitive to plant parameters under high biomass conditions, particularly vegetation water content, but the sensitivity to surface characteristics, particularly to soil moisture, is also appreciable.     

Genesis of tin and tantalum mineralization in pegmatites from the Bynoe area,Pine Creek Geosyncline,Northern Territory

The tin‐ and tantalum‐bearing pegmatites of the Bynoe area are located in the western Pine Creek Geosyncline. They are emplaced within psammopelitic rocks in the contact aureole of the Two Sisters Granite. The latter is a Palaeoproterozoic, fractionated, granite with S‐type characteristics and comprises a syn‐ to late‐orogenic, variably foliated, medium‐grained biotite granite and a post‐orogenic, coarse‐grained biotite‐muscovite granite. The pegmatites comprise a border zone of fine grained muscovite + quartz followed inward by a wall zone of coarse grained muscovite + quartz which is in turn followed by an intermediate zone of quartz + feldspar + muscovite. A core zone of massive quartz is present in some occurrences. Feldspars in the intermediate zone are almost completely altered to kaolinite. This zone contains the bulk of cassiterite, tantalite and columbite mineralization. Fluid inclusions in pegmatitic quartz indicate that early Type A (CO₂ + H₂O ± CH₄) inclusions were trapped at the H₂O‐CO₂ solvus at P～100 MPa, T～300°C (range 240–328°C) and salinity ～6 wt% eq NaCl. Pressure‐salinity corrected temperatures on Type B (H₂O + ～20% vapour), C (H₂O + < 15% vapour) and D (H₂O + halite + vapour) inclusions also fall within the range of Type A inclusions. Oxygen and hydrogen isotope data show that kaolin was either formed in isotopic equilibrium with meteoric waters or subsequent to its formation, from hydrothermal fluid, underwent isotopic exchange with meteoric waters. Fluid inclusion waters from core zone quartz show enrichment in deuterium suggesting metamorphic influence. Isotope values on muscovite are consistent with a magmatic origin. It is suggested that the pegmatites were derived from the post‐orogenic phase of the Two Sisters Granite. Precipitation of cassiterite took place at about 300°C from an aqueous fluid largely as a result of increase in pH due to feldspar alteration.     

Understanding nesting ecology and behaviour of green marine turtles at Setiu,Terengganu, Malaysia

In this paper, we emphasize the importance of understanding the nesting ecology and nesting behaviour of green marine turtles (Chelonia mydas). Data were collected from 2007 until 2012 from nesting beaches at Setiu Terengganu, Malaysia. We focused on one of the beaches, Telaga Papan, based on data collected in 2012. We recorded the distribution of nesting areas, the emergence hour and the correlation between successful nesting attempts and false crawls. Telaga Papan had a significantly higher distribution of green marine turtle nesting compared with the other five beaches (ANOVA, F_5,42 = 8.874, P < 0.01, mean = 36.750 ± 3.727). The highest number of successful nesting attempts was recorded in 2012 (mean = 28.714). A majority of the species landed between 22:00 and 23:59 h (25%). There was a strong correlation between successful nesting attempts and false crawls (r_s = 0.883, P = 0.02). Based on these findings on the nesting ecology and nesting behaviour of green marine turtles, we suggest that scientific research, strict monitoring, awareness programs and policy implementation should be carried out proactively. Such activities are necessary to reduce the anthropogenic pressures at the nesting beaches as well as to ensure more successful nesting attempts of green marine turtles in Setiu.