东海OBN和三维DAS-VSP数据的联合采集与处理方法研究

张少华,苟 量,余 刚,,刘海波,张 昊,曹中林,陈沅忠,,何光明,吴俊军,,王熙明,王艳华,

(1.中国石油集团东方地球物理勘探有限责任公司,河北涿州072751;2.中油奥博(成都)科技有限公司,四川成都611731;3.上海石油天然气有限公司,上海200041)

近年来,光纤传感技术已经应用于地面地震数据、海洋地震数据、井中地震数据和井-地联合地震数据的采集,推动了光纤传感技术在地球物理特别是地震数据采集中的应用。光纤传感技术是一项革命性的新技术,光纤因其体积小、不带电、分布式、高密度、多参量、耐高温、高压、全井段接收和低成本等特征,必将带来井下、海洋和陆地地球物理技术的一场革命。井中分布式光纤声波传感(distributed acoustic sensing,DAS)技术已广泛应用于井中VSP数据采集、水力压裂微地震监测和精准工程监测,可实现油气井全生命周期监测和管理。分布式光纤传感技术在油气资源勘探开发领域的规模化推广应用,已经从井中延伸到了陆地和海洋;从井下单分量测量拓展到了井下和陆地三分量测量(三分量分布式声波传感铠装光缆);从单井单参数测量发展到了多井多参数同步测量,调制解调仪器也从单通道单参数发展到了多通道多参数复合调制解调系统。分布式光纤声波传感器采用独特的分布式光纤探测技术,对沿光纤传输路径上的空间分布和随时间变化的信息进行测量或监控。该技术利用光纤本身作为传感器进行信号采集,在井中地震勘探、井中-地面立体勘探、油气藏长期动态监测和微地震监测等方面被广泛应用,凭借质量轻、高灵敏度、高空间分辨率、全井段、高效率、低成本等优势,近年来得到快速发展,成为贯穿油气井全生命周期的新兴油气藏勘探、开发和监测技术。

在井中地震勘探领域,井中DAS-VSP数据采集越来越被认为是用于采集井下垂直地震剖面数据的三分量地震检波器阵列的可行替代方案[1]。2012年底,某石油公司在墨西哥湾的深水工区内进行OBN数据采集时,在套管外预先安装了铠装光缆的两口井中同步联合采集了三维DAS-VSP数据,然后进行了三维DAS-VSP上行反射波数据的逆时偏移处理和多次波的逆时偏移处理,对盐丘下的油气储层构造进行了有效成像[2-3]。2015年,某公司在墨西哥湾的深水联合区块内进行时移三维OBN数据采集时,在油气管柱外预先安装了铠装光缆的5口井中同步联合采集了时移三维DAS-VSP数据,时移三维DAS-VSP数据的处理结果被用于监测该深水油田的注水驱油效果[4]。海上OBN和三维DAS-VSP数据联采的相关文献[5]为我们提供了有关海上三维DAS-VSP数据应用非常有价值的信息。WU等[6]在DAS-VSP数据采集井周围获得的高质量三维DAS-VSP数据成像结果也给了我们信心和决策支持,使我们能够同时进行新规划的东海高密度OBN数据和海洋三维DAS-VSP数据的井-海联合采集。

井-地或井-海联合地震勘探技术是陆地或海洋地震和井中垂直地震剖面技术相结合形成的一项立体地震勘探方法,实现了地面三维地震数据与井中地震数据的同步采集及同步处理,达到地面地震与井中地震优势互补的目的。VSP采集是将三分量检波器或铠装光缆置于井中,记录地表震源激发的地震波的直达和反射信号。相对于地面地震勘探,VSP的优点是具有准确的时深关系、井旁地层信息精确,并且波场信息丰富且直观,VSP观测到的反射波传播路径短、受近地表低降速带和环境噪声影响均较小,因此,VSP资料具有高分辨率和高信噪比等优势。同时,井下检波器更靠近储层,利用零井源距VSP和井-地或井-海联采的DAS-VSP资料可以得到准确的时深关系、地层速度、反褶积算子、球面扩散补偿因子TAR值、吸收衰减因子Q值和井眼周围的各向异性参数。这些参数可用于提升地表三维地震数据的处理效果,提高地震资料的保真度、分辨率和成像精度。井驱地面地震资料处理包括速度模型标定与修正、静校正、反褶积、去除多次波、高频恢复、各向异性偏移、Q补偿或Q偏移等。由于常规井下三分量检波器阵列的级数有限(不超过100级),无法在一次地面震源的激发中采集全井段的三维VSP数据,因此,难以利用井下三分量检波器阵列实现高效率低成本的全井段井-地或井-海三维地震数据的联采作业。利用套管内外或油气管柱内外布设的全井段铠装光缆,可以方便快捷地实现全井段VSP数据采集和井-地或井-海三维地震数据联合采集作业。井-地联合勘探提升了地震数据成像结果的精度和质量,提高了目的层反射波信噪比和分辨率,有利于识别地下特殊地质体,精细研究井周围构造、储层及油层的变化特征,从而提高地震数据对地质目标的认识、描述和表征能力[7]。

