低渗透裂缝性油气藏非稳态窜流因子研究

深入分析了常定窜流因子和现有的非稳态因子对于低渗透基质-裂缝窜流的不适应问题,根据对23 个低渗岩心渗流曲线的拟合分析结果,引入了变渗透率系数,并在此基础上推导了非线性扩散方程,也就是低渗透基质-裂缝系统间窜流的控制微分方程;经过无量纲化处理,引入了动边界条件,分别利用积分方法和矩方法导出了非线性扩散方程在前后两个阶段的近似解析解,并在此基础上构建了非稳态窜流因子的表达形式;新的窜流因子可以有效应用于低渗透基质-裂缝系统间的窜流计算,有限元数值模拟的结果验证了新因子在计算非线性非稳态窜流方面的准确性与可靠性.      

低渗透多孔介质中的非线性渗流理论简

文中论述了低渗透性多孔介质中非线性渗流理论的几个问题,阐明了渗流流体的性质,指出了多孔介质对流体通过的选择性,提出了新的非线性渗流方程,用实验资料对其进行了验证,分析了该方程演变功能,表明它可以描述各种渗流规律.该方程的各项参数都可从实验中直接得到,应用方便,并且参数的物理意义明确.     

结构可靠性仿真方法研究简

可靠性仿真是对复杂结构系统进行可靠性分析计算最为有效的方法,引起了研究者越来越广泛的关注. 综合目前结构可靠性仿真方法研究现状,介绍了蒙特卡罗（Monte-Carlo,MC）法、极限状态方程重构法、随机有限元法（stochastic finite element method,SFEM）进行结构可靠性仿真计算的特点及发展现状,归纳了固体火箭发动机（solid rocket motor,SRM）结构可靠性的特点及当前进行仿真的方法;并在此基础上分析了目前结构可靠性仿真方法及对SRM 结构可靠性进行仿真分析计算时存在的不足及需要进一步研究的内容.     

幂律流体圆柱绕流的格子波尔兹曼模拟简

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明：插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.     

串列风力机尾流干扰的研究简

广义致动盘方法是通过引入体积力代替叶片的致动盘技术与三维Navier-Stokes 方程相结合来获得风力机周围流场信息的一种方法. 该方法避免了花费大量网格与计算资源去求解风力机叶片的附面层,从而可以把更多的网格与计算资源用于风力机尾流流场的模拟,非常适合用于风力机尾流流场的研究. 以NH1500风力机为计算模型,将常规CFD (computational fluid dynamics) 方法与广义致动盘方法计算所得的叶片载荷分布进行比较,以验证广义致动盘方法的可行性. 然后使用广义致动盘方法对风场中串列风力机进行数值模拟,研究串列风力机之间间距变化时,上游风力机产生的尾流对下游风力机的干扰影响.     

新型涡轮流量计的理论模型简

基于动量力矩定理和流体边界层理论,建立了一新型涡轮流量计的理论模型. 利用该理论模型可以分析叶轮几何参数对流量计计量性能的影响. 以一DN50 涡轮流量计为例,利用本文所建立的理论模型对该流量计的仪表系数进行了计算,并在一体积管标准装置对该流量计的仪表系数进行标定. 结果表明,计算值和实验值较为吻合,误差在±;3.5% 以内,验证了模型的有效性,从而为仪表的结构优化设计提供了理论依据.     

浮泥对非线性波的运动响应简

将浮泥作为高黏度牛顿流体,与水体构成互不相混的两相流体,基于垂向二维RANS 方程(Reynolds-averaged Navier-Stokes)、SST 湍流模式(shear-stress transport),利用流体体积法(volume of fluid,VOF),结合数值造波、消波技术模拟波浪,建立了波泥耦合运动的数学模型,研究了线性波、Stokes波和椭余波作用下,浮泥界面的运动响应. 发现在各种波浪条件下,界面位移与表面位移间有相近的波形,两者间存在相位差,其值随周期和波高的变化不大;界面波高基本正比于表面波高,它随周期的增长呈单调增大的趋势.     

基于有限变形理论的数值流形方法研究简

