高黏粒含量新吹填淤泥加固新技术室内研发

新吹填淤泥处于“稀泥汤”状态,工程特性极差,几乎无承载力,现有真空预压技术加固效果不理想,地基有效加固深度小,土体强度增长有限。为此,提出了“真空预压联合化学加固”的新思路。室内采用3种环保型化学加固剂设计了3种真空预压联合化学加固方案（VPCT1、VPCT2、VPCT3）和1种纯化学加固方案（CT）,与现行真空预压加固方案（VPT）进行了对比试验研究。研究结果表明：(1) 添加有效的化学加固剂后,新吹填淤泥在抽真空前较短的静置时间内就基本完成了自重沉积过程;(2) 化学加固剂选取得当能有效提高土体的加固效果;(3) 宏观上,VPCT3方案加固效果较均匀,加固后土体中黏粒含量明显降低、湿密度最大,孔隙比降低幅度明显,无侧限抗压强度也最大,为66.7 kPa（为VPT方案的3.5倍）;(4) 微观上,VPCT3方案加固后,土颗粒之间的孔隙总面积、孔隙平均周长、总孔隙数、孔隙比、孔隙率等最小,单元体之间有良好的定向度、排序性最好,连结方式以面-面为主导,相互之间接触紧密,团聚现象最明显。因此,相对于现行的真空预压技术而言,VPCT3方案土体加固效果非常明显,也较均匀,可以进一步开展现场试验研究,探讨其付诸于工程实践的可行性。     

基于Bootstrap方法的岩土体参数联合分布模型识别

小样本容量岩土体参数最优联合概率分布模型的识别是一个富有挑战性的问题。基于Bootstrap提出了小样本容量岩土体参数最优边缘分布函数和最优Copula函数识别方法。简要介绍了岩土体参数联合概率分布函数构造的Copula方法,采用AIC准则识别最优的边缘分布函数和Copula函数。将识别结果表示为不同备选边缘分布函数和Copula函数为最优边缘分布和最优Copula的权重系数集合,以基桩荷载-位移双曲线参数试验数据为例证明了所提方法的有效性。结果表明：基于小样本容量岩土体参数试验数据估计的样本均值、标准差和相关系数具有较大的离散性,这种离散性进一步导致了统计量AIC值存在较大变异性。提出的基于Bootstrap的最优边缘分布函数和最优Copula函数识别方法不仅可以有效地考虑统计量AIC值的变异性,而且能够综合地反映不同备选边缘分布函数和Copula函数为最优边缘分布和最优Copula函数的概率,为小样本容量岩土体参数最优边缘分布函数和最优Copula函数的识别提供了一条有效的途径。     

距离法定量评价储层地质模型的不确定性

针对目前常用的储层地质模型不确定性评价方法存在的主要问题,提出了先用距离函数计算模型之间的差异,再以差异的大小来判断不确定性大小的方法.以WZ油田西区为例,采用相控物性参数建模技术,利用顺序高斯模拟方法建立渗透率的三维模型.对各种度量差异的距离函数进行对比研究和分析,结果显示曼哈顿距离函数和欧氏距离函数能较好地刻画模型之间的差异.选用欧氏距离函数计算模型之间的差异,其原理是先计算每两个模型之间相对应的每一网格节点的渗透率值差的平方和,然后取平方根,得到一个表征各模型之间差异的矩阵.根据该矩阵可得到各个模型之间的差异程度,差异越大,不确定性就越大.最后通过对比模型过井剖面图分析结果与距离矩阵分析结果,说明了本方法的正确性,结果显示该方法能有效评价随机模拟生成的储层地质模型的不确定性.     

南北地震带壳幔结构成像研究

<正>对地壳和上地幔的速度结构进行成像是研究地球深部结构及动力学问题的一个有效途径。南北地震带地区是青藏高原与相邻块体的过渡带,这里的地质构造十分复杂,地震活动频繁。在南北地震带地区进行壳幔结构成像研究,有利于提高对青藏高原及其与相邻块体相     

利用TanDEM-X生成DEM的精度评定

目前,许多学者对TanDEM-X生成DEM开展了一些研究,其研究成果也显示了TanDEM-X生成高精度DEM的可行性。为了验证TanDEM-X/TerraSAR-X干涉生成的DEM能否满足测图要求,需要对其进行精度评价和分析。相对于C波段的ERS、ASAR和L波段的ALOS,X波段的高分辨率TerraSAR影像干涉条纹更密集,解缠更加困难。针对这一问题,本文设计了一种低分辨率SRTM辅助高分辨率的X波段的TerraSAR干涉相位解缠方案,提高了解缠的效率和精度。同时,本文提出了一种基于协方差函数的方法对TDX/TSX DEM进行精度分析和评价。该方法通过对各个距离上的协方差值进行拟合,消除了高程误差异常对InSAR DEM精度评价的影响,可以更加客观真实地反映DEM的精度。实例分析结果表明：采用协方差函数方法来评价DEM的精度是可行的,对于试验研究区域,TDX/TSX干涉生成的DEM总体精度为1.42 m,能够满足1：10 000测图要求,为我国空白地区的测图提供了有利条件。     

亚塑性模型不同积分算法的数值实现

亚塑性模型是以Jaumann应力率张量及变形率张量描述的一种率型本构关系,本构关系在非线性有限元分析计算中具有关键作用,解决用应变增量求解应力响应的问题需要一个时间积分过程。针对亚塑性本构模型发展了自适应隐式和显式两种不同的积分算法,给出了误差控制的方法,同时推导了自适应隐式积分算法所需的一致切向模量,并采用了两个不同的单元,利用ABAQUS平台比较了两种积分算法的数值模拟结果。为了实现从ABAQUS/Standard到ABAQUS/Explicit的过渡,开发了UMAT-VUMAT接口,从而可以使已有的UMAT子程序用于大变形动力问题分析。算例分析证明了研究结果的正确性。     

