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

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

封面图片说明

封面图为一组油气井地下渗流场的模拟图(带箭头的线为流线,其他为等值线,颜色深浅表示压力的大小,圆圈的核心为井点).油气井井组的渗流是油气田开发方案制定及方案调整的重要依据,无论是对常规油气藏还是对非常规油气藏都有着重要的意义.从油气藏井组的渗流场图可以看出:流体的流动方向、压力分布的不均匀程度以及不同井产量对整个井组渗流场的影响.从渗流场图还可以确定不流动区域的位置,为调整井的钻井提供可靠的依据.(图文供稿:刘曰武,中国科学院力学研究所)     

LARGE-STRAIN CONSOLIDATION OF FOUNDATION WITH SAND DRAINS BASED ON HANSBO'S FLOW 1)

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

矩形页岩气藏压裂直井解析解产能评价

为了准确评价矩形渗流区域内页岩气藏压裂直井的产能,加强气井生产过程中的科学性管理。通过运用Fick扩散定律和达西渗流规律,建立页岩气藏压裂直井基质扩散方程和裂缝渗流方程。结合时间叠加原理,通过Green函数给出矩形封闭边界和矩形定压边界下无限导流垂直裂缝井不稳定渗流的解析解。研究结果明确了页岩气藏矩形边界压裂直井的渗流特征及不同扩散特征下的压力响应,提出了矩形页岩气藏压裂直井产能评价方法,达到了快速评价及准确预测气井产能的目的。矩形页岩气藏的流动能分为7个阶段。通过将本文提出的矩形封闭边界模型与实际生产数据拟合,得到单井产量平均拟合误差为11.6%,压力平均拟合误差为8.3%,累计产气量平均拟合误差为7.12%。     

凝析油气渗流理论研究进展

本文综述了凝析油气渗流研究的现状。简要介绍了凝析油气藏特征。开采方案及其试井分析方法,着重讨论了凝析油气渗流的基本方程、解析解法和数值模拟。在此基础上,提出了该理论研究的发展前景。      

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

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

????????????????????????????о?

基于射孔完井水平裂缝中支撑剂均匀充填，考虑充填层强度、油层和流体 物性以及裂缝的渗流条件，研究了地应力及渗流引起的附加应力联合作用的支撑剂充填层渗 流模型和力学模型，并进行现场应用. 研究表明，所建立的模型对压裂油井支撑剂回采出砂 状况分析具有较高的符合率，油井总体符合率、不出砂井符合率和出砂井符合率分 别为84.66\%, 84.18\%和88.10\%, 为大庆油田压裂油井防砂措施实施方案制定提 供了决策支持，有利于提高措施有效率.     

低渗透煤层气藏中气-水两相不稳定渗流动态分析

针对低渗透煤层渗流问题,考虑启动压力梯度及其引起的动边界和动边界内吸附气解吸作用的渗流模型研究目前仅限于单相流,而更符合实际的气-水两相渗流动边界模型未见报道.本文综合考虑了煤层吸附气的解吸作用、气-水两相渗流、非达西渗流、地层应力敏感等影响因素,进行了低渗透煤层的气-水两相渗流模型研究.采用了试井技术中的"分相处理"方法,修正了两相渗流的综合压缩系数和流度,并基于含气饱和度呈线性递减分布的假设,建立了煤层气藏的气-水两相渗流耦合模型.该数学模型不仅可以描述由于低渗透煤层中渗流存在启动压力梯度而产生的可表征煤层有效动用范围随时间变化的移动边界,还可以描述煤层有效动用范围内吸附气的解吸现象以及吸附气解吸作用所引起的煤层含气饱和度的上升;为了提高模型精度,控制方程还保留了二次压力梯度项.采用了稳定的全隐式有限差分方法进行了模型的数值求解,并验证了数值计算方法的正确性,获得了模型关于瞬时井底压力与压力导数响应的双对数特征曲线,由此分析了各渗流参数的敏感性影响.本文研究结果可为低渗透煤层气藏开发的气-水两相流试井技术提供渗流力学的理论基础.     

各向异性双重介质垂直裂缝井两相流体渗流

