The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.     

A new approach to estimate canopy evaporation and canopy interception capacity from evapotranspiration and sap flow measurements during and following wetting

Canopy interception and its evaporation into the atmosphere during irrigation or a rainfall event are important in irrigation scheduling, but are challenging to estimate using conventional methods. This study introduces a new approach to estimate the canopy interception from measurements of actual total evapotranspiration (ET) using eddy covariance and estimation of the transpiration from measurements of sap flow. The measurements were conducted over a small‐scale sprinkler‐irrigated cotton field before, during and after sprinkler irrigation. Evaporation and sap flow dynamics during irrigation show that the total ET during irrigation increased significantly because of the evaporation of free intercepted water while transpiration was suppressed almost completely. The difference between actual ET and transpiration (sap flow) during and immediately following irrigation (post irrigation) represents the total canopy evaporation while the canopy interception capacity was calculated as the difference between actual ET and transpiration (sap flow) during drying (post irrigation) following cessation of the irrigation. The canopy evaporation of cotton canopy was calculated as 0.8 mm, and the interception capacity was estimated to be 0.31 mm of water. The measurement uncertainty in both the non‐dimensional ET and non‐dimensional sap flow was shown to be very low. Copyright © 2015 John Wiley & Sons, Ltd.     

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.     

运用ENVI实现SPOT-5卫星影像的RVI提取

运用ENVI软件处理长沙地区的SPOT-5影像,将全色影像和多波段影像进行融合,分析融合后的影像,再选择感兴趣区域进行影像裁剪;采用ENVI的波段运算(Band3/Band2)提取该地区比值植被指数并进行分析。结果表明:比值植被指数对植被覆盖率较高的区域有非常高的敏感度;在植被覆盖率小于50%的时候,其敏感度明显降低,在绿色植被覆盖区域的比值植被指数远大于1,在裸露地表、建筑物、水体等没有植被覆盖区域的比值植被指数则在1附近,而比值植被指数大于2的则出现在高覆盖健康绿色植被上。     

A note on “Origin of components in Chilean thermal waters” by

Risacher et al. (2011) have presented voluminous data on thermal waters of Chile that is quite appreciable; however, their work still falls short on several counts. The most notable shortcoming of the work is the presentation and treatment of data. The interpretations are based on questionable premises (viz., extent of seawater intrusion) and considerations (viz., using average chemical composition of rock types for geochemical modeling, considering volcanic rocks as the only rock type in contact with the geothermal fluids) and assumptions not always substantiated by facts (with proper references) could have been corroborated. Use of Cl/Br ratio for discrimination purpose is unconvincing, considering the uncertainty in the measurement of low (<1 mg/L) Br concentration.     

延边开山屯花岗闪长岩-石英闪长岩体年代学、地球化学特征及对古太平洋板块俯冲作用时限的制约

近年来古太平洋构造域的构造演化备受学者关注。本文报道的延边开山屯地区花岗闪长岩-石英闪长岩体LAICP-MS U-Pb年龄表明其形成时间为早侏罗世早期(198±1Ma),所采样品可根据Zr/Hf值分为高Zr/Hf值组花岗闪长岩和低Zr/Hf值组石英闪长岩。高Zr/Hf值组花岗闪长岩起源深度浅,富集Rb、Th、U、K等大离子亲石元素(LILEs),贫Nb、Ta、Ti等高场强元素(HFSEs),具壳源岩浆的特点。低Zr/Hf值组为壳源岩浆与来自深部的亏损地幔岩浆混合而成,岩石亏损Nb、Ta、Zr、Hf、Ti等高场强元素,具有典型的弧型岩浆地球化学特征。岩体中存在细粒闪长质包体,镜下可见针柱状磷灰石。开山屯岩体属钙碱性系列岩石,结合前人资料,认为其与该地区同时代火成岩组成北-东向分布的早侏罗世活动大陆边缘型火成岩带,而位于该带西侧的小兴安岭-张广才岭地区存在同时代弧后拉张带,两者构成典型的大陆弧与弧后拉张带模型,共同揭示了早侏罗世早期古太平洋板块对东北地区的俯冲作用。     

黄海南部小黄鱼洄游分布及种群结构特征

