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.     

Yield Forecasting of Spring Maize Using Remote Sensing and Crop Modeling in Faisalabad-Punjab Pakistan

Real time, accurate and reliable estimation of maize yield is valuable to policy makers in decision making. The current study was planned for yield estimation of spring maize using remote sensing and crop modeling. In crop modeling, the CERES-Maize model was calibrated and evaluated with the field experiment data and after calibration and evaluation, this model was used to forecast maize yield. A Field survey of 64 farm was also conducted in Faisalabad to collect data on initial field conditions and crop management data. These data were used to forecast maize yield using crop model at farmers’ field. While in remote sensing, peak season Landsat 8 images were classified for landcover classification using machine learning algorithm. After classification, time series normalized difference vegetation index (NDVI) and land surface temperature (LST) of the surveyed 64 farms were calculated. Principle component analysis were run to correlate the indicators with maize yield. The selected LSTs and NDVIs were used to develop yield forecasting equations using least absolute shrinkage and selection operator (LASSO) regression. Calibrated and evaluated results of CERES-Maize showed the mean absolute % error (MAPE) of 0.35–6.71% for all recorded variables. In remote sensing all machine learning algorithms showed the accuracy greater the 90%, however support vector machine (SVM-radial basis) showed the higher accuracy of 97%, that was used for classification of maize area. The accuracy of area estimated through SVM-radial basis was 91%, when validated with crop reporting service. Yield forecasting results of crop model were precise with RMSE of 255 kg ha^?1, while remote sensing showed the RMSE of 397 kg ha^?1. Overall strength of relationship between estimated and actual grain yields were good with R² of 0.94 in both techniques. For regional yield forecasting remote sensing could be used due greater advantages of less input dataset and if focus is to assess specific stress, and interaction of plant genetics to soil and environmental conditions than crop model is very useful tool.     

Timescales of partial melting in the Himalayan middle crust: insight from the Leo Pargil dome, northwest India

The Leo Pargil dome (LPD) in northwest India exposes an interconnected network of pre-, syn-, and post-kinematic leucogranite dikes and sills that pervasively intrude amphibolite-facies metapelites of the mid-crustal Greater Himalayan sequence. Leucogranite bodies range from thin (5-cm-wide) locally derived sills to thick (2-m-wide) crosscutting dikes extending at least 100 m. Three-dimensional exposures elucidate crosscutting relations between different phases of melt injection and crystallization. Combined laser ablation inductively coupled plasma mass spectrometry U–Th/Pb geochronology and trace element analysis on well-characterized monazite grains from nineteen representative leucogranites yields a large, internally consistent data set of approximately 700 U–Th/Pb and 400 trace element analyses. Grain-scale variations in age correlate with trace element distributions and indicate semi-continuous crystallization of monazite from 30 to 18 Ma. The youngest U–Th/Pb ages in a given sample are consistent with the outcrop-scale crosscutting relations, whereas older ages within individual samples record inheritance from partially crystallized melt and source metapelites. U–Th/Pb isotopic and trace element data are incorporated into a model of melting within the LPD that involves (1) steady-state equilibrium batch melting of compositionally homogeneous metapelitic sources; (2) pulses of increased melt mobility lasting 1–2 m.y. resulting in segregation of melt from its source and amalgamation into mixed magmas; and (3) rapid emplacement and final crystallization of leucogranite bodies. Melt systems in the LPD evolved from locally derived, in situ melt in migmatitic source rocks into a vast network of dikes and sills in the overlying non-migmatitic host rocks.     

Treatment of Hydrodynamic Loading for Random Sea State

- In order to employ cost effective frequency domain analysis for off-shore structures treatment of hydrodynamic loading is essential. Drag and inertia dominated, resonating and antiresonating cases under random sea states are analyzed to highlight the implications and relative merits of four salient linearization techniques.     

Free Vibration Analysis of Symmetrically Laminated Composite Plates on Elastic Foundation and Coupled with Stationary Fluid

