水驱油微观孔隙内流场分析

摘要目前大部分油田已经开始进入高含水期，关于如何提高原油采收率变得更加艰难。运用计算机软件模拟的方法来研究水驱油过程的运移规律是一种简单可行的方法。文章选用等温条件下广义的牛顿流体方程来建立水在微观孔隙空间内流动时的流动方程；并且运用计算流体动力学软件Polyflow对假设的模型进行计算；得到了流场的应力分布，速度分布；分析了原油粘度、流道宽度、油膜高度、压力梯度、界面张力对油膜周围应力、速度分布的影响。并且简要分析了原油粘度、流道宽度、油膜高度、压力梯度、界面张力对油膜变形的影响。结果表明：压力梯度、油膜高度以及油水之间的界面张力对驱替效果起着决定性的作用，在其他条件相同的条件下，油膜高度越大时，油膜顶部速度越大，应力越大，变形也越大；压力梯度越大时，油膜顶部速度越大，应力越大，变形也越大；界面张力越小时，油膜顶部速度越大，应力越大，变形也越大。由此可以推断，压力梯度越大，油膜高度越大，油水之间界面张力越小时，驱油效率越高。关键词：水驱油；计算流体力学；速度分布；应力分布；变形AbstractSincemostcurrentoilfieldisintothehighwatercutperiod,itiscriticalandhardtocontinuetoimprovetherecoveryofcrudeoil.Tostudywaterfloodingmigrationthroughthecomputersoftwaresimulationisafeasiblemethod.ThisarticlechoosesgeneralizedisothermalNewtonfluidconstitutiveequationstoestablishflowsflowequationofthewaterinthemicroporespace.ThecomputationalfluiddynamicssoftwarePolyflowisusingtomodelnumericalcalculation.Thismodelisgivenfortheboundaryconditions.Wegottheflowfield,thevelocitydistributionofthestressdistribution,analyzedtheinfluenceofcrudeoilviscosity,oilfilmheight,passagewidthandinterfacialtension,pressuregradientonandspeedthestressdistribution,simplyanalyzedtheinfluenceofcrudeoilviscosity,oilfilmheight,passagewidthandinterfacialtension,pressuregradientonoilfilmdeformation.Theresultshowthat,thespeedreachmaximumintheoilfilmplace,themaximumspeedoftheimportislessthanthatofexport,thisiscausedbytheinfluenceoftheflowofoilfilmontheexport.Underthesamecondition,thegreaterthepressuregradientis,themoreobviousoilfilmdrivingeffectis.Thenarrowertheflowrunnerofoilis,themoreobvioustheeffectis.Thehigheroilfilmheightis,themoreobvioustheeffectis.Thesmallertheinterfacialtensionbetweenoilandwateris,themoreobvioustheeffectis.Thelowercrudeoilviscosityis,themoreobvioustheeffectis.Thegreaterthepressuregradientis,thegreatertheextremestressis.Thenarrowertheflowrunnerofoilis,hegreatertheextremestressis.Thehigheroilfilmheightis,hegreatertheextremestressis.Thesmallertheinterfacialtensionbetweenoilandwateris,thegreatertheextremestressis.Thelowercrudeoilviscosityis,hegreatertheextremestressis.Pressuregradient,theheightofoilfilmandtheinterfacialtensionbetweenoilandwaterhaveamajorinfluenceondeformation.Keywords:waterflooding;computationalfluiddynamics;distributionofspeed;distributionofstress;deformation第1章概述1.1水驱油微观孔隙内流场分析的目的及意义随着经济的快速发展,中国对石油的依赖日益增长,这一趋势在这个世纪将继续下去。不断促进对外石油贸易的同时,必须大力加强有效开采国内石油资源,以确保更多的石油生产来满足中国的重要战略目标的需要。由于目前大多数油田已进入高含水阶段,如何进一步提高石油采收率成为至关重要的问题。通过软件模拟方法研究水驱油的方法是一种可行的方法。使用计算流体动力学软件,运用计算机超强的数值运算能力进行模拟,从而进一步研究石油开采运移过程中的一些规律。因此,水驱油物理过程流场对石油工程师提高原油采收率具有重要意义。石油开采是不可重复的,更确切的了解一个地区的地下水驱油的流场情况对石油工程师是非常重要的,同时也对石油生产具有重要的指导意义,通过研究的目的可以掌握水驱油过程的运移规律,压力分布,并通过比较不同地层条件情况,深化研究水驱油的流场情况,帮助石油工程师决策石油生产,提高生产效率,节约生产成本。1.2水驱油微观孔隙内流场分析的背景石油的开采过程分为自然方法采油(一次采油)、水驱油(二次采油)、增强水驱油(三次采油)。在1940年以前，油田开发的技术水平较低，主要依靠消耗天然能量开采。一般采收率为5%到10%，称之为一次采油。渗流理论发展到一定程度后，达西定律逐渐地运用到了流体在多孔介质中的渗流，表明油井的产量与油井的压力梯度成正比。这使人们意识到一次采油过程中低原油采收率的主要原因是油层能量的衰竭，并在这个基础上提出了通过注入人工能量来保持油层压力从而开发油田的二次采油法。但是经过二次采油后仍有60%—70%剩余油残留在地下采不出来。通过多种方法预测显示我国油田的平均采收率只有34.2%，有部分油田甚至只有20%—25%。因此，多年来，国内外的石油工作者针对如何提高原油采收率这个问题进行了大量研究，从而提出了三次采油法。即通过增强水驱油方法来提高原油的采收率，主要包括化学法、混相法、热力法和微生物法等。由于作用原理不同，化学法又可进一步分为聚合物驱、碱驱、表面活性剂驱、碱一聚合物复合驱、碱一表面活性剂一聚合物复合驱以及表面活性剂一碱一聚合物复合驱。由于所用的混相剂不同，混相法分为溶剂混相驱、烃混相驱、CO2混相驱和N2混相驱。化学法和混相法主要用来开采轻质油，而粘度大于50mPa·S的重质油的开采主要采用热力法，热力法又可分为蒸汽吞吐、蒸汽驱和火烧油层。泡沫驱又属于化学驱，但是通常又用于混相驱和蒸汽驱作为流度控制技术。微生物采油法根据注入工艺又可以分为微生物单井吞吐和微生物驱。另外，注浓硫酸工艺以及用微波、电磁波、磁场等进行井底带处理工艺等提高石油采收率方法也都在一些文献中有报道，但目前尚处于探索阶段，需要进一步研究和发展。1.2.1碱水驱油在石油工业中，传统的注水方式提高原油采收率的效果不明显，水的粘度较低，有效渗透率高，导致流度比偏大，水突破较早。同时，由于水和油的化学性质不同，使得两相不混容，导致界面张力高、原油采收率低。因此，人们上尝试了各种化学的、物理的及热采的方法，并且根据油藏和油层中流体的特性选择其中一种或几种方法结合使用。在所用各种化学试剂中，已经考虑到将碱加入水驱的水中，用于清除孔壁中的油，使孔隙介质中机械圈闭的油滴流动。由于操作简单，价格经济实惠，在提高原油采收率方面，碱水驱将继续受到关注。另外，经研究，单纯的冷碱水驱只能采出大约57%的地层油，但是若使用碱蒸气驱则可以将采收率提高到80%。但是，热碱不仅会在酸性油中产生表面活性剂，而且还导致地层原油的活度降低。同时由于碱能够与地层岩石发生反应，导致储层损害，所以碱驱油有被淘汰的可能。1.2.2聚合物驱油聚合物驱是指在注入地层的水中加入水溶性的高相对分子量的聚合物，增加水相的粘度，降低水相的渗透率，改善油水粘度比，从而扩大水驱油的体积波及系数，从而提高原油的采收率。当油水粘度比很大时，采出的液体中含水率上升速度很快，油水粘度比较小时，采出液中含水率的上升速度将逐渐缓慢。聚合物驱油剂的一个重要特点就是它能够在石油开发的控制条件下，通过较低的聚合物浓度来保持高的粘度。目前油田常见的聚合物是聚丙烯酞胺以及生物聚合物。但水通常为假塑性流体，粘度随剪切速度的变化而变化，而且很有可能堵塞地层通道，所以在很多的油田并不适合使用。1.2.3三元复合驱三元复合驱简称ASP，是有碱、聚合物和表面活性剂组成。各组分之间相互补充相互增效从而达到提高采收率的目的，研究表明，这种方式能够使采收率达到80%，因此成为当前研究的热点。尽管增强水驱的方法的驱替效率较高，但是由于其成本也较高，而且由于大多数油田的开采年限较短，所以在实际的油田中，使用普通的水驱油开采方式的比较多，因此本文的模拟系统也是基于普通水驱油的基础上进行的。据预测，在21世纪，原油需求总数量为2500到2600亿吨，按照目前油藏开发的技术和措施年均只能提供380亿吨。因此，必须将采收率提高到65%一70%才能满足需求。提高油气产量的途径主要有两种，即使用高效的生产技术和增加石油地质储量。当前提高原油采收率已经成为石油工程师们面临的最重要的挑战之一。首先较低的波及系数和洗油效率使得油藏中有2/3的地质储量采不出来。另外随着油田开采时间增加，地层中的水形成了通路，使得驱油效率快速下降。油田开发分析法、物理模拟实验法及数值模拟法是油田研究高含水期微观水驱剩余油分布及挖潜措施的重要方法。油田开发分析法主要是通过研究油田的动、静态资料来掌握油田的变化规律及开发状况。物理模拟实验方法是通过理论计算与实际资料对比，找出油田开发中的主要矛盾及开发与管理中存在的主要问题，从而有针对性地采取措施。由于三维油藏系统的复杂性，同时同一油藏又不能从油藏原始条件下重新进行二次注水，所有的不利因素导致我们弄不清开发时各种变量所产生的影响。运用模拟和数学方法来计算地下注水动态只能得出被大大简化了的情况的计算结果。应用物理模拟技术，可以较客观地描述地下流体的运动和变化规律，并定量地描述剩余油饱和度。1.3国内文献综述杨清彦、兰玉波、贾忠伟、张洪兴应用微观仿真模型进行水驱油实验，并且通过微观驱油动态的彩色图像量化处理系统对实验的过程和结果进行分析，并且从定量和定性两个角度对水驱油过程进行了分析。水驱油时，由于润湿性不同，注入水在微观空隙中有不同的渗流机理。在亲油岩石的孔隙中，注入水沿着孔隙中心突入岩石；在亲水岩石孔隙中，注入水附着在岩石空隙壁上爬入。由于储层岩石具有非均质性，剩余油在空隙中又可以分为分散型剩余油和连片状剩余油两类。具体形式有簇状、角状、柱状、膜状和“孤岛”状等分布形态。另外不同的化学驱替对于不同的剩余油形式的作用效果不同。李星民，殷广明，殷茵从多孔介质孔隙分布特征与水驱油微观机理出发，针对具单峰对数正态分布特征的孔隙结构，给出了一种可以定量描述油水在原始状态下以及水驱过程中的微观分布及其变化的方法。该方法可较真实地揭示出水驱油过程中波及系数与油水微观分布及其变化的关系，并可以推广应用到具有双峰甚至多峰分布特征的孔隙结构中[1~4]。何顺利，李中锋，门成全，杨文新研究制作了能够真实反映出岩石微观孔隙结构的仿真模型，并且运用了图像分析技术和微观渗流物理实验模拟技术对水驱油过程及残余油的分布进行观察，模拟了不同驱替速度、孔隙结构、粘度比等环境下的水驱油过程。建立了用分形维数定量表征残余油及空间分布的测定方法。研究表明，驱替速度越大，形成的残余油量越少；原油粘度越大，残余油量越多[5~7]。姜汉桥，陈民锋二人模拟储层岩心的一次驱油过程和二次驱油过程，建立了油水两相三维准静态孔隙模型；通过将试验结果与模型计算对比，验证了该模型的有效性。利用建立的孔隙网络模型，研究岩心在亲水条件下，储层微观孔隙结构参数对水驱油驱替特征的影响。研究表明，水驱油的驱替效果随着储层岩心的孔喉比变小而变好，随着空隙配位数和形状因子增大而变好，随着残余油饱和度越的变小而变好，因此形状因子对水驱油效率影响很大[8~11]。第2章计算流体动力学基础2.1计算流体动力学概述2.1.1什么是计算流体动力学计算流体力学或计算流体动力学，英文ComputationalFluidDynamics，简称CFD,是用电子计算机和离散化的数值方法对流体力学问题进行数值模拟和分析的一个分支。CFD的基本原理可以概括为：把原来在时间域及空间域上连续的物理量的场，如压力场和速度场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。CFD可以认为是流动的基本方程(质量守恒方程、能量守恒方程、动量守恒方程)控制下对流动的数值模拟。通过数值模拟，我们可以得到非常复杂问题的流场内各个位置上的基本物理量(如速度、温度、压力、浓度等)的分布，以及这些物理量随着时间变化的变化情况，确定漩涡分布特征、空化特征及脱流区等。CFD方法与传统的理论分析方法、实验测量方法构成了研究流体流动问题的一个完整体系。理论分析方法的优点在于其所得的结果具有普遍性，各种影响因素清晰可见，是指导实验研究和验证新的数值计算方法的理论基础。然而，它通常需要对计算的对象进行抽象和简化，才有可能得出理论解。对于非线性的情况，只有少数流动才能给出解析结果。实验测量方法所得到的实验结果是真实可信的，它是理论分析方法和数值方法的基础，其重要性不容小觑。然而，实验往往受到模型尺寸、人身安全、流场扰动和测量精度的限制，有时可能很难通过试验方法得到结果。此外，实验还会遇到经费投入、人力和物力的巨大消耗及周期长等许多困难。而CFD方法恰好克服了前面两种方法的弱点，在计算机上实现一个特定的计算，就好象在计算机上做一次物理实验。例如，机翼的绕流，通过计算机并将其结果在屏幕上显示，就可以看到流场的各种细节；如激波的运动、强度，涡的生成与传播，流动的分离、表面的压力分布、受力大小及其随时间的变化等。数值模拟可以形象地再现流动情景，与做实验没有什么区别。2.1.2计算流体动力学的特点CFD的长处是适应性强、应用面广。首先，流动问题的控制方程一般是非线性的，自变量多，计算域的几何形状和边界条件复杂，很难求得解析解，而用CFD方法则有可能找出满足工程需要的数值解；其次，可利用计算机进行各种数值实验，例如，选择不同流动参数进行物理方程中各项有效性和敏感性试验，从而进行方案比较。再者，它不受物理模型和实验模型的限制，省钱省时，有较多的灵活性，能给出详细完整的资料，很容易模拟特殊尺寸、高温、有毒等真实条件和实验中只能接近而无法达到的理想条件。CFD也存在一定的局限性。首先数值解法是一种离散近似的计算方法，依赖于物理上合理、数学上适用、适合于在计算机上进行计算的离散的有限的数学模型，且最终结果不能提供任何形式的解析表达式，只是有限个离散点上的数值解，并有一定的计算误差；第二，它不像物理模型实验一开始就能给出流动现象并定性描述，往往需要由原体观测或物理模型试验提供某些流动参数，并需要对建立的数学模型进行验证；第三，程序的编制及资料的收集、整理与正确利用，在很大程度上依赖于经验与技巧。此外，因数值处理方法等原因有可能导致计算结果的不真实，例如产生数值黏性和频散等伪物理效应。当然，某些缺点和局限性可通过某种方式克服或弥补。另外，CFD因涉及大量的数值计算需要较高的计算机软硬件配置。2.1.3计算流体动力学的应用领域近年来，CFD方法有了飞速的发展，一些经典流体力学的一些近似计算法和图解法已经被替代；过去的一些典型教学试验现在完全可以借助CFD手段在计算机上实现。所有涉及流体流动、热交换、分子输运等现象的问题，都可以通过计算流体力学的方法进行分析和模拟。CFD不仅作为一个研究工具，而且还作为设计工具在水利工程、土木工程、环境工程、食品工程、海洋结构工程、工业制造等领域发挥作用。典型的应用场合及相关的工程问题包括：（1）泵、水轮机和风机等流体机械内部的流体流动；（2）飞机和航天飞机等飞行器的设计；（3）汽车流线外形对性能的影响；（4）洪水波及河口潮流计算；（5）风载荷对高层建筑物稳定性及结构性能的影响；（6）温室及室内的空气流动及环境分析；（7）电子元器件的冷却；（8）换热器性能分析及换热器片形状的选取；（9）河流中污染物的扩散；（10）汽车尾气对街道环境的污染；（11）食品中细菌的运移。对这些问题的处理，过去主要借助于基本的理论分析和大量的物理模型实验，而现在大多采用CFD的方式加以分析和解决，CFD技术现已发展到完全可以分析三维黏性湍流及漩涡运动等复杂问题的程度。

2.2Polyflow软件概述2.2.1Polyflow软件简介POLYFLOW是基于有限元法的计算机流体动力学(CFD)软件，专用于黏弹性材料的流动模拟，1982年由比利时Louvia大学研发，1988年成立了公司，1997年被世界著名的流体分析软件公司FLUENT公司收购。POLYFLOW具有强大的解决非牛顿流体以及非线性问题的能力，可以解决聚合物、食品、玻璃等加工过程中遇到的等温/非等温、二维/三维、稳态/非稳态的流动问题。并且POLYFLOW软件的升级和非牛顿流体力学数值模拟技术的发展密切相关，POLYFLOW软件继承了当今最新、最完善而适当的运算法则，它的本构方程也是基于该领域最新的研究成果。与其他CFD软件相比，POLYFLOW具有如下优点：(1)POLYFLOW采用变形网格、接触算法、网格重叠技术，提高了计算的准确性；(2)POLYFLOW与GAMBIT、IDEAS、PATRAN有良好的数据接口，可以使用他们生成网格，同样支持多种网格类型，比如三角形、四边形、四面体、五面体、六面体、组合网格等；(3)POLYFLOW在模拟复杂的流变特性流体或者黏弹性流体时，主要有广义牛顿模型、屈服应力模型和黏弹性模型3种模型。正因为如此，POLYFLOW使用与聚合物材料的挤出成型(包括直接挤压、反向挤压、混合挤压等)、吹塑成型、模压成型、热成型、涂复成型等各种成型工艺的模拟计算，还可应用在热传输、纺丝、共混、渗流、玻璃熔炉中的流动、轮胎面的制造与设计等。2.2.2Polyflow软件的结构POLYFLOW软件包括POLYMAT、GAMBIT、POLYSTAT、POLYDATA、FLUENT/Post(FIELVIEW)5个主要模块，它们由一个主控程序POLYMAN来执行。其中POLYMAT模块用于创建黏弹性材料库，对实验数据根据用户指定的模型进行拟合。GAMBIT是前处理模块，用于创建划分网格、几何模型、设定边界和子域。POLYDATA是数据模块，用来设定物理模型、计算类型、边界条件、材料参数、求解方程和数值参数等。POLYFLOW实际上是一个求解器，求解POLYDATA设置完成后输出的数据文件，生成后处理模块可以读取的结果文件。FLUENT/Post是后处理模块，可以把POLYFLOW求解生的结果文件通过图像的形式显示出来，以便用户有效地观察和分析流动计算结果。POLYSTAT模块则用于分析加工流场的共混效果。由于Polyflow软件独有的特点，本文将采用此软件进行流场计算。

