CN109812237B - 一种满足异形井眼清洁的钻井液排量确定方法 - Google Patents

Abstract

本发明提供了一种满足异形井眼清洁的钻井液排量确定方法，属于石油天然气勘探开发钻井领域。该方法所述方法首先基于钻井基础参数建立异形井眼几何模型，设定初始的排量，然后通过数值计算得到稳定状态下井眼各个位置的岩屑浓度分布，通过岩屑浓度分布判断是否满足井眼清洁效果，并判断钻井液排量是否在合理排量范围内，如果有一个为否，则调整钻井液流变性，然后重复上述过程，直到钻井液排量满足井眼清洁效果且在合理排量范围内，此时的钻井液排量即为最终确定的钻井液排量。利用本发明可以克服现有经验法或基于规则井眼的环空水力学法的不足，解决钻井现场仅凭经验或环空水力学法在异形井眼中调整钻井液排量精度不足的问题。

Description

一种满足异形井眼清洁的钻井液排量确定方法

技术领域

本发明属于石油天然气勘探开发钻井领域，具体涉及一种满足异形井眼清洁的钻井液排量确定方法。

背景技术

钻井液是钻井过程中的循环携带岩屑的流体介质。井眼清洁程度是指钻井时井筒中岩屑携带的状态。钻井中的岩屑携带状况直接关系到钻井速度及井下安全等多个方面，是十分重要的钻井参数。钻井过程中，如果岩屑不能及时携带出井眼，可能造成下钻困难、下钻不到位以及卡钻等井下复杂情况。因此井眼清洁成为能否安全、快速钻进的首要问题。

实践表明，井眼的形状和钻井液的排量、流变参数对钻井中的岩屑携带有着显著的影响。在实际钻井过程中，井眼常常出现扩径的异形现象，在异形井眼的扩径位置附近，钻井液的流动状态和岩屑携带情况变得更加复杂，容易产生岩屑堆积的现象。

根据异形井眼中井眼清洁指标的要求，可以对钻井液排量和性能调整做出合理判断，利用现场小型试验，优化钻井液流变参数，调整钻井液排量，可以改善岩屑携带效果，有效增强钻井安全并降低钻井成本。

目前，现场采用的钻井液排量选取方法主要有：(1)仅凭经验调整排量和钻井液流变参数；(2)在井眼规则的假设下，基于环空水力学法计算不同钻井液排量和流变参数下的岩屑携带效果，根据岩屑携带效果是否在合理范围确定是否调整钻井液排量和流变参数。

目前的方法存在着一些不足，主要体现在：(1)仅凭经验调整排量和钻井液流变参数盲目性较大；(2)假定的规则井眼形状可能与实际井眼形状偏差较大，未采用最符合实际井眼形状计算岩屑携带效果，计算精度低，无法准确判断岩屑在井眼中的堆积位置。以上的不足可能导致钻井液排量和流变性能维护调控不当，浪费钻井液材料，增加钻井成本，极易引起井下复杂情况的发生。

发明内容

本发明的目的在于解决上述现有技术中存在的难题，提供一种满足异形井眼清洁的钻井液排量确定方法，解决钻井现场仅凭经验或基于规则井眼的调整钻井液排量的问题。在满足异形井眼清洁要求的前提下，实现现场钻井液排量的合理调控，确定钻井液维护处理措施，满足安全高效钻井的要求并降低钻井成本。

本发明是通过以下技术方案实现的：

一种满足异形井眼清洁的钻井液排量确定方法，所述方法首先基于钻井基础参数建立异形井眼几何模型，设定初始的排量，然后通过数值计算得到稳定状态下井眼各个位置的岩屑浓度分布，通过岩屑浓度分布判断是否满足井眼清洁效果，并判断钻井液排量是否在合理排量范围内，如果有一个为否，则调整钻井液流变性，然后重复上述过程，直到钻井液排量满足井眼清洁效果且在合理排量范围内，此时的钻井液排量即为最终确定的钻井液排量。

所述方法包括：

(1)基于钻井基础参数，建立异形井眼几何模型，以N倍钻具直径的长度为纵向范围，以井眼扩大的最大值为横向范围，确定xy二维计算区域；

(2)根据异形井眼几何模型划分网格；

(3)从振动筛前的高架槽内取钻井液，用旋转粘度计测量钻井液在不同转速Φ下的读数θ，将旋转粘度计的读数转变为相应的一组剪切速率γ和剪切应力τ；根据剪切速率和剪切应力，回归拟合4种流变模型；

(4)，计算各流变模型拟合的相关系数，确定相关系数最大的流变模型为拟合效果最佳的钻井液流变流型；

(5)设定排量Q，设置钻井液流型的入口和出口的边界条件；

(6)设定初始条件：按照循环开始前环空可能出现的最高岩屑浓度设定岩屑浓度α₂的初始值；

(7)基于连续方程、动量方程和能量方程计算得到的每个网格的两相混合流速

