各向异性渗透率油藏数值模拟

建立三维三相各向异性油藏渗流数值模拟模型并应用于油藏开发.模型包含完全各向异性渗透率张量,采用十九点有限差分格式和全隐式迭代方法求解.可以描述不同区域不同主方向的各向异性渗透率分布,准确计算各向异性油藏内流体运动过程、剩余油分布和压力分布.     

各向异性油藏水平井多角度人工裂缝线性单元计算方法

考虑各向异性油藏渗透率张量的表征,利用Green函数和拉普拉斯变换建立裂缝流动的一维单元,裂缝的流量分布采用节点线性插值,裂缝内的流动处理为线性积分,耦合地层与人工裂缝的流动,建立有限导流裂缝井底压力的求解方法.结果表明:多裂缝压裂水平井存在压裂裂缝线性流、地层线性流、系统径向流3种流动形态,压裂裂缝条数越多,相同的生产时间,无因次井底压降越小;裂缝条数对流动影响明显.随着裂缝条数的增加,压降变化减小;裂缝长度和导流能力有相似的变化.人工裂缝与井筒角度越大,产能越大,当裂缝垂直于井筒时,产量最大;地层最大渗透率方向垂直于人工裂缝时产量最大,平行裂缝时产量最小.当人工裂缝垂直于井筒,并同时垂直于地层最大渗透率方向时,达到最大产量值.     

考虑渗透率张量的非均质油藏有限元数值模拟方法

针对具有混合边界的非均质油藏,考虑全张量形式的渗透率,建立弹性微可压缩单相流体不稳定渗流问题的数学模型.根据变分原理,将压力微分方程的边值问题转化为泛函的极值问题,建立渗流模型的有限元方程.针对典型的均质和非均质渗流问题进行模拟计算,得到油藏内压力动态分布曲线,并分析曲线特征.研究表明,有限元法计算精度很高,适用于求解利用渗透率张量表征的非均质油藏渗流问题.为非均质油藏的开发和精细油藏数值模拟提供了理论依据.     

裂缝性油藏水平井产能计算方法

根据裂缝性油藏水平井渗流特征,将渗流区域等效为近井区域和远井区域.远井区域视为连续均匀介质.在近井区域,基于势叠加、镜像反映原理和微元线汇思想,建立考虑天然裂缝和水平井生产段共同渗流的耦合模型.组合两个流动区域得到裂缝性油藏水平井产能评价模型,并进行求解.实例计算结果表明,由于考虑了裂缝之间的干扰,用建立的模型计算的水平井产能小于用理论公式计算的水平井的产能,与实际产能比较接近;随着裂缝密度增加,产量逐渐增加,但增加的幅度降低.     

特低渗透油藏面积井网见水时间计算

特低渗透油藏在水驱开发过程中常表现出渗透率各向异性和流体非达西渗流特征,运用非达西渗流公式及流线积分法,得到特低渗透油藏五点井网、反九点井网和菱形反九点井网在油水两相非活塞式驱替条件下的油井见水时间.计算长庆鄂尔多斯盆地某特低渗透油藏井网见水时间,结果与实际动态符合较好.分析油水黏度比、渗透率各向异性及启动压力梯度对油井见水时间的影响,结果表明：油水黏度比越大,油井的见水时间越早;当菱形反九点井网长轴方向井距与短轴方向井距之比与渗透率各向异性强度大致相等时,长轴方向与短轴方向上的角井能实现均衡驱替;启动压力梯度延缓了油井的见水时间,且生产井距越大,启动压力梯度的影响越显著.     

基于离散裂缝模型的低渗透油藏开发数值模拟

基于低渗透油藏的人工压裂开采的离散裂缝模型,考虑基质岩块中非达西渗流模式,对此类油藏的注水开发进行数值模拟研究.首先对裂缝进行降维处理,建立相应的数学模型;考虑裂缝的复杂几何形状,基于非结构化网格应用控制体积有限差分法进行数值求解,通过数例验证方法的正确性;给出两个注水开发数值算例,分析不同非达西渗流模式对开发效果的影响.     

低渗透油藏压裂水平井井筒与油藏耦合的非稳态模型

推导低渗透各向异性油藏压裂水平井井筒与盒式油藏耦合的非稳态模型,并给出求解方法.模型考虑摩阻和加速度的影响,并可以使用不同类型的约束条件.实例计算表明,压裂水平井的流动可分为非稳态阶段和拟稳态阶段.在非稳态阶段,各条裂缝的产量相差不大,总产量随着裂缝条数的增加呈线性增加;在拟稳态阶段,两端裂缝产量高于中部裂缝的产量.受摩阻和加速度压降的影响,空间上位置对称的裂缝在流量上呈现不对称性.井筒内压力损失的存在将使水平井的产量降低,并使井筒内的压力分布不均匀.在最小井底流压的基础上固定流量时,裂缝条数越多,稳产期也越长.     

水平井体积压裂缝网表征及流动耦合模型

针对岩石脆性系数高且发育天然裂缝的储层,提出表征水平井体积压裂形成裂缝网络的三种基本模式,并将渗流过程划分为油藏流动和缝网内部流动.在此基础上,利用势叠加原理导出油藏流动控制方程,利用有限差分方法建立缝网内部有限导流等式;其次,采用星三角变换法处理人工缝与天然缝的交汇流动;最后,耦合两部分流动矩阵方程得到水平井体积压裂缝网渗流数学模型.该模型表明:当水平井改造段长度一定时,压裂段数与段内分簇数是决定产能的最主要因素,其次是人工裂缝半长和人工缝导流能力,而天然裂缝密度和导流能力对产量影响较小.实例应用表明,实际产油量与模型计算值一致,误差较小.     

部分射开压裂直井渗流压力动态分析

建立符合非均质油藏部分射开地层实际的物理模型和三维各向异性矩形油藏的不稳定渗流数学模型,考虑不渗透顶、底边界和定压顶、底边界等边界条件的组合,通过无因次量纲变换、拉普拉斯变换、傅里叶余弦变换和分离变量等方法,得到拉普拉斯域的解析解,利用斯蒂芬森数值反演方法,得出实数域的压力数值解.绘制压力动态曲线,并进行敏感性分析.计算结果与数值模拟基本吻合,证实方法的可靠性.敏感性分析表明:压力动态曲线可分为早期线性流、中期径向流、晚期球形流、边界控制流四个流动期.裂缝长度主要影响早期线性流,渗透率各向异性主要影响中期径向流,储层射开程度和裂缝方位主要影响晚期球形流,边界条件和油藏宽度主要影响边界控制流.该方法可以确定最优射开程度、垂向渗透率等参数,为油藏工程分析和压裂工艺设计提供指导.     

变形双重介质分形气藏非线性渗流理论模型及数值研究

以Warren-Root模型为基础,引入分形参数和压缩系数,考虑压力对具有分形特征的渗透率和孔隙度的影响,建立变形双重介质分形气藏渗流数学模型.应用高收敛和高精度的有限元法求解这类含有第二类边界条件的数学模型.利用压力曲线,压力导数,参数对压力的影响率详细分析了部分参数对压力的影响及在变形双重介质分形气藏渗流数学模型中所具有的独特意义.     

有限导流压裂定向井耦合流动模型

考虑人工裂缝、斜井筒和地层中的耦合流动,将压裂定向井的流动分为地层向裂缝渗流、裂缝内流动和斜井筒变质量管流.采用镜像反映和势叠加原理,建立地层中势分布和渗流数学模型.应用边值理论,将人工裂缝面离散为裂缝网格,并与斜井筒变质量管流模型进行耦合,建立有限导流压裂定向井耦合流动数学模型,并形成相应的迭代算法.应用分析表明:该模型与Prats图版法相比,对于垂直缝吻合度较高;裂缝导流能力和井斜角对产量及裂缝面压力分布影响较大.     

基于不可逆热力学原理的油藏多孔介质耦合作用模型

油藏多孔介质渗流过程存在多场耦合作用,目前的一些研究方法和分析模型具有一定局限性。本文将不可逆热力学原理应用于油藏复杂渗流过程,在局域平衡和连续性假设基础上,提出了基于表征单元体(Representative Elementary Volume,简称REV)的元胞描述方法;推导出了多孔介质渗流过程考虑导热和扩散作用的耦合方程;基于所建模型,分析了稠油热采过程的熵产变化特性。稠油热采过程除包括热传导和扩散耦合作用外,还包括黏滞流动和化学反应等不可逆过程,可根据本文模型扩展。     

基于模拟有限差分法的水驱油藏渗透率时变数值模拟

针对高含水油田注入水的长期冲刷导致储层物性参数逐渐变化从而影响油田开发及剩余油分布预测的问题,提出一种基于模拟有限差分法的水驱油藏渗透率时变模拟方法。首先拟合现场油井试井渗透率,建立通过冲刷时间表征渗透率的时变模型,然后对模型采用模拟有限差分法进行离散并采用IMPES方法对其求解,最后分析不同渗透率倍数对剩余油、压力、流线分布和开发效果的影响。结果表明：考虑渗透率时变后,水窜加剧、波及系数减小、开发效果变差。与传统的油藏描述方法相比,该方法符合水驱油藏的实际情况,能够为水驱油田后期深度挖潜提供帮助。     

分形低渗透缝洞型油藏非线性渗流规律

引入分形理论,建立考虑低速非达西效应的分形三重介质缝洞型油藏数学模型,通过Laplace变换及Stehfest数值反演方法求出井底压力,借助Matlab编程绘制压力动态曲线,划分渗流阶段,分析渗流规律,进行非线性参数敏感性分析.最后结合实际算例,验证模型的正确性.结果表明：分形三重介质油藏渗流过程分为早期纯井储,过渡流,缝洞窜流,拟径向流,基质与溶洞、裂缝窜流及总体径向流6个渗流阶段;低速非达西效应对渗流的影响随时间的推移逐渐增大;启动压力梯度越大,总径向流阶段压力动态曲线上翘幅度越大;分形系数影响整个渗流过程,随着分形系数的增大,裂缝迂曲程度随之增大,致使渗流阻力增加,引起压力动态曲线整体上移幅度增大.     