本文介绍了东海某OBN数据采集工区利用井下铠装光缆同步采集三维DAS-VSP数据,以及三维DAS-VSP数据的成像处理及初步解释结果。利用井-海联采的井中三维DAS-VSP数据,可以获得精确的时深关系、地层速度、反褶积算子、TAR值、Q值和井眼周围的各向异性参数。这些参数可用于提高地表三维地震数据的处理。井驱地面地震资料处理包括速度模型标定与修正、静校正、反褶积、多次波处理、高频恢复、各向异性偏移、Q补偿或Q偏移等。通过井-地或井-海联合地震勘探,可以大幅度提高地震数据的成像质量,提高目的层反射波的信噪比和分辨率,最终提高数据表征地质目标的能力。

1 工区地质背景及问题分析

1.1 地质背景

平湖油气田位于东海陆架盆地西湖凹陷西侧平湖断裂构造带中部,西侧以平湖主断层为界。平湖主断层呈北北东走向,延伸长度大于100km,对平湖油气田和平湖组沉积的形成起着重要控制作用。平湖油气田由放鹤亭、八角亭、望湖亭、双照亭等8个大小构造组成,其中放鹤亭、八角亭、望湖亭构造为主要构造,均有钻井控制,并试获工业油气流,即筹划开发的区块。

平湖油藏的特征主要为背斜型底水块状油藏。气藏主要为断块背斜构造和岩性控制的层状凝析气藏,受构造和岩性控制的高压层状凝析气藏。

西湖凹陷南北长约550km,东西宽约110km,面积约6×104km2,主要沉积层为第三系地层,可预测沉积厚度大于12000m。西湖凹陷构造区总体分为西部斜坡区、中央洼陷区、东部断阶区3大区带,目前主要油气田和井位分布在西部斜坡区和中央洼陷区的中央反转构造带、西次凹。

1.2 存在的问题和地质目标

平湖油气田在2007年采集了三维OBC数据,但是此OBC数据的质量和成果资料的品质制约着油藏评价和油田开发工作的进一步深入。评价、开发面临的主要问题有:①目的层埋藏较深(3300m以上),三维OBC数据的中深层成像品质较差,难以满足构造精细解释、优势储层识别的地质需求;②储层横向变化较大,储层预测难度大;③地震数据成像精度差,深部资料信噪比较低;④原始地震资料噪声干扰较为严重,主体构造存在地震资料空白区。因此,决定在该区进行新一轮的三维OBN地震勘探,以满足构造精细解释、优势储层识别和预测的需要。

VSP技术是对地面地震的有效补充,可以提取地层地球物理参数用来提高地面地震数据处理的精度。因此,VSP井地联合勘探技术在油气田勘探开发中具有良好的应用前景。在对井况数据分析后,甲方决定在B5和B7两口井中开展三维DAS-VSP和三维OBN井-海联采作业。

本次三维OBN勘探的地质目标是:①求取准确的纵波速度资料,包括平均速度和层速度等;②准确标定各地震地质反射层;③提供地层吸收衰减参数Q、真振幅恢复因子TAR、反褶积参数;④进行海上OBN资料多次波分析;⑤精确提取用于偏移的各向异性参数、井驱参数,预期处理分辨率提高5～10Hz;⑥获取井旁高精度高分辨率成像数据。

2 井-海三维OBN和三维DAS-VSP数据联合采集

目标勘探区域位于中国东海平湖油气田。工区海底较为平坦,水深为80～100m。工区内有2个作业平台区,分别为八角亭、放鹤亭主平台。为了获取高质量的OBN地震数据,用于深层油气藏精细刻画、岩性油气藏物性预测、裂缝性油气藏裂缝分布、在浅浮雕结构中准确评价剩余油气藏构造,对数据的采集要求为:①目的层地震反射能量强,信噪比和分辨率满足地质目标要求;②提高深部数据的能量和信噪比;③优化平台区周边OBN勘探设计,减少地震数据空白区。图1显示了三维OBN和三维DAS-VSP数据井-海联合采集的炮点分布情况。