的值绘制钻井液异形井眼处的流速分布图，根据每个网格的固相流动速度差计算得到的每个网格的岩屑浓度α₂的值绘制钻井液异形井眼处的岩屑浓度分布图；其中，α₂＝α_2，原+v_2，流进-v_2，流出，其中α_2，原是迭代上一步的岩屑浓度，v_2，流进是每个网格岩屑的流进速度，v_2，流出是每个网格岩屑的流出速度；

(8)如果步骤(7)中得到的每个网格的岩屑浓度值α₂值都小于预设的岩屑浓度最大允许值α_2，max，则判定步骤(5)中的排量Q值满足井眼清洁要求，转入步骤(9)，否则返回步骤(5)；

(9)判断满足井眼清洁要求的排量Q是否在合理排量范围内，如果否，则添加能够改变钻井液的流变性的钻井液助剂，调整钻井液流变性，然后返回步骤(3)；如果是，则将其设定为钻井液排量。

所述步骤(1)中的N取5-20。

所述步骤(1)中的钻井基础参数包括：

钻头直径D_b，单位为m；钻杆外径D_p，单位为m；钻井液密度ρ_m，单位为kg/m³；岩屑密度ρ_d，单位为kg/m³；岩屑等效直径Dd，单位为m；异形井眼井径扩大率、井径扩大处初始岩屑浓度，单位为％。

所述步骤(2)中划分网格时将井眼的轴对称问题简化为2维平面计算问题，并对异形井眼位置处的网格进行加密。

所述步骤(3)中的4种流变模型为：

宾汉流型：τ＝τ_YP+μ_PV×γ；

幂律流型：τ＝K×γⁿ；

卡森模式：τ^1/2＝τc^1/2+η_∞ ^1/2×γ^1/2；

赫-巴流型：τ＝τ_HB+K×γⁿ

其中，τ_YP、μ_PV、K、n、τc、η_∞、τ_HB分别通过数据拟合计算得到，K为稠度系数,单位是Pa·sⁿ，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)基于钻井基础参数(包括钻头直径D_b(m)、钻杆外径D_p(m)、钻井液密度ρ_m(kg/m³)、岩屑密度ρ_d(kg/m³)、岩屑等效直径Dd(m)、异形井眼井径扩大率、井径扩大处初始岩屑浓度(％))，利用现有的Fluent软件建立异形井眼的几何模型，以10倍钻具直径的长度为纵向范围，以井眼扩大的最大值为横向范围，确定计算区域；

(2)根据井眼几何模型划分网格，为提高计算精度并节约计算时间，将井眼的轴对称问题简化为2维平面计算问题，并对异形井眼位置加密网格；

(3)从振动筛前的高架槽内取钻井液样，用旋转粘度计测量钻井液在不同转速Φ下的读数θ，即θ₃，θ₆，θ₁₀₀，θ₂₀₀，θ₃₀₀，θ₆₀₀(角标对应旋转粘度计转速)。将旋转粘度计读数转变为相应的一组剪切速率数组γ([5.11,10.22,170.33,511,340.67,1022])和一组剪切应力τ数组([θ₃×0.511，θ₆×0.511，θ₁₀₀×0.511，θ₂₀₀×0.511，θ₃₀₀×0.511，θ₆₀₀×0.511])；

(4)根据剪切速率和剪切应力，回归拟合4种流变模式，并计算各模型拟合的相关系数：宾汉流型(τ＝τ_YP+μ_PV×γ)、幂律流型(τ＝K×γⁿ)、卡森模式(τ^1/2＝τc^1/2+η_∞ ^1/2×γ^1/2))和赫-巴流型(τ＝τ_HB+K×γⁿ)，这4个公式中，除开剪切速率γ和剪切应力τ是已知一组数据，τ_YP、μ_PV、K、n、τc、η_∞、τ_HB分别通过数据拟合计算得到。式中：K—稠度系数,Pa·sⁿ，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/π/(D_b ²-D_p ²)；出口边界条件设为流出速度边界条件(outflow)，计算方法为：4Q/π/(D_b ²-D_p ²)；井壁和钻具外壁处常采用无滑移边界条件，设定边界处的流动速度为0；

(6)设定初始条件：按照循环开始前环空可能出现的最高岩屑浓度设定初始岩屑浓度α_2，0，以这个浓度值作为岩屑浓度初始值进行后续迭代计算，每迭代计算一次更新一次环空中的岩屑浓度分布，分析岩屑沉降堆积情况；

(7)基于计算流体力学的固液两相流的数值方法(即下面的连续方程、动量方程和能量方程)(其中，

为两相混合流速，α_k从第k相的浓度，α₂为岩屑的浓度)，根据数值计算得到的每个网格的

