考虑桩体内径向渗流的复合地基固结解

考虑了桩体内径竖向组合渗流和涂抹区土体水平渗透系数的两种变化模式,得到了此类固结问题的控制方程.采用新的初始条件,得到了控制方程的解答,给出了桩周土体和桩体内的平均孔压.在此基础上,把桩体的固结压缩引人对复合地基总平均固结度的定义中,分别给出了复合地基按应力和按应变定义的总平均固结度.结果表明:按应力和按应变定义的复合地基的总平均固结度相等;考虑桩体内径向渗流比不考虑桩体内径向渗流的复合地基固结度要小,而且随着桩间距的增大,两者之间的差值逐渐减小;考虑涂抹区水平渗透系数逐渐变化比假定涂抹区水平渗透系数保持不变的复合地基的总平均固结度大.     

考虑Hansbo渗流的砂井地基大变形固结问题

基于Gibson大变形固结理论,通过引入Hansbo渗流模型和软黏土非线性变形特性,同时考虑土层的沉积效应,分别建立了以孔隙比e和超孔压u表示的砂井地基大变形固结方程。通过与已有研究成果的对比,验证了本文方法的可靠性和砂井地基大变形固结方程两种描述方式的等效性。通过FlexPDE得到方程的数值解,在此基础上,研究了砂井地基大、小应变固结理论与Barron固结解的差异性,探讨了Hansbo渗流参数m和I1对砂井地基固结的影响,最后对比分析了砂井地基轴对称固结、径向固结和竖向一维固结的关系。研究结果表明:考虑Hansbo渗流的砂井地基非线性大变形固结模型的固结速率最慢,且在固结后期,Hansbo渗流时大、小应变固结与Barron固结的平均固结度基本趋近;随着Hansbo渗流参数m和I1的逐渐增大,砂井地基的固结速率逐渐降低;随着砂井影响半径的增大,砂井地基轴对称固结与径向固结的差异性也越来越大,且在固结早期会出现一维竖向固结速率高于轴对称固结和径向固结的现象。     

考虑介质微观非均质特征的乳状液渗流模型

油水乳化渗流是三元复合驱和热力采油过程中常见的现象,地层介质的微观孔隙结构特征对乳状液流动有着重要影响.现有描述乳状液渗流的理论模型都属于确定性方法,只能反映出介质孔隙结构的体积平均效果.当介质内部微观非均质性相比其尺寸不能被忽略时,采用确定性方法描述会与实验结果存在偏差.基于连续时间随机游走理论建立了描述乳状液渗流的随机理论模型.该模型引入反映液滴微观运动特征的跃迁时间和跃迁位移两个概率分布函数来反映多孔介质微观非均质特征.研究结果表明该模型能很好地刻画实验曲线中出现的与介质尺度相关的拖尾现象,可作为更一般的过滤模型.     

非达西渗流拟启动压力梯度推算

岩体启动压力梯度的大小是石油开采,工程防渗处理等方面需要考虑的影响因素。非线性渗流现象已被广泛所认同,低渗透性岩渗流存在启动压力梯度,它基本代表了流体产生渗流时的压力梯度大小。试验对试样进水端水压进行长期稳定控制,出水端采用精确测定流出水体积变化量的方法进行流量测定; 当渗出端水体积的变化量与时间呈直线关系时,即认为该时段渗透特征符合达西定律。通过对具有不同渗透性的软弱岩系统的试验测试,根据实测渗透系数与压力梯度的正相关关系,提出了拟启动压力梯度推算的"三次平均法",即:首先统一调整实测渗透系数值,然后确定绝对渗透系数值,再根据该值确定V0值, 3次取平均推算拟启动压力梯度的方法。分析表明:拟启动压力梯度随围压升高而增大,随岩石的绝对渗透系数降低而增大; 在渗透压差增大过程测得的拟启动压力梯度大于降低过程的拟启动压力梯度,这些都与理论分析结果相吻合。实例说明拟启动压力梯度推算"三次平均法"较为合理,易于推广使用。这为通过室内试验确定工程岩体启动渗流压力梯度的大小,提出了一种新的思路。     

