文聚驱井采出液流动规律：多维度解析与工程应用

文聚驱井采出液流动规律：多维度解析与工程应用一、引言1.1研究背景与意义在全球能源需求持续增长的大背景下，石油作为最重要的能源之一，其稳定供应对国家经济发展和社会稳定至关重要。随着石油勘探开发的不断深入，传统采油工艺在面对复杂油藏条件时，愈发难以满足日益增长的能源需求，提高采收率成为石油工业发展的关键任务。聚合物驱油技术作为一种重要的提高采收率方法，在注水开发油田中得到了广泛的推广与应用。其原理是向油层注入聚合物，利用聚合物增加水相粘度，并因聚合物的滞留引起油层渗透率下降。这一技术有效控制了水淹层段中水相流度，改善了水油流度比，提高了水淹层段的实际驱油效率；同时降低了高渗透率的水淹层段中流体总流度，缩小了高低层段间水线推进速度差，调整了吸水剖面，进而提高了实际波及系数。以国内的胜利油田为例，采用聚合物驱油技术后，采收率从20%大幅提高到了70%以上，生产效益显著提升。然而，聚合物驱采油井在见聚后，采出液性质发生显著变化。随着采出液中聚合物溶液含量的不同，采出液从牛顿流体转变为非牛顿流体，甚至呈现出弹性效应。这种性质的改变导致聚驱井产生了一系列较水驱井更为严重的生产问题。其中，杆管偏磨问题尤为突出，不仅增加了设备维护成本和更换频率，还影响了油井的正常生产时效，降低了原油产量。研究聚驱井采出液的流动规律，对于深入理解聚驱井生产过程中的物理现象、解决生产问题以及进一步提高采收率具有重要意义。从解决生产问题角度来看，通过掌握采出液流动规律，可以精准分析杆管偏磨等问题产生的机理。例如，有杆抽油系统在下冲程中，泵端集中轴向压力和轴向分布力的综合作用使抽油杆受压，加之偏心环空内聚驱采出液因弹性对抽油杆产生侧向力，这两个因素是导致聚驱井杆管偏磨的主要原因。明确这些机理后，就能够针对性地提出有效的解决措施，如优化抽油杆结构、调整采出液性质等，从而减少杆管偏磨，降低设备损耗，提高油井生产的稳定性和经济性。从提高采收率角度而言，深入了解采出液在井筒和油藏中的流动规律，有助于优化聚合物驱油方案。可以根据不同油藏条件和采出液特性，合理调整聚合物的注入量、注入方式以及驱油工艺参数，提高聚合物在油藏中的波及效率和驱油效果，进一步挖掘油藏潜力，提高原油采收率，满足全球对石油能源不断增长的需求，保障国家能源安全。1.2研究现状与进展聚合物驱油技术作为提高采收率的重要手段，自20世纪50年代提出后，便受到了国内外学者的广泛关注。国外在聚合物驱油技术的基础研究和应用方面起步较早，取得了一系列重要成果。美国早在1964年就开始了聚合物驱油的现场试验，随后在得克萨斯州、俄克拉何马州等多个油田进行了大规模应用，通过优化聚合物的注入参数和驱油工艺，有效提高了原油采收率。前苏联也在聚合物驱油技术方面进行了大量研究和实践，在西西伯利亚等油田成功应用该技术，实现了原油产量的稳定增长。国内对聚合物驱油技术的研究始于20世纪70年代，经过多年的技术攻关和现场试验，在理论研究和实际应用方面都取得了显著进展。大庆油田作为国内聚合物驱油技术应用的典范，从1996年开始大规模推广聚合物驱油技术，目前已成为油田提高采收率的主导技术之一。通过不断优化聚合物的配方和注入工艺，大庆油田聚合物驱油的平均采收率提高了10%-15%，为保障国家能源安全做出了重要贡献。胜利油田、辽河油田等也在聚合物驱油技术方面进行了大量实践，取得了良好的应用效果。在聚驱井采出液流动规律的研究方面，国内外学者主要从聚合物溶液的流变性、环空内流体的流动以及机采井的杆管偏磨问题等几个关键方向展开深入探索。聚合物溶液的流变性是研究聚驱井采出液流动规律的基础。国外学者如Bird等在粘弹性流体的本构方程研究方面取得了重要成果，提出了Maxwell模型、Oldroyd-B模型等多种描述粘弹性流体特性的本构方程，为聚合物溶液流变性的研究提供了重要的理论基础。国内学者在聚合物溶液流变性研究方面也取得了丰硕成果。中国石油大学的黄善波通过实验研究了聚合物溶液在不同剪切速率、温度和浓度条件下的粘度变化规律，发现聚合物溶液的粘度随着剪切速率的增加而降低，呈现出典型的剪切变稀特性；同时，温度升高会导致聚合物溶液的粘度下降，而浓度增加则会使粘度增大。大庆油田的常瑞清应用德国HAAK公司RS100流变仪及数显密度计，对现场抽油机井产出液以及配制聚合物溶液开展了室内实验研究，测量了油井产出液和聚合物溶液在不同剪切速率下粘度随时间的变化、产出液粘度随含水的变化、不同含聚浓度井产出液的粘度变化、配制聚合物溶液密度变化情况，实验结果为抽油机井的优化设计提供了依据。环空内流体的流动是聚驱井采出液流动规律研究的重要内容。国外学者如Hanks等采用数值模拟方法研究了同心环空中幂律流体的非定常流动，分析了环空间距、流体性质等因素对流动的影响。国内学者在这方面也进行了深入研究。黄善波采用固有函数法从理论上分别对同心环空中牛顿流体、Maxwell型粘弹性流体在周期性压力梯度作用下的振荡流、内壁面沿轴向做往复运动的环空内的非定常流动进行了研究；采用SIMPLE数值方法分析了同心环空内幂律流体因入口速度脉动产生的振荡流和内壁面沿轴向做往复运动的环空内的非定常流，得到了各种流体在两种情形下环空内非定常流动的速度、平均速度、壁面处应力或摩擦系数等的变化规律，分析和讨论了环空间距、流体性质、振荡频率和振幅等因素对流动的影响。机采井的杆管偏磨问题是聚驱井生产中面临的一个严重问题，国内外学者对此进行了大量研究。国外学者通过实验和数值模拟等方法，分析了杆管偏磨的原因和影响因素，提出了一些防治措施，如优化抽油杆结构、采用耐磨材料等。国内学者在杆管偏磨问题的研究方面也取得了重要成果。黄善波分析表明，有杆抽油系统在下冲程中泵端集中轴向压力和轴向分布力的综合作用使抽油杆处于受压状态而有可能失稳弯曲，以及偏心环空内聚驱采出液因弹性而对抽油杆产生的侧向力是聚驱井发生杆管偏磨的两个主要原因，其中抽油杆在因弹性而产生的侧向力作用下的偏磨是聚驱井所特有的。采用Maxwell模型计算了采出液在偏心环空内流动时的侧向力，指出聚合物溶液的增稠特性使聚驱井的中和点上移、弯曲载荷增大、临界弯曲载荷减小，随抽油杆的增长，侧向力很容易超过杆管发生偏磨时的载荷。尽管国内外在聚驱井采出液流动规律研究方面取得了一定成果，但仍存在一些不足之处。现有研究大多集中在单一因素对采出液流动规律的影响，而实际生产中采出液的流动受到多种因素的综合作用，对多因素耦合作用下采出液流动规律的研究还相对较少。在采出液流变性的研究中，对于复杂油藏条件下聚合物溶液的粘弹性特性及其对流动的影响，还需要进一步深入研究。此外，目前的研究成果在实际生产中的应用还存在一定的局限性，如何将理论研究成果更好地转化为实际生产中的有效措施，仍有待进一步探索。1.3研究内容与方法1.3.1研究内容聚驱井采出液的制备及性能测试：根据实际油藏条件和聚合物驱油工艺，在实验室中制备不同聚合物浓度、不同油水比例的聚驱井采出液模拟样品。运用先进的实验仪器，如旋转流变仪、密度计等，精确测量采出液的粘度、密度、粘弹性等关键性能参数，深入分析聚合物浓度、温度、剪切速率等因素对采出液性能的影响规律。聚驱井采出液在管道输送过程中的流动规律研究：构建聚驱井采出液在管道中流动的物理模型和数学模型，综合运用理论分析、数值模拟和实验研究等方法，深入探究采出液在不同管径、流速、温度等条件下的流动特性。重点研究采出液的流型转变、压力损失、速度分布等关键流动参数的变化规律，分析聚合物的粘弹性对流动稳定性的影响，揭示聚驱井采出液在管道输送过程中的复杂流动机制。聚驱井采出液在岩石孔隙中的渗透力学分析：利用微观可视化实验技术和数值模拟方法，研究聚驱井采出液在岩石孔隙中的渗流特性。分析岩石孔隙结构、渗透率、润湿性等因素对采出液渗流的影响，建立考虑聚合物吸附、滞留和岩石孔隙变形的渗流力学模型，深入探讨采出液在岩石孔隙中的微观渗流机理，为提高油藏采收率提供理论依据。聚驱井采出液在不同油层条件下的应用研究：结合实际油藏地质数据，对不同油层条件下的聚驱井采出液进行数值模拟和现场试验研究。分析油层厚度、渗透率分布、原油粘度等因素对聚驱井采出液流动和驱油效果的影响，优化聚合物驱油方案，提出适合不同油层条件的聚驱井采出液应用工艺和技术措施，提高聚驱井的生产效率和原油采收率。1.3.2研究方法实验室模拟测试：在实验室环境中，严格模拟聚驱井的实际工况，制备各种不同参数的采出液样品。利用高精度的实验仪器，如流变仪、密度计、界面张力仪等，对采出液的粘度、密度、粘弹性、界面张力等物理性质进行精确测量。通过改变实验条件，如温度、压力、剪切速率等，系统研究各因素对采出液性能的影响规律，为后续的理论分析和数值模拟提供可靠的实验数据支持。数值模拟方法：采用计算流体力学（CFD）等先进的数值模拟技术，建立聚驱井采出液在井筒、管道和油藏中的流动模型。运用有限元法、有限体积法等数值计算方法，对模型进行求解，模拟采出液在不同条件下的流动过程，得到采出液的速度场、压力场、浓度场等详细信息。通过与实验室模拟结果进行对比验证，不断优化数值模型，提高模拟结果的准确性和可靠性，为聚驱井采出液流动规律的研究提供高效、便捷的研究手段。野外实验法：选择具有代表性的聚驱井现场，开展野外实验研究。在实际生产条件下，对聚驱井采出液的流动参数进行实时监测和数据采集，如井口压力、流量、含水率等。通过对现场数据的分析，深入了解采出液在实际油藏环境中的流动特性和变化规律，验证实验室模拟和数值模拟的结果，为聚驱井采出液流动规律的研究提供真实、可靠的现场数据支持，同时为聚驱井采出液的实际应用提供技术指导。二、文聚驱井采出液特性分析2.1采出液的组成与性质文聚驱井采出液是一种成分复杂的多相混合物，主要由聚合物、原油、水以及少量的固体颗粒等组成。其中，聚合物作为提高采收率的关键添加剂，其在采出液中的含量和性质对采出液的整体特性起着至关重要的作用。在文聚驱油过程中，常用的聚合物为部分水解聚丙烯酰胺（HPAM），它具有长链高分子结构，在水中能够通过分子链的伸展和缠绕形成网状结构，从而增加水相的粘度。原油在采出液中以油滴的形式存在，其含量和性质因油藏条件的不同而有所差异。原油的粘度、密度、组成成分等因素会影响采出液的流动性和油水分离特性。一般来说，高粘度的原油会使采出液的整体粘度增加，导致流动阻力增大；而原油中的轻质组分和重质组分的比例也会影响采出液的稳定性和界面性质。水是采出液的主要连续相，其矿化度、pH值等性质对采出液的性质也有重要影响。高矿化度的水会使聚合物分子链发生卷曲，降低聚合物的增粘效果；而pH值的变化则可能影响聚合物的水解程度和稳定性，进而改变采出液的流变性。与传统水驱采出液相比，文聚驱井采出液最显著的特点是其非牛顿流体特性。牛顿流体在流动过程中，剪切应力与剪切速率呈线性关系，其粘度不随剪切速率的变化而改变，遵循牛顿粘性定律，数学表达式为\tau=\mu\dot{\gamma}，其中\tau为剪切应力，\mu为动力粘度，\dot{\gamma}为剪切速率。然而，文聚驱井采出液中的聚合物分子链在流动过程中会发生取向和变形，导致其剪切应力与剪切速率之间呈现非线性关系，粘度随剪切速率的变化而变化，不满足牛顿粘性定律，属于非牛顿流体。通过旋转流变仪对文聚驱井采出液进行流变测试，得到其流变曲线（如图1所示）。从图中可以看出，在低剪切速率范围内，采出液的粘度较高，且随剪切速率的增加变化较小；当剪切速率增大到一定程度后，粘度迅速下降，呈现出明显的剪切变稀特性。这种剪切变稀特性是由于聚合物分子链在高剪切速率下发生解缠结和取向，导致分子间的相互作用减弱，从而使粘度降低。【此处插入图1：文聚驱井采出液流变曲线】文聚驱井采出液还表现出明显的弹性效应。当采出液受到剪切作用时，不仅会产生粘性流动，还会储存一部分弹性变形能，在剪切作用停止后，这部分弹性变形能会释放出来，使流体发生弹性回复。这种弹性效应在许多流动现象中都有体现，如爬杆效应（韦森堡效应）和射流胀大现象。在爬杆效应实验中，将一根旋转的圆柱放入文聚驱井采出液中，采出液会沿圆柱向上爬升，而牛顿流体则会在离心力作用下沿圆柱向下流动。这是因为采出液的弹性使其在旋转圆柱表面产生了法向应力，从而推动流体向上爬升。在射流胀大实验中，当采出液从毛细管中射出时，射流的直径会大于毛细管的直径，这是由于采出液在毛细管内流动时储存的弹性变形能在射出毛细管后释放，导致射流发生膨胀。文聚驱井采出液的弹性效应可以用第一法向应力差N_1=\sigma_{11}-\sigma_{22}和第二法向应力差N_2=\sigma_{22}-\sigma_{33}来描述，其中\sigma_{ij}为应力张量的分量。第一法向应力差通常为正值，且远大于第二法向应力差，它是导致爬杆效应和射流胀大现象的主要原因。通过实验测量采出液在不同剪切速率下的第一法向应力差和第二法向应力差，可以定量地研究其弹性效应。采出液的弹性效应与其聚合物含量、分子量以及溶液的浓度等因素密切相关。一般来说，聚合物含量越高、分子量越大，采出液的弹性效应越明显。此外，温度、剪切历史等因素也会对采出液的弹性产生影响。温度升高会使聚合物分子链的热运动加剧，分子间的相互作用减弱，从而降低采出液的弹性；而剪切历史则会使聚合物分子链发生降解和取向，改变其微观结构，进而影响采出液的弹性性能。2.2采出液的流变性研究流变性是指流体在外力作用下发生变形和流动的特性，对于文聚驱井采出液而言，其流变性研究至关重要，因为它直接影响着采出液在井筒、管道以及油藏孔隙中的流动行为。为了深入探究文聚驱井采出液的流变性，本研究采用了先进的旋转流变仪进行实验，并结合理论分析方法，系统地研究了采出液的粘度、非牛顿指数、稠度系数等流变参数随浓度、温度、剪切速率等因素的变化规律。实验选用安东帕MCR302型旋转流变仪，该仪器具有高精度、宽量程的特点，能够精确测量不同条件下采出液的流变参数。实验过程中，首先制备一系列不同聚合物浓度的采出液样品，聚合物浓度范围设定为500mg/L-2500mg/L，以模拟实际生产中可能遇到的不同情况。同时，控制温度在25℃-65℃之间变化，以研究温度对采出液流变性的影响。在测量粘度时，采用同心圆筒测量系统，将采出液样品注入外筒与内筒之间的环形间隙中。通过设定不同的剪切速率，从0.1s⁻¹逐渐增加到1000s⁻¹，测量采出液在不同剪切速率下的剪切应力，根据粘度的定义\eta=\frac{\tau}{\dot{\gamma}}（其中\eta为粘度，\tau为剪切应力，\dot{\gamma}为剪切速率），计算得到采出液的粘度。实验结果表明，文聚驱井采出液的粘度随聚合物浓度的增加而显著增大（如图2所示）。当聚合物浓度从500mg/L增加到2500mg/L时，采出液在低剪切速率下（如0.1s⁻¹）的粘度从10mPa・s左右增大到100mPa・s以上。这是因为随着聚合物浓度的升高，聚合物分子链之间的相互作用增强，形成了更为致密的网状结构，阻碍了流体的流动，从而导致粘度增大。【此处插入图2：不同聚合物浓度下采出液粘度随剪切速率的变化曲线】采出液的粘度还呈现出明显的剪切变稀特性，即随着剪切速率的增加，粘度逐渐降低。在低剪切速率范围内，聚合物分子链相互缠绕，形成较为稳定的结构，此时粘度较大；当剪切速率增大时，分子链在剪切力的作用下逐渐解缠结并取向，分子间的相互作用减弱，粘度随之降低。例如，对于聚合物浓度为1500mg/L的采出液，在剪切速率为0.1s⁻¹时，粘度约为50mPa・s，而当剪切速率增加到1000s⁻¹时，粘度降至10mPa・s以下。温度对采出液粘度的影响也十分显著（如图3所示）。随着温度的升高，采出液的粘度逐渐降低。当温度从25℃升高到65℃时，聚合物浓度为1000mg/L的采出液在剪切速率为1s⁻¹时的粘度从30mPa・s左右降至15mPa・s左右。这是由于温度升高使聚合物分子链的热运动加剧，分子间的相互作用力减弱，分子链的柔韧性增加，更容易发生解缠结和取向，从而导致粘度降低。【此处插入图3：不同温度下采出液粘度随剪切速率的变化曲线】非牛顿指数n和稠度系数K是描述非牛顿流体流变特性的重要参数，它们可以通过幂律模型\tau=K\dot{\gamma}^n（其中\tau为剪切应力，\dot{\gamma}为剪切速率，K为稠度系数，n为非牛顿指数）来确定。通过对实验数据进行拟合，可以得到不同条件下采出液的非牛顿指数和稠度系数。实验结果表明，文聚驱井采出液的非牛顿指数n均小于1，表明其属于假塑性流体，即具有剪切变稀特性。随着聚合物浓度的增加，非牛顿指数n逐渐减小，这意味着采出液的剪切变稀特性更加明显。例如，当聚合物浓度为500mg/L时，非牛顿指数n约为0.8，而当聚合物浓度增加到2500mg/L时，n减小至0.6左右。这是因为高浓度的聚合物分子链之间的相互作用更强，在剪切作用下分子链的解缠结和取向更加容易，导致剪切变稀特性增强。稠度系数K则随着聚合物浓度的增加而增大，反映了聚合物浓度对采出液粘性的影响。当聚合物浓度从500mg/L增加到2500mg/L时，稠度系数K从0.5Pa・sⁿ左右增大到2.0Pa・sⁿ以上。这表明随着聚合物浓度的升高，采出液的粘性增大，流动阻力增加。温度对非牛顿指数n和稠度系数K也有一定的影响。随着温度的升高，非牛顿指数n略有增大，说明温度升高使采出液的剪切变稀特性有所减弱。这是因为温度升高使聚合物分子链的热运动增强，分子链的柔韧性增加，在一定程度上抑制了分子链在剪切作用下的解缠结和取向，从而使剪切变稀特性减弱。而稠度系数K则随着温度的升高而减小，表明温度升高降低了采出液的粘性。这是由于温度升高导致分子间的相互作用力减弱，使采出液的粘性降低。三、聚驱井采出液在同心环空内的非定常流分析3.1牛顿流体同心环空内的非定常流3.1.1振荡流分析为了深入研究牛顿流体在同心环空内的振荡流特性，构建相应的物理及数学模型。考虑一个由半径为R_1的内管和半径为R_2的外管组成的同心环形通道，其中充满牛顿流体。假设流体不可压缩，且流动为层流状态。在通道入口处施加一个周期性变化的压力梯度G(t)=G_0\sin(\omegat)，其中G_0为压力梯度的幅值，\omega为振荡频率，t为时间。基于上述物理模型，根据牛顿流体的本构方程\tau=\mu\frac{du}{dr}（其中\tau为剪切应力，\mu为动力粘度，u为速度，r为径向坐标）以及连续性方程\frac{\partialu}{\partialr}+\frac{u}{r}=0和动量方程\rho\frac{\partialu}{\partialt}=-\frac{\partialp}{\partialz}+\mu(\frac{\partial^2u}{\partialr^2}+\frac{1}{r}\frac{\partialu}{\partialr})（其中\rho为流体密度，p为压力，z为轴向坐标），建立数学模型。将压力梯度G(t)代入动量方程中，并结合边界条件：在r=R_1处，u=0；在r=R_2处，u=0，得到如下定解问题：\begin{cases}\rho\frac{\partialu}{\partialt}=G_0\sin(\omegat)+\mu(\frac{\partial^2u}{\partialr^2}+\frac{1}{r}\frac{\partialu}{\partialr})&(R_1<r<R_2,t>0)\\u(R_1,t)=0&(t>0)\\u(R_2,t)=0&(t>0)\\u(r,0)=0&(R_1<r<R_2)\end{cases}为求解上述数学模型，采用分离变量法。设u(r,t)=R(r)T(t)，代入方程中可得：\frac{\rho}{T}\frac{dT}{dt}=G_0\sin(\omegat)+\frac{\mu}{R}(\frac{d^2R}{dr^2}+\frac{1}{r}\frac{dR}{dr})等式左边仅与时间t有关，右边仅与径向坐标r有关，因此两边应等于一个常数，设为-\lambda^2。由此得到两个常微分方程：\begin{cases}\frac{dT}{dt}+\frac{\lambda^2T}{\rho}t=\frac{G_0}{\rho}\sin(\omegat)\\\frac{d^2R}{dr^2}+\frac{1}{r}\frac{dR}{dr}+\frac{\lambda^2}{\mu}R=0\end{cases}对于第一个常微分方程，利用一阶线性非齐次常微分方程的求解方法，其通解为：T(t)=e^{-\frac{\lambda^2t}{\rho}}\left(\int\frac{G_0}{\rho}\sin(\omegat)e^{\frac{\lambda^2t}{\rho}}dt+C_1\right)通过积分运算可得：T(t)=e^{-\frac{\lambda^2t}{\rho}}\left(\frac{G_0}{\rho}\frac{\omegae^{\frac{\lambda^2t}{\rho}}}{\lambda^4+\omega^2\rho^2}(\lambda^2\sin(\omegat)-\omega\rho\cos(\omegat))+C_1\right)对于第二个常微分方程，它是零阶贝塞尔方程，其通解为：R(r)=C_2J_0(\frac{\lambdar}{\sqrt{\mu}})+C_3Y_0(\frac{\lambdar}{\sqrt{\mu}})其中J_0和Y_0分别为零阶第一类和第二类贝塞尔函数。根据边界条件u(R_1,t)=0和u(R_2,t)=0，可得：\begin{cases}C_2J_0(\frac{\lambdaR_1}{\sqrt{\mu}})+C_3Y_0(\frac{\lambdaR_1}{\sqrt{\mu}})=0\\C_2J_0(\frac{\lambdaR_2}{\sqrt{\mu}})+C_3Y_0(\frac{\lambdaR_2}{\sqrt{\mu}})=0\end{cases}为使C_2和C_3不全为零，系数行列式必须为零，即：\begin{vmatrix}J_0(\frac{\lambdaR_1}{\sqrt{\mu}})&Y_0(\frac{\lambdaR_1}{\sqrt{\mu}})\\J_0(\frac{\lambdaR_2}{\sqrt{\mu}})&Y_0(\frac{\lambdaR_2}{\sqrt{\mu}})\end{vmatrix}=0由此可确定\lambda的值，进而得到速度分布u(r,t)的表达式。得到速度分布后，可进一步计算其他参数。平均速度\overline{u}可通过对速度在环空横截面上积分并除以环空横截面积得到：\overline{u}=\frac{1}{\pi(R_2^2-R_1^2)}\int_{R_1}^{R_2}2\piru(r,t)dr壁面处应力\tau_w可根据牛顿流体本构方程在壁面处的取值计算：\tau_w|_{r=R_1}=\mu\left(\frac{\partialu}{\partialr}\right)_{r=R_1}\tau_w|_{r=R_2}=\mu\left(\frac{\partialu}{\partialr}\right)_{r=R_2}通过数值计算，分析不同参数对振荡流的影响。当环空间距\DeltaR=R_2-R_1增大时，速度分布的变化较为明显（如图4所示）。在相同的压力梯度幅值和振荡频率下，环空间距增大，流体的流动阻力减小，速度幅值增大。这是因为环空间距增大，流体的流通面积增大，在相同流量下，流速相对减小，从而导致速度梯度减小，粘性力对流体的阻碍作用减弱，使得速度幅值增大。【此处插入图4：不同环空间距下速度随时间的变化曲线】振荡频率\omega对速度分布也有显著影响（如图5所示）。随着振荡频率的增加，速度幅值减小，且速度变化的周期变短。这是因为振荡频率增加，流体在单位时间内受到的压力变化次数增多，流体来不及充分响应压力的变化，导致速度幅值减小。同时，由于压力变化周期变短，速度变化的周期也相应变短。【此处插入图5：不同振荡频率下速度随时间的变化曲线】流体性质（如动力粘度\mu）的改变同样会影响振荡流特性（如图6所示）。当动力粘度增大时，速度幅值明显减小。这是因为动力粘度增大，流体的粘性增强，内部摩擦力增大，阻碍了流体的流动，使得速度幅值减小。【此处插入图6：不同动力粘度下速度随时间的变化曲线】3.1.2内壁面做往复运动的环空内非定常流建立物理与数学模型来研究牛顿流体在内壁面做往复运动的环空内的非定常流。物理模型仍为同心环形通道，内管半径为R_1，外管半径为R_2，充满牛顿流体。不同的是，内管内壁面以速度u_w(t)=u_0\sin(\omegat)沿轴向做往复运动，其中u_0为运动幅值，\omega为运动频率。基于此物理模型，建立数学模型。连续性方程和动量方程与振荡流分析中的相同，边界条件变为：在r=R_1处，u=u_0\sin(\omegat)；在r=R_2处，u=0。由此得到定解问题：\begin{cases}\rho\frac{\partialu}{\partialt}=-\frac{\partialp}{\partialz}+\mu(\frac{\partial^2u}{\partialr^2}+\frac{1}{r}\frac{\partialu}{\partialr})&(R_1<r<R_2,t>0)\\u(R_1,t)=u_0\sin(\omegat)&(t>0)\\u(R_2,t)=0&(t>0)\\u(r,0)=0&(R_1<r<R_2)\end{cases}同样采用分离变量法求解该数学模型。设u(r,t)=R(r)T(t)，代入方程得到两个常微分方程：\begin{cases}\frac{dT}{dt}+\frac{\lambda^2T}{\rho}t=0\\\frac{d^2R}{dr^2}+\frac{1}{r}\frac{dR}{dr}+\frac{\lambda^2}{\mu}R=0\end{cases}对于第一个常微分方程，其通解为：T(t)=C_4e^{-\frac{\lambda^2t}{\rho}}对于第二个常微分方程，通解为：R(r)=C_5J_0(\frac{\lambdar}{\sqrt{\mu}})+C_6Y_0(\frac{\lambdar}{\sqrt{\mu}})根据边界条件u(R_1,t)=u_0\sin(\omegat)和u(R_2,t)=0，确定C_4、C_5和C_6的值，从而得到速度分布u(r,t)的表达式。对计算结果进行分析，讨论不同因素对流动的影响。当内壁面运动幅值u_0增大时，环空内流体的速度幅值也随之增大（如图7所示）。这是因为内壁面运动幅值增大，直接带动与内壁面接触的流体层运动速度增大，通过粘性力的传递，使得整个环空内流体的速度幅值都增大。【此处插入图7：不同内壁面运动幅值下速度随时间的变化曲线】运动频率\omega的变化对速度分布有重要影响（如图8所示）。随着运动频率增加，速度幅值减小，且速度变化更加剧烈。这是因为运动频率增加，流体在短时间内受到内壁面的加速和减速作用更加频繁，使得流体难以保持较高的速度，速度幅值减小。同时，由于频繁的速度变化，速度分布的波动更加明显。【此处插入图8：不同运动频率下速度随时间的变化曲线】环空间距\DeltaR对流动也有一定影响（如图9所示）。环空间距增大，速度幅值略有增大，但变化相对较小。这是因为环空间距增大，虽然流体的流通面积增大，流动阻力减小，但由于内壁面运动对流体的带动作用主要集中在内壁面附近，随着环空间距增大，这种带动作用的影响范围相对减小，因此速度幅值的增大并不明显。【此处插入图9：不同环空间距下速度随时间的变化曲线】3.2幂律流体同心环空内的非定常流3.2.1非定常流模型建立构建幂律流体在同心环空内非定常流的物理模型，考虑一个由半径为R_1的内管和半径为R_2的外管组成的同心环形通道，其中充满幂律流体。假设流体不可压缩，且流动为层流状态。幂律流体的本构方程为\tau=K(\dot{\gamma})^{n-1}\dot{\gamma}，其中\tau为剪切应力张量，K为稠度系数，n为幂律指数，\dot{\gamma}为剪切速率张量，\dot{\gamma}=\frac{1}{2}(\nabla\vec{u}+(\nabla\vec{u})^T)，\vec{u}为速度矢量。守恒型的控制方程包括连续性方程\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{u})=0，其中\rho为流体密度，以及动量方程\frac{\partial(\rho\vec{u})}{\partialt}+\nabla\cdot(\rho\vec{u}\vec{u})=-\nablap+\nabla\cdot\boldsymbol{\tau}，其中p为压力。定解条件如下：在初始时刻t=0时，速度\vec{u}(r,0)=0；在r=R_1处，若内管做轴向往复运动，速度u(R_1,t)=u_0\sin(\omegat)，若研究进口流速正弦变化的振荡流，此处速度条件根据具体入口条件确定；在r=R_2处，速度u(R_2,t)=0。3.2.2数值计算与结果分析采用有限体积法对控制方程进行离散化，将计算区域划分为一系列的控制体积，在每个控制体积上对控制方程进行积分，将偏微分方程转化为代数方程。离散后的方程通过迭代求解，常用的迭代算法有SIMPLE算法（Semi-ImplicitMethodforPressure-LinkedEquations）及其改进算法，如SIMPLER算法（SIMPLERevised）、SIMPLEC算法（SIMPLEConsistent）等。以具体算例进行计算，设定同心环空内管半径R_1=0.01m，外管半径R_2=0.02m，幂律流体的稠度系数K=0.5Pa\cdots^n，幂律指数n=0.8。当进口流速以正弦规律变化时，设进口流速u_{in}(t)=u_{max}\sin(\omegat)，其中u_{max}=0.5m/s，\omega=2\pirad/s。得到速度分布、平均速度和壁面处应力随时间的变化规律（如图10-12所示）。从速度分布（图10）可以看出，在一个周期内，速度沿径向呈现出非均匀分布，靠近内管处速度变化较为剧烈，这是因为进口流速的变化对内管附近流体的影响更为直接。平均速度（图11）也随时间做周期性变化，其幅值与进口流速的幅值和频率密切相关。壁面处应力（图12）同样呈现周期性变化，且内管壁面处的应力幅值大于外管壁面处，这是由于内管附近的速度梯度较大。【此处插入图10：进口流速正弦变化时速度随时间和径向位置的变化】【此处插入图11：进口流速正弦变化时平均速度随时间的变化】【此处插入图12：进口流速正弦变化时壁面处应力随时间的变化】当内壁面做轴向往复运动时，设内管内壁面速度u_w(t)=u_0\sin(\omegat)，其中u_0=0.3m/s，\omega=3\pirad/s。分析速度分布、平均速度和壁面处应力的变化情况（如图13-15所示）。速度分布（图13）显示，在内壁面运动的带动下，靠近内管的流体速度明显增大，且速度分布随时间发生周期性变化。平均速度（图14）也随时间做周期性波动，其波动范围与内壁面运动的幅值和频率有关。壁面处应力（图15）在一个周期内呈现出明显的变化，内管壁面处的应力变化更为显著，这是由于内壁面的运动直接作用于内管附近的流体，导致此处的剪切应力变化较大。【此处插入图13：内壁面轴向往复运动时速度随时间和径向位置的变化】【此处插入图14：内壁面轴向往复运动时平均速度随时间的变化】【此处插入图15：内壁面轴向往复运动时壁面处应力随时间的变化】通过对不同条件下的计算结果分析可知，环空间距、流体性质（稠度系数K和幂律指数n）、振荡频率\omega和振幅等因素对幂律流体在同心环空内的非定常流有显著影响。环空间距增大，流体的流通面积增大，速度分布会变得相对均匀，平均速度增大，壁面处应力减小；稠度系数增大，流体的粘性增强，流动阻力增大，速度幅值减小，壁面处应力增大；幂律指数n的变化会改变流体的非牛顿特性，当n减小时，流体的剪切变稀特性更加明显，速度分布和壁面处应力的变化也会相应改变；振荡频率增大，流体在单位时间内受到的扰动次数增多，速度幅值减小，平均速度的波动频率增加；振幅增大，速度和壁面处应力的幅值也随之增大。3.3粘弹性流体同心环空内的非定常流3.3.1基本方程与起动流分析建立粘弹性流体在同心环空内非定常流的基本方程，考虑一个由半径为R_1的内管和半径为R_2的外管组成的同心环形通道，其中充满粘弹性流体。假设流体不可压缩，且流动为层流状态。粘弹性流体的本构方程采用上随体Maxwell模型，其表达式为\tau_{ij}+\lambda\frac{\delta\tau_{ij}}{\deltat}=\mu_0A_{ij}，其中\tau_{ij}为偏应力张量，\lambda为松弛时间，\mu_0为零剪切黏度，A_{ij}为一阶Rivlin-Erickson张量，\frac{\delta\tau_{ij}}{\deltat}为上随体导数，其表达式为\frac{\delta\tau_{ij}}{\deltat}=\frac{\partial\tau_{ij}}{\partialt}+u_k\frac{\partial\tau_{ij}}{\partialx_k}-\tau_{kj}\frac{\partialu_i}{\partialx_k}-\tau_{ik}\frac{\partialu_j}{\partialx_k}。守恒型的控制方程包括连续性方程\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{u})=0，其中\rho为流体密度，以及动量方程\frac{\partial(\rho\vec{u})}{\partialt}+\nabla\cdot(\rho\vec{u}\vec{u})=-\nablap+\nabla\cdot\boldsymbol{\tau}，其中p为压力。在起动流分析中，假设在初始时刻t=0时，流体静止，即\vec{u}(r,0)=0。在t>0时，突然施加一个恒定的压力梯度G=-\frac{\partialp}{\partialz}，同时内壁面以速度u_w(t)=u_0\sin(\omegat)沿轴向做往复运动，其中u_0为运动幅值，\omega为运动频率。边界条件为：在r=R_1处，u=u_0\sin(\omegat)；在r=R_2处，u=0。为求解上述数学模型，采用固有函数法。首先对控制方程进行无量纲化处理，定义无量纲参数：\bar{r}=\frac{r}{R_2}，\bar{t}=\frac{\mu_0t}{\rhoR_2^2}，\bar{u}=\frac{uR_2}{\mu_0G}，Ha=\frac{\lambda\mu_0}{\rhoR_2^2}（Ha为Hartree数，表示流体弹性大小），A=\frac{u_0}{R_2G}（A为内壁面往复运动的无量纲振幅），Rev=\frac{\omega\rhoR_2^2}{\mu_0}（Rev为振荡雷诺数）。将无量纲参数代入控制方程和边界条件中，得到无量纲化的控制方程和边界条件。然后采用分离变量法，设\bar{u}(\bar{r},\bar{t})=\bar{v}(\bar{r})\bar{T}(\bar{t})，代入无量纲化的控制方程中，得到关于\bar{v}(\bar{r})和\bar{T}(\bar{t})的常微分方程。对于\bar{v}(\bar{r})的常微分方程，其解为\bar{v}(\bar{r})=C_1J_0(\beta\bar{r})+C_2Y_0(\beta\bar{r})，其中J_0和Y_0分别为零阶第一类和第二类贝塞尔函数，\beta为待定常数。根据边界条件\bar{u}(\bar{R}_1,\bar{t})=A\sin(Rev\bar{t})和\bar{u}(\bar{R}_2,\bar{t})=0，可以确定C_1、C_2和\beta的值。对于\bar{T}(\bar{t})的常微分方程，其解为\bar{T}(\bar{t})=C_3e^{-\frac{\beta^2\bar{t}}{Ha}}+C_4。将\bar{v}(\bar{r})和\bar{T}(\bar{t})的解代入\bar{u}(\bar{r},\bar{t})=\bar{v}(\bar{r})\bar{T}(\bar{t})中，得到速度分布\bar{u}(\bar{r},\bar{t})的表达式。计算结果表明，在起动阶段，流体速度随时间逐渐增大，最终达到稳定状态。弹性效应使得流体速度分布更加均匀，且速度幅值相对较小。随着Hartree数的增大，即流体弹性增强，速度达到稳定状态所需的时间变长，且速度分布的均匀性更加明显。例如，当Hartree数从0.1增加到0.5时，在相同时间内，速度分布的最大值与最小值之差减小了约30%，表明速度分布更加均匀。内壁面运动幅值和频率对速度分布也有显著影响。运动幅值增大，速度幅值相应增大；运动频率增大，速度变化更加频繁，且速度幅值略有减小。当运动幅值从0.1增加到0.3时，速度幅值增大了约50%；当运动频率从1增加到3时，速度幅值减小了约20%。3.3.2振荡流与内壁面往复运动非定常流分析在振荡流分析中，假设在入口处施加一个周期性变化的压力梯度G(t)=G_0\sin(\omegat)，其中G_0为压力梯度的幅值，\omega为振荡频率。边界条件为：在r=R_1处，u=0；在r=R_2处，u=0。对于此数学模型，同样采用固有函数法求解。通过分离变量，将控制方程转化为关于径向函数和时间函数的常微分方程。径向函数的解由贝塞尔函数表示，根据边界条件确定贝塞尔函数的系数和特征值。时间函数的解则根据压力梯度的变化形式求解。最终得到速度分布u(r,t)的表达式为：u(r,t)=\sum_{n=1}^{\infty}A_nJ_0(\beta_nr)\left[B_n\cos(\omega_nt)+C_n\sin(\omega_nt)\right]其中A_n、B_n和C_n为待定系数，\beta_n为特征值，\omega_n为与振荡频率相关的参数。计算结果显示，速度随时间做周期性变化，且在一个周期内，速度分布呈现出复杂的形态。弹性效应使得速度分布的峰值减小，分布更加平缓。例如，当弹性参数（如松弛时间\lambda）增大时，速度分布的峰值降低了约25%，同时速度分布在环空内更加均匀。振荡频率对速度分布影响显著，随着振荡频率增加，速度变化的周期变短，速度幅值减小。当振荡频率从1Hz增加到3Hz时，速度幅值减小了约30%。在内壁面做轴向往复运动的非定常流分析中，假设内壁面以速度u_w(t)=u_0\sin(\omegat)沿轴向做往复运动，其中u_0为运动幅值，\omega为运动频率。边界条件为：在r=R_1处，u=u_0\sin(\omegat)；在r=R_2处，u=0。采用与振荡流类似的求解方法，得到速度分布u(r,t)的表达式。分析结果表明，内壁面的往复运动对环空内流体的速度分布产生了明显的影响。靠近内壁面处，速度变化较为剧烈，且速度幅值较大；随着径向距离的增加，速度变化逐渐平缓，幅值逐渐减小。弹性效应使得流体对内壁面运动的响应发生改变，速度分布的波动减小。当弹性增强时，速度分布的波动幅度减小了约20%。运动幅值和频率的增加都会使速度幅值增大，但运动频率的增加会使速度变化更加频繁，导致速度分布的波动加剧。当运动幅值从0.1m/s增加到0.3m/s时，速度幅值增大了约60%；当运动频率从2Hz增加到4Hz时，速度分布的波动幅度增大了约35%。四、聚驱井采出液在偏心环空内的非定常流分析4.1幂律流体偏心环空内的非定常流4.1.1数学模型建立在研究幂律流体在偏心环空内的非定常流时，构建准确的物理模型是基础。考虑由半径为R_1的内管和半径为R_2的外管组成的偏心环形通道，内管与外管的偏心距为e。通道内充满幂律流体，假设流体不可压缩，且流动为层流状态。基于上述物理模型，给出相应的控制方程和定解条件。幂律流体的本构方程为\tau=K(\dot{\gamma})^{n-1}\dot{\gamma}，其中\tau为剪切应力张量，K为稠度系数，n为幂律指数，\dot{\gamma}为剪切速率张量，\dot{\gamma}=\frac{1}{2}(\nabla\vec{u}+(\nabla\vec{u})^T)，\vec{u}为速度矢量。守恒型的控制方程包括连续性方程\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{u})=0，其中\rho为流体密度，以及动量方程\frac{\partial(\rho\vec{u})}{\partialt}+\nabla\cdot(\rho\vec{u}\vec{u})=-\nablap+\nabla\cdot\boldsymbol{\tau}，其中p为压力。定解条件如下：在初始时刻t=0时，速度\vec{u}(r,0)=0；在r=R_1处，若内管做轴向往复运动，速度u(R_1,t)=u_0\sin(\omegat)，若研究进口流速正弦变化的振荡流，此处速度条件根据具体入口条件确定；在r=R_2处，速度u(R_2,t)=0。为了更方便地求解偏心环空内的流动问题，引入双极坐标系(\xi,\eta,z)。在双极坐标系下，坐标变换关系为x=\frac{a\sinh\eta}{\cosh\eta-\cos\xi}，y=\frac{a\sin\xi}{\cosh\eta-\cos\xi}，z=z，其中a为与偏心距相关的常数。将控制方程和定解条件在双极坐标系下进行转换。速度矢量\vec{u}在双极坐标系下表示为\vec{u}=u_{\xi}\vec{e}_{\xi}+u_{\eta}\vec{e}_{\eta}+u_{z}\vec{e}_{z}，其中\vec{e}_{\xi}、\vec{e}_{\eta}和\vec{e}_{z}分别为双极坐标系下的单位矢量。连续性方程在双极坐标系下的形式为：\frac{\partial(\rhoh_{\xi}h_{\eta}u_{\xi})}{\partial\xi}+\frac{\partial(\rhoh_{\xi}h_{\eta}u_{\eta})}{\partial\eta}+\frac{\partial(\rhoh_{\xi}h_{\eta}u_{z})}{\partialz}=0其中h_{\xi}=\frac{a}{\cosh\eta-\cos\xi}，h_{\eta}=\frac{a}{\cosh\eta-\cos\xi}为拉梅系数。动量方程在双极坐标系下的形式较为复杂，以z方向动量方程为例：\begin{align*}&\rho\left(\frac{\partialu_{z}}{\partialt}+u_{\xi}\frac{\partialu_{z}}{\partial\xi}+u_{\eta}\frac{\partialu_{z}}{\partial\eta}+u_{z}\frac{\partialu_{z}}{\partialz}\right)\\=&-\frac{\partialp}{\partialz}+\frac{1}{h_{\xi}h_{\eta}}\left[\frac{\partial}{\partial\xi}\left(\frac{h_{\eta}}{h_{\xi}}\tau_{\xiz}\right)+\frac{\partial}{\partial\eta}\left(\frac{h_{\xi}}{h_{\eta}}\tau_{\etaz}\right)+\frac{\partial}{\partialz}\left(h_{\xi}h_{\eta}\tau_{zz}\right)\right]\end{align*}其中\tau_{ij}为应力张量在双极坐标系下的分量，可根据幂律流体本构方程在双极坐标系下的形式计算得到。边界条件在双极坐标系下也相应转换。在\xi=\xi_1（对应内管表面）处，u_{z}(\xi_1,\eta,t)=u_0\sin(\omegat)；在\xi=\xi_2（对应外管表面）处，u_{z}(\xi_2,\eta,t)=0。通过以上步骤，建立了双极坐标系下幂律流体在偏心环空内非定常流的数学模型。4.1.2数值求解与结果分析在建立数学模型后，采用有限体积法对其进行数值求解。有限体积法的核心思想是将计算区域划分为一系列控制体积，在每个控制体积上对控制方程进行积分，将偏微分方程转化为代数方程，进而求解。在偏心环空的计算区域划分中，根据双极坐标系的特点，合理布置网格节点，确保网格的正交性和疏密分布能够准确捕捉流场的变化。在数值求解过程中，对控制方程进行离散化处理。以动量方程为例，将其在每个控制体积上进行积分，利用高斯散度定理将体积分转化为面积分，得到离散化的动量方程。离散化后的方程通过迭代求解，常用的迭代算法有SIMPLE算法（Semi-ImplicitMethodforPressure-LinkedEquations）及其改进算法，如SIMPLER算法（SIMPLERevised）、SIMPLEC算法（SIMPLEConsistent）等。这些算法通过引入压力修正方程，实现速度和压力的耦合求解，确保计算结果的准确性和收敛性。以具体算例进行计算分析。设定偏心环空内管半径R_1=0.01m，外管半径R_2=0.02m，偏心距e=0.005m，幂律流体的稠度系数K=0.5Pa\cdots^n，幂律指数n=0.8。当进口流速以正弦规律变化时，设进口流速u_{in}(t)=u_{max}\sin(\omegat)，其中u_{max}=0.5m/s，\omega=2\pirad/s。得到速度分布、平均速度和壁面处应力随时间的变化规律（如图16-18所示）。从速度分布（图16）可以看出，由于偏心环空的结构特点，速度分布呈现出明显的非对称性。在靠近内管的一侧，速度变化较为剧烈，且速度值相对较大；而在靠近外管的一侧，速度变化相对平缓，速度值较小。这是因为偏心环空内，内管附近的流体受到内管的影响更大，流速更快，而外管附近的流体受到的影响较小，流速较慢。平均速度（图17）也随时间做周期性变化，其幅值与进口流速的幅值和频率密切相关。壁面处应力（图18）同样呈现周期性变化，且内管壁面处的应力幅值大于外管壁面处，这是由于内管附近的速度梯度较大，根据幂律流体本构方程，剪切应力与速度梯度相关，因此内管壁面处的应力更大。【此处插入图16：进口流速正弦变化时速度随时间和径向位置的变化】【此处插入图17：进口流速正弦变化时平均速度随时间的变化】【此处插入图18：进口流速正弦变化时壁面处应力随时间的变化】当内壁面做轴向往复运动时，设内管内壁面速度u_w(t)=u_0\sin(\omegat)，其中u_0=0.3m/s，\omega=3\pirad/s。分析速度分布、平均速度和壁面处应力的变化情况（如图19-21所示）。速度分布（图19）显示，在内壁面运动的带动下，靠近内管的流体速度明显增大，且速度分布随时间发生周期性变化。在一个周期内，速度先增大后减小，呈现出与内壁面运动相似的周期性。平均速度（图20）也随时间做周期性波动，其波动范围与内壁面运动的幅值和频率有关。壁面处应力（图21）在一个周期内呈现出明显的变化，内管壁面处的应力变化更为显著，这是由于内壁面的运动直接作用于内管附近的流体，导致此处的剪切应力变化较大。【此处插入图19：内壁面轴向往复运动时速度随时间和径向位置的变化】【此处插入图20：内壁面轴向往复运动时平均速度随时间的变化】【此处插入图21：内壁面轴向往复运动时壁面处应力随时间的变化】通过对不同条件下的计算结果分析可知，环空间距、流体性质（稠度系数K和幂律指数n）、振荡频率\omega和振幅等因素对幂律流体在偏心环空内的非定常流有显著影响。环空间距增大，流体的流通面积增大，速度分布会变得相对均匀，平均速度增大，壁面处应力减小；稠度系数增大，流体的粘性增强，流动阻力增大，速度幅值减小，壁面处应力增大；幂律指数n的变化会改变流体的非牛顿特性，当n减小时，流体的剪切变稀特性更加明显，速度分布和壁面处应力的变化也会相应改变；振荡频率增大，流体在单位时间内受到的扰动次数增多，速度幅值减小，平均速度的波动频率增加；振幅增大，速度和壁面处应力的幅值也随之增大。4.2粘弹性流体偏心环空内的非定常流4.2.1数学模型与数值求解为深入探究粘弹性流体在偏心环空内的非定常流特性，建立准确的物理及数学模型至关重要。考虑一个由半径为R_1的内管和半径为R_2的外管组成的偏心环形通道，内管与外管的偏心距为e，通道内充满粘弹性流体。假设流体不可压缩，且流动为层流状态。粘弹性流体的本构方程采用上随体Maxwell模型，其表达式为\tau_{ij}+\lambda\frac{\delta\tau_{ij}}{\deltat}=\mu_0A_{ij}。其中\tau_{ij}为偏应力张量，\lambda为松弛时间，\mu_0为零剪切黏度，A_{ij}为一阶Rivlin-Erickson张量，\frac{\delta\tau_{ij}}{\deltat}为上随体导数，其表达式为\frac{\delta\tau_{ij}}{\deltat}=\frac{\partial\tau_{ij}}{\partialt}+u_k\frac{\partial\tau_{ij}}{\partialx_k}-\tau_{kj}\frac{\partialu_i}{\partialx_k}-\tau_{ik}\frac{\partialu_j}{\partialx_k}。守恒型的控制方程包括连续性方程\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{u})=0，其中\rho为流体密度，以及动量方程\frac{\partial(\rho\vec{u})}{\partialt}+\nabla\cdot(\rho\vec{u}\vec{u})=-\nablap+\nabla\cdot\boldsymbol{\tau}，其中p为压力。定解条件如下：在初始时刻t=0时，速度\vec{u}(r,0)=0；在r=R_1处，若内管做轴向往复运动，速度u(R_1,t)=u_0\sin(\omegat)，若研究进口流速正弦变化的振荡流，此处速度条件根据具体入口条件确定；在r=R_2处，速度u(R_2,t)=0。为简化求解过程，引入双极坐标系(\xi,\eta,z)，在双极坐标系下，坐标变换关系为x=\frac{a\sinh\eta}{\cosh\eta-\cos\xi}，y=\frac{a\sin\xi}{\cosh\eta-\cos\xi}，z=z，其中a为与偏心距相关的常数。将控制方程和定解条件在双极坐标系下进行转换。速度矢量\vec{u}在双极坐标系下表示为\vec{u}=u_{\xi}\vec{e}_{\xi}+u_{\eta}\vec{e}_{\eta}+u_{z}\vec{e}_{z}，其中\vec{e}_{\xi}、\vec{e}_{\eta}和\vec{e}_{z}分别为双极坐标系下的单位矢量。连续性方程在双极坐标系下的形式为：\frac{\partial(\rhoh_{\xi}h_{\eta}u_{\xi})}{\partial\xi}+\frac{\partial(\rhoh_{\xi}h_{\eta}u_{\eta})}{\partial\eta}+\frac{\partial(\rhoh_{\xi}h_{\eta}u_{z})}{\partialz}=0其中h_{\xi}=\frac{a}{\cosh\eta-\cos\xi}，h_{\eta}=\frac{a}{\cosh\eta-\cos\xi}为拉梅系数。动量方程在双极坐标系下的形式较为复杂，以z方向动量方程为例：\begin{align*}&\rho\left(\frac{\partialu_{z}}{\partialt}+u_{\xi}\frac{\partialu_{z}}{\partial\xi}+u_{\eta}\frac{\partialu_{z}}{\partial\eta}+u_{z}\frac{\partialu_{z}}{\partialz}\right)\\=&-\frac{\partialp}{\partialz}+\frac{1}{h_{\xi}h_{\eta}}\left[\frac{\partial}{\partial\xi}\left(\frac{h_{\eta}}{h_{\xi}}\tau_{\xiz}\right)+\frac{\partial}{\partial\eta}\left(\frac{h_{\xi}}{h_{\eta}}\tau_{\etaz}\right)+\frac{\partial}{\partialz}\left(h_{\xi}h_{\eta}\tau_{zz}\right)\right]\end{align*}其中\tau_{ij}为应力张量在双极坐标系下的分量，可根据上随体Maxwell本构方程在双极坐标系下的形式计算得到。边界条件在双极坐标系下也相应转换。在\xi=\xi_1（对应内管表面）处，u_{z}(\xi_1,\eta,t)=u_0\sin(\omegat)；在\xi=\xi_2（对应外管表面）处，u_{z}(\xi_2,\eta,t)=0。在数值求解方面，采用有限体积法对控制方程进行离散化处理。将计算区域划分为一系列控制体积，在每个控制体积上对控制方程进行积分，将偏微分方程转化为代数方程。离散化后的方程通过迭代求解，常用的迭代算法有SIMPLE算法（Semi-ImplicitMethodforPressure-LinkedEquations）及其改进算法，如SIMPLER算法（SIMPLERevised）、SIMPLEC算法（SIMPLEConsistent）等。在迭代过程中，通过不断调整速度和压力的迭代值，使得离散化方程满足收敛条件，从而得到稳定的数值解。为确保数值计算的准确性和可靠性，对计算结果进行网格无关性检验。通过逐步加密网格，对比不同网格数量下的计算结果。当网格数量增加到一定程度后，计算结果不再发生明显变化，此时认为计算结果与网格无关，所采用的网格数量满足计算精度要求。同时，对计算结果进行误差分析，通过与理论解或实验数据进行对比，评估计算结果的误差范围。若误差在可接受范围内，则表明数值计算方法和结果是可靠的；若误差较大，则需要进一步优化数值计算方法或调整计算参数，以提高计算结果的准确性。4.2.2计算结果与影响因素分析通过数值计算，得到了粘弹性流体在偏心环空内非定常流的详细结果，包括速度分布、压力分布、应力分布等，并对这些结果进行了深入分析，探讨了不同因素对流动的影响。当压力梯度呈正弦变化时，即G(t)=G_0\sin(\omegat)，其中G_0为压力梯度的幅值，\omega为振荡频率。速度分布呈现出复杂的周期性变化（如图22所示）。在一个周期内，速度先随着压力梯度的增大而增大，达到峰值后又随着压力梯度的减小而减小。靠近内管处的速度变化较为剧烈，这是因为内管附近的流体受到压力变化的影响更为直接。弹性效应使得速度分布更加均匀，与牛顿流体相比，粘弹性流体的速度峰值相对较小。例如，在相同的压力梯度条件下，牛顿流体的速度峰值为u_{max1}，而粘弹性流体的速度峰值为u_{max2}，u_{max2}约为u_{max1}的80%。【此处插入图22：压力梯度正弦变化时速度随时间和径向位置的变化】平均速度也随时间做周期性变化（如图23所示），其幅值与压力梯度的幅值和振荡频率密切相关。压力梯度幅值增大，平均速度的幅值也随之增大；振荡频率增大，平均速度的变化周期变短，幅值略有减小。当压力梯度幅值从G_{01}增大到G_{02}时，平均速度的幅值增大了约30%；当振荡频率从\omega_1增加到\omega_2时，平均速度的幅值减小了约15%。【此处插入图23：压力梯度正弦变化时平均速度随时间的变化】壁面处应力同样呈现周期性变化（如图24所示），内管壁面处的应力幅值大于外管壁面处。这是由于内管附近的速度梯度较大，根据粘弹性流体本构方程，剪切应力与速度梯度相关，因此内管壁面处的应力更大。弹性效应使得壁面处应力的变化相对平缓，与牛顿流体相比，粘弹性流体壁面处应力的波动范围较小。例如，在相同的流动条件下，牛顿流体壁面处应力的波动范围为\Delta\tau_1，而粘弹性流体壁面处应力的波动范围为\Delta\tau_2，\Delta\tau_2约为\Delta\tau_1的70%。【此处插入图24：压力梯度正弦变化时壁面处应力随时间的变化】当内壁面沿轴向做往复运动时，设内管内壁面速度u_w(t)=u_0\sin(\omegat)，其中u_0为运动幅值，\omega为运动频率。速度分布在靠近内壁面处变化剧烈，随着径向距离的增加，速度变化逐渐平缓（如图25所示）。在内壁面运动的带动下，靠近内管的流体速度明显增大，且速度分布随时间发生周期性变化。在一个周期内，速度先增大后减小，呈现出与内壁面运动相似的周期性。弹性效应使得流体对内壁面运动的响应发生改变，速度分布的波动减小。当弹性增强时，速度分布的波动幅度减小了约20%。【此处插入图25：内壁面轴向往复运动时速度随时间和径向位置的变化】平均速度也随时间做周期性波动（如图26所示），其波动范围与内壁面运动的幅值和频率有关。内壁面运动幅值增大，平均速度的幅值相应增大；运动频率增大，平均速度的波动频率增加，幅值略有减小。当内壁面运动幅值从u_{01}增加到u_{02}时，平均速度的幅值增大了约40%；当运动频率从\omega_1增加到\omega_2时，平均速度的幅值减小了约18%。【此处插入图26：内壁面轴向往复运动时平均速度随时间的变化】壁面处应力在一个周期内呈现出明显的变化（如图27所示），内管壁面处的应力变化更为显著。这是由于内壁面的运动直接作用于内管附近的流体，导致此处的剪切应力变化较大。弹性效应使得壁面处应力的变化相对稳定，与牛顿流体相比，粘弹性流体壁面处应力的波动范围较小。例如，在相同的运动条件下，牛顿流体壁面处应力的波动范围为\Delta\tau_3，而粘弹性流体壁面处应力的波动范围为\Delta\tau_4，\Delta\tau_4约为\Delta\tau_3的65%。【此处插入图27：内壁面轴向往复运动时壁面处应力随时间的变化】影响粘弹性流体在偏心环空内非定常流的因素众多，包括环空间距、流体性质（松弛时间\lambda、零剪切黏度\mu_0）、振荡频率\omega和振幅等。环空间距增大，流体的流通面积增大，速度分布会变得相对均匀，平均速度增大，壁面处应力减小；松弛时间增大，流体的弹性增强，速度分布更加均匀，壁面处应力的波动范围减小；零剪切黏度增大，流体的粘性增强，流动阻力增大，速度幅值减小，壁面处应力增大；振荡频率增大，流体在单位时间内受到的扰动次数增多，速度幅值减小，平均速度的波动频率增加；振幅增大，速度和壁面处应力的幅值也随之增大。五、聚驱机采井杆管偏磨的机理分析5.1聚驱井中轴向力作用下的杆管偏磨5.1.1采出液宏观流变性与井筒流动影响聚驱采出液的宏观流变性是影响井筒内流动的关键因素，而其粘度特性在其中起着核心作用。聚驱采出液中含有聚合物，使其呈现出复杂的流变特性。通过大量实验研究发现，聚驱采出液的粘度与聚合物浓度密切相关。当聚合物浓度较低时，采出液的粘度相对较小；随着聚合物浓度的逐渐增加，采出液的粘度显著增大。例如，在某聚驱油田的实验中，当聚合物浓度从500mg/L增加到1500mg/L时，采出液在低剪切速率下的粘度从5mPa・s左右迅速增大到30mPa・s以上。这是因为聚合物分子在溶液中形成了网状结构，随着浓度的增加，分子间的相互作用增强，导致粘度大幅上升。采出液的粘度还受温度影响显著。温度升高时，聚合物分子的热运动加剧，分子间的相互作用力减弱，使得采出液的粘度降低。以某聚合物驱采出液为例，在25℃时，其粘度为20mPa・s，当温度升高到50℃时，粘度下降至10mPa・s左右。这种粘度随温度的变化特性对井筒内流动产生重要影响，在不同温度条件下，采出液的流动阻力和流动状态会发生改变。剪切速率对聚驱采出液的粘度也有明显影响，呈现出剪切变稀特性。当剪切速率较低时，聚合物分子链相互缠绕，采出液粘度较高；随着剪切速率的增大，分子链逐渐解缠结并取向，粘度随之降低。例如，在剪切速率从0.1s⁻¹增加到10s⁻¹的过程中，某聚驱采出液的粘度从40mPa・s降至15mPa・s。这种剪切变稀特性使得采出液在井筒内的流动特性更加复杂，不同部位的流速和粘度分布会因剪切速率的变化而不同。聚合物溶液的宏观流变性不仅体现在粘度特性上，还表现出粘弹性。聚合物分子链的长链结构使其在流动过程中能够储存和释放弹性变形能。在低剪切速率下，聚合物溶液主要表现为粘性，随着剪切速率的增加，弹性效应逐渐显现。例如，在一些实验中，当剪切速率达到一定值时，聚合物溶液会出现明显的弹性现象，如爬杆效应和射流胀大现象。这种粘弹性对井筒内流动的影响不可忽视，它会改变流体的速度分布和压力分布，进而影响杆管的受力情况。为了深入研究采出液性质变化对井筒内流动的影响，建立准确的物理和数学模型至关重要。物理模型考虑了井筒的实际结构，包括油管、抽油杆以及其间的环形空间，同时考虑了采出液的多相组成和流变特性。数学模型基于流体力学的基本方程，如连续性方程、动量方程和能量方程，结合聚驱采出液的本构方程进行建立。对于粘弹性聚驱采出液，常采用上随体Maxwell模型作为本构方程，以准确描述其复杂的流变行为。通过数值模拟方法求解建立的数学模型，得到了一系列关于井筒内流动的重要结果。在速度分布方面，由于采出液的非牛顿特性和井筒结构的影响，速度分布呈现出非均匀性。靠近抽油杆和油管壁面处的速度梯度较大，而在环空中心区域速度相对较为均匀。例如，在某模拟计算中，靠近抽油杆壁面处的速度为0.5m/s，而在环空中心区域速度达到1.0m/s。这种速度分布的不均匀性会导致采出液对杆管的作用力不均匀，从而增加杆管偏磨的风险。压力分布也受到采出液性质和井筒结构的影响。在井筒内，压力沿轴向逐渐降低，且在不同径向位置存在压力差。这是由于采出液的粘性和流动阻力导致的压力损失，以及