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

K. Huang  Z. Zhang  Ahmad Ghassemi 《国际地质力学数值与分析法杂志》2013,37(13):2021-2038

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.  相似文献   

Water waves generated by instantaneous disturbances at the bed of a sloping beach     

Rumpa Chakraborty 《地球物理与天体物理流体动力学》2013,107(5):481-496

The two-dimensional problem of the generation of water waves due to instantaneous disturbances prescribed at the bed of a beach sloping at an arbitrary angle is studied here. It is formulated in terms of an initial-boundary-value problem for the velocity potential describing the motion in the fluid region assuming the linear theory. Using the Laplace transform in time and the Mellin transform in distance, the problem is reduced to solving a difference equation whose method of solution is of considerable importance in the literature. The form of the free surface is obtained in terms of a multiple infinite integral that is evaluated by the method of steepest-descent. For some prescribed forms of the disturbance at the bed of the beach, the free surface is depicted in a number of figures for different beach angles. It is observed that as the beach angle decreases, the maximum wave height increases, which is plausible.  相似文献   

GIS‐evaluation of two slope‐calculation methods regarding their suitability in slope analysis using high‐precision LiDAR digital elevation models     

M. Irfan Ashraf  Zhengyong Zhao  Charles P.‐A. Bourque  Fan‐Rui Meng 《水文研究》2012,26(8):1119-1133

  相似文献   

A geochronological and petrological study of anatectic paragneiss and associated granite dykes from the Day Nui Con Voi metamorphic core complex,North Vietnam: constraints on the timing of metamorphism within the Red River shear zone     

R. M. Palin  M. P. Searle  D. J. Waters  R. R. Parrish  N. M. W. Roberts  M. S. A. Horstwood  M.‐W. Yeh  S.‐L. Chung  T. T. Anh 《Journal of Metamorphic Geology》2013,31(4):359-387

The Red River shear zone (RRSZ) is a major left‐lateral strike‐slip shear zone, containing a ductilely deformed metamorphic core bounded by brittle strike‐slip and normal faults, which stretches for >1000 km from Tibet through Yunnan and North Vietnam to the South China Sea. The RRSZ exposes four high‐grade metamorphic core complexes along its length. Various lithologies from the southernmost core complex, the Day Nui Con Voi (DNCV), North Vietnam, provide new constraints on the tectonic and metamorphic evolution of this region prior to and following the initial India–Asia collision. Analysis of a weakly deformed anatectic paragneiss using P–T pseudosections constructed in the MnO–Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (MnNCKFMASHTO) system provides prograde, peak and retrograde metamorphic conditions, and in situ U–Th–Pb geochronology of metamorphic monazite yields texturally controlled age constraints. Tertiary metamorphism and deformation, overprinting earlier Triassic metamorphism associated with the Indosinian orogeny and possible Cretaceous metamorphism, are characterized by peak metamorphic conditions of ~805 °C and ~8.5 kbar between c. 38 and 34 Ma. Exhumation occurred along a steep retrograde P–T path with final melt crystallizing at the solidus at ≥~5.5 kbar at ~790 °C. Further exhumation at ~640–700 °C and ~4–5 kbar at c. 31 Ma occurred at subsolidus conditions. U–Pb geochronological analysis of monazite from a strongly deformed pre‐kinematic granite dyke from the flank of the DNCV provides further evidence for exhumation at this time. Magmatic grains suggest initial emplacement at 66.0 ± 1.0 Ma prior to the India–Asia collision, whereas grains with metamorphic characteristics indicate later growth at 30.6 ± 0.4 Ma. Monazite grains from a cross‐cutting post‐kinematic dyke within the core of the DNCV antiform provide a minimum age constraint of 25.2 ± 1.4 Ma for the termination of fabric development. A separate and significant episode of monazite growth at c. 83–69 Ma is suggested to be the result of fluid‐assisted recrystallization following the emplacement of magmatic units.  相似文献   

Snow cover and snowmelt of an extensive High Arctic wetland: spatial and temporal seasonal patterns     

Jane Assini 《水文科学杂志》2013,58(4):738-755

Abstract

This study examined the end-of-winter snow storage, its distribution and the spatial and temporal melt patterns of a large, low gradient wetland at Polar Bear Pass, Bathurst Island, Nunavut, Canada. The project utilized a combination of field observations and a physically-based snowmelt model. Topography and wind were the major controls on snow distribution in the region, and snow was routinely scoured from the hilltop regions and deposited into hillslopes and valleys. Timing and duration of snowmelt at Polar Bear Pass were similar in 2008 and 2009. The snowmelt was initiated by an increase in air temperature and net radiation receipt. Inter-annual variability in spatial snowmelt patterns was evident at Polar Bear Pass and was attributed to a non-uniform snow cover distribution and local microclimate conditions. In situ field studies and modelling remain important in High Arctic regions for assessing wetland water budgets and runoff, in addition to model parameterization and validation of satellite imagery.

