页岩气藏的双重介质-离散裂缝模型

考虑页岩气藏开发中渗流的多尺度效应,提出了一个基于裂缝-孔隙双重介质的离散裂缝模型．在该模型中,基质、天然裂缝和人工压裂裂缝采用各自控制方程独立计算,不同介质之间通过流量交换相互关联．为分析模型可靠性,分别和基于渗透率粗化及压裂裂缝导流能力无穷大的模型对比．数值算例显示,伴随着网格细分,该模型与精确渗透率粗化模型具有相同计算精度,两者收敛速度均较快,但该模型易推广到多相流动问题,而等压模型对产量将有所高估．研究了地质参数和工艺参数对气井产量的影响规律．计算结果表明天然裂缝渗透率及基质孔隙扩散系数对产气速率有着重要影响,产气速率伴随着人工压裂裂缝导流能力、长度以及数目的增加而增加,但是增加幅度会逐步趋缓．     

一种利用岩块特征长度估算二维渗吸量的新方法

对于存在大量亲水疏油岩块的天然裂缝性油藏,自发渗吸是基质岩块和裂缝系统油水交换的重要驱动机制.对于一维渗吸问题,前期存在自相似解析解,后期可以用近似解析解描述,从而可以很方便地估算出岩块剩余油饱和度和岩块出油速率.本文研究了二维矩形岩块在不同边界条件下渗吸过程的近似解析解,通过引入特征长度,将二维岩块渗吸问题等效为一维,从而可以有效利用一维近似解析解来估算出二维岩块渗吸过程中剩余油饱和度和岩块出油速率随时间的变化.本文提出了一个新的混合特征长度,在渗吸前期取岩块面积与岩块开放边界长度之比;在渗吸后期,特征长度根据前期特征长度和Ma等提出的特征长度按饱和度插值确定.通过数值检验比较了几种常见的特征长度取法和本文提出的混合特征长度取法,结果表明新的混合特征长度总体性能优于其他已有的特征长度取法,并且适用于本文列举的所有边界条件.     

一种基于边界积分方程的隐式翼型反设计算法

本文由边界元方法出发,将适用于单连通空间Laplace问题的边界积分方程推广到带环量的多连通空间中,并对离散边界积分方程中的矩阵元积分式解析化,以避免在翼型尾缘处尖点附近直接利用数值积分计算矩阵元导致的数值振荡,对于以翼型表面压力分布为收敛目标的反设计问题,利用Newton- Raphson迭代求解满足该目标压力的非线...     

可变形儒可夫斯基翼型亚音速非定常气动力的研究

对于翼面变形速度远小于来流速度情况下的儒可夫斯翼型亚音速绕流问题,通过仿射变换将可压缩流动转换成不可压缩流动,将解析解和离散涡方法相结合计算变形机翼的不可压缩流动速度场,再利用逆变换得到变形机翼的亚音速流动速度场,进而分析非定常气动力特性,建立变形机翼的准定常升力系数和非定常附加升力系数在可压缩和不可压缩两种状态下的简单近似对应关系。计算结果显示变形机翼的非定常气动升力近似等于准定常计算结果叠加上虚拟质量力导致的非定常附加升力,该非定常附加升力随翼型变形速率呈线性关系,由机翼当前时刻飞行姿态、翼型及其变形速率确定,与具体变形历史过程无关。低来流马赫数时虚拟质量力导致的非定常效应显著,高亚音速流动时准定常升力起主导作用。同时还分析了不同马赫数下机翼往复变形过程中升力的变化特性,指出尽管高亚音速变形机翼的气动升力近似等于准定常气动升力,但不能忽视非定常附加升力的影响,非定常附加升力将导致完成往复变形需要外界输入正比于Ma∞/[(1-Ma2∞)]的功。     

复杂边界稠油油藏蒸汽注采过程自适应网格法的研究

