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.     

Revised runoff curve number for runoff prediction in the Loess Plateau of China

The soil conservation service (now Natural Resources Conservation Service) Curve Number (SCS-CN), one of the most commonly used methods for surface runoff prediction. The runoff calculated by this method was very sensitive to CN values. In this study, CN values were calculated by both arithmetic mean (CN_C) and least square fit method (CN_F) using observed rainfall-runoff data from 43 sites in the Loess Plateau region, which are considerably different from the CN₂ values obtained from the USDA-SCS handbook table (CN_T). The results showed that using CN_C instead of CN_T for each watershed produce little improvement, while replacing CN_T with CN_F improves the performance of the original SCS-CN method, but still performs poorly in most study sites. This is mainly due to the SCS-CN method using a constant CN value and discounting of the temporal variation in rainfall-runoff process. Therefore, three factors—soil moisture, rainfall depth and intensity—affecting the surface runoff variability are considered to reflect the variation of CN in each watershed, and a new CN value was developed. The reliability of the proposed method was tested with data from 38 watersheds, and then applied to the remaining five typical watersheds using the optimized parameters. The results indicated that the proposed method, which boosted the model efficiencies to 81.83% and 74.23% during calibration and validation cases, respectively, performed better than the original SCS-CN and the Shi and Wang (2020b) method, a modified SCS-CN method based on tabulated CN value. Thus, the proposed method incorporating the influence of the temporal variability of soil moisture, rainfall depth, and intensity factors suggests an accurate runoff prediction for general applications under different hydrological and climatic conditions on the Loess Plateau region.     

Paleogene of the Huanghekou Sag in the Bohai Bay Basin,NE China: deposition–erosion response to a slope break system of rift lacustrine basins

The sequence architecture and depositional systems of the Paleogene lacustrine rift succession in the Huanghekou Sag, Bohai Bay Basin, NE China were investigated based on seismic profiles, combined with well log and core data. Four second‐order or composite sequences and seven third‐order sequences were identified. The depositional systems identified in the basin include: fan delta, braid delta, meander fluvial delta, lacustrine and sublacustrine fan. Identification of the slope break was conducted combining the interpretation of faults of each sequence and the identification of syndepositional faults, based on the subdivision of sequence stratigraphy and analysis of depositional systems. Multiple geomorphologic units were recognized in the Paleogene of the Huanghekou Sag including faults, flexures, depositional slope break belts, ditch‐valleys and sub‐uplifts in the central sag. Using genetic division principles and taking into consideration tectonic features of the Paleogene of the Huanghekou Sag, the study area was divided into the Northern Steep Slope/Fault Slope Break System, the Southern Gentle Slope Break System and T10 Tectonic Slope Break System/T10 Tectonic Belt. Responses of slope break systems to deposition–erosion are shown as: (1) basin marginal slope break is the boundary of the eroded area and provenance area; (2) ditch‐valley formed by different kinds of slope break belts is a good transport bypass for source materials; (3) shape of the slope break belt of the slope break system controls sediments types; (4) the ditch‐valley and sub‐sag of a slope break system is an unloading area for sediments; and (5) due to their different origins, association characteristics and developing patterns, the Paleogene slope break belt systems in the Huanghekou Sag show different controls on depositional systems. The Northern Fault Slope Break system controls the deposition of a fan delta‐lacustrine‐subaqueous fan, the Southern Gentle Slope Break system controls the deposition of a fluvial–deltaic–shallow lacustrine and sublacustrine fan, and the T10 Tectonic Slope Break System controls the deposition of shallow lacustrine beach bar sandbodies. The existence of a slope break system is a necessary but not a sufficient condition for studying sandbody development. The formation of effective sandbodies along the slope break depends on the reasonable coupling of effective provenance, necessary association patterns of slope break belt, adequate unloading space and creation of definite accommodation space. Copyright © 2013 John Wiley & Sons, Ltd.     

顾及形态特征和语义信息的双线道路识别方法

双线道路识别与提取是城市路网综合的关键。提出了一种城市双线道路提取的方法。在构建好的道路网眼的基础上,综合考虑双线道路网眼形态特征以及构成网眼路段之间的语义相似度,设计了识别双线道路网眼的综合指标。通过判断网眼综合指标确定最终的双线道路网眼,最后根据识别出的双线道路网眼提取出双线道路。实验表明,该方法能有效提取出城市道路数据中的双线道路。     

X射线衍射技术应用于宝石鉴定-合成及晶体结构研究进展

因宝石检测常附带无损要求、宝石评价具科学性和不确定性等特点,X射线衍射(XRD)在宝石学这门矿物学的分支学科中的应用具有特殊性,其应用领域主要集中在三个方面。①宝石鉴定领域:对单晶宝石和多晶宝石进行物相鉴定、宝石产地特征研究以及宝石矿物多型的种类划分等。例如根据XRD图谱,可将品种繁多的微晶石英隐晶质变种玉髓中不同的SiO₂质矿物相精确区分;通过黏土矿物成分及其含量可将鸡血石中"地"的种类进行划分。②宝石矿物的晶体结构研究:对成分复杂的单晶宝石进行晶体化学分析、获取矿物类质同象替代情况、利用结晶度划分宝石品质。例如晶格参数中的c/a比值可将绿柱石的Al八面体类质同象替代和Be四面体类质同象替代进行区分;软玉中通过XRD数据计算出的结晶度与软玉的品质关系密切。③宝石学研究技术的开发:对宝石的优化处理过程进行监测,在宝石合成过程中检验矿物晶形大小、结晶质量及内应力等,以有机宝石的晶体结构为基础的人工养殖技术研究。例如尖晶石在热处理过程中,Mg离子由四面体中的T位迁移至八面体中的M位,导致M-O和T-O键长随温度变化,并反映在尖晶石的色调中;在托帕石合成过程中,XRD数据显示770~800℃时形成含氟托帕石(与天然托帕石结构相同),而在950℃和1000℃出现含氟托帕石离解成刚玉和莫来石。目前,我国在宝石学的XRD应用方面起步较晚,国外在宝石鉴定、优化与合成方面的应用相对成熟。本文认为,该技术与扫描电镜、质子探针等分析技术的联用仍有很大发展空间。     

Biogenic volatile compound emissions from a temperate forest,China: model simulation

Emissions of biogenic volatile organic compounds (BVOC) were measured using a relaxed eddy accumulation (REA) technique on an above-canopy tower in a temperate forest (Changbai Mountain, Jilin province, China) during the 2010 and 2011 summer seasons. Solar global radiation and photosynthetically active radiation (PAR) were also measured. Based on PAR energy dynamic balance, an empirical BVOC emission and PAR transfer model was developed that includes the processes of BVOC emissions and PAR transfer above the canopy level, including PAR absorption and consumption, and scattering by gases, liquids, and particles (GLPs). Simulated emissions of isoprene and monoterpenes were in agreement with observations. The averages of the relative estimator biases for the flux were 39.3 % for isoprene, and 27.1 % for monoterpenes in the 2010 and 2011 growing seasons, with NMSE (normalized mean square error) values of 0.133 and 0.101, respectively. The observed and simulated mean diurnal variations of isoprene and monoterpenes in the 2010 and 2011 growing seasons were evaluated for the validation of the empirical model. Under observed atmospheric conditions, the sensitivity analysis showed that emissions of isoprene and monoterpenes were more sensitive to changes in PAR than to water vapor content or to the magnitude of the scattering factor. The emissions of isoprene and monoterpenes in the 2010 and 2011 growing seasons (from June to September) were estimated using this empirical model along with hourly observational data, with mean hourly emissions of 1.71 and 1.55 mg m^?2 h^?1 for isoprene, and 0.48 and 0.47 mg m^?2 h^?1 for monoterpenes in 2010 and 2011, respectively. As formaldehyde (HCHO) is considered as the main oxidation product of isoprene and monoterpenes, it is necessary to investigate the link between HCHO and BVOC emissions. GOME-2 HCHO vertical column densities (VCDs) can be used to estimate BVOC emission fluxes in the Changbai Mountain temperate forest.     

相对湿度和PM2.5浓度对大气能见度的影响研究:基于小时资料的多站对比分析

为了全面分析浙江省不同区域能见度变化基本特征及影响机理,基于杭州、宁波、温州3个国家基本气象站2013-2014年逐时能见度观测资料,比较分析了3市能见度变化的基本特征。发现3市不同等级能见度出现频率基本一致,随着能见度等级的提高,出现频率逐渐降低;从能见度的日变化来看,07时（北京时）前后最低,之后缓慢上升,14-15时达到最高,随后逐渐下降;全年有两个能见度较低时段,分别出现在12月-次年2月和5-6月;总体而言,宁波能见度最优,杭州和温州大体相当。功率谱分析结果表明,3市能见度均有显著的日周期,高频波段呈现出多个显著谱峰,低频波段存在若干显著谱峰。进一步开展机理分析,发现相对湿度和PM_2.5浓度是调制大气能见度的关键因子,相对湿度增大、PM_2.5浓度升高导致能见度降低。在同一相对湿度等级下,初始阶段能见度随PM_2.5浓度的升高迅速降低,到达“拐点”之后降低速率趋于缓慢。在同一PM_2.5浓度水平下,相对湿度越大,能见度越低,说明水汽对能见度也有重要影响。基于相对湿度和PM_2.5浓度两个因子,采用非线性拟合方案构建了大气能见度定量统计模型,总体而言模型拟合效果较好。最后针对研究中存在的不足和未来值得进一步发掘的科学问题进行了讨论。     

电涌保护器级间退耦元件的选型分析

针对电涌保护器（SPD）级间能量配合试验中配合失效的问题,主要利用波的传输理论与基尔霍夫定律进行分析,结合JB/T 10618-2006（组合式电涌保护器）的规定假设一套组合型SPD,通过基础理论分析与精确计算得出：当被保护装置属于低频（直流）设备时,在不影响信号强度的情况下,退耦元件采用电感或电阻型元件均可;当被保护装置属于高频设备时,退耦元件只能选用低电感量的电阻型退耦元件.退耦元件的参数值均可按照波的传输理论与基尔霍夫定律进行理论计算.这在组合型SPD退耦元件的选型过程中,具有重要的指导意义.