图1 三维OBN和三维DAS-VSP数据井-海联合采集炮点分布(a)和OBN气枪震源条带激发采集阶段计划(b)

表1为三维DAS-VSP观测系统设计参数。只对井源距4000m矩形范围内的OBN炮点进行接收,激发参数与三维OBN地震数据采集范围内的OBN激发参数相同。主要激发参数包括震动频率、时长等,分别以B5和B7井为中心。

表1 三维DAS-VSP数据采集观测系统设计参数

表2是井下铠装光缆接收参数,表3是各井三维DAS-VSP数据采集的接收炮数。图2是本三维OBN数据采集项目实际气枪震源分布和全覆盖次数。

表2 井下铠装光缆接收参数

表3 OBN和2台DAS仪器接收炮数

图2 实际三维OBN数据采集项目气枪震源分布(a)和全覆盖次数(b)

对于用铠装光缆进行DAS-VSP数据的采集,井身结构和井轨迹对于套管内铠装光缆的布设和后续的三维DAS-VSP数据处理十分重要。图3显示了B7井和B5井的井下套管结构和铠装光缆(蓝色线)布设示意。两口井均为斜井并有多层套管。B7井水深88m,完井井深4270m,生产封隔器安装在3595m处。B5井水深81.2m,完井井深4270m,目的层最大井斜11.4°～36.4°。

图3 B7井(a)和B5井(b)的井下套管结构以及套管内铠装光缆(蓝色实线)布设示意

三维DAS-VSP和三维OBN数据的联合采集使用相同的激发震源在井下和海底同时获取信息。在B5井和B7井的套管内部署了两条带有多根耐高温单模光纤的自耦合铠装光缆。B5井和B7井的井下铠装光缆长度分别为3085m和3537m。每条铠装光缆连接到每个平台上的两台uDAS®调制解调仪器,以同时记录两套三维DAS-VSP数据,并将两套数据叠加在一起以进行增强处理。图4是用于三维DAS-VSP数据采集的uDAS®调制解调仪器、数据存储硬盘阵列和井口旁的铠装光缆卷筒的实物照片。从海底到钻孔底部的井下铠装光缆进行全井段三维DAS-VSP数据接收,输出数据的级间距为1m。铠装光缆在B5井内记录了62473次海面气枪震源的激发信号,在B7井内记录了48411次海面气枪震源的激发信号。在不影响OBN数据采集的情况下,24小时连续记录三维DAS-VSP数据,并在三维DAS-VSP数据采集完成后根据每个气枪震源激发的GPS时间对uDAS®调制解调仪器连续记录的三维DAS-VSP数据进行切分。

图4 三维DAS-VSP数据采集所使用的uDAS®调制解调系统(a)、数据存储硬盘阵列(b)和井口旁的铠装光缆卷筒(c)的实物照片

3 三维DAS-VSP数据处理

常规三维VSP数据成像处理步骤包括:观测系统定义、预处理、初至拾取、静校正、振幅补偿、反褶积、上下行波场分离、速度分析和成像处理[8]。此次井-海联采的三维DAS-VSP数据,利用组合去噪处理、子波反褶积、矢量波场分离、各向异性速度建模、地震偏移速度场校正、上行反射波成像、共接收点道集优化成像、下行多次波成像等VSP数据处理方法,对三维DAS-VSP数据进行了成像处理。具体三维DAS-VSP数据处理的主要步骤为:①振铃噪声压制,压制振铃噪声和异常幅度干扰,为上下行波场分离奠定基础;②信号反褶积,利用下行波场得到的子波进行气泡压制和子波零相位处理;③波场分离,将上行波场和下行波场分开,提高用于反射波成像的上行波场数据的信噪比,进行目标波场分析或成像;④浅水区多次波压制,压制水层多次波,突出地层真实反射特征;⑤三维DAS-VSP数据偏移,利用角度域高斯束偏移方法以及倾角约束的积分法偏移对三维DAS-VSP数据进行偏移处理,得到高精度的三维DAS-VSP数据成像结果。