将自适应网格法推广到复杂边界稠油油藏的蒸汽注采过程,针对复杂边界附近的网格提出相应的粗化算法。首先,在实施自适应网格算法之前对边界上最精细网格上的计算参数进行预处理以提高计算精度,然后,再利用同样的预处理方法对自适应网格法中边界处的各层次粗网格的渗透率进行粗化。在建立动态AMR网格系统的网格粗化准则中,仅采用油藏温度和各相饱和度的空间变化作为控制阈值,这样边界区域在相变锋面未到达时将自适应地采用粗网格进行计算。数值算例显示边界附近自适应地采用粗网格进行计算并不影响油藏数值模拟的计算精度,自适应网格法在保持计算精度的同时,大幅度提高了计算速度。     

含泥质夹层油藏蒸汽注采过程的自适应网格算法研究

根据泥质夹层的低渗特性及空间分布，本文提出了一种含泥质夹层油藏网格渗透率的粗化计算方法，并在此基础上，将自适应网格算法应用于含泥质夹层油藏的数值模拟，提升其计算效率．在计算过程中，网格的动态划分仅依据流体物理量的变化，泥质夹层区域不全部采用细网格，仅针对流动锋面处的泥质夹层采用细网格，其余泥质夹层处采用不同程度的粗网格．相较于传统算法，网格数大幅下降．数值算例表明，自适应网格算法的计算结果精度与全精细网格一致，能够准确模拟出泥质夹层对于流体的阻碍作用，同时计算效率得到大幅提升，约为全精细网格算法的3~7 倍．     

纤维缠绕复合材料压力容器渐进损伤分析

利用大型有限元软件ANSYS的参数化设计语言(APDL)建立了纤维缠绕复合材料压力容器的有限元模型,该模型真实地反应了复合材料压力容器封头处纤维缠绕层的厚度以及纤维缠绕角度沿子午线不断变化的情况,即充分体现该压力容器结构的真实性.针对建立的有限元模型,对复合材料压力容器在一定内压下的应变进行了分析,将分析的结果与试验结果比较,验证该模型的准确性.在此基础上进行了复合材料压力容器的渐进损伤分析,获得了复合材料压力容器在外载荷逐渐增加的情况下,复合材料缠绕层逐层失效破坏的详细信息,为复合材料压力容器的设计提供一定的依据.     

非均匀介质中流动和传热问题的有限分析法

数值求解非均匀介质中的输运问题广泛应用于科学计算和工程领域．介质的强非均匀性给相关问题的准确求解带来极大的困难．近年来，本课题组将有限分析法拓展到该领域，建立了非均匀介质中输运问题的有限分析法．该算法基于网格奇点邻域内类拉普拉斯方程局部解析解构建，算法具有很高的精度，且不依赖于介质的非均匀性强度．不管相邻网格传导率差异如何，仅需对原始网格进行很少地细分就可以获得非常准确的计算结果，因此与其他传统数值算法相比，可以大幅提高计算精度和效率．该算法可广泛应用于求解非均匀多孔介质中的渗流、复合材料中的热传导及电场分布等问题．     

导弹水下点火推力峰值问题的数值研究

通过对燃烧室内燃气采用常微分控制方程,对喷管-燃气泡流场采用准一维模型,数值研究了导弹水下点火初期推力峰值与燃气泡初始体积、喷管尾盖打开压强、发射深度等参数的关系。对水流场采用势流理论并通过边界元方法求解,水流场和燃气泡流场由交界面上的压力连续和速度连续条件相互耦合;利用准一维气流场模型计算得到的结果表明导弹在水下点火后瞬间所受推力将急剧增大,与利用球形气泡模型得到的结果一致。计算结果还表明推力峰值随着燃气泡初始体积的增加迅速降低,随着喷管尾盖打开压强的增加而增大,随尾盖质量的增加而增大;在一定范围内发射深度的变化对推力峰值影响较小,计算结果可为工程设计部门提供参考。     

多连通区域的Dirichlet空间上的复合算子

对于定义在有界连通区域上的Dirichlet空间,我们着重讨论了其上复合算子的有界性、可逆性、Fredholm性及K-理论性质.