Free vibration analysis of symmetrically laminated composite plates resting on Pasternak elastic support and coupled with an ideal, incompressible and inviscid fluid is the objective of the present work. The fluid domain is considered to be infinite in the length direction but bounded in the depth and width directions. In order to derive the eigenvalue equation, Rayleigh-Ritz method is applied for the fluid-plate-foundation system. The efficiency of the method is proved by comparison studies with those reported in the open literature. At the end, parametric studies are carried out to examine the impact of different parameters on the natural frequencies.     

Photogrammetric Capabilities of the Kodak DC40, DCS420 and DCS460 Digital Cameras

The aim of this study is to assess the impact of sensor size and resolution of different digital camera sensors upon the accuracy and precision of three dimensional data derived by photogrammetry. Kodak DC40, DCS420 and DCS460 digital cameras were used to produce digital images of retroreflective targets in a 4 m three dimensional test field. The image locations of the targeted points were automatically measured using an off the shelf image processing software package. Two different sub-pixel measurement approaches were examined: centre of gravity and weighted mean. From the automated sub-pixel measurement of the targeted points, results indicate that the high resolution DCS460 camera produces optimum results using either the weighted mean or centre of gravity approaches. Although this was perhaps expected, the far lower resolution DC40 camera performed better than was originally anticipated, suggesting potential for the cheaper DC40 for many applications.     

Development of Tide–Surge Interaction Model for the Coastal Region of Bangladesh

In this paper, a two-dimensional, vertically integrated hydrodynamic model is developed taking into account entrained air bubbles during storm surges as well as incorporating inverted barometer, and river and land dynamics. The model is specifically designed for the coastal region of Bangladesh. A nested scheme method with a fine mesh scheme (FMS), capable of incorporating the complex coastline and all major offshore islands accurately, nested into a coarse mesh scheme (CMS) covering up to 15° N latitude in the Bay of Bengal is used. To incorporate the small and big offshore islands in the Meghna estuarine region with its complex coastline accurately, a very fine mesh scheme (VFMS) is again nested into the FMS. Along the northeast corner of the VFMS, the Meghna river discharge is taken into account. The coastal and island boundaries are approximated through proper stair steps. The model equations are solved by a semi-implicit finite difference technique using a staggered C-grid. A stable appropriate tidal condition over the model domain is generated by applying tidal forcing with the four major tidal constituents M ₂, S ₂, K ₁, and O ₁ along the southern open boundary of the CMS. This tidal regime is introduced as the initial state of the sea for nonlinear interaction of tide and surge. The model is applied to simulate water levels due to the interaction of tide and surge associated with the cyclones April 1991 and Aila at different coastal and island locations along the coast of Bangladesh. The results are found to be quite satisfactory with root mean square error of ～0.50 m as calculated for both the storm events. Tests of sensitivities on water levels are carried out for air bubbles, offshore islands, river discharge, inverse barometer, and grid resolution. The presence of air bubbles increases simulated water levels a little bit in our model, and the contribution of air bubbles in increasing water level is found around 2 %. Further, water levels are found to be influenced by offshore islands, river discharge, inverse barometer as well as grid resolution.     

Household vulnerability and resilience in flood hazards from disaster-prone areas of Punjab,Pakistan

Natural Hazards - An environmental variation has caused Pakistan an alarming portrait of vulnerability in flood disasters. The government has focused on a number of realistic actions, heartening...     

An application of three-dimensional electrical resistivity imaging for the detection of an underground waste-water system

A 3D electrical resistivity imaging survey is presented in this paper. The objective was to investigate an underground wastewater system at the University of Malaya, Malaysia. Apparent resistivity data were collected along ten parallel lines using a Wenner-Schlumberger configuration; electrode cables were oriented in the x-direction with 3 m spacing. Roll-along measurements using a line spacing of 3 m were carried out covering a grid of 20 × 10 electrodes. All data sets were merged into a single data file in order to perform a 3D inversion. Two different 3D least squares algorithms, based on the robust inversion method and the smoothness-constrained technique, were used for the inversion of the apparent resistivity data. Both the horizontal and vertical extents of the anomalous zones found by inversion are displayed. The results indicate the superiority of the robust inversion method over the smoothness-constrained technique at this site. The results are in sufficient accordance with previously known information about the investigation area. The results show that 3D electrical resistivity imaging surveys, in combination with an appropriate 3D inversion method, can be highly useful for engineering and archaeological investigations as well as for environmental applications.