环境敏感区废弃钻井液无害化处理新技术

王勇; 鄢捷年; 吴江; 李志勇; 徐生江; 韦火云

doi:10.3969/j.issn.1001-0890.2014.02.013

环境敏感区废弃钻井液无害化处理新技术

1.
石油工程教育部重点实验室（中国石油大学（北京））,北京 102249;
2.
中国石油集团渤海钻探工程有限公司,天津 300457;
3.
中海石油（中国）有限公司湛江分公司,广东湛江 524057

基金项目:

国家自然科学“复杂油气井钻井与完井基础研究”（编号：51221003）和“低渗裂缝性砂岩凝析气藏损害机理及有效保护开发一体化新方法研究”（编号：51374225）、北京市高等学校青年英才计划项目“全过程欠平衡钻井耐高温高强度冻胶阀研究”（编号：YETP0671）、中国石油大学（北京）科研“全过程欠平衡钻井耐高温高强度冻胶阀研究”（编号：KYJJ2012-02-15）资助。

详细信息

作者简介:
王勇（1971—），男，河北南宫人，1993年毕业于江汉石油学院钻井工程专业，在读博士研究生，高级工程师，主要从事钻井技术研究及相关管理工作。

中图分类号: TE254
计量
- 文章访问数: 2814
- HTML全文浏览量: 96
- PDF下载量: 3807
出版历程
- 收稿日期: 2013-05-26
- 修回日期: 2014-01-11
- 刊出日期: 1899-12-31

New Technology of Wasted Drilling Fluid Treatment in Environmentally Sensitive Areas

1.
MOE Key Laboratory of Petroleum Engineering(China University of Petroleum(Beijing)), Beijing, 102249, China;
2.
CNPC Bohai Drilling Engineering Company Ltd., Tianjin, 300457, China;
3.
Zhanjiang Branch, CNOOC, Zhanjiang, Guangdong, 524057, China

摘要

摘要: 针对环境敏感区的废弃钻井液无害化处理难题，在确定脱稳剂 SD-131、促凝剂 SD-132、分散剂 SD-133等废弃钻井液无害化处理剂加量的基础上，采用单因素试验方法及原子分光光度计、红外测油仪等测试手段，以处理后固体物浸出液中石油类和重金属含量为评价指标，研发了环境敏感区废弃钻井液处理新技术。废弃钻井液经无害化处理后，浸出液色度和悬浮物含量低，未出现返浆现象；同时，浸出液的化学需氧量为23~56 mg/L，总铬含量为0.86~1.12 mg/L，Cr6+含量为0.22~0.40 mg/L，石油类含量为2.67~4.23 mg/L，达到了废弃钻井液无害化处理的目的。在温米油田环境敏感区进行的现场试验表明，废弃钻井液无害化处理后，自然干化23~30 d，仍具有一定的孔隙度和类似土壤的特征，可承载推土机进行井场恢复。现场试验结果证明，该技术可实现对废弃钻井液中污染物的有效转化与封固。
- 环境敏感区 /
- 废钻井液 /
- 废物处理 /
- 污染物分析 /
- 温米油田
Abstract: In view of the disposal of wasted drilling fluid in environmentally sensitive areas,a wasted drilling fluid innocuous treatment formula has been developed on the basis of determining the dosage of destablizer SD-131,coagulation accelerator SD-132,and dispersant SD-133.After innocuous treatment,the infiltrate fluid of drilling fluid has low chromaticity and low content of suspended materials.Different materials and elements in the infiltrate fluid are indicated by COD:23-56 mg/L,total Cr:0.86-1.12 mg/L,Cr6+ 0.22-0.40 mg/L,and hydrocarbon:2.67-4.23 mg/L,which meet the requirement of innocuous treatment of wasted drilling fluid.The test of the technology in environmentally sensitive area of Wenmi Oilfield shows the drilling fluid still has certain porosity and is similar to soil 20-30 d of drying after treatment,and can bear bulldozer to restore well site.Field test shows the technology can effectively transform and solidify pollutants in drilling fluid.
- environmental sensitive area /
- waste drilling fluid /
- waste disposal /
- pollutant analysis /
- Wenmi Oilfield

HTML全文

