气体钻井声波超前测距方法与数值模拟

简旭; 李皋; 王军; 韩旭; 黄兵; 王松涛

doi:10.11911/syztjs.2022016

气体钻井声波超前测距方法与数值模拟

简旭^{1, 2,},
李皋²,
王军¹,
韩旭³,
黄兵¹,
王松涛¹

1.
中国石油集团川庆钻探工程有限公司川东钻探公司，重庆 401147
2.
油气藏地质及开发工程国家重点实验室（西南石油大学），四川成都 610500
3.
中国石油西南油气田分公司蜀南气矿，四川泸州 646000

基金项目: 国家自然科学基金重点项目“气体钻井安全监测的前兆预警关键传感器研究”（编号：61731016）资助

详细信息

作者简介:
简旭（1991—），男，河南南阳人，2014年毕业于长江大学勘查技术与工程专业，2020年获西南石油大学油气井工程专业博士学位，工程师，主要从事气体钻井、随钻声波探测方面的研究。E-mail：jxsky-blue@foxmail.com

中图分类号: P631.8⁺1
计量
- 文章访问数: 317
- HTML全文浏览量: 153
- PDF下载量: 48
出版历程
- 收稿日期: 2021-04-13
- 修回日期: 2022-01-14
- 录用日期: 2022-03-08
- 网络出版日期: 2022-04-07
- 刊出日期: 2022-06-08

Acoustic Advance Ranging Method in Gas Drilling and Its Numerical Simulation

JIAN Xu^{1, 2,},
LI Gao²,
WANG Jun¹,
HAN Xu³,
HUANG Bing¹,
WANG Songtao¹

1.
Chuandong Drilling Company, CNPC Chuanqing Drilling Engineering Company Limited, Chongqing, 401147, China
2.
State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu, Sichuan, 610500, China
3.
Southern Sichuan Gas District, PetroChina Southwest Oil & Gasfield Company, Luzhou, Sichuan, 646000, China

摘要

摘要:
气体钻井的安全性受储层地质条件影响明显，在钻井过程中，受气体所产生噪声影响，地层反射信号识别难度大。为提高探测准确度，降低气体钻井过程中钻头前方岩性界面的不确定性，设计了近钻头随钻冲击震源短节，基于地震波自激自收原理，提出了一种适用于气体钻井的声波超前测距方法。通过冲击试验分析聚四氟乙烯对振动尾波的衰减效果，通过数值模拟验证气体钻井声波超前测距方法的可行性，分析探测距离对反射波振幅的影响。冲击试验表明，使用聚四氟乙烯后空心圆柱上振动尾波衰减显著。数值模拟结果表明，PML边界条件可有效消除模型边界反射的干扰，利用反射纵波信号到达时间所计算测距的误差为1 m。随着探测距离增大，反射纵波相对强度逐渐增大，反射波振幅显著降低。研究结果验证了近钻头冲击震源探测钻头前方岩性界面的可行性，为气体钻井随钻声波超前探测技术的发展提供了新思路。
- 气体钻井 /
- 近钻头 /
- 冲击震源 /
- 超前探测 /
- 岩性界面
Abstract:
The safety of gas drilling is greatly influenced by local geological conditions of reservoirs. It is difficult to identify formation reflection signals in gas drilling due to the influence of noise generated from gas. Therefore, a near-bit impact source sub while drilling was designed, and a method of acoustic advance ranging for gas drilling was proposed on the basis of the self-excitation and self-reception principle of seismic waves. In this way, measurement accuracy could be raised, and the uncertainty of lithological interfaces ahead of the bit could be reduced. Specifically, the attenuation effect of polytetrafluoroethylene (PTFE) on vibration coda waves was analyzed through an impact test; the feasibility of the acoustic advance ranging method in gas drilling was verified through numerical simulations, and the influence of detection distance on the amplitude of reflected waves was analyzed. The impact test revealed that the coda waves of vibration waves on the hollow cylinder were significantly attenuated upon the use of PTFE. Numerical simulations indicated that the boundary conditions of a perfectly matched layer (PML) can effectively eliminate the interference of boundary reflection of the model, and the ranging error calculated from the arrival time of reflected P-waves was 1 m. In addition, the relative intensity of reflected P-waves shows a gradual increase, and the amplitude of reflected waves decreases greatly as the increase in detection distance. The research results verify the feasibility of detecting lithological interfaces ahead of the bit by near-bit impact sources and provide a new idea for the development of acoustic advance detection technology while drilling in gas drilling.
- gas drilling /
- near-bit /
- impact source /
- advance detection /
- lithological interface