原有数值流形方法通过积累每一时步的小变形而得到结构最终的大变形,然而,当结构发生大变形、大转动时往往产生较大计算误差. 针对该问题,从动量守恒方程以及应力边界条件的积分弱形式出发,引入流形方法的插值函数,建立了基于有限变形理论的数值流形方法. 通过对比改进前后流形方法的计算迭代格式,指出了原有流形方法计算大变形问题时的误差来源. 最后,通过大变形悬臂梁和旋转块体算例对有限变形流形方法进行了验证. 数值结果表明,改进后的流形方法能够很好地处理大变形大转动问题,消除了转动所带来的计算误差,其计算结果与解析解及ABAQUS 软件求得的数值解相吻合.     

非均布随机参数结构非线性响应的概率密度演化简

首先考察了概率密度演化理论中的点演化和群演化与概率空间剖分的关系. 继而,讨论了点集筛选的基本准则. 在此基础上推广了点集偏差的概念,对非均匀、非正态的一般多维分布,提出了广义F 偏差（GF 偏差）的概念,避免了偏差计算的NP 难解问题. 探索了GF 偏差与EF 偏差的关系. 以GF 偏差最小化为准则,建议了概率空间最优剖分与点集重整的新策略. 结果表明,上述方法能够处理包含多达数10 个随机变量的结构动力响应概率密度演化分析问题. 最后,指出了需要进一步研究的问题.     

低稀释度条件下乙烯点火特性的激波管研究简

在激波管中利用反射激波后高温环境加热燃料,以燃料点火过程中氢氧自由基特征发射光谱强度的急剧变化作为点火发生的标志,在温度范围800 1 650K,压力0.2MPa, 0.7MPa, 1.2MPa,化学当量比为0.5, 1, 2,O₂ 浓度为空气含量20% 的条件下,进行了C₂H₄/O₂/Ar 混合气在低稀释度条件下点火特性的实验研究. 获得了乙烯点火延时随温度、压力、化学当量比、燃料以及氧化剂浓度等参数变化的拟合关系式. 对乙烯点火转爆轰现象进行了初步观察,考察了初始温度对乙烯点火特性以及点火转爆轰的影响.     

裂缝性低渗透油藏流-固耦合理论与数值模拟

根据裂缝性低渗油藏的储层特征,建立适合裂缝性砂岩油藏渗流的等效连续介质模型。将渗流力学与弹塑性力学相结合,建立裂缝性低渗透油藏的流-固耦合渗流数学模型,并给出其数值解.通过数值模拟对一实际井网开发过程中孔隙度、渗透率的变化以及开发指标进行计算,并和刚性模型以及双重介质模型的计算结果进行了分析比较.     

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

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

低渗油层压裂水平井两相流研究

依据压裂水平井不同流动区域的流动规律, 将压裂水平井的渗流分为裂缝中的高速非达西流动区、裂缝控制影响的椭圆渗流区、远离裂缝的基质非达西渗流区, 考虑启动压力梯度的影响, 对压裂水平井两相渗流进行了分析, 得到了低渗透油层压裂水平井的产量公式. 研究结果表明, 裂缝的导流能力越大, 压裂水平井的产量越高. 但随着开采时间的增加, 其产量递减幅度越大; 压裂裂缝长度越小, 压裂水平井的初始产量越高. 但随着生产时间的推移, 压裂裂缝的长度越大, 产量的递减幅度越小; 中间裂缝长, 两翼裂缝短的情况下, 压裂水平井的产量最高.     

低渗岩心驱替过程中的压力分布特征研究

考虑低渗岩心的应力敏感性,通过实验和数值计算,研究了一维水驱油过程中,压力的分布及变化规律,并对其影响因素进行了分析. 结果表明：压力分布由于受到应力敏感性的影响而显示出明显的非线性特征,在前缘面处存在一个明显的分界点;应力敏感性程度、黏度比对压力分布的影响较大.     

页岩储层体积压裂的地应力变化研究