斜井、大位移井和水平井等复杂结构井钻井过程中，钻柱在井眼中会出现偏心旋转，影响了环空钻井液的流动和循环压耗，进而使井底钻井液当量循环密度（ECD）发生变化，而准确预测环空摩阻压降是复杂结构井控压钻井的重要理论基础^[1-2]。钻井液环空摩阻压降常规计算方法假设钻柱在井眼内居中，且未考虑钻柱旋转，因而无法直接计算复杂结构井的环空摩阻压降。为此，国内外学者开展了偏心环空摩阻压降计算理论和试验研究^[3-25]，建立了各种偏心环空流动摩阻压降计算模型，但均存在一定局限性，如解析模型中只考虑了偏心，没有考虑钻柱旋转，尚未建立符合现场工况的偏心环空摩阻压降计算模型。此外，环空层流摩阻压降的数值本身并不高，且偏心和旋转对环空摩阻压降的影响一般小于20%，因此，试验研究对压力计精度和操作精度的要求更加苛刻，导致文献公布的试验数据不系统，规律认识不一致。所以，基于试验数据建立的回归模型适用性差，现场应用时存在计算结果不合理现象。为此，笔者应用数值模拟方法，研究了正常钻井参数条件下，转速和偏心对复杂结构井典型环空（ϕ127.0 mm钻杆和ϕ215.9 mm井眼）摩阻压降的影响，进而建立了偏心旋转环空无因次摩阻压降计算模型，可以精细描述环空压力场和准确计算ECD，为控压钻井水力参数优化提供指导。

1. 三维偏心环空流动数值模拟

1.1 数值模拟参数

复杂结构井的生产井段井径大多为215.9 mm，通常为控压钻井的主要井段。该井段钻进时，钻柱转速一般不超过120 r/min，钻井液排量不超过40 L/s。ϕ127.0 mm钻杆一般应用ϕ168.0 mm接头，考虑井眼10%的扩径率，钻杆接头紧贴井壁情况下的钻柱偏心度为62.89%，偏心度计算公式可表示为：

$E=\frac{\delta }{{r}_{\mathrm{w}}-{r}_{\mathrm{d}}}$

(1)

式中：E为偏心度； $\delta$ 为偏心距，即沿井筒轴线视图中钻杆和井筒2个圆心间的距离，m； ${r}_{\mathrm{w}}$ 为井筒半径，m； ${r}_{\mathrm{d}}$ 为钻杆半径，m。

数值模拟时，假设钻井液性能和排量不变，改变钻柱转速和偏心度，主要模拟参数设定为：钻井液流性指数n为0.615，稠度系数K为0.518 Pa·sⁿ，钻井液密度1 300 kg/m³、入口速度1.50 m/s，钻柱偏心度0～67.42%，转速0～114.65 r/min。由于偏心度和转速模拟矩阵划分得越细密，建立的回归模型预测结果精度越高，但模型也越复杂。综合考虑这2方面因素和偏心度对摩阻压降梯度的影响规律，不均匀设置了12个偏心度：0（偏心距0），22.47%（偏心距10 mm），26.97%（偏心距12 mm），31.46%（偏心距14 mm），35.96%（偏心距16 mm），40.45%（偏心距18 mm），44.94%（偏心距20 mm），49.44%（偏心距22 mm），53.93%（偏心距24 mm），58.43%（偏心距26 mm），62.92%（偏心距28 mm）和67.42%（偏心距30 mm）；均匀设置了7个转速：0，19.11，38.22，57.32，76.43，95.54和114.65 r/min，构建了84个数值模拟点。

1.2 物理模型及网格划分

为消除入口段对模拟结果的影响，设定环空长度为5 m（大于20倍环空直径），重力方向为X轴，环空轴线方向为Z轴，应用Design Modeler建立了12个偏心度的环空流道几何模型，并进行了结构化网格划分以加速收敛。以22.47%偏心度模型为例，网格划分结果如图1所示，共划分了640 000个控制体网格，673 680个节点。网格最大扭曲度为0.28，最小扭曲度为1.54×10⁻²，平均扭曲度为0.13，网格质量优良，满足FLUENT对网格质量的要求，有利于计算收敛。

图 1 22.47%偏心度下偏心环空物理模型与网格划分

Figure 1. Physical model and meshing of eccentric annulus with the eccentricity of 22.47%

下载: 全尺寸图片幻灯片

1.3 数值模拟方法有效性分析及参数设置

