0 引言
目前,随着我国含油气盆地勘探程度的不断提高,勘探难度逐渐加大,复杂构造、复杂岩性、深层/超深层油气藏已逐渐成为主要勘探对象[1 ] ,勘探的地质目标具有“隐、碎、薄、小”的特点,如复杂小断块、致密碎屑岩等[2 -3 ] ,对地震偏移成像技术提出了更高的要求。井中地震具有资料信噪比高、储层识别能力强的优点。此外,它能够适应检波器的任意排列,可减少面元数据体丢失,理论上能够对井周的“隐、碎、薄、小”等储集体实现高精度成像。
井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震。随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] 。但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震。对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] 。在石油地震勘探领域,井中地震技术的研究相对较少。
在实际生产中,钻井的设计受到生产要求或客观因素影响,井型除常规的直井外还存在斜井、曲斜井、水平井等多种情况。与常规的直井相比水平井在非常规油气资源高效勘探和开发上有着重要的作用[17 ] 。而在海洋勘探中,由于作业平台的空间有限,其井型一般都为有角度的斜井或空变的斜井[18 ] 。由于井中地震的激发震源位于井中,不同井型使得震源所处空间位置有较大差异,造成地震波的传播路径存在显著不同,这势必会对最终的成像结果产生影响。因此在应用井中地震技术前,需详细分析井型对成像结果的影响,进而指导实际生产中激发井段的选择。
考虑到应用井中地震技术的目的是对深层复杂构造进行精细成像,对比常用的Kirchhoff积分偏移、单程波偏移、逆时偏移等成像算法,逆时偏移算法具有没有角度限制、适合复杂介质成像的优势。此外,由于震源位于地下,能量微弱受地层吸收衰减作用影响严重,地震波在传播过程中会伴随着能量衰减和波形畸变现象,影响偏移剖面精度。综合以上因素,对于井中地震,应对其进行Q 补偿逆时偏移成像,以而获得高保真度的偏移成像剖面。
早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果。而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力。目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] 。总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点。为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] 。其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果。但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题。基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM。
综上所述,井中地震是一种可以实现对井周“隐、碎、薄、小”地质目标进行精细成像的物探方法。在实际生产中,井型不同对成像结果存在影响,然而,当前井型对于偏移成像的影响具体如何,还没有定性或定量的认识。本文使用黏声逆时偏移方法进行井中地震成像,并基于Zhu等[30 ] 提出的分数阶黏声波方程进行数值模拟,通过在理论模型中进行直井、斜井、曲斜井、水平井与多井的成像试算,探究井型对偏移质量的影响。数值实验结果阐明了井型与成像剖面间的影响关系,为井中地震采集系统设计提供了理论支撑。
1 井中黏声逆时偏移基本原理
1.1 井中地震特点
井中地震技术是一种井中激发,地面接收的特殊物探方法。由于井中地震观测系统特殊,导致其与地面地震存在较大差异,地面和井中地震观测系统如图1 所示。从图中可见,与地面地震相比,由于井中震源位于地下且靠近勘探目标,故其激发的地震波减少了一次低降速带损耗,规避了面波的产生,有助于采集高分辨率及信噪比的地震信号;同时由于井中地震在地面上不存在震源干扰,因此,在地表可进行大范围检波器布置,有利于进行高密度采集;且可以采用规则化观测系统在地面布置检波器,以降低面元数据体丢失,便于在采集区域获得更丰富的有效信息。
图1
图1
井中地震(a)与地面地震(b)观测系统示意
Fig.1
Schematic diagram of 2D borehole seismic(a) and ground seismic(b) observation system
井中地震由于震源位于井中,因此当震源位于不同井型时,地震波的传播路径也不尽相同。不同井型地震波的传播路径如图2 所示。通过对图2 a、b、c的对比发现,直井与曲斜井和水平井的地震波传播路径主要差异在于偏移距的变化,由于波的传播路径不同故其偏移距也不同。综合以上分析可以发现,井型的变化主要影响了地震波的传播路径,进一步对偏移距产生了影响。考虑到井中地震的偏移距大小与反射点覆盖范围密切相关,因此在理论上井型的变化会对最终成像剖面的成像范围产生影响。
图2
图2
井中地震地震波的传播路径示意
a—直井路径示意;b—曲斜井路径示意;c—水平井路径示意
Fig.2
Propagation path of seismic wave in well
a—schematic diagram of vertical well path;b—schematic diagram of curved inclined shaft path;c—schematic diagram of horizontal well path
尽管井中地震技术具有诸多优点,但其同样存在受地层衰减效应影响大,采集的地震数据层间多次波干扰严重、覆盖范围有限的缺点。总体来说,在理论上井中地震技术是一种适用于井周区域小构造精细勘探的特殊物探方法。
1.2 分数阶黏声波方程
在地震勘探领域,介质的衰减特性通常被认为与频率呈线性关系。因此,在地震资料频带范围内,可以假定Q 与频率无关,即常Q 模型。基于此假设,Kjartansson[28 ] 推导了常Q 模型的速度相位频散和衰减关系:
(1) c p = c 0 ω ω 0 γ
(2) α = t a n π γ 2 ω c p 。
式中:ω 0 是参考角频率;ω 是角频率;c 0 是参考角频率对应的声波速度;c p 是声波速度;γ = t a n - 1 ( 1 / Q ) / π ,是无量纲参数;α 是衰减因子。基于上述常Q 频散关系,Zhu等[30 ] 推导了时间域分数阶黏声波方程:
(3) 1 c 2 ∂ 2 p ∂ t 2 = η - ∇ 2 γ + 1 p + τ ∂ ∂ t - ∇ 2 γ + 1 2 p
(4) η = - c 0 2 γ ω 0 - 2 γ c o s ( π γ ) , τ = - c 0 2 γ - 1 ω 0 - 2 γ s i n ( π γ ) ,
式中:p 是地震波场;∇ 2 是Laplace 算子;c 2 = c 0 2 c o s 2 ( π γ / 2 ) 。
对于方程的稳定性问题,本文通过对方程(3)的左侧进行广义傅里叶变换,并对二阶时导数进行中心差分。下一时刻的波场p n + 1 可以由当前时刻的波场p 和前一时刻的波场p n - 1 表示如下:
(5) p n + 1 = 2 p n - p n - 1 + Δ t 2 c 0 2 η k 2 γ + 2 p n + Δ t c 0 2 τ k 2 γ + 1 p n - p n - 1
(6) p n + 1 p n = 2 + a Δ t 2 + b Δ t - 1 - b Δ t 1 0 p n p n - 1
(7) a = c 0 2 γ η k 2 γ + 2 , b = c 0 2 τ k 2 γ + 1 。
本文使用特征值法分析介质波动方程的稳定条件,方程若保持稳定则其矩阵的特征值必须小于或等于1,故通过对式(6)的特征值求取可得,分数阶方程的稳定性条件为:
(8) Δ t ≤ - b - - 2 a a = t a n ( π γ ) c 0 k + 2 ω 0 γ ( c 0 k ) γ + 1 c o s ( π γ ) 。
为了保证解对所有波都是稳定的,c 0 = c m a x 是最大速度,k = π / Δ x 是奈奎斯特空间波数,此时γ = t a n - 1 1 / Q m i n / π 。
为了便于实现基于上述方程的黏声逆时偏移,本文使用统一的黏声波方程进行波场的正向和反向延拓,具体如下:
(9) 1 c 2 ∂ 2 p ∂ t 2 = ∇ 2 p + β 1 [ η ( - ∇ 2 ) γ + 1 - ∇ 2 ] p + β 2 τ ∂ ∂ t ( - ∇ 2 ) γ + 1 2 p
方程(9)右端第二项和第三项分别为表示相位频散和振幅衰减。其中,β 1 可以取1和0,β 2 可以取1、0和-1。对于相位频散项,β 1 取1表示包含相位频散,β 1 取0表示无相位频散。对于振幅衰减项,β 2 取1表示振幅衰减,β 2 取0表示无振幅衰减,β 2 取-1表示振幅衰减补偿。
基于上述分析,Q 补偿逆时偏移需包含相位频散项,而且需对振幅衰减进行补偿,也即β 1 = 1 ,β 2 = - 1 ,则方程(9)可重写为:
(10) 1 c 2 ∂ 2 p ∂ t 2 = ∇ 2 p + [ η ( - ∇ 2 ) γ + 1 - ∇ 2 ] p - τ ∂ ∂ t ( - ∇ 2 ) γ + 1 2 p 。
本文利用方程(10)进行波场的正向和反向延拓,校正地震波传播过程中的振幅衰减和相位频散。
1.3 Q 补偿逆时偏移原理
叠前逆时偏移分为3个步骤:①震源波场正传;②检波点波场反传;③应用成像条件成像。在波场模拟的过程中本文选用Liu等[40 ] 的混合吸收边界作为边界条件,该方法将计算区域分为内部区、过渡区和边界区。在进行波场模拟时,内部区中用双程波方程求解,边界区用单程波方程求解,过渡区波场由单程波和双程波的波场加权得到;在成像条件的选取上,本文选用了在成像过程中同时利用了震源波场和检波点波场的信息,可以对所有波进行成像的互相关成像条件。互相关成像条件可表示为:
(11) I ( x ) = ∫ 0 T S ( x , t ) R ( x , t ) d t
式中:I ( x ) 表示成像结果;S ( x , t ) 和R ( x , t ) 分别表示t 时刻成像空间x 处的震源波场和检波点波场;x = ( x , y , z ) 表示成像点空间位置矢量;T 为最大记录时间。
图3 为井中地震波场正向延拓及逆时偏移原理示意。其中图3 a、b分别表示无衰减和衰减介质波场正向延拓过程。对于衰减介质,波场正向延拓过程中振幅的衰减情况可以表示为:
(12) R ^ ( x , t ) = R ( x , t ) e - α L D e - α L U
图3
图3
井中地震波场延拓及逆时偏移原理示意
a—无衰减介质波场正向延拓;b—衰减介质波场正向延拓;c—声波逆时偏移;d—Q 补偿逆时偏移
Fig.3
Principles of wavefield propagation and reverse time migration in boreholes seismic
a—forward continuation of wave field in non attenuating medium;b—forward continuation of wave field in attenuating medium;c—acoustic inverse time migration;d—Q compensation inverse time migration
由于实际地层是一种黏弹性介质,当地震波在黏弹性介质中传播时,波场的振幅和相位会发生变化,最终影响成像剖面。图3 c、d分别表示声波逆时偏移和Q 补偿逆时偏移的基本原理。对于衰减介质,需在震源波场正向延拓和检波点波场反向延拓过程中进行补偿,基于图3 b所示的地震波场衰减规律,Q 补偿逆时偏移过程中,任意成像点处的震源波场和检波点波场可分别使用算子e + α L D 和e + α L U 进行衰减补偿,具体可以表示为:
(13) S C ( x , t ) = S ( x , t ) e + α L D
(14) R C ( x , t ) = R ^ ( x , t ) e + α L U
利用补偿后的震源波场S C ( x , t ) 和检波点波场R C ( x , t ) ,应用零延迟互相关条件实现成像,具体为:
(15) I C ( x ) = ∫ 0 T S C ( x , t ) R C ( x , t ) d t
结合式(12)、(14)和(15),由零延迟互相关成像条件得到的补偿后图像I C ( x ) 理论上等于声波逆时偏移得到的图像I ( x ) ,即介质的吸收衰减效应得到完全补偿。
当不考虑地层的黏滞性时波场正向延拓和逆时延拓在参数选取合理的情况下都是稳定的。但当地震波在黏滞性介质中传播时,由于波场正传过程是一个能量衰减过程,故该过程稳定。但在逆时延拓的过程中,黏滞波动方程会对震源和检波点波场进行补偿,使其能量呈指数型增大,与此同时波场中高频噪声也被呈指数放大,因此出现高频不稳定现象。为保证波场的稳定[30 ] ,通过在波数域设置低通滤波器将高频噪声去除。与此同时,本文考虑到在应用互相关成像条件进行逆时偏移成像时,成像过程中存在低波数噪声干扰现象。因此本文通过在成像剖面应用Laplace滤波的方式来压制低波数噪声。
在偏移成像过程中由于井中地震特殊的观测方式,导致其成像范围有限,有效成像范围近似为“拱形”。图4 a为水平层状介质中,井中地震射线路径示意,从图中可以发现,对于井中地震数据,不同采样点对应的反射点位置在水平方向和深度方向均有变化,反射点轨迹成“拱形”。随深度增加,反射点位置逐渐远离震源,到达最深位置时,反射点与震源相距最大,水平距离为偏移距的一半。由此可见,当震源深度一定时偏移距越大,反射点水平方向覆盖范围越大;图4 b为不同激发深度覆盖范围分布,可以看到当震源处于不同深度时,其覆盖范围也有所差异。当检波器排列固定时,震源深度越大,其覆盖范围越小,最大覆盖范围为一半炮检距,此时震源位于井口[13 ] 。
图4
图4
井中地震射线路径(a)及不同激发深度覆盖范围(b)示意
Fig.4
Seismic ray path(a) and coverage map with different excitation depth(b) in well
2 数值模拟
2.1 井中地震数值模拟
本文选用砂泥岩薄互层模型应用声波和黏声波方程进行数值试验,图5 为所用模型,大小为242×600,深度方向采样间隔为5 m,水平方向采样间隔为10 m。采用主频为20 Hz的雷克子波作为震源,记录时长为5 s,采样间隔为0.5 ms,共设置242个检波器均匀分布于地表,道间距为10 m。为了压制震源引起的强照明现象,本文采用设置衰减因子的方式对震源处能量进行压制,而由于井口距检波器位置较近,衰减因子的影响导致成像剖面井口处存在一定的能量不均匀现象。
图5
图5
砂泥岩薄互层速度模型(a)和Q 模型(b)
Fig.5
Veolocity model(a) and Q model(b) of sand-shale thin interbed
图6 为井中地震逆时偏移剖面,激发井段长度为1 000 m,炮数20炮,炮间隔50 m,由水平方向1 200 m深度300 m处自上而下依次激发,图中红色竖线代表井位,图6 a~c分别为声波偏移剖面,黏声未补偿偏移剖面,黏声补偿后偏移剖面。本文选用不考虑衰减情况的声波逆时偏移结果(图6 a)作为参考剖面,通过对比未补偿和补偿后偏移结果可见,成像剖面与理论预期相符。图6 b为考虑了黏弹性介质的衰减影响,但未进行补偿的逆时偏移剖面,与参考剖面相比,该剖面能量明显减弱,且仅可对浅层进行成像,深度2km以下的薄互层区域难以成像,整体成像效果极差。图6 c为补偿后的逆时偏移剖面,与图6 b对比可以看出,补偿后的剖面消除了地层衰减作用对成像剖面质量的影响,解决了成像剖面能量弱的问题,同时对模型的浅、中、深层均进行了良好的刻画,结果与参考剖面相符。图7 为与图6 对应的薄互层模型局部偏移剖面,通过图7 b、c的对比可以看出:地层的吸收衰减效应对薄互层影响严重,当未进行衰减补偿时,薄互层处能量弱,层位难以识别;补偿后的剖面能量得到恢复,层位刻画清晰与参考剖面(声波局部偏移剖面)结果相符。
图6
图6
薄互层模型逆时偏移剖面
a—声波偏移剖面;b—黏声未补偿偏移剖面;c—黏声补偿后偏移剖面
Fig.6
Reverse time migration profile of thin interbedded model
a—acoustic migration profile;b—viscoelastic uncompensated migration profile;c—migration profile after viscoelastic compensation
图7
图7
薄互层模型逆时偏移剖面局部放大
a—声波偏移剖面;b—黏声未补偿偏移剖面;c—黏声补偿后偏移剖面
Fig.7
Partial enlarged view of reverse time migration profile of thin interbedded model
a—acoustic migration profile;b—viscoelastic uncompensated migration profile;c—migration profile after viscoelastic compensation
通过对比图6 a~c和图7 a~c可得出以下结论:地层的吸收衰减效应会严重影响偏移剖面的成像质量,使目的层难以识别;采用稳定Q 补偿的井中逆时偏移成像方法能够有效地补偿介质吸收衰减对成像结果造成的影响,改善成像质量。
2.2 井型影响因素分析
由于井中震源所处位置随井位一同变化,考虑到实际生产中存在直井、斜井、曲斜井、水平井等多种情况。因此,本文仍选用砂泥岩薄互层模型基于Q 补偿黏声波方程,通过水平井、斜井和多井情况的数值模拟试验以探究当震源处于不同井型情况下对井中地震逆时偏移成像的影响。图8 为补偿后的井中地震逆时偏移剖面, 由于井中观测系统覆盖区域呈“拱形”,成像范围有限,因此本文在成像剖面中将无法成像的区域切除,仅保留有效成像范围,并用黑色曲线圈出。剖面中红色竖线代表井位,黑色曲线代表其成像范围,本文将该成像剖面作为参考剖面用以检测井型变化对井中地震成像结果的影响。
图8
图8
薄互层模型逆时偏移参考剖面
Fig.8
Reverse time migration reference profile of thin interbedded model
为探究当存在水平激发段时对偏移结果的影响,选用一个简单的井位布置形式,即在直井后增加一段水平井段,模拟实际生产中的激发井存在水平激发井段的情况。图9 为补偿后的井中地震偏移剖面,直井井段中震源布置方式与图6 中相同,在水平井段中震源炮间隔50 m,由左向右依次激发。图9 a为水平井段长度为500 m的偏移剖面,炮数10炮;图9 b为水平井段长度为1000 m的偏移剖面,炮数20炮。通过与参考剖面的对比可见,当存在水平井段时,偏移剖面照明范围有一定提升,但增长有限;通过图9 a、b成像剖面对比可以看出,当水平井段长度增加时,其照明范围也随之扩大,但对薄互层处影响有限。当存在水平激发井段时,成像效果改善有限,仅在照明范围方面有一定变化。考虑实际生产成本及震源放置的难易程度问题,本文认为在井中地震资料的采集过程中,水平井段并不适合作为激发井段。
图9
图9
薄互层模型水平井逆时偏移剖面
a—水平井段500 m;b—水平井段1 000 m
Fig.9
Reverse time migration profile of horizontal well with thin interbed model
a—horizontal well section with 500 m;b—horizontal well section with 1 000 m
为探究倾斜井对偏移成像的影响,本文设计了右侧倾斜井和左侧倾斜井两组数值模拟实验,每组实验均包含3个测井倾斜程度变化的数值测试,激发方式为自上而下依次等间隔激发,炮数均为20炮。
为探究右倾程度对成像效果的影响,分别对斜率为0.5、1以及曲斜井进行数值模拟,得到图10 a、b、c所示的右倾斜井补偿后偏移剖面。其中图10 a为激发井段斜率为0.5的斜井,竖直方向1 000 m,水平方向500 m;图10 b为激发井段斜率为1的斜井,竖直方向1 000 m,水平方向1 000 m;图10 c为曲斜井,激发井段竖直方向1 300 m,水平方向1 000 m。通过与参考剖面的对比可以看出,当存在右侧倾斜激发井段时,成像范围有一定提升,在深层薄互层处成像范围提升极其明显;通过图10 a、b、c的成像剖面对比可以看出,随着激发井段右倾程度增加时,模型右侧的横向探测范围也随之增加,且范围增长明显。
图10
图10
薄互层模型右斜井逆时偏移剖面
a—斜率为0.5偏移剖面;b—斜率为1偏移剖面;c—右斜井偏移剖面
Fig.10
Reverse time migration profile of right inclined well in thin interbedded model
a—offset profile with slope of 0.5;b—offset profile with slope of 1;c—offset profile of right inclined shaft
为探究左倾程度对成像效果的影响,分别对斜率为-0.5、-1以及曲斜井进行数值模拟,得到图11 a、b、c所示的左倾斜井补偿后偏移剖面。其中图11 a激发井段斜率为-0.5的斜井,竖直方向1 000 m,水平方向500 m;图11 b激发井段斜率为-1的斜井,竖直方向1 000 m,水平方向1 000 m;图11 c为曲斜井,激发井段竖直方向1 300 m,水平方向1 000 m。将该模型与参考剖面对比可以看出,当存在左侧倾斜激发井段时,成像范围有一定提升,但整体提升不明显;通过图11 a、b、c的成像剖面对比可以看出,随着激发井段左倾程度增加时,模型左侧的横向探测范围也随之增加,但范围增长有限。
图11
图11
薄互层模型左斜井逆时偏移剖面
a—斜率为-0.5偏移剖面;b—斜率为-1偏移剖面;c—左斜井偏移剖面
Fig.11
Reverse time migration profile of left inclined well in thin interbedded model
a—offset profile with slope of -0.5;b—offset profile with slope of -1;c—offset profile of left inclined shaft
综合两组倾斜井实验的数值模拟结果,发现斜井激发时相当于在参考井横向上增加了炮点,模型上部分的平层在偏移剖面中成像范围扩大;对模型下部分的薄互层构造成像范围的刻画也有所改变,当激发井段方向与地层倾斜方向相同时,受偏移距的影响,成像范围基本没有变化。当激发井段方向与地层倾斜方向相反时,成像范围扩大明显。
考虑到实际生产中存在的多井联采情况,本文设计双井和三井数值实验,以探究存在多激发井时的偏移成像效果。图12 为补偿后的井中地震偏移剖面,激发井段中震源布置方式与图6 相同,震源起始深度均为300 m。图12 a为双井逆时偏移剖面,两口井分别位于水平方向800 m和1600 m处;图12 b为三井逆时偏移剖面,3口井分别位于水平方向600 m、1 200 m和1 800 m处。通过与参考剖面的成像范围对比可见,当进行多井激发时模型整体的横向探测范围有明显扩大;通过双井和三井的偏移剖面对比可见,三井的成像范围明显大于双井,且深层薄互层构造相对清晰。
图12
图12
薄互层模型多井逆时偏移剖面
a—双井偏移剖面;b—三井偏移剖面
Fig.12
Reverse time migration profile in thin interbed model with multiple wells
a—double well migration profile;b—mitsui migration profile
3 结论
本文使用黏声逆时偏移进行井中地震成像,通过砂泥岩薄互层模型进行数值模拟,对直井、斜井、曲斜井和水平井进行不同井型的数值模拟实验。综合成像结果,通过对不同井型下的成像剖面对比分析可得出以下结论:
1)稳定Q 补偿逆时偏移方法在井中地震技术中有良好的适用性,能够有效地补偿介质吸收衰减造成的影响,并在后续的逆时偏移过程中极大改善了成像的质量。
2)通过薄互层模型的数值模拟结果,发现当激发井段变化时,其成像范围也会随之改变。当激发井段存在倾斜或水平井段时,由于横向上炮点增加,横向探测范围有不同程度提高;当激发井段为斜井,方向与地层倾斜方向相反时,成像范围扩大及其明显。
3)模型测试表明,当存在多口激发井时,井中地震横向探测范围有明显的提升,随激发井数量的增加其探测范围也随之增加,且深层薄互层构造的刻画精度也有一定程度的提高。
参考文献
View Option
[1]
何海清 , 范土芝 , 郭绪杰 , 等 . 中国石油 “十三五” 油气勘探重大成果与 “十四五” 发展战略
[J]. 中国石油勘探 , 2021 , 26 (1 ): 17 -30 .
DOI:10.3969/j.issn.1672-7703.2021.01.002
[本文引用: 1]
“十三五”期间,中国石油天然气股份有限公司(简称中国石油)油气勘探工作在地质理论和技术创新、 油气勘探突破发现和增储等方面取得众多成果。文章在梳理总结“十三五”期间油气勘探成果、理论技术进展的基础 上,提出了“十四五”油气勘探发展思路与战略。“十三五”以来,中国石油坚持创新和资源战略,突出高效勘探, 加大新区新领域风险勘探,加强重点盆地、重点领域集中勘探,强化基础研究、整体研究和领域目标研究;创新和完 善了凹陷区砾岩成藏、浅水三角洲大面积岩性成藏、深层海相碳酸盐岩“三古”控藏、走滑断裂“断溶体”控藏、前 陆冲断带深层成藏、页岩油“四要素”成藏富集、页岩气“三控”成藏富集、深层火山岩“三元”控藏等地质理论与 认识;攻关形成了“双高”和“双复杂”三维地震采集处理解释、深层—超深层钻完井、深层测试、长水平井压裂技 术等为代表的勘探(评价)技术。在四川、塔里木、准噶尔、鄂尔多斯、松辽、渤海湾等大盆地岩性—地层、海相碳 酸盐岩、前陆冲断带、页岩油气、成熟探区、火山岩领域取得了34 项重大勘探成果,油气储量实现历史性长期高位 增长,为石油产量稳中上升、天然气产量快速增长奠定了资源基础。“十四五”期间,中国石油将科学谋划油气勘探 的战略方向和战略重点,推动六大盆地和五大领域常规油气持续规模增储及非常规资源快速发展,加快落实油气规模 增储大场面,积极准备战略接替区和重大接替领域,夯实高质量发展资源基础。
He H Q , Fan T Z , Guo X J , et al . PetroChina: Major achievements in oil and gas exploration during the 13th Five-Year Plan period and development strategy for the 14th Five-Year Plan
[J]. China Petroleum Exploration , 2021 , 26 (1 ): 17 -30 .
[本文引用: 1]
[2]
何登发 , 李德生 , 童晓光 , 等 . 中国沉积盆地油气立体综合勘探论
[J]. 石油与天然气地质 , 2021 , 42 (2 ):265 -284 .
[本文引用: 1]
He D F , Li D S , Tong X G , et al . Integrated 3D hydrocarbon exploration in sedimentary basins of China
[J]. Oil & Gas Geology , 2021 , 42 (2 ): 265 -284 .
[本文引用: 1]
[3]
杨勤勇 , 杨江峰 , 王咸彬 , 等 . 中国石化物探技术新进展及发展方向思考
[J]. 中国石油勘探 , 2021 , 26 (1 ): 121 -130 .
DOI:10.3969/j.issn.1672-7703.2021.01.010
[本文引用: 1]
物探技术作为油气勘探不可或缺的重要手段,为中国油气储量和产量增长提供了技术保障。经过数十年 的发展,中国石油化工股份有限公司(简称中国石化)物探技术形成了具有自身特色的地震资料采集、处理和解释技 术系列。“十三五”以来,在高密度地震勘探、可控震源高效采集、深层—超深层碳酸盐岩储层成像与描述、致密碎 屑岩“甜点”识别、非常规一体化和软硬件研发等技术方面取得了重要进展,有力支撑了中国石化重点领域油气勘探 开发。围绕当前中国石化主要勘探开发领域,详细分析东部断陷盆地碎屑岩、深层—超深层碳酸盐岩、中西部致密碎 屑岩、非常规、火成岩、山前带、海域七大领域面临的地质问题,明确了具有针对性的物探技术需求,结合国内外物 探行业技术发展趋势,提出了未来中国石化在支撑当前技术攻关、引领未来技术研发、“卡脖子”软件装备开发3 个 方面的发展方向,以推动中国石化物探技术的快速发展,对促进中国物探行业高质量发展具有重要意义。
Yang Q Y , Yang J F , Wang X B , et al . Sinopec:Progress and development direction of geophysical prospecting technology
[J]. China Petroleum Exploration , 2021 , 26 (1 ): 121 -130 .
[本文引用: 1]
[4]
Weatherby B B . Method of making sub-surface determinations
[P]. US , US2062151 ,1936 -11 -24.
[本文引用: 1]
[5]
Deily F H , Dareing D W , Paff G H , et al . Downhole measurements of drill string forces and motions
[J]. Journal of Engineering for Industry , 1968 , 90 (2 ):217 -225 .
DOI:10.1115/1.3604617
URL
[本文引用: 1]
[6]
Squire W D , Alsup J M . Linear signal processing and Ultrasonic transversal filters
[J]. IEEE Transactions on Microwave Theory & Techniques , 1969 , 17 (11 ):1020 -1040 .
DOI:10.1109/TMTT.1969.1127091
URL
[本文引用: 1]
[7]
Haldorsen J , Miller D E , Walsh J J . Walk-away VSP using drill noise as a source
[J]. Geophysics , 1995 , 60 (4 ):978 .
DOI:10.1190/1.1443863
URL
[本文引用: 1]
[8]
杨微 . 随钻地震信号检测方法研究 [D]. 北京 : 中国地震局地球物理研究所 , 2007 .
[本文引用: 1]
Yang W . Single detection of the drill bit seismic wave whlie drilling [D] .Beijing : Institute of Geophysics ,China Earthquake Administration, 2007 .
[本文引用: 1]
[9]
吕海川 , 朱伟伦 , 贾衡天 , 等 . 随钻VSP测量中地震波场的数值模拟
[J]. 石油机械 , 2017 , 45 (2 ):10 -12 ,44.
[本文引用: 2]
Lyu H C , Zhu W L , Jia H T , et al . Numerical simulation of seismic wave field in VSP-WD
[J]. China Petroleum Machinery , 2017 , 45 (2 ):10 -12 ,44.
[本文引用: 2]
[10]
Liang Z H . Wavefield processing of reverse VSP data
[J]. SEG Technical Program Expanded Abstracts , 1991 , 10 (1 ):1646 .
[本文引用: 1]
[11]
胡建平 . 变偏移距VSP射线追踪模型
[J]. 西安工程学院学报 , 1998 , 20 (S1 ):10 -13 .
[本文引用: 1]
Hu J P . Walkaway VSP ray tracing model
[J] .Journal of Earth Sciences and Environment , 1998 , 20 (S1 ):10 -13 .
[本文引用: 1]
[12]
朱龙生 . 多方位角逆VSP层析成像 [D]. 西安 : 长安大学 , 2003 .
[本文引用: 1]
Zhu L S . Multi-azimuth Inverse VSP tomography [D] .Xi'an : Chang’an University , 2003 .
[本文引用: 1]
[13]
胡明顺 . 煤层气RVSP地震勘探成像方法研究 [D]. 北京 : 中国矿业大学 , 2013 .
[本文引用: 2]
Hu M S . Study on RVSP Seismic Imaging for Coalbed Methane Exploration [D] .Beijing : China University of Mining and Technology , 2013 .
[本文引用: 2]
[14]
金红娣 , 潘冬明 , 杨光 . RVSP等效地面处理方法研究
[J]. 地球物理学进展 , 2015 , 30 (2 ):641 -649 .
[本文引用: 1]
Jin H D , Pan D M , Yang G . Study on equivalent surface data processing method in RVSP
[J] .Progress in Geophysics , 2015 , 30 (2 ):641 -649 .
[本文引用: 1]
[15]
张辉 . 碳酸岩裸露区煤田RVSP勘探技术研究与应用 [D]. 北京 : 中国矿业大学 , 2018 .
[本文引用: 1]
Zhang H . Research and application of RVSP exploration technology in Carbonate exposed coalfield [D] .Beijing : China University of Mining and Technology , 2018 .
[本文引用: 1]
[16]
Hu M S , Pan D M , Zhou F B , et al . Multi-hole joint acquisition of a 3D-RVSP in a karst area:Case study in the Wulunshan Coal Field,China
[J]. Appl. Geophys. , 2020 , 17 :37 -53 .
DOI:10.1007/s11770-020-0808-8
URL
[本文引用: 1]
[17]
邹才能 , 张国生 , 杨智 , 等 . 非常规油气概念, 特征, 潜力及技术——兼论非常规油气地质学
[J]. 石油勘探与开发 , 2013 , 40 (4 ): 385 -399 ,454.
[本文引用: 1]
Zou C N , Zhang G S , Yang Z , et al . Geological concepts, characteristics, resource potential and key techniques of unconventional hydrocarbon:On unconventional petroleum geology
[J]. Petroleum Exploration and Development , 2013 , 40 (4 ): 385 -399 ,454.
[本文引用: 1]
[18]
范廷恩 , 余连勇 , 杨飞龙 , 等 . 斜井 VSP 高斯射线束正演方法
[J]. 中国海上油气 , 2014 , 26 (5 ):30 -35 .
[本文引用: 1]
Fan T E , Yu L Y , Yang F L , et al . A method of Gaussian beam forward modeling in deviated-well VSP
[J]. China Offshore Oil and Gas , 2014 , 26 (5 ):30 -35 .
[本文引用: 1]
[19]
刘财 , 冯晅 , 张瑾 . 稳定的迭代法反Q滤波
[J]. 石油地球物理勘探 , 2013 , 48 (6 ):890 -895 .
[本文引用: 1]
Liu C , Feng X , Zhang J . A stable inverse Qfiltering using the iterative filtering method
[J]. OGP , 2013 , 48 (6 ):890 -895 .
[本文引用: 1]
[20]
吴吉忠 , 杨晓利 , 龙洋 . 一种稳定高效的等效Q值反Q滤波算法及应用
[J]. 石油地球物理勘探 , 2016 , 51 (1 ): 63 -70 .
[本文引用: 1]
Wu J Z , Yang X L , Long Y . A robust approach of inverse Q filtering with equivalent Q
[J]. OGP , 2016 , 51 (1 ):63 -70 .
[本文引用: 1]
[21]
Hale D . Q-adaptive deconvolution
[J]. SEG Technical Program Expanded Abstracts , 1982 :82 -83 .
[本文引用: 1]
[22]
Hargreaves N D . Similarity and the inverse Q filter:Some simple algorithms for inverse Q filtering
[J]. Geophysics , 1992 , 57 (7 ): 944 -947 .
DOI:10.1190/1.1443307
URL
[本文引用: 1]
[23]
Deng F , McMechan G A . Viscoelastic true-amplitude prestack reverse-time depth migration
[J]. Geophysics , 2008 , 73 (4 ):S143 -S155.
DOI:10.1190/1.2938083
URL
[本文引用: 1]
[24]
Dutta G , Schuster G T . Attenuation compensation for least-squares reverse time migration using the viscoacoustic-wave equation
[J]. Geophysics , 2014 , 79 (6 ): S251 -S262.
DOI:10.1190/geo2013-0414.1
URL
[本文引用: 1]
[25]
Bai J , Chen G , Yingst D , et al . Attenuation compensation in viscoacoustic reverse time migration
[J]. SEG Technical Program Expanded Abstracts , 2013 :3825 -3830 .
[本文引用: 1]
[26]
Tian K , Huang J , Bu C , et al . Viscoacoustic reverse time migration by adding a regularization term
[C]// New Orleans:2015 SEG Annual Meeting , 2015 :4127 -4131 .
[本文引用: 1]
[27]
田坤 , 张学涛 , 李国磊 . 添加正则化项的黏声逆时偏移成像方法研究
[J]. CT 理论与应用研究 , 2017 , 26 (6 ):669 -677 .
[本文引用: 1]
Tian K , Zhang X T , Li G L . Viscoacoustic reverse time migration by adding a regularization term
[J]. Computerized Tomography Theory and Applications , 2017 , 26 (6 ):669 -677 .
[本文引用: 1]
[28]
Kjartansson E . Constant Q-wave propagation and attenuation
[J]. Journal of Geophysical Research:Solid Earth , 1979 , 84 (B9 ):4737 -4748 .
[本文引用: 2]
[29]
Zhang Y , Zhang P , Zhang H . Compensating for visco-acoustic effects in reverse-time migration
[M]// SEG Technical Program Expanded Abstracts , 2010 :3160 -3164 .
[本文引用: 1]
[30]
Zhu T , Harris J M . Modeling acoustic wave propagation in heterogeneous attenuating media using decoupled fractional Laplacians
[J]. Geophysics , 2014 , 79 (3 ):S165 -S174.
[本文引用: 4]
[31]
吴玉 , 符力耘 , 陈高祥 . 基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移
[J]. 地球物理学报 , 2017 , 60 (4 ): 1527 -1537 .
[本文引用: 1]
Wu Y , Fu L Y , Chen G X . Forward modeling and reverse time migration of viscoacoustic media using decoupled fractional Laplacians
[J]. Chinese J. Geophys. , 2017 , 60 (4 ):1527 -1537 .
[本文引用: 1]
[32]
Zhu T , Harris J M . Improved seismic image by Q-compensated reverse time migration:Application to crosswell field data, west Texas
[J]. Geophysics , 2015 , 80 (2 ):B61 -B67.
DOI:10.1190/geo2014-0463.1
URL
[本文引用: 1]
[33]
罗文山 , 陈汉明 , 王成祥 , 等 . 时间域黏滞波动方程及其数值模拟新方法
[J]. 石油地球物理勘探 , 2016 , 51 (4 ):707 -713 .
[本文引用: 1]
Luo W S , Chen H M , Wang C X , et al . A novel time-domain viscoacoustic wave equation and its numerical simulation
[J]. OGP , 2016 , 51 (4 ):707 -713 .
[本文引用: 1]
[34]
Hu W , Zhou T , Ning J . An efficient Q-RTM algorithm based on local differentiation operators
[M]// SEG Technical Program Expanded Abstracts , 2016 :4183 -4187 .
[本文引用: 1]
[35]
Li Q , Zhou H , Zhang Q , et al . Efficient reverse time migration based on fractional Laplacian viscoacoustic wave equation
[J]. Geophysical Journal International , 2016 , 204 (1 ):488 -504 .
DOI:10.1093/gji/ggv456
URL
[本文引用: 1]
[36]
Sun J , Zhu T . Strategies for stable attenuation compensation in reverse-time migration
[J]. Geophysical Prospecting , 2018 , 66 (3 ):498 -511 .
DOI:10.1111/1365-2478.12579
URL
[本文引用: 1]
[38]
冀国强 , 石颖 . 正则化形式的稳定粘声逆时偏移成像方法
[J]. 石油物探 , 2020 , 59 (3 ):374 -381 .
[本文引用: 1]
Ji G Q , Shi Y . Stable and regularized visco-acoustic reverse time migration
[J]. Geophysical Prospecting for Petroleum , 2020 , 59 (3 ):374 -381 .
[本文引用: 1]
[39]
陈汉明 , 汪燚林 , 周辉 . 一阶速度—压力常分数阶黏滞声波方程及其数值模拟
[J]. 石油地球物理勘探 , 2020 , 55 (2 ): 302 -310 .
[本文引用: 1]
Chen H M , Wang Y L , Zhou H . A novel constant fractional-order Laplacians viscoacoustic wave equation and its numerical simulation method
[J]. OGP , 2020 , 55 (2 ): 302 -310 .
[本文引用: 1]
[40]
Liu Y , Sen M K . A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation
[J]. Geophysics , 2010 , 75 (2 ): A1 -A6.
DOI:10.1190/1.3295447
URL
[本文引用: 1]
中国石油 “十三五” 油气勘探重大成果与 “十四五” 发展战略
1
2021
... 目前,随着我国含油气盆地勘探程度的不断提高,勘探难度逐渐加大,复杂构造、复杂岩性、深层/超深层油气藏已逐渐成为主要勘探对象[1 ] ,勘探的地质目标具有“隐、碎、薄、小”的特点,如复杂小断块、致密碎屑岩等[2 -3 ] ,对地震偏移成像技术提出了更高的要求.井中地震具有资料信噪比高、储层识别能力强的优点.此外,它能够适应检波器的任意排列,可减少面元数据体丢失,理论上能够对井周的“隐、碎、薄、小”等储集体实现高精度成像. ...
中国石油 “十三五” 油气勘探重大成果与 “十四五” 发展战略
1
2021
... 目前,随着我国含油气盆地勘探程度的不断提高,勘探难度逐渐加大,复杂构造、复杂岩性、深层/超深层油气藏已逐渐成为主要勘探对象[1 ] ,勘探的地质目标具有“隐、碎、薄、小”的特点,如复杂小断块、致密碎屑岩等[2 -3 ] ,对地震偏移成像技术提出了更高的要求.井中地震具有资料信噪比高、储层识别能力强的优点.此外,它能够适应检波器的任意排列,可减少面元数据体丢失,理论上能够对井周的“隐、碎、薄、小”等储集体实现高精度成像. ...
中国沉积盆地油气立体综合勘探论
1
2021
... 目前,随着我国含油气盆地勘探程度的不断提高,勘探难度逐渐加大,复杂构造、复杂岩性、深层/超深层油气藏已逐渐成为主要勘探对象[1 ] ,勘探的地质目标具有“隐、碎、薄、小”的特点,如复杂小断块、致密碎屑岩等[2 -3 ] ,对地震偏移成像技术提出了更高的要求.井中地震具有资料信噪比高、储层识别能力强的优点.此外,它能够适应检波器的任意排列,可减少面元数据体丢失,理论上能够对井周的“隐、碎、薄、小”等储集体实现高精度成像. ...
中国沉积盆地油气立体综合勘探论
1
2021
... 目前,随着我国含油气盆地勘探程度的不断提高,勘探难度逐渐加大,复杂构造、复杂岩性、深层/超深层油气藏已逐渐成为主要勘探对象[1 ] ,勘探的地质目标具有“隐、碎、薄、小”的特点,如复杂小断块、致密碎屑岩等[2 -3 ] ,对地震偏移成像技术提出了更高的要求.井中地震具有资料信噪比高、储层识别能力强的优点.此外,它能够适应检波器的任意排列,可减少面元数据体丢失,理论上能够对井周的“隐、碎、薄、小”等储集体实现高精度成像. ...
中国石化物探技术新进展及发展方向思考
1
2021
... 目前,随着我国含油气盆地勘探程度的不断提高,勘探难度逐渐加大,复杂构造、复杂岩性、深层/超深层油气藏已逐渐成为主要勘探对象[1 ] ,勘探的地质目标具有“隐、碎、薄、小”的特点,如复杂小断块、致密碎屑岩等[2 -3 ] ,对地震偏移成像技术提出了更高的要求.井中地震具有资料信噪比高、储层识别能力强的优点.此外,它能够适应检波器的任意排列,可减少面元数据体丢失,理论上能够对井周的“隐、碎、薄、小”等储集体实现高精度成像. ...
中国石化物探技术新进展及发展方向思考
1
2021
... 目前,随着我国含油气盆地勘探程度的不断提高,勘探难度逐渐加大,复杂构造、复杂岩性、深层/超深层油气藏已逐渐成为主要勘探对象[1 ] ,勘探的地质目标具有“隐、碎、薄、小”的特点,如复杂小断块、致密碎屑岩等[2 -3 ] ,对地震偏移成像技术提出了更高的要求.井中地震具有资料信噪比高、储层识别能力强的优点.此外,它能够适应检波器的任意排列,可减少面元数据体丢失,理论上能够对井周的“隐、碎、薄、小”等储集体实现高精度成像. ...
Method of making sub-surface determinations
1
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
Downhole measurements of drill string forces and motions
1
1968
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
Linear signal processing and Ultrasonic transversal filters
1
1969
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
Walk-away VSP using drill noise as a source
1
1995
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
1
2007
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
1
2007
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
随钻VSP测量中地震波场的数值模拟
2
2017
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
... [9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ].在石油地震勘探领域,井中地震技术的研究相对较少. ...
随钻VSP测量中地震波场的数值模拟
2
2017
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
... [9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ].在石油地震勘探领域,井中地震技术的研究相对较少. ...
Wavefield processing of reverse VSP data
1
1991
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
变偏移距VSP射线追踪模型
1
1998
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
变偏移距VSP射线追踪模型
1
1998
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
1
2003
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
1
2003
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
2
2013
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
... 在偏移成像过程中由于井中地震特殊的观测方式,导致其成像范围有限,有效成像范围近似为“拱形”.图4 a为水平层状介质中,井中地震射线路径示意,从图中可以发现,对于井中地震数据,不同采样点对应的反射点位置在水平方向和深度方向均有变化,反射点轨迹成“拱形”.随深度增加,反射点位置逐渐远离震源,到达最深位置时,反射点与震源相距最大,水平距离为偏移距的一半.由此可见,当震源深度一定时偏移距越大,反射点水平方向覆盖范围越大;图4 b为不同激发深度覆盖范围分布,可以看到当震源处于不同深度时,其覆盖范围也有所差异.当检波器排列固定时,震源深度越大,其覆盖范围越小,最大覆盖范围为一半炮检距,此时震源位于井口[13 ] . ...
2
2013
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
... 在偏移成像过程中由于井中地震特殊的观测方式,导致其成像范围有限,有效成像范围近似为“拱形”.图4 a为水平层状介质中,井中地震射线路径示意,从图中可以发现,对于井中地震数据,不同采样点对应的反射点位置在水平方向和深度方向均有变化,反射点轨迹成“拱形”.随深度增加,反射点位置逐渐远离震源,到达最深位置时,反射点与震源相距最大,水平距离为偏移距的一半.由此可见,当震源深度一定时偏移距越大,反射点水平方向覆盖范围越大;图4 b为不同激发深度覆盖范围分布,可以看到当震源处于不同深度时,其覆盖范围也有所差异.当检波器排列固定时,震源深度越大,其覆盖范围越小,最大覆盖范围为一半炮检距,此时震源位于井口[13 ] . ...
RVSP等效地面处理方法研究
1
2015
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
RVSP等效地面处理方法研究
1
2015
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
1
2018
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
1
2018
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
Multi-hole joint acquisition of a 3D-RVSP in a karst area:Case study in the Wulunshan Coal Field,China
1
2020
... 井中地震技术按照震源的类型可分为随钻井中地震和爆震井中地震.随钻井中地震,即依靠随钻地震的思想,利用钻头的振动作为震源进行井地观测[4 ] ,基于该思想国内外很多学者对旋转钻头和随钻井中地震做了相关研究[5 ⇓ ⇓ ⇓ -9 ] .但随钻地震存在着探测深度有限的弊端,不利于深层复杂构造区域的地震勘探,由此发展了井中激发震源的爆震井中地震.对于爆震井中地震,由于其需要在保护井筒的同时激发出较强的能量,对震源的要求很高,因此相关研究较少且主要集中在浅层层析成像、煤层井中地震研究及工程物探领域[9 ⇓ ⇓ ⇓ ⇓ ⇓ ⇓ -16 ] .在石油地震勘探领域,井中地震技术的研究相对较少. ...
非常规油气概念, 特征, 潜力及技术——兼论非常规油气地质学
1
2013
... 在实际生产中,钻井的设计受到生产要求或客观因素影响,井型除常规的直井外还存在斜井、曲斜井、水平井等多种情况.与常规的直井相比水平井在非常规油气资源高效勘探和开发上有着重要的作用[17 ] .而在海洋勘探中,由于作业平台的空间有限,其井型一般都为有角度的斜井或空变的斜井[18 ] .由于井中地震的激发震源位于井中,不同井型使得震源所处空间位置有较大差异,造成地震波的传播路径存在显著不同,这势必会对最终的成像结果产生影响.因此在应用井中地震技术前,需详细分析井型对成像结果的影响,进而指导实际生产中激发井段的选择. ...
非常规油气概念, 特征, 潜力及技术——兼论非常规油气地质学
1
2013
... 在实际生产中,钻井的设计受到生产要求或客观因素影响,井型除常规的直井外还存在斜井、曲斜井、水平井等多种情况.与常规的直井相比水平井在非常规油气资源高效勘探和开发上有着重要的作用[17 ] .而在海洋勘探中,由于作业平台的空间有限,其井型一般都为有角度的斜井或空变的斜井[18 ] .由于井中地震的激发震源位于井中,不同井型使得震源所处空间位置有较大差异,造成地震波的传播路径存在显著不同,这势必会对最终的成像结果产生影响.因此在应用井中地震技术前,需详细分析井型对成像结果的影响,进而指导实际生产中激发井段的选择. ...
斜井 VSP 高斯射线束正演方法
1
2014
... 在实际生产中,钻井的设计受到生产要求或客观因素影响,井型除常规的直井外还存在斜井、曲斜井、水平井等多种情况.与常规的直井相比水平井在非常规油气资源高效勘探和开发上有着重要的作用[17 ] .而在海洋勘探中,由于作业平台的空间有限,其井型一般都为有角度的斜井或空变的斜井[18 ] .由于井中地震的激发震源位于井中,不同井型使得震源所处空间位置有较大差异,造成地震波的传播路径存在显著不同,这势必会对最终的成像结果产生影响.因此在应用井中地震技术前,需详细分析井型对成像结果的影响,进而指导实际生产中激发井段的选择. ...
斜井 VSP 高斯射线束正演方法
1
2014
... 在实际生产中,钻井的设计受到生产要求或客观因素影响,井型除常规的直井外还存在斜井、曲斜井、水平井等多种情况.与常规的直井相比水平井在非常规油气资源高效勘探和开发上有着重要的作用[17 ] .而在海洋勘探中,由于作业平台的空间有限,其井型一般都为有角度的斜井或空变的斜井[18 ] .由于井中地震的激发震源位于井中,不同井型使得震源所处空间位置有较大差异,造成地震波的传播路径存在显著不同,这势必会对最终的成像结果产生影响.因此在应用井中地震技术前,需详细分析井型对成像结果的影响,进而指导实际生产中激发井段的选择. ...
稳定的迭代法反Q滤波
1
2013
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
稳定的迭代法反Q滤波
1
2013
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
一种稳定高效的等效Q值反Q滤波算法及应用
1
2016
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
一种稳定高效的等效Q值反Q滤波算法及应用
1
2016
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Q-adaptive deconvolution
1
1982
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Similarity and the inverse Q filter:Some simple algorithms for inverse Q filtering
1
1992
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Viscoelastic true-amplitude prestack reverse-time depth migration
1
2008
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Attenuation compensation for least-squares reverse time migration using the viscoacoustic-wave equation
1
2014
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Attenuation compensation in viscoacoustic reverse time migration
1
2013
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Viscoacoustic reverse time migration by adding a regularization term
1
2015
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
添加正则化项的黏声逆时偏移成像方法研究
1
2017
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
添加正则化项的黏声逆时偏移成像方法研究
1
2017
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Constant Q-wave propagation and attenuation
2
1979
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
... 在地震勘探领域,介质的衰减特性通常被认为与频率呈线性关系.因此,在地震资料频带范围内,可以假定Q 与频率无关,即常Q 模型.基于此假设,Kjartansson[28 ] 推导了常Q 模型的速度相位频散和衰减关系: ...
Compensating for visco-acoustic effects in reverse-time migration
1
2010
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Modeling acoustic wave propagation in heterogeneous attenuating media using decoupled fractional Laplacians
4
2014
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
... 综上所述,井中地震是一种可以实现对井周“隐、碎、薄、小”地质目标进行精细成像的物探方法.在实际生产中,井型不同对成像结果存在影响,然而,当前井型对于偏移成像的影响具体如何,还没有定性或定量的认识.本文使用黏声逆时偏移方法进行井中地震成像,并基于Zhu等[30 ] 提出的分数阶黏声波方程进行数值模拟,通过在理论模型中进行直井、斜井、曲斜井、水平井与多井的成像试算,探究井型对偏移质量的影响.数值实验结果阐明了井型与成像剖面间的影响关系,为井中地震采集系统设计提供了理论支撑. ...
... 式中: ω 0 是参考角频率; ω 是角频率; c 0 是参考角频率对应的声波速度; c p 是声波速度; γ = t a n - 1 ( 1 / Q ) / π ,是无量纲参数; α 是衰减因子.基于上述常Q 频散关系,Zhu等[30 ] 推导了时间域分数阶黏声波方程: ...
... 当不考虑地层的黏滞性时波场正向延拓和逆时延拓在参数选取合理的情况下都是稳定的.但当地震波在黏滞性介质中传播时,由于波场正传过程是一个能量衰减过程,故该过程稳定.但在逆时延拓的过程中,黏滞波动方程会对震源和检波点波场进行补偿,使其能量呈指数型增大,与此同时波场中高频噪声也被呈指数放大,因此出现高频不稳定现象.为保证波场的稳定[30 ] ,通过在波数域设置低通滤波器将高频噪声去除.与此同时,本文考虑到在应用互相关成像条件进行逆时偏移成像时,成像过程中存在低波数噪声干扰现象.因此本文通过在成像剖面应用Laplace滤波的方式来压制低波数噪声. ...
基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移
1
2017
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移
1
2017
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Improved seismic image by Q-compensated reverse time migration:Application to crosswell field data, west Texas
1
2015
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
时间域黏滞波动方程及其数值模拟新方法
1
2016
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
时间域黏滞波动方程及其数值模拟新方法
1
2016
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
An efficient Q-RTM algorithm based on local differentiation operators
1
2016
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Efficient reverse time migration based on fractional Laplacian viscoacoustic wave equation
1
2016
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
Strategies for stable attenuation compensation in reverse-time migration
1
2018
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
A stable and efficient approach of Q reverse time migration
1
2018
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
正则化形式的稳定粘声逆时偏移成像方法
1
2020
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
正则化形式的稳定粘声逆时偏移成像方法
1
2020
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
一阶速度—压力常分数阶黏滞声波方程及其数值模拟
1
2020
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
一阶速度—压力常分数阶黏滞声波方程及其数值模拟
1
2020
... 早期的衰减补偿是对地震记录应用反Q 滤波达到波场放大的目的[19 -20 ] ,主要用于Kirchhoff偏移和单程波偏移中[21 -22 ] ,但仅是对振幅相位进行逐道补偿,没有考虑地震波传播路径问题,因此无法达到合理的补偿效果.而基于黏滞波动方程的Q 补偿偏移可以在波场延拓过程中进行补偿,更具有说服力.目前基于标准线性固体(SLS)模型的方程被广泛使用,用以描述衰减过程[23 -24 ] ,许多学者对其进行了相关研究并最终实现了稳定的Q 补偿逆时偏移(Q-RTM)[25 ⇓ -27 ] .总体来说基于SLS模型的方程,具有计算速度快、便于实现的优点,但同时也存在衰减项和频散项相互耦合的缺点.为解决这一弊端,诸多学者发展了以常Q 模型为基础的方程[28 -29 ] .其中,Zhu等[30 ] 、吴玉等[31 ] 基于平面波频散关系提出了解耦的分数阶黏声波方程,该方程显式分离了衰减项和频散项,同时对振幅和相位进行了校正,具有较高的精度,且在West Texas实际数据中[32 ] 进行了应用,并取得了良好的效果.但其同样存在需要利用伪谱法求解Laplace算子,计算量大,衰减和频散项解耦不彻底的问题.基于该方程国内外诸多学者进行了一系列研究[33 ⇓ ⇓ ⇓ ⇓ ⇓ -39 ] ,总的来说解耦的分数阶Laplace 方程具有较高的计算精度,便于实现稳定的Q -RTM. ...
A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation
1
2010
... 叠前逆时偏移分为3个步骤:①震源波场正传;②检波点波场反传;③应用成像条件成像.在波场模拟的过程中本文选用Liu等[40 ] 的混合吸收边界作为边界条件,该方法将计算区域分为内部区、过渡区和边界区.在进行波场模拟时,内部区中用双程波方程求解,边界区用单程波方程求解,过渡区波场由单程波和双程波的波场加权得到;在成像条件的选取上,本文选用了在成像过程中同时利用了震源波场和检波点波场的信息,可以对所有波进行成像的互相关成像条件.互相关成像条件可表示为: ...