建立了各向异性双重介质中垂直裂缝井两相渗流数学模型,用有限差分法求得了其解,进行了算例分析。得到了含水饱和度沿径向、最大及最小渗透率方向的分布和水驱油前缘等饱和度面随时间的变化,分别讨论了各向异性、吸渗作用对这种介质中两相流体渗流的影响。结果表明：渗透率较在的方向,前缘推进速度较快,油井见水较早;吸渗使前缘推进速度减慢,使油井见水晚。注水开发这种油藏时,见水前注入率不能太大以充分发挥吸渗作用,存在最佳注入率。这对于水力压裂垂直裂缝井注水开发碳酸盐岩裂缝性油藏有指导意义。     

??????????????????????о?

根据岩心实验提出了两参数的非线性渗流方程,并在此基础上结合Buckly-Leverett 方程建立了非线性非混相驱替方程. 通过对该方程求解分析可知,与线性达西渗流相比,非 线性渗流的驱油相率较差;同一含水饱和度下的含水率更高;而储层中各点的压力梯度更大, 导致所需注采压差更大. 说明有效的开发低渗透油藏需要比中高渗油藏更小的井距和更高的 注采压力体系.     

复杂边界条件下渗流场流线分布研究

流线分布研究已受到油藏工程师们的普遍关注。本文从无限大油藏稳态渗流场基本解出发,结合边界元方法求解出复杂边界条件下稳态渗流场的势分布,以此为基础提出了流线场的生成办法,并给出了应用实例。本文方法的优点在于求解过程中将问题的维数降低了一维,减少了计算量;计算精度较高且具有一定的普遍应用性;适应于求解任意形状的包括定压、定流量或混合边界在内的组合边界问题。实例表明,利用本文方法产生的流线分布因能够较为直观地反映出油藏流体在注采井间的运动轨迹,为优化井网和注入方案提供了重要依据。     

钻井-固井过程软泥岩井筒结构完整性分析

在建立套管-水泥环-软泥岩地层井筒结构力学模型的基础上,基于多物理过程数值方法,研究了钻固井过程中近井区地层与固井水泥环的塑性损伤形变与应力演化规律。结果表明:钻井成孔过程和渗流-应力耦合效应为软泥岩地层井筒结构失稳的重要机制,固井水泥硬化过程可进一步促进近井区地层的塑性损伤,但程度有限;固井过程中,水泥环内产生的环向拉应力与径向压应力是导致水泥环开裂失效的主要原因;采用低弹性模量、高泊松比固井水泥浆体系可改善固井水泥环受力形态,有利于维持井筒结构在后期生产作业过程中的完整性。     

污染区对部分打开井井底压力响应的影响

钻井造成的污染区对部分打开井的影响不同于裸眼井. 为了分析污染区对部分打开井井底压力响应的影响,建立了一种部分打开井的二维轴对称渗流模型,模型考虑了真实的污染区以及储层渗透率各向异性特征.利用有限元数值方法对模型进行求解,获得了部分打开井的井底压力响应及储层压力分布. 根据压力响应及压力分布特征,将部分打开井的压力响应过程划分为5 个流动阶段,其中早期局部径向流动和椭球流动是该类井最典型的特征. 对污染区的影响分析表明:传统方法中的表皮系数S 并不等于污染区引起的机械表皮系数S_d;无量纲井筒储存系数不能与机械表皮系数组合. 修正了传统方法中部分打开井的井底压力公式,验证了部分打开井的总表皮计算公式,为该类井的井底压力响应解释及产能预测提供理论指导.      

复合地基群井固结模型理论研究

针对复合地基固结问题,考虑井阻和土体径竖向渗流的影响,建立以"桩–土–桩"为单元的复合地基群井计算模型,采用解析解法,推导出瞬时载荷作用下复合地基群井固结解析解,讨论了复合地基固结特点和孔压等的变化规律。研究表明,复合地基群井固结模型为复合地基固结理论的研究提供了一种新思路,在井径比大于5时,基于中心桩、边桩流量相等假设,复合地基群井固结理论所求固结度和传统固结理论所求固结度相差很小,可忽略不计;桩体渗透系数、内外扰动区土体径向渗透系数、土体竖向渗透系数越大,复合地基固结速率越快;堆载越大,最终达到稳定时的沉降量和竖向应变越大;地基表层土体的孔压消散速率和竖向应变速率大于地基深部土体。     

A CONSOLIDATION MODEL OF GROUP WELLS ON COMPOSITE FOUNDATION1)

