Research on Downhole Electromagnetic Repeater Transmission Characteristics in Ultra High Resistivity Gypsum-Salt Layers
-
摘要: 为解决超高阻盐膏层随钻电磁信号由于急速衰减易导致传输中断的问题，建立了超高阻盐膏层电磁中继转发模型，运用空间和阶次分别细化的自适应HP有限元算法，获得了超高阻地层环境下不同频率载波电磁中继信号衰减速度和传输距离的关系，形成了电磁中继器安装位置的预测方法，与常规有限元算法相比，该方法具有收敛速度和计算精度方面的优势。伊朗AGHA JARI区块AJ214井现场试验证明，应用所建模型计算的盐膏层电磁信号衰减速度与实际应用效果相符，中继器安装位置的预测准确率超过90%。研究表明，利用建立的传输模型可以准确预测中继器安装位置，充分发挥中继器转发能力，确保电磁信号可以在超高阻盐膏层连续传输，为电磁随钻测量技术的现场应用提供了理论参考。Abstract: In order to solve the problem that electromagnetic telemetry signal easily interrupts transmission in ultra high resistivity gypsum-salt layer due to fast attenuation,an adaptive HP finite element method (FEM),which could automatically select space or order thinning,was used to simulate the electromagnetic relay in ultra high resistivity gypsum-salt layers.The falling velocity and transmission distance of electromagnetic signals in ultra-high resistivity gypsum-salt layer under different carrier frequencies were obtained.A prediction method of relay installation position was formed,and compared with conventional finite element calculation.It had advantages in convergence rate and calculation accuracy.The the algorithm was verified correct in Well AJ214 of AGHA JARI Block,Iran.The calculation result of signal attenuation rate tallied well with the practical application,and the prediction accuracy of repeaters’ installation position was over 90%.The research showed that the transmission model could accurately predict the installation position of repeaters,and ensure the continuous transmission of electromagnetic signals across the ultra high resistivity gypsum-salt layer.The model can serve as a theoretical reference and foundation for EM-MWD applications.
-
随着我国油气资源勘探开发逐渐向深水、深地迈进，油气钻井过程中的压力体系愈来愈复杂，负压力窗口地层越来越多，如南海陵水区块、新疆顺北地区及准噶尔盆地等存在较多的负压力窗口地层，给油气钻井作业带来严重的技术挑战[1–4]。控压钻井技术是解决该问题的重要技术手段，但控压下套管前需制定下套管方案和确定补偿回压，而计算下套管准备阶段、下套管过程中及钻井液循环过程中的井筒压力是制定下套管方案和确定补偿回压的关键。
目前，国内外针对下套管过程中井筒压力的计算开展了大量研究，提出了稳态法以及瞬态法波动压力计算模型。20世纪70年代前，波动压力计算均是基于稳态法理论进行计算的。G. E. Cannon[5]最早对波动压力开展了研究；G. S. Ormsby[6]给出了不同流态下井内波动压力的计算方法；E. H. Clark[7]基于一维稳定流动建立了波动压力计算模型，并分析了钻井液排量、管柱运动速度对波动压力的影响程度；J. A. Burkhardt[8]利用现场井内压力实测数据，分析了井内波动压力的变化规律，认为井内波动压力是由钻井液的流动阻力、惯性力和静切力综合作用产生的；F. J. Schuh等人[9-11]先后提出了宾汉模式、幂律模式和Roberston-Stiff模式下的波动压力计算方法；王锦昌[12]认为建立波动压力模型时，要考虑稳定器安放位置对环空间隙的影响；彭齐等人[13]以槽流模型为基础，分别讨论了层流、紊流状态下的波动压力，建立了基于钻柱运动的稳态井筒波动压力计算模型；He Shiming等人[14]建立了基于回归方程的波动压力模型，模型的计算误差为±3%。以上研究主要分析了钻井液密度、稳定器数量、环空间隙及钻柱偏心度等对波动压力的影响。
20世纪70年代中期至今，随着计算机的发展，波动压力瞬态计算法得到了普及。A. Lubinski等人[15]考虑流体的压缩性和流道管壁的膨胀性及流体的惯性，建立了波动压力计算模型；M. Lal[16]预测了井眼内各点的最大波动压力和抽汲压力随时间的变化规律；Zhong Bing等人[17]提出用混合隐式特征线法求解不同流道瞬态波动压力的方法，得出了不同情况下瞬态波动压力的数值解；王超等人[18]根据不稳定流动理论，采用特征线求解方法，分析了压力波在井内的传递与转化过程；Zhang Feifei等人[19]提出了深水钻井近井筒瞬态波动压力、井筒压力与地应力场振荡耦合的井筒稳定性分析方法；吴鹏程等人[20]考虑液–固两相介质，采用数值模拟法分析了影响起下钻瞬态波动压力的因素；谭东[21]分析了井底压力与溢流的关系；屈俊波[22]引入动态摩阻，分析了钻具加速度对波动压力的影响规律。然而，以上波动压力计算模型均是基于井筒内流体为均质流体建立的，无法实现负压力窗口地层的防漏压稳。为此，笔者以井筒内多密度梯度钻井液为基础，建立了下套管过程中的液柱结构及波动压力计算模型，以乐探1井为例，模拟了起钻后、下套管过程中和钻井液循环等3个阶段井筒内的液柱结构及压力分布情况，形成了负压力窗口油气井控压下套管井筒压力控制技术。
1. 负压力窗口下套管技术难点
负压力窗口下套管过程中的复杂性及主要技术难点包括以下几点：
1）采用单一密度的钻井液无法实现负压力窗口的防漏压稳，需采用多密度梯度的钻井液液柱结构。
2）下套管以及钻井液循环过程中，井筒内多密度梯度钻井液液柱结构的变化复杂。
3）井筒内液柱的结构在下套管过程中时刻变化，下套管过程中的井筒压力计算模型复杂，难以确定合理的套管下入速度。
综上所述，下套管过程中井筒内液柱结构及井筒压力计算是解决负压力窗口下套管技术难点的关键。
2. 下套管井筒内液柱结构计算模型
2.1 完钻后井筒内液柱结构的计算模型
对于负压力窗口井段，由于起钻后井筒内无法施加回压，因此采用单一密度的钻井液无法实现负压力窗口地层的防漏压稳，需要在起钻过程中根据压力分布情况循环不同密度的钻井液。
假设某裸眼井段有n套压力体系，对应孔隙压力及漏失压力的当量密度分别为（ρp1，ρL1）、（ρp2，ρL2）、（ρp3，ρL3）、……、（ρpn，ρLn），所在井深分别为h1′，h2′，h3′，…，hn′，完钻时井筒内钻井液密度为ρ0，如图1（a）所示。由于井底漏失压力当量密度ρL1小于上部孔隙压力当量密度，因此，采用由下而上的方法设计多密度梯度钻井液液柱结构：
1）起钻前，由于需要用加重钻井液压稳上部地层，因此，为保证不压漏井底地层，需要用密度为ρ1的低密度钻井液替换部分井底的钻井液，如图1（b）所示；
2）起钻时，钻杆上提至一定高度后，密度为ρ1低密度钻井液的高度变为h1，为压稳上部地层，循环密度为ρ2的加重钻井液，如图1（c）所示；
3）同理，根据不同的压力体系，设计上提钻杆高度和循环密度为ρi的加重钻井液，如图1（c）所示；
4）钻杆全部起出后，无法施加回压，井筒内多密度梯度钻井液的静液柱压力起到防漏压稳的作用，设计的多密度梯度钻井液密度与使用长度分别为（ρ1，h1）、（ρ2，h2）、（ρ3，h3）、……、（ρm， hm），如图1（d）所示。
设计密度为ρi钻井液与使用长度hi时，需满足以下条件：
ppj⩽ (1) 式中：k为设计的第k种密度的钻井液；ρi为第i段钻井液的密度，kg/cm3；hi为第i段钻井液的长度，m；ppj为第j段的孔隙压力，MPa；pLj为第j段的漏失压力，MPa。
2.2 下套管过程中井筒内液柱结构的计算模型
随着套管下入，套管排开井内液体的体积不断变化，加之套管下至一定深度后，需要循环钻井液，对于井身结构复杂的油气井，下套管过程中井筒内液柱结构的计算尤其困难，为了实现下套管过程井筒内液柱结构的准确计算，采用如下的计算方法。
1）假设井筒内多密度梯度钻井液体积由下到上分别为V0，V1，V2，……，Vm，送入钻杆与套管的长度分别为LD和LC，两者之和为井深，下套管过程分为s个阶段，每个阶段的套管下入速度分别为v1，v2，v3，…，vs，对应的下入井深分别为H1，H2，H3，…，Hs，如图2所示。
2）计算任一时刻下入管柱在井内的体积Vp，与井筒内k～m种密度钻井液体积之和Vf比较，若Vp>Vf，则第m，m−1，m−2，……，k种密度钻井液被替换出井筒。
{V}_{\mathrm{f}}={\sum }_{i=k}^{m}{V}_{i} (2) 式中：Vi为第i种钻井液的体积，m3；Vf井筒内 k～m种密度钻井液体积之和，m3。
3）计算此时在井筒内的第0，1，2，……，k−1种密度钻井液的长度，分套管引鞋以上及以下2部分分别计算各段液柱的结构，套管引鞋以上环空中液柱结构的计算流程如图3所示。
![]()
图 3 下套管过程中环空液柱结构的计算流程Figure 3. Calculation flow chart of annular slurry column structure while running casing4）管柱在第n个阶段下到井深Hn后，循环密度为ρcp的钻井液，循环结束后，继续下套管至第n+1个阶段的井深Hn+1，重复过程2）～3），即可求出下套管过程中井筒内的液柱结构。
3. 下套管井筒压力计算模型
3.1 假设条件
建立下套管井筒压力计算模型时，进行以下假设：1）钻井液的流动均视为一维不稳定流动；2）认为井筒中流体的流动为绝热流动，与地层无热量交换；3）认为井筒内不同密度钻井液之间为等截面接触。
3.2 模型描述
负压力窗口下套管时，以套管引鞋为分界线，将井筒分为上下2部分计算其压力。根据下套管作业实际工况，分为2个阶段建立下套管过程中的井筒压力计算模型：1）控压下套管井筒压力计算模型；2）循环钻井液过程井筒压力计算模型。
3.2.1 控压下套管井筒压力计算模型
1）顶替钻杆未入井时，尾管无法实现控压，此时的井筒压力p1等于套管引鞋上端静液柱压力phu、套管引鞋下端静液柱压力phd、波动压力po之和。可以根据下式计算出该段地层安全密度窗口下的套管下入速度。
{p_1} = {p_{\text{o}}}({v_{\text{c}}}) + {p_{\text{hu}}} + {p_{\text{hd}}} < {p_{\text{L}}} (3) 2）顶替钻杆入井后，可以通过施加回压压稳地层。此时套管下入过程中的井筒压力p1包括4部分：套管引鞋下端的静液柱压力phd、套管引鞋上端的静液柱压力phu、管柱下入引起的波动压力po、井口施加的回压pb。为了保证套管能安全下入到井内，要求严格控制井筒压力在安全密度窗口范围内，物理模型如图4（a）所示。
{p_1} = {p_{\text{o}}}({v_{\text{c}}}) + {p_{\text{hu}}} + {p_{\text{hd}}} + {p_{\text{b}}} < {p_{\text{L}}} (4) 式中：p1为下套管过程中的井筒压力，MPa；phu为套管引鞋上端静液柱压力，MPa；phd为套管引鞋下端的静液柱压力，MPa；po为波动压力，MPa；vc为套管下入速度，m/s；pb为井口施加的回压，MPa；pL为漏失压力，MPa。
通过计算套管引鞋上下端的静液柱压力、波动压力，进而找到井筒压力与套管下入速度的关系，并根据安全密度窗口确定套管合理的下入速度与环空回压。
3.2.2 循环钻井液过程井筒压力计算模型
循环钻井液（灌浆）过程中的井筒压力p2包括井筒内钻井液的静液柱压力和灌浆过程流体上返时的摩阻压降pf2部分，物理模型如图4（b）所示。
{p_2} = {p_{\text{hu}}} + {p_{\text{hd}}} + {p_{\text{f}}} < {p_{\text{L}}} (5) 式中：p2为循环钻井液过程中的井筒压力，MPa；pf为钻井液循环摩阻压降，MPa。
3.3 模型求解
3.3.1 静液柱压力
参考2.2节中关于下套管过程中井筒内液柱结构的计算方法，可以求得井筒内静液柱压力ph。
{p_{\text{h}}} = \sum\limits_{n = 0}^{k - 1} {{\rho _n}g} {h_n} + \sum\limits_{l = k - 1}^m {{\rho _l}g} {h_l} (6) 式中：ph为井筒内静液柱压力，MPa；ρn为套管引鞋上部环空第n段钻井液的密度，kg/cm3；hn为环空第n段液柱长度，m；ρl为套管引鞋下部第l段钻井液的密度，kg/cm3；hl为套管引鞋下部第l段液柱长度，m。
3.3.2 摩阻压降
循环钻井液时，钻井液由套管注入，从环空返出。以幂律流体为例，环空中钻井液上返时的摩阻压降为：
p_{\text{f}} = \frac{{4KL}}{{{D_{{\text{hy}}}}}}{\left(\frac{{2n + 1}}{{3n}}\frac{{12v_\text{r}}}{{{D_{{\text{hy}}}}}}\right)^n} (7) v_{\text{r}}= \frac{{4Q_{\text{r}}}}{{{\text{π }}(D_{\text{h}}^{\text{2}} - D_{{\text{co}}}^{\text{2}})}} (8) 式中：K为返出钻井液的稠度系数，Pa·sn；L为返出钻井液上返距离，m；n为返出钻井液的流性指数；Dhy为环空水力直径，m；vr为返出钻井液平均流速，m/s；Dh为环空直径，m；Dco为套管外径，m；Qr为返出钻井液的排量，L/s。
3.3.3 波动压力
基于井内不稳定流动理论，利用波动压力瞬态计算法计算波动压力时只考虑传递到井口和井底处的压力变化。
1）控制方程。根据质量守恒和动量守恒定律，得到井筒内一维不稳定流动的控制方程：
\left\{ \begin{gathered} \frac{Q}{S}\frac{{\partial p }}{{\partial z}} + \frac{{\partial p}}{{\partial t}} + \frac{{\rho {c^2}}}{S}\frac{{\partial Q}}{{\partial z}} = 0 \\ \frac{{\partial p}}{{\partial z}} + \frac{\rho }{S}\frac{{\partial Q}}{{\partial t}} + \frac{{\rho Q}}{{{S^2}}}\frac{{\partial Q}}{{\partial z}} + {p_{\text{f}}} = 0 \\ \end{gathered} \right. (9) \, 其中\qquad\qquad\quad \;c= \frac{1}{{\sqrt {\rho (\alpha + \beta )} }} \qquad\quad (10) 式中：Q为钻井液流量，L/s；c为压力波传播速度，m/s；S为流道截面积，m2；t为套管下放时间，s；z为套管长度，m；β为流道膨胀系数；α为钻井液压缩系数；ρ为钻井液密度，kg/cm3。
2）方程求解。采用特征线法求解控制方程组，首先根据修正的Lister方法建立控制方程组式（9）的特征方程，然后利用有限差分法求解该特征方程，进而得到控制方程组的解。控制方程组式（9）的特征方程为：
\left\{ \begin{gathered} \pm {\text{d}}p + \frac{{\rho c}}{S}{\text{d}}Q + c{p_{\text{f}}}{\text{d}}t = 0 \\[-2pt] \frac{{{\text{d}}z}}{{{\text{d}}t}} = v_\text{r} \pm c \\ \end{gathered} \right. (11) 根据式（11）的第2个方程，分别绘制斜率为arctan
\left( \dfrac{1}{v_\text{r}+c} \right) 且与空间坐标轴相交于A点的前向特征线及斜率为arctan\left( \dfrac{1}{v_\text{r}-c}\right) 且与空间坐标轴相交于B点的后向特征线，2条特征线的交点为点C，如图5所示。当选择的空间步长∆z和时间步长∆t足够小时，可以把特征线段AC和BC看作直线时，则AC和BC线上的点满足以下条件：\left\{ \begin{gathered} {\text{d}}p + \frac{{\rho c}}{S}{\text{d}}Q + c{p_{\text{f}}}{\text{d}}t = 0 \\[-1pt] - {\text{d}}p + \frac{{\rho c}}{S}{\text{d}}Q + c{p_{\text{f}}}{\text{d}}t = 0 \\ \end{gathered} \right. (12) 线段AC满足式（12）中的第1个方程，线段BC满足式（12）中的第2个方程，C点是特征线方程的交点，所以C点满足式（12）的2个方程。式（12）的有限差分格式可以写为：
\left\{ \begin{gathered} {p_C} - {p_A} + \frac{{\rho c}}{S}({Q_C} - {Q_A}) + c\Delta t{p_{{\text{f}}A}} = 0 \\ - {p_C} + {p_B} + \frac{{\rho c}}{S}({Q_C} - {Q_B}) + c\Delta t{p_{{\text{f}}B}} = 0 \\ \end{gathered} \right. (13) 式中：pA，pB和pC分别为A点、B点和C点处的瞬态波动压力，MPa；QA，QB和QC分别为A点、B点和C点处的钻井液流量，L/s；pfA和pfB分别为A点和B点处的摩阻压降，MPa。
若已知pA，QA，pfA，pB，QB和pfB，即可求得流道中任意点处在任意时刻的瞬态压力pC和流量QC。
{Q_C} = \frac{{{p_A} + {p_B} - \dfrac{{\rho c}}{S}({Q_B} - {Q_A}) - c\Delta t({p_{{\text{f}}A}} - {p_{{\text{f}}B}})}}{{\dfrac{{2\rho c}}{S}}} (14) {p_C} = \frac{{{p_A} + {p_B} - \dfrac{{\rho c}}{S}({Q_B} - {Q_A}) - c\Delta t({p_{{\text{f}}A}} - {p_{{\text{f}}B}})}}{2} (15) 3）定解条件。实际下套管作业中，管柱在井内的运动造成井内流体在2个流道中流动，即环空流道和井底裸眼流道。管柱运动引起的波动压力在套管引鞋处产生，由套管引鞋处分别向井口和井底衰减，环空流道的井口边界与大气相连通，所以环空流道井口边界压力为0，且2个流道在交汇点（即套管引鞋）处的压力相等。套管引鞋处的流体受到管柱的顶替作用，环空流道的流体向上流动，井底裸眼流道的流体向下流动，不考虑钻井液在井底的渗透，所以井底裸眼流道的边界流量为0，且2个流道在交汇点（即套管引鞋）处的流量之和为管柱顶替流量。
a. 环空流道的井口边界压力为0，流量由前向特征线方程求得。
\left\{ \begin{gathered} {p_C} = 0 \\ {Q_C} = \frac{{{p_A} + \dfrac{{\rho {c_A}}}{{{S_A}}}{Q_A} - {c_A}{p_{{\text{f}}A}}\Delta t}}{{\dfrac{{\rho {c_A}}}{{{S_A}}}}} \\ \end{gathered} \right. (16) b. 井底流道井底边界流量为0，压力由前向特征线方程求得。
\left\{ \begin{gathered} {Q_C} = 0 \\ {p_C} = {p_A} + \frac{{\rho {c_A}}}{{{S_A}}}{Q_A} - {c_A}{p_{{\text{f}}A}}\Delta t \\ \end{gathered} \right. (17) c. 交汇点即套管引鞋处管柱顶替流量为2个流道流量之和。
\left\{ \begin{gathered} {Q_\text{p}} = {Q_{C1}} + {Q_{C2}} \\ {p_C} = {p_\text{a}} + \frac{{\rho {c_\text{a}}}}{{{S_\text{a}}}}({Q_{C1}} - {Q_\text{a}}) + {c_\text{a}}{p_{{\text{fa}}}}\Delta t \\ {p_C} = {p_\text{d}} + \frac{{\rho {c_\text{d}}}}{{{S_\text{d}}}}({Q_{C2}} - {Q_\text{d}}) + {c_\text{d}}{p_{{\text{fd}}}}\Delta t \\ \end{gathered} \right. (18) 式中：cA，ca和cd分别为A点、环空流道和井底流道的压力波传播速度，m/s；SA，Sa和Sd分别为A点、环空和井底的流道截面积，m2；QC1和QC2分别为C点环空和井底流道的流量，m3/s；Qp为下入套管排开流量，m3/s；pfa和pfd分别为环空和井底的摩阻压降，MPa。
联立式（14）和式（15），代入定解条件，可得到下套管过程中井筒各位置处的压力和流量，具体求解流程如图6所示。根据图6中的波动压力求解流程，利用C#编程语言开发了波动压力计算程序，实现了负压力窗口下套管过程中波动压力的计算。
4. 计算实例
在上述下套管井筒液柱结构及井筒波动压力计算模型的基础上，以乐探1井ϕ139.7 mm油层尾管下入过程为研究对象，进行实例计算。乐探1井是新疆油田的一口风险探井，位于准噶尔盆地南缘冲断带霍玛吐背斜带吐谷鲁背斜，完钻井深7 050.00 m，该井存在4个漏层，分别在井深6 152.00，6 674.00，6 947.00和7 028.00 m处，井深6 152.00 m处的当量循环密度为2.445 kg/L，井深6 674.00 m处的当量循环密度为2.446 kg/L，井深6 947.00 m处的当量循环密度为2.432 kg/L，井深7 028.00 m处的当量循环密度为2.432 kg/L。
由控压钻井数据（见表1）可知，5 620.00～5 746.00和5 746.00～6 013.00 m井段的地层压力大于井深6 674.00，6 947.00和7 028.00 m等处的漏失压力，符合负压力窗口地层特点，是负压力窗口井的典型代表。
表 1 乐探1井裸眼各层位地层压力Table 1. Formation pressure of open-hole layers in Well Letan-1井段/m 层位 压力系数 实际钻井液密度/
（kg·L−1）5 620～5 746 K1h 2.49～2.54 2.39～2.41 5 746～6 013 K1h-K1q 2.45～2.48 2.34～2.40 6 013～6 160 K1q 2.42～2.45 2.55～2.59 6 160～6 808 K1q 2.40～2.42 2.56～2.57 6 808～7 050 J3k 2.34～2.39 2.34～2.36 通过模拟下套管过程中井筒内的液柱结构以及井筒压力变化情况，可以给出下套管过程钻井液循环方案，推荐不同阶段套管下入速度及控压值。
4.1 完钻后井筒内的液柱结构
针对负压力窗口地层，将井筒内的液柱设计成多密度梯度钻井液液柱结构，根据地层压力体系调整起钻过程中注入钻井液的密度，地层压力偏高井段采用重浆压稳，地层压力偏低井段采用轻浆或重浆与轻浆搭配使用，以实现防漏压稳。利用完钻后井筒内液柱结构设计模型，结合乐探1井裸眼段各层位地层压力体系，设计了乐探1井完钻后井筒内的液柱结构，结果如表2所示。
表 2 乐探1井完钻后井筒内的液柱结构Table 2. Slurry column structure in Well Letan-1 after drilling井深/m 钻井液密度/
（kg·L−1）段长/m 液柱压力/
MPa当量循环密度/
（kg·L−1）1 523.00 2.50 1 523 37.352 2.500 3 560.00 2.50 2 037 49.957 2.500 5 852.00 2.38 2 292 53.513 2.453 5 900.00 2.38 48 1.121 2.452 6 152.00 2.35 252 5.809 2.448 6 300.00 2.35 48 3.412 2.446 6 674.00 2.10 374 7.705 2.427 6 947.00 2.10 273 5.624 2.414 7 028.00 2.10 81 1.669 2.410 7 050.00 2.10 22 0.453 2.409 4.2 下套管过程中的钻井液循环方案
套管下至一定深度后，需要根据密度窗口调整井筒内钻井液液柱结构。基于起钻时井筒内形成的多密度梯度钻井液液柱结构，考虑套管下至一定深度时井筒内部分钻井液被排出井外引起井筒内液柱结构的变化情况，设计循环钻井液的密度、用量和排量。采用下套管过程井筒内液柱结构的计算方法，结合乐探1井的地层压力体系，分别设计了3种不同密度梯度的钻井液循环方案（见表3）。
表 3 乐探1井套管下入不同深度处的钻井液循环方案Table 3. Drilling fluid circulation schemes at different depths of casing running in Well Letan-1套管下深/m 钻井液循环方案1 钻井液循环方案2 钻井液循环方案3 3 560.00 2.43 kg/L+2.32 kg/L 2.45 kg/L+2.35 kg/L 2.47 kg/L+2.38 kg/L 5 900.00 2.32 kg/L 2.35 kg/L 2.38 kg/L 7 050.00 2.32 kg/L 2.30 kg/L 2.28 kg/L 利用3.1.2节中的井筒压力计算流程，计算3种钻井液循环方案下井筒当量循环密度的分布，结果见图7。由图7可以看出：采用钻井液循环方案1时，6 700.00～7 100.00 m井段的当量循环密度大于漏层的当量循环密度；采用钻井液循环方案3时，5 900.00～6 600.00 m井段的当量循环密度大于漏层的当量循环密度；采用钻井液循环方案2时，全井段的当量循环密度都处于安全密度窗口范围内，因此，选用循环钻井液方案2。对于该方案，若送入钻杆已进入井筒，可通过设计合理的排量和控压值实现井筒压力控制，以套管下至井深5 900.00 m时为例，通过计算得到密度2.35 kg/L的钻井液分别以0.1，0.2，0.3和0.4 m3/min 4种排量循环时的井口控压值分别为4.65，3.90，3.27和2.64 MPa。通过改变套管下入速度调节下套管过程中的井筒压力，可以实现负压力窗口地层套管的安全下入。不同井段套管下入速度的方案如表4所示。
![]()
图 7 不同循环钻井液方案下井筒当量循环密度的分布Figure 7. Equivalent circulating density distribution of wellbore under different drilling fluid circulation schemes表 4 不同井段的套管下入速度方案Table 4. Casing running speed schemes for different well sections井段m 方案A/（m·s−1） 方案B/（m·s−1） 方案C/（m·s−1） 0～1 523.00 0.160 0.160 0.160 1 523.00～3 560.00 0.160 0.160 0.160 3 560.00～5 900.00 0.145 0.145 0.145 5 900.00～6 674.00 0.137 0.137 0.137 6 674.00～7 050.00 0.124 0.110 0.137 利用3.2节的控压下套管井筒压力模型，计算分析了以上3种套管下入方案的井筒当量循环密度分布情况（见图8）。由图8可以看出：采用套管下入速度方案A和B时，下套管过程中井筒当量循环密度均处于安全密度窗口范围内；采用套管下入速度方案C时， 6 800.00～7 050.00 m井段的当量循环密度大于漏层的当量密度。考虑下套管效率，选用套管下入速度方案A。同时计算了该方案下的井筒波动压力和环空回压，结果如表5所示。
![]()
图 8 套管不同下放速度方案的ECD变化Figure 8. Equivalent circulating density changes under different casing running speed schemes表 5 采用套管下入速度方案A时不同井段的波动压力和环空回压Table 5. Fluctuating pressure and annular back pressure at different well sections under casing running speed scheme A井段/m 套管下入速度/
（m·s−1）波动压力/
MPa环空回压/
MPa0～1 523.00 0.160 0.32 1 523.00～3 560.00 0.160 0.77 1.641 3 560.00～5 900.00 0.145 1.74 3.427 5 900.00～6 674.00 0.137 1.70 4.041 6 674.00～7 050.00 0.124 1.80 4.457 乐探1井采用以上设计的套管下入速度方案（见表5）和循环钻井液方案2（见表3），进行控压下套管作业，套管顺利下至设计位置，没有出现井漏与溢流情况，为控压固井安全施工奠定了基础。
5. 结 论
1）以一维不稳定流动理论为基础，结合井筒内多密度梯度钻井液的液柱结构，建立了控压下套管过程中井筒内液柱结构以及井筒压力的计算模型，通过利用其计算套管下至不同井深处时的井筒压力分布，找到井筒压力与套管下放速度的关系，从而确定套管合理的下入速度。
2）根据提出的井筒压力计算方法，模拟计算了乐探1井下套管及循环钻井液过程中液柱结构的变化情况，推荐了下套管过程中钻井液循环方案和不同井段的套管下入速度及井口控压值。结果表明：采用多密度梯度钻井液液柱结构可以实现负压力窗口地层的防漏压稳；采用优化的钻井液循环方案和套管下入速度方案下套管，当量循环密度均在安全密度窗口范围之内，不会出现井漏与溢流情况。
3）建立的控压下套管过程中井筒内液柱结构以及井筒压力计算模型，为负压力窗口油气井下套管方案的制定提供了理论依据。
-
[1] 闫宏亮,石文龙,李琳.随钻测量信息传输方式的发展现状综述研究[J].重庆科技学院学报(自然科学版),2015,17(6):69-72,83. YAN Hongliang,SHI Wenlong,LI Lin.Research on the status and future of the MWD data transmission[J].Journal of Chongqing University of Science and Technology (Natural Sciences Edition),2015,17(6):69-72,83. [2] 袁鹏斌,余荣华,欧阳志英.无线随钻测量信息传输的现状与问题[J].焊管,2010,33(10):65-69. YUAN Pengbin,YU Ronghua,OUYANG Zhiying,et al.The current status and problems of wireless signal transmission in measurement while drilling[J].Welded Pipe and Tube,2010,33(10):65-69. [3] 马哲,杨锦舟,赵金海.无线随钻测量技术的应用现状与发展趋势[J].石油钻探技术,2007,35(6):112-115. MA Zhe,YANG Jinzhou,ZHAO Jinhai,et al.Status quo and development trend of MWD technique[J].Petroleum Drilling Techniques,2007,35(6):112-115. [4] 汤明,何世明,邢景宝,等.大牛地气田DP14水平井氮气泡沫钻井实践与认识[J].天然气工业,2010,30(3):74-76. TANG Ming,HE Shiming,XING Jingbao,et al.Practices and knowledge from nitrogen foam drilling at the Well DP-14 in the Daniudi Gas Field[J].Natural Gas Industry,2010,30(3):74-76. [5] 呼石磊,鄢泰宁,李晓.地层对电磁随钻测量信号的影响研究[J].煤炭科学技术,2011,39(9):114-117. HU Shilei,YAN Taining,LI Xiao.Study on strata affected to electromagnetic with drilling measuring signal[J].Coal Science and Technology,2011,39(9):114-117. [6] 刘修善,侯绪田,涂玉林,等.电磁随钻测量技术现状及发展趋势[J].石油钻探技术,2006,34(5):4-9. LIU Xiushan,HOU Xutian,TU Yulin,et al.Developments of electromagnetic measurement while drilling[J].Petroleum Drilling Techniques,2006,34(5):4-9. [7] 李林,弓志谦,王磊,等.扩展EM-MWD传输深度及提高可靠性的方法与对策[J].钻采工艺,2010,33(4):25-27,38. LI Lin,GONG Zhiqian,WANG Lei,et al.Methods of improving EM-MWD transmissive depth and reliability[J].Drilling Production Technology,2010,33(4):25-27,38. [8] 邵养涛,姚爱国,张明光.电磁波随钻遥测技术在钻井中的应用与发展[J].煤田地质与勘探,2007,35(3):77-80. SHAO Yangtao,YAO Aiguo,ZHANG Mingguang.Application and development of electro-magnetic telemetry in drilling operation[J].Coal Geology Exploration,2007,35(3):77-80. [9] 胡长翠,张明友,张琴,等.井下测试数据无线传输技术探讨[J].钻采工艺,2011,34(1):48-51. HU Changcui,ZHANG Mingyou,ZHANG Qin,et al.Research on wireless telemetry technology of downhole test data[J].Drilling Production Technology,2011,34(1):48-51. [10] 胡越发,杨春国,高炳堂.井下电磁中继传输技术研究及应用[J].科技导报,2015,33(15):66-71. HU Yuefa,YANG Chunguo,GAO Bingtang.Research and application of downhole electromagnetic relay transmission technology[J].Science Technology Review,2015,33(15):66-71. [11] 郑俊华,宗艳波,钱德儒,等.井下环境模拟试验装置研制[J].科学技术与工程,2017,17(8):21-25. ZHENG Junhua,ZONG Yanbo,QIAN Deru,et al.Development of downhole environment simulation test equipment[J].Science Technology and Engineering,2017,17(8):21-25. [12] 李勇华,杨锦舟,杨震,等.随钻电阻率地层边界响应特征分析及应用[J].石油钻探技术,2016,44(6):111-116. LI Yonghua,YANG Jinzhou,YANG Zhen,et al.The analysis and application of formation interface response characteristics of the resistivity LWD tool[J].Petroleum Drilling Techniques,2016,44(6):111-116. [13] 李翠,高德利,刘庆龙,等.邻井随钻电磁测距防碰计算方法研究[J].石油钻探技术,2016,44(5):52-59. LI Cui,GAO Deli,LIU Qinglong,et al.A method of calculating of avoiding collisions with adjacent wells using electromagnetic ranging surveying while drilling tools[J].Petroleum Drilling Techniques,2016,44(5):52-59. [14] 解茜草,赵志峰.随钻电磁波测井响应时域有限差分数值模拟[J].计量与测试技术,2014,41(8):1-3. XIE Xicao,ZHAO Zhifeng.A study on modeling of electromagnetic wave propagation resistivity logging tool based on the FDTD method[J].Metrology Measurement Technique,2014,41(8):1-3. [15] 陈晓晖,刘得军,刘得芳,等.随钻电阻率测量系统响应模拟算法探讨[J].电子测量技术,2009,32(6):21-26. CHEN Xiaohui,LIU Dejun,LIU Defang,et al.Discussion on methods of simulation for resistivity logging-while-drilling (LWD) measurements[J].Electronic Measurement Technology,2009,32(6):21-26. [16] 胡松,王晓畅,孔强夫.水平井随钻电磁波电阻率数值模拟[J].科学技术与工程,2017,17(14):59-66. HU Song,WANG Xiaochang,KONG Qiangfu.Numerical simulation of LWD resistivity in horizontal wells[J].Science Technology and Engineering,2017,17(14):59-66. [17] 杨震,肖红兵,李翠.随钻方位电磁波仪器测量精度对电阻率及界面预测影响分析[J].石油钻探技术,2017,45(4):115-120. YANG Zhen,XIAO Hongbing,LI Cui.Impacts of accuracy of azimuthal electromagnetic logging-while-drilling on resistivity and interface prediction[J].Petroleum Drilling Techniques,2017,45(4):115-120. [18] 倪卫宁,张晓彬,万勇,等.随钻方位电磁波电阻率测井仪分段组合线圈系设计[J].石油钻探技术,2017,45(2):115-120. NI Weining,ZHANG Xiaobin,WAN Yong,et al.The design of the coil system in LWD tools based on azimuthal electromagnetic-wave resistivity combined with sections[J].Petroleum Drilling Techniques,2017,45(2):115-120. -
期刊类型引用(4)
1. 田志宾，张嗣祺，杨庚佳，罗小兵. 超高温井下热管理系统优化设计与实验研究. 华中科技大学学报(自然科学版). 2025(04): 150-156 . 
百度学术
2. 田志宾，彭嘉乐，鄢星宇，魏赞庆，杨庚佳，罗小兵. 测井仪被动式热管理系统室温冷却研究. 石油钻探技术. 2024(01): 146-154 . 本站查看
3. 顾玉洋，夏竹君，罗鹏，王勇，李世举. MRCT-HT旋转井壁取心技术在南海东部A井中的应用. 石油和化工设备. 2024(04): 110-112+103 . 百度学术
4. 张波涛，罗鸣，张万栋，刘峰，孙艳军. 海上超高温高压井取心工具关键性能评价及应用研究. 中国科技论文. 2024(07): 760-768 . 百度学术
其他类型引用(0)
下载:
计量
- 文章访问数: 4759
- HTML全文浏览量: 75
- PDF下载量: 7347
- 被引次数: 4
