反映动态启动压力梯度的低渗透油藏渗流模型

基于流体边界层理论,建立了一种描述在低渗透油藏流体非线性渗流特征的新模型,获得了动态启动压力梯 度的连续变化函数。应用考虑动态启动压力梯度的非线性渗流模型,建立了一维单相流体流动的数学模型,并应用全 隐式方法对离散方程进行了线性化处理;对比研究了线性渗流模型、拟线性和反映动态启动压力梯度的非线性渗流模 型下地层压力的分布特征及流体非线性渗流系数对动边界扩展的影响。数值结果表明：受动态启动压力梯度的影响 在近井区域地层压力下降幅度更大,在远井区域的地层压力越高,地层压力下降范围更小;非线性系数数值越小,流体 流动的非线性程度越强,动边界的扩展速度越慢。     

隧道渗流的数值方法

用边界元方法，求解了各向异性介质的渗流场问题，计算了双孔洞隧道周围渗流区域的水压分布，并分析了渗流区域不同边界条件对水压分布的影响。     

大坝渗流的边界元分析

本文采用边界元方法计算了二维正交各向异性大坝的渗流问题，得到大坝内部压力分布，分析了在给定水深情况下坝体的稳定性，为实际的工程应用提供了理论依据．     

低渗油藏中含启动压力梯度水平井生产动态

水平井渗流为三维流动，远比垂井的渗流复杂，因而水平井的压力产量公式往往非常复杂，为了简化，进行了一些理论分析，在没有边界影响的情况下，水平井生产时在地层中形成的等压面是一个旋转椭球面。     

高压深层裂缝性油气藏中不定常渗流的数值研究

根据实验结果和裂缝性压力敏感地层的特性,建立了适合于不同外边界条件下裂缝性压力敏感地层的不定常渗流模型．推导出了裂缝性压力敏感地层的不定常渗流模型的有限元方程．在对区域进行非结构化网格自动划分的基础上,利用有限元计算方法计算出了无限大地层、圆形有界地层以及线性组合边界地层的裂缝性压力敏感地层的井底压力的理论曲线,并对理论曲线特征进行了分析和对比．确定了压力敏感系数以及有效井径组合参数对理论曲线的影响,得到了确定压力敏感油气藏的地层参数的方法．通过某气田高温高压深层气井的实例分析,得到了较好的分析结果,同时也验证该文所建立的裂缝性压力敏感地层的不定常渗流模型的正确性．     

煤层单相水非稳态渗流模型研究

为了得到排水过程中的非稳态渗流特征,通过相应的坐标变换,建立同时考虑各向异性及应力敏感现象的单相水非稳态渗流模型,采用近似解法并根据物质平衡原理,给出内边界定产、外边界无限大条件下的煤储层压力分布及动边界变化规律.算例分析表明:煤层排水过程中,压降漏斗以椭圆形式向外拓展,应力敏感性越强,井底流压和地层压力下降越快;初始渗透率和应力敏感系数的各向异性共同控制泄流面积的大小及形态.     

水平裂缝水平井的压力动态分析方法

根据渗流力学原理，提出了一种分析计算水平井伴有有限导流水平裂缝瞬变压力的方法，即在计算机软件中通过系统模拟与参数拟合的方法动态分析与研究水平裂缝水平井的压力。假设裂缝内流动的是一维流动，地层内流动是三维流动，让两个流场中解出的压力表达式在裂缝边界上相等以求出流量分布，最后求出地层的瞬态压力分布。通过计算比较了不同长宽比的水平裂缝对瞬态压力的影响，比较不同裂缝传导率对瞬态压力的影响。从导数曲线上可以看出代表有限导流水平裂缝特征的两个流动段，裂缝的不稳定流动过渡段和地层内的垂直径向流动段。     

曲线拟合坐标技术在复杂边界有压渗流数值计算中的应用

建立了曲线拟合坐标系下的坝下复杂边界渗流的数值计算模型，给出了边界条件和求解方法。计算结果表明，用该方法求得的各种物理量分布符合物理规律；与有限元法相比，该方法无需解大型方程组，节省计算内存和时间，是一种计算复杂边界渗流的有效方法。     

考虑动边界影响的低渗透油藏反九点水力压裂注采井网数值模拟研究