原始三维DAS-VSP数据质量分析结果见图5。包括主频率、信噪比、能量和背景噪声水平。结果表明:①从近到远的单次激发原始数据能量强、信噪比高、主频高、频带宽,低频能量丰富,初至清晰、连续且容易拾取,整体数据质量好;②三维DAS-VSP数据存在不同程度的背景噪声、缆波干扰(振铃噪声),部分井段因海水涌动不可避免存在干扰;③三维DAS-VSP数据波场丰富,多次波较为发育;④铠装光缆多次下井采集的三维DAS-VSP数据无差异,说明采集设备的稳定性和参数设置没有异常;⑤全井高密度三维DAS-VSP数据满足井中地震参数计算和井周三维DAS-VSP精确成像要求。

图5 B7井(a)和B5井(b)的原始三维DAS-VSP数据质量分析结果(主频、信噪比、能量和背景噪声水平)

DAS-VSP勘探具有高密度、全井段覆盖等优点,但由于部分井段铠装光缆与套管壁耦合不良,光缆的缆波(振铃噪声)较大。平湖三维DAS-VSP数据振铃噪声去除前、后结果对比如图6所示。从图6a可以看出,振铃噪声的视速度在时域上接近于有效信号,在频谱中表现为周期性脉冲信号,在自相关谱中呈周期性。传统的去噪方法难以有效去除此类噪声。

图6 三维DAS-VSP数据振铃噪声去除效果对比a 原始单炮记录; b 振铃噪声去除后的单炮记录; c 去除的振铃噪声

本文提出了一种基于τ-p变换的振铃噪声衰减技术。τ-p变换是近几年发展起来的一种新的处理方法,是平面波场分离的方法。它是在炮检距轴线上应用线性时差和振幅求和实现波场的平面波分解,这种方法也叫做倾斜叠加。在国外,τ-p变换已被广泛应用于地震资料处理的各个领域,形成了一个τ-p变换处理系统。例如在τ-p变换域分解干涉,收敛直线性波,进行动校正,叠加,重采样,道内插以及双曲速度滤波等。τ-p域对与视速度接近的噪声有压制效果,这是它的优势,常用于多次波压制等。首先,利用τ-p变换将三维DAS-VSP数据变换到τ-p域,然后,利用线性预测方法对周期性振铃噪声进行预测,最后,从原始数据中减去预测出来的振铃噪声。

与地面地震数据相比,VSP数据的波场更为丰富,不仅有上行反射波,还有下行反射波,且其照明范围更广。一般下行反射波包括自由表面多次波、层间多次波等。如图7所示,当有两层反射界面和一个VSP测量井时,S为地面激发震源,R为井下接收检波器。直达波、一次反射波、自由表面多次波和层间多次波分别沿绿线、蓝线、红线和粉色线传播。一次反射波的传播路径是从地面震源点S到反射层的反射点再到井下接收点R。自由表面多次波的传播路径是从地面震源点S到反射层的反射点,再到自由表面(地表)的反射点,然后到井下接收点R。井下接收检波器R不仅可以接收下伏地层的自由表面多次波,还可以接收上覆地层的自由表面多次波。对比反射点的位置可以看出,自由表面多个反射点在地下浅部波阻抗界面上的横向展布较宽,在地下深部波阻抗界面上的横向展布随深度的增加而变窄。经过下行反射波成像处理后,可以实现井下接收检波器R点上方的一次反射波和自由表面多次波成像,再加上井下接收检波器R点下方的上行反射波的成像处理,三维DAS-VSP数据的上行波和多次波成像可以使横向成像范围不断拓宽。

图7 VSP测量的各种井中地震波场传播路径示意

图6b显示了该数据中去除振铃噪声的效果。采用τ-p变换的振铃噪声衰减技术后,振铃噪声得到了有效的衰减和压制,去除后的噪声数据(图6c)里已经看不到有效的反射波信息。

自由表面多次波逆时偏移法已广泛应用于三维DAS-VSP数据成像[9-12]。FERGUSON等[13]从标量波动方程推导出单程波场延拓的傅里叶积分方程,给出了角频域的显式表达式。基于一次反射波单程叠前深度偏移方法,结合自由表面多次波的传播特性,只需修改波场延拓法即可实现多次波成像。自由表面多次波单程波场连续成像方法可描述为以下几个步骤。

1) 将VSP数据的Z分量重新排列到共检波点道集(CRG)域数据中。

2) 将震源子波向前延拓到表面。在井中的接收点R处设置震源子波。震源子波的波场以步长Δz逐层向上延拓,从z到z-Δz的延拓公式为:

ψs(x,z+Δz,ω)=ψs(x,z,ω)e-ikzΔz

