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法确定阻尼比,二者确定的动剪切模量和阻尼比均最接近平均值。     

Landslide susceptibility zonation method based on C5.0 decision tree and K-means cluster algorithms to improve the efficiency of risk management

Machine learning algorithms are an important measure with which to perform landslide susceptibility assessments,but most studies use GIS-based classification methods to conduct susceptibility zonation.This study presents a machine learning approach based on the C5.0 decision tree(DT)model and the K-means cluster algorithm to produce a regional landslide susceptibility map.Yanchang County,a typical landslide-prone area located in northwestern China,was taken as the area of interest to introduce the proposed application procedure.A landslide inventory containing 82 landslides was prepared and subse-quently randomly partitioned into two subsets:training data(70％landslide pixels)and validation data(30％landslide pixels).Fourteen landslide influencing factors were considered in the input dataset and were used to calculate the landslide occurrence probability based on the C5.0 decision tree model.Susceptibility zonation was implemented according to the cut-off values calculated by the K-means clus-ter algorithm.The validation results of the model performance analysis showed that the AUC(area under the receiver operating characteristic(ROC)curve)of the proposed model was the highest,reaching 0.88,compared with traditional models(support vector machine(SVM)=0.85,Bayesian network(BN)=0.81,frequency ratio(FR)=0.75,weight of evidence(WOE)=0.76).The landslide frequency ratio and fre-quency density of the high susceptibility zones were 6.76/km2 and 0.88/km2,respectively,which were much higher than those of the low susceptibility zones.The top 20％interval of landslide occurrence probability contained 89％of the historical landslides but only accounted for 10.3％of the total area.Our results indicate that the distribution of high susceptibility zones was more focused without contain-ing more"stable"pixels.Therefore,the obtained susceptibility map is suitable for application to landslide risk management practices.