莺歌海盆地有着丰富的天然气资源^[1]，同时也存在异常高温高压地层^[2-6]。东方X区块位于莺歌海盆地的莺歌海凹陷南部，具有良好的成藏地质条件，但该区块构造复杂、岩性预测难、存在异常高温高压地层，并且安全密度窗口窄，造成钻井作业风险高，给钻井、地质工程带来了极大挑战。

针对异常高压地层的钻探风险，近年来主要利用井筒电缆VSP资料预测地层压力^[7-10]，或利用随钻声波测井资料进行地层压力预测和监测^[11-12]，或将电缆VSP资料和随钻声波测井资料相结合监测地层压力 ^[13]。但是，利用电缆VSP资料预测地层压力存在以下不足：1）进行一次电缆VSP测量需要3～5 d时间，作业时间长；2）利用电缆VSP资料预测目的层地层压力的精度与测量时井底距目的层的距离密切相关，需要尽量靠近目的层才能获得较高的预测精度。随钻地震技术将随钻地震资料和其他随钻测井数据实时结合起来，实时更新时深关系，标定层位，随钻监测钻头前方地层压力。与电缆VSP相比，其具有以下技术优势：1）时效性。随钻地震是在接钻柱间隙进行数据采集，不额外占用平台时间。2）实时性。随钻地震数据通过钻井液脉冲遥测系统实时传输到地面和解释中心。3）重复性。随钻地震在同一深度或关键层位可重复测量，从而保障数据和监测结果的准确性。4）实用性。随钻地震可以边钻边测，而不是靠近目的层才开始测量。

鉴于随钻地震技术具有上述优势，东方X区块的第一口探井DF-X1井在钻井过程中应用了随钻地震技术，不仅获得了高质量的时深关系，准确预测了关键标志层位T30、T31和储层黄流组二段A1砂体的深度，而且在钻入高压储层之前准确预测了高压储层的孔隙压力系数和破裂压力系数，并利用该技术优化了钻井参数，成功将ϕ244.5 mm套管下至高压储层上部的泥岩层，确保了高压储层ϕ212.7 mm井段有可控的安全密度窗口，保证了该井的顺利完钻。

1.   随钻地震数据的采集

随钻地震技术利用随钻地震仪器测量地震波从地表传播到井下接收器之间的时间，并记录4分量波形数据以提高解释精度。采用气枪震源在地表或者海面激发地震波，钻具组合中的随钻地震仪器在井眼中进行数据采集。该技术最大的特点是具有实时性，能够提供实时时深关系、速度和波场数据，并且是在接钻杆的间隙进行震源激发和数据采集，不额外占用平台时间。由于地面地震速度信息的不确定性，钻前的时深关系模型往往存在较大的不确定性。随钻地震seismicVISION测量仪器可减小这种不确定性，在直井利用该仪器可预测钻头前方的层位和地层压力，为钻井决策提供依据，降低钻井风险。

seismicVISION测量仪器内有4个传感器、1个井下处理器和1块内存。4个传感器为3个正交分量检波器和1个水听器。DF-X1井采用的seismicVISION测量仪器耐温150 ℃，耐压172.25 MPa。随钻地震数据采集如图1所示。其工作流程为：在连接钻杆的过程中，钻井作业和循环作业均停止，气枪震源激发，每一道数据震源激发10～15次，seismicVISION测量仪器实时记录地震数据，地震数据在井下完成自动叠加和初至拾取。时长约500 ms的地震数据通过钻井液脉冲遥测系统传输至地面，然后通过网络将地震数据传送到数据处理中心，处理人员对初至时间拾取进行快速质量控制，获得时深关系。此外对实时数据进行波场处理，获得走廊叠加道。使用上述时深关系将时间域的地面地震剖面转换到深度域，从而获得目的层位的深度信息，进而了解钻头当前位置与目的层之间的距离。实时作业过程中，每接收到一个深度的地震数据，都要完成上述处理过程。

图  1  随钻地震数据采集示意

Figure  1.  Schematic of seismic data acquisition while drilling

下载: 全尺寸图片 幻灯片

2.   随钻地震数据的处理与利用

2.1   实时时深关系和层位标定

随钻地震数据的处理方法与电缆VSP基本相同，处理流程如图2所示。随钻地震数据处理可以分为常规处理和特殊处理2部分。其中，常规处理主要包括数据叠加、初至提取、频谱分析、静校正、振幅恢复、波场分离、反褶积和走廊叠加。不同于电缆VSP，随钻地震的初至时间拾取可以在井下自动完成并传至地面，解释人员进行现场质控后得到初至时间，通过几何校正将其校正为垂直入射时间，得到时深关系。

图  2  随钻地震数据处理流程

Figure  2.  Process flow of seismic data while drilling

下载: 全尺寸图片 幻灯片