HTML全文

气体钻井具有提高机械钻速和保护储层的优势，但储层地质条件对其安全性有较大影响^[1]。随钻超前探测技术可提前预测到钻头前方的岩性界面、断层、溶洞和异常高压等，降低钻头前方地层信息的不确定性^[2-3]，保证气体钻井的安全性。目前的随钻超前探测技术主要包括随钻地震前探技术、近钻头电磁波前视技术、钻头随钻地震技术^[4]。随钻地震前探技术利用近钻头处的检波器接收地表可控震源激发的直达波信号和钻头前方岩性界面产生的反射波信号，但该技术分辨率低，地震波的能量和频率被地层显著衰减，影响了钻头前方岩性界面的探测精度^[5-9]。为提高随钻前探技术的分辨率，斯伦贝谢和哈里伯顿公司研制了电磁波前视工具，实时获取钻头前方地层电阻率的变化，探测直井或低角度斜井中地层边界^[10-13]。但采用该技术进行远距离探测时，需要大功率储能装置，且气体介质对电磁波的衰减作用更强。钻头随钻地震技术利用牙轮钻头的轴向振动，实现钻头前方探测^[14]。但牙轮钻头信号强度弱，地面噪音对信号的干扰大。为了克服牙轮钻头震源的不足，哈里伯顿公司研制了扫频水力脉冲工具^[15]，该工具可在近钻头处激发宽频地震波信号，但连续冲击振动波信号会影响地层反射波信号的识别。

综合分析上述超前探测技术的不足，为降低地层对高频信号的衰减，提出了一种适用于气体钻井的随钻声波超前测距方法，并以此为基础，设计了无需井下供电的随钻冲击震源短节，以适应远距离探测。通过开展冲击试验，分析了聚四氟乙烯对空心圆柱振动尾波的衰减效果，并通过数值模拟验证了该测距方法的可行性。

1. 随钻声波超前测距方法

1.1 测距原理

在钻柱上进行震源激发和接收，受噪声等多因素影响，地层反射波信号识别难度大。因此，为避免噪声等的影响，选择钻柱停止钻进且保持循环通气的钻杆连接期间进行声源激发。采用地震波自激自收方式探测钻头前方地层信息，在近钻头处安装声源激发器和检波器，利用检波器接收钻头前方岩性界面的反射波，获取地层反射波的到达时间，在同一地层钻进一定距离后，再次激发声波获取地层反射波的到达时间。根据钻进距离以及反射波的时间差，计算钻头到前方岩性界面的距离。

当钻头前方存在岩性界面时，冲击震源短节激发地震波信号，如图1所示（E代表冲击震源，R代表检波器），反射波的到达时间为 $\Delta {t_1}$ ；同一地层继续钻进距离S₀后，冲击震源再次激发同样的地震波信号，反射波的到达时间为 $\Delta {t_2}$ 。

图 1 近钻头声波测距示意

Figure 1. Near-bit acoustic ranging

下载: 全尺寸图片幻灯片

钻头到前方岩性界面的距离S可表示为：

$S = \frac{{{v_{\text{p}}}\left( {\Delta {t_2} - \dfrac{{2{S_1}}}{{{v_{{\text{d}}}}}} - \dfrac{{{S_2}}}{{{v_{{\text{d}}}}}}} \right)}}{{2\cos \varphi }}$

(1)

${v_{\text{p}}} = \frac{{2{S_0}\cos \varphi }}{{\Delta {t_1} - \Delta {t_2}}}$

(2)

式中：S为钻头到岩性界面的距离，m；S₁为冲击震源短节到检波器的距离，m；S₂为检波器到钻头端部的距离，m；v_d为机械钻速，m/s；v_p为地层钻速，m/s； $\varphi$ 为岩性界面倾角，（°）；S₀为钻头钻进的距离，m； $\Delta {t_1}$ 和 $\Delta {t_2}$ 为反射波的到达时间，s。

将式（2）代入式（1）可得:

$S = \frac{{{S_0}\left( {\Delta {t_2} - \dfrac{{2{S_1}}}{{{v_{{\text{d}}}}}} - \dfrac{{{S_2}}}{{{v_{{\text{d}}}}}}} \right)}}{{\Delta {t_1} - \Delta {t_2}}}$

(3)

当待钻地层速度、岩性界面倾角未知时，可通过式（3）求得钻头至前方岩性界面的距离。

1.2 随钻冲击震源设计