Monthly fishery survey data of the small yellow croaker Larimichthys polyactis in the southern Yellow Sea from2003 and 2013 were employed to evaluate the variation in the resource distribution and biological characteristics(especially body length and sex ratio) in the population on a decadal scale.The results indicated that the small yellow croaker migrated from the Shawai fishing ground to the Dasha fishing ground in spring and was mainly distributed in the central and western parts of the Dasha fishing ground in April and May.Larimichthys polyactis in the Dasha fishing ground migrated eastward to offshore wintering grounds in autumn and reached the central Dasha fishing ground in October and November.The small yellow croaker entered the western waters of the Shawai fishing ground in winter.A large number of age 0+ fish occurred in the Shawai and Jiangwai fishing grounds in October of 2003 and 2013.The body lengths of the spawning stock and wintering stock in 2013 were larger than those in 2003,and the monthly sex ratios(female to male) were significantly less than 1 in both years.The monthly distribution of this fish in the southern Yellow Sea was consistent with a previous finding that "the stock migrated between the wintering grounds in the west of Jeju Island and the Lüsi spawning grounds" but tended to move more northward,with the spawning grounds extending outward.In the past decade,body length variation experienced a decline after an increase,rather than a steady decrease.The sex ratio in the single-stick stow net showed a tendency to increase over the decade,but was either less than or more than 1 depending on the fishing gear;therefore,further studies should be conducted to determine the sex structure.     

苏里格大型致密砂岩气田储层结构与水平井提高采收率对策

在分析苏里格气田辫状河体系及储层结构特征的基础上,提出了剖面储量集中度的概念,建立了厚层块状型、垂向叠置泛连通型、分散局部连通型、分散孤立型4种砂体分布模式,探讨了不同储层结构下的水平井采出程度,提出了水平井提高采收率技术对策。研究结果表明,辫状河沉积体系复合有效砂体由于“阻流带”的存在,直井动用不完善,水平井能克服“阻流带”的影响,提高层内储量动用程度;但由于砂体多层状分散分布,水平井开发会导致纵向含气层系储量动用不充分,影响层间采出程度。对于剖面储量集中度高的厚层块状型、垂向叠置泛连通型储层,采用水平井整体开发,Ⅰ+Ⅱ类井比例达70%以上,可显著提高储量动用程度和采收率。对于剖面储量集中度低的分散局部连通型储层,采用直井井网开发后进行甜点式优选水平井井位加密部署,可提高采收率10%以上。     

基于人工鱼群算法和SCA等效理论的井中横波速度预测方法

在建立页岩岩石物理模型的基础上,根据等效自相容近似(SCA)岩石物理模型,构建出岩石的纵波速度、横波速度与岩石密度、组分和孔隙度等的定量关系,得出使理论纵波速度和实际纵波速度最接近的孔隙纵横比,进而将该孔隙纵横比作为约束条件来实现横波速度预测。反演算法利用人工鱼群算法来计算最佳孔隙纵横比,并将预测的横波速度与实际测得的横波速度对比,证明了人工鱼群算法的有效性。     

确定砂土动剪切模量和阻尼比的方法对比

确定土体动剪切模量的常用方法有规范法、Kumar法和自相关函数法,确定相应阻尼比的方法有规范法、Das and Luo法、Kokusho法、Kumar法和互相关函数法,为了分析不同方法所产生差异,实现定量化对比分析,笔者以福建标准砂（粒径为0.5~1.0 mm）为研究对象,采用不排水的应力控制动三轴试验,探讨不同的确定土体动剪切模量和阻尼比方法的差异性,并给出了不同土体条件建议选用的方法。结果表明：1）3种方法确定动剪切模量的结果有一定的差异,随剪应变的增大结果的差异逐渐增大,有效围压对结果的差异性有所影响,当剪应变为4×10^-3,有效围压为100 kPa时,3种方法差异显著,相对误差最大接近20%;2）而5种方法确定阻尼比的结果差异显著,随着剪应变的增大,5种方法确定的阻尼比相对误差大体上均在迅速减小,只有规范法在有效围压为100 kPa时,其相对误差有较小的增大趋势;5种方法中,Kumar法确定的阻尼比最接近平均阻尼比,互相关函数法远高于平均阻尼比,Das and Luo法和Kokusho法确定的阻尼比基本一致但低于平均阻尼比。建议以后的工程应用中,加载方式为应力控制时,可采用自相关函数法确定动剪切模量,采用Kumar法确定阻尼比,二者确定的动剪切模量和阻尼比均最接近平均值。