面炮叠前深度偏移与控制照明技术

基于波场延拓的叠前深度偏移是实现复杂构造地质体成像的最可靠方法,但存在着计算量大、对观测系统适应性差等缺点。面炮偏移是波动方程实现精确叠前成像的另一类方法,具有较高的计算效率,不存在偏移孔径问题,而且可以通过控制照明方法,解决平面波在目标区域的能量补偿问题。本文采用面炮成像技术进行叠前深度偏移成像,通过对面炮震源下行波场的质量控制以及射线参数的个数和范围的选取,以达到最佳的成像效果。采用不同深度点上的控制照明技术,较大地提高了目标地层的成像精度。数据实验表明面炮成像技术是一种快速有效的方法,其成像精度与单平方根算子的共炮点道集偏移和双平方根算子的共中心点道集偏移相当,但在计算速度上要快得多,而且易于并行计算。     

叠前地震数据特征波场分解、偏移成像与层析反演

本文提出了一套叠前地震数据稀疏表达(特征波场合成)、深度偏移成像和层析成像的处理流程.不同于传统的变换域中的数据稀疏表达理论, 本文利用局部平面波的传播方向(慢度矢量), 在中心炮检点处同时进行波束合成, 从而将地震数据投影到局部平面波域(高维空间)中.由于波束合成后的地震数据描述了局部平面波的方向特征, 因此称之为特征波场.然而波束合成算法需要估计局部平面波的慢度矢量.当地震数据受噪声干扰时, 难以在常规τ-p谱中自动估计局部平面波的射线参数(慢度矢量).本文提出了基于反演理论的特征波场合成方法, 可以同时反演局部平面波及其传播方向, 从而提高特征波合成的自动化程度并保持方法的稳健性.通过特征波场合成, 可以将地震数据分解为单独的震相(波形).这样的数据可以直接用来成像及反演.在局部平面波域中, 由于局部平面波的入射与出射射线参数已知, 传统的Kirchhoff叠前深度偏移(PSDM)和高斯束/控制束PSDM可以实现从"沿等时面的画弧"到"向反射点(段)的直接投影"的转变, 叠前偏移的效率以及成像质量可以同时提高.此外, 特征波场与地下反射点(段)的一对一映射关系使得叠前深度偏移与层析成像融为一体, 可以极大地提高速度反演的效率.数值试验证明了特征波场合成、叠前深度成像以及层析反演的有效性.     

表驱三维角度域Kirchhoff叠前时间偏移成像方法

针对传统Kirchhoff叠前时间偏移方法的一些不足,以及振幅随入射角、方位角变化(AVA/AVAZ)分析的需要,本文提出一种基于射线理论的三维叠前时间偏移角度域成像方法.它通过横向均匀介质中稳健的射线追踪建立单程波走时和传播角度的数值表,然后在此基础上估算反射波双程走时以及在界面处传播的方位角和入射角,最后基于脉冲响应叠加原理获得三维构造图像和方位\|角度域共成像点道集.与传统方法不同之处在于,上述过程均考虑了地震波在垂向变速介质中的射线弯曲效应和三维传播特征,有利于准确提取随入射角和方位角变化的振幅和时差信息.理论模型合成数据和实际地震资料测试结果展示了方法的优越性与实用性.     

基于共聚焦点道集的叠前深度偏移

共聚焦点(CFP)偏移技术是一种基于等时原理,将Kirchhoff积分法的一步偏移分两步聚焦（即激发聚焦和检波聚焦）来完成的叠前地震成像方法.该方法借助于逆时聚集算子和共聚焦点道集来实现叠前偏移成像.基于共聚焦点道集的叠前深度偏移是把基于共炮集的深度偏移的算法引入到CFP技术上来,基于波场延拓的理论来实现偏移成像,该方法首先生成共聚焦点道集,然后基于面炮合成的理论合成聚焦震源,最终通过相关成像来实现叠前偏移成像.该方法选取较少的聚焦点就可以实现对于地下构造的偏移成像,和炮域波动方程偏移相比,其计算效率得到了提高.通过模型试算和实际资料的试处理,验证了该方法在实现叠前深度偏移成像上的正确性和有效性.     