基于子区间法的隐式功能函数非概率可靠性方法研究

从统计资料不足的角度出发,研究一种岩土工程隐式功能函数非概率可靠性方法,作为概率可靠性方法的有益补充。首先利用系列区间扩展函数降低区间相关性,建立基于子区间法的隐式函数区间解的迭代程序;然后采用试验设计的方法,从全面设计的角度建立统计优化模型搜索最优解,达到一个目标区间内获得一个简单计算方案的目的;最后通过对功能函数中应力函数的区间解设置置信区间,形成有限计算代价下的偏向于绝对安全的隐式功能函数非概率可靠性模型。通过岩土工程中最小支护阻力函数演示该方法的具体操作程序。工程算例结果表明,由此得到的最小支护阻力上限值与蒙特卡洛模拟法最大值的相对误差为0.83%,偏向于绝对安全。     

喀斯特山区不同岩性分区下的景观生态质量差异——以罗甸县为例 

以贵州省罗甸县为例,获取岩性及2000、2005、2010和2015年4期景观格局分布等数据,从景观稳定程度和干扰程度视角对比分析了不同岩性分区的人为活动对景观生态质量变化的影响模式。研究表明:（1）利用景观稳定程度和干扰程度进行评价,能较好地反映不同岩性分区的景观生态质量差异;（2）非喀斯特区域的景观稳定程度高于亚喀斯特区域,典型喀斯特区域的景观稳定程度最低;（3）亚喀斯特区域的景观干扰程度高于非喀斯特区域,典型喀斯特区域受到的干扰最小;（4）非喀斯特区域的景观生态质量最好,亚喀斯特区域次之,典型喀斯特区域的景观生态质量较差;（5）2000—2015年,不同岩性分区下的景观生态质量变化呈现出不同趋势。典型喀斯特区域4期的景观生态质量分别为0.30、0.18、-0.24、-0.18,2000—2010年景观生态质量有所下降,但2010—2015年景观生态质量处于上升状态;亚喀斯特区域4期的景观生态质量分别0.55、0.17、0.32、-0.51,2000-2005年景观生态质量下降,而2005—2010年景观生态质量上升,2010-2015年景观生态质量又处于下降状态;非喀斯特区域4期的景观生态质量分别为0.66、0.38、0.27、-0.30,2000-2015年景观生态质量一直处于下降状态。受岩性及人类活动影响的亚喀斯特区域景观生态质量更为敏感。因此,在生态环境治理措施中应根据岩性的差异进行分区对待,更需关注亚喀斯特区域的景观生态质量。      

基于主(微)量元素分析法的细粒沉积岩沉积环境判别——以大民屯凹陷安福屯地区沈352井沙四下亚段为例

为明确大民屯凹陷安福屯地区沈352井沙四下亚段细粒沉积岩的沉积环境及演化特征,利用地球化学方法,根据岩石中w(V)/w(V+Ni)判断氧化还原环境的变化,根据w(Sr)/w(Ba)、w(V/)w(Ni)、w(Mg)/w(Ca)判断沉积环境中盐度的变化,根据w(Zr)/w(Al)及w(Al2O3)、w(K2O)、w(TiO2)、w(P2O5)、w(V)判断水体深度的变化,根据w(Sr)/w(Cu)、c(为几种元素的比值)判断沉积环境的古气候变化,分析研究区沉积环境特征及变化规律。结果表明:处于陆相沉积环境的沈352井沙四下亚段由深到浅分为三部分,即在3 280~3 352m处为还原—弱还原环境,盐度为半咸水—咸水,水体较深,属于半潮湿—半干燥气候;在3 248~3 280m处为弱还原环境,盐度为咸水,水体较浅,属于干热气候;在3 169~3 248m处为还原环境,盐度为咸水—半咸水,水体较深,属于温湿气候;其中粒屑发育段(3 328~3 333m)较为特殊,为氧化环境,盐度为微咸水,水体较浅,属于干热气候。该研究对分析研究区细粒沉积岩形成时的沉积环境及勘探潜力具有一定参考意义。     

基于粘弹性球体地球模型的震后位移与重力变化计算软件

在考虑地球曲率、成层结构、可压缩性与自重的前提下,Tanaka等［1-2］提出一套较为完备的粘弹性球体位错理论,可计算全球任意位置由地震产生的同震与震后形变（含位移、重力变化、大地水准面变化）。Gao等［3］给出了与上述理论相匹配、界面友好的计算软件,能计算30个震后时间点对应的震后形变。本研究针对Gao等的软件进行改进,可计算震后任意时间点对应的震后形变。新软件由3个部分组成：1）与32个震后时间点相对应的32套离散格林函数数值框架;2）格林函数插值计算程序,可针对上述32套格林函数数值框架进行插值运算,输出任意震后时间点对应的格林函数数值结果;3）积分计算程序,调用上述格林函数数值结果,计算任意类型地震在地表任意位置产生的同震与震后形变。一般情况下,使用者只需按要求准备辅助文件,提供发震断层模型和观测站位置信息,以及震中周围地区地幔粘滞性因子,先后运行格林函数插值计算程序和积分计算程序,即可计算出目标地震在地表任意位置产生的同震与震后形变。本文基于粘弹球体位错理论与弹性球体位错理论,分别计算2011年日本MW9.0地震引起的远场同震位移,2套结果的高度一致性证明了新程序的正确性。最后,介绍需要注意的若干事项,便于使用者掌握该软件。