考虑Hansbo渗流的砂井地基大变形固结问题 1)

基于Gibson大变形固结理论,通过引入Hansbo渗流模型和软黏土非线性变形特性,同时考虑土层的沉积效应,分别建立了以孔隙比$e$和超孔压$u$表示的砂井地基大变形固结方程。通过与已有研究成果的对比,验证了本文方法的可靠性和砂井地基大变形固结方程两种描述方式的等效性。通过FlexPDE得到方程的数值解,在此基础上,研究了砂井地基大、小应变固结理论与Barron固结解的差异性,探讨了Hansbo渗流参数$m$和$I_1$对砂井地基固结的影响,最后对比分析了砂井地基轴对称固结、径向固结和竖向一维固结的关系。研究结果表明：考虑Hansbo渗流的砂井地基非线性大变形固结模型的固结速率最慢,且在固结后期,Hansbo渗流时大、小应变固结与Barron固结的平均固结度基本趋近;随着Hansbo渗流参数$m$和$I_1$的逐渐增大,砂井地基的固结速率逐渐降低;随着砂井影响半径的增大,砂井地基轴对称固结与径向固结的差异性也越来越大,且在固结早期会出现一维竖向固结速率高于轴对称固结和径向固结的现象。     

考虑次固结的软土地基一维非线性固结理论

将软黏土变形分为有效应力变化引起的变形和次固结引起的变形,推导了软土新型应力应变关系。然后结合Davis固结理论,建立了考虑次固结的一维非线性固结控制方程,并对其进行解析求解。通过与数值方法对比,验证了解析解的可靠性。在此基础上,分析了次固结对软土地基固结沉降的影响。结果表明:考虑次固结的孔压消散速度与固结速度较不考虑次固结的慢;忽略次固结将低估软黏土地基的工后沉降。     

循环荷载作用下考虑Hansbo渗流的饱和黏土固结分析

针对传统一维Terzaghi固结理论的局限性,本文引入了考虑时间效应的UH (unified hardening)本构模型,将Hansbo渗流模型耦合到UH本构模型固结方程中,建立了几种典型循环荷载作用下基于UH模型和考虑Hansbo渗流模型的固结方程。采用FlexPDE软件进行了数值分析,将计算结果与已有文献结果对比,验证了该算法的可靠性。在此基础上,研究了饱和黏土在几种典型循环荷载作用下的渗流固结特性,分析了各模型参数对饱和黏土固结特性的影响。结果表明：在循环荷载作用下,地基的平均固结度、沉降量等始终处于循环状态;Hansbo渗流模型参数对平均固结度的影响相对于Darcy渗流参数影响程度较大,其中UH模型中次固结系数对固结过程影响程度较大。另外,Hansbo渗流参数及土的回弹指数和渗透指数等对固结过程的影响主要体现在中期,固结过程后期将处于稳定循环状态。     

饱和土地基的非轴对称固结

固结问题是土力学的基础课题之一，现有的研究主要针对饱和土的平面固结和轴对称固结问题。本文应用Fourier展开、Laplace变换和Hankel变换求解了各向同性饱和土的Blot固结方程；讨论了半空间各向同性饱和土的非轴对称固结。借助数值计算，本文的方法可以用于分析各向同性饱和土的非轴对称固结问题。     

采空区流场非达西渗流一种新的迭代算法

对Bachmat非线性渗流方程求解问题，结合采空区问题的力学特点，提出基于变渗透率达西 (Darcy)渗流求解非线性渗流的迭代算法. 建立了冒落采空区非线性漏风流态的有限元数 值模型，通过对复杂边界采空区的漏风渗流的计算，得到与实际流态更接近的流动规律(风 压等值线和流函数线)和采空区漏风强度分布(速度场). 结果表明，迭代算法是振荡性收 敛的，收敛速度快；与Darcy渗流相比，采空区非线性渗流速度场趋向于平缓，在工作面附 近30\,m范围内最大风速降低到原来的0.61倍，其余大部分区域速度略有增加. 迭代方法满 足工程要求.     

