泸州区块深层页岩气水平井优快钻井技术

李涛; 杨哲; 徐卫强; 阳强; 杨兆亮

doi:10.11911/syztjs.2022036

泸州区块深层页岩气水平井优快钻井技术

1.
中国石油西南油气田分公司工程技术研究院, 四川成都610017
2.
中国石油西南油气田分公司开发事业部, 四川成都 610066

基金项目: 中国石油集团重大现场试验项目“深层页岩气水平井优快钻完井技术现场试验”（编号：2019F-31-03）和中国石油西南油气田分公司科技项目“川渝地区深层页岩气水平井钻井降温技术研究与应用” （编号：20200302-09）联合资助

详细信息

作者简介:
李涛（1993—），男，四川成都人，2015年毕业于西南石油大学石油工程专业，工程师，主要从事钻井完井技术研究与钻井故障及复杂管理工作。E-mail：litao01@petrochina.com.cn

中图分类号: TE243⁺.1
计量
- 文章访问数: 576
- HTML全文浏览量: 183
- PDF下载量: 242
出版历程
- 收稿日期: 2022-03-12
- 修回日期: 2022-08-23
- 网络出版日期: 2022-11-06
- 刊出日期: 2023-01-31

Optimized and Fast Drilling Technology for Deep Shale Gas Horizontal Wells in Luzhou Block

1.
Engineering Technology Research Institute, PetroChina Southwest Oil & Gasfield Company, Chengdu, Sichuan, 610017, China
2.
Development Division, PetroChina Southwest Oil & Gasfield Company, Chengdu, Sichuan, 610066, China

摘要

摘要:
泸州区块深层页岩气储层龙马溪组页岩埋藏深、温度高，钻井过程中面临井漏频发、难钻地层机械钻速低、旋转导向仪器高温下的故障率高和钻井周期长等技术难点。为此，在分析钻遇地层特征和钻井技术难点的基础上，进行了井身结构和井眼轨道优化、高效钻头优选、堵漏技术措施优化、油基钻井液地面降温技术和欠平衡钻井技术研究，形成了泸州区块深层页岩气水平井优快钻井技术。该技术在泸州区块34口页岩气水平井进行了应用，应用井平均完钻井深5 760 m，平均水平段长1 890 m，平均机械钻速由5.61 m/h提高至7.03 m/h，创造了造斜段+水平段一趟钻完钻、单趟钻进尺2 330 m等纪录。泸州区块深层页岩气水平井优快钻井技术支撑了该区块深层页岩气安全钻井和钻井提速提效，对其他地区深层页岩气水平井钻井提速具有一定的借鉴价值。
- 深层页岩气 /
- 水平井 /
- 钻井提速 /
- 井漏 /
- 钻头优选 /
- 钻井液降温 /
- 欠平衡钻井 /
- 泸州区块
Abstract:
The Longmaxi shale gas reservoirs in Luzhou are deeply buried and have a high temperature. During the drilling process, technical problems such as frequent lost circulation, low rate of penetration (ROP), high failure rate of steering instruments due to high temperature, and long drilling period are frequently encountered. Therefore, on the basis of analyzing the actual formation conditions and technical challenges, the research on optimization of casing program and wellbore trajectory, plugging treatment improvement for lost circulation, selection of efficient bits, surface cooling of oil-based drilling fluid and underbalanced drilling technologies were conducted. As a result, an optimized and fast drilling technology for deep shale horizontal well in Luzhou Block was developed. This technology was applied to 34 shale gas horizontal wells in Luzhou Block. Average drilling depth of the wells was 5 760 m and the average length of horizontal sections was 1 890 m. In addition, the average ROP was increased from 5.61 m/h to 7.03 m/h. The inclined and horizontal sections were completed in one trip and the drilling footage of a single trip reached 2 330 m, which created a new record. The optimized and fast drilling technology for deep shale gas horizontal wells in Luzhou Block has ensured the safe drilling of deep shale gas, improved ROP and efficiency. It can provide a technical reference for ROP improvement of deep shale gas horizontal wells in other regions.
- deep shale gas /
- horizontal well /
- ROP improvement /
- lost circulation /
- bit selection /
- cooling of drilling fluid /
- underbalanced drilling /
- Luzhou Block

HTML全文

很多陆上探井的完钻井深超过8 000 m，如塔里木油田轮探1井和满深3井的完钻井深分别为8 882和8 010 m；顺北油田鹰1井、蓬1井和顺北5-5H井的完钻井深分别为8 588，8 450和8 520 m^[1-2]。随着井深增加，地层温度快速上升，井深8 000 m处的初始地层温度已超过180 ℃。井下高温对保持钻井液和固井水泥浆性能、防止钻具腐蚀产生不利影响。因此，国内外学者基于数值模型和解析模型预测和分析钻进中的井筒温度分布特征，为钻前优化钻井液等工作液的性能、优选钻具和测量工具提供依据。

其中，通过数值模型求解，即基于能量守恒原理，考虑井筒–地层各控制区域在径向上和轴向上的热交换机理，建立井筒–地层瞬态传热模型（井筒–地层传热数值模型），并应用有限差分法和有限元法对该模型求解^[3]。截至目前，相关学者已开展了直井和水平井井筒温度敏感因素分析^[4-5]，结合不同钻井完井工况，论述了溢流、非循环时间及水泥浆水化热对井筒温度分布的影响^[6-9]。为了降低数值模型的求解难度，A. R. Hasan 等人^[10-11]应用无因次时间温度，将井筒–地层间热交换表征为拟稳态传热过程，钻柱内和环空中的流体为稳态传热过程，建立了井筒拟稳态传热模型（井筒–地层传热解析模型）。基于此方法，唐林等人^[12]研究了循环时间、注入排量对超深井钻井井筒温度的影响，窦亮彬等人^[13]分析了影响注CO₂井筒温度的因素，Nian Yongle等人^[14]评价了废弃油井的地热资源潜能，李勇等人^[15-16]考察了漏失工况下井筒温度分布以及气井井筒温度对水泥环应力的影响特征。

文献调研可知^[17-20]，数值模型和解析模型广泛应用于钻井、开采中井筒温度分布特征分析，为高温井井筒流体流动安全与压力控制提供了重要的分析方法。井下温度需要精细化评价，如水泥浆温度差值为5 ℃时，其稠化时间有1 h左右的误差，若井下温度预测过低，可能导致水泥浆在套管内提前凝固。同时，井下循环温度预测过低会缩短定向工具使用寿命，增加钻井成本。此外，准确预测井下温度，可为研发、优选钻井液、水泥浆、完井液及压裂液等工作液提供依据^[21]。然而，目前尚未深入分析井筒–地层瞬态传热模型和井筒拟稳态传热模型计算精度及影响其计算精度的因素，对传热模型的应用条件与范围还不清楚。为此，笔者基于微分单元能量守恒原理，建立了钻井流体循环期间井筒–地层传热数值模型和井筒–地层传热解析模型，并对其进行了求解，然后以顺北油田某超深井井身结构、钻井参数及测温数据为基础，分析了2种模型的井筒温度计算精度及其影响因素，以期为科学实时评价超深井井筒温度提供依据。

1. 井筒–地层传热物理模型

钻进时，钻井流体循环过程中井筒–地层各控制区域（包括钻柱内流体、钻柱壁、环空流体、井壁及近井壁地层）为耦合体，一个控制区域热量的变化会影响其他区域热量的改变。因此，井筒–地层各控制区域发生复杂的热交换方式，会使钻柱内和环空流体温度随钻进时间变化。

研究井筒流体温度变化时，需建立井筒–地层瞬态传热模型并进行求解，以获得钻柱内和环空流动流体的温度变化特征。钻进时钻柱内与环空流体温度的变化情况（物理模型）如图1所示（图1中：q_p(z)为钻柱内流体向下流入微元体的热量，J；Δz为微元体长度，m；q_p(z+Δz)为钻柱内流体向下流出微元体的热量，J；q_pf为单位长度钻柱内流体循环摩阻产生的热量，J；q_pa为钻柱内流体与钻柱内壁面交换的热量，J；q_a(z+Δz)为环空内流体向上流入微元体的热量，J；q_a(z)为环空内流体向上流出微元体的热量，J；q_af为单位长度环空内流体循环摩阻产生的热量，J；q_Fa为从井壁传入环空的热量，J；q_ap为从环空传入钻柱的热量，J；t_p为钻柱内流体温度，℃；t_w为钻柱壁温度，℃；t_a为环空流体温度，℃；t_c1为井壁温度，℃；t_i，j为地层中空间网格第i行第j列的温度，℃）。

图 1 钻进中钻柱内与环空流体温度的变化

Figure 1. Fluid temperature variations in the drill string and annulus during drilling

下载: 全尺寸图片幻灯片

钻井液在地面以温度t_in从钻柱内向下流动，钻柱内流动传热控制单元的热量变化包括：1）向下流入的热量q_p(z)与流出热量之差q_p(z+Δz)；2）钻柱内壁以热对流方式传入钻柱内流体的热量q_pa；3）流体循环摩阻产生的热量q_pf。流体向下流经钻头喷嘴后流入环空，在井底钻柱内流体的温度与环空流体的温度相同。在环空向上流动过程中，控制单元热量变化包括：1）向上流入的热量q_a(z+Δz)与流出热量之差q_a(z)；2）钻柱外壁以热对流方式传入钻柱壁的热量q_ap；3）井壁以热对流方式传入环空的热量q_Fa；4）流体循环摩阻产生的热量q_af。流体经一个循环周期从环空返至地面。

2. 井筒–地层传热数值模型

基于井筒–地层各控制区域的能量守恒原理，分别建立了钻柱内流体、钻柱壁、环空流体及近井壁地层温度数值模型（统称为井筒–地层传热数值模型），并采用全隐式有限差分法求解了该模型。

2.1 钻柱内流体传热模型

结合图1所示物理模型反映的传热机理，建立钻柱内传热模型^[6]。在时间dt内，钻柱内钻井液向下流动携带进入单元体的热量为：

${\rm{d}}{q_{\text{1}}} = {q_{{\rm{p}}(z)}} - {q_{{\rm{p}}(z + \Delta z)}} = {\rho _{\text{m}}}Q{c_{\text{m}}}\left| {{t_{{\rm{p}}(z)}} - {t_{{\rm{p}}(z + \Delta z)}}} \right|{\rm{d}}t$

(1)

在时间dt内，钻柱内钻井液与钻柱内壁在径向上对流换热产生的热量为：

${\rm{d}}{q_{\text{2}}} = 2{\text{π}}{r_{\text{0}}}{h_{{\text{ci}}}}({t_{\text{w}}} - {t_{\text{p}}})\Delta z{\rm{d}}t$

(2)

令单位长度的功率为Q_p，则外界对单元体做功产生的热量为：

${\rm{d}}W = {Q_{\text{p}}}\Delta z{\rm{d}}t$

(3)

在时间dt内，钻柱内单元体流体内能的变化为：

${\rm{d}}E = {\rho _{\text{m}}}{c_{\text{m}}}{\text{π}}{r_{\text{0}}}\Delta z\frac{{\partial {t_{\text{p}}}}}{{\partial t}}{\rm{d}}t$

(4)

基于能量守恒，对式（1）—式（4）整理可得：

$\begin{split} {\rho _{\text{m}}}Q{c_{\text{m}}}&\left| {{{t}_{{\rm{p}}(z)}}- {{t}_{{\rm{p}}(z + \Delta z)}}} \right|{\rm{d}}t + 2{\text{π}}{r_{\text{0}}}{h_{{\text{ci}}}}({t_{\text{w}}} - {t_{\text{p}}})\Delta z{\rm{d}}t + \\ &\quad {Q_{\text{p}}}\Delta z{\rm{d}}t = {\rho _{\text{m}}}{c_{\text{m}}}{\text{π}}r_{\text{0}}^2\Delta z\frac{{\partial {t_{\text{p}}}}}{{\partial t}}{\rm{d}}t \end{split}$

(5)

等式两边同时除以 $\Delta z{\rm{d}}t$ ，则有：

$\begin{split} - {\rho _{\text{m}}}Q{c_{\text{m}}}& \frac{{\left| {{t_{{\rm{p}}(z + \Delta z)}} - {t_{{\rm{p}}(z)}}} \right|}}{{\Delta z}} + 2{\text{π}}{r_{\text{0}}}{h_{{\text{ci}}}}({t_{\text{w}}} - {t_{\text{p}}}) + {Q_{\text{p}}} = \\ &\qquad\qquad {\rho _{\text{m}}}{c_{\text{m}}}{\text{π}}r_{\text{0}}^2\frac{{\partial {t_{\text{p}}}}}{{\partial t}} \end{split}$

(6)

由微分基本定义， $\displaystyle\frac{{\left| {{t_{{\rm{p}}(z + \Delta z)}} - {t_{{\rm{p}}(z)}}} \right|}}{{\Delta z}} \approx \frac{{\partial {t_{\text{p}}}}}{{\partial z}}$ ，整理可得：

${Q_{\text{p}}} - {\rho _{\text{m}}}Q{c_{\text{m}}}\frac{{\partial {t_{\text{p}}}}}{{\partial z}} - 2{\text{π}}{r_{\text{0}}}{h_{{\text{ci}}}}({t_{\text{p}}} - {t_{\text{w}}}) = {\rho _{\text{m}}}{c_{\text{m}}}{\text{π}}r_{\text{0}}^2\frac{{\partial {t_{\text{p}}}}}{{\partial t}}$

(7)

式中： $\;{\rho _{\text{m}}}$ 为钻井液密度，kg/m³； ${c_{\text{m}}}$ 为钻井液比热容，J/（kg·℃）； ${h_{{\text{ci}}}}$ 为钻柱内壁面对流换热系数，W/（m²·℃）； $Q$ 为钻井液排量，m³/s； ${Q_{\text{p}}}$ 为钻柱内流体热源项，W/m； ${r_{\text{0}}}$ 为钻柱内径，m；t为时间，s。

2.2 钻柱壁传热模型

钻柱壁温度与钻柱内钻井液和环空钻井液的热对流相关。钻柱壁热量变化因素包括：轴向上，钻柱热传导产生的热量；径向上，在钻柱外壁和钻柱内部由热对流传入钻柱壁和钻柱内流体的热量；单位时间内单元体内能的变化。用数学模型可表示为：

${\lambda _{\text{w}}}\frac{{{\partial ^2}{t_{\text{w}}}}}{{\partial {z^2}}} + \frac{{2{r_{\text{0}}}{h_{{\text{ci}}}}}}{{r_{\text{1}}^2 - r_{\text{0}}^2}}({t_{\text{p}}} - {t_{\text{w}}}) - \frac{{2{r_{\text{1}}}{h_{{\text{co}}}}}}{{r_{\text{1}}^2 - r_{\text{0}}^2}}({t_{\text{w}}} - {t_{\text{a}}}) = {\rho _{\text{w}}}{c_{\text{w}}}\frac{{\partial {t_{\text{w}}}}}{{\partial t}}$

(8)

式中： ${\rho _{\text{w}}}$ 为钻柱密度，kg/m³； ${c_{\text{w}}}$ 为钻柱比热容，J/（kg·℃）； ${\lambda _{\text{w}}}$ 为钻柱导热系数，W/（m·℃）； ${h_{{\text{co}}}}$ 为钻柱外壁面对流换热系数，W/（m²·℃）； ${r_{\text{1}}}$ 为钻柱外径，m。

2.3 环空流体传热模型

影响环空钻井液温度的因素包括：1）钻井液沿井筒向上流动时携带的热量；2）井壁和钻柱外壁发生热对流获得的热量；3）外界对单元体流体所做的功；4）钻井液自身内能的变化。据此，环空流体传热数学模型可表示为：

$\begin{split} {\rho _{\text{m}}}Q{c_{\text{m}}}\frac{{\partial {t_{\text{a}}}}}{{\partial z}} + 2{\text{π}} & {r_{\text{1}}}{h_{{\text{co}}}}({t_{\text{w}}} - {t_{\text{a}}}) - 2{\text{π}}{r_{\text{2}}}{h_{\text{b}}}({t_{\text{a}}} - {t_{{\text{c1}}}}) + {Q_{\text{a}}} = \\ & {\rho _{\text{m}}}{c_{\text{m}}}{\text{π}}(r_{\text{2}}^2 - r_{\text{1}}^2)\frac{{\partial {t_{\text{a}}}}}{{\partial t}} \\[-15pt] \end{split}$

(9)

式中： $h{}_{\text{b}}$ 为井壁对流换热系数，W/（m²·℃）； ${r_{\text{2}}}$ 为井壁半径，m； ${Q_{\text{a}}}$ 为环空内流体热源项，W/m。

2.4 井壁传热模型

井壁由套管段和裸眼段（岩石）构成，则影响井壁温度的因素有：轴向上，热传导获得的热量；径向上，与内层水泥环/岩石发生热传导、与环空钻井液产生对流换热；自身内能的变化。据此，井壁传热数学模型可表示为：

${\lambda _{{\text{c1}}}}\frac{{{\partial ^2}{t_{{\text{c1}}}}}}{{\partial {z^2}}} + \frac{{2{r_{\text{2}}}{h_{\text{b}}}({t_{\text{a}}} - {t_{{\text{c1}}}})}}{{r_{\text{3}}^2 - r_{\text{2}}^2}} - \frac{{2{\lambda _{{\text{o1}}}}({t_{{\text{c1}}}} - {t_{{\text{c2}}}})}}{{r_{\text{3}}^2 - r_{\text{2}}^2}} = \rho {}_{{\text{c1}}}{c_{{\text{c1}}}}\frac{{\partial {t_{{\text{c1}}}}}}{{\partial t}}$

(10)

式中： ${\lambda _{{\text{c1}}}}$ 为套管/岩石的导热系数，W/（m·℃）； ${\lambda _{{\text{o1}}}}$ 为水泥环/岩石的导热系数，W/（m·℃）； ${t_{{\text{c2}}}}$ 为井水泥环/岩石温度，℃； $\rho {}_{{\text{c1}}}$ 为套管/岩石的密度，kg/m³； ${c_{{\text{c1}}}}$ 为套管/岩石的比热容，J/（kg·℃）； ${r_{\text{3}}}$ 为套管外径，m。

2.5 近井壁地层传热模型

近井壁各控制区域包括多层套管、水泥环及地层，其热交换可看作单一圆筒壁在轴向和径向上的热传导。据此，传热模型可表示为：

$\frac{{{\partial ^2}{t_{{\rm{f}}k}}}}{{\partial {r^2}}} + \frac{1}{{{r_k}}}\frac{{\partial {t_{{\rm{f}}k}}}}{{\partial r}} + \frac{{{\partial ^2}{t_{{\rm{f}}k}}}}{{\partial {z^2}}} = \frac{{{\rho _{{\rm{f}}k}}{c_{{\rm{f}}k}}}}{{{\lambda _{{\rm{f}}k}}}}\frac{{\partial {t_{{\rm{f}}k}}}}{{\partial t}}$

(11)

式中： ${t_{{{\rm{f}}_k}}}$ 为套管、水泥环或岩石的温度，℃； ${\rho _{{\rm{f}}k}}$ 为套管、水泥环或岩石的密度，kg/m³； ${c_{{{\rm{f}}_k}}}$ 为套管、水泥环或岩石的比热容，J/（kg·℃）； ${\lambda _{{{\rm{f}}_k}}}$ 为套管、水泥环或岩石的导热系数，W/（m·℃）；k为套管层级大小，通常4 ≤ k ≤ 11。

2.6 井筒–地层传热数值模型求解

上述井筒–地层传热数值模型涉及多个偏微分方程组。对此，采用全隐式有限差分法对该模型进行求解：

${\alpha _{ij}}t_{i - 1,j}^{n + 1} + {\beta _{ij}}t_{i,j}^{n + 1} + {\gamma _{ij}}t_{i + 1,j}^{n + 1} + {\delta _{ij}}t_{i,j - 1}^{n + 1} + {\lambda _{ij}}t_{i,j + 1}^{n + 1} = {c_{ij}}$

(12)

把所有节点的有限差分方程用矩阵形式表示：

${\boldsymbol{AT}}={\boldsymbol{C}}$

(13)

式中： ${\alpha _{ij}}$ ， ${\beta _{ij}}$ ， ${\gamma _{ij}}$ ， ${\delta _{ij}}$ 和 ${\lambda _{ij}}$ 均为矩阵系数；i, j分别为径向和轴向上的空间网格节点；n为时间节点。

结合式（12）的离散形式，将上述偏微分方程在时间上和空间上离散，并合并空间和时间节点相同的温度项，构建微分方程的离散格式。以式（7）为例，具体偏微分方程的离散格式为：

$\begin{split} \frac{{{\rho _{\text{m}}}Q{c_{\text{m}}}}}{{\Delta {z_j}}} &t_{{\text{0}},j - 1}^{n + 1} - \left( {\frac{{{\rho _{\text{m}}}Q{c_{\text{m}}}}}{{\Delta {z_j}}} + 2\text{π} {r_{\text{0}}}{h_{{\rm{c}}i}} + \frac{{\text{π} r_{\text{0}}^2{\rho _{\text{m}}}{c_{\text{m}}}}}{{\Delta t}}} \right) t_{{\text{0}},j}^{n + 1}+ \\ & 2\text{π} {r_{\text{0}}}{h_{{\rm{c}}i}} t_{{\text{1}},j}^{n + 1} = - {Q_{\text{p}}} - \frac{{\text{π} r_{\text{0}}^2{\rho _{\text{m}}}{c_{\text{m}}}}}{{\Delta t}} t_{{\text{0}},j}^n \\[-15pt] \end{split}$

(14)

式中：下标0，1分别代表钻柱内流体和钻柱壁。

将井筒–地层各区域的瞬态传热模型及边界条件均按照式（14）离散，构成多个微分方程组（如式（13）），然后应用全隐式有限差分法对方程组进行求解。

3. 井筒–地层传热解析模型

基于微分单元体能量平衡原理，建立了钻柱内和环空流体热量平衡方程（统称为井筒–地层传热解析模型），但建立模型时简化了地层和井筒间的热交换。建立后，对该解析模型进行了求解。

3.1 环空流体传热模型

从地层远处传递到井筒的总热量q_f可表示为^[10]：

${\text{ }}{q_{\text{f}}} = \frac{{2\text{π} {\lambda _{\text{f}}}}}{{{\rho _{\text{m}}}Q{t_{\text{D}}}}}({t_{\rm{i}}} - {t_{{\text{c1}}}}){\text{ }}$

(15)

式中： ${\lambda _{\text{f}}}$ 为地层导热系数，W/（m·℃）；t_i为初始地层温度，℃ ； ${t_{\text{D}}}$ 为无因次温度，计算方法可参考文献[11]。

根据傅里叶热传导定律，地层与环空流体交换的热量q_a可表示为^[12]：

${\text{ }}{q_{\text{a}}} = \frac{{2\text{π} {r_{\text{2}}}{U_{\text{a}}}}}{{{\rho _{\text{m}}}Q}}({t_{{\text{c1}}}} - {t_{\text{a}}})\Delta z{\text{ }}$

(16)

式中， ${U_{\text{a}}}$ 为钻井流体循环时环空内流体与地层之间的综合传热系数，W/（m·℃）。

联立式（15）和式（16）消除 ${t_{{\text{c1}}}}$ ，可得：

${\text{ }}{q_{\text{a}}} = \frac{{{c_{\text{m}}}}}{A}({t_{\rm{i}}} - {t_{\text{a}}})\Delta z$

(17)

$\,\text{其中} \qquad\qquad A = \frac{{{c_{\text{m}}}Q{\rho _{\text{m}}} ({\lambda _{\text{f}}} + {r_{\text{2}}}{U_{\text{a}}}{t_{\text{D}}})}}{{2\text{π} {r_{\text{2}}}{U_{\text{a}}}{\lambda _{\text{f}}}}}$

(18)

3.2 钻柱内流体传热模型

根据能量守恒原理，钻柱内流体与环空内流体交换的热量q_p可表示为：

${\text{ }}{q_{\text{p}}} = \frac{{{c_{\text{m}}}}}{B}({t_{\text{a}}} - {t_{\text{p}}})\Delta z$

(19)

$\,\text{其中} \qquad\qquad\qquad\quad B = \frac{{{c_{\text{m}}}Q{\rho _{\text{m}}}}}{{2\text{π} {r_{\text{0}}}{U_{\text{p}}}}} \qquad\quad\qquad$

(20)

式中， ${U_{\text{p}}}$ 为钻柱内流体与环空内流体间的综合传热系数，W/（m·℃）。

3.3 井筒–地层传热解析模型求解

基于钻柱内和环空流体微分单元能量守恒原理，对井筒–地层传热解析模型进行求解：

${t_{\text{p}}} = {\beta _{\text{1}}}{{\rm{e}}^{{k_{\text{1}}}z}} + {\beta _{\text{2}}}{{\rm{e}}^{{k_{\text{2}}}z}}{\text{ + }}{G_{\text{t}}}\left( {z - B} \right) + {t_{{\text{mf}}}}$

(21)

${t_{\text{a}}} = {\beta _{\text{1}}}(B{k_{\text{1}}} + 1){{\rm{e}}^{{k_{\text{1}}}z}} + {\beta _{\text{2}}}(B{k_{\text{2}}} + 1){{\rm{e}}^{{k_{\text{2}}}z}} + {G_{\text{t}}}z + {t_{{\text{mf}}}}$

(22)

$\,\text{其中} \qquad\qquad\qquad {k_{\text{1}}} = \frac{{1 + \sqrt {\frac{{4A}}B + 1} }}{{2A}}$

(23)

${k_{\text{2}}} = \frac{{1 - \sqrt {\frac{{4A}}B + 1} }}{{2A}}$

(24)

${\beta _{\text{1}}} = \varphi - \frac{{B{G_{\text{t}}} + \varphi B{k_{\text{1}}}{{\rm{e}}^{{k_{\text{1}}}H}} - {t_{\text{b}}}}}{{B({k_{\text{1}}}{{\rm{e}}^{{k_{\text{1}}}H}} - {k_{\text{2}}}{{\rm{e}}^{{k_{\text{2}}}H}})}}$

(25)

${\beta _{\text{2}}} = \frac{{B{G_{\text{t}}} + \varphi B{k_{\text{1}}}{{\rm{e}}^{{k_{\text{1}}}H}} - {t_{\text{b}}}}}{{B({k_{\text{1}}}{{\rm{e}}^{{k_{\text{1}}}H}} - {k_{\text{2}}}{{\rm{e}}^{{k_{\text{2}}}H}})}}{\text{ }}$

(26)

$\varphi ={t}_{\text{in}}－{t}_{\text{mf}}+B{G}_{\text{t}}$

(27)

${t_{{\text{mf}}}} = A{t_{{\text{fa}}}} + \left( {A + B} \right){t_{{\text{fp}}}} + {t_{{\text{in}}}} + {t_{\text{s}}}$

(28)

式中： ${G_{\text{t}}}$ 为地温梯度，℃/100m；H为井深，m； ${t_{\text{b}}}$ 为钻头压降产生的温度，℃； ${t_{{\text{fp}}}}$ 和 ${t_{{\text{fa}}}}$ 分别为钻柱内和环空内流体流动摩阻产生的温度，℃； ${t_{{\text{in}}}}$ 和 ${t_{\text{s}}}$ 分别为入口温度和地表温度，℃。

4. 算例分析

以顺北油田某超深井为例，利用该井井身结构与钻井参数进行了计算分析。

该井用ϕ444.5 、ϕ311.1和ϕ215.9 mm钻头分别钻至井深800，3 065及6 348 m，下入ϕ339.7、ϕ250.8及ϕ177.8 mm套管。在井深6 099 m处造斜，井眼曲率为7.46°/30 m，靶点A和靶点B对应井深分别为6 353和6 759 m，井斜角为30.35°。6 348～6 750 m井段采用ϕ149.2 mm钻头钻进，钻进中钻井液排量为16 L/s，密度为1.35 kg/L，塑性黏度为22 mPa·s，动切力为8 Pa，比热容为2 600 J/（kg·℃）、导热系数为0.32 W/（m·℃）。钻具主要由ϕ88.9 mm加重钻杆和ϕ101.6 mm钻杆组成。

以井深6 750 m处钻井参数为参考，地表温度与流体入口温度分别为10和30 ℃，地温梯度为2.467 ℃/100m，利用井筒–地层传热数值模型和解析模型计算钻柱内和环空流体的温度分布，结果见图2和图3。

图 2 数值模型计算的钻柱内与环空流体的温度分布

Figure 2. Fluid temperature distributions in the drill string and annulus calculated by numerical model

下载: 全尺寸图片幻灯片

图 3 解析模型计算的钻柱内与环空流体的温度分布

Figure 3. Fluid temperature distributions in the drill string and annulus calculated by analytical model

下载: 全尺寸图片幻灯片

从图2和图3可以看出，2种模型计算出的环空流体温度均高于钻柱内流体温度，表明钻井液从钻柱内向下流动过程中不断被环空钻井液加热，距井底越近，温度越高。

对比图2和图3可知，环空流体和钻柱内流体的温度差，数值解小于解析解，2种模型计算结果的差值分别为5 ℃和13 ℃左右。导致该结果的主要原因是：建立解析模型时忽略了钻柱内外壁对热对流交换的影响，采用综合传热系数表示环空流体与钻柱内流体的热交换（如式（20）中U_p），钻进中，钻柱内外钻井液与壁面进行强烈的对流换热，交换的热量远高于热传导方式下的热量；若仅用综合传热系数U_p来表述径向上钻柱内外壁面的热交换，减弱了热交换效率，使环空流体热量不能有效传递至钻柱内。因此，解析模型下的环空流体温度要远高于钻柱内流体温度。

未钻进时，井筒–地层达到热力学平衡状态，钻柱内与环空流体的温度与初始地层温度接近。钻进中井筒–地层为耦合体，流体循环不断将下部井段地层热量携带到地面，使下部井段流体温度低于初始地层温度，而上部井段高于初始地层温度。因此，井筒流体与近井壁地层温度间互相影响，即井筒–地层传热为瞬态过程。将图2和图3中温度分布与初始地层温度分布进行对比可知，解析模型得到的钻柱内与环空流体的温度均远低于初始地层温度。导致该现象的原因是：建立解析模型时对瞬态传热进行了简化，采用无因次温度t_D表示地层-井筒间热量传递的拟稳态过程，见式（15）。同时，用综合传热系数U_a表示流体循环时井壁与环空流体间的热交换效率，见式（16）。2种因素的共同作用，使解析模型计算出的钻柱内与环空流体的温度低于数值模型。

根据MWD实测数据，钻进10 h后环空温度为156.5 ℃。基于此，利用数值模型和解析模型计算环空流体温度，结果见图4所示。

图 4 数值模型和解析模型计算的环空流体温度

Figure 4. Fluid temperature in the annulus calculated bynumericaland analytical models

下载: 全尺寸图片幻灯片

从图4可以看出，数值模型和解析模型计算出的井底温度分别为158.79和145.64 ℃，两者相差13.15 ℃，与实测温度相比，2种模型的计算误差分别为1.46%和6.94%，表明数值模型的计算精度远高于解析模型。因此，进行钻进中井筒温度预测与评价研究时，推荐使用数值模型。

为了解近井壁地层温度扰动特征，利用数值模型计算了近井壁地层温度与环空流体温度之差的分布情况（见图5），分析流体循环对近井壁温度的影响。

图 5 近井壁温度与环空流体温度差分布

Figure 5. Distribution of the difference between the temperature near well walls and that of annular fluid

下载: 全尺寸图片幻灯片

从图5可以看出：井深5 100 m以深，近井壁地层温度高于环空流体温度，表明循环钻井液时携带了下部井段近井壁地层热量；井深5 100 m以浅，近井壁地层温度低于环空流体温度，表明钻井液在环空向上流动过程中携带下部地层热量并不断加热上部地层，使环空流体温度高于初始地层温度（见图2）。同时，距井壁越近，环空流体温度与地层温度的差越小；反之，距井壁越远，环空流体温度与地层温度的差越大。距井壁1.07和2.84 m处的地层温度与环空流体温度之差几乎相等，表明2种距离的地层温度也相等，也说明钻进中钻井液循环对近井壁地层的扰动距离在1.07 m左右。井深6 348 m处，近井壁地层温度与环空流体温度之差发生波动，由于技术套管下至该处，套管与地层的导热系数（分别为43.2和2.05 W/（m·℃））不同引起该现象。因此，计算井筒–地层瞬态温度时，需要考虑井身结构对近井壁地层传热效率的影响，使计算结果更能深入揭示流体循环条件下井筒–地层各控制区域的传热机理与温度演变特性。

5. 结　论

1）基于微分单元能量平衡原理，建立了钻进流体循环时的井筒–地层传热数值模型和井筒–地层传热解析模型，分别应用全隐式有限差分法和解析法对2种模型进行了求解。计算结果表明，下部井段环空流体温度低于初始地层温度，上部井段环空流体温度高于初始地层温度。

2）解析模型应用无因次时间函数简化了远处地层传入近井壁的热交换方式，采用综合传热系数表征近井壁地层传入环空的热量。同时，解析模型应用综合传热系数表示环空流体与钻柱内流体间的热交换，忽略了钻柱内外壁热对流的热交换作用，导致解析模型计算出的钻柱内与环空流体的温度低于数值模型。

3）近井壁地层温度数值模型能有效解释钻进中环空流体与近井壁温度演变特性，结合实测结果可知，数值模型和解析模型的计算误差分别为1.46%和6.94%，两者计算结果之差为13.15 ℃。推荐采用数值模型进行井筒瞬态温度预测。

图 1 泸203井区某平台井井眼轨道优化前后对比

Figure 1. Comparison of wellbore trajectory before and after optimization of a platform well in Well Block Lu 203

下载: 全尺寸图片幻灯片

表 1 泸州区块深层页岩气水平井各开次中完井深、层位优化前后对比

Table 1 Comparisons of completion depth and layer before and after optimization of different sections in deep shale gas horizontal wells in Luzhou Block

开次	井眼尺寸/ mm	套管尺寸/ mm	调整情况	中完井深/ m	中完层位
一开	660.4	508.0	优化前	50	沙溪庙组
一开	660.4	508.0	优化后	150	沙溪庙组
二开	444.5	339.7	优化前	1 100	须家河组顶
二开	444.5	339.7	优化后	1 200	须家河组中部
三开	311.1	250.8	优化前	2 900	韩家店组顶
三开	311.1	250.8	优化后	2 800	栖霞组顶

下载: 导出CSV

表 2 泸州区块2019—2021年钻井指标变化情况

Table 2 Variations of drilling index in Luzhou Block from 2019 to 2021

井眼直径/ mm	年份	钻井周期/d	机械钻速/ （m·h⁻¹）	平均钻进趟数
444.5	2019	18.4	7.28	4.4
	2020	17.6	7.52	4.0
	2021	12.5	10.21	3.5
311.1	2019	35.4	5.79	6.5
	2020	38.8	5.22	6.3
	2021	29.9	6.63	5.2
215.9	2019	59.1	4.72	14.0
	2020	58.8	5.31	10.3
	2021	49.5	6.53	8.7

下载: 导出CSV

参考文献(18)

[1]	冉启华,赵晗. 威远页岩气水平井钻井关键技术及发展方向[J]. 钻采工艺,2020,43(4):12–15. doi: 10.3969/J.ISSN.1006-768X.2020.04.04 RAN Qihua, ZHAO Han. Key technology and development direction of horizontal well drilling in Weiyuan shale gas reservoir[J]. Drilling & Production Technology, 2020, 43(4): 12–15. doi: 10.3969/J.ISSN.1006-768X.2020.04.04
[2]	樊好福,臧艳彬,张金成,等. 深层页岩气钻井技术难点与对策[J]. 钻采工艺,2019,42(3):20–23. doi: 10.3969/J.ISSN.1006-768X.2019.03.06 FAN Haofu, ZANG Yanbin, ZHANG Jincheng, et al. Technical difficulties and countermeasures of deep shale gas drilling[J]. Drilling & Production Technology, 2019, 42(3): 20–23. doi: 10.3969/J.ISSN.1006-768X.2019.03.06
[3]	刘厚彬,崔帅,孟英峰,等. 深层脆性页岩水平井井壁崩落失稳研究[J]. 断块油气田,2021,28(3):323–328. LIU Houbin, CUI Shuai, MENG Yingfeng, et al. Study on wellbore caving and instability of horizontal well in deep brittle shale[J]. Fault-Block Oil & Gas Field, 2021, 28(3): 323–328.
[4]	宋保健,孙凯,乐守群,等. 涪陵页岩气田钻井提速难点与对策分析[J]. 钻采工艺,2019,42(4):9–12. doi: 10.3969/J.ISSN.1006-768X.2019.04.03 SONG Baojian, SUN Kai, LE Shouqun, et al. Drilling acceleration challenges at Fuling Shale Gas Field and solutions[J]. Drilling & Production Technology, 2019, 42(4): 9–12. doi: 10.3969/J.ISSN.1006-768X.2019.04.03
[5]	石崇东,王万庆,史配铭,等. 盐池区块深层页岩气水平井钻井关键技术研究[J]. 石油钻探技术,2021,49(6):23–28. doi: 10.11911/syztjs.2021007 SHI Chongdong, WANG Wanqing, SHI Peiming, et al. Research on key drilling technology for horizontal wells in the deep shale gas reservoirs in Yanchi Block[J]. Petroleum Drilling Techniques, 2021, 49(6): 23–28. doi: 10.11911/syztjs.2021007
[6]	谭玮,房舟,罗成波,等. 致密高应力水平井优快钻完井难点与对策[J]. 断块油气田,2020,27(5):653–656. TAN Wei, FANG Zhou, LUO Chengbo, et al. Difficulties and solutions of optimized fast drilling and completion of horizontal well in tight and high stressed formation[J]. Fault-Block Oil & Gas Field, 2020, 27(5): 653–656.
[7]	席传明,史玉才,张楠,等. 吉木萨尔页岩油水平井JHW00421井钻完井关键技术[J]. 石油钻采工艺,2020,42(6):673–678. XI Chuanming, SHI Yucai, ZHANG Nan, et al. Key technologies for the drilling and completion of shale oil horizontal well JHW00421 in Jimusaer[J]. Oil Drilling & Production Technology, 2020, 42(6): 673–678.
[8]	赵洪波,朱迪斯,黄正,等. 南华北盆地亳州—阜阳地区页岩气钻井技术[J]. 石油钻采工艺,2020,42(6):679–683. ZHAO Hongbo, ZHU Disi, HUANG Zheng, et al. Shale gas drilling technologies used in Bozhou-Fuyang area of the southern north China Basin[J]. Oil Drilling & Production Technology, 2020, 42(6): 679–683.
[9]	王建龙,冯冠雄,刘学松,等. 长宁页岩气超长水平段水平井钻井完井关键技术[J]. 石油钻探技术,2020,48(5):9–14. doi: 10.11911/syztjs.2020086 WANG Jianlong, FENG Guanxiong, LIU Xuesong, et al. Key technology for drilling and completion of shale gas horizontal wells with ultra-long horizontal sections in Changning Block[J]. Petroleum Drilling Techniques, 2020, 48(5): 9–14. doi: 10.11911/syztjs.2020086
[10]	李传武,兰凯,杜小松,等. 川南页岩气水平井钻井技术难点与对策[J]. 石油钻探技术,2020,48(3):16–21. doi: 10.11911/syztjs.2020055 LI Chuanwu, LAN Kai, DU Xiaosong, et al. Difficulties and countermeasures in horizontal well drilling for shale gas in southern Sichuan[J]. Petroleum Drilling Techniques, 2020, 48(3): 16–21. doi: 10.11911/syztjs.2020055
[11]	刘伟,毛莉君,欧阳伟. 页岩气油基交联固化堵漏技术研究与应用[J]. 钻井液与完井液,2021,38(2):207–211. doi: 10.3969/j.issn.1001-5620.2021.02.013 LIU Wei, MAO Lijun, OUYANG Wei. Research on plugging technology of shale gas oil-based crosslinking solidification[J]. Drilling Fluid & Completion Fluid, 2021, 38(2): 207–211. doi: 10.3969/j.issn.1001-5620.2021.02.013
[12]	郑述权,谢祥锋,罗良仪,等. 四川盆地深层页岩气水平井优快钻井技术:以泸203井为例[J]. 天然气工业,2019,39(7):88–93. doi: 10.3787/j.issn.1000-0976.2019.07.011 ZHENG Shuquan, XIE Xiangfeng, LUO Liangyi, et al. Fast and efficient drilling technologies for deep shale gas horizontal wells in the Sichuan Basin: A case study of Well Lu 203[J]. Natural Gas Industry, 2019, 39(7): 88–93. doi: 10.3787/j.issn.1000-0976.2019.07.011
[13]	王建龙,于志强,苑卓,等. 四川盆地泸州区块深层页岩气水平井钻井关键技术[J]. 石油钻探技术,2021,49(6):17–22. doi: 10.11911/syztjs.2021062 WANG Jianlong, YU Zhiqiang, YUAN Zhuo, et al. Key technologies for deep shale gas horizontal well drilling in Luzhou Block of Sichuan Basin[J]. Petroleum Drilling Techniques, 2021, 49(6): 17–22. doi: 10.11911/syztjs.2021062
[14]	田逢军,王运功,唐斌,等. 长庆油田陇东地区页岩油大偏移距三维水平井钻井技术[J]. 石油钻探技术,2021,49(4):34–38. doi: 10.11911/syztjs.2021079 TIAN Fengjun, WANG Yungong, TANG Bin, et al. Drilling technology for long-offset 3D horizontal shale oil wells in the Longdong area of the Changqing Oilfield[J]. Petroleum Drilling Techniques, 2021, 49(4): 34–38. doi: 10.11911/syztjs.2021079
[15]	王涛,刘锋报,罗威,等. 塔里木油田防漏堵漏技术进展与发展建议[J]. 石油钻探技术,2021,49(1):28–33. doi: 10.11911/syztjs.2020080 WANG Tao, LIU Fengbao, LUO Wei, et al. The technical advance and development suggestions for leakage prevention and plugging technologies in the Tarim Oilfield[J]. Petroleum Drilling Techniques, 2021, 49(1): 28–33. doi: 10.11911/syztjs.2020080
[16]	柳伟荣,倪华峰,王学枫,等. 长庆油田陇东地区页岩油超长水平段水平井钻井技术[J]. 石油钻探技术,2020,48(1):9–14. doi: 10.11911/syztjs.2020029 LIU Weirong, NI Huafeng, WANG Xuefeng, et al. Shale oil horizontal drilling technology with super-long horizontal laterals in the Longdong region of the Changqing Oilfield[J]. Petroleum Drilling Techniques, 2020, 48(1): 9–14. doi: 10.11911/syztjs.2020029
[17]	孙永兴,贾利春. 国内 3000 m长水平段水平井钻井实例与认识[J]. 石油钻采工艺,2020,42(4):393–401. SUN Yongxing, JIA Lichun. Cases and understandings on the drilling of horizontal well with horizontal section of 3000 m long in China[J]. Oil Drilling & Production Technology, 2020, 42(4): 393–401.
[18]	光新军,叶海超,蒋海军. 北美页岩油气长水平段水平井钻井实践与启示[J]. 石油钻采工艺,2021,43(1):1–6. doi: 10.13639/j.odpt.2021.01.001 GUANG Xinjun, YE Haichao, JIANG Haijun. Drilling practice of shale oil & gas horizontal wells with long horizontal section in the North America and its enlightenment[J]. Oil Drilling & Production Technology, 2021, 43(1): 1–6. doi: 10.13639/j.odpt.2021.01.001

施引文献(19)

期刊类型引用(15)

1.	刘清友，赵建国，方世纪. 全液压连续油管牵引器工作机理研究. 钻采工艺. 2025(01): 216-227 . 百度学术
2.	曾凡辉，胡大淦，郭建春，张柟乔，郑彬涛，白小嵩，陈掌星. 基于钻录井资料的页岩气水平井分段多簇差异化压裂参数智能优化. 天然气工业. 2025(02): 84-94 . 百度学术
3.	唐弘程. 海上调整井优快钻井技术管理研究. 中国石油和化工标准与质量. 2024(08): 177-179 . 百度学术
4.	宋先知，姚学喆，许争鸣，周蒙蒙，王庆辰. 超深井控温钻井隔热涂层参数影响机制研究. 石油钻探技术. 2024(02): 126-135 . 本站查看
5.	陈建国，汪伟，都伟超. 渝西大安区块超深层页岩气水平井钻井提速关键技术研究. 钻探工程. 2024(04): 154-162 . 百度学术
6.	冯秀鲁. 油田深层油藏钻井提速技术研究及应用. 西部探矿工程. 2024(08): 68-70+74 . 百度学术
7.	程童，黎波，张中，黄兰，张庆，卢海兵，李晓晨，陈小龙，巩建平. 渝西区块深层页岩气井钻头优选研究. 四川地质学报. 2024(S1): 15-20 . 百度学术
8.	Yaoran Wei，Yongcun Feng，Zhenlai Tan，Tianyu Yang，Shuang Yan，Xiaorong Li，Jingen Deng. Simultaneously improving ROP and maintaining wellbore stability in shale gas well:A case study of Luzhou shale gas reservoirs. Rock Mechanics Bulletin. 2024(03): 91-103 . 必应学术
9.	欧翔，谭凯，周楚翔. 深层钻井堵漏材料的研究现状与发展思考. 材料导报. 2024(S2): 615-620 . 百度学术
10.	赵文彬，宋文豪. 永川黄202区块页岩气水平井钻井提速对策分析. 天然气勘探与开发. 2023(02): 127-134 . 百度学术
11.	温景东. 页岩气开采安全管理现状及改进策略. 石化技术. 2023(07): 218-220 . 百度学术
12.	张文平，许争鸣，吕泽昊，赵雯. 深层页岩欠平衡钻井气液固三相瞬态流动传热模型研究. 石油钻探技术. 2023(05): 96-105 . 本站查看
13.	倪维军，杨国昊，翟喜桐，马龙飞. 延安气田富县区域下古生界水平井优快钻井技术. 石油工业技术监督. 2023(12): 44-48 . 百度学术
14.	罗双平，刘青，韩巍，谭宇龙，周焱. 基于井筒工程数据的钻井提速评价分析方法——以高石梯—磨溪区块为例. 天然气技术与经济. 2023(06): 21-28 . 百度学术
15.	朱明明，孙欢，孙艳，丛成，侍德益，贾继国. 陇东致密油区域恶性出水漏层堵漏技术. 石油钻探技术. 2023(06): 50-56 . 本站查看

其他类型引用(4)

资源附件(0)

图(1) / 表(2)

计量

文章访问数: 576
HTML全文浏览量: 183
PDF下载量: 242
被引次数: 19

1. 井筒–地层传热物理模型
2. 井筒–地层传热数值模型
2.1 钻柱内流体传热模型
2.2 钻柱壁传热模型
2.3 环空流体传热模型
2.4 井壁传热模型
2.5 近井壁地层传热模型
2.6 井筒–地层传热数值模型求解
3. 井筒–地层传热解析模型
3.1 环空流体传热模型
3.2 钻柱内流体传热模型
3.3 井筒–地层传热解析模型求解
4. 算例分析
5. 结　论

1. 井筒–地层传热物理模型
2. 井筒–地层传热数值模型
2.1 钻柱内流体传热模型
2.2 钻柱壁传热模型
2.3 环空流体传热模型
2.4 井壁传热模型
2.5 近井壁地层传热模型
2.6 井筒–地层传热数值模型求解
3. 井筒–地层传热解析模型
3.1 环空流体传热模型
3.2 钻柱内流体传热模型
3.3 井筒–地层传热解析模型求解
4. 算例分析
5. 结　论

参考文献(18)

施引文献(19)

资源附件(0)

泸州区块深层页岩气水平井优快钻井技术

作者简介: 李涛（1993—），男，四川成都人，2015年毕业于西南石油大学石油工程专业，工程师，主要从事钻井完井技术研究与钻井故障及复杂管理工作。E-mail：litao01@petrochina.com.cn

计量

出版历程

Optimized and Fast Drilling Technology for Deep Shale Gas Horizontal Wells in Luzhou Block

1. 井筒–地层传热物理模型

2. 井筒–地层传热数值模型

2.1 钻柱内流体传热模型

2.2 钻柱壁传热模型

2.3 环空流体传热模型

2.4 井壁传热模型

2.5 近井壁地层传热模型

2.6 井筒–地层传热数值模型求解

3. 井筒–地层传热解析模型

3.1 环空流体传热模型

3.2 钻柱内流体传热模型

3.3 井筒–地层传热解析模型求解

4. 算例分析

5. 结 论

期刊类型引用(15)

其他类型引用(4)

计量

出版历程

目录

1. 井筒–地层传热物理模型

2. 井筒–地层传热数值模型

2.1 钻柱内流体传热模型

2.2 钻柱壁传热模型

2.3 环空流体传热模型

2.4 井壁传热模型

2.5 近井壁地层传热模型

2.6 井筒–地层传热数值模型求解

3. 井筒–地层传热解析模型

3.1 环空流体传热模型

3.2 钻柱内流体传热模型

3.3 井筒–地层传热解析模型求解

4. 算例分析

5. 结 论

作者简介:
李涛（1993—），男，四川成都人，2015年毕业于西南石油大学石油工程专业，工程师，主要从事钻井完井技术研究与钻井故障及复杂管理工作。E-mail：litao01@petrochina.com.cn

5. 结　论

5. 结　论