A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

基于ROMS和4DVAR的沿轨与网格化SSH数据同化效果评价

Remote sensing products are significant in the data assimilation of an ocean model. Considering the resolution and space coverage of different remote sensing data, two types of sea surface height(SSH) product are employed in the assimilation, including the gridded products from AVISO and the original along-track observations used in the generation. To explore their impact on the assimilation results, an experiment focus on the South China Sea(SCS) is conducted based on the Regional Ocean Modeling System(ROMS) and the four-dimensional variational data assimilation(4 DVAR) technology. The comparison with EN4 data set and Argo profile indicates that, the along-track SSH assimilation result presents to be more accurate than the gridded SSH assimilation, because some noises may have been introduced in the merging process. Moreover, the mesoscale eddy detection capability of the assimilation results is analyzed by a vector geometry–based algorithm. It is verified that, the assimilation of the gridded SSH shows superiority in describing the eddy's characteristics, since the complete structure of the ocean surface has been reconstructed by the original data merging.     

南海风暴潮对风场和外围海水敏感性研究

为探索来自西北太平洋台风风暴潮与南海局地生成台风风暴潮不同,本研究在假设两种台风气象条件相同情况下,研究随台风而来外围海水所形成增水对南海沿岸的影响。以0814"黑格比"强台风风暴潮为基础,使用ROMS(regional ocean modeling system)模式进行数值模拟并通过设计对比试验方法进行研究,研究发现在台风登陆时引起的增水最大,最大增水出现在台风路径右侧,其中在沿岸区域,外围海水形成增水约占总增水10%,且大约3 h后出现增水回震现象。同时,设计对比试验,研究来自西北太平洋台风风暴潮对台风路径、台风强度、台风移动速度和流入角等气象条件敏感性,并获得与前人一致的结果。     

基于“DEM-NDVI-土地覆盖分类”的天山博格达自然遗产地山地垂直带提取与变化分析

本文基于Landsat影像数据获取天山博格达自然遗产地土地覆盖分类,结合归一化植被指数（NDVI）和数字高程模型（DEM）构建“DEM-NDVI-土地覆盖分类”散点图分析研究区植被受海拔和坡向的水热空间变化影响的分布特征,通过概率统计分析提取博格达遗产地山地垂直带,并结合研究区的气温、降水数据和NDVI变化特征分析垂直带变化的原因。研究结果表明：① 本文利用“DEM-NDVI-土地覆盖分类”散点图,揭示了研究区1989年和2016年的NDVI值和分类类别随着海拔上升的变化特征,其中NDVI值随着海拔上升呈现“倒U形”变化,而不同分类类别在一定的海拔区间内呈现出聚集效应,且不同分类类别有明显的高程界限。② 1989年和2016年博格达遗产地山地垂直带分带上限分别为：1278 m和1185 m（温带荒漠草原带）、1784 m和1759 m（山地草原带）、2706 m和2730 m（山地针叶林带）、3272 m和3293 m（高山草甸带）、3636 m和3690 m（高山垫状植被带）。③ 博格达遗产地1989年和2016年山地垂直带受区域气温升高和降雨增加的影响有较为明显的改变,其中温带荒漠草原带最为敏感,其上限变化最大,向下收缩93 m;山地针叶林带的分布范围则向两侧扩张49 m;山地草甸带带宽基本保持不变,但整体上移了约20 m;冰雪带则受到全球气候变暖的影响向上退缩54 m。     

面向稀疏散布数据集的时空Kriging优化