基于Rayleigh积分的聚焦束估算层状介质的三维地震观测系统分辨率

通过对聚焦束计算过程的简化,用Rayleigh积分实现了三维层状介质的三维观测系统聚焦束的计算.(1)用目标点处激发的地表响应代替CFP道集的形成;(2)用单频谐波计算代替全频率范围的计算;(3)用插入虚拟界面的方法实现了起伏界面的Rayleigh积分波场外推.此三点改进实现了在微机上进行三维层状介质的聚焦束计算,从而估算出三维观测系统的分辨率.此项研究对于三维观测系统的优化设计具有指导作用.     

复杂构造三维地震观测系统设计的共聚焦分辨率分析

三维地震观测系统共聚焦分辨率特性分析突破传统以点论证为基础的观测系统分辨率分析方法,面向地质目标定量预测三维观测系统地震成像的空间分辨率和振幅精度.基于Fourier有限差分（FFD）大步长波场延拓和Born-Kirchhoff小步长波场插值递推方法,本文介绍了一种复杂介质条件下三维地震观测系统共聚焦分辨率特性快速分析方法.对给定的速度模型,该方法能够分析拟采用的三维地震观测系统设计方案对复杂构造的成像分辨率与AVP属性,从而为进一步的偏移成像与储层分析提供保证.最后本文以SEG/EAGE三维盐丘模型为例设计满覆盖为16次的三维地震观测系统,并实施三维共聚焦分辨率特性分析.     

基于Gabor框架合成小波束叠前深度偏移

利用小波束域波场分解和传播在空间及方向上的双重局域性,通过Fourier传播算子进行深度外推波场,完成小波束的传播,最终应用成像条件求取成像值.合成小波束和小束源记录,进行部分小束源叠前偏移成像的数值试验,获得的成像结果表明,小波束叠前深度成像是切实有效的.     

基于拟线性Born近似的叠前深度偏移方法

在频率波数域和频率空间域实现了一种基于拟线性Born近似的叠前深度偏移方法,并在2-D空间进行了Marmousi模型炮集数据的处理.通过与Split-StepFourier、Phase-Screen和稳定的Born近似叠前深度偏移等方法比较,认为基于拟线性Born近似的叠前深度偏移方法不仅在效果上要优于前三者,而且还能更好地处理速度横向变化.在散射波场计算中,使用了一个更稳定的散射波场计算公式,扩大了拟线性Born近似的应用范围,使基于拟线性Born近似的叠前深度偏移方法能够适应更强的横向速度变化.     

基于拟线性Born近似的叠前深度偏移方法

在频率波数域和频率空间域实现了一种基于拟线性Born近似的叠前深度偏移方法，并在2-D空间进行了Marmousi模型炮集数据的处理.通过与Split-StepFourier、Phase-Screen和稳定的Born近似叠前深度偏移等方法比较，认为基于拟线性Born近似的叠前深度偏移方法不仅在效果上要优于前三者，而且还能更好地处理速度横向变化.在散射波场计算中，使用了一个更稳定的散射波场计算公式，扩大了拟线性Born近似的应用范围，使基于拟线性Born近似的叠前深度偏移方法能够适应更强的横向速度变化.     

双复杂条件下的波动方程叠前深度偏移

复杂地表条件下的地震勘探越来越被人们所关注.双复杂条件下的叠前深度偏移方法是解决复杂地表条件和复杂地质构造成像的有效手段."波场上延"法能实现由非水平观测界面开始的偏移过程,解决复杂地表对地下构造成像的影响.复杂理论模型的试算以及实际资料处理表明,"波场上延"方法较好地克服了起伏地形对地下构造成像的影响,取得了令人满意的效果,实现了波动方程基准面校正和深度成像的有机结合.     

Diffraction problems of 3D seismic imaging