(1)

式中:ψs(x,z,ω)为z深度角频域的震源波场;kz为z方向的波数。这里,波场延拓算子可以是分步傅里叶算子、傅里叶有限差分算子或广义屏幕算子等。

3) 反方向继续延拓一步。CRG波场根据公式(2)从表面以步长Δz逐层向下延拓。

ψr(x,z+Δz,ω)=ψr(x,z,ω)eikzΔz

(2)

式中:ψr(x,z,ω)为自由表面多次波在深度z处角频域的波场;kz为z方向的波数。

4) 继续向前延拓一步。第2步的震源子波波场根据公式(3)从表面逐层以步长Δz向下扩展。

ψr(x,z-Δz,ω)=ψr(x,z,ω)eikzΔz

(3)

式中:ψr(x,z,ω)为自由表面多次波在深度z处角频域的波场;kz为z方向的波数。

5) 提取成像值。将第3步和第4步扩展的波场进行互相关,取零时成像值。重复步骤3),4)和5),直到达到最大成像深度。

三维DAS-VSP数据下行多次反射波成像处理技术流程如图8所示。

图8 三维DAS-VSP数据下行多次反射波成像处理技术流程

由于B5和B7井均为斜井,震源距平台500m,三维DAS-VSP上行反射波成像孔径和覆盖面积沿井斜方向非常有限,几乎无法获取井底储层附近的构造成像(图9和图10)。采用下行多次波进行三维DAS-VSP数据成像,可以大幅度增加成像孔径和覆盖面积。根据三维VSP下行多次波的特点,利用三维VSP下行多次波成像技术,大大扩展了三维DAS-VSP数据的成像范围,提高了三维DAS-VSP数据的整体成像质量和水平[14]。

图9 B5井三维DAS-VSP数据成像比较a 上行反射波成像; b 下行多次反射波成像

图10 B7井三维DAS-VSP数据成像比较a 上行反射波成像; b 下行多次反射波成像

下行多次波成像分辨率较低,因此引入下行多次波Q偏移来减小地层吸收衰减系数的影响,提高多次波成像分辨率。由图9和图10可以看出,B5井和B7井的三维DAS-VSP数据成像结果均有两种,一种是基于上行反射波数据的叠前深度偏移成像,另一种是基于下行多次波数据的叠前深度偏移成像。最终的偏移速度模型采用的是利用三维DAS-VSP数据计算出的各向异性系数更新后的各向异性速度模型[15]。图9a和图10a中的蓝色实线是沿B5和B7井的井轨迹采集的三维DAS-VSP数据的上行反射波数据的偏移成像范围。图9b和图10b是利用三维DAS-VSP数据的下行多次反射波数据的偏移成像,其结果有效地拓展了三维DAS-VSP数据的成像范围。

图11为B7井平湖三维DAS-VSP下行多次反射波成像结果的立体展示。该三维成像数据体可用于指导三维构造的精细解释并生成三维属性数据体,用于井周围的储层表征。由于下行多次反射波成像扩大了覆盖范围,该成像结果可用来获得井周围详细的储层构造图并追踪储层中的流体分布。

图11 B7井平湖三维DAS-VSP下行多次反射波成像结果的三维立体显示

图12为B7井周围早期的三维OBC数据成像、新采集的三维DAS-VSP上行波成像和三维DAS-VSP下行多次反射波成像对比。图12中的蓝色虚线是三维DAS-VSP上行反射波数据成像覆盖区域。图12c中的红色虚线是三维DAS-VSP下行多次反射波成像的覆盖区域,横向成像范围明显得到扩大。当这个三维DAS-VSP成像数据镶嵌到三维OBN数据成像数据体中时,大多数反射界面都匹配良好。

图12 B7井周围早期三维OBC数据成像(a)、新采集的三维DAS-VSP上行反射波数据成像(b)和三维DAS-VSP下行多次反射波数据成像(c)结果

图13显示了B5井周围早期OBC数据成像、新采集的三维OBN数据成像和三维DAS-VSP下行多次反射波成像的比较结果。与早期三维OBC数据成像结果相比,三维DAS-VSP下行多次反射波数据的成像(图13c)显著地改进了早期OBC数据的成像质量。在新采集的OBN数据成像和三维DAS-VSP数据成像中,与早期三维OBC数据成像相比,储层边界的刻画变得更加容易和准确。

4 地质成果解释