应用FLUENT进行三维偏心环空流动数值模拟分析时，需要确定的主要问题有：流态模型（层流或湍流模型）、离散格式、算法、非牛顿流体最大黏度及边界条件。笔者基于可以获取准确解析解的同心环空不旋转工况，确定上述参数。流速为1.50 m/s时，应用摩阻压降解析计算方法，得到雷诺数为2 300，层流摩阻压降计算结果（894.86 Pa/m）略大于紊流摩阻压降计算结果（853.00 Pa/m），故数值模拟研究时选择层流模型。为提高计算精度，选择二阶迎风格式，应用SIMPLE算法加速收敛。入口边界设定为速度入口（1.50 m/s），出口边界为压力出口，钻杆壁面和环空壁面为固体静止壁面，壁面处无滑移；在模拟钻柱旋转时，钻杆壁面处设定为对应的旋转角速度。调整钻井液黏度为最大值0.055 Pa·s时，数值模拟结果为896.91 Pa/m，与解析结果（894.86 Pa/m）相对误差为0.23%，故应用上述条件进行偏心旋转工况下的数值模拟分析。

2. 数值模拟结果分析

2.1 偏心环空流动特征分析

以偏心度53.93%和转速114.65 r/min的情况为例，在模型轴线4 m处建立一个截面，该界面处轴向、切向速度剖面如图2、图3所示。为进一步分析偏心和旋转对流速的影响，在该截面处建立一条自窄间隙到宽间隙的高边方向线，获取并绘制偏心度53.93%、转速114.65 r/min，偏心度53.93%、转速19.11 r/min，偏心度22.47%、转速114.65 r/min和偏心度22.47%、转速19.11 r/min等4种工况下的轴向和切向流速沿该线的分布，结果见图4和图5。为方便分析，定义X轴正方向为井筒高边方向(自钻杆下方窄环空间隙指向钻杆上方的宽环空间隙)，Y轴为截面中垂直于X轴的方向。

图 2 偏心环空轴向流速剖面

Figure 2. Axial velocity profile of eccentric annulus

下载: 全尺寸图片幻灯片

图 3 偏心环空切向流速剖面

Figure 3. Tangential velocity profile of eccentric annulus

下载: 全尺寸图片幻灯片

图 4 轴向流速沿高边方向线分布

Figure 4. Axial velocity distribution along the high-sidedirection line

下载: 全尺寸图片幻灯片

图 5 切向流速沿高边方向线分布

Figure 5. Tangential velocity distribution along the high-side direction line

下载: 全尺寸图片幻灯片

从图2可以看出：1）钻杆上方环空中间区带的轴向流速最高，钻杆下方窄环空区带的轴向流速最低；2）钻杆顺时针旋转时，Y轴右侧环空的轴向流速高于左侧环空；3）钻杆壁面和井壁处的轴向流速为0，与数值模拟中2个壁面处的无滑移边界条件一致。

从图3可以看出：1）距离钻杆壁面越近，钻井液的切向流速越大；2）钻杆壁面下部区域的钻井液切向流速最大，上部宽环空区带的切向流速最小。

从图4可以看出：1）钻杆上方宽环空区带和钻杆下方窄环空区带的钻井液轴向流速均基本呈抛物线分布，壁面处轴线流速为0，中间位置轴向流速最大；2）偏心度越大，宽环空区带钻井液轴向流速越大，窄环空区带钻井液轴向流速越小；3）偏心度越大，转速对窄环空区带的钻井液轴向流速分布影响越大；4）相同偏心度条件下，不同转速下的宽环空区带钻井液轴向流速分布基本重合，即转速对宽环空区带的钻井液轴向流速影响极小。

从图5可以看出：1）窄环空区带的钻井液切向流速大于宽环空区带；2）偏心度越大，宽环空区带钻井液切向流速越低，窄环空区带钻井液切向流速越大，这2个区域钻井液切向流速差异也越大；3）转速比偏心度对钻井液切向流速的影响更大。

钻柱偏心和旋转对环空钻井液流速分布的影响复杂，难以直接得到二者对摩阻压降的影响规律。但从岩屑运移的角度来说，钻柱偏心度越大，转速对窄环空区带钻井液轴向流速和切向流速的影响越大，而钻杆下方恰恰是岩屑堆积的重点区域，因此可以推断，提高转速可显著改善大斜度井段岩屑的运移效果。

2.2 摩阻压降梯度数值模拟计算结果

流体在流道内的流动通常分为流动发展阶段和稳定阶段，为了准确获取摩阻压降梯度数值模拟结果，沿模型中钻柱上方流域中心点建立一条轴线，分析压力沿该轴线的变化情况，以确定数值模拟的流动稳定阶段。钻井液静压力沿轴线变化的模拟结果见图6。

图 6 钻井液静压力沿宽流域中心线的分布

Figure 6. Static drilling fluid pressure distribution along the center line of wide flow domain

下载: 全尺寸图片幻灯片

从图6可以看出，流道入口附近压降梯度较大，距流道入口1.50 m之后，钻井液静压力与Z轴位置呈理想线性关系，这表明流道入口1.50 m内为流动发展阶段，之后进入流动稳定阶段。基于此，在模型轴线2.00和5.00 m处建立2个截面，分析2个截面之间的压差，求取摩阻压降梯度。不同偏心度和钻杆转速下的摩阻压降梯度计算结果见表1。

表 1 偏心环空摩阻压降梯度数值模拟计算结果

Table 1. Numerical computation results of frictional pressure drop gradient in eccentric annulus

偏心度，%	不同钻柱转速下的摩阻压降梯度/（Pa·m⁻¹）
偏心度，%	0 r/min	19.11 r/min	38.22 r/min	57.32 r/min	76.43 r/min	95.54 r/min	114.65 r/min
0	898.91	898.69	898.73	898.80	898.91	899.03	899.18
22.47	921.99	919.85	916.67	915.89	915.74	914.72	912.90
26.97	930.93	927.63	922.68	921.31	921.12	919.76	917.43
31.46	939.65	934.85	927.73	925.76	925.82	924.38	921.85
35.96	946.11	939.65	930.41	928.23	929.19	928.37	926.25
40.45	947.04	939.42	929.21	927.89	931.11	932.09	931.20
44.94	938.21	931.38	923.14	924.87	932.38	936.60	937.83
49.44	916.31	914.30	912.82	920.49	934.65	943.11	946.98
53.93	882.02	890.43	901.17	917.22	939.97	952.41	958.70
58.43	841.34	865.33	892.13	916.90	948.91	964.11	971.14
62.92	800.91	842.91	886.76	917.93	957.17	974.83	977.08
67.42	763.60	820.64	879.73	913.85	956.12	977.72	968.23

下载: 导出CSV

| 显示表格

2.3 偏心度对偏心环空摩阻压降梯度的影响

基于表1数据，绘制不同钻杆转速下环空摩阻压降梯度与偏心度的关系曲线，如图7所示（图中红色加粗直线为根据解析式计算出的同心不旋转摩阻压降梯度）。从图7可以看出：1）不同转速下，偏心度对环空摩阻压降梯度的影响规律不同，且没有单调关系；2）大多数情况下，偏心旋转环空摩阻压降梯度大于同心不旋转环空摩阻压降梯度，但在偏心度大于53.93%且转速小于38 r/min时，偏心旋转环空摩阻压降梯度低于同心不旋转环空摩阻压降梯度，也就是说，大斜度井段若采取滑动钻进方式且不考虑岩屑堆积造成的流动截面积减小，环空摩阻压降梯度会低于旋转钻进工况；3）不同转速下，偏心环空摩阻压降梯度与偏心度的关系曲线在偏心度约45%处存在一个交点，偏心度小于45%时，随偏心度增大，偏心环空摩阻压降梯度均略有增大；4）偏心度大于45%且转速高于57 r/min时，偏心环空摩阻压降梯度随偏心度增大而增大；5）偏心度大于45%且转速小于57 r/min时，偏心环空摩阻压降梯度随偏心度增大而减小。

图 7 偏心环空摩阻压降梯度与钻柱偏心度的关系曲线

Figure 7. Variation of frictional pressure drop gradient ineccentric annulus with eccentricity of drill string

下载: 全尺寸图片幻灯片

2.4 转速对偏心环空摩阻压降梯度的影响

基于表1数据，绘制不同钻杆偏心度下环空摩阻压降梯度与转速的关系曲线（见图8）。从图8可以看出：1）钻杆偏心度小于53.93%时，转速对环空摩阻压降梯度的影响较小，随着转速增大，偏心环空摩阻压降梯度略有降低；2）钻杆偏心度大于53.93%时，偏心环空摩阻压降梯度随转速增大而增大，但转速大于60 r/min后，转速对偏心环空摩阻压降梯度的影响减弱，也就是说，低转速下偏心度对摩阻压降梯度的影响远大于高转速情况。

图 8 偏心环空摩阻压降梯度与钻柱转速的关系曲线

Figure 8. Variation frictional pressure drop gradient in eccentric annulus with rotational speed of drill string

下载: 全尺寸图片幻灯片

3. 偏心环空摩阻压降梯度预测模型的建立

将表1中的摩阻压降梯度数据除以同心不旋转工况下的摩阻压降梯度，得到无因次偏心环空摩阻压降梯度（见表2）。由于转速和偏心度对环空摩阻压降梯度的影响规律复杂，不具有单调特征，若强行建立一个综合考虑转速和偏心度的二元关系模型R=f（N，E），则预测精度会比较低。为此，按照偏心度和转速进行分类，分别建立偏心环空摩阻压降梯度预测模型。

表 2 无因次偏心环空摩阻压降梯度数值模拟结果

Table 2. Numerical simulation results of dimensionless eccentric frictional pressure drop gradient

偏心度，%	不同钻柱转速下的无因次偏心环空摩阻压降梯度
偏心度，%	0 r/min	19.11 r/min	38.22 r/min	57.32 r/min	76.43 r/min	95.54 r/min	114.65 r/min
0	1.0000	0.9998	0.9998	0.9999	1.0000	1.0001	1.0003
22.47	1.0257	1.0233	1.0198	1.0189	1.0187	1.0176	1.0156
26.97	1.0356	1.0320	1.0264	1.0249	1.0247	1.0232	1.0206
31.46	1.0453	1.0400	1.0321	1.0299	1.0299	1.0283	1.0255
35.96	1.0525	1.0453	1.0351	1.0326	1.0337	1.0328	1.0304
40.45	1.0535	1.0451	1.0337	1.0322	1.0358	1.0369	1.0359
44.94	1.0437	1.0361	1.0270	1.0289	1.0372	1.0419	1.0433
49.44	1.0194	1.0171	1.0155	1.0240	1.0398	1.0492	1.0535
53.93	0.9812	0.9906	1.0025	1.0204	1.0457	1.0595	1.0665
58.43	0.9360	0.9626	0.9925	1.0200	1.0556	1.0725	1.0804
62.92	0.8910	0.9377	0.9865	1.0212	1.0648	1.0845	1.0870
67.42	0.8495	0.9129	0.9787	1.0166	1.0636	1.0877	1.0771

下载: 导出CSV

| 显示表格

建立心环空摩阻压降梯度预测模型时，定义了无因次旋转因子 ${R}_{N}$ 和无因次偏心因子 ${R}_{E}$ ：

${R}_{N}=\dfrac{{\Delta p}_{{\rm{f}},N}}{\Delta {p}_{{\rm{f}},0}}$

(2)

${R}_{E}=\dfrac{{\Delta p}_{{\rm{f}},E}}{\Delta {p}_{{\rm{f}},0}}$

(3)

式中：R_N为无因次旋转因子，即相同偏心度下某转速时的摩阻压降梯度与基准摩阻压降梯度（即转速为0、偏心度为0时的摩阻压降梯度）的比值；R_E为无因次偏心因子，即相同转速下某偏心度时的摩阻压降梯度与基准摩阻压降梯度的比值； $\Delta {p}_{{{\rm{f}}},0}$ 为同心不旋转工况下的幂律流体摩阻压降梯度，Pa/m； $\Delta {p}_{{{\rm{f}}},N}$ 为相同偏心度不同转速时的摩阻压降梯度，Pa/m； $\Delta {p}_{{{\rm{f}}},E}$ 相同转速不同偏心度时的摩阻压降梯度，Pa/m。

将2类偏心环空摩阻压降梯度预测模型的预测结果与数值模拟结果进行比较发现，偏心度分类预测模型的残差平方和为3.61×10⁻⁴，而转速分类预测模型的残差平方和为2.25×10⁻³，前者的预测精度比后者高一个数量级，所以最终确定将表2中的数据进行回归，建立偏心度分类环空摩阻压降梯度预测模型。

1）E=0时的预测模型为：

${R}_{N}=7.254\times 1{0}^{-8}{N}^{2}-3.991\times 1{0}^{-6}N+0.999\,\,9$

(4)

2）E=22.47%时的预测模型为：

${R}_{N}=3.965\times 1{0}^{-7}{N}^{2}-2.369\times 1{0}^{-4}N+1.025\,\,0$

(5)

3）E=26.97%时的预测模型为：

${R}_{N}=8.114\times 1{0}^{-7}{N}^{2}-1.247\times 1{0}^{-4}N+1.035\,\,0$

(6)

4）E=31.46%时的预测模型为：

${R}_{N}=1.588\times 1{0}^{-6}{N}^{2}-2.405\times 1{0}^{-4}N+1.045\,\,0$

(7)

5）E=35.96%时的预测模型为：

${R}_{N}=2.539\times 1{0}^{-6}{N}^{2}-3.645\times 1{0}^{-4}N+1.052\,\,0$

(8)

6）E=40.45%时的预测模型为：

${R}_{N}=3.579\times 1{0}^{-6}{N}^{2}-4.356\times 1{0}^{-4}N+1.053\,\,0$

(9)

7）E=44.94%时的预测模型为：

$\begin{aligned}{R}_{N}=&-6.572\times 1{0}^{-8}{N}^{3}+1.544\times 1{0}^{-5}{N}^{2}-\\ &9.16\times 1{0}^{-4}N+1.045\,\,0 \end{aligned}$

(10)

8）E=49.44%时的预测模型为：

$\begin{aligned}{R}_{N}=&-8.842\times 1{0}^{-8}{N}^{3}+1.855\times 1{0}^{-5}{N}^{2}-\\ &5.72\times 1{0}^{-4}N+1.021\,\,0 \end{aligned}$

(11)

9）E=53.93%时的预测模型为：

$\begin{aligned}{R}_{N}=&-1.067\times 1{0}^{-7}{N}^{3}+1.875\times 1{0}^{-5}{N}^{2}-\\ &9.066\times 1{0}^{-6}N+0.982\,\,2 \end{aligned}$

(12)

10）E=58.43%时的预测模型为：

$\begin{aligned}{R}_{N}=&-1.141\times 1{0}^{-7}{N}^{3}+1.496\times 1{0}^{-5}{N}^{2}+\\ &1.039\times 1{0}^{-3}N+0.936\,\,6\end{aligned}$

(13)

11）E=62.92%时的预测模型为：

$\begin{aligned}{R}_{N}=&-1.158\times 1{0}^{-7}{N}^{3}+8.551\times 1{0}^{-6}{N}^{2}+\\ &2.252\times 1{0}^{-3}N+0.891\,\,6 \end{aligned}$

(14)

12）E=67.42%时的预测模型为：

$\begin{aligned}{R}_{N}=&-1.277\times 1{0}^{-7}{N}^{3}+3.688\times 1{0}^{-6}{N}^{2}+\\ &3.251\times 1{0}^{-3}N+0.850\,\,3 \end{aligned}$

(15)

式中：N为转速，r/min。

利用上述预测模型计算偏心环空摩阻压降梯度的主要步骤为：

1）基于钻井液流变性和流道尺寸，根据解析式^[26]计算同心不旋转工况下的幂律流体摩阻压降梯度 $\Delta {p}_{{{\rm{f}}},0}$ ；

2）根据计算点所在井段特征（增斜段、降斜段还是稳斜段），以及井斜角、造斜率、钻柱尺寸、井眼尺寸和稳定器尺寸等参数，依据钻柱屈曲理论^[17,27-28]，确定偏心度E；

3）查找偏心度E在12个偏心度的区间[E₁，E₂]，选取相应的预测模型计算偏心度E₁和E₂时的无因次旋转因子 ${R}_{{N},1}$ 和 ${R}_{{N},2}$ ，应用插值法计算偏心度E时的无因次旋转因子 ${R}_{{N}}$ ；

4）根据式（2）计算偏心环空摩阻压降梯度 ${\Delta p}_{{{\rm{f}}},{N}}$ 。

4. 模型预测结果验证

南海某水平井水平段的钻具组合为：ϕ215.9 mm钻头+ ϕ171.0 mm旋转导向工具+ϕ206.0 mm倒划眼稳定器+ϕ165.1 mm液压震击器+ϕ127.0 mm加重钻杆8根+ϕ127.0 mm钻杆。钻井液密度为1 280 kg/m³，稠度系数为0.72 Pa·sⁿ，流性指数为0.66；钻井液排量为0.033 m³/s，钻柱转速为110 r/min。

基于上述数据，首先，应用Haciislamoglu方法^[17]计算了偏心度为20%、30%、40%和50%时文中提出的无因次摩阻压降梯度系数 ${R}_{\mathrm{l}\mathrm{a}\mathrm{m}}$ （对应文中的R_N），分别为0.94、0.98、0.98和0.98，均小于1.0，与现场工况不符，进一步印证了基于有限室内试验数据建立的回归模型的实用性不强。之后，应用建立的偏心环空摩阻压降梯度预测模型，以10 m为距离步长，计算了该井钻井液循环工况下ϕ215.9 mm井段偏心环空摩阻压降。为方便与现场PWD实测数据进行对比，按照已有公式计算得到了该井段的当量循环密度，计算值与实测值对比见图9，平均相对误差为0.45%，验证了模型的准确性。

图 9 南海某水平井ϕ215.9 mm井段当量循环密度计算值与实测值对比

Figure 9. Comparison between calculated and measured values of ECD in ϕ215.9 mm section of a horizontal well in the South China Sea

下载: 全尺寸图片幻灯片

钻井液循环工况下某井深处当量循环密度的计算公式为^[29]：

$\rho_{{\rm{ce}}}=\dfrac{{p}_{{\rm{h}}}+{p}_{{\rm{f}}}}{{D}_{{\rm{v}}}}$

(16)

$\text{其中}\qquad\qquad\qquad {p}_{{\rm{f}}}=\sum _{i=1}^{N}{\Delta {L}_{i}\Delta p}_{{\rm{f}},i}$

(17)

式中：ρ_ce为当量循环密度，kg/m³； ${p}_{{\rm{h}}}$ 为静液柱压力，Pa； ${p}_{{\rm{f}}}$ 为摩阻压降，根据井眼轨迹分段计算并累计求和，Pa； ${D}_{{\rm{v}}}$ 为目标点垂深，m。

5. 结论与认识

1）三维环空流动数值模拟结果显示：偏心度越大，转速对窄环空区带的钻井液轴向流速分布影响越大，但转速对宽环空区带的钻井液轴向流速影响极小；窄环空区带的钻井液切向流速大于宽环空区带，偏心度越大，宽环空区带钻井液切向流速越低，窄环空区带的钻井液切向流速越大，这2个区域的钻井液切向流速差异也越大，转速比偏心度对钻井液切向流速的影响更大。

2）模拟计算结果显示，偏心环空摩阻压降梯度与偏心度、钻柱转速之间的关系复杂，不存在单调关系，难以建立综合考虑转速和偏心度的具有较高预测精度的二元关系模型；而按照偏心度和转速进行分类，分别建立偏心环空摩阻压降梯度预测模型，更有利于工程应用。

3）对比预测模型的预测结果与数值模拟结果发现，偏心度分类预测模型的预测精度比转速分类预测模型的预测精度高一个数量级，为此，通过对数值模拟数据回归分析，建立了偏心度分类环空摩阻压降梯度预测模型。

4）应用偏心度分类环空摩阻压降梯度预测模型，计算了南海某水平井的ϕ215.9 mm井段当量循环密度，并与PWD测试结果进行了对比，平均相对误差为0.45%，验证了模型的准确性，可用于精细描述钻井环空循环压力场和准确计算当量循环密度。

参考文献(12)

[1]	李爱英.废钻井液污染分析及处理方法的探讨[J].油气田环境保护,1998,8(2):15-17. Li Aiying.Discussion of waste drilling fluid contamination analysis and processing methods[J].Environmental Protection of Oil Gas Fields,1998,8(2):15-17.
[2]	苏勤,何青水,张辉,等.国外陆上钻井废弃物处理技术[J].石油钻探技术,2010,38(5):106-110. Su Qin,He Qingshui,Zhang Hui.Foreign onshore drilling waste treatment technology[J].Petroleum Drilling Techniques,2010,38(5):106-110.
[3]	张兴儒,吴振烈.油气田环境保护[M].北京:石油工业出版社,1994:15-21. Zhang Xingru,Wu Zhenlie.Environmental protection of oil gas fields[M].Beijing:Petroleum Industry Press,1994:15-21.
[4]	薛玉志,马云谦,李公让,等.海上废弃钻井液处理研究[J].石油钻探技术,2008,36(5):19-23. Xue Yuzhi,Ma Yunqian,Li Gongrang,et al.Offshore waste drilling fluid treatment[J].Petroleum Drilling Techniques,2008,36(5):19-23.
[5]	赵雄虎,王风春.废弃钻井液处理研究进展[J].钻井液与完井液,2004,21(2):43-48. Zhao Xionghu,Wang Fengchun.Research development of waste drilling fluids disposal[J].Drilling Fluid Completion Fluid,2004,21(2):43-48.
[6]	朱墨,夏斌,张进.废钻井液对环境污染及固化处理室内研究[J].钻井液与完井液,1995,12(3):8-14. Zhu Mo,Xia Bin,Zhang Jin.Laboratory study on the innocuous treatment of waste drilling fluid[J].Drilling Fluid Completion Fluid,1995,12(3):8-14.
[7]	乔东宇,张兴国,王云,等.新疆油田废弃钻井液固化处理技术[J].新疆石油天然气,2011,7(4):82-85. Qiao Dongyu,Zhang Xingguo,Wang Yun,et al.Waste drilling fluid treating technology in Xinjiang Oilfield[J].Xinjiang Oil Gas,2011,7(4):82-85.
[8]	汤超,邓皓,王蓉沙,等.废弃钻井液处理技术[J].石油化工腐蚀与防护,2010,27(2):21-24. Tang Chao,Deng Hao,Wang Rongsha,et al.Research on treatment technologies of waste drilling fluid[J].Corrosion Protection in Petrochemical Industry,2010,27(2):21-24.
[9]	Getliff J M,Silverstone M P,Shearman A K,et al.Waste management and disposal of cuttings and drilling fluid waste resulting from the drilling and completion of wells to produce orinoco very heavy oil in Eastern Venezuela.SPE 46600,1998.
[10]	Ozumba C I,Benebo T E T.Waste recycling initiatives in an exploration company in Nigeria.SPE 73841,2002.
[11]	张强,刘亮,党宁海.废弃钻井液处理技术研究[J].大观周刊,2012(40):114-115. Zhang Qiang,Liu Liang,Dang Ninghai.Research on treatment technologies of waste drilling fluid[J].Daguan Weekly,2012(40):114-115.
[12]	许海彬.钻井废泥浆和固化材制成的泵房:中国,2466286[P].2001-12-19. Xu Haibin.Pumping stations and curing materials made of mud:CN,2466286[P].2001-12-19.

施引文献(6)

期刊类型引用(6)

1.	沈吉阳，孙晨皓，王建民，张木楠，邹树江，彭勇. 基于钻柱动力学大位移井钻具接头可靠性研究. 石油矿场机械. 2024(02): 1-9 . 百度学术
2.	陈雨飞，安锦涛，张辉，李军，周英操，路宗羽. 考虑钻杆接头的小井眼环空压耗计算模型. 钻井液与完井液. 2024(03): 296-304 . 百度学术
3.	王彪，李军，杨宏伟，詹家豪，张更，龙震宇. 基于工程参数变化趋势的溢流早期智能检测方法. 石油钻探技术. 2024(05): 145-153 . 本站查看
4.	赵金，赵星，蔡鹏，彭齐，饶嘉骐. 连续管多孔喷嘴射流工具研制与清洗参数优化. 石油钻探技术. 2023(03): 83-89 . 本站查看
5.	李文拓，罗鸣，黄洪林，李军，肖平. 高温高压小井眼水平井环空ECD综合计算模型. 石油钻采工艺. 2023(03): 259-268 . 百度学术
6.	龚勇，蒋凯，闫许峰，康建平，林黎明. 钻具旋转时同心窄间隙环空流动的数值模拟. 钻探工程. 2023(06): 56-62 . 百度学术

其他类型引用(0)

资源附件(0)

计量

文章访问数: 2814
HTML全文浏览量: 96
PDF下载量: 3807
被引次数: 6

1. 三维偏心环空流动数值模拟
1.1 数值模拟参数
1.2 物理模型及网格划分
1.3 数值模拟方法有效性分析及参数设置
2. 数值模拟结果分析
2.1 偏心环空流动特征分析
2.2 摩阻压降梯度数值模拟计算结果
2.3 偏心度对偏心环空摩阻压降梯度的影响
2.4 转速对偏心环空摩阻压降梯度的影响
3. 偏心环空摩阻压降梯度预测模型的建立
4. 模型预测结果验证
5. 结论与认识

1. 三维偏心环空流动数值模拟
1.1 数值模拟参数
1.2 物理模型及网格划分
1.3 数值模拟方法有效性分析及参数设置
2. 数值模拟结果分析
2.1 偏心环空流动特征分析
2.2 摩阻压降梯度数值模拟计算结果
2.3 偏心度对偏心环空摩阻压降梯度的影响
2.4 转速对偏心环空摩阻压降梯度的影响
3. 偏心环空摩阻压降梯度预测模型的建立
4. 模型预测结果验证
5. 结论与认识

参考文献(12)

施引文献

资源附件(0)

环境敏感区废弃钻井液无害化处理新技术

作者简介: 王勇（1971—），男，河北南宫人，1993年毕业于江汉石油学院钻井工程专业，在读博士研究生，高级工程师，主要从事钻井技术研究及相关管理工作。

计量

出版历程