Migration is essential to seismic imaging. It is carried out by backward extrapolation of the wavefield registered on the observation surface. The quality of images depends on the accuracy of the wavefield reconstruction at interior subsurface points. From the theory based on the exact solution of the scalar wave equation it is known that, for accurate wave extrapolation, data must be obtained from an infinite observation surface. Limiting of migration apertures, which is inevitable in practice, leads to artefacts in extrapolated fields. The distortion they cause in 2D and 3D imaging is different. In 2D migration, the artefacts known as truncation effects are much weaker than the signals being extrapolated and for this reason attract no special attention. In 3D migration, diffractions caused by an aperture edge are stronger and may create serious problems. For a circular aperture, their amplitudes are comparable to the amplitudes of the signals themselves. The study of aperture diffractions is intended to help in the search for ways of either suppressing them efficiently or deliberately utilizing them in order to improve imaging.
In optics, diffractions by an aperture play a constructive role in image making. This research shows that the same may take place in seismic imaging.     

控制照明与面向目标成像的观测系统设计

常规地震观测系统设计方法基于地下水平层 状介质的假设,通常不能适应复杂构造情况。我们 从控制照明的思想出发,提出了一种面向目标成像 的地震观测系统设计方法,该方法需要一个由初步 地震解释得到的速度模型。利用单程傅立叶有限差 分波场传播算子将目标层的平面源延拓到地表,通 过分析从目标层延拓到地表的波场能量的分布,可 以确定目标层成像所需要的炮点或者检波点的位置。 利用SEG-EAGE盐丘模型数值试算结果表明,该方法 用于设计面向目标成像的特定地震采集系统。     

基于BS-MUSIC的定向钻孔雷达三维成像算法

定向钻孔雷达由位于井中的一个发射天线和分布于圆环上的呈等角度分布的四个接收天线组成,能确定井旁地质目标的深度和方位.然而,目标方位确定和三维成像的算法问题一直是该技术面临的挑战.本文将空间谱估计中的BS-MUSIC算法应用于定向钻孔雷达目标方位识别算法中,在此基础上提出并实现了利用四个接收天线的信号对井周地质体进行三维空间成像的算法.使用合成数据对以上算法的效果进行验证,并和通常的MUSIC算法进行比较,结果表明,MUSIC算法和基于波束空间的MUSIC算法均可以在方位确定和三维成像方面取得非常好的效果,但计算速度方面,BS-MUSIC算法相比于MUSIC算法有了非常明显的提升.

     

Semi-automatic fault/fracture interpretation based on seismic geometry analysis

Fault and fracture interpretation is a fundamental but essential tool for subsurface structure mapping and modelling from 3D seismic data. The existing methods for semi-automatic/automatic fault picking are primarily based on seismic discontinuity analysis that evaluates the lateral changes in seismic waveform and/or amplitude, which is limited by its low resolution on subtle faults and fractures without apparent vertical displacements in seismic images. This study presents an innovative workflow for computer-aided fault/fracture interpretation based on seismic geometry analysis. First, the seismic curvature and flexure attributes are estimated for highlighting both the major faults and the subtle fractures in a seismic volume. Then, fault probability is estimated from the curvature and flexure volumes for differentiation between the potential faults and non-faulting features in the geometric attributes. Finally, the seeded fault picking is implemented for interpreting the target faults and fractures guided by the knowledge of interpreters to avoid misinterpretation and artefacts in the presence of faulting complexities as well as coherent seismic noises. Applications to two 3D seismic volumes from the Netherlands North Sea and the offshore New Zealand demonstrate the added values of the proposed method in imaging and picking the subtle faults and fractures that are often overlooked in the conventional seismic discontinuity analysis and the following fault-interpretation procedures.     

三维地震与地面微地震联合校正方法

由于地面微地震监测台站布设在地表,会受到地表起伏、低降速带厚度和速度变化的影响,降低了微地震事件的识别准确度和定位精度,限制了地面微地震监测技术在复杂地表地区的应用.因此,将三维地震勘探技术的思路引入到地面微地震监测中,提出了三维地震与地面微地震联合校正方法,将油气勘探和开发技术更加紧密地结合在一起.根据三维地震数据和低降速带测量数据,通过约束层析反演方法建立精确的近地表速度模型,将地面微地震台站从起伏地表校正到高速层中的平滑基准面上,有效消除复杂近地表的影响.其次,根据射孔数据和声波测井速度信息,通过非线性反演方法建立最优速度模型,由于已经消除复杂近地表的影响,在进行速度模型优化时不需要考虑近地表的影响,因而建立的速度模型更加准确.最后,在精确速度模型的基础上,通过互相关方法求取剩余静校正量,进一步消除了复杂近地表和速度模型近似误差的影响.三维地震与地面微地震联合校正方法采用逐步校正的思路,能够有效消除复杂近地表的影响,提高微地震数据的品质和速度模型的精确度,保证了微地震事件的定位精度,具有良好的应用前景.

     