Editor Z.W. Kundzewicz

Citation Assini, J. and Young, K.L., 2012. Snow cover and snowmelt of an extensive High Arctic wetland: spatial and temporal seasonal patterns. Hydrological Sciences Journal, 57 (4), 738–755.  相似文献   

Potential Influence of a Developing La Niña on the Sea-Ice Reduction in the Barents–Kara Seas     

Rui Luo  Peng Zhang  Juan Dou 《大气与海洋》2013,51(3):182-194

Abstract

The potential influence of a developing La Niña on Arctic sea-ice annual variability is investigated using both observational data and an atmospheric general circulation model. It is found that during the developing phase of an eastern Pacific (EP) La Niña event in June, July, and August (JJA) and September, October, and November (SON), the sea-ice concentration (SIC) over the Barents–Kara Seas declines more than 15%. The local atmospheric circulation pattern associated with the EP La Niña is characterized as a weak decrease in geopotential height over the Barents–Kara Seas, combined with an anticyclone in the North Atlantic. The corresponding southerly winds push warm waters northward into the key sea-ice reduction region and directly accelerate sea-ice melt. Meanwhile, the abundant moisture contained in the lower troposphere is transported into the Arctic region by winds resulting from the local barotropic structure. The humid atmosphere contributes to both net shortwave and longwave radiation and thus indirectly accelerates the decline in sea ice. Simulations by the European Centre Hamburg Model, version 5.4, are forced by observed sea surface temperature anomalies associated with EP La Niña events. The results of the simulations capture the North Atlantic anticyclone and reproduce the moisture transport, which supports the premise that an EP La Niña plays a crucial role in sea-ice reduction over the Barents–Kara sector from the perspective of atmospheric circulation and net surface heat flux.  相似文献   

Mosaic desert pavement influences water infiltration and vegetation distribution on fluvial fan surfaces     

Dejin Wang  Wenzhi Zhao  Hong Zhou  Weicheng Luo  Hu Liu 《水文研究》2021,35(9):e14373

Desert pavements (DPs) are critical for maintaining ecological stability and promoting near-surface hydrological cycling in arid regions. However, few studies have focused on eco-hydrological processes of DPs in the ecological systems of fluvial fans. Although DP surfaces appear to be barren and flat, we found that the surfaces are characterized by surface mosaic patterns of desert pavement (mosaic DP) and bare ground (mosaic BG). We investigated the effects of mosaic DP on water infiltration and vegetation distribution at six sites in fluvial fans (one on a hillside and five within the sectors of fans) along a southwest belt transect in northern Linze County, in the central Hexi Corridor (China). We found significant differences in mosaic DP between the hillside and sector sites in terms of pavement thickness and vesicular horizon thickness (Av thickness), particle composition, and bulk density, although significant differences were absent for mass soil water content, gravel coverage, and surface gravel size. The mosaic DP inhibited water infiltration by the pavement layer, where the sorptivity (S), initial infiltration rate (i_int), steady-state infiltration rate (i_sat) and infiltration time (T) averaged 1.19 cm/min^-0.5, 0.64 cm/min, 0.13 cm/min and 12.76 min, respectively. Where the pavement layer was scalped, the S, i_int, and i_sat increased by 0.27 cm/min^-0.5, 0.52 cm/min, and 0.40 cm/min, respectively, and the T reduced by 7.42 min. Water infiltration was mainly controlled by the pavement layer thickness (+), Av thickness (−), surface gravel coverage (−), fine earth (+) and fine gravel (−) in the pavement layer. The DP surfaces only had a sparse covering of shrubs, but an abundance of herbs. Few shrubs were present on the mosaic DP, but a greater number of shrubs and herbs grew on the mosaic BG. It can be concluded that DPs can maintain vegetation stability for different surface mosaic patterns. This study deepens our understanding of the eco-hydrological cycle of DP landscapes in arid regions.  相似文献   

A comparison of fibre‐optic distributed temperature sensing to traditional methods of evaluating groundwater inflow to streams     

Martin A. Briggs  Laura K. Lautz  Jeffrey M. McKenzie 《水文研究》2012,26(9):1277-1290