2.2.3Polyflow软件的求解过程图2-1给出了Polyflow的计算框图，其求解过程分为以下几步：否是否是建立控制方程确立初始条件和边界条件划分计算网格，生成计算节点建立离散方程离散初始条件和边界条件给定求解控制参数求解离散方程显示和输出结果解收敛否？图2-1Polyflow的计算框图2.3流体动力学控制方程流体流动要受物理守恒定律的支配，基本的守恒定律包括：质量守恒定律、动量守恒定律、能量守恒定律。如果流动处于湍流状态，系统还要遵守附加的湍流输运方程[12~16]。

2.3.1质量守恒方程任何流体流动都必须遵守质量守恒定律。该定律可以表述为：单位时间内流体微元体中质量的增加量等于同一时间间隔内流入该微元体的净质量[17~19]。通过这个定律，可以得出质量守恒方程：QUOTE∂ρ∂t+∂(ρu)∂x+代入矢量符号QUOTEdiva=∂ax∂x+∂ay∂y(2-2)式中ρ——密度，kg/m3；T——时间，s；u、v、w——速度矢量u在x、y、z方向的分量，m/s。式(2-2)为瞬态三维可压流体的质量守恒方程。若流体不可压，密度ρ为常数，则其变为：(2-3)若流动处于稳态，则密度ρ不随时间变化，则式(2-3)变为：QUOTE∂(ρu)∂x+∂(pv)∂2.3.2动量守恒方程动量守恒定律也是任何流动系统都必须满足的基本定律。该定律可表述为：微元体中流体的动量对时间的变化率等于外界作用在这个微元体上的各种力的和。根据这一定律，可导出x、y、z三个方向的动量守恒方程：(2-5)QUOTE∂(ρu)∂t+divρvu(2-7)式中p——流体微元体上的压力，Pa；τxx、τxy、τxz——因分子黏性作用而产生的作用在微元体表面上的黏性应力τ的分量，Pa；Fx、Fy、Fz——微元体上的体力，N。若体力只有重力且z轴竖直向上，则Fx=0、Fy=0、Fz=–ρg。2.3.3能量守恒方程能量守恒定律是包含有热交换的流动系统必须满足的基本定律。该定律可描述为：微元体中能量的增加率等于进入微元体的净热流量加上体力与面力对微元体所做的功[17~19]。流体的能量E通常是内能i，动能K＝（u2+v2+w2）/2和势能P三项之和，我们可针对总能量E建立能量守恒方程。但是，这样得到的能量守恒方程并不是很好用，一般是从中扣除动能的变化，从而得到关于内能i的守恒方程。而我们知道，内能i与温度T之间存在一定关系，即i=cpT，这样，可得到以温度T为变量的能量守恒方程：QUOTE∂(ρT)∂t+divρu式中Cp——比热容，J/(kg·K)；T——温度，K；K——流体的传热系数，W/(m2·K)；St——流体的内热源及由于黏性作用流体机械能转化为热能的部分，J。2.3.4组分质量守恒方程对于一个确定的系统而言，组分质量守恒定律可表述为：系统内某种化学组分对时间的变化率，等于系统界面净扩散流量与通过化学反应产生的该组分的生产率之和。根据组分质量守恒定律，可写出组分s的组分质量守恒方程：QUOTE∂(ρcx)∂式中Cs——组分的体积浓度，L/L；ρCs——该组分的质量浓度，kg/L；Ds——该组分的扩散系数，m2/s；Ss——系统内部单位时间内单位体积通过化学反应产生的该组分的质量，kg。各组分质量守恒方程之和就是连续性方程ΣSs=0。因此，如果共有z个组分，那么只有z－1个独立的组分质量守恒方程。将组分守恒方程各项展开，上式可改写为：QUOTE∂ρcs∂t+∂ρcsu∂组分质量守恒方程简称组分方程。一种组分的质量守恒方程事实上是一个浓度传输方程。因此，组分方程在有些情况下称浓度传输方程。2.3.5控制方程通用形式为了便于对各种控制方程进行分析，并用同一程序对各控制方程进行求解，建立各基本控制方程的通用形式是比较方便的。用φ表示通用变量，则此通用方程表示为：QUOTE∂(ρϕ)∂t+divρu其展开形式为：QUOTE∂(ρϕ)∂t+∂(ρ式中ϕ——通用变量；u、v、w——求解变量；Γ——广义扩散系数；S——广义源项。对于特定的方程，ϕ、Γ、S具有特定的形式，表2-1给出了三个符号与各特定方程的对应关系。表2-1通用控制方程中各符号的具体形式ϕΓS连续方程100动量方程uiμ能量方程TST组分方程cxDxSs

第3章两相界面理论两相之间分界称界面，包括气—液，气—固，液—液，液—固，固—固等界面。其中液—气界面，固—气界面，也称表面。界面有一定厚度(几个分子到几十个分子的厚度)，并非几何面。在此区域内，可称为界面相(性质与相邻两相内部截然不同)。物质表面层分子与内部分子所具有情况(能量、作用力)不相同的。物质内部的分子其周围均有同类分子包围着，平均地看，所受用围分子的引力是球形对称的，各个方向的力被此互相抵消，合力为零，所以物质内部的分子可以任意移动，而不消耗功。而处在表面上的分子则不同，一方面受到体相分子的作用，另一方面又受到性质不同的另一相中物质分子的作用，使得表面层分子受到不对称的作用力。如将内部分子移至表面，必须对所移动的分子施加外力、做功，表面层的分子受到指向物体内部并垂直于表面的作用力(合力)，使物体表面有自动缩小的趋势，表面能量较物体内部大，因而产生各种表面现象[20~23]。3.1界面张力由于表面层的分子受到的作用力是指向物体内部并垂直于表面，使物体表面有自动缩小的趋势，如图3-1所示。对于液体，其内部的分子对表面层中分子吸引力，使表面层中的分子恒受到指向液体内部的拉力。图3-1图3-1分子间的作用力在两相界面存在着一种张力，它垂直于表面的单位线段上，并指向使表面紧缩的方向并与之相切。作用在单位边界线上的力就称为表面张力γ，单位是N/m。对于两相界面来说，表面张力只是自由表面能的一种表示方法，实际上两相界面上并不存在着什么“张力”。只有在三相周界上，该表面能才呈现出来表面张力的作用，即只有在三相周界上，表面能才以“张力”的形式表现出来。对于表面张力，作用点在表面上，若为弯曲面则在切面，使表面积收缩体现了方向；其大小是单位长度上的力。影响界面张力的因素有物质的本性，还有温度、压力等。3.1.1物质的本性两相接触时，其界面层自由表面能的大小直接与两相分子的性质有关；两相间分子的极性差越大，相间分子力场的不平衡就越严重，表面能也就越大。非极性物质与非极性物质之间、极性物质与极性物质彼此更容易吸引，处于界面上的分子之间差异更小，界面张力也更小。3.1.2温度的影响大多数液体物质的表面张力随温度的上升呈线性下降趋势。但在临界温度以前(距临界温度30K以内)，有明显偏差。一般线性经验关系最简的形式是：(3-1)式中b——温度系数，K-1；T——热力学温度，K。此外还有经验关系式：QUOTEγ=γ0(1-TTc)式中Tc——液体物质的临界温度，K；N——经验系数。对与有机体n的平均值为1.21。当温度为Tc时，表面张力趋于零。还可以用多项式来表达γ-T之间的函数关系。例如水的表面张力和温度的关系为：γ=75.796-0.145t-0.00024t2(3-3)式中t——温度，℃。该式适用的温度范围为10～60℃。3.1.3由热力学基本关系式dG=–SdT+Vdp+γdA，在T一定时，根据全微分的性质，有：QUOTE(∂γ∂p)A，T=(因为表面层物质密度低于液体体相密度，所以一般dV为正，即表面张力随压力增加而增加。但实际情况相反，表面张力随压力增加而减小。通常某种液体的表面张力为该液体与含有液体的蒸汽的空气相接触时的值。因为在一定温度下液体的蒸汽压是定值，因此只能靠加惰性气体等方法或改变气相中空气的压力来改变气相的压力。气相压力增加，气相中物质在液体中的溶解度增加，并可能产生吸附，会使表面张力下降。压力的影响与气体的溶解及吸附的影响相比，后者更甚。所以，液体表面张力一般随气相压力增加而降低。如20℃时，水在0.098MPa压力下表面张力为72.8mN/m，在9.8MPa压力下为66.4mN/m[24]。3.2弯曲液面的附加压力3.2.1弯曲表面上的附加压力在一个大的容器中，静止的液体表面一般是一个平面，但在某些特殊情况下，由于液体和固体间的相互润湿，液体会沿固体表面延展，使得液气表面是一个弯曲表面(如图3-2)[25]。图3-2图3-2弯曲液面的附加压力由于表面张力的作用，弯曲液面上的表面张力不是水平面而是沿界面处与表面相切。对于凸面，表面张力将有一指向液体内部的合力，液体内部压力大于外部压力，使它受到一个附加压力。用ps表示附加压力，ps=p内－p外。凹面正好相反，液体内部的压力小于外面的压力，也受到一个附加压力。附加压力总是指向液面的曲率中心。使凹面一侧的压力比凸面一侧的高[26]。1805年Young-Laplace导出了附加压力与曲率半径之间的关系式：QUOTEps=σ1R1+1R式中Ps——曲面的附加压力，Pa；σ——两相间界面张力，N/m；R1，R2——任意曲面的两个主曲率半径，m。根据数学上规定，凸面的曲率半径取正值，凹面的曲率半径取负值。所以，凸面的附加压力指向液体，凹面的附加压力指向气体。当曲面为球面时，两个曲率半径等同，此时拉普拉斯方程变为：QUOTEps=2σr(3-6式中r——球面半径，m。对于平面，两曲率半径R1，R2为无限大，Ps＝0。3.2.2如果将一根毛细管插入润湿相液体中，则管内液气界面为凹形，那么液体受到一个附加向上的压力，使湿相液面上升一定的高度；反之，如果把毛细管插入到非润湿相中，.则管内液体界面呈凸形，液体受到一个向下的附加压力，使非润湿相液面下降一定的高度。这种在毛细管中产生的液面上升或下降的曲面附加压力称为毛细管压力，简称毛管力。由于附加压力所引起的毛细管中液体上升或下降的现象就称为毛细现象(如图3-3)[27]。图3-3图3-3由于毛细管力引起的液体升降由力的平衡得到：QUOTEps=Δρgh(该式表明：液柱上升(或下降)的高度直接与毛管力ps有关，毛管力越大，液柱上升(或下降)越高。对实际油层来说，毛管力的作用是不容忽视的。当岩石亲水时，毛管力为正，h也为正，水面就会上升，如将亲水岩心浸泡在油中，就会发现水能在毛管力作用下自动进入岩心，驱出了岩心中的非湿相油这一过程称为吸吮过程或自吸过程；反之表现出岩心不能自动吸水，毛管力成为阻力，若要使水进入岩心，则必须施加一个外力克服毛管力，才能使水驱油，此为驱替过程。毛管数是一个无量纲量，代表粘性力与界面张力的比值。用Ca表示。QUOTECa=μUγ(3-8)式中μ——水的粘度，Pa·s；U——水的特征速度，m/s；Γ——两相的界面张力，N/m。粘性力的作用使油膜变形，而界面张力使油膜恢复原来的形状，阻止变形。毛管数就反映了这两种力的力量对比，当毛管数较大时，粘性力占主导地位，油膜变形，变形到一定程度时，就有可能断裂；当毛管数较小时，界面张力占主导地位，油膜恢复原来的形状，不会破裂。经实验证明：如果雷诺数Re较小时，在稳定剪切流中，即使粘性力再大，也不会使液滴破裂。但是，一般情况下，随着毛管数的增大，油膜的变形就越大，从而可以增大洗油效率[28~30]。3.3润湿液体润湿与不润湿固体，统称为润湿现象。润湿现象是固体与液体接触处的表面现象，对同一种液体能润湿某些固体，但是对另一固体则不能润湿。地下油藏也是如此。对于某些亲油岩石，油附着在表面，而对于一些亲水的岩石，油则以球形或椭球形在岩石表面存在。正因为如此，针对不同情况就需要采用不同的方案进行驱油。对于亲水岩石表面的油，可以在岩石表面滚动，与岩石的附着力小，所以很容易被驱走；而对于亲油岩石表面的油，以油膜形式附着在岩石上，很难被驱走。这种亲油岩石表面的油膜就是本文要研究的问题。水驱已经不能将其驱走，所以采用聚合物驱，利用水的高的粘弹特性，将这些油膜逐渐剥离、断裂，形成油滴被溶液带走。如将一小块岩石放入水银中，可以发现水银并不在岩块表面上铺展。相反将岩块放入水中，岩石表面就被水所润湿，此时，岩石表面与水接触，设它们的界面张力为σws。在岩块投入水之前，岩块处于空气中时，设其表面张力为σgs。投入水的前后，岩块表面积A没有发生改变，但表面能变化了，设表面能的改变为ΔU，则：QUOTEΔU=σws⋅A-σgs⋅通常σgs>σws，表面能的改变ΔU<0为负值，说明表面能降低。第4章数值计算及结果分析4.1模型的描述目前，经水驱仍然存在于地下的残余油还占有很大的比例。这些残余油可划分为以下几类：(1)油滴；(2)油柱；(3)油膜；(4)盲端油；(5)簇状油本文主要考虑亲油岩石表面的残余油膜。由于油膜的尺寸为微米级的，重力场的影响就忽略不计。通过使用计算流体软件Polyflow，对其流场进行分析，进而得到油膜的受力情况。4.2速度场模拟图4-1~图4-10分别给出了不同压力梯度下、不同流道宽度、不同油膜高度、不同界面张力、不同原油粘度的流场速度等值线和速度分布云图。4.2.1压力梯度对于速度的影响图4-1图4-1不同压力梯度情况下的速度等值线图4-2图4-2不同压力梯度情况下速度分布云图图中由上到下分别为压力梯度为0.0002MPa、0.002MPa、0.02MPa的速度等值线和速度云图。由图4-1和图4-2可知，随着压力梯度的增加，流道中流体的流动速度明显增加。压力梯度为0.02MPa/m时，油膜顶部出现明显的流动，说明驱油效果明显，并且压力梯度越大越明显。4.2.2流道宽度对于速度的影响图4-3图4-3不同流道宽度情况下的速度等值线图4-4图4-4不同流道宽度情况下的速度分布云图图中由上到下分别为流道宽度为15μm、20μm、30μm速度等值线和速度云图。由图4-3和图4-4可知，流道宽度为15μm时，油膜顶部出现明显的流动，说明驱油效果比较明显，并且流道越窄越明显。4.2.3油膜高度对于速度的影响图4-5图4-5不同油膜高度情况下的速度等值线图4-6图4-6不同油膜高度情况下的速度分布云图图中由上到下分别为油膜高度为3μm、4μm、5μm速度等值线和速度云图。由图4-5和4-6可知，油膜高度为5μm时，油膜顶部出现明显的流动，说明驱油效果比较明显，并且油膜越高越明显。4.2.4界面张力对于速度的影响图4-7图4-7不同界面张力情况下的速度等值线图4-8图4-8不同界面张力情况下的速度分布云图图中由上到下分别为界面张力为0.001mN/m、0.01mN/m、0.1mN/m速度等值线和速度云图。由图4-7和4-8可知，界面张力为0.001mN/m时，油膜顶部出现明显的流动，说明驱油效果比较明显，并且界面张力越小越明显。4.2.5油膜粘度对于速度的影响图4-9图4-9不同油膜粘度情况下的速度等值线图4-10图4-10不同油膜粘度情况下的速度分布云图图中由上到下分别为油膜粘度为0.015Pa·s、0.02Pa·s、0.03Pa·s的速度等值线和速度分布云图。由图4-9和图4-10可知，原油粘度为0.015Pa·s时，油膜顶部出现明显的流动，说明驱油效果比较明显，并且原油粘度越小越明显。比较分析图4-1~图4-10可以得出以下结论：在油膜处速度达到最大值，进口最大速度小于出口最大速度，这是由于出口流动受到油膜的影响。在其他情况一样的条件下，压力梯度越大，油膜被驱动效果越明显；流道越窄，油膜被驱动效果越明显；油膜高度越高，油膜被驱动效果越明显；原油与水的界面张力小，油膜被驱动的效果越明显；原油粘度越低，油膜被驱动的效果越明显。4.3应力场模拟图4-11~图4-25分别给出了不同压力梯度下、不同流道宽度、不同油膜高度、不同界面张力、不同油膜粘度流场的第一法向应力、第二法向应力、切向应力

4.3.1压力梯度对于应力的影响图4-11图4-11不同压力梯度情况下的第一法向应力图4-12图4-12不同压力梯度情况下的第二法向应力图中由上到下分别为压力梯度为0.0002MPa、0.002MPa、0.02MPa的τ11，τ12，τ22的等值线图。从图4-11～图4-13中可以看出：随着压力梯度的增加，流场中τ11，τ12，τ22的极值都增大。即压力梯度越大，流场中的应力极值就越大，但是并不是随着压力梯度的成倍增加而成倍增加。图4-13图4-13不同压力梯度情况下的切向应力4.3.2流道宽度对于应力的影响图4-14图4-14不同流道宽度情况下的第一法向应力图中由上到下分别为流道宽度为15μm、20μm、30μm的τ11，τ12，τ22的等值线图从图4-14～图4-16中可以看出：随着流道宽度的增加，流场中τ11，τ12，τ22的极值都增大。即油膜宽度越大，流场中的应力极值就越大。图4-15图4-15不同流道宽度情况下的第二法向应力图4-16图4-16不同流道宽度情况下的切向应力

4.3.3油膜高度对于应力的影响图4-17图4-17不同油膜高度情况下的第一法向应力图4-18图4-18不同油膜高度情况下的第二法向应力图中由上到下分别为油膜高度为3μm、4μm、5μm的τ11，τ12，τ22的等值线图。从图4-17～图4-19中可以看出：随着油膜高度的增大，流场中τ11，τ12，τ22的极值都增大。即油膜突起高度越大，流场中的应力极值就越大。图4-19图4-19不同油膜高度情况下的切向应力4.3.4界面张力对于应力的影响图4-20图4-20不同界面张力情况下的第一法向应力图中由上到下分别为界面张力为0.001mN/m、0.01mN/m、0.1mN/m的τ11，τ12，τ22的等值线图。从图4-20～图4-22中可以看出：随着界面张力的增大，流场中τ11，τ12，τ22的极值都增大。即原油与水之间的界面张力越大，流场中的应力极值就越大。图4-21图4-21不同界面张力情况下的第二法向应力图4-22图4-22不同界面张力情况下的切向应力4.3.5油膜粘度对于应力的影响图中由上到下分别为油膜粘度为0.015Pa·s、0.02Pa·s、0.03Pa·s的τ11，τ12，τ22的等值线图。从图4-23～图4-25中可以看出：随着原油粘度的增大，流场中τ11，τ12，τ22的极值都增大。即原油粘度越大，流场中的应力极值就越大。图4-23图4-23不同原油粘度情况下的第一法向应力图4-24图4-24不同原油粘度情况下的第二法向应力图4-25图4-25不同原油粘度情况下的切向应力

4.4油膜变形模拟图4-26~图4-30分别给出了不同压力梯度下、不同流道宽度、不同油膜高度、不同界面张力、不同油膜粘度的变形情况4.4.1压力梯度对于油膜变形的影响图4-26图4-26不同压力梯度情况下的油膜变形图中由上到下分别为压力梯度为0.02MPa、0.002MPa、0.0002MPa的油膜变形图4.4.2流道宽度对于变形的影响图4-27图4-27不同流道宽度情况下的油膜变形图中由上到下分别为流道宽度为15μm、20μm、30μm的油膜变形图。

4.4.3油膜高度对于变形的影响图4-28图4-28不同油膜高度情况下的油膜变形图中由上到下分别为油膜高度为3μm、4μm、5μm的油膜变形图。4.4.4界面张力对于油膜变形的影响图中由上到下分别为界面张力为0.001mN/m、0.01mN/m、0.1mN/m的油膜变形图。图4-29图4-29不同界面张力情况下的油膜变形

4.4.5原油粘度对于油膜变形的影响图4-30图4-30不同原油粘度情况下的油膜变形图中由上到下分别为油膜粘度为0.015Pa·s、0.02Pa·s、0.03Pa·s的油膜变形图。由图4-26~图4-30可知，压力梯度、油膜高度以及油水之间的界面张力对于油膜的变形起着决定性的作用，在其他条件相同的条件下，压力梯度越大油滴，变形也越大，油膜高度越大，变形越大，界面张力越小，变形越大。结论本文采用牛顿流体本构模型来描述驱替液的流变性，建立驱替液在微观孔隙空间中流动的物理模型和数学模型，使用POLYFLOW软件进行数值分析，在不同压力梯度，不同流道宽度和油膜高度，不同界面张力，不同原油粘度情况下，求解出流动的速度场、应力场以及应变场得到以下结论：(1)油膜顶部的流动速度与压力梯度、流道宽度和油膜高度、界面张力、原油粘度有关。压力梯度越大，流道宽度越小，油膜高度越大，界面张力越小，原油粘度越小时，油膜顶部的流动速度越大；(2)作用在残余油上的法向应力、切向应力与压力梯度、界面张力、油膜高度和流道宽度、原油粘度有关，压力梯度越大，油膜高度越大，流道宽度越宽，原油粘度越大时，法向应力与切向应力越大。为进一步分析小油滴剥离油膜时的条件以及受力奠定基础；(3)压力梯度、油膜高度以及油水之间的界面张力对于油膜的变形起着决定性的作用，在其他条件相同的条件下，压力梯度越大油滴，变形也越大，油膜高度越大，变形越大，界面张力越小，变形越大。由此可以推断，压力梯度越大，油膜高度越大，油水之间界面张力越小时，驱油效率越高。

参考文献[1]隋新光.大庆油田低渗透储层动用条件及渗流特征研究[D].北京:中国石油大学,2000.[2]朱九成.岩石孔隙结构分形模型及随机网络仿真[J].石油学报,1995,19(6):46-51.[3]李玉彬.微焦点X-CT成象技术研究岩石微观特征[J].油气采收率,2000,7(4):50-52.[4]汪周华,吐依洪江,郭平,等.凝析气藏水驱机理研究[J].西南石油,2006,28(6):36-39.[5]王域辉,廖淑华.分形与石油[M].北京:石油工业出版社,1994:33-142.[6]张济忠.分形[M].北京:清华大学出版社,1995:1-66.[7]ZHANGLi-qiang,JIYou-liang,MAWen-jie,etal.Characteristicsoffra