考虑低渗透地层中因启动压力梯度存在而产生动边界的影响,利用CMG三维三相黑油油藏数值模拟器对开发低渗透油藏的正方形反九点水力压裂注采井网进行了精确的数值模拟。数值模拟结果表明:相对于达西渗流模式等压面足够光滑、整个油藏区块均参与流动,低渗透油藏启动压力梯度为0.2 MPa/m下的等压面(油藏初始压力)空间分布有较大间断性,在远离注水井和生产井的部分区域仍保持了油藏的初始状态,该区域还未参与到油藏的渗流系统中;地层中存在随时间变化的动边界,会降低油田注水的油藏波及体积和原油采出程度。由饱和度的三维空间分布图可以看出:达西渗流模式下的压裂注水效果较好,油水分布比较集中,界面清晰,更接近于活塞式水驱油;而启动压力梯度的存在会使压裂注水井驱油效果变差,油水存在较大的油水混合渗流区域,油水分布比较分散,属于典型的非活塞式水驱油。由绘制的油田开发特征曲线可以得出:启动压力梯度越大,累计产油量越小,累计产油量随时间上升越慢,生产井含水上升越快,日产油量越低;启动压力梯度对低渗透油藏压裂注水开发效果影响显著。     

用边界元素法计算直流电测井响应

从三维稳恒电流场问题的边界积分方程的基本解出发，导出了轴对称条件下的边界元方程，得到了普通电阻率测井和侧向测井响应计算的实用算法，并计算了几种简单测井响应问题，结果表明，该方法可以用于复杂条件的直流电测井响应的计算。     

基于SBFEM的坝-库水相互作用分析

将W o lf和Song提出的比例边界有限元法(sca led boundary fin ite e lem en t m ethod,简称SBFEM)应用于坝-库水相互作用分析.在库水不可压缩假定的前提下推导了坝面动水压力与附加质量矩阵的基本方程.二维重力坝和三维拱坝坝面动水压力的算例表明,与有限元法比较,SBFEM计算精度明显提高,同时由于具有降维的特点,计算工作量也在很大程度上有所节约.还分析了库底和水面坡度变化以及水库边界形状适当变化对动水压力分布和数值大小的影响.结果表明,一般情况下水库边界的变化对动水压力的分布形状基本上不发生影响,对动水压力数值大小也影响不大.     

流动场管内阴极保护的电位分布计算

应用边界元方法计算了管内介质处于流动状态时，受到阴极保护的管道内壁的瞬时电位分布和稳定电位分布，并在实验室内进行了小型的流动场客内阴极保护实验，计算与实测结果较吻合，证实了边界元方法应用于管内阴极保护问题的有效性和准确度。     

采空区沿空通风巷道边界漏风分布的迭代计算

为了提高计算采空区漏风流动状态的有限元计算精度,提出用迭代方法处理沿空通风巷道(或工作面)边界条件,给出了迭代格式和计算方法,针对工作面的侧漏风问题提出了工作面漏风率的概念,结合算例,从理论上证明沿程漏风巷道的风压分布(1类边界条件)是非线性的,指出在漏风量占总供风量的比例不是很大的情况下,工作面风压分布基本近似于线性分布;迭代收敛的快慢与漏风量大小有关.     

三重介质油藏压力动态分析的边界元法

油藏形状对井底压力动态响应有着显著的影响,为了研究任意形状三重介质油藏的压力动态,本文首先建立了由基岩系统、裂缝系统和溶洞系统组成的任意形状的三重介质油藏试井解释模型.利用Laplace变换和边界元方法获得油藏内任意点的压力,进而由杜哈美原理得到了考虑井筒储存和表皮效应的井底压力.讨论了介质间窜流、弹性储容比、复杂边界以及油藏内存在不渗透区域等因素对压力响应的影响.分析表明,溶洞—裂缝窜流系数越大,半对数压力曲线上第一个过渡段出现得越早,双对数压力导数曲线上第一个下凹部分出现得越早;基岩—裂缝窜流系数和溶洞—裂缝窜流系数类似,只是其影响压力响应曲线的第二个过渡段;裂缝弹性储容比的大小控制着第一个过渡段持续时间的长短,表现为压力曲线上第一个过渡段台阶的宽窄和压力导数曲线上第一个凹兜的深浅;油藏的边界条件对油藏的压力动态存在明显影响,与定压边界条件相比,分段定压边界条件明显延迟了压力导数曲线下掉的时间;不渗透区域的存在使得压力导数曲线上第二个凹兜变浅,晚期径向流段抬升,偏离了0.5水平线.     

