镇泾油田致密砂岩钻井泥浆侵入正、反演计算

钻井泥浆侵入到地层中会影响储层的饱和度和地层电阻率,通过油、水两相渗流模型和混合流体模型正演出泥浆侵入后的地层饱和度和电阻率分布,反演过程对阵列感应采用最优化方法来完成,并对混合流体正演模型和反演模型进行了修正,修正后的模型计算结果更准确。镇泾油田致密砂岩泥浆侵入深度和电阻率的正、反演计算结果表明,泥浆的侵入深度与钻井泥浆滤失量、地层孔隙度和持水率有关,泥浆侵入深度的复杂性与地层孔隙结构的复杂性是一致的。     

基于Fredlund & Xing模型的渗流分析在川东红层梯田滑坡中的应用

四川省巴中市通江县新场镇七家沟村二社水头上滑坡为川东典型红层梯田滑坡,区内类似斜坡覆盖范围广泛,但针对该地区此类滑坡的渗流分析研究较少。在对滑坡区进行大量野外调查、勘查、资料收集的基础上,综合运用卫星遥感、无人机航拍、机载LiDAR、室内试验等技术手段,基于Fredlund&Xing土水特征曲线数学模型,采用Geo-Studio中SEEP/W模块进行渗流分析,将不同时间阶段分析结果与SLOPE/W模块耦合,继而得到稳定系数与降雨、时间之间的动态变化关系,揭示了滑坡的变形过程和形成机理。研究显示：（1）连续降雨促使坡体地下水位升高,稳定性降低,最终导致水头上滑坡整体失稳破坏;（2）滑坡变形过程、地下水出露特征与Geo-Studio计算结果基本吻合,说明基于Fredlund&Xing数值模拟在类似地区能提供较为准确的模拟结果;（3）开垦梯田（水田）会降低坡体稳定性。研究结果可为四川山区类似滑坡灾害进行隐患排查和主动防范提供理论支撑,为防灾减灾提供参考依据。     

不同水力荷载路径下煤体微观渗流特征及宏观破坏研究

为探究煤体在不同水力荷载路径下的微观渗流特征,通过CT三维重建技术,建立了基于煤体微观结构的双介质渗流模型,并设计了恒定速度加载（LP1）、常规脉动加载（LP2）和强加载缓卸载加载（LP3）3种水力荷载路径,进行了在不同循环荷载路径下的煤体注水渗流数值模拟试验。另外利用自行搭建的煤体注水装置,探究了3种荷载路径下的宏观损伤情况。结果表明：联通裂隙占孔裂隙结构的73.49%,是影响煤体渗流的主要因素,而孔径为9～23μm的孔隙是孤立孔隙中的主要部分,数量和体积占比均超过了50%;煤体内部孔裂隙结构和荷载路径对渗流速度和压力的分布特征具有重要影响。LP1沿注入方向平均渗流速度的波动较大,LP3则抑制效果明显,且LP3路径较LP1和LP2随时间有较大跨度。LP3的荷载路径相对于LP1具有脉动荷载的疲劳损伤效果,并且造成的损伤程度强于LP2。该研究可为煤体微观结构研究和煤层注水技术参数优化提供方向。     

考虑潜水面的基坑二维稳态渗流场解析解

对考虑潜水面悬挂式止水帷幕支护下的基坑二维稳态渗流场进行了解析研究,并给出一种求解潜水面的解析方法。根据对称性取基坑半截面,并根据边界的连续条件将其分为3个规则区域,使用分离变量法分别将3个区域的水头分布表示为级数解的形式,结合区域之间的连续性条件及级数正交条件得到各区域渗流场的显式解,根据潜水面所满足的总水头等于位置水头的条件确定潜水面。将解析解方法计算结果与室内试验、有限元分析结果进行比对分析,结果验证了解析解的正确性,解析方法相较于有限元数值方法具有更高效的计算效率。对潜水面位置进行参数分析,发现止水帷幕插入深度、基坑宽度及深度对潜水面位置有不可忽视的影响。随着过水断面厚度的增加,潜水面位置逐渐降低,不考虑极端状态的情况下可认为帷幕底部至不透水层顶面的距离与潜水面在止水帷幕上的位置呈线性关系;随着基坑尺寸的增加,潜水面位置呈下降趋势,基坑内侧半宽度对潜水面的影响明显小于帷幕底部至不透水层顶面的距离,且随着基坑内侧半宽度的增加影响越来越小;潜水面位置随着基坑深度的增大而降低,基坑深度对潜水面位置的影响相对较大。     