地质成果包括B5和B7井周边的储层构造精细解释、储层边界追踪和流体分布预测等。图14 是基于精确VSP时深关系约束的井震地层识别和标定。图14a显示了与VSP剖面和OBN剖面相对应的合成记录和走廊叠加;图14b是基于VSP数据和合成记录的走廊叠加的地震地层的识别与标定。这个识别过程是在井轨迹处对DAS-VSP数据的成像剖面和OBN数据成像剖面的井震关系标定,其良好的标定结果使我们有信心使用三维DAS-VSP数据进行井周围详细地质构造解释和储层预测。

图14 基于精确VSP时深关系约束的联合井震地层标定a 与VSP剖面和OBN剖面相对应的合成记录和走廊叠加(井-震关系); b 基于VSP数据和合成记录的走廊叠加与地震地层的识别和标定

工区早期OBC地震资料与人工合成记录的标定效果不是很好。基于多井零井源距数据准确的时深关系约束,有效提高了标定精度。通过横向井间比对,证实了地震层位识别的一致性,为工区地震层位解释奠定了基础。

图15是3D DAS-VSP数据成像和OBN数据成像的时间切片以及曲率属性切片(2800ms)对比。从时间切片上(图15a)可以看出,3D DAS-VSP成像在北部B5井区西侧地层反射内部清晰,东侧地层特征与OBN数据基本一致;南部B7井区,西侧地层细节更为丰富。从曲率属性切片对比(图15b)看,构造特征整体相近,均能清晰刻画北部主断裂形态,局部小断层存在一定差异;南侧B7井区3D VSP资料信噪比高,能够更为准确刻画小断裂的发育情况。该资料为精细储层预测、有利断块精细刻画奠定了基础,为油气藏开发后期加密井或侧钻井部署提供可靠的资料。

图15 3D DAS-VSP数据成像和OBN数据成像的2800ms时间切片(a)以及曲率属性切片(b)对比

图16显示了多井三维DAS-VSP联合数据反演剖面及平湖组砂体沿四井连井剖面的分布情况。对于反演剖面,利用三维DAS-VSP数据反演振幅特性和波形特征,进而分析预测砂体的横向厚度及分布。反演剖面中的红色储层表示最高的油气饱和度,黄色表示较低的油气饱和度。砂体的三维横向厚度分布以及砂体内油气饱和度预测对于工区优化提高油气采收率开发方案以提高该成熟油气田的产量具有非常重要的价值。

图16 多井联合三维DAS-VSP数据反演剖面(a)和平湖组砂体沿四井连井剖面(b)的分布情况

5 结论

DAS-VSP测量越来越被认为是用于井中地震数据采集的井下三分量地震检波器阵列的可行替代方案。井地联合地震勘探技术是地面地震与井中地震数据同时采集相结合形成的三维立体地震勘探方法。利用与地面地震资料同时采集的零偏移距VSP数据,可以获得准确的时深关系、地层速度、球面扩散补偿因子TAR值、吸收衰减因子Q值和各向异性参数。这些参数可以应用于地面地震数据的井驱处理,显著提高地震数据的保幅特性、分辨率和成像精度。在东海某工区进行的三维OBN数据采集的同时用井下铠装光缆同步采集了三维DAS-VSP数据,并开展了三维DAS-VSP数据的成像处理和综合解释,其结果表明,与早期三维OBC数据成像结果相比,本次新采集的OBN数据成像结果和三维DAS-VSP下行多次反射波成像结果显示出显著的成像质量改进,目标地层和储层边界的追踪变得更容易和更加可靠。使用三维DAS-VSP数据可以帮助我们进行更为详细的地质构造解释和井周围的储层构造成像。基于多井三维DAS-VSP数据准确的时深关系约束,有效提高了地层和储层的标定精度。通过井间横向比对,证实了地震层位识别的一致性,为研究区地震层位解释奠定了基础。

利用多井三维DAS-VSP联合数据的属性反演剖面及储层砂体沿连井剖面的振幅特性和波形特征,可以分析预测工区内井周围砂体的横向厚度变化和砂体的油气饱和度,进而规划提高油气采收率的优化开发方案以提高该成熟油气田的产量。

致谢：感谢上海石油天然气有限公司允许发布数据和BGP与中油奥博(成都)科技有限公司为本项目在数据采集和处理阶段提供支持和指导,以及BGP“井中DAS-VSP关键技术研究”(06-01-2021)科研项目和中国石油研究项目“井下光纤智能监测项目(2020F-44)”的大力支持。