非牛顿流体在非均质油藏渗流压力场实验

在非均质油藏模型上进行非牛顿流体流动物理模拟实验，对比研究水驱、聚合物驱和交联聚合物对提高石油采收率的影响．通过布置高精度的压差传感器测量不同驱替过程模型中的渗流压力场的动态变化，成胶后的交联聚合物封堵了高渗条区，改变了油藏内流体流动方向，驱替出低渗区内油，提高了采收率．     

复杂边界非均质渗流场流线分布研究

建立了考虑源(汇)影响的含有不渗透区域复杂边界条件下非均质油藏稳定渗流的数学模型。利用扰动边界元方法求解数学模型,获得了地层中任意一点的压力公式.在此基础上,提出了流线场的生成方法。绘制了考虑非均质性、复杂边界和不渗透区域影响的流线分布图,并分析了流线分布的特征。通过分析表明,渗流场的非均质性和不渗透区域的存在都对流线分布存在较大的影响。利用本文方法产生的流线分布图能够较为直观地反映出油藏流体在注采井间的运动轨迹,为优化井网和注入方案提供了重要依据。     

动边界双重介质油藏低速非达西渗流试井模型

裂缝性油藏中基质岩块的渗透率一般很低,大量岩心测试实验证实在基质岩块内的液体渗流和在一定含水饱和度下的气体渗流将偏离达西渗流,往往出现低速非达西渗流,表现出启动压力梯度以及流体流动边界不断向外扩展等特殊现象。本文充分考虑启动压力梯度与动边界的影响,建立了微可压缩双重介质油藏低速非达西渗流的试井数学模型,对时间和空间变量...     

考虑固结的基坑开挖回弹问题有限层求解

在综述前人有关基坑变形研究成果的基础上,指出开挖回弹量直接影响到基坑的稳定性和建筑物的后期沉降,但相关研究很少考虑土体在开挖过程中的固结影响.现引入固结有限层法,通过增加开挖部分的耦合矩阵和开挖部分的渗透矩阵,建立三维横观各向同性土体考虑固结的开挖回弹问题的有限层求解格式,并编制了相应的计算程序,以分析开挖过程中坑底土体的动态固结过程和负孔压的消散规律,研究固结对坑底土体回弹的影响程度.结果表明坑底土体的孔压随时间逐步消散,从负孔压慢慢消散到零,且符合Mandel-Cryer效应,坑底的回弹量随时间不断增加,且中心点回弹量较大,边缘点回弹量相对较小,总体成中心对称的倒锅底形.     

用DQM求解成层地基一维非线性固结问题

由于多层地基的一维非线性固结问题求解的复杂性,其解析解很难求得。本文基于Davis和Raymond一维非线性固结理论,利用DQM(Differential Quadrature Method)导了初始有效应力沿深度变化、任意边界条件、任意荷载作用下成层地基一维非线性固结的统一表达式,求得了孔压、有效应力和平均固结度的解答。通过解的收敛性分析讨论了DQM解的有效性。由于DQM解对于固结间题各种复杂条件具有统一的矩阵表达式,更便于编程计算和工程应用。最后,用本文解答对三层地基一维非线性固结问题进行了讨论。     

径向积分边界元法确定非均质土石坝渗流自由面

浸润面位置的确定是无压渗流分析中非常重要的问题,本文将径向积分边界单元法应用到渗流问题中,通过直接将区域积分转化为边界积分的技术,克服了传统边界元解决该类问题的缺陷,编写了迭代计算程序,并用二、三维算例证明了算法的有效性。     

基于分形理论的煤层气非达西流动分析