时空Kriging法通过将变异函数向时空域进行扩展得到时空变异函数,有效地利用时空邻近的采样点综合进行插值,由于时空稀疏散布数据集具有单一时刻下样本点数量少以及时空分布不规律的特点,难以满足使用时空Kriging插值法的基本条件,导致插值精度不高,据此本文提出了优化方法：通过多时段叠置拟合空间变异函数的方法,综合利用时空邻域内的采样点以解决单一时刻下空间邻域内数量不足情况;控制时间变异对空间变异函数拟合的误差影响;采用积合式模型构建时空变异函数进行插值。最后使用Argo海温数据进行插值实验,在相同条件下与时空Kriging法以及时空权重法的交叉验证结果对比得出,该方法在保证拟合所需采样点数量要求的同时,有效削减了一般时空Kriging法中时间变异对空间变异函数拟合结果的干扰,插值结果的绝对误差均值从0.5降低至0.2以内,稳定性进一步增强,改善了时空Kriging法在稀疏散布数据条件下精度上的不足。     

PRELIMINARY STUDY ON THE FLOODING AND DROUGHT CLAMITY DURING PAST 1500 YEARS IN THE HAI′AN REGION，JIANGSU PROVINCE

Collection and arrangement of the historical records of climatic changes and environment evolution,espectial-lyin the aspect of calamities,are made on the history documents of past 1500 years about Haiˊan region,Jiangsu Province.There existed two obvious flooding-drought frequently-occurring periods:one was from 1550 AD to 1850 AD and another was 100 AD to 1200AD.The period of 1550 AD to1850 AD is interrupted by two relatively arid and cold climatic periods:one was from 1630 AD to 1700 AD and another was 1750 AD to 1820 AD.The main characteristic of the calamity periods is that they occurred by turns,and sometimes,both drought and flooding occurred in the same year.The instability of the climatic changes in the Little Ice Age may be the main reason of the frequently-occurring flooding and drought in Haiˊan region.Research results also show that the frequently-occurring periods of flooding and drought is in close relationship with the solar activity,and therefore,occurrence of the flooding and drought may be in relation with the intensity of the solar activity.This hypothesis may need further study in the future.     

青藏高原隆升的非线性动态有限元仿真研究

根据青藏高原的地质特征建立分析模型，采用3维动态有限元方法，在计算仿真板块速度场的基础上，计算在青藏高原的隆升过程中该地区地壳岩石的等效应力和位移随时间的变化，计算仿真得到的速度场与1998年GPS观测的速度场吻合良好；与过去一贯的假设相反，计算结果反映出地壳应力场不是静态的，而是此起彼伏，不断变化的，应力值最大且变化最剧烈的地区在克什米尔地区、鄂尔多斯地区和鲜水河－小江断裂带，与地震多发区域吻合。     

大同—阳高Ms6.1地震前兆综合加权信息熵的变化

利用“地震预报计算机专家系统”的思想对大同 -阳高 Ms6 .1地震前每一前兆异常事件进行综合评估 ,以每一异常的最可能发震时间来计算发震概率 ,利用地震前兆综合加权信息熵研究了系统熵值与地震的关系。对华北地区的地震前兆综合加权信息熵研究表明 ,在大同 -阳高 Ms6 .1地震前 ,信息熵出现了明显的减熵有序变化     

对青藏东北缘现今块体划分、运动及变形的初步研究

利用2维非连续变形分析方法（DDA），以位移代替围压作为边界约束力，研究青藏东北缘现今块体划分及其运动变形。根据该地区地质构造及地震活动，以GPS点测量位移作为模拟结果约束点，得出了较合理的块体划分模型和随时间演化的主应变分布图，并把应变高值区与近几年来发生的5级以上地震作对比，得出了研究区内地震危险性可能较大的区域。另外，对模拟的甘青块体与阿拉善块体的边缘带断裂左旋运动做了大概计算。     

抽油机悬点投影测量仪

油田抽油机的抽油杆幌动幅度过大，是引起抽油杆断裂的一个重要因素，检测这种幌动幅度是防止抽油杆断裂的一种有效手段。设计了一种基于面阵CCD和普通半导体激光器（LD）测量这种幌动的悬点投影测量方法，通过数字卷积滤波，达到了范围为0－40mm，误差＜0.2mm的技术检测指标。