VTI介质隐式有限差分平面波偏移

这里研究了VTI（Vertical Transverse Isotropy）介质隐式有限差分（IFD）波场外推算子和VTI介质平面波偏移。在文中设计了VTI介质IFD波场外推算子,将Taylor分析法求得的差分系数,作为初始解,用非线性优化方法迭代求得算子系数。将设计的外推算子频散曲线和偏移脉冲响应与理论解进行了对比分析,验证出外推算子有较高的精度。从各向同性介质平面波偏移理论出发,结合VTI介质IFD波场外推算子,将平面波偏移推广到VTI介质中。对Hess VTI标准模型拟合数据进行偏移,偏移结果验证了平面波偏移在VTI介质中的有效性,以及在加强特定陡倾界面成像方面的优势。并且,与常规VTI介质叠前深度偏移相比,VTI介质平面波偏移在保证成像质量的前提下,大大减少了计算量,为各向异性偏移成像研究提供了新的技术思路。     

井地联合提取VTI介质各向异性参数

在新疆塔里木盆地某地区的宽方位三维地震数据中,VTI介质火成岩的各向异性和垂向速度变化引起非双曲线时差。通过对实际资料的系统分析,我们不能利用DMO反演法,移动震源VSP反演法和多偏移距VSP初至旅行时反演法提取Thomsen参数用于各向异性地震资料的处理。因此本文提出一个新的方法:运用地震资料和零偏移距VSP井求取火成岩的各向异性参数δ;在比较精确确定火成岩地层深度和厚度的情况下,利用偏移距VSP井,求取了各向异性参数ε,从而获得火成岩的瞬时各向异性参数η;充分考虑VTI介质垂向速度变化对时差曲线的影响,获得了各向异性参数ηeff。最后,利用所提取的各向异性参数,进行VTI介质各向异性速度分析和成像,成像质量得到明显改善。     

共偏移距可分表示法叠前深度偏移

基于双平方根方程的共偏移距可分表示法叠前深度偏移用于复杂介质成像, 该方法在中点-偏移距坐标中同时向下延拓炮点和检波点波场, 实现方式采用正反傅立叶变换.构造的双平方根方程波场延拓算子能够使波数域变量与空间(速度) 域变量分离, 波数域内进行相移计算, 在空间域对因介质横向变速引起的时移作修正.地震数据偏移不需要逐炮计算, 具有较高的效率.对Marmousi模型数据的偏移成像结果显示, 该方法较好地成像强横向变速介质中的复杂构造.      

基于双相介质弹性波波场分离的逆时偏移方法

随着地震勘探的发展,以单相介质模型为基础的声波逆时偏移和弹性波逆时偏移很难比较精确地描述地下介质和满足复杂的工程要求,尤其是在处理地下含流体多孔介质时,所以需要引入双相介质模型。首先基于Biot理论双相介质一阶速度-应力弹性波方程,利用散度场和旋度场的波场分离理论将固相和流相中的纵横波场解耦分离,并由波场快照图得出固相与流相中的纵波波场与横波波场均完成分离,而快纵波与慢纵波相互耦合均存在于纵波波场中。采用高阶交错网格有限差分方法,结合PML吸收边界条件,推导了具有波场分离的双相介质弹性波正传与反传的偏移延拓算子,结合互相关成像条件,可以得到双相介质弹性波固相与流相中的各分量偏移成像剖面。最后对模型使用该双相介质弹性波波场分离的逆时偏移方法进行偏移成像试算并取得成功,效果较常规双相介质弹性波逆时偏移成像效果有了明显提升。     

基于互相关成像条件的VTI介质保幅高斯束逆时偏移

保幅高斯束逆时偏移是一种兼具高斯束偏移灵活、高效和逆时偏移高精度特性的深度域成像算法,不仅能够消除偏移过程中由于入射角变化和几何扩散引起的振幅误差,而且具有面向目标成像的能力。其以Kirchhoff积分为基础,采用高斯束的加权积分来构建格林函数,进而实现正反向波场的延拓,最后通过相应的成像条件进行成像。这里从基于相速度的各向异性射线追踪算法出发,采用互相关成像条件,将保幅高斯束逆时偏移推广到VTI介质。随后通过模型试算对方法的正确性和有效性进行了验证。     

基于解耦传播的波场分解方法在VTI介质弹性波逆时偏移中的应用