循环剪切作用下三维粗糙裂隙非线性渗流特性数值模拟研究

研究了循环剪切过程中剪切位移和循环次数对三维粗糙裂隙非线性水力特性的影响。首先,基于天然粗糙裂隙面的高程数据生成三维裂隙面,并通过前人的研究结果确定循环剪切过程中的法向位移变化值,建立了16个开度空间模型,从而得到了不同剪切位移和循环次数下接触面积和开度分布情况,发现当循环次数一定时,随着剪切位移的增加,原本吻合的上下表面发生错动,平均开度和各向异性增大,接触面积减小;剪切位移一定时,随着循环次数的增加,上下表面起伏在剪应力作用下被磨平,接触面积增大,平均开度减小。将模型导入数值模拟软件进行计算得到流量数据,结果表明,水力梯度随流量变化规律符合Forchheimer定律。循环1次时,随着剪切位移从4 mm增加到10 mm,线性Fochheimer系数减小了88.3%,非线性Fochheimer系数减小了95.2%;循环2次时,随着剪切位移从4 mm增加到10 mm,线性Fochheimer系数减小了95.4%,非线性Fochheimer系数减小了99.7%,两系数减小的幅度随循环次数增大而增大。进一步分析发现,剪切位移的增加对裂隙渗透性有促进作用,循环次数增加对渗透性起抑制作用。此外,...     

考虑渗流的饱和粉土地层盾构开挖面稳定分析

土压平衡盾构在穿越高水位钱塘江粉土地层时,地下水与密封舱之间的高水压差会产生过大的指向开挖面的渗透力。钱塘江粉土黏聚力低,开挖面自立性差,且土体模量小,过大的渗透力会使土体产生更大的变形,导致开挖面失稳。数值模拟中合理选取土的本构关系与模型参数,对计算结果的可靠性十分重要。为了研究渗流条件下开挖面稳定性问题,基于室内试验,对钱塘江粉土的摩尔-库仑本构模型和摩尔-库仑应变软化本构模型参数进行研究。研究发现,钱塘江粉土的剪应变是轴向应变的1.5倍;根据三轴卸荷试验得到应变软化内摩擦角和剪应变的关系,当剪应变为4.7%时内摩擦角达到峰值为25°,剪应变为7.4%时内摩擦角达到残余值为21°,剪应变在4.7%～7.4%之间时,内摩擦角线性减小。采用有限差分软件建立39组三维分析模型,通过与离心机试验的极限支护压力和地表沉降对比发现,相同水位情况下,摩尔-库仑应变软化本构模型计算的极限支护压力与离心机试验的最小误差为3.1%,最大地表沉降的误差为5.67%。证明了在钱塘江粉土地层中使用摩尔-库仑应变软化本构模型的可行性及模型参数的可靠性。     

岩土基坑开挖过程中的周围土体渗流-应力耦合模型

传统渗流-应力耦合模型仅考虑土体的各向同性特征,忽略了工程中应力与渗流的非线性影响,导致计算结果与实际情况偏差较大,由此,设计一种岩土基坑开挖过程中的周围土体渗流-应力耦合模型。搭建测试装置,利用瞬态试验法测定岩土基坑周围土体渗透率,分析岩土基坑周围土体基质块微观结构,得到土体变形方程并求解;从宏观角度分析土体应力,得到应力边界的面力结构,构建土体渗流-应力耦合模型。根据实际情况设定工程边界条件,选择某建筑工程中岩土基坑作为分析实例,测试结果表明,该模型计算结果与实际情况的误差低于0.2 m,可以为岩土基坑开挖提供参考。     