页岩储层属于致密超低渗透储层,需改造形成复杂缝网才有经济产能.体积压裂是页岩储层增产改造的主要措施,而地应力场特别是水平主应力差值是体积压裂的关键控制因素. 理论研究表明:(1)当初始两向水平主应力差较小时,容易形成缝网,反之不易产生缝网;(2)人工裂缝的形成能够改变地层初始应力场. 因此应在前人研究的基础上优化设计压裂方式,以克服和翻转初始水平主应力差值,产生体积缝网.基于此,建立了页岩气藏水平井体积压裂数值模型,模型中采用多孔介质流固耦合单元模拟页岩基质的行为,采用带有孔压的"cohesive"单元描述水力裂缝的性质,模型对"Texas Two-Step" 压裂方法进行了数值模拟,模拟结果得到了压裂过程中地层应力场的分布及其变化,模拟结果和解析公式计算结果吻合良好.模拟结果表明:(1)裂缝的产生减弱了地层应力场的各向异性;(2 对于低水平应力差页岩储层,采用"Texas Two-Step"压裂方法可以产生缝网. 对于采用"Texas Two-Step"压裂方法无法产生缝网的高应力差页岩储层,提出了三次应力"共振" 和四次应力"共振" 压裂方法并进行了数值模拟,模拟结果得到了压裂过程中页岩储层应力场的分布及其变化,得到了缝网形成的区域,模拟结果表明:(1)对于高应力差页岩储层,采用"Texas Two-Step" 压裂方法无法产生缝网;(2)对于高应力差页岩储层,三次应力"共振" 和四次应力"共振"压裂方法是有效的体积压裂缝网形成的方法.      

Y型微通道两相流内部流动特性刘

利用显微粒子图像测速技术、高速度数码显微系统及数值模拟方法研究了Y 型微通道内液滴的形成. 主要考虑了Y 型角度（45°,90°,135°,180°）、两相流量大小等因素的影响. 发现在挤压机制中,Y 型微通道内分散相液滴的形成主要受到来自连续相的剪切作用,Y 型角度越小,分散相所受到的剪切作用越大. 在液滴生成过程中,连续相速度剖面呈非对称抛物线型分布. 当Y 型角度小于180°时,角度的变化对液滴直径大小影响较小,但角度的减小会加快液滴的生成时间. 当Y 型角度为180°时,生成的液滴体积最大且生成时间最长. 毛细数对液滴直径和生成时间的变化同时产生影响,连续相毛细数的增大使得连续相在两相交汇位置处对分散相的作用力更集中,导致分散相更易破裂.     

页岩气藏压裂水平井试井分析

页岩气藏资源丰富,开发潜力巨大,已成为目前研究的热点.与常规气藏相比,页岩气藏运移机制复杂,流动模式呈非线性,有必要考虑页岩气的吸附解吸,天然微裂缝的应力敏感性,人工裂缝内的非达西流等非线性因素对压裂水平井压力响应的影响. 基于双重介质和离散裂缝混合模型,分别采用Langmuir等温吸附方程描述吸附解吸,渗透率指数模型描述应力敏感,Forchheimer方程描述非达西效应,建立页岩气藏压裂水平井数值试井模型. 运用伽辽金有限元法对模型进行求解.根据试井特征曲线,划分流动阶段,着重分析非线性因素对压力响应的影响.结果表明:页岩气藏压裂水平井存在压裂裂缝线性流、压裂裂缝径向流、地层线性流、系统径向流及封闭边界影响5 种流动阶段.吸附解吸的影响发生窜流之后,Langmuir吸附体积增大,拟压力导数曲线凹槽更加明显,系统径向流出现时间与压力波传播到边界时间均延迟;天然裂缝系统的应力敏感性主要影响试井曲线的晚期段,拟压力和拟压力导数曲线均表现为上翘,应力敏感效应越强,上翘幅度越大;高速非达西效应对早期段影响较大,非达西效应越强,拟压力降幅度越大,试井曲线上翘.与解析解的对比以及矿场实例验证了模型的正确性与适用性.      

EVALUATION ON STRESS INTENSITY FACTOR OF CRACK UNDER DYNAMIC LOAD USING NUMERICAL MANIFOLD METHOD

相较于传统有限元,数值流形方法（numerical manifold method, NMM） 的一个显著优点是在处理裂纹问题时网格无需与裂纹重合,这就方便了岩体破坏过程的模拟. 基于包含裂尖增强函数的NMM,采用Newmark 隐式动力学算法进行时间积分,重点研究了动力载荷条件下裂纹动态应力强度因子（dynamic stress intensity factor,DSIF） 的求解方法. 针对典型的线弹性动力裂纹问题,给出了NMM 的数值算例. 结果表明NMM 能够准确计算动载荷条件下裂纹的DSIF,并且具有较好的收敛性.