直接三维有限元求解转轮内部流场周期性边界条件的处理

本文重点介绍在以变分原理为基础的直接三维有限元求解转轮内部流场时,周期性边界条件的处理过程,以及所用的计算程序,并对 FB-8 型风机做了计算,得到了满意的速度场和压力场。     

W型通风冒落采空区流场数值模拟计算

针对冒落采空区通风比较复杂的问题，基于非均质多孔介质漏风渗流方程建立了冒落采空区漏风流态的有限元数值模型．结合现场实例，对W形通风复杂边界形式的采空区风流规律进行了求解。从理论上描绘了W形通风采空区的风压分布等值线和流函数线（分布解），该计算方案对W形通风复杂边界条件仍能满足流网正交；求解方便迅速，可操作性好。指出W形通风采空区明显降低工作面两端风压羞，缩小漏风范围的效果。     

动水压力的边界元方法

本文介绍了用边界单元法分析动水压力的原理,并且对于二维问题给出了常量边界元的支配方程以及有关系数矩阵的计算公式。实例计算的结果表明,边界单元法具有工作量少、精度高、适应性强的优点。文中还讨论了当坝体上游面倾斜、库底倾斜或水库为无限远时对作用在坝体上游面上动水压力的影响。     

考虑应力敏感和边界层效应的页岩油油藏渗流模型研究

由于页岩油储层孔隙结构复杂,原油渗流过程中存在明显的多尺度效应,而储层中存在应力敏感效应和边界层效应对原油渗流特征影响使得其流动机理更加复杂。基于毛细管模型,结合边界层和应力敏感理论,建立了页岩油储层渗流新模型并进行了验证。根据新模型分析了应力敏感效应、边界层效应对页岩油渗流特征的影响。研究结果表明：页岩油新渗流模型与渗流实验结果拟合性好,具有较高可靠性;边界层效应和应力敏感效应是影响潜江页岩油储层渗流特征的主要因素;对于含有微裂缝和溶蚀孔的页岩油储层,应力敏感效应是影响储层渗流特征的主控因素;随着压力的梯度的增加,影响渗流特征的因素由边界层和应力敏感效应的综合作用所决定。通过该研究,加深了影响页岩油油藏渗流特征相关因素的认识,为合理开发页岩油油藏提供了理论依据。     

基于三角形剖分的复杂GPR模型有限元法正演模拟

针对基于矩形网格剖分的时域有限差分法(FDTD)和有限单元法(FEM),对于物性参数分布复杂或几何特征不规则的模型适应性差的问题,从雷达波所满足的Maxwell方程出发,推导探地雷达(GPR)有限元波动方程,通过采用三角形网格剖分和线性插值基函数,在满足时间步长与空间网格差分稳定性前提下,应用Galerkin有限单元法求解GPR波波动方程;同时为消除FEM进行GPR正演模拟时来自截断边界处的超强反射,采用透射边界条件把GPR波在截断边界处的反射波透射出去,进而压制来自截断边界处的反射波。然后,编制GPR有限元正演模拟的Matlab程序。应用该程序分别对起伏分界面、"V"字形2个复杂地电模型进行FEM正演模拟,得到基于三角形网格剖分的FEM正演模拟GPR剖面图,并把该正演模拟剖面图与常规的基于矩形剖分的FEM正演模拟剖面图进行对比,结果表明:基于三角形剖分的FEM对于复杂GPR模型的物性参数分界面拟合更好,其模拟所得的正演剖面与实际模型更相符,具有更高的模拟精度,更有利于指导雷达剖面的数据解译。     

深水中复合材料椭圆柱壳声辐射的数值分析

用有限元／边界元法分析深水中复合材料椭圆柱壳的声辐射．首先利用有限元分析壳体的受迫振动，然后提取壳体表面（耦合面）数据，作为声场分析的边界条件，并计及流-固耦合作用，从而可用边界元计算壳体的辐射功率以及壳体表面、近场和远场声压，结果表明系统的耦合振动模态对于其辐射声场具有重要影响．