走廊叠加处理后，通过对比走廊叠加道和地面地震剖面，对主要标志层位进行标定，以确定走廊叠加道和地面地震之间的时间漂移量。然后，再对比地面地震提取的时深关系，更新钻头前方地层的预测深度，确定高压层位置。

2.2   波阻抗反演获取地层层速度

随钻地震数据的特殊处理，主要是进行贝叶斯波阻抗反演得到钻头前方地层的层速度^[14]。基于地震褶积模型，地震记录可以看作是震源子波和地层界面反射系数的褶积：

式中：t为时间，s；S(t)为地震记录，本文采用随钻地震数据处理后得到的上行波场走廊叠加结果；W(t)为地震子波，可以从随钻地震数据处理后的下行波场中提取；R(t)为地层界面的反射系数。

按照水平层状地层模型，地层界面的反射系数和波阻抗的关系可以表示为：

式中：i指第i个地层界面；Z_i+1为第i+1层的波阻抗，kg/（s·m²）；Z_i为第i层的波阻抗，kg/（s·m²）；R_i为第i个界面的反射系数。

根据式（1）得到反射系数序列，再根据式（2）得到波阻抗剖面。波阻抗为地震波在地层中的传播速度和地层密度的乘积，其表达式为：

式中：ρ_i为地层密度，kg/m³；v_i为地层的层速度，m/s。

反演得到波阻抗剖面后，再输入地层的密度就可以得到地层的层速度剖面，包括钻头前方未钻地层的层速度。未钻地层的密度采用邻井相同地层的实测密度或根据研究区域的经验直接输入密度常值。获得钻头前方地层的层速度后，就可以用其来预测钻头前方地层的压力。

2.3   随钻地震地层压力监测

地层压力监测主要是监测地层的孔隙压力和破裂压力。地层孔隙压力通常是指地层孔隙中流体的压力，若地层中的流体压力因某种地质条件的作用高于正常静水压力，则称之为地层超压或地层异常高压。地层破裂压力是指钻井过程中使井壁岩石产生破裂形成裂缝的钻井液液柱压力。

有关地层孔隙压力和破裂压力的计算，目前已有很多经验方法^[15-19]。东方X区块目的层及其上覆岩层主要为砂泥岩地层，可采用Eaton法^[20]计算该区块地层的孔隙压力，计算公式为：

式中：p_p为地层孔隙压力，Pa；σ_v为上覆岩层压力，Pa；p_n为正常静水压力，Pa；Δt为实测地层声波时差，μs/m；Δt_n为地层正常压实情况下的声波时差，μs/m。

上覆岩层压力通过对上覆海水和地层密度积分得到^[21]：

式中：ρ_w为海水密度，通常为（1.03～1.06）×10³ kg/m³；g为重力加速度，m²/s；h为海水层厚度，m；ρ_i为第i层的地层密度，kg/m³；ΔD_i为深度采样间隔，m。

地层破裂压力采用Eaton关系式^[22]计算，其计算公式为：

式中：p_f为地层破裂压力，Pa；µ为泊松比。

在缺少实测数据的情况下，可以利用泊松比随深度变化的关系式计算泊松比：

式中：${D_{{\rm{v,bml}}}}$为以泥线为基准的垂深，m；A₁，A₂和A₃为模型参数。

随钻压力监测工作流程主要包括3方面（见图3）：1）钻前压力预测。基于邻井的测井数据及钻井、地质和地震资料，建立钻前压力预测模型。2）随钻压力监测。利用随钻地震数据、随钻测量数据、钻井和录井等资料，实时监测当前地层和钻头前方地层的孔隙压力与破裂压力，结合井下的当量循环密度（ECD）和井况，预防钻井风险。3）钻后评估。根据实测压力数据标定地层压力预测结果，改进、完善地层压力预测模型。

图  3  地层压力监测工作流程

Figure  3.  Workflow of formation pressure monitoring

下载: 全尺寸图片 幻灯片

钻井过程中，根据随钻地震提供的实时时深关系标定层位，利用随钻地震数据，声波、伽马和电阻率测井曲线及录井数据等，实时监测地层的孔隙压力和破裂压力。对于钻头前方未钻地层，利用随钻地震获得地层的层速度计算孔隙压力和破裂压力，实时提供前方地层压力的变化趋势，优化钻井液的密度，确保钻井安全。

3.   随钻地震资料的应用

DF-X1预探井位于莺歌海盆地莺歌海凹陷的东方X构造，钻探目的是落实黄流组二段海底扇含气性和天然气资源量，主要目的层为黄流组二段A1砂体。图4所示为DF-X1井过井地面地震剖面，主要地震标志界面有T30、T31和T40。A1砂体为海底扇水道砂，海底扇整体处于“泥包砂”沉积环境中。东方X区块的目的层温度为185～190 ℃，预测地层孔隙压力系数在2.1以上，是典型的高温高压风险区。

图  4  DF-X1井过井地面地震剖面

Figure  4.  Surface seismic profile of Well DF-X1

下载: 全尺寸图片 幻灯片

3.1   随钻地震数据采集和处理

根据DF-X1井实际钻井情况，震源采用3个G-Gun空气枪组合，位于海平面下6.80 m处，单个气枪压力15.16～16.19 MPa。井下使用seismicVISION测量仪器，采集4分量（水听器+3分量检波器）地震数据，道间距为14.40 m。DF-X1井中进行了3次随钻地震数据采集（见表1），共计采集96道地震数据，测量井段共计1 004.16 m。随钻地震第一次数据采集是在ϕ339.7 mm套管内进行的，第二次和第三次数据采集是在ϕ311.1 mm裸眼井中进行的。

表  1  DF-X1井随钻地震各趟采集参数

Table  1.  Parameters acquisition of each run in Well DF-X1 by seismic while drilling

采集次数	采集道数	起始深度/m	结束深度/m	井筒类型	道间距/m
1	22	2 767.06	3 067.20	套管	14.40
2	30	3 081.17	3 381.12	裸眼	14.40
3	44	3 369.60	3 771.22	裸眼	14.40

下载: 导出CSV 

| 显示表格

钻井液脉冲遥测系统将地震波形数据和自动拾取的初至时间从井下传输至地面，然后通过平台加密网络将数据传输到处理中心。对接收到的数据进行检查，控制其质量，并审核井下自动拾取的初至时间获得时深关系。图5所示为3次随钻地震实时波形数据。由图5可知，初至起跳干脆（拾取的初至时间为红色实线）。

图  5  DF-X1井随钻地震实时数据

Figure  5.  Real-time seismic data while drilling in Well DF-X1

下载: 全尺寸图片 幻灯片

3.2   地层深度预测

实时作业过程中，每接收到一个深度的地震数据，都要进行上述处理过程。下面详细介绍第三次随钻地震数据的处理情况。图6所示为第三次实时时深关系和钻前时深关系的对比情况。从图6可以看出：在深度较浅时，相同时间下钻前预测模型预测的深度与随钻地震数据的深度有差别，随着深度增大二者逐渐靠近。

图  6  随钻地震实时时深关系与钻前及邻井时深关系对比

Figure  6.  Comparison of the real-time time-depth relationship among seismic while drilling, the pre-drilling and that in adjacent wells

下载: 全尺寸图片 幻灯片

对实时随钻地震数据进行波场处理，获得走廊叠加道。通过对比走廊叠加道和地面地震，获得地面地震双程时的漂移量。DF-X1井的走廊叠加道和地面地震存在17 ms的双程时时差。对从时间域地震数据上的层位双程时做相应漂移，漂移之后二者基本对应。根据时深关系转换到深度域，对目的层深度进行随钻预测和标定，得到钻头当前位置与目的层之间的距离。图7所示为第三次实时随钻地震数据处理结果，从左到右分别为双程时的反褶积后上行波场、走廊叠加道、增益走廊叠加道及走廊叠加道和地面地震的对比，主要层位在随钻地震走廊叠加道上都明显可见。图7中TD_now标记为第三次随钻地震采集时的井底位置。

图  7  第三次随钻地震数据实时处理结果

Figure  7.  Real-time processing results of the seismic data while drilling in the third run

下载: 全尺寸图片 幻灯片

实际钻井过程中不断调整处理储层上部的标志层参数，来提高随钻地震预测精度。黄流组二段A1砂体上部，主要地震标志层位有4个，分别为T30、T30_C、T30_E和T30_H。利用3次实时随钻地震数据对这4个标志层的深度进行预测，预测结果和最终实钻深度见表2。由于T30标志层上获得的地震道数据较少，预测误差为13.00 m。随着地震道数据增加，T30_C，T30_E和T30_H等3个标志层的预测误差分别为5.00，4.00和1.00 m。随着继续钻进，对目的层深度进行更新，钻头与目的层距离小于100 m时，发现钻遇砂岩的层速度增大，且砂岩气显示很好。利用第三次随钻地震数据预测的T31层位深度与实钻深度的误差为3.00 m，A1砂体顶预测深度与实钻深度的误差为6.00 m。

表  2  随钻地震地层深度实时预测结果

Table  2.  Real-time formation depth prediction results of seismic while drilling

地震层位	随钻地震预测深度/m	实钻深度/m	误差/m
T30	3 158.00	3 145.00	+13.00
T30_C	3 461.00	3 466.00	–5.00
T30_E	3 584.00	3 580.00	+4.00
T30_H	3 641.00	3 642.00	–1.00
T31	3 787.00	3 784.00	+3.00
A1砂顶	3 825.00	3 819.00	+6.00