There are several methods for determining the spatial distribution and magnitude of groundwater inputs to streams. We compared the results of conventional methods [dye dilution gauging, acoustic Doppler velocimeter (ADV) differential gauging, and geochemical end‐member mixing] to distributed temperature sensing (DTS) using a fibre‐optic cable installed along 900 m of Ninemile Creek in Syracuse, New York, USA, during low‐flow conditions (discharge of 1·4 m³ s^?1). With the exception of differential gauging, all methods identified a focused, contaminated groundwater inflow and produced similar groundwater discharge estimates for that point, with a mean of 66·8 l s^?1 between all methods although the precision of these estimates varied. ADV discharge measurement accuracy was reduced by non‐ideal conditions and failed to identify, much less quantify, the modest groundwater input, which was only 5% of total stream flow. These results indicate ambient tracers, such as heat and geochemical mixing, can yield spatially and quantitatively refined estimates of relatively modest groundwater inflow even in large rivers. DTS heat tracing, in particular, provided the finest spatial characterization of groundwater inflow, and may be more universally applicable than geochemical methods, for which a distinct and consistent groundwater end member may be more difficult to identify. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Garnet crystal plasticity in the continental crust,new example from south Madagascar     

J.‐E. MARTELAT  K. MALAMOUD  P. CORDIER  B. RANDRIANASOLO  K. SCHULMANN  J.‐M. LARDEAUX 《Journal of Metamorphic Geology》2012,30(4):435-452

Garnet (10 vol.%; pyrope contents 34–44 mol.%) hosted in quartzofeldspathic rocks within a large vertical shear zone of south Madagascar shows a strong grain‐size reduction (from a few cm to ～300 μm). Electron back‐scattered diffraction, transmission electron microscopy and scanning electron microscope imaging coupled with quantitative analysis of digitized images (PolyLX software) have been used in order to understand the deformation mechanisms associated with this grain‐size evolution. The garnet grain‐size reduction trend has been summarized in a typological evolution (from Type I to Type IV). Type I, the original porphyroblasts, form cm‐sized elongated grains that crystallized upon multiple nucleation and coalescence following biotite breakdown: biotite + sillimanite + quartz = garnet + alkali feldspar + rutile + melt. These large garnet grains contain quartz ribbons and sillimanite inclusions. Type I garnet is sheared along preferential planes (sillimanite layers, quartz ribbons and/or suitably oriented garnet crystallographic planes) producing highly elongated Type II garnet grains marked by a single crystallographic orientation. Further deformation leads to the development of a crystallographic misorientation, subgrains and new grains resulting in Type III garnet. Associated grain‐size reduction occurs via subgrain rotation recrystallization accompanied by fast diffusion‐assisted dislocation glide. This plastic deformation of garnet is associated with efficient recovery as shown by the very low dislocation densities (10¹⁰ m^?3 or lower). The rounded Type III garnet experiences rigid body rotation in fine‐grained matrix. In the highly deformed samples, the deformation mechanisms in garnet are grain‐size‐ and shape‐dependent: dislocation creep is dominant for the few large grains left (>1 mm; Type II garnet), rigid body rotation is typical for the smaller rounded grains (300 μm or less; Type III garnet) whereas diffusion creep may affect more elliptic garnet (Type IV garnet). The P–T conditions of garnet plasticity in the continental crust (≥950 °C; 11 kbar) have been identified using two‐feldspar thermometry and GASP conventional barometry. The garnet microstructural and deformation mechanisms evolution, coupled with grain‐size decrease in a fine‐grained steady‐state microstructure of quartz, alkali feldspar and plagioclase, suggests a separate mechanical evolution of garnet with respect to felsic minerals within the shear zone.  相似文献   

Rate and historical change of direct recharge from precipitation constrained by unsaturated zone profiles of chloride and oxygen‐18 in dry river bed of North China Plain     

Ruiqiang Yuan  Xianfang Song  Dongmei Han  Yinghua Zhang  Liang Zhang  Bing Zhang  Xiting Long  Yilei Yu 《水文研究》2012,26(9):1291-1301

Quantifying of direct recharge derived from precipitation is crucial for assessing sustainability of well‐irrigated agriculture. In the North China Plain, the land use is dominated by groundwater‐irrigated farmland where the direct recharge derived from precipitation and irrigation. To characterize the mean rate and historical variance of direct recharge derived from precipitation, unsaturated zone profiles of chloride and δ¹⁸O in the dry river bed of the Beiyishui River were employed. The results show that archival time scale of the profile covers the duration from 1980 to 2002 (corresponding to depths from 5 to 2 m) which is indicated by matching the δ¹⁸O peaks in the isotope profile with the aridity indexes gained by instrumental records of annual precipitation and annual potential evaporation. Using the chloride mass balance method, the mean rate of the direct recharge corresponding to the archival time scale is estimated to be 3·8 ± 0·8 mm year^?1, which accounts for about 0·7% of the long‐term average annual precipitation. Further, the direct recharge rates vary from 2·1 to 6·8 mm year^?1 since 1980. Despite the subhumid climate, the estimate of recharge rates is in line with other findings in semiarid regions. The low rate of direct recharge is considered as a result of the relative dry climate in recent decades. In dry river bed, unsaturated zone profiles of chloride and δ¹⁸O combined with instrumental records could offer valuable information about the direct recharge derived from precipitation during droughts. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献