CN109812237B - 一种满足异形井眼清洁的钻井液排量确定方法 - Google Patents
一种满足异形井眼清洁的钻井液排量确定方法 Download PDFInfo
- Publication number
- CN109812237B CN109812237B CN201711170726.XA CN201711170726A CN109812237B CN 109812237 B CN109812237 B CN 109812237B CN 201711170726 A CN201711170726 A CN 201711170726A CN 109812237 B CN109812237 B CN 109812237B
- Authority
- CN
- China
- Prior art keywords
- drilling fluid
- discharge capacity
- rock debris
- drilling
- special
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Earth Drilling (AREA)
Abstract
本发明提供了一种满足异形井眼清洁的钻井液排量确定方法，属于石油天然气勘探开发钻井领域。该方法所述方法首先基于钻井基础参数建立异形井眼几何模型，设定初始的排量，然后通过数值计算得到稳定状态下井眼各个位置的岩屑浓度分布，通过岩屑浓度分布判断是否满足井眼清洁效果，并判断钻井液排量是否在合理排量范围内，如果有一个为否，则调整钻井液流变性，然后重复上述过程，直到钻井液排量满足井眼清洁效果且在合理排量范围内，此时的钻井液排量即为最终确定的钻井液排量。利用本发明可以克服现有经验法或基于规则井眼的环空水力学法的不足，解决钻井现场仅凭经验或环空水力学法在异形井眼中调整钻井液排量精度不足的问题。
Description
技术领域
本发明属于石油天然气勘探开发钻井领域，具体涉及一种满足异形井眼清洁的钻井液排量确定方法。
背景技术
钻井液是钻井过程中的循环携带岩屑的流体介质。井眼清洁程度是指钻井时井筒中岩屑携带的状态。钻井中的岩屑携带状况直接关系到钻井速度及井下安全等多个方面，是十分重要的钻井参数。钻井过程中，如果岩屑不能及时携带出井眼，可能造成下钻困难、下钻不到位以及卡钻等井下复杂情况。因此井眼清洁成为能否安全、快速钻进的首要问题。
实践表明，井眼的形状和钻井液的排量、流变参数对钻井中的岩屑携带有着显著的影响。在实际钻井过程中，井眼常常出现扩径的异形现象，在异形井眼的扩径位置附近，钻井液的流动状态和岩屑携带情况变得更加复杂，容易产生岩屑堆积的现象。
根据异形井眼中井眼清洁指标的要求，可以对钻井液排量和性能调整做出合理判断，利用现场小型试验，优化钻井液流变参数，调整钻井液排量，可以改善岩屑携带效果，有效增强钻井安全并降低钻井成本。
目前，现场采用的钻井液排量选取方法主要有：(1)仅凭经验调整排量和钻井液流变参数；(2)在井眼规则的假设下，基于环空水力学法计算不同钻井液排量和流变参数下的岩屑携带效果，根据岩屑携带效果是否在合理范围确定是否调整钻井液排量和流变参数。
目前的方法存在着一些不足，主要体现在：(1)仅凭经验调整排量和钻井液流变参数盲目性较大；(2)假定的规则井眼形状可能与实际井眼形状偏差较大，未采用最符合实际井眼形状计算岩屑携带效果，计算精度低，无法准确判断岩屑在井眼中的堆积位置。以上的不足可能导致钻井液排量和流变性能维护调控不当，浪费钻井液材料，增加钻井成本，极易引起井下复杂情况的发生。
发明内容
本发明的目的在于解决上述现有技术中存在的难题，提供一种满足异形井眼清洁的钻井液排量确定方法，解决钻井现场仅凭经验或基于规则井眼的调整钻井液排量的问题。在满足异形井眼清洁要求的前提下，实现现场钻井液排量的合理调控，确定钻井液维护处理措施，满足安全高效钻井的要求并降低钻井成本。
本发明是通过以下技术方案实现的：
一种满足异形井眼清洁的钻井液排量确定方法，所述方法首先基于钻井基础参数建立异形井眼几何模型，设定初始的排量，然后通过数值计算得到稳定状态下井眼各个位置的岩屑浓度分布，通过岩屑浓度分布判断是否满足井眼清洁效果，并判断钻井液排量是否在合理排量范围内，如果有一个为否，则调整钻井液流变性，然后重复上述过程，直到钻井液排量满足井眼清洁效果且在合理排量范围内，此时的钻井液排量即为最终确定的钻井液排量。
所述方法包括：
(1)基于钻井基础参数，建立异形井眼几何模型，以N倍钻具直径的长度为纵向范围，以井眼扩大的最大值为横向范围，确定xy二维计算区域；
(2)根据异形井眼几何模型划分网格；
(3)从振动筛前的高架槽内取钻井液，用旋转粘度计测量钻井液在不同转速Φ下的读数θ，将旋转粘度计的读数转变为相应的一组剪切速率γ和剪切应力τ；根据剪切速率和剪切应力，回归拟合4种流变模型；
(4)，计算各流变模型拟合的相关系数，确定相关系数最大的流变模型为拟合效果最佳的钻井液流变流型；
(5)设定排量Q，设置钻井液流型的入口和出口的边界条件；
(6)设定初始条件：按照循环开始前环空可能出现的最高岩屑浓度设定岩屑浓度α2的初始值；
(7)基于连续方程、动量方程和能量方程计算得到的每个网格的两相混合流速
的值绘制钻井液异形井眼处的流速分布图，根据每个网格的固相流动速度差计算得到的每个网格的岩屑浓度α2的值绘制钻井液异形井眼处的岩屑浓度分布图；其中，α2＝α2，原+v2，流进-v2，流出，其中α2，原是迭代上一步的岩屑浓度，v2，流进是每个网格岩屑的流进速度，v2，流出是每个网格岩屑的流出速度；
(8)如果步骤(7)中得到的每个网格的岩屑浓度值α2值都小于预设的岩屑浓度最大允许值α2，max，则判定步骤(5)中的排量Q值满足井眼清洁要求，转入步骤(9)，否则返回步骤(5)；
(9)判断满足井眼清洁要求的排量Q是否在合理排量范围内，如果否，则添加能够改变钻井液的流变性的钻井液助剂，调整钻井液流变性，然后返回步骤(3)；如果是，则将其设定为钻井液排量。
所述步骤(1)中的N取5-20。
所述步骤(1)中的钻井基础参数包括：
钻头直径Db，单位为m；钻杆外径Dp，单位为m；钻井液密度ρm，单位为kg/m3；岩屑密度ρd，单位为kg/m3；岩屑等效直径Dd，单位为m；异形井眼井径扩大率、井径扩大处初始岩屑浓度，单位为％。
所述步骤(2)中划分网格时将井眼的轴对称问题简化为2维平面计算问题，并对异形井眼位置处的网格进行加密。
所述步骤(3)中的4种流变模型为：
宾汉流型：τ＝τYP+μPV×γ；
幂律流型：τ＝K×γn；
卡森模式：τ1/2＝τc1/2+η∞ 1/2×γ1/2；
赫-巴流型：τ＝τHB+K×γn
其中，τYP、μPV、K、n、τc、η∞、τHB分别通过数据拟合计算得到，K为稠度系数,单位是Pa·sn，n为流性指数，γ为剪切速率,单位是s-1，η∞是卡森粘度,单位是mPa·s，τ是剪切应力,单位是Pa，τc是卡森屈服值,单位是Pa，τYP是屈服值,单位是Pa，μPV是塑性粘度,单位是mPa·s，τHB是赫巴切力。
所述步骤(5)中的设置钻井液流型的入口和出口的边界条件是这样实现的：
将入口的边界条件取为流入速度边界条件；
将出口的边界条件取为流出速度边界条件；
井壁和钻具外壁处采用无滑移边界条件，设定边界处的流动速度为0。
所述步骤(7)中的流速分布图的绘制区域为步骤(1)中设置的xy二维计算区域，流速为z值，用颜色深度表示；
所述步骤(7)中的岩屑浓度分布图的绘制区域为步骤(1)设置的xy二维计算区域，岩屑浓度为z值，用颜色深度表示。
所述步骤(9)中的合理排量范围为：0～50L/s。
与现有技术相比，本发明的有益效果是：利用本发明，可以克服现有经验法或基于规则井眼的环空水力学法的不足，解决钻井现场仅凭经验或环空水力学法在异形井眼中调整钻井液排量精度不足的问题。本发明可以在满足异形井眼清洁要求的前提下，实现现场钻井液排量和流变参数的合理调控，达到安全高效钻井的目的。
附图说明
图1本发明方法的步骤框图
图2实施例1中的异形井眼计算网格图
图3实施例1中的钻井液卡森流变模型拟合结果图
图4实施例1中的排量10L/s时环空流速分布图
图5实施例1中的排量10L/s时稳定循环状态下的岩屑浓度分布图
图6实施例1中的排量40L/s时流速分布图
图7实施例1中的排量40L/s时稳定循环状态下的岩屑浓度分布图
图8实施例2中的钻井液赫巴流变模型拟合结果图
图9实施例2中的排量10L/s时流速分布图
图10实施例2中的排量10L/s时稳定循环状态下的岩屑浓度分布图
图11实施例2中的40排量40L/s时流速分布图
图12实施例2中的排量40L/s时稳定循环状态下的岩屑浓度分布图
图13实施例3中的异形井眼计算网格图
图14实施例3中的钻井液卡森流变模型拟合结果图
图15实施例3中的10L/s排量下流速分布图
图16实施例3中的10L/s排量下的清除岩屑过程中的岩屑浓度分布图
图17实施例3中的10L/s排量下的稳定循环时的岩屑浓度最终分布图
图18实施例3中的40L/s排量下流速分布图
图19实施例3中的40L/s排量下的清除岩屑过程中的岩屑浓度分布图
图20实施例3中的40L/s排量下的稳定循环时的岩屑浓度最终分布图。
具体实施方式
下面结合附图对本发明作进一步详细描述：
实际钻井中常常遇到井眼扩径使得原来规则的圆筒状井眼编程不规则的异形井眼。本方法根据实际钻井过程中遇到的异形井眼的形状建立井眼几何计算区域，并基于当前使用的钻井液密度、流变参数、岩屑大小和密度参数建立井眼计算流体力学固液两相流模型包括计算区域内的两相流体控制方程、初始条件，通过调整钻井液的排量改变计算区域的流入、流出边界条件，再经过数值计算得到稳定状态下井眼异形位置的岩屑浓度分布及岩屑浓度最大值，根据岩屑浓度最大值是否小于规定的限值判断井眼清洁效果是否满足要求，如果不满足则通过修改排量大小再进行计算，直到计算得到的岩屑浓度最大值小于规定的限值为止，最终得到的钻井液排量值为确定的满足异形井眼清洁的钻井液排量值。
如图1所示，本发明的一种满足异形井眼清洁的钻井液排量确定方法步骤如下：
(1)基于钻井基础参数(包括钻头直径Db(m)、钻杆外径Dp(m)、钻井液密度ρm(kg/m3)、岩屑密度ρd(kg/m3)、岩屑等效直径Dd(m)、异形井眼井径扩大率、井径扩大处初始岩屑浓度(％))，利用现有的Fluent软件建立异形井眼的几何模型，以10倍钻具直径的长度为纵向范围，以井眼扩大的最大值为横向范围，确定计算区域；
(2)根据井眼几何模型划分网格，为提高计算精度并节约计算时间，将井眼的轴对称问题简化为2维平面计算问题，并对异形井眼位置加密网格；
(3)从振动筛前的高架槽内取钻井液样，用旋转粘度计测量钻井液在不同转速Φ下的读数θ，即θ3，θ6，θ100，θ200，θ300，θ600(角标对应旋转粘度计转速)。将旋转粘度计读数转变为相应的一组剪切速率数组γ([5.11,10.22,170.33,511,340.67,1022])和一组剪切应力τ数组([θ3×0.511，θ6×0.511，θ100×0.511，θ200×0.511，θ300×0.511，θ600×0.511])；
(4)根据剪切速率和剪切应力，回归拟合4种流变模式，并计算各模型拟合的相关系数：宾汉流型(τ＝τYP+μPV×γ)、幂律流型(τ＝K×γn)、卡森模式(τ1/2＝τc1/2+η∞ 1/2×γ1/2))和赫-巴流型(τ＝τHB+K×γn)，这4个公式中，除开剪切速率γ和剪切应力τ是已知一组数据，τYP、μPV、K、n、τc、η∞、τHB分别通过数据拟合计算得到。式中：K—稠度系数,Pa·sn，n—流性指数，γ—剪切速率,s-1；；η∞—卡森粘度(接近钻头喷嘴粘度),mPa·s；τ—剪切应力,Pa；τc—卡森屈服值,Pa，τYP-屈服值(动切力),Pa；μPV—塑性粘度,mPa·s；τHB—赫巴切力；根据数值回归拟合的方法计算流变模型的相关系数大小，确定相关系数最大的流变模式为拟合效果最佳的钻井液流型；获得的钻井液粘度参数(τYP、μPV；K、n；τc、η∞；τHB、K、n)用在步骤7中固液两相混合物粘度的计算中；；
(5)设定边界条件：设定钻井液的排量Q(初始可设为额定最大排量，例如50L/s)，设置模型的入口和出口的边界条件。将入口的边界条件取为流入速度边界条件，计算方法为：4Q/π/(Db 2-Dp 2)；出口边界条件设为流出速度边界条件(outflow)，计算方法为：4Q/π/(Db 2-Dp 2)；井壁和钻具外壁处常采用无滑移边界条件，设定边界处的流动速度为0；
(6)设定初始条件：按照循环开始前环空可能出现的最高岩屑浓度设定初始岩屑浓度α2，0，以这个浓度值作为岩屑浓度初始值进行后续迭代计算，每迭代计算一次更新一次环空中的岩屑浓度分布，分析岩屑沉降堆积情况；
(7)基于计算流体力学的固液两相流的数值方法(即下面的连续方程、动量方程和能量方程)(其中，
为两相混合流速，αk从第k相的浓度，α2为岩屑的浓度)，根据数值计算得到的每个网格的
值绘制钻井液异形井眼处的流速分布图(绘制区域为步骤(1)中设置的xy二维计算区域，流速为z值，用颜色深度表示)，根据数值计算中每个网格的固相流动速度差计算得到的每个网格的α2值(α2＝α2，原+v2，流进-v2，流出，其中α2，原是迭代上一步的岩屑浓度，v2，流进是每个网格岩屑的流进速度，v2，流出是每个网格岩屑的流出速度)绘制钻井液异形井眼处的岩屑浓度分布图(绘制区域为步骤(1)设置的xy二维计算区域，岩屑浓度为z值，用颜色深度表示)。
固液两相流数值计算的控制方程如下：
①连续方程：
其中：
ρk为第k相的密度，αk为第k相的密度体积分数，
为第k相的速度；
②动量方程，由固液两相的动量叠加而成：
其中μm是固液两相混合物粘度，
③能量方程：
其中：
Ek为第k相的机械能，hk为第k的重力势能，p为压力，ρk为第k相的密度，vk是第k相的速度；
(8)如果步骤(7)中数值计算得到的每个网格的岩屑浓度值α2值都小于预设的岩屑浓度最大允许值α2，max(例如3％)，则步骤(5)中的排量Q值满足岩屑输送要求，转入步骤(9)，否则返回步骤(5)，重新进行(5)-(7)的计算，直到满足每个网格的岩屑浓度值α2值都小于预设的岩屑浓度最大允许值α2，max为止；
(9)判断满足井眼清洁的排量Q是否在合理排量范围内，所述合理排量范围为0～50L/s；如果不在合理排量范围内，则添加能够改变钻井液的流变性的钻井液助剂，调整钻井液流变性，然后返回步骤(3)；如果排量Q在合理排量范围，则将其设定为钻井排量。
实施例1：
(1)获取钻井基础参数如表1所示：
| 钻头直径D<sub>b</sub>，m | 0.216 |
| 钻杆外径Dp,m | 0.127 |
| 钻井液密度ρ<sub>m,</sub>，kg/m<sup>3</sup> | 1750 |
| 岩屑密度ρ<sub>d</sub>,kg/m<sup>3</sup> | 2700 |
| 岩屑等效直径Dd,m | 0.005 |
| 井径扩大率 | 100％ |
| 井径扩大处初始岩屑浓度 | 10％ |
表1
(2)建立异形井眼的几何模型，划分网格，并在异形井眼处加密网格，如图2所示。
(3)在高架槽处取钻井液样，用旋转粘度计测量一组剪切速率和剪切应力，如表2所示。
表2
根据剪切速率和剪切应力，回归拟合宾汉流型流变模式，并计算各模型的相关系数：(τ＝τYP+μPV×γ)、幂律流型(τ＝K×γn)、卡森模式
和赫-巴流型(τ＝τHB+Kγn)，如表3所示。
| 流变模式类型 | 拟合的公式 | 相关系数 |
| 宾汉流型 | τ＝4.41+0.032*γ | 0.975 |
| 幂律流型 | τ＝1.14*γ^0.47 | 0.965 |
| 卡森模式 | τ<sup>1/2</sup>＝1.91<sup>1/2</sup>+0.02<sup>1/2</sup>*γ<sup>1/2</sup> | 0.994 |
| 赫-巴流型 | τ＝1.90+0.25*γ^0.70 | 0.998 |
表3
(4)确定相关系数最大的赫-巴模型作为计算用的流变类型，如图3所示，并为模型中的流变参数赋值；
(5)根据设定的排量范围(10L/s-40L/s)，设置模型的入口和出口的边界条件。先将入口的边界条件取为10L/s。基于计算流体力学的固液两相流方法计算钻井液流动和岩屑分布，如图4和图5所示，分析得到稳定循环状态下流速图可知在异形井眼的扩径处，钻井液的流动存在涡流，流速相对较低；分析岩屑浓度的分布图可知异形井眼扩径处岩屑浓度较高，岩屑浓度在8％以上，最高岩屑浓度值22.2％，判断容易产生岩屑堆积的位置为井径扩大处的下边沿；岩屑携带效果不理想。
(6)同样在(3)的钻井液流变条件下，将入口的边界条件取为40L/s。分析得到稳定循环状态下异形井眼的扩径处钻井液的流动仍存在涡流，如图6和图7所示，流速相对较低；分析岩屑浓度的分布图可知异形井眼扩径处岩屑浓度较高，岩屑浓度在7％以上，最高岩屑浓度值13.4％，判断容易产生岩屑堆积的位置为井径扩大处的下边沿；岩屑携带效果不理想。因为不存在满足井眼清洁要求的排量范围，因此需要调整钻井液的流变性能。
实施例2
钻井基础参数和异形井眼形状与实施例1相同，在实施例1中的钻井液流变条件下，不存在满足井眼清洁要求的排量范围，在实施例2中通过添加钻井液助剂，调整钻井液流变性。调整后的钻井液流变性如表4所示.
| 旋转粘度计转速，r/min | 600 | 300 | 200 | 100 | 6 | 3 |
| 剪切速率，/s | 1022 | 511 | 340.7 | 170.3 | 10.2 | 5.11 |
| 旋转粘度计读数θ | 172 | 128 | 99 | 71 | 36 | 27 |
| 剪切应力τ，Pa | 87.9 | 65.4 | 50.6 | 36.3 | 18.4 | 13.8 |
表4
根据剪切速率和剪切应力，回归拟合4种流变模式，确定相关系数最大的赫-巴模型作为计算用的流变类型，如表5和图8所示。
| 流变模式类型 | 拟合的公式 | 相关系数 |
| 宾汉流型 | τ＝20.74+0.072*γ | 0.949 |
| 幂律流型 | τ＝8.01*γ^0.32 | 0.956 |
| 卡森模式 | τ<sup>1/2</sup>＝12.55<sup>1/2</sup>+0.036<sup>1/2</sup>*γ<sup>1/2</sup> | 0.989 |
| 赫-巴流型 | τ＝10.35+1.73*γ^0.547 | 0.994 |
表5
(2)将钻井液流变条件下将入口的边界条件取为10L/s。分析得到稳定循环状态下异形井眼的扩径处钻井液的流动为向上的层流，如图9和图10所示；分析岩屑浓度的分布图可知异形井眼扩径处大部分区域岩屑浓度在1％左右，最高岩屑浓度值2.45％，判断岩屑携带效果理想。
(3)将钻井液流变条件下将入口的边界条件取为40L/s。分析得到稳定循环状态下异形井眼的扩径处钻井液的流动为涡流状态，如图11和图12所示；分析岩屑浓度的分布图可知异形井眼扩径处大部分区域岩屑浓度在3％左右，最高岩屑浓度值4.09％，判断岩屑携带效果不如相同流变条件下的10L/s排量。
(4)正常维护钻井液，保持(1)中稳定的流变参数，并优选10L/s的排量钻进。
实施例3：
(1)获取钻井基础参数如下：
| 钻头直径D<sub>b</sub>，m | 0.216 |
| 钻杆外径Dp，m | 0.127 |
| 钻井液密度ρ<sub>m，</sub>，kg/m<sup>3</sup> | 1800 |
| 岩屑密度ρ<sub>d</sub>，kg/m<sup>3</sup> | 2700 |
| 岩屑等效直径Dd，m | 0.005 |
| 井径扩大率 | 79％ |
| 井径扩大处初始岩屑浓度 | 10％ |
表6
(2)建立异形井眼的几何模型，划分网格，并在异形井眼处加密网格，如图13所示。
(3)在高架槽处取钻井液样，用旋转粘度计测量一组剪切速率和剪切应力，如表7所示。
表7
根据剪切速率和剪切应力，回归拟合4种流变模式，并计算各模型的相关系数：宾汉流型(τ＝τYP+μPV×γ)、幂律流型(τ＝K×γn)、卡森模式
和赫-巴流型(τ＝τHB+Kγn)，如表8。
| 流变模式类型 | 拟合的公式 | 相关系数 |
| 宾汉流型 | τ＝20.687+0.0767*γ | 0.951 |
| 幂律流型 | τ＝7.767*γ^0.339 | 0.957 |
| 卡森模式 | τ<sup>1/2</sup>＝12.18<sup>1/2</sup>+0.039<sup>1/2</sup>*γ<sup>1/2</sup> | 0.990 |
| 赫-巴流型 | τ＝10.73+1.534*γ^0.571 | 0.996 |
表8
(4)确定相关系数最大的赫-巴模型作为计算用的流变类型，如图14，并为模型中的流变参数赋值；
(5)根据设定的排量范围(10L/s-40L/s)，设置模型的入口和出口的边界条件。先将入口的边界条件取为10L/s。基于计算流体力学的固液两相流方法计算钻井液流动和岩屑分布，分析得到稳定循环状态下流速图可知在异形井眼的扩径处，钻井液的流动为层流，流速相对较低；分析岩屑浓度的分布图可知异形井眼扩径处岩屑浓度在1％左右，与稠浆中的岩屑浓度相同，表明异形井眼处的岩屑被完全携带出，因此岩屑携带效果理想，如图15、图16和图17所示，。
(6)将入口的边界条件取为40L/s。分析得到稳定循环状态下异形井眼的扩径处钻井液的流动存在涡流，流速相对较低；分析岩屑浓度的分布图可知异形井眼扩径处岩屑浓度较高，岩屑浓度在5％以上，最高岩屑浓度值9.15％，判断容易产生岩屑堆积的位置为井径扩大处的中部和下边沿；岩屑携带效果不理想。因此40L/s的排量效果不如10L/s的排量，因此选择10L/s的排量作为携带岩屑的排量，如图18、图19和图20所示。
本发明主要用于石油天然气勘探开发领域，在复杂井眼条件下满足井眼清洁的要求对钻井液和钻井参数进行优化。在钻井中，井眼清洁效果直接关系到井眼净化、钻井速度及井下安全等多个方面，是十分重要的钻井参数。在实际钻井过程中，井眼形状常常出现扩径的异形现象。在异形井眼的位置附近，钻井液的流动状态和岩屑携带情况变得更加复杂，岩屑不容易被携带出井眼，容易产生岩屑堆积的现象，可能造成下钻困难、下钻不到位以及卡钻等井下复杂情况。因此异形井眼内的井眼清洁效果成为能否安全、快速钻进的重要问题。利用本发明，可以解决钻井时计算异形井眼中井眼清洁效果的问题。本发明可以克服现有经验法或环空流体力学法的不足，解决钻井现场仅凭经验或基于规则井眼的调整钻井液排量和流变参数的问题。本发明可以在满足异形井眼清洁要求的前提下，实现现场钻井液排量和流变参数的合理调控，达到安全高效钻井的目的。本发明可用于石油勘探开发领域，但不受该实施范围的限制，亦可用于地质勘探、水文水井、地热等领域的井眼清洁效果评价。
上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述具体实施方式所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。
Claims (4)
1.一种满足异形井眼清洁的钻井液排量确定方法，其特征在于：所述方法首先基于钻井基础参数建立异形井眼几何模型，设定初始的排量，然后通过数值计算得到稳定状态下井眼各个位置的岩屑浓度分布，通过岩屑浓度分布判断是否满足井眼清洁效果，并判断钻井液排量是否在合理排量范围内，如果有一个为否，则调整钻井液流变性，然后重复上述过程，直到钻井液排量满足井眼清洁效果且在合理排量范围内，此时的钻井液排量即为最终确定的钻井液排量；
所述方法包括：
(1)基于钻井基础参数，建立异形井眼几何模型，以N倍钻具直径的长度为纵向范围，以井眼扩大的最大值为横向范围，确定xy二维计算区域；
(2)根据异形井眼几何模型划分网格；
(3)从振动筛前的高架槽内取钻井液，用旋转粘度计测量钻井液在不同转速Φ下的读数θ，将旋转粘度计的读数转变为相应的一组剪切速率γ和剪切应力τ；根据剪切速率和剪切应力，回归拟合4种流变模型；
(4)，计算各流变模型拟合的相关系数，确定相关系数最大的流变模型为拟合效果最佳的钻井液流变流型；
(5)设定排量Q，设置钻井液流型的入口和出口的边界条件；
(6)设定初始条件：按照循环开始前环空可能出现的最高岩屑浓度设定岩屑浓度α2的初始值；
(7)基于连续方程、动量方程和能量方程计算得到的每个网格的两相混合流速
的值绘制钻井液异形井眼处的流速分布图，根据每个网格的固相流动速度差计算得到的每个网格的岩屑浓度α2的值绘制钻井液异形井眼处的岩屑浓度分布图；其中，α2＝α2，原+v2，流进-v2，流出，其中α2，原是迭代上一步的岩屑浓度，v2，流进是每个网格岩屑的流进速度，v2，流出是每个网格岩屑的流出速度；
(8)如果步骤(7)中得到的每个网格的岩屑浓度值α2值都小于预设的岩屑浓度最大允许值α2，max，则判定步骤(5)中的排量Q值满足井眼清洁要求，转入步骤(9)，否则返回步骤(5)；
(9)判断满足井眼清洁要求的排量Q是否在合理排量范围内，如果否，则添加能够改变钻井液的流变性的钻井液助剂，调整钻井液流变性，然后返回步骤(3)；如果是，则将其设定为钻井液排量；
所述步骤(1)中的N取5-20；
所述步骤(2)中划分网格时将井眼的轴对称问题简化为2维平面计算问题，并对异形井眼位置处的网格进行加密；
所述步骤(3)中的4种流变模型为：
宾汉流型：τ＝τYP+μPV×γ；
幂律流型：τ＝K×γn；
卡森模式：τ1/2＝τc1/2+η∞ 1/2×γ1/2；
赫-巴流型：τ＝τHB+K×γn
其中，τYP、μPV、K、n、τc、η∞、τHB分别通过数据拟合计算得到，K为稠度系数,单位是Pa·sn，n为流性指数，γ为剪切速率,单位是s-1，η∞是卡森粘度,单位是mPa·s，τ是剪切应力,单位是Pa，τc是卡森屈服值,单位是Pa，τYP是屈服值,单位是Pa，μPV是塑性粘度,单位是mPa·s，τHB是赫巴切力；
所述步骤(5)中的设置钻井液流型的入口和出口的边界条件是这样实现的：
将入口的边界条件取为流入速度边界条件；
将出口的边界条件取为流出速度边界条件；
井壁和钻具外壁处采用无滑移边界条件，设定边界处的流动速度为0。
2.根据权利要求1所述的满足异形井眼清洁的钻井液排量确定方法，其特征在于：所述步骤(1)中的钻井基础参数包括：
钻头直径Db，单位为m；钻杆外径Dp，单位为m；钻井液密度ρm，单位为kg/m3；岩屑密度ρd，单位为kg/m3；岩屑等效直径Dd，单位为m；异形井眼井径扩大率、井径扩大处初始岩屑浓度，单位为％。
3.根据权利要求2所述的满足异形井眼清洁的钻井液排量确定方法，其特征在于：所述步骤(7)中的流速分布图的绘制区域为步骤(1)中设置的xy二维计算区域，流速为z值，用颜色深度表示；
所述步骤(7)中的岩屑浓度分布图的绘制区域为步骤(1)设置的xy二维计算区域，岩屑浓度为z值，用颜色深度表示。
4.根据权利要求3所述的满足异形井眼清洁的钻井液排量确定方法，其特征在于：所述步骤(9)中的合理排量范围为：0～50L/s。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711170726.XA CN109812237B (zh) | 2017-11-22 | 2017-11-22 | 一种满足异形井眼清洁的钻井液排量确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711170726.XA CN109812237B (zh) | 2017-11-22 | 2017-11-22 | 一种满足异形井眼清洁的钻井液排量确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109812237A CN109812237A (zh) | 2019-05-28 |
| CN109812237B true CN109812237B (zh) | 2021-01-29 |
Family
ID=66601227
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711170726.XA Active CN109812237B (zh) | 2017-11-22 | 2017-11-22 | 一种满足异形井眼清洁的钻井液排量确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109812237B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110230479A (zh) * | 2019-06-12 | 2019-09-13 | 中国海洋石油集团有限公司 | 一种大斜度井油基钻井液井眼清洁的控制方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| BR112012007730A2 (pt) * | 2009-10-06 | 2016-08-23 | Prad Res & Dev Ltd | planejamento e monitoramento de testes de formação |
| CN203025064U (zh) * | 2012-12-18 | 2013-06-26 | 中国石油化工股份有限公司 | 钻井液性能测试装置 |
| CN106469228B (zh) * | 2015-08-14 | 2019-09-20 | 中国石油化工股份有限公司 | 基于在线测试与地质参数的热采井吸汽剖面解释方法 |
| CN105257242B (zh) * | 2015-11-20 | 2017-07-25 | 东北石油大学 | 偏心环空注水泥附加量确定方法 |
| CN105952400B (zh) * | 2016-04-26 | 2018-07-06 | 西南石油大学 | 一种环空井眼清洁实时监测方法 |
-
2017
- 2017-11-22 CN CN201711170726.XA patent/CN109812237B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109812237A (zh) | 2019-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lichun et al. | Drilling fluid loss model and loss dynamic behavior in fractured formations | |
| CN101684727B (zh) | 超深井稠油掺稀比例确定的优化方法及其掺稀混配器 | |
| CN106437609B (zh) | 一种高温高压超深井全过程塞流防漏固井设计方法 | |
| Szanyi et al. | Near-wellbore modeling of a horizontal well with computational fluid dynamics | |
| CN103590824A (zh) | 经过多段压裂改造后的致密气藏水平井的产能计算方法 | |
| CN109812236B (zh) | 一种确定异形井眼中的井眼清洁效果的方法 | |
| Fisher et al. | Numerical modelling of cake formation and fluid loss from non-Newtonian muds during drilling using eccentric/concentric drill strings with/without rotation | |
| Ferroudji et al. | 3D numerical and experimental modelling of multiphase flow through an annular geometry applied for cuttings transport | |
| CN109812237B (zh) | 一种满足异形井眼清洁的钻井液排量确定方法 | |
| Liu et al. | A new solution to enhance cuttings transport in mining drilling by using pulse jet mill technique | |
| CN108756830B (zh) | 一种水平井砾石充填模拟计算方法 | |
| Noah | Optimizing drilling fluid properties and flow rates for effective hole cleaning at high-angle and horizontal wells | |
| Kragset et al. | Effect of buoyancy and inertia on viscoplastic fluid: Fluid displacement in an inclined eccentric annulus with an irregular section | |
| Egbue | Improving the Efficiency of Transportation of Cuttings in Wellbore-Experimental Investigations on Critical Rolling and Lifting Velocities | |
| CN110593856B (zh) | 一种固井安全作业密度窗口测定方法 | |
| Roustaei | Yield stress fluid flows in uneven geometries: applications to the oil & gas industry | |
| CN109815511B (zh) | 一种满足异形井眼清洁的钻井液密度确定方法 | |
| Dutra et al. | Analysis of interface between Newtonian and non-Newtonian fluids inside annular eccentric tubes | |
| Hussain et al. | Analysis of yield-power-law fluid flow in irregular eccentric annuli | |
| Ali | A parametric study of cutting transport in vertical and horizontal well using computational fluid dynamics (CFD) | |
| Zulqarnain et al. | The effect of constant and variable eccentricity on the spacer performance during primary well cementing | |
| Ferroudji et al. | Effects of the inner pipe rotation and rheological parameters on the axial and tangential velocity profiles and pressure drop of yield Power-Law fluid in eccentric annulus | |
| CN114077812B (zh) | 一种用于设计钻井液流变参数的方法 | |
| Wang et al. | Numerical simulation of jet excavation in conductor jetting operations | |
| Braghini et al. | Effect of rheology on flow displacement during cementing process in oil wells |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |