一种二阶锥线性化方法在上限有限元中的应用研究

对于平面应变条件下岩土稳定性分析,基于线性规划的上限有限元需对常用的摩尔-库仑屈服准则形成的二阶锥约束进行线性化,直接地处理方法是以外接多边形替代锥体投影形成的圆域。为了提高线性化精度往往需直接增加外接多边形边数,从而造成线性规划模型中决策变量包含大量的塑性乘子变量,使计算难度大为增加甚至变得不可行。为此,引入Ben-Tal和Nemirovsky提出的一种二阶锥线性化方法,并将其嵌入到自编的上限有限元程序。经算例分析发现,该法与外接多边形线性化方法所获计算结果相互印证,且其通过适量的增加决策变量和等式约束数目,能保证摩尔-库仑屈服准则线性化精度,同时形成的线性规划规模更小,可望应用于基于线性规划模型的上限有限元中。     

基于微米CT页岩微裂缝表征方法研究

页岩储层中的微裂缝是页岩气的储集空间,也是页岩气运移、散失的重要通道。通过微米CT扫描、场发射扫描电镜(FE- SEM)和岩芯观测等手段,明确了焦石坝地区JYA井龙马溪组页岩微裂缝结构特征,对比了微裂缝在不同条件下(地表或地层)的张开度及其对页岩渗透率的贡献,初步探讨了微裂缝对页岩纳米孔隙保存的影响。结果表明：① 基于微米CT可将页岩样品分为基质组分、高密度矿物(如黄铁矿)及裂缝三个部分,其体积比分别为97. 88%、0. 88%和1. 24%;② 微裂缝在空间上主要呈孤立状、细长条状或锯齿弯曲状,其长度分布范围为3. 4~299μm,以3. 4~20μm为主,平均为12. 44μm;③ 微裂缝在地层中的张开度约为2. 2~30. 33μm,平均为3. 97μm,相对于地表张开度(平均为6. 24μm)降低了36. 38%。尽管这部分微裂缝体积仅占页岩样品总体积的0. 7%,却提供了页岩储层中70%以上的渗透率。笔者认为JYA井龙马溪组页岩内部大量发育的网状结构微裂缝是页岩气沿着“孔隙—微裂缝—宏观裂缝”路径发生强烈“解吸—运移—逸散”作用的重要通道,进而导致孔隙因流体压力降低而压缩变形和闭合,不利于纳米孔隙的保存,这对于研究微裂缝在页岩储层中主要为储集空间或逸散通道具有重要的现实意义。     

支持向量机在储层厚度预测和计算中的应用

本文通过对油田储层结构的分析,运用支持向量机的理论和方法,建立了用于预测和计算储层厚度的支持向量机回归模型,并对该模型从参数变化范围、核函数选择、误差评价的标准等多方面进行了探讨,找出了建立储层厚度预测模型的一种有效方法,通过对实际储层厚度的预测,证明该方法在预测和计算储层厚度中具有较高的参考价值.     

基于Web服务的多源Argo数据共享与可视化应用研究

本文分析了当前海洋Argo观测数据在数据和模型共享与应用方面存在的不足,结合Web Service技术,提出在分布式网络环境下将海洋Argo数据与应用模型封装为Web服务的方式,实现多源海洋Argo数据集成与可视化应用的方法。基于该方法构建了"多源海洋Argo数据动态集成与远程可视化系统平台"。通过该平台的构建,为用户提供海洋Argo数据的透明的Web应用,实现了基于网络对多源海洋Argo数据的集成与交互式可视化表达。     

经纬仪横轴误差对竖直角观测的影响

通过对经纬仪横轴本身不水平和由于竖轴不竖直导致横轴不水平2方面分析横轴倾斜对竖直角观测的影响,推出相应的公式.并通过实验数据得出,通过水准管安置经纬仪,由于横轴倾斜角非常小,不必考虑横轴倾斜对竖直角观测的影响.     

影响地震活动性因素的非均匀细胞自动机模拟研究

基于细胞自动机模型,构建了由81×81个细胞单元组成非均匀二维单断层介质样本,研究了断层结构细观非均匀性以及相关模型参数变化对模拟输出结果和强度分布特征的影响。研究结果表明:结构非均匀性是影响地震活动性的控制性参量,随着非均匀程度的增强,地震事件数量增多,震级分布变得越来越均匀,破裂演化模式由相对"脆性"向"塑性"变化;介质结构的初始条件和应力加载的随机性对地震活动演化过程基本不产生影响,此外,应力降分配系数减小使b值主体线性段向大震级延伸;随着局部摩擦损耗系数增大,模拟地震事件的数量有所增加,大震级事件数量和震级衰减增强b,值主体线性区间变窄,使断层的变形破坏方式由相对"脆性"向"塑性"变化。这些对理解孕震过程的复杂性是有帮助的。     

新疆巴楚和塔什库尔干地震台下方地壳结构研究

利用位于新疆塔里木盆地的巴楚地震台（BCH）和西昆仑山区的塔什库尔干地震台（TAG）多年记录的高质量远震波形数据,应用接收函数H-k叠加方法研究了台站下方的地壳结构。研究结果显示,BCH和TAG地震台下方地壳厚度差异分明,区域地壳厚度与地形有很好的对应关系;两个台都具有高的地壳平均波速比Vp／Vs值。BCH台下方的地壳厚度为44km、地壳内的平均波速比为1．849,在该台地壳中部21km处存在清晰的间断面,该间断面内存在低的平均P波波速和作／K值,该间断面的深度与邻近巴楚台的伽师震源区的精确定位的震源深度下界面一致,可能为地壳内的脆．韧转换带。TAG台下方的地壳厚度为69km,地壳内的平均波速比为1．847,较厚的地壳和高波速比可能表明该台下地壳的物质易于发生塑性流动,为地壳的变形和增厚创造了条件。     

基于线性规划模型的极限分析上限有限元的实现

极限分析上限有限元法常利用三角形常应变率单元和摩尔-库仑屈服函数线性化的方法,以形成较易求解的线性规划模型。然而为满足计算精度,需引入大量的优化变量,增加了计算和存储的难度。为此,基于线性规划模型,利用MATLAB编制极限分析上限有限元程序,针对线性规划模型等式约束矩阵的高度稀疏性的特点,以稀疏矩阵的方式存储,从而在一定程度上解决了上述问题,使得上限有限元法能处理较大规模的岩土工程稳定性问题。以条形基础地基承载力课题为例进行算例分析,验证该方法的有效性,同时讨论了模型网格单元和塑性乘子数目对计算结果精度的影响。     

超浅覆大断面暗挖隧道下穿富水河道施工风险分析及控制研究

以广－珠铁路江门隧道下穿富水河道段为依托,针对其超浅覆、大断面、富水等特点,对超浅覆大断面暗挖隧道下穿富水河道的风险源进行系统分析,阐述风险产生的原因及造成的危害。在对可能采用的预加固手段及开挖方案进行初步比选后,采用三维数值模拟手段进一步量化比选,提出在地表隔水措施基础上,采用水平旋喷与大管棚复合超前支护并结合三台阶法开挖方案。研究表明,该方案能有效止水并加固地层,有效规避塌方风险,保证隧道施工安全的同时,加快施工进度及降低造价,可为类似工程提供参考。