根据上述思路设计了一种适用于气体钻井的随钻冲击震源短节，剖面如图2所示。上接头与短节本体采用花键接头，可为钻柱传递扭矩。冲击震源短节内使用聚四氟乙烯圆筒衰减钻柱振动波的尾波^[16]。高性能磁钢可吸住滑块，防止滑块二次冲击，降低二次冲击波对地层反射波的干扰。

图 2 冲击震源短节结构示意

1. 上接头；2. 花键接头；3. 台阶；4. 内套筒；5. 锥形螺纹a；6. 弹簧；7. 外套筒；8. 滑块；9. 高性能磁钢；10. 短节本体；11. 下接头；12. 聚四氟乙烯圆筒A；13. 传动螺纹；14. 变径接箍；15. 聚四氟乙烯圆筒B；16. 锥形螺纹b

Figure 2. Structure of the impact source sub

下载: 全尺寸图片幻灯片

冲击震源短节激发振动波时，上提上接头，逆时针小扭矩旋转上接头使内套筒传动螺纹与滑块传动螺纹脱离，弹簧推动滑块冲击高性能磁钢，并吸附在高性能磁钢上。震源短节完成振动波激发后，下放上接头使内套筒传动螺纹部分嵌入滑块传动螺纹，上提上接头一定距离，顺时针旋转使内套筒传动螺纹完全嵌入滑块传动螺纹，最后下放上筒体使上筒体与中筒体的花键接头配合。

2. 冲击试验

为模拟钻柱上冲击振动波的特征，将速度传感器固定在空心圆柱表面，开展冲击试验。将圆柱形304不锈钢放置于空心圆柱顶端，使用锤子敲击；随后在空心圆柱顶端依次放置圆柱形聚四氟乙烯、304不锈钢，使用相同锤子敲击，分析得到振动波形和频谱（见图3），使用聚四氟乙烯后所接收首波为雷克子波，首波后尾波被显著衰减，振动波主频为425.7 Hz，该振动波形有利于地层反射波信号识别。

图 3 使用聚四氟乙烯所得振动波形和频谱

Figure 3. Vibration waveforms and spectra received using PTFE

下载: 全尺寸图片幻灯片

3. 近钻头反射波场数值模拟

3.1 三维黏弹性波动方程有限差分

为分析地层对声场的耗散效应，选择Kelvin-Voight模型模拟近钻头井眼声场的传播特征^[17-18]。交错网格中波场分量在时间和空间中的位置如图4所示（其中，空间步长 $\Delta x = \Delta y = \Delta z$ ，时间步长为 $\Delta t$ ）。

图 4 交错网格中的波场分量

Figure 4. Wavefield components in staggered grids

下载: 全尺寸图片幻灯片

不同波场分量在时间和空间上的差分格式如下：

$\begin{array}{c} \sigma _{mm}^n(i,j,k) = {\sigma _{mm}}\left[ {\left(i + \dfrac{1}{2}\right)\Delta x,\left(j + \dfrac{1}{2}\right)\Delta y,\left(k + \dfrac{1}{2}\right)\Delta z;n\Delta t} \right] \qquad (m=x,y,z) \\ \end{array}$

(4)

$\tau _{xy}^n(i,j,k) = {\tau _{xy}}\left[ {i\Delta x,j\Delta y,\left(k + \frac{1}{2}\right)\Delta z;n\Delta t} \right]$

(5)

$\tau _{xz}^n\left( {i,j,k} \right) = {\tau _{xz}}\left[ {i\Delta x,\left( {j + \frac{1}{2}} \right)\Delta y,k\Delta z;n\Delta t} \right]$

(6)

$\tau _{yz}^n(i,j,k) = {\tau _{yz}}\left[ {\left( {i + \frac{1}{2}} \right)\Delta x,j\Delta y,k\Delta z;n\Delta t} \right]$

(7)

$v_x^{n + {1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}\left( {i,j,k} \right) = {v_x}\left[ {i\Delta x,\left( {j + \frac{1}{2}} \right)\Delta y,\left( {k + \frac{1}{2}} \right)\Delta z;\left( {n + \frac{1}{2}} \right)\Delta t} \right]$

(8)

$v_y^{n + {1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}\left( {i,j,k} \right) = {v_y}\left[ {\left( {i + \frac{1}{2}} \right)\Delta x,j\Delta y,\left( {k + \frac{1}{2}} \right)\Delta z;\left( {n + \frac{1}{2}} \right)\Delta t} \right]$

(9)

$v_z^{n + {1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}\left( {i,j,k} \right) = {v_z}\left[ {\left( {i + \frac{1}{2}} \right)\Delta x,\left( {j + \frac{1}{2}} \right)\Delta y,k\Delta z;\left( {n + \frac{1}{2}} \right)\Delta t} \right]$

(10)

式中： ${v_x}$ ， ${v_y}$ 和 ${v_z}$ 分别为x，y和z方向的质点速度分量，m/s； ${\sigma _{xx}}$ ， ${\sigma _{yy}}$ 和 ${\sigma _{zz}}$ 分别为x，y和z方向的正应力分量，MPa； ${\tau _{xy}}$ ， ${\tau _{yz}}$ 和 ${\tau _{zx}}$ 分别为x，y和z方向的切应力分量，MPa。

式（8）在时间2阶空间4阶的差分格式为：

$\begin{array}{c} v_x^{n + 1/2}\left( {i,j + \dfrac{1}{2},k + \dfrac{1}{2}} \right) = v_x^{n - 1/2}\left( {i,j + \dfrac{1}{2},k + \dfrac{1}{2}} \right) {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} + \dfrac{{\Delta t}}{{{\rho _{i,j{\text{ + }}\dfrac{{\text{1}}}{{\text{2}}},k{\text{ + }}\dfrac{{\text{1}}}{{\text{2}}}}}\Delta x}}\displaystyle\sum\limits_{{\text{N}} = 1}^2 {{C_{\text{N}}}} \left[ \sigma _{xx}^n\left( {i + \dfrac{{2n - 1}}{2},j + \dfrac{1}{2},k + \dfrac{1}{2}} \right) - \right. \\ \left. \sigma _{xx}^n\left( {i - \dfrac{{2n - 1}}{2},j + \dfrac{1}{2},k + \dfrac{1}{2}} \right) \right] {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} + \dfrac{{\Delta t}}{{{\rho _{i,j{\text{ + }}\dfrac{{\text{1}}}{{\text{2}}},k{\text{ + }}\dfrac{{\text{1}}}{{\text{2}}}}}\Delta y}}\displaystyle\sum\limits_{{\text{N}} = 1}^2 {{C_{\text{N}}}} \left[ \tau _{xy}^n\left( {i,j + n,k + \dfrac{1}{2}} \right) -\right. \left. \tau _{xy}^n\left( {i,j + 1 - n,k + \dfrac{1}{2}} \right) \right] \hfill \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} + \dfrac{{\Delta t}}{{{\rho _{i,j{\text{ + }}\dfrac{{\text{1}}}{{\text{2}}}{\text{,}}k{\text{ + }}\dfrac{{\text{1}}}{{\text{2}}}}}\Delta z}}\displaystyle\sum\limits_{{\text{N}} = 1}^2 {{C_{\text{N}}}} \left[ \tau _{xz}^n\left( {i,j + \dfrac{1}{2},k + n} \right) - \right. \left. \tau _{xz}^n\left( {i,j + \dfrac{1}{2},k + 1 - n} \right) \right] \hfill \\ \end{array}$

(11)

3.2 稳定性条件和边界吸收条件

进行有限差分数值模拟时，为避免产生较大数值频散，时间步长、空间步长以及物理参数的选取均需满足差分计算的稳定性要求，三维有限差分格式的稳定性条件为^[17]：

$\Delta t \leqslant \dfrac {1}{{v_{\text{p}}}\sqrt {\Delta {x^{ - 2}} + \Delta {y^{ - 2}} + \Delta {z^{ - 2}}}}$

(12)

式中， $\Delta x$ ， $\Delta y$ 和 $\Delta z$ 分别为x，y和z方向的网格步长。

为消除人工边界反射对数值模拟结果的影响，采用完全匹配层（PML）吸收边界条件（见图5），当声波从计算区域传播到PML层时，声波能量以指数形式衰减。

图 5 三维PML吸收边界示意

Figure 5. Three-dimensional (3D) PML absorption boundary

下载: 全尺寸图片幻灯片

为分析三维PML边界条件的吸收效果，利用三维数值模拟程序模拟声波在砂岩地层的传播特征。基于冲击试验结果，选取雷克子波为声源。数值模拟模型网格大小300×300×300，空间步长为0.01 m，时间步长为1×10⁻⁵ s，声源主频为500 Hz，边界层厚度为1.50 m。砂岩纵波波速为3 563 m/s，横波波速为2 354 m/s，密度为2.65 g/cm³，纵波品质因子为90，横波品质因子为60。使用三维PML边界条件后，声波在砂岩中的三维波场快照如图6所示。当声波传播到模型边界时，未出现界面反射现象，表明其消除模型界面反射的效果较好。

图 6 使用PML边界条件后的三维波场快照

Figure 6. 3D wavefield snapshots after using PML boundary conditions

下载: 全尺寸图片幻灯片

3.3 可行性验证

为了验证井眼有限差分模拟的准确性，计算了轴对称情况下钻柱内接收到的波形，并与实轴积分法的结果进行对比。采用声源主频为500 Hz的雷克子波，空间步长为0.01 m，时间步长为1×10⁻⁵ s。声波模拟的介质参数见表1，钻柱纵向振动时纵波的品质因子为4～15^[19]，取钻铤纵波品质因子为10，通过纵横波品质因子之间的关系求取横波品质因子^[20]。对比有限差分法（FDTD）和实轴积分法（RAI）计算出的波形（见图7），发现有限差分计算出的波形存在频散，导致尾波比较显著，但两者的钻铤模式波和钻铤界面反射波吻合较好，验证了有限差分数值模拟算法的可行性。

表 1 计算模型中不同介质的参数

Table 1. Parameters of different media in the calculation model

介质类型	纵波速度/（m·s⁻¹）	横波速度/（m·s⁻¹）	密度/（kg·m⁻³）	纵波品质因子	横波品质因子	外半径/m
钻铤	5 860	3 130	7 850	10	4	0.09
空气	380	0	20			0.12
砂岩	3 563	2 354	2 650	90	60	1.00
泥岩	2 600	1 700	2 200	60	35	1.00

下载: 导出CSV

| 显示表格

图 7 有限差分和实轴积分波形对比

Figure 7. Comparison of finite-difference and real-axis integral waveforms

下载: 全尺寸图片幻灯片

3.4 不同时刻声波传播特征

目前气体钻井所钻遇地层主要为砂泥岩地层，因此以砂岩和泥岩为介质进行数值模拟，介质参数见表1。三维井眼模型如图8所示，钻铤与井周地层之间的环空为气体。声源位于钻铤内，在钻铤底部依次布置接收器，钻铤底端与反射界面之间的距离为5 m。

图 8 三维井眼模型

Figure 8. 3D borehole model

下载: 全尺寸图片幻灯片

为了分析声波在三维井眼模型中的传播特征，记录了X-Z平面内不同时刻的波场快照（见图9）。三维井眼模型的空间步长为0.02 m，时间步长为1×10⁻⁵ s，声源主频为500 Hz，钻铤与前方岩性界面间的距离为10 m。径向上钻铤波经环空气体透射到地层，转换为地层纵波和横波；轴向上钻铤波直接透射到地层，形成地层纵波和横波，轴向纵波速度大于径向纵波速度。当地层纵波遇到前方岩性界面时，形成反射纵波和横波，声波振幅被明显衰减。岩性界面产生的反射纵波最先透射到钻铤内，被钻铤内的检波器接收到，因此可利用纵波探测前方岩性界面。随着传播距离增加，纵波振幅逐渐降低，导致钻铤内接收到的反射纵波振幅减弱。

图 9 不同时刻的波场快照

Figure 9. Wavefield snapshots at different moments

下载: 全尺寸图片幻灯片

3.5 探测距离对反射波振幅影响

利用三维声波数值模拟模型分析了探测距离对反射波振幅的影响，将钻铤底界面与岩性界面之间的距离从5 m增加到30 m，距离间隔为5 m。当探测距离为5和15 m时，与直达波振幅相比，波形图中反射波振幅非常弱，需要分离出直达波信号才能识别出反射波（见图10）。随探测距离增加，反射波到达时间增加，振幅降低。因此在远探距离探测中，反射波振幅要强于背景噪音，才有利于提取出弱地层反射波信号。当探测距离为5和15 m时，钻铤底界面与岩性界面之间的垂直距离分别为6.07 m和16.08 m，与实际探测距离之间的误差为1 m。当待钻岩石波速未知时，利用式（3）求得钻铤底界面与岩性界面之间的垂直距离为6.07 m，与已知探测距离条件下一致，验证了该测距公式的准确性。

图 10 钻铤内距钻铤底界面2 m处所接收波形

Figure 10. Waveforms received at places 2 m away from bottom interface within a drill collar

下载: 全尺寸图片幻灯片

当探测距离为5和15 m时，使用PML边界条件可以很好地消除模型边界反射的干扰（见图11），岩性界面处的反射波由地层传入钻铤，随探测距离增加，反射波的振幅显著降低，与时域波形中反射波振幅的变化趋势一致。

图 11 不同探测距离的波场快照

Figure 11. Wavefield snapshots at different detection distances

下载: 全尺寸图片幻灯片

拟合激发源处反射波相对强度与探测距离的关系曲线可得

${A_{\text{R}}} = 20\lg \frac{{{{A}_{\text{d}}}}}{{{{A}_{\text{r}}}}}$

(13)

式中：A_R为直达波与反射波之间的相对强度，dB；A_d为直达波的峰值幅度，g；A_r为反射波的峰值幅度，g。

根据钻铤内距钻铤底界面2 m处接收到的波形，分别读取不同探测距离下的直达波振幅和地层反射波振幅，利用式（13）求取直达波振幅与反射波振幅之间的相对强度（见图12），随着探测距离增加，反射波的相对强度逐渐增大。当探测距离大于20 m时，反射波相对强度的增加趋势变缓，随着探测距离增大，反射波振幅逐渐减小，且呈指数衰减。

图 12 探测距离与相对强度之间的关系

Figure 12. Relationship between detection distance and relative intensity

下载: 全尺寸图片幻灯片

4. 结论与建议

1）钻头前方存在异常地质体时会产生反射波信号，需降低井下噪音和钻柱振动尾波对地层反射波信号的干扰。将聚四氟乙烯应用在近钻头冲击震源短节上，可显著衰减空心圆柱的振动尾波，提高识别弱地层反射波信号的能力。

2）Kelvin-Voight模型和有限差分方法可用于模拟低频声源下三维井眼中声波的传播特征，获取钻头到前方岩性界面的距离；随着探测距离增大，反射纵波与直达波之间的相对强度逐渐增大，识别反射波信号的难度增大。

3）与气体钻井条件相比，常规钻井液钻井可在停钻接立柱期间停止钻井液循环，其井下噪音干扰少，有利于地层反射波信号的识别，因此后续可尝试将该声波超前测距方法用于常规钻井液钻井，研制相应样机在地面开展超前探测试验，进一步验证模型，改进工具结构。

图 1 近钻头声波测距示意

Figure 1. Near-bit acoustic ranging

下载: 全尺寸图片幻灯片

图 2 冲击震源短节结构示意

Figure 2. Structure of the impact source sub

下载: 全尺寸图片幻灯片

图 3 使用聚四氟乙烯所得振动波形和频谱

Figure 3. Vibration waveforms and spectra received using PTFE

下载: 全尺寸图片幻灯片

图 4 交错网格中的波场分量

Figure 4. Wavefield components in staggered grids

下载: 全尺寸图片幻灯片

图 5 三维PML吸收边界示意

Figure 5. Three-dimensional (3D) PML absorption boundary

下载: 全尺寸图片幻灯片

图 6 使用PML边界条件后的三维波场快照

Figure 6. 3D wavefield snapshots after using PML boundary conditions

下载: 全尺寸图片幻灯片

图 7 有限差分和实轴积分波形对比

Figure 7. Comparison of finite-difference and real-axis integral waveforms

下载: 全尺寸图片幻灯片

图 8 三维井眼模型

Figure 8. 3D borehole model

下载: 全尺寸图片幻灯片

图 9 不同时刻的波场快照

Figure 9. Wavefield snapshots at different moments

下载: 全尺寸图片幻灯片

图 10 钻铤内距钻铤底界面2 m处所接收波形

Figure 10. Waveforms received at places 2 m away from bottom interface within a drill collar

下载: 全尺寸图片幻灯片

图 11 不同探测距离的波场快照

Figure 11. Wavefield snapshots at different detection distances

下载: 全尺寸图片幻灯片

图 12 探测距离与相对强度之间的关系

Figure 12. Relationship between detection distance and relative intensity

下载: 全尺寸图片幻灯片

表 1 计算模型中不同介质的参数

Table 1 Parameters of different media in the calculation model

介质类型	纵波速度/（m·s⁻¹）	横波速度/（m·s⁻¹）	密度/（kg·m⁻³）	纵波品质因子	横波品质因子	外半径/m
钻铤	5 860	3 130	7 850	10	4	0.09
空气	380	0	20			0.12
砂岩	3 563	2 354	2 650	90	60	1.00
泥岩	2 600	1 700	2 200	60	35	1.00

下载: 导出CSV

参考文献(20)

[1]	李皋,李诚,孟英峰,等. 气体钻井随钻安全风险识别与监控[J]. 天然气工业,2015,35(7):66–72. doi: 10.3787/j.issn.1000-0976.2015.07.010 LI Gao, LI Cheng, MENG Yingfeng, et al. While-drilling safety risk identification and monitoring in air drilling[J]. Natural Gas Industry, 2015, 35(7): 66–72. doi: 10.3787/j.issn.1000-0976.2015.07.010
[2]	ESMERSOY C,HAWTHORN A,李辉,等. 降低钻头前的不确定性[J]. 国外测井技术,2016,37(3):61–68. ESMERSOY C, HAWTHORN A, LI Hui, et al. To reduce the uncertainty of the drill bit[J]. World Well Logging Technology, 2016, 37(3): 61–68.
[3]	杨书博,乔文孝,赵琪琪,等. 随钻前视声波测井钻头前方声场特征研究[J]. 石油钻探技术,2021,49(2):113–120. doi: 10.11911/syztjs.2021020 YANG Shubo, QIAO Wenxiao, ZHAO Qiqi, et al. The characteristics of the acoustic field ahead of the bit in“look-ahead”acoustic logging while drilling[J]. Petroleum Drilling Techniques, 2021, 49(2): 113–120. doi: 10.11911/syztjs.2021020
[4]	苏义脑,徐义,盛利民,等. 随钻地震技术研究进展[J]. 石油钻采工艺,2010,32(5):1–7. doi: 10.3969/j.issn.1000-7393.2010.05.001 SU Yinao, XU Yi, SHENG Limin, et al. Review on study progress of seismic while drilling technology[J]. Oil Drilling & Production Technology, 2010, 32(5): 1–7. doi: 10.3969/j.issn.1000-7393.2010.05.001
[5]	ESMERSOY C, UNDERHILL W, HAWTHORN A. Seismic measurement while drilling: Conventional borehole seismic on LWD[R]. SPWLA-2001-RR, 2001.
[6]	EL TOKHY M, FAHMY A, EL MARAKEBY H, et al. Seismic while-drilling, real-data reposition the well in the seismic volume and improve the accuracy of the depth prognosis of the target events[R]. SPE 176763, 2015.
[7]	史鸿祥,李辉,郑多明,等. 基于随钻地震测井的地震导向钻井技术:以塔里木油田哈拉哈塘区块缝洞型储集体为例[J]. 石油勘探与开发,2016,43(4):662–668. doi: 10.1016/S1876-3804(16)30085-4 SHI Hongxiang, LI Hui, ZHENG Duoming, et al. Seismic guided drilling technique based on seismic while drilling (SWD): a case study of fracture-cave reservoirs of Halahatang block, Tarim Oilfield, NW China[J]. Petroleum Exploration and Development, 2016, 43(4): 662–668. doi: 10.1016/S1876-3804(16)30085-4
[8]	药晓江,卢华涛,尚捷,等. 随钻测井仪流道转换器优化设计与数值分析[J]. 石油钻探技术,2021,49(5):121–126. doi: 10.11911/syztjs.2021069 YAO Xiaojiang, LU Huatao, SHANG Jie, et al. Optimization design and numerical analysis of flow passage converters in LWD tools[J]. Petroleum Drilling Techniques, 2021, 49(5): 121–126. doi: 10.11911/syztjs.2021069
[9]	高永德,刘鹏,杜超,等. 随钻地震技术在莺歌海盆地高温高压地层钻井中的应用[J]. 石油钻探技术,2020,48(4):63–71. doi: 10.11911/syztjs.2020049 GAO Yongde, LIU Peng, DU Chao, et al. The application of seismic while drilling in high temperature, high pressure reservoirs of the Yinggehai Basin[J]. Petroleum Drilling Techniques, 2020, 48(4): 63–71. doi: 10.11911/syztjs.2020049
[10]	CONSTABLE M V, ANTONSEN F, STALHEIM S O, et al. Looking ahead of the bit while drilling: from vision to reality[J]. Petrophysics, 2016, 57(5): 426–446.
[11]	HARTMANN A, GOREK M, FULDA C, et al. Early bed boundary detection while drilling-testing and application of a new resistivity device[R]. IPTC-12063-MS, 2008.
[12]	WU H H, GOLLA C, PARKER T, et al. A new ultra-deep azimuthal electromagnetic LWD sensor for reservoir insight[R]. SPWLA-2018-X, 2018.
[13]	黄明泉,杨震. 随钻超深电磁波仪器探测深度及响应特征模拟[J]. 石油钻探技术,2020,48(1):114–119. doi: 10.11911/syztjs.2019132 HUANG Mingquan, YANG Zhen. Simulation to determine depth of detection and response characteristics while drilling of an ultra-deep electromagnetic wave instrument[J]. Petroleum Drilling Techniques, 2020, 48(1): 114–119. doi: 10.11911/syztjs.2019132
[14]	ANCHLIYA A. A review of seismic while drilling (SWD) techniques: A journey from 1986 to 2005[R]. SPE 100352, 2006.
[15]	KOLLE J J, THEIMER K. Seismic-while-drilling using a swept impulse source[R]. SPE 92114, 2005.
[16]	JIAN Xu, LI Hongtao, LI Gao, et al. Lithological interface detection using an impact source[J]. Shock and Vibration, 2020, 2020: 4189419.
[17]	LIU Qinghuo, SCHOEN E, DAUBE F, et al. A three‐dimensional finite difference simulation of sonic logging[J]. The Journal of the Acoustical Society of America, 1996, 100(1): 72–79. doi: 10.1121/1.415869
[18]	CARCIONE J M, POLETTO F, GEI D. 3-D wave simulation in anelastic media using the Kelvin-Voigt constitutive equation[J]. Journal of Computational Physics, 2004, 196(1): 282–297. doi: 10.1016/j.jcp.2003.10.024
[19]	ANGONA F A. Drill string vibration attenuation and its effect on a surface oscillator drilling system[J]. Journal of Engineering for Industry, 1965, 87(2): 110–114. doi: 10.1115/1.3670771
[20]	WANG Yun, LU Jun, SHI Ying, et al. PS-wave Q estimation based on the P-wave Q values[J]. Journal of Geophysics and Engineering, 2009, 6(4): 386–389. doi: 10.1088/1742-2132/6/4/006

施引文献(10)

期刊类型引用(6)

1.	于瑞丰，刁斌斌，高德利，刘喆，张富磊. 救援井在穿越井段提高测量精度的计算方法研究. 钻采工艺. 2024(04): 132-139 . 百度学术
2.	张来斌，谢仁军，殷启帅. 深水油气开采风险评估及安全控制技术进展与发展建议. 石油钻探技术. 2023(04): 55-65 . 本站查看
3.	蒲文学，范光第，朱建建，赵国山. 磁性随钻测斜仪所需无磁钻具长度及影响因素研究. 石油钻探技术. 2022(04): 129-134 . 本站查看
4.	王瑞芹，李翠，李浩然，刘阳，范兰兰. 目标井口通电救援井连通探测系统测距导向算法研究. 内蒙古石油化工. 2022(08): 10-14 . 百度学术
5.	吴艳，余文艳，刘永辉，晋新伟，周岩，朱宽亮. 主动磁测距技术在裸眼落鱼井中的重入试验. 钻采工艺. 2022(06): 14-19 . 百度学术
6.	于瑞丰，刁斌斌，高德利. 基于邻井距离测量误差的救援井磁测距工具优选方法. 石油钻探技术. 2021(06): 118-124 . 本站查看

其他类型引用(4)

资源附件(0)

图(12) / 表(1)

计量

文章访问数: 317
HTML全文浏览量: 153
PDF下载量: 48
被引次数: 10

1. 随钻声波超前测距方法
1.1 测距原理
1.2 随钻冲击震源设计
2. 冲击试验
3. 近钻头反射波场数值模拟
3.1 三维黏弹性波动方程有限差分
3.2 稳定性条件和边界吸收条件
3.3 可行性验证
3.4 不同时刻声波传播特征
3.5 探测距离对反射波振幅影响
4. 结论与建议

1. 随钻声波超前测距方法
1.1 测距原理
1.2 随钻冲击震源设计
2. 冲击试验
3. 近钻头反射波场数值模拟
3.1 三维黏弹性波动方程有限差分
3.2 稳定性条件和边界吸收条件
3.3 可行性验证
3.4 不同时刻声波传播特征
3.5 探测距离对反射波振幅影响
4. 结论与建议

参考文献(20)

施引文献(10)

资源附件(0)

气体钻井声波超前测距方法与数值模拟

作者简介: 简旭（1991—），男，河南南阳人，2014年毕业于长江大学勘查技术与工程专业，2020年获西南石油大学油气井工程专业博士学位，工程师，主要从事气体钻井、随钻声波探测方面的研究。E-mail：jxsky-blue@foxmail.com

计量

出版历程