由于弹性波逆时偏移更符合实际情况,而且转换波的成像结果有更高的分辨率,因此弹性波逆时偏移越来越受到重视。然而,弹性波逆时偏移需要多波多分量数据,为了减少成像结果中的串扰假象,在逆时偏移过程中进行P波和S波分解就变得非常必要。结合基于向量的激发振幅成像条件,我们把基于解耦传播的波场分解方法应用到弹性波逆时偏移中,并对比了其在各向同性介质和横向各向同性（VTI）介质中的应用效果。结果说明,该方法可以在各向同性介质中完全分解P波和S波,并保留波场的向量信息。虽然在VTI介质中有较小的分解残余,但是该分解残余不会在逆时偏移结果中产生明显的串扰;因此,这种波场分解方法可以应用于各向同性介质和VTI介质弹性波逆时偏移。该方法是在时间空间域实现的,可以在波场传播过程中直接对P波和S波进行分离,应用方便,计算效率高。与利用Helmholtz分解的弹性波逆时偏移相比,该方法避免了在PS波成像结果中的极性反转问题。在复杂Hess VTI模型的逆时偏移结果中,高速岩体和断层的成像清晰,甚至是两个低速薄夹层也能较好成像;这说明该方法对复杂介质具有较好的适应性。PS波成像结果中的各向异性体成像清晰,说明各向异性介质弹性波逆时偏移可以对传统逆时偏移不能很好成像的构造进行成像。     

煤矿采区三维地震叠前深度偏移技术研究

影响波动方程叠前深度偏移成像精度的主要因素除了要求勘探区覆盖次数均匀、单炮记录信噪比高外,还必须建立准确的偏移速度场及选择合适的偏移孔径,其中速度深度模型的建立原则取决于偏移后共反射点道集同相轴是否呈水平形态,而偏移孔径的选择依据则是目的层倾角的大小。在地层倾角变化剧烈的吉新煤矿区,深度域层速度场采用三次速度迭代求得,并确定偏移孔径为600m×600m。应用表明：叠前深度偏移比叠后时间偏移处理在复杂构造区成像准确、断点空间位置清晰,且在解决复杂地质构造成像问题的同时还能够提高资料的信噪比和分辨率,是解决复杂构造成像的有效工具。     

基于P＋S波震源的弹性波叠前逆时偏移方法

基于弹性波动方程,推导了二维非均匀介质条件下纵横波解耦且无转换波的一阶速度一应力弹性波方程,使得弹性波在延拓过程中实现矢量偏移保证偏移质量,且不产生其他类型的转换干扰波。从激发时间成像条件以及反射系数理论出发,并与纵横波解耦且无转换波的一阶速度一应力弹性波方程结合,给出了适用于P＋s波震源模拟数据的成像条件以及偏移方法。文中给出了两个典型模型的理论记录的偏移结果,试验表明,笔者提出的方法能够对P＋s波震源的产生的波场准确归位,提高资料处理精度。     

多信息多技术深度速度建模技术在复杂构造中的应用研究

中国西部等复杂区的油气勘探潜力较大,但地质条件复杂,造成常规偏移成像的精度不够,给解释带来了诸多困难。根据其地质特点,围绕高精度成像的需求,开展高精度速度建模技术研究,形成了配套的技术策略,为其进一步改善复杂区的成像精度,落实勘探目标奠定了基础。在该配套技术中:基于真地表的速度建模技术可以通过真地表数据准备来解决现有大圆滑偏移基准面对成像波场的改变、初至层析近地表建模技术的应用弥补了反射建模技术在近地表应用的不足、多信息约束速度建模可提高速度模型的可靠性、各向异性假设的成像技术更加接近复杂构造地下介质的假设,这些技术的配套应用可以有效提高复杂区成像位置的精度和准确度。在实际资料中的应用充分说明该技术的有效性。     

基于程函方程的克希霍夫全平面域偏移成像方法

本文分别以直射线和弯曲射线初至波层析反演技术重建的速度为偏移速度,采用波前扩展外推方式求解程函方程,求取了成像区域内各节点初至波时间场,应用克希霍夫积分偏移技术对复杂模型叠前数据进行全平面偏移成像,成像结果表明:弯曲射线层析反演方法获得的速度场较直射线层析结果对速度异常体的刻画更准确,以弯曲射线层析反演速度场作为偏移速度可以获得较好的偏移成像结果。     

利用高阶交错网格有限差分法数值模拟VTI介质井孔声场

横向各向同性（TI）介质是岩石地球物理中常见的一种现象,研究其井孔声场传播特征对声波测井理论以及为声波测井解释提供依据具有重要意义。针对具有垂直对称轴的横向各向同性（VTI）介质,根据柱坐标系条件下的弹性波波动方程,推导了速度-应力交错有限差分公式,采用时间二阶、空间十阶的交错有限差分算法对VTI介质中的井孔声场进行数值模拟。给出了在均匀介质中井孔声场不同时刻的波场快照,以及不同各向异性系数的VTI介质中的波场快照,计算了井轴上声源激发出的声波全波列波形。结果表明,在其他条件不变的条件下,VTI地层的各向异性系数的增大对横波的传播影响不大,但会使得纵波在纵向上的传播速度相对变小,径向上变化不大。各向异性系数的增大会使声波测井全波列首波信号时差变大,声波幅度略变小。      

纵波VTI介质各向异性参数的求取

在新疆塔里木盆地某地区的宽方位三维地震数据中,通过速度分析和不同炮检距范围数据的叠加剖面的比较,发现了各向同性基于双曲线时差曲线的NMO动校正时,远炮检距数据存在比较大的剩余动校正量。分析表明,是由于VTI介质火成岩的各向异性和垂向速度变化引起非双曲线时差。应用纵波短排列地震数据求取的均方根动校正速度,分别以扫描法和旅行时法得到了各自相互验证的各向异性参数。最后,利用所提取的各向异性参数,进行VTI介质各向异性速度分析、动校正和成像,剩余动校正得到消除,成像质量得到明显改善。     

VTI介质多波速度与各向异性系数求取及应用

针对具有垂直对称轴的横向各向同性(VTI)介质,研究了多波速度与地层各向异性系数之间的关系。采用适应大偏移距以及强各向异性介质的双平方根方程,进行速度和各向异性系数分析,并将批量计算与交互解释相结合,实现了高精度各向异性介质多波速度与各向异性系数提取,并形成了一套利用多波地震资料求取地震波纵、横波速度及各向异性系数的方法和软件。所求取的参数可用于多波地震资料处理,提高多波资料的成像质量,也可为地层地质解释提供参考。实际资料的应用效果证明了该方法和软件的正确性及有效性。     

基于P-SV波分离的VTI介质射线追踪方法

近年来,VTI介质参数建模和成像技术发展迅速,其中VTI介质射线追踪是基础,其精度和稳定性直接影响建模和成像的效果。VTI介质射线追踪大多通过前人提出的频散关系推导出程函方程和射线方程,该方法主要基于纵横波相互独立的假设,采用令横波速度等于零的手段,在某些特殊情况下精度较低、稳定性较差。针对其缺点提出一种新的射线追踪方法,该方法基于P-SV波分离得到的q P波频散关系,从而推导出程函方程和射线方程,提高了精度和稳定性。最后通过模型试算验证该方法的正确性和稳定性。     

弱各向异性介质地震波传播特征分析

随着油气田勘探精度要求越来越高,各向异性已成为不得不考虑的重要因素,因此对地震波在各向异性介质中传播特征的研究具有一定的价值。利用Thomsen提出的弱各向异性理论,通过理论模型数值计算,分析了Thomsen各向异性参数对弱各向异性VTI介质中三种体波的影响。分析结果表明:①参数ε、δ分别影响qP波大角度、中小角度传播的各向异性;②参数(ε-δ)影响qSV波中等角度传播时的各向异性;③参数γ影响qSH波大角度传播的各向异性,且群速度的各向异性要比相速度的各向异性更明显,并会出现复杂的交会区。通过这些研究分析,进一步认识了弱各向异性介质相速度及群速度随各向异性参数变化的规律。     

Angle-domain common-image gathers from anisotropic Gaussian beam migration and its application to anisotropy-induced imaging errors analysis

An approach for extracting angle-domain common-image gathers (ADCIGs) from anisotropic Gaussian beam prestack depth migration (GB-PSDM) is presented in this paper. The propagation angle is calculated in the process of migration using the real-value traveltime information of Gaussian beam. Based on the above, we further investigate the effects of anisotropy on GB-PSDM, where the corresponding ADCIGs are extracted to assess the quality of migration images. The test results of the VTI syncline model and the TTI thrust sheet model show that anisotropic parameters ε, δ, and tilt angle ??, have a great influence on the accuracy of the migrated image in anisotropic media, and ignoring any one of them will cause obvious imaging errors. The anisotropic GB-PSDM with the true anisotropic parameters can obtain more accurate seismic images of subsurface structures in anisotropic media.     

各向异性上覆层对裂缝储层地震AVAz影响的模拟及分析

为认识和研究上覆层为VTI时,垂直裂缝发育储层的地震反射AVAz响应特性,提出了相适应的正演模拟方法。该方法类比Shoenberg传播矩阵法中水平慢度分量的表达形式,引入了用入射角、方位角和VTI介质中相速度来表达的水平慢度分量作为输入,简化了各向异性反射率法的框架。利用该方法模拟和分析了各向异性层状模型的地震反射响应随方位的变化,讨论了入射角、上覆层各向异性程度对目标层顶底界面反射波AVAz特征的影响。结果表明,采用适当大的观测入射角进行AVAz分析更为合理,上覆层为VTI时会减弱目标层的AVAz特征。测井尺度理论模型的数值模拟说明了该方法的实用性。     

倾斜层位偏移规律研究

层状介质模型是对实际地质情况较为准确的近似。而在实际的数据处理过程中,以等效均匀介质代替层状介质,波传播的几何路径、介质速度的差别会使对倾斜界面位置的判断偏离实际情况。针对设计的理论模型和从测井数据提取的实际介质模型,笔者运用几何射线理论和波场数值模拟技术,定量计算并提取了层状介质模型和等效均匀介质模型两种情况下反射点在时间和空间上的偏移量。结果证明,等效均匀介质模型近似会造成对倾斜界面位置的判断不准确; 准确的速度模型是偏移归位研究的关键,应尽量使用由地震测井资料求得的瞬时速度。