值绘制钻井液异形井眼处的流速分布图(绘制区域为步骤(1)中设置的xy二维计算区域，流速为z值，用颜色深度表示)，根据数值计算中每个网格的固相流动速度差计算得到的每个网格的α₂值(α₂＝α_2，原+v_2，流进-v_2，流出，其中α_2，原是迭代上一步的岩屑浓度，v_2，流进是每个网格岩屑的流进速度，v_2，流出是每个网格岩屑的流出速度)绘制钻井液异形井眼处的岩屑浓度分布图(绘制区域为步骤(1)设置的xy二维计算区域，岩屑浓度为z值，用颜色深度表示)。

固液两相流数值计算的控制方程如下：

①连续方程：

其中：

ρ_k为第k相的密度，α_k为第k相的密度体积分数，

为第k相的速度；

②动量方程，由固液两相的动量叠加而成：

其中μ_m是固液两相混合物粘度，

是两相流中的第k相(固相)的漂移速度：

③能量方程：

其中：

E_k为第k相的机械能，h_k为第k的重力势能，p为压力，ρ_k为第k相的密度，v_k是第k相的速度；

(8)如果步骤(7)中数值计算得到的每个网格的岩屑浓度值α₂值都小于预设的岩屑浓度最大允许值α_2，max(例如3％)，则步骤(5)中的排量Q值满足岩屑输送要求，转入步骤(9)，否则返回步骤(5)，重新进行(5)-(7)的计算，直到满足每个网格的岩屑浓度值α₂值都小于预设的岩屑浓度最大允许值α_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×γⁿ)、卡森模式

和赫-巴流型(τ＝τ_HB+Kγⁿ)，如表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×γⁿ)、卡森模式

和赫-巴流型(τ＝τ_HB+Kγⁿ)，如表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)设定初始条件：按照循环开始前环空可能出现的最高岩屑浓度设定岩屑浓度α₂的初始值；

(7)基于连续方程、动量方程和能量方程计算得到的每个网格的两相混合流速

的值绘制钻井液异形井眼处的流速分布图，根据每个网格的固相流动速度差计算得到的每个网格的岩屑浓度α₂的值绘制钻井液异形井眼处的岩屑浓度分布图；其中，α₂＝α_2，原+v_2，流进-v_2，流出，其中α_2，原是迭代上一步的岩屑浓度，v_2，流进是每个网格岩屑的流进速度，v_2，流出是每个网格岩屑的流出速度；

(8)如果步骤(7)中得到的每个网格的岩屑浓度值α₂值都小于预设的岩屑浓度最大允许值α_2，max，则判定步骤(5)中的排量Q值满足井眼清洁要求，转入步骤(9)，否则返回步骤(5)；

(9)判断满足井眼清洁要求的排量Q是否在合理排量范围内，如果否，则添加能够改变钻井液的流变性的钻井液助剂，调整钻井液流变性，然后返回步骤(3)；如果是，则将其设定为钻井液排量；

所述步骤(1)中的N取5-20；

所述步骤(2)中划分网格时将井眼的轴对称问题简化为2维平面计算问题，并对异形井眼位置处的网格进行加密；

所述步骤(3)中的4种流变模型为：

宾汉流型：τ＝τ_YP+μ_PV×γ；

幂律流型：τ＝K×γⁿ；

卡森模式：τ^1/2＝τc^1/2+η_∞ ^1/2×γ^1/2；

赫-巴流型：τ＝τ_HB+K×γⁿ

其中，τ_YP、μ_PV、K、n、τc、η_∞、τ_HB分别通过数据拟合计算得到，K为稠度系数,单位是Pa·sⁿ，n为流性指数，γ为剪切速率,单位是s^-1，η_∞是卡森粘度,单位是mPa·s，τ是剪切应力,单位是Pa，τc是卡森屈服值,单位是Pa，τ_YP是屈服值,单位是Pa，μ_PV是塑性粘度,单位是mPa·s，τ_HB是赫巴切力；

所述步骤(5)中的设置钻井液流型的入口和出口的边界条件是这样实现的：

将入口的边界条件取为流入速度边界条件；

将出口的边界条件取为流出速度边界条件；

井壁和钻具外壁处采用无滑移边界条件，设定边界处的流动速度为0。

2.根据权利要求1所述的满足异形井眼清洁的钻井液排量确定方法，其特征在于：所述步骤(1)中的钻井基础参数包括：

钻头直径D_b，单位为m；钻杆外径D_p，单位为m；钻井液密度ρ_m，单位为kg/m³；岩屑密度ρ_d，单位为kg/m³；岩屑等效直径Dd，单位为m；异形井眼井径扩大率、井径扩大处初始岩屑浓度，单位为％。

3.根据权利要求2所述的满足异形井眼清洁的钻井液排量确定方法，其特征在于：所述步骤(7)中的流速分布图的绘制区域为步骤(1)中设置的xy二维计算区域，流速为z值，用颜色深度表示；

所述步骤(7)中的岩屑浓度分布图的绘制区域为步骤(1)设置的xy二维计算区域，岩屑浓度为z值，用颜色深度表示。

4.根据权利要求3所述的满足异形井眼清洁的钻井液排量确定方法，其特征在于：所述步骤(9)中的合理排量范围为：0～50L/s。

Images