在考虑了煤层的分形特征和启动压力梯度影响的基础上,建立了无限大煤层中气体低速非达西流动的数学模型,并求得了量纲为一的井底压力的Laplace空间解析解,并根据数值求解结果绘制了典型的井底压力动态曲线。研究结果表明:在定产量生产的情况下,分形维数和量纲为一的启动压力梯度对早期井底压力动态无显著影响;在生产的中后期,由于受二者的影响,压力导数曲线上的径向流水平直线段消失;量纲为一的井储系数的影响主要表现在续流阶段,而吸附系数则主要影响煤基质向裂隙扩散的过渡阶段。     

基于Hansbo渗流的圆柱土样Biot固结分析

诸多黏性土渗流试验表明,在低水力梯度下渗流会出现明显偏离Darcy定律的现象.为了分析渗流的非Darcy特性对固结过程的影响,引入Hansbo渗流方程描述圆柱土样内的渗流,重新推导轴对称条件下的Biot固结方程,并给出方程的Crank--Nicolson有限差分格式.通过与Darcy渗流条件下轴对称Biot固结方程解析解的对比,验证计算方法的有效性.然后分析Hansbo模型参数对圆柱土样固结过程的影响.计算结果表明:与Darcy渗流相比,Hansbo渗流会延缓圆柱土样的固结过程.随着Hansbo渗流参数m或I1的增大,在固结前期,Mandel--Cryer效应会更加显著,即孔隙水压力峰值将提高,且达到该峰值的时间会延迟;在固结中后期,孔压消散滞后的现象也更加明显.不过,Hansbo渗流对位移的影响很小.     

渗透各向异性可压缩地基固结的平面应变分析

对多层地基的平面应变固结问题进行了研究,并同时考虑了土体的渗透各向异性和孔隙流体的可压缩性.从平而应变Biot固结的控制方程出发,对时间t,坐标z和x进行Laplace和Fourier变换,建立了地基表面(z=O)和任意深度z处的基本量在Laplauce-Fourier变换域内的传递矩阵关系.利用传递矩阵法,结合土层连续条件和边界条件,并应用Laplace-Fourier逆变换技术,推导出渗透各向异性可压缩多层地基平面应变固结的理论解.基于该解,编制了计算程序,并进行了数值计算.讨论了土体的渗透各向异性、孔隙流体的可压缩性以及地基的分层特性对地基同结的影响,分析结果表明:土体的渗透各向异性、孔隙流体的可压缩性,以及地基的分层特性对地基的固结行为有着重要的影响.     

解析法求解成层渗透各向异性地基Biot固结轴对称问题

采用解析法研究成层渗透各向异性地基,该法从渗透各向异性Ｂｉｏｔ固结轴对称问题的基本方程（静力平衡方程,物理方程及渗透连续方程）出发,利用Ｌａｐｌａｃｅ～Ｈａｎｋｅｌ变换及有关矩阵理论等,得到Ｂｉｏｔ固结基本量不同深度之间的传递矩阵。利用传递矩阵,边界条件以及Ｌａｐｌａｃｅ～Ｈａｎｋａｌ逆变换技术可求解多层渗透各向异性地基体系,采用更为有效的Ｆ．Ｄｕｒｂｉｎ的方法实现Ｌａｐｌａｃｅ逆变换。编制了计算     

考虑在地基固结和流变的上部结构和地基基础的共同作用问题

本文提出一种能同时考虑地基的瞬时变形,团结变形和流变的地基基础与上部结构共同作用的一种新的计算方法,在这种方法中,土的计算考虑了Biot固结理论和流变,上部结构利用子结构法进行凝聚可得边界刚度矩阵和边界力,利用位移协讯条件则可以得到共同作用的控制方程,对上述方程进行Laplace变换,可以得到Laplace变换域内的控制方程,在Laplace变换域内对上述方程进行数值求解并进行相应的Laplace逆变换则可以得到时间境内任意时刻的解。