油气井流体力学教学课件PPT固液两相流动 91页

'第三章　环空固液两相流动颗粒沉降速度牛顿流体非牛顿流体钻井液的携带能力Moore关系式Chien模式Walker和Mayes关系式岩屑传输率两层稳定运移模型岩屑运移机理两层运移模型岩屑床厚度经验模式实验研究主要内容1 第三章　环空固液两相流动两层不稳定运移模型模型假设与时间相关方程的建立差分离散计算结果Chien模式Walker和Mayes关系式三层稳定运移模型几何模型数学模型的建立传输方式和通用方程模拟结果主要内容2 第三章　环空固液两相流动第一节颗粒沉降速度1.1牛顿流体重力W、浮力Fbo和粘性阻力F平衡粘性阻力F与球形颗粒的沉降速度之间的关系(stokes)则颗粒沉速为：适用于雷诺数小于0.1的流动（工程单位）3 第三章　环空固液两相流动第一节颗粒沉降速度1.1牛顿流体对于雷诺数大于0.1的流动，引入摩擦系数的概念：颗粒沉降速度公式可用于雷诺数小于0.1的情况。4 第三章　环空固液两相流动第一节颗粒沉降速度1.1牛顿流体形状圆球度球体1.0八面体0.85立方体0.81圆柱体h=r/15h=r/10h=r/3h=rh=2rh=3rh=10rh=20r0.250.320.590.830.870.960.690.585 第三章　环空固液两相流动第一节颗粒沉降速度1.1牛顿流体例1钻井液的密度为8.33lbm/gal，粘度为1cp，含沙量约为1%（体积分数）。如果钻井液停止循环30分钟，将有多少砂砾沉降到井底？假定砂砾的平均粒径为0.0025in，圆球度为0.81，比重为2.6,沙砾沉降物的孔隙度为0.4。解：6 第三章　环空固液两相流动第一节颗粒沉降速度1.1牛顿流体7 第三章　环空固液两相流动第一节颗粒沉降速度1.1牛顿流体例1钻井液的密度为8.33lbm/gal，粘度为1cp，含沙量约为1%（体积分数）。如果钻井液停止循环30分钟，将有多少砂砾沉降到井底？假定砂砾的平均粒径为0.0025in，圆球度为0.81，比重为2.6,沙砾沉降物的孔隙度为0.4。解：=如果钻井液停止循环30min，井筒下方沉降的沙砾高度为303ft。如果那么井底沙砾沉降物的高度为：8 第三章　环空固液两相流动第一节颗粒沉降速度1.2非牛顿流体9 第三章　环空固液两相流动第一节颗粒沉降速度1.2非牛顿流体微粒直径必须大于例2泥浆密度为9lbm/gal，胶凝作用力为5lbm/100sqft，沙砾比重为2.6。试计算能在钻井液中悬浮的沙砾的最大直径。10 第三章　环空固液两相流动第一节颗粒沉降速度1.4Moore关系式非牛顿流体牛顿流体11 第三章　环空固液两相流动第一节颗粒沉降速度1.4Moore关系式当流核雷诺数≤3时，流体处于层流状态，对于过渡流雷诺数，12 第三章　环空固液两相流动第一节颗粒沉降速度1.5Chien模式对于聚合物钻井液：对于粘土在水中悬浮问题，推荐使用塑性粘度代替表观粘度。在颗粒雷诺数小于100时：13 第三章　环空固液两相流动第一节颗粒沉降速度1.6Walker和Mayes关系式Walker和Mayes提出了用磨擦系数来求滑脱速度的关系式，适用于圆片。颗粒雷诺数大于100时，认为流动模式是紊流，并假定f=1.12，则：14 第三章　环空固液两相流动第一节颗粒沉降速度1.6Walker和Mayes关系式雷诺数小于100时则用于确定颗粒雷诺数的表观粘度由下式得到：15与=刻度盘读数×1.066=电机速度×1.03 第二节岩屑传输率岩屑向地面运移的速度为流体流动速度和颗粒沉降速度之差。传输速度：颗粒相对于地面移动的速度传输率为传输速度与平均环空返速之比：传输率为正值时，岩屑被携带到地面。当颗粒沉降速度为零时，岩屑平均速率与平均环层速率相等，岩屑传输率为整数。随着颗粒沉降速度的增加，传输率减小，岩屑在环空向地面传输途中浓度增大。所以，岩屑传输率是钻井液携岩能力的量度。第三章　环空固液两相流动16 第二节岩屑传输率预测岩屑传输率的各种方法对比第三章　环空固液两相流动17 第二节岩屑传输率各种方法预测岩屑传输效率的误差直方图第三章　环空固液两相流动18 第二节岩屑传输率Moore和Chien关系式预测下滑速度的误差直方图第三章　环空固液两相流动19 通过提高泥浆泵排量可以降低环空泥浆平均密度，从而提高岩屑传输效率。但是，随着泥浆泵排量增加到达一定值后，井底压力会随着流量的增加而增大，这是因为在环空产生了过大的磨擦压力损失。同时，由于钻柱中的过大摩擦压力损失，钻头处的喷射冲击力开始随着流量的增大而降低。因此，存在一个最优化流量，在该流量下，钻进每英尺理论费用最小。第二节岩屑传输率第三章　环空固液两相流动20 第二节岩屑传输率岩屑传输优化结果实例第三章　环空固液两相流动21 第二节岩屑传输率泥浆粘度对优化岩屑传输效率和降低每英尺费用的影响实例调整流体性质（粘度、密度）及增加环空返速可增大岩屑传输效率。在多数情况下，使用低粘度流体可实现每英尺的理论费用。第三章　环空固液两相流动22 第二节岩屑传输率例3.3计算传输效率，岩屑直径0.25in，比重2.6(21.6-lbm/gal)。粘土钻井液密度为9.01bm/gal，其环空返速为120ft/min(2.0ft/s)，环空尺寸为10×5in。分别用Moore、Chien、Walker和Mayes关系式求解。从旋转粘度计上读得参数如右表。假定岩屑直径及厚度均约为0.25in。解：Moore模式稠度系数K和流性指数N根据300转和600转的粘度计读数求得：第三章　环空固液两相流动转子转速标度盘读数(rpm)(degree)3206331001320022300306005023 例3.3计算传输效率，岩屑直径0.25in，比重2.6(21.6-lbm/gal)。粘土钻井液密度为9.01bm/gal，其环空返速为120ft/min(2.0ft/s)，环空尺寸为10×5in。分别用Moore、Chien、Walker和Mayes关系式求解。从旋转粘度计上读得参数如上表。假定岩屑直径及厚度均约为0.25in。第二节岩屑传输率环空返速为120ft/min(2ft/s)时的表观牛顿粘度为：第三章　环空固液两相流动24 例3.3计算传输效率，岩屑直径0.25in，比重2.6(21.6-lbm/gal)。粘土钻井液密度为9.01bm/gal，其环空返速为120ft/min(2.0ft/s)，环空尺寸为10×5in。分别用Moore、Chien、Walker和Mayes关系式求解。从旋转粘度计上读得参数如上表。假定岩屑直径及厚度均约为0.25in。第二节岩屑传输率假定是过渡态雷诺数（过渡流型），下滑速度为雷诺数由下式求得：第三章　环空固液两相流动25 例3.3计算传输效率，岩屑直径0.25in，比重2.6(21.6-lbm/gal)。粘土钻井液密度为9.01bm/gal，其环空返速为120ft/min(2.0ft/s)，环空尺寸为10×5in。分别用Moore、Chien、Walker和Mayes关系式求解。从旋转粘度计上读得参数如上表。假定岩屑直径及厚度均约为0.25in。第二节岩屑传输率既然雷诺数介于3～300之间，则求解正确。故在环空流速2.0ft/s、岩屑下滑速度为0.49ft/s条件下的岩屑传输效率为：第三章　环空固液两相流动26 例3.3计算传输效率，岩屑直径0.25in，比重2.6(21.6-lbm/gal)。粘土钻井液密度为9.01bm/gal，其环空返速为120ft/min(2.0ft/s)，环空尺寸为10×5in。分别用Moore、Chien、Walker和Mayes关系式求解。从旋转粘度计上读得参数如上表。假定岩屑直径及厚度均约为0.25in。第二节岩屑传输率Chien模式对粘土∕水泥浆，Chien推荐用塑性粘度代替表观粘度。用范氏粘度计的300转和600转读数得到：第三章　环空固液两相流动27 第二节岩屑传输率假定为过渡态，下滑速度可由公式（3.11）求得：例3.3计算传输效率，岩屑直径0.25in，比重2.6(21.6-lbm/gal)。粘土钻井液密度为9.01bm/gal，其环空返速为120ft/min(2.0ft/s)，环空尺寸为10×5in。分别用Moore、Chien、Walker和Mayes关系式求解。从旋转粘度计上读得参数如上表。假定岩屑直径及厚度均约为0.25in。第三章　环空固液两相流动28 第二节岩屑传输率例3.3计算传输效率，岩屑直径0.25in，比重2.6(21.6-lbm/gal)。粘土钻井液密度为9.01bm/gal，其环空返速为120ft/min(2.0ft/s)，环空尺寸为10×5in。分别用Moore、Chien、Walker和Mayes关系式求解。从旋转粘度计上读得参数如上表。假定岩屑直径及厚度均约为0.25in。则此速度下的颗粒雷诺数为：既然这个雷诺数小于100，则结果正确，传输效率由下式给出：第三章　环空固液两相流动29 第二节岩屑传输率例3.3计算传输效率，岩屑直径0.25in，比重2.6(21.6-lbm/gal)。粘土钻井液密度为9.01bm/gal，其环空返速为120ft/min(2.0ft/s)，环空尺寸为10×5in。分别用Moore、Chien、Walker和Mayes关系式求解。从旋转粘度计上读得参数如上表。假定岩屑直径及厚度均约为0.25in。Walker-Mayer模式在静止液体中颗粒下沉的剪切力估算：14lbm/100sqft的剪切应力对应粘度计读数为：第三章　环空固液两相流动转子转速标度盘读数(rpm)(degree)3206331001320022300306005030 第二节岩屑传输率例3.3计算传输效率，岩屑直径0.25in，比重2.6(21.6-lbm/gal)。粘土钻井液密度为9.01bm/gal，其环空返速为120ft/min(2.0ft/s)，环空尺寸为10×5in。分别用Moore、Chien、Walker和Mayes关系式求解。从旋转粘度计上读得参数如上表。假定岩屑直径及厚度均约为0.25in。假定流型为过渡态，下滑速度由下式求得：第三章　环空固液两相流动31 第二节岩屑传输率例3.3计算传输效率，岩屑直径0.25in，比重2.6(21.6-lbm/gal)。粘土钻井液密度为9.01bm/gal，其环空返速为120ft/min(2.0ft/s)，环空尺寸为10×5in。分别用Moore、Chien、Walker和Mayes关系式求解。从旋转粘度计上读得参数如上表。假定岩屑直径及厚度均约为0.25in。此下滑速度下的颗粒雷诺数为：岩屑传输效率由下式给出：第三章　环空固液两相流动32 第三章　环空固液两相流动研究概况33 第三章　环空固液两相流动34 第三章　环空固液两相流动35 3.1两层稳定模型假设在一定环空返速下，因倾斜或水平环空内的快速沉降而使绝大多数岩屑颗粒聚集在环空底部形成岩屑床；只有少数颗粒在液流中作悬浮运动。假定岩屑床的厚度均匀，且其中岩屑体积浓度为52%。假定岩屑床面以上颗粒作悬浮运动，不产生滚动、翻转、滑脱。固液体系介质均不可压缩，且固液两相无质量交换。岩屑颗粒具有相同的直径和圆球度。钻杆偏心但不旋转。两层模型示意图第三章　环空固液两相流动第三节两层稳定运移模型36 连续性方程：运动方程：岩屑扩散方程：受阻颗粒沉速：第三章　环空固液两相流动第三节两层稳定运移模型3.1两层稳定模型37 水动力摩擦力：岩屑床中流体流动压降：第三章　环空固液两相流动第三节两层稳定运移模型3.1两层稳定模型38 理论岩屑床厚度与环空返速间的关系图在相同环空返速下，环空压降与井斜角间的关系图第三章　环空固液两相流动第三节两层稳定运移模型3.1两层稳定模型39 第三章　环空固液两相流动第三节两层稳定运移模型3.2岩屑床厚度经验模式40 模型假设：(1)底层为岩屑床(岩屑由于重力的作用而沉积在环空低侧)。此层中的岩屑体积浓度假定为52%；(2)顶层为悬浮层，岩屑颗粒的径向扩散瞬时达到平衡；(3)在底层与顶层中，固液两相之间均无滑脱，固液两相之间无质量交换；(4)固液体系均为不可压缩介质，固体颗粒具有相同的直径和圆球度；(5)底层与顶层的界面高度不变，横截面上的流体压力视为静压分布；(6)钻柱偏心但不旋转。模型示意图：第三章　环空固液两相流动第四节两层不稳定运移模型4.1模型假设41 连续性方程运动方程辅助方程几何关系式、各种剪切应力、摩擦系数、以及各层密度的关系式悬浮层颗粒浓度的求法同前。边界条件第三章　环空固液两相流动第四节两层不稳定运移模型4.2模型方程42 低粘度的钻井液在井眼中达到稳态大约需1200s左右，而高粘度的钻井液在2000s后还在冲蚀岩屑床。从图中也可以看出，高粘度的钻井液在稳态时所形成的岩屑床比低粘度的钻井液所形成的岩屑床低一些。时间距离钻头长度岩屑床高度时间距离钻头长度岩屑床高度第三章　环空固液两相流动第四节两层不稳定运移模型4.3计算结果43 低速下只有部分岩屑被除去。随着钻井液流速的增加(右图)，岩屑床可以完全被除去。从图中可以看出岩屑的沙丘传输机制。在6000s后，井段的大部分岩屑均被除去。因此，在一定流速条件下，长时间循环钻井液可以完全除去岩屑床。时间距离钻头长度岩屑床高度时间距离钻头长度岩屑床高度第三章　环空固液两相流动第四节两层不稳定运移模型4.3计算结果44 由于研究和计算的问题较为复杂，因此对模型作了如下假设：（1）钻井液为不可压缩流体，除悬浮层外，不考虑其它各层内固液两相之间的滑动；（2）钻井液和岩屑物性连续；（3）悬浮层岩屑浓度相对很小，且符合扩散定律；（4）均匀层岩屑体积浓度为55%，分散层浓度为均匀层浓度的0.8倍。第三章　环空固液两相流动45 连续性方程动量方程第三章　环空固液两相流动46 动量方程连续性方程■■■■■■■■模型求解有限差分：使用交错网格和一阶迎风格式差分方程组第三章　环空固液两相流动第三章　环空固液两相流动47 （1）划分网格，给定初始值和边界条件；（2）求解动量方程预测步，得到；（3）将代入动量方程修正步，得到的和方程；（4）将以上方程代入连续性方程预测步，求解出、和；（5）将压力代入动量方程修正步，得到速度；（6）将速度代入基本连续性方程，求解出正确的和；（7）求解下一时间步至结束。由以上方程和求解过程可以看出，每一步均是求解比较简单的一次线性方程或方程组，无须迭代计算，因此该方法可以有效的提高计算效率，节省计算时间。SETS数值求解SETS方法，英文全称为StabilityEnhancingTwoStep方法，包括预测和修正两步，在每一时间步无需迭代，因此能够有效提高计算效率。第三章　环空固液两相流动48 现有不稳定模型三层模型层数两层三层岩屑运移方式悬浮、层移悬浮、滑动、滚动、层移应用范围单一井段全井段钻杆旋转不考虑考虑泥浆流变模式幂律幂律、宾汉、卡森、赫巴适应性两层、单层三层、拟三层、两层、单层求解SIMPLE方法（非线性）SETS方法（线性）井底边界静态动态实钻模拟不能能第三章　环空固液两相流动49 以大港ES-H6大位移井为例，其井身结构如图所示。第三章　环空固液两相流动50 大位移井钻井实践表明，216mm（8-1/2in）井段虽然也是大斜度井段，但由于井眼尺寸小，泵排量可以满足要求，因此该井段一般不会出现携岩问题。而311mm（12-1/4in）井段同样是大斜度井段，但由于环空尺寸大，泵排量一般达不到携岩要求，因此该井段钻进时携岩问题最为突出。因此，本研究重点模拟311mm井眼岩屑运移情况。使用139.7mm钻杆，泥浆密度为1.1g/cm3，岩屑密度为2.5g/cm3，岩屑粒径为5mm，泥浆性能参数为n=0.5，K=0.8Pa﹒sn，钻杆转速默认为130r/min，ROP为50m/h，偏心度为0.5，默认排量为60L/s。311mm井段从测深1350m连续钻进至4050m，然后停钻洗井，不考虑倒划眼、起下钻等过程。第三章　环空固液两相流动51 311mm井段从测深1350m连续钻进至4050m，然后停钻洗井。Q=60L/s，钻进用时54h，洗井用时20h。第三章　环空固液两相流动52 岩屑床高度分布变化图中绿线表示钻进过程，红线表示停钻洗井过程。Q=50L/s钻进过程，钻头逐渐后移洗井过程，岩屑床逐渐被带出第三章　环空固液两相流动53 岩屑床高度分布变化Q=70L/s第三章　环空固液两相流动54 岩屑床高度分布变化Q=90L/s第三章　环空固液两相流动55 井底ECD变化钻进过程洗井过程第三章　环空固液两相流动56 钻杆旋转影响Q=60L/s钻进时洗井时第三章　环空固液两相流动57 岩屑床高度工程计算模型敏感性排序：排量、泥浆有效粘度、井眼内径和钻杆外径、井斜角、偏心度、钻杆钻速、钻头直径、机械钻速、泥浆密度、岩屑直径、岩屑密度。第三章　环空固液两相流动58 实测ECD计算理论无岩屑床ECD理论有岩屑床ECD井眼清洁的影响因素第三章　环空固液两相流动59 第三章　环空固液两相流动60 实测ECD计算理论无岩屑床ECD研究和实践表明，大位移井因其大斜度段和水平段长，易形成岩屑床，当环空中存在岩屑床的时候，即使岩屑床厚度很小，环空压耗也会有大幅度增加，因此可以利用环空压耗监测井眼清洁程度。工程上常将ECD作为一个重要工作参数进行监测，因此可以将立管压力转换为ECD，基于岩屑床动态运移工程计算模型，通过实测ECD、理论无岩屑床和存在一定高度岩屑床时ECD的对比分析，可以对环空平均岩屑床高度进行定量评价，实时监测井眼清洁程度。实测ECD理论无岩屑床ECD理论有岩屑床ECD第三章　环空固液两相流动61 第三章　环空固液两相流动62 第三章　环空固液两相流动63 第三章　环空固液两相流动64 第三章　环空固液两相流动65 （1）17-1/2”井段钻至1342m，井斜角从0°变化到77.88°，下13-3/8”套管至1341m，固井。采用海水聚合物泥浆体系，密度为1.05-1.12g/cm3，粘度数据为38/24/19/13/6/2。（2）12-1/4”钻至4026m，井斜角从77.88°变化到87.78°，下9-5/8”套管至4024m，固井。采用Ultradril泥浆体系，密度为1.10-1.20g/cm3，流变参数为51/36/29/20/7/5。（3）8-1/2”井段钻至4590m，井斜角从87.78°变化到88.89°，下7”筛管至4047m。采用Ultradril泥浆体系，密度为1.05-1.15g/cm3，流变参数为55/37/30/21/7/5。以大港油区ES-H6井为例进行分析，其井身结构与泥浆数据为：实例分析66 171/2in井段——岩屑床高度分析实例分析67 171/2in井段——最小排量分析159.0L/S1.14m/s实例分析68 171/2in井段——ECD分析计算表明，岩屑床对ECD的影响很大，由于大斜度段长，现有泵无法达到理论最小排量，环空必然形成一定的岩屑床，因此应严格控制ROP、排量等，以控制大斜度井段岩屑床。实例分析69 171/2in井段——岩屑床高度钻杆转速是井眼清洁的有效和主要控制手段，应尽可能旋转钻杆，破坏岩屑床。ROP减小，环空岩屑床高度降低。实例分析70 171/2in井段——加大洗井频率计算表明，排量为60l/s，ROP为50m/h，洗井间隔为300m时，大斜度井段环空将形成约16%～18%的岩屑床，洗井间隔为150m时，大斜度井段环空将形成约14%～15%的岩屑床，岩屑床高度降低。实例分析71 171/2in井段——循环洗井时间以环空残留1%高度岩屑床为标准，当排量为60L/s，钻杆转速为130rpm时，洗井间隔为150m时，每次洗井所需时间约为60-150min，洗井时间随测深增加而增加。实例分析72 121/4in井段——岩屑床高度计算实例分析73 121/4in井段——最小排量分析80.0L/S1.06m/s实例分析74 121/4in井段——ECD实例分析75 121/4in井段——循环压耗实例分析76 121/4in井段——加大洗井频率计算表明，排量为60l/s，ROP为30m/h，洗井间隔为300m时，大斜度井段环空将形成约11%～17%的岩屑床，洗井间隔为150m时，大斜度井段环空岩屑床将控制在9%～13%。实例分析77 121/4in井段——循环洗井时间计算以环空残留1%高度岩屑床为标准，当排量为60L/s，钻杆转速为130rpm，洗井间隔为150m时，每次洗井所需时间约为70-140min，洗井时间随测深增加而增加。钻杆旋转能有效减少洗井时间，因此循环洗井时应尽可能的保持钻杆高速旋转。实例分析78 81/2in井段——岩屑床高度计算实例分析79 81/2in井段——最小排量分析37.0L/S1.06m/s实例分析80 81/2in井段——ECD实例分析81 81/2in井段——循环压耗实例分析82 81/2in井段——岩屑床高度实例分析83 81/2in井段——加大洗井频率实例分析84 12-1/4”井段ECD分布图12-1/4”井段从1340m钻至4026m，如图所示，实测ECD与无岩屑床时ECD基本完全重合，整个井段钻进过程顺利完成。以大港油区ES-H6井为例进行分析，计算数据同前。实例分析85 8-1/2”井段ECD分布图8-1/2”井段从4026m钻至4590m，由于环空尺寸进一步减小，环空返速进一步增加，该井段携岩较为理想。实际排量基本达到了理论最小排量，井眼清洁良好，无岩屑床存在，实测ECD与无岩屑床时ECD基本完全重合，整个井段钻进过程顺利完成。实例分析86 庄海Es-L1井12-1/4”井段ECD分布图庄海Es-L3井12-1/4”井段ECD图实例分析87 庄海Es-H2井12-1/4”井段ECD图庄海Es-H5井12-1/4”井段ECD图实例分析88 张海34-17井12-1/4”井段ECD图庄海8Ng-H4井12-1/4”井段ECD图实例分析89 庄海8Es-H2井身结构参数及泥浆数据：171/2”（444.5mm）井段钻至1356m，井斜角从0°变化到76.74°，下13-3/8”套管至1354m，固井。采用海水聚合物泥浆体系，密度为1.05-1.10g/cm3，粘度数据为38/24/19/13/6/2。平均机械钻速32m/h。（牙轮钻头）出现问题：在定向井段，由于井眼比较大，环空返速低携岩困难，发生井眼不清洁现象，采用间断打入稠塞70L/s大排量循环的措施，才解决了泥浆携岩的问题。实例分析90 庄海8Es-H2（软件预测）实例分析91'