二维非均匀多孔介质中不可压两相驱替的有限分析算法

为提高油藏数值模拟算法的计算效率,在求解单向稳态渗流的有限分析算法基础上,构建二维非均匀多孔介质中不可压两相渗流的有限分析算法.算法中,网格界面上的平均渗透率不是简单地取为相邻网格渗透率的调和平均值,而是通过奇点邻域解析解积分求得.相比于传统的数值算法,有限分析算法随着网格的加密,能够很快地收敛(仅需将原始网格细分至2×2或3×3),并且其计算精度和收敛性不依赖于介质的非均匀强度,从而计算效率得到提高.     

含平行裂缝储层中地震波频散、衰减、及频变各向异性*

裂缝厚度(张开度)对裂缝储层有效渗透率等起关键作用因而探测裂缝厚度对裂缝储层的勘探开发具有重要意义。针对此,该文主要介绍地震波在含平行裂缝储层中的频散、衰减及频变各向异性特征,及其与裂缝厚度的关系。重点探讨了流体运动与散射两种主要机制,并介绍了最新的理论模型研究进展。研究表明,裂缝厚度对流体运动与散射机制均会产生重大影响,主要影响频率范围为流体运动的高频区域和散射的低频区域。该频率范围常在地震频带内,故在实际地震勘探中应考虑裂缝厚度的影响,其同时为地震探测裂缝厚度提供了理论基础。最后,还探讨了流体运动与散射机制的耦合作用,并给出了相应的理论模型。结果表明,两种机制相互耦合时,地震波的频散衰减特征会发生明显变化,应予以考虑。     

使用混合网格计算非达西渗流

针对垂直裂缝井的特殊流动模式,从非达西定律出发,建立二维平面的非达西渗流方程.通过建立一组无量纲量,最终得到无量纲的渗流方程及其定解条件.假定外边界为圆形,用PEBI网格及混合网格对求解区域进行网格划分,用有限差分法对无量纲的方程进行离散,最终得到垂直裂缝井的井底压力数值解.根据此数值解并考虑井筒存储和表皮因子的影响,得到真实垂直裂缝井的井底压力.对计算结果的分析表明,使用混合网格求解非达西渗流井底压力相当准确,该方法也适用于水平井等更复杂井型及复杂边界的问题求解.     

考虑水力压裂缝和天然裂缝动态闭合的三维离散缝网数值模拟

基于嵌入式离散裂缝模型, 提出一种可在三维空间中考虑应力状态影响的裂缝动态闭合表征方法。将任意方向裂缝的开度和渗透率考虑为作用在裂缝平面法向有效应力的函数, 同时用裂缝传导率变化表征支撑剂充填的水力压裂缝与被开启的天然裂缝由于油藏开发过程中地层流体压力下降而发生的动态闭合行为。研究表明: 致密油藏开发以缝控储量为主。对压裂水平井进行产能评价时, 裂缝动态闭合会导致产能的部分损失, 其影响不可忽略; 水力压裂缝的支撑剂材料属性及天然裂缝的刚度是其中的主控因素。因此需要增大支撑剂的浓度、粒径大小并改善支撑剂的性质, 在最大程度上降低裂缝闭合对生产的不利影响。     

基于多尺度量纲统一原则的油藏孔隙分形输运特性研究

油藏孔隙分布具有分形特征,如何将分形描述方法应用到油藏多孔介质内流体的输运特性研究过程中,目前的研究方法存在一定局限,主要是没有考虑量纲统一原则和尺度效应等。本文基于量纲统一原则,利用分形描述方法,对油藏多孔介质内输运过程渗透率和渗流流量进行了推导,得到了输运特性参数的单位尺度转换系数,据此可以实现基于微观分形特征描述来分析宏观输运过程。本文建立的分形输运描述方程能够实现不同尺度层级间的转换,得到的宏观特性参数结果更接近实测值。基于本文所建模型,利用某油藏数据进行了分析验证,得到了输运流量及渗透率变化规律与实际油藏的分布规律吻合。     

缝洞型油藏致密基质灰岩的压力敏感性规律的NMR研究

相比于常规储层,缝洞型油藏灰岩储具有集体结构复杂、缝洞尺度差异大、裂缝起主要沟通作用等特点,使用常规方法研究其流动特性和流体分布规律异常困难．该文采用自行研制的核磁共振（NMR）在线驱替装置对塔河奥陶系油藏15块致密基质灰岩岩芯进行了压力敏感性实验,获得了高压及流动状态下岩芯的横向弛豫时间（T₂）谱变化情况．在此基础上分析了灰岩在压敏实验过程中微裂缝尺寸、孔隙度及渗透率的变化规律,并与常规测试分析结果进行了对比．结果表明：NMR在线检测技术可以定量测量岩芯孔隙及微裂缝的变化,为岩芯渗透率变化规律提供微观解释．