Seismic hazard of the Northern Apennines based on 3D seismic sources

Seismic hazard has been computed for the Northern Apennines in northern Italy based on a new seismogenic zonation. This zonation considers inclined (dipping) planes as seismogenic sources, defined on the basis of all the seismotectonic information available so far. Although these geometries are extremely rough because they simplify with a few inclined elements the totality of faults constituting a source, this model mimics the tectonic style better than that based on horizontal planes. Nevertheless, for a comparison between the new ground motions obtained and those available in the literature, the plane version of the zonation has been developed, where horizontal areas (the standard seismogenic zones), representing the surficial projection of the inclined planes, are used as seismogenic sources.     

基于稀疏约束反演的三维混采数据分离

混合震源采集技术相对于传统的地震数据采集,在极大提高采集效率的同时引入了混叠噪声,很大程度上影响了成像结果的精度.二维混采数据中,我们通常利用混叠噪声在非共炮域呈非相干分布这一特点来压制混叠噪声,从而实现混合震源数据分离.相对于二维混采数据,三维混采数据具有数据量巨大,构建混合震源算子困难,混合度的增加引入了高强度混叠噪声的特点.针对上述问题,本文采用稀疏约束反演方法在Radon域实现混采数据分离,混叠噪声强度比较大的情况下,稀疏约束反演方法能够得到更高精度的分离结果;利用震源激发的GPS时间通过长记录的方式在共接收点道集对上一次迭代分离结果做混合、伪分离,实现了单个共接收点道集自身混合、伪分离,避免了对整个数据做运算,同时不需要构建混合震源算子.通过模拟数据和实际数据计算来验证上述方法的适用性.

     

2D多尺度非线性地震速度成像

将遗传算法和单纯形算法相结合,得到了一种高效、健全的2D混合地震走时反演方法.把速度场划分为不同的空间尺度,定义网格节点上的速度作为待反演参数,采用双三次样条函数模型参数化,正问题采用有限差分走时计算方法,反问题采用多尺度混合反演方法.首先在较大的空间尺度内反演,然后减小空间尺度,将大尺度的反演结果作为次一级尺度反问题的初始模型,再进行混合反演,如此类推逐次逼近全局最优解.一个低速度异常体的数值模拟试验和抗走时扰动试验表明该方法是有效和健全的.我们将该方法应用到青藏高原东北缘阿尼玛卿缝合带东段上部地壳速度结构研究中,并与前人的成果进行了对比.     

三维共偏移距叠前地震数据反假频射线束形成偏移成像

对稀疏/非规则采样或者低信噪比数据,射线束提取困难并伴随有假频产生,对叠加剖面和道集造成严重干扰.为了提升射线束偏移在稀疏和低信噪比地震数据采集中的成像效果,本文提出基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法.射线束偏移首先将地震数据划分为超道集,经过部分NMO后转化为以射线束中心定义的共偏移距数据,倾斜叠加和反假频操作均在局部共中心点坐标上实现.时间域倾斜叠加是对地震数据的时移累加操作,三角低通滤波同样可以在时间域完成,在对地震数据进行因果和反因果积分后,亦为地震数据的时移累加.因此,三角低通滤波与倾斜叠加可在时间域结合同时完成,避免了频域滤波的正反傅里叶变换.本文在反假频公式中加入权重系数,用以对反假频的程度进行控制,达到分辨率和噪声压制的最佳折衷.以某海上三维实际数据为例,文中展示了反假频射线束形成对偏移叠加剖面和共成像点偏移距道集中的噪声进行了有效压制.

     

基于DMVPN接入地震行业网络设计

介绍山东省地震行业网络建设中县级地震台网使用的动态多点VPN接入方式,对中心站点和分支站点的配置进行说明,分析DMVPN相较于其他VPN的优势。作为一种经济有效的地震行业网络接入方式,DMVPN具有较强的推广、应用价值。