针对复合地基固结问题,考虑井阻和土体径竖向渗流的影响,建立 以"桩-土-桩"为单元的复合地基群井计算模型,采用解析解法,推导出瞬时载荷作用下复合地基群井固结解析解,讨论了复合地基固结特点和孔压等的变化规律。研究表明,复合地基群井固结模型为复合地基固结理论的研究提供了一种新思路,在井径比大于5 时,基于中心桩、边桩流量相等假设,复合地基群井固结理论所求固结度和传统固结理论所求固结度相差很小,可忽略不计;桩体渗透系数、内外扰动区土体径向渗透系数、土体竖向渗透系数越大,复合地基固结速率越快;堆载越大,最终达到稳定时的沉降量和竖向应变越大;地基表层土体的孔压消散速率和竖向应变速率大于地基深部土体。     

变形三重介质低渗透油藏垂直裂缝井的非线性流动分析

不仅考虑低渗透油藏具有启动压力梯度的渗流特征,还考虑应力敏感地层中介质的变形;发展了Cinco-Ley H.提出的有限导流垂直裂缝井双线性流理论,将流体在垂直裂缝与地层中形成的流动划分为两个区域—垂直裂缝中的线性流区域和变形三重介质低渗透油藏中的非线性流区域;由此建立了变形三重介质低渗透油藏有限导流垂直裂缝井的非线性流...     

低渗透多孔介质渗流动边界模型的解析与数值解

考虑启动压力梯度的低渗透多孔介质非达西渗流模型属于强非线性动边界问题, 分别利用相似变量变换方法和基于空间坐标变换的有限差分方法, 对内边界变压力情况下、考虑启动压力梯度的一维低渗透多孔介质非达西渗流动边界模型进行了精确解析与数值求解研究. 研究结果表明:该动边界模型存在唯一的精确解析解, 且所求得的精确解析解可严格验证数值解的正确性;且当启动压力梯度值趋于零时, 非达西渗流动边界模型的精确解析解将退化为达西渗流情况下的精确解析解. 由求解结果作出的非零无因次启动压力梯度下的地层压力分布曲线表现出紧支性特点, 其与达西渗流模型的有显著不同. 因此, 研究低渗透多孔介质中非稳态渗流问题时, 应该考虑动边界的影响. 研究内容完善了低渗透多孔介质的非达西渗流力学理论, 为低渗透油气藏开发的试井解释与油藏数值模拟技术提供了理论基础.      

复合含水层地区深基坑降水三维渗流场数值模拟——以上海环球金融中心基坑降水为例

长江三角洲分布有巨厚、松散的第四纪沉积层,地下含水系统为一复合含水层系统,深基坑降水一般采用非完整井,且由于深基坑周围连续挡水墙难以深入含水层底板,所以其地下水渗流场变化极其复杂,具有明显的三维流特征。本文以上海环球金融中心深基坑降水为例,采用数值模拟方法,模拟了在多层含水层复合存在、含水层最深底板埋深达149m、基坑周围挡水连续墙埋深达34m、基坑内地下水位降深达26m的情况下的地下水三维渗流场特征,并以此为据确定出了该基坑降水的最优化方案。     

抽水引起的含水层水平应变──地裂缝活动新机理

基于Darcy -Gersevanov广义渗流关系式 ,本文首先导出了单井定流量抽水引起的泰斯承压含水层水平运动和应变解析表达式 ,并据此提出了抽水诱发地裂缝活动的含水层水平应变新机理 ,其主要观点为 :受井栅的约束阻挡作用 ,抽水活动将在含水层内部形成近井径向挤压区和远井径向拉张区两个应变区域 ,随着抽水时间的增长 ,近井径向挤压区将逐渐向外扩展 ;当地裂缝及其下伏断层或裂隙位于抽水活动引起的径向拉张应变区域时 ,地裂缝将出现加速活动 ;反之 ,当地裂缝处于抽水活动引起的径向挤压应变区域时 ,地裂缝活动反而受到抑制。上述地裂缝活动机理通过大同机车工厂现场抽水实验得到了初步证实。     

渗流力学研究的现状和发展趋势

渗流力学在能源、环境、水利、岩土、交通、生物等工程领域有广泛的应用，渗流力学经过约一个半世纪的发展，已经积累了相当多的成果。随着现代科学的发展以及生产实践需求的不断提高，渗流力学仍然有进一步完善和发展的广阔空间。本文简述了水利、环境、能源工程和生物学中的渗流力学问题，扼要概括了渗流力学理论研究的现状，并指出渗流力学在多孔介质描述、裂缝型介质渗流、多相多组分渗流、物理化学渗流、非线性渗流、非饱和渗流、微观渗流、渗流模拟等理论及相关方法和测试技术等方面的发展趋势。