考虑库水升降和滑带弱化作用的岸坡启滑机制分析

以三峡库区凉水井滑坡为研究对象,采用理论分析与数值模拟的研究手段,构建了滑带强度弱化模型,提出了渗流驱动的滑坡启滑判据,应用Geo-studio有限元程序分析了库水不同升降速率对滑坡稳定性的影响,揭示了库水升降作用下岸坡渗流场演变规律和渗流驱动下的启滑机制。研究表明：（1）渗透压力与渗透时间的变化是滑带土强度弱化的关键因素,弱化到临界强度时,在渗流驱动下发生压剪破坏而启滑,由局部向整体以渐进模式破坏失稳;（2）库水升降过程中,坡体内孔隙水压力滞后性较明显,水位升降速率会影响坡体地下水响应时程,升降速率越大,孔隙水压力变化越大,渗流驱动力越大,滑坡稳定性变化越快,越趋近于渐进破坏;（3）库水位从175 m降到145 m,凉水井滑坡滑面法向应力最大降低了38.19%,剪应力最大降低了22.20%,有效法向应力最大落差为168.64 kPa,抗剪强度最大落差为63.45 kPa。以上分析结论与规律可为涉水滑坡启滑机制分析、库岸山体滑坡失稳研究及应急防治工程等提供科学依据和理论方法。     

基于四参数随机生长法重构土体的格子玻尔兹曼细观渗流研究

多孔介质模型的重构问题是土体细观渗流机理研究的基础和关键。由四参数随机生长法(QSGS)构建土体模型,采用格子玻尔兹曼方法(LBM),通过MATLAB自编程序研究重构土在不同条件下的细观渗流机理。结果表明:随模型尺寸增大,孔隙连通程度显著提高,300×300格点大小的模型连通孔隙率增长幅度(34.38%)最大,继续扩大模型尺寸发现增加不明显;流体粒子在孔隙连通性好、孔径大的区域,会形成主渗流通道,且存在指进效应,孔道中间流速最大,可达0.0324,越靠近孔壁流速越小;大孔隙率土的流速比小孔隙率土大,而低孔隙率土中的流速相比大孔隙土更稳定;LBM模拟渗透率与经典K-C公式计算结果对比发现,孔隙率越高计算渗透率越准确(n=0.78,误差为10.22%);土颗粒越小,渗流孔道越细窄、分布越密集,对应的速度场分布更为均匀,同时流速也更小。该研究成果能较好地揭示重构土的细观渗流机理,也可为现有细观土体孔隙流研究提供一定借鉴。     

岩石弹塑性应力-渗流-损伤耦合模型研究(Ⅱ):参数反演及数值模拟

海底隧道或近海富水区隧道围岩长期处于地下水环境中,围岩稳定性受渗流场的影响较为明显,由于渗流场作用使得应力场、损伤场发生变化,而围岩应力场、损伤场的变化又对渗流场产生反作用,三场耦合效应十分显著。针对耦合模型中参数多、确定难度大的问题,进行耦合模型中损伤参数的反演。采用基于岩石弹塑性应力-渗流-损伤耦合模型所编制的计算程序和智能位移反分析程序,对大连地铁海事大学试验线路过河段隧道施工过程中的围岩稳定性进行数值计算。根据现场监测位移采用耦合模型进行损伤参数反演,其中耦合计算中采用应力场与渗流场分别迭代求解的间接耦合方法进行有限元计算,利用反演的参数对隧道围岩应力场、渗流场、损伤场分布规律及衬砌结构的受力特征进行了分析。研究结果表明:利用位移反分析法得到的围岩力学参数进行类似地质条件的隧道围岩数值分析是可行的,进而可以预测围岩的变形破坏模式,判断围岩的稳定性。与此同时,通过数值计算可知,地下水的渗流作用对近海隧道的围岩变形有一定的影响,增加了围岩的应力、位移,从围岩-支护结构共同作用原理考虑,进行隧道支护结构设计是应该考虑三场耦合效应的,计算结果可以指导隧道防排水施工质量的改进与提高,为近海富水区隧道开挖设计提供一定的理论参考。