Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

非等间距GM(1,1)模型的总体最小二乘算法及其病态问题

在非等间距GM(1,1)模型中，系数矩阵中有无误差的常数项和有误差的随机项，并且系数矩阵与观测向量误差同源，即系数矩阵与观测向量中有相同的元素存在，这些相同元素应该有相同的改正数，为此本文推导了一种适合非等间距GM(1,1)模型求解的总体最小二乘算法。同时，考虑到非等间距GM(1,1)模型中存在病态问题时影响总体最小二乘计算结果的稳定性，提出对系数矩阵常数列乘以某一常数的方法，以改善病态问题。     

内蒙古鲁新井田水文地质条件研究

矿产资源开发导致了地下水失衡,地下水失衡又给矿产开发造成了极大的安全隐患,近年来,矿产资源开发与地下水环境保护之间的矛盾愈发凸显。通过对内蒙古鲁新井田典型的水文地质条件进行分析研究,分析采矿导致地下水失衡机理,深入剖析矿井开采充水条件及矿山开发对地下水环境的影响,合理提出了促进矿产开发与地下水保护相互协调的对策建议,为实现"采矿保水"协调统一提供了基础地质依据。为类似地区矿产开发过程中遵循自然规律,趋利避害,保障生产安全,保护地下水环境安全,实现资源绿色开发有较好的指导借鉴作用。     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

基于灯光数据的京津冀城市多标度异速分析

城镇体系的时空演化是一个复杂的动力学过程,具有无尺度性,其特征可以采用异速标度指数定量描述。论文以京津冀为例,基于1992—2013年夜间灯光数据的校正结果,提取城镇夜间灯光信息,采用多标度异速分析方法,分析京津冀城镇体系的相对发展特征及其空间分异格局。研究发现：①京津冀大城市的增长空间局限日益显著。虽然大城市如京津唐的绝对发展水平高,但接近环境承载量的极限,22 a间增长快速的是大城市交接地带的小城市(如三河、迁安、廊坊等)。②相对增长速度快的城市从东北到西南呈现带状分布——政府选择发展的“雄安新区”就在这个轴带上。由此可以得出两个结论：其一,京津冀城市演化的趋势相对均衡化,空间结构异质性逐步减弱,整个城镇体系的发展趋于相对平衡状态;其二,大、小城市“两极”贯通或许是增强城镇体系功能的有效途径。京津冀的空间格局优化不能仅仅考虑主要的、规模大的中心城市,而应该基于等级结构协调整个城镇网络体系。     

自然资源一体化政务服务系统及数据融合建设——以南京为例

自然资源部制定的信息化建设总体方案中，要求构建“互联网+自然资源政务服务”应用体系，实现“服务事项标准统一、整体联动、业务协同，自然资源政务服务和共享开放能力全面提高”。一体化政务服务系统作为政务服务应用体系的重点建设内容，也是落实自然资源管理“两统一”职责的重要支撑，如何整合各类数据资源、发挥系统作用，成为亟待研究的问题。本文在分析自然资源业务管理需求后，提出了一体化政务服务系统及数据融合的建设思路。     

基于纹理特征的围填海SAR图像分水岭分割方法

本文提出了一种基于纹理特征的围填海SAR图像分水岭分割方法,首先对机载MiniSAR图像进行灰度共生矩阵纹理滤波,获得纹理特征图像,再对纹理特征图像进行分水岭算法分割,将获得的形态学重建图像进行门限阈值分割,得到最后的二值化分割结果。该方法一方面通过调整灰度共生矩阵纹理滤波的窗口大小,抑制了斑点噪声的影响;另一方面,利用分水岭算法对边缘模糊杂乱图像的优势,提高了围填海信息提取的准确性。实验结果表明,本方法对高分辨率SAR图像围填海监测图像的分割效果良好。     

海岸破波带内悬沙空间分布试验研究和解析表达式

关于海岸破波带内悬沙浓度水平和垂向分布的研究对于计算海岸输沙率和地形演变具有重要意义。本研究进行了规则波、波群和不规则波三种波浪情况破波带内悬沙浓度的水平和垂向分布的试验测量。试验在大尺度波浪水槽进行,接近实际海岸波况尺度。给出了破波带内多断面悬沙垂向分布的细致测量结果,并以此为基础给出了预报实际海岸破波带内悬沙浓度水平和垂向分布解析表达式,讨论了形成这些分布的物理原因和不同波况、不同破波带区域对分布的影响。     

300 BC-1900 AD无定河流域城镇时空格局演变

通过分析和整理从战国中晚期（约300 BC）至清末（约1900 AD）无定河流域历代县级及以上城址的位置、兴废年代数据,结合行政区划沿革、经济社会发展、政权更迭等资料,分析流域城镇格局的时空演变过程。研究表明：① 受气候周期性波动影响,无定河流域城镇的兴起与衰废具有明显的周期性特征,城镇几何中心的移动轨迹具有明显的“西北—东南”向潮汐性运动特征;② 城址存续年限普遍较短,具有明显的阶段性特征,流域内曾存在过4个阶段性中心城镇,中心城镇移动的方向和过程与城镇几何中心的潮汐性移动过程趋势一致;③ 流域城镇空间格局可分为3种类型,即秦汉与隋唐时期的沿河流谷地分布,宋、明两代的沿边境线与长城分布,以及元、清两代集中于流域下游分布。流域城镇未来的发展布局应重视气候变化对城镇分布的长期影响,关注区域环境的脆弱性,合理安排城镇体系发展规模与布局。