下载: 导出CSV 

| 显示表格

3.3   随钻压力监测

受采集环境及仪器的限制，随钻地震数据信噪比低于电缆VSP，数据低频成分低至10 Hz左右，因此通常不做反演处理，只根据初至时间提取获得时深关系。DF-X1井随钻地震数据实时处理过程中，为降低钻遇异常高压层的风险，尝试对随钻地震数据进行反演处理，来获得钻头前方地层的层速度。

由于随钻地震资料的频谱是带限的，因此波阻抗反演过程中利用随钻声波数据和录井获得的地层密度数据进行约束，反演结果如图8所示。图8中从左到右分别为走廊叠加道、利用反演得到的阻抗计算得到的合成记录、走廊叠加道和合成记录的差、反射系数、波阻抗及反演得到的地层层速度。走廊叠加道和合成记录的差值越小，表明反演结果越好。从图8可以看出，在2.465 s处走廊叠加道和合成记录的差值较大，反演层速度明显增大，实测层速度没有增大，这主要是由走廊叠加道上同一时间的波峰造成，因此判定这个异常反演结果与随钻地震波场处理有关，而非存在高速薄层。

图  8  随钻地震波阻抗反演结果

Figure  8.  Results of seismic wave impedance inversion while drilling

下载: 全尺寸图片 幻灯片

利用时深关系，将时域地层层速度转换到深度域，同时与随钻自然伽马、声波和电阻率测井数据、当量循环密度及气测录井资料结合监测地层压力，预测钻头前方地层的孔隙压力。图9所示为利用随钻地震数据对ϕ311.1 mm井段和钻头前方地层孔隙压力和破裂压力的监测结果。从图9可以看出，钻至井深3 751.00 m时，利用随钻地震数据预测的地层层速度在井深3 770.00 m以深明显降低，表明地层孔隙压力升高。预测钻头前方地层孔隙压力系数为2.20～2.22，结合随钻地震层位预测结果，预测T31层在井深3 787.00 m处，A1砂体顶部在井深3 825.00 m处，于是决定在钻过T31层后，在A1砂体上部泥岩层中完，最后ϕ311.1 mm井段按照随钻监测结果钻至井深3 798.00 m中完。

图  9  随钻地震监测地层压力监测结果

Figure  9.  Formation pressure monitoring results of seismic while drilling

下载: 全尺寸图片 幻灯片

ϕ244.5 mm套管准确下到高压储层A1砂体上方的泥岩中，使ϕ212.7 mm井段在后续钻井过程中有足够的安全压力窗口，确保了储层井段顺利完钻。完钻后进行了电缆声波测井，并测量了地层孔隙压力，结果如图10所示。从图10可以看出，随钻地震时差与电缆声波时差整体匹配得很好。随钻地震预测A1砂体地层的孔隙压力系数为2.20～2.22，实测地层孔隙压力系数为2.27，预测精度达到3.0%。ϕ244.5 mm套管鞋处的地层破裂压力系数预测结果是2.32，实际地漏测试结果是2.34，预测精度达到1.0%。

图  10  随钻地震储层段地层压力监测结果与实测结果对比

Figure  10.  Comparison of the formation pressure monitoring results of seismic while drilling and that of the actually measured in the reservoir section

下载: 全尺寸图片 幻灯片

DF-X1井钻井过程中，利用随钻地震数据准确预测了地层界面深度，关键层位T31预测精度为3.00 m，A1砂体深度预测精度为6.00 m。A1砂体地层孔隙压力系数预测误差为0.05，预测精度达到3.0%；地层破裂压力系数预测误差为0.02，预测精度达到1.0%。随钻地震监测到地层孔隙压力系数在井深3 770.00 m以深上升，并结合层位预测结果在实际钻井过程中确定了ϕ311.1 mm井段的中完井深，成功将ϕ244.5 mm套管下至高压储层A1砂体上方的泥岩层中，确保了ϕ212.7 mm井段的安全压力窗口，保证了该井安全顺利完钻。

4.   结　论

1）应用随钻地震技术，成功预测了莺歌海盆地关键的地质层位和异常高压层，提高了钻井过程中地层层位的预测精度，增强了随钻地层压力监测的时效性和准确性。

2）随钻地震技术可以在不增加钻时的情况下，重复对目标地层进行数据采集，获得钻头前方地层的层速度，为随钻压力预监测提供了可靠的数据基础，提高了随钻压力监测精度，降低了海上高温高压地层的钻探工程风险。

3）现场应用表明，综合利用随钻地震数据、随钻测井数据、钻井参数、录井气测及电测数据等，可以监测地层压力和预测异常高压地层，对钻井风险进行提示，可以保障高温高压地层钻井安全顺利进行。