局部指数标架小波束在复杂介质方向照明分析中的应用

提出了用局部指数标架小波束进行角度域分解的方法,解决了局部余弦基和局部正弦基缺乏单一方向性的问题.局部指数标架由局部余弦基和局部正弦基线性组合而成,是冗余度为2的紧标架.利用局部余弦变换和局部正弦变换的快速算法,能使基于局部指数标架进行方向照明分析的计算效率较常用的局部倾斜叠加和Gabor-Daubechies标架等方法具有更为明显的优势.通过计算二维SEG/EAGE模型和SIGSBEE模型的方向照明图以及采集系统倾角响应图证实了本文方法的有效性.该方法的高效性使其在三维模型的方向照明分析和大规模的工业应用中具有广阔前景.     

一种起伏地表条件下的照明补偿方法

起伏的地表条件限制了采集孔径范围并造成深层地震照明不足,为改善该类地区的成像质量,本文提出了一种起伏地表条件下的照明补偿方法.首先,基于小波束波场延拓算子和逐步累加的外推方法在波场延拓过程中解决起伏地表面的影响,并引入空间滤波函数压制虚拟层内的偏移噪音;其次,利用局部指数标架对上、下行波场分解,得到局部角度域成像和照明补偿因子.再次,利用计算出的成像值和照明补偿因子,在局部倾角域完成照明补偿.SEG起伏地表模型测试证明了本方法的有效性,深层构造照明度明显加强,不同角度成像振幅更加均衡,该技术为提高起伏地表地区的成像品质提供了新的手段.     

基于局部余弦基小波束的观测系统沉降法叠前深度偏移

将局部余弦基小波束波场分解、传播与观测系统沉降法叠前深度偏移相结合,推导了源-检波器观测系统沉降法传播算子.本算法中,先对频率域的共点源和共点检波器道集做局部余弦小波束分解,然后分别沿共小波束源和共小波束检波器在深度方向延拓得到下一层波场.每个深度的波场,都等效于把源和检波器放在该层后所能接收到的地震记录,每点的像值由炮点和检波点重合时的零时刻波场值给出.通过二维SEG/EAGE盐丘模型和Marmousi模型的偏移成像结果验证该方法理论推导的正确性.另外,结果显示该方法继承了小波束域波场延拓在速度扰动较大情况下波传播及成像精度高的优点.     

散射系数矩阵的计算及其应用

利用GaborDaubechies（G?D）小波束域波场分解和传播在空间和方向上的双重局域性，提出了基于G?D小波束域叠前深度偏移进行角度域成像和计算局部散射系数矩阵的方法. 以简单分层模型为例，对不同探测系统的局部散射系数矩阵分布特征进行分析. 分析结果表明，在一定的探测系统几何布局下，由本文方法得到的局部散射系数矩阵能够较真实地反映局部结构的散射(或反射)特性. 通过局部散射系数矩阵进一步外推具有不同速度反差的水平界面随角度变化的反射系数，并估计界面的空间位置和倾角等说明局部散射系数矩阵的潜在应用.     

基于Gabor-Daubechies小波束叠前深度偏移的角度域共成像道集

传统炮检距域共像集（CIG）在复杂介质中因波传播的多路径而存在反射体位置不确定的问题. 角度域CIG由于克服了这一缺陷而逐步成为速度分析、AVA以及振幅保真偏移成像等研究的主要手段. 以波动理论为基础的地震偏移成像方法的发展为获得高质量的角度域CIG提供了可靠的实现途径. 其中,基于波场局域化分解和传播的小波束域波场延拓和偏移成像方法,因其波场分解基本函数和传播算子在空间和方向上的双重局域特性,而成为角度相关分析研究的有效工具. 本文在采用Gabor Daubechies框架分解的小波束叠前角度域偏移成像基础上,利用不同的叠加方法由局部角度域像矩阵得到了反射角域CIG（CRAIG）和倾角域CIG（CDAIG）. 以SEG EAGE二维盐体模型为例,通过对CRAIG和CDAIG的对比,探讨了这两种角度域CIG的特点及其在地震偏移成像中的潜在应用.     

散射系数矩阵的计算及其应用

利用Gabor Daubechies（G? D）小波束域波场分解和传播在空间和方向上的双重局域性，提出了基于G? D小波束域叠前深度偏移进行角度域成像和计算局部散射系数矩阵的方法. 以简单分层模型为例，对不同探测系统的局部散射系数矩阵分布特征进行分析. 分析结果表明，在一定的探测系统几何布局下，由本文方法得到的局部散射系数矩阵能够较真实地反映局部结构的散射(或反射)特性. 通过局部散射系数矩阵进一步外推具有不同速度反差的水平界面随角度变化的反射系数，并估计界面的空间位置和倾角等说明局部散射系数矩阵的潜在应用.     

柯西波场及在角度域成像条件中的应用

角度域成像道集是叠前深度偏移的重要输出结果,它是偏移速度分析、各向异性分析和AVA分析等研究工作的基础.目前存在的角度域成像道集的生成方法受计算效率或角度分辨率的影响,仍然满足不了实际生产的要求.角度域成像道集的生成方法可以大致分为直接法和间接法两大类.在直接方法中,波矢量方向计算和局部平面波分解是两个最重要的内容,它们共同决定角度域成像条件的实现效率和角度域成像道集的质量.为了完善现有的角度域成像道集生成方法,本文提出一种新的波矢量计算方法和局部平面波分解方法.本文先用波动方程任意时刻的柯西条件构造一个只含原波场负频率成分的柯西波场,然后根据柯西波场在时间波数域的振幅谱计算波场的波矢量方向.该方法仅在需要计算波矢量方向的时刻合成柯西波场,不需要增加额外的数据读写操作,是一种高效的波矢量计算方法.本文还以柯西波场为基础提出一种高效的局部平面波分解方法,保证角度域成像条件的实现效率.结合柯西波场和局部平面波分解方法,本文最后给出一种新的角度域成像方法.文中最后的数值实验证明该方法得到的角度域成像道集具有理想的角度分辨率,可以反映地下构造的角度照明情况.     

2D共炮时间域高斯波束偏移

针对传统射线方法在奇异区成像精度不高,而2D频率域高斯波束叠前深度偏移需要计算成像点处每个频率的格林函数,影响计算效率的问题,本文通过使用复走时代替实走时,改变频率域下成像公式的积分顺序,给出了在时间域下进行高斯波束偏移的方法和计算公式.本文使用复杂数值模型验证了2D时间域高斯波束叠前偏移方法的正确性,并同传统射线偏移成像结果做了对比.对比结果表明时间域高斯波束偏移在成像精度上优于传统射线偏移.     

时-空局域化地震波传播方法:Dreamlet叠前深度偏移

提出了一种在时间和空间上完全局域化的波场分解和传播算法─dreamlet偏移方法.Dreamlet是一种脉冲-小波束形式的波场分解原子,它利用多维局部分解变换,把时空域波场映射到局部时间-频率-空间-波数相空间,并用局部相空间的传播算子(dreamlet算子)沿深度延拓.本文利用多维局部余弦变换实现dreamlet算法,分解后的波场系数和传播算子不仅有很好的稀疏性,且均为实数,也即波的传播和成像过程完全在实数域实现.文中推导了局部余弦基dreamlet波场分解和传播算子理论公式并将其应用于叠前深度偏移.在dreamlet相空间波的传播过程为稀疏矩阵相乘,而且延拓后的地表数据波场的有效时间长度随深度的增加不断减小,从而可以减少需要传播的波场系数.二维SEG/EAGE盐丘和SIGSBEE模型算例验证了理论推导的正确性,成像结果显示该方法在横向速度变化剧烈情况下有很好的精度.     

局域化相空间中的VSP偏移成像方法(英文)

对VSP资料进行偏移成像可提高井附近地下结构的成像分辨率。本文给出了一种基于局域化相空间波场分解的VSP偏移成像方法。此方法采用了基于Gabor-Daubechies紧标架的延拓算子(G-D延拓算子)及其高频渐近形式对相空间波场进行延拓;基于局部平面假设,提出了一种局部角度域相关成像条件。合成和实际VSP资料的偏移成像结果表明,在满足渐近展开的条件下利用G-D延拓算子的高频近似式能够有效的减少计算时间;采用局部角度域相关成像条件能够在不增加计算量的同时,有效减弱VSP成像剖面上的偏移假象。     

单向传播子的局部平面波分解以及在定向照明分析和目标特征成像中的应用

本文基于波场的局部平面波分解，得到空间和方向都具有局部性的局部平面波小波束，并由此对Fourier有限差分传播算子进行局部平面波分解.这种用小波束分解的单向传播算子，可被用来对采集系统进行定向照明分析和采集倾角响应估算，以分析采集系统布局和上覆结构对面向目标的叠前偏移成像质量的影响.进而，利用局部平面波成像矩阵及其成像相册，根据目标结构的主倾角方向，抽取并叠加部分局部平面波像，可以进行目标结构定向化特征成像.最后，对具有速度强横向变化上覆结构的陡倾角目标断层的复杂二维SEG-EAGE盐丘模型数据，进行采集系统目标照明分析和目标结构定向化特征成像.     

小波束叠前深度偏移和照明:基于Marmousi模型的测试

小波束源有很强的局部性和方向性,它可以很容易的实现局部照明和偏移,除此之外,它的偏移效果也要好于常规的偏移方法。本文介绍了基于框架理论的小波束叠前深度偏移的基本原理。我们解释了基于Gaussian函数的G-D框架。小波束分解提供了局部空间和方向上的信息。我们在小波域中,分别应用矩形窗和Gaussian窗合成了小波束源和小波束记录,并基于傅立叶有限差分算子,对marmousi模型进行了叠前深度偏移试算。通过对单个,以及多个小波束源偏移结果对比分析,验证了Gaussian小波束叠前深度偏移方法的有效性。     

面向目标的小束源照明和成像

结合小束源和Fourier传播子,应用Fourier传播子进行面向目标的小束源照明和成像.小束源的合成通过小波束变换中的小束函数获得,而Fourier传播子进行波场外推,完成小束源的波传播.小束源具有空间位置和方向的双局部特性,对于照明和成像,它更加灵活,并有更多的控制方法.通过照明分析,选择面向目标的有效照明小束源,进行部分小束源偏移,可以提供更好的成像质量和计算效率.作为数值试验,我们选用Fourier有限差分传播子,对 Marmousi模型和二维SEG-EAGE盐丘模型的数据,试验了小束源对目标结构的方向照明和成像.获得的结果表明,结合小束源和Fourier传播子进行面向目标的照明和成像是切实有效的.     

Seismic data compression by an adaptive local cosine/sine transform and its effects on migration

The local cosine/sine basis is a localized version of the cosine/sine basis with a window function which can have arbitrary smoothness. It has orthogonality and good time and frequency localization properties. The adaptive local cosine/sine basis is a best‐basis obtained from an overabundant library of cosine/sine packets based on a cost‐functional. We propose a 2D semi‐adaptive (time‐adaptive or space‐adaptive) local cosine transform (referred to as a 2D semi‐ALCT) and apply it to the SEG–EAEG salt model synthetic data set for compression. From the numerical results, we see that most of the important features of the data set can be well preserved even in the high compression ratio (CR=40:1) case. Using reconstructed data from the highly compressed ALCT coefficients (CR=40:1) for migration, we can still obtain a high‐quality image including subsalt structures. Furthermore, we find that the window partition, generated by the 2D semi‐ALCT, is well adapted to the characteristics of the seismic data set, and the compression capability of the 2D semi‐ALCT is greater than that of the 2D uniform local cosine transform (2D ULCT). We find also that a (32, 32) or (32, 64) minimum (time, space) window size can generate the best compression results for the SEG–EAEG salt data set.     

Prestack Migration/Imaging Using Synthetic Beam Sources and Plane Sources

The principle and procedures of beam-source and plane-source synthesis for wave-equation based imaging/migration are introduced and tested with Gabor-Daubechies and local-cosine beamlet prestack migration. Plane-source synthesis and imaging are treated as a special case of approximation to the beam-source synthesis and imaging. Beam sources with properly selected windows provide both space and direction localizations for illumination and are very flexible for target-oriented migration using wave-theory based methods. On the other hand, plane-source imaging can increase the computation efficiency of depth migration significantly, up to 5-10 times for the SEG-EAGE 2D salt model and Marmousi model, while maintains the image quality similar or comparable to that of point-source (common-shot) imaging/migration.     

基于Dreamlet变换的地震数据压缩理论与方法

为达到更加有效地表示地震数据的目的,仅仅将地震数据当作普通的图像数据处理是远远不够的,地震数据中蕴含的地震波的运动学特性也应作为重要因素而被考虑到.本文讨论了利用Dreamlet变换方法实现地震数据压缩的方法,并针对地震数据本身所蕴含的频散关系特性进一步提出了多尺度Dreamlet变换压缩方法.Dreamlet变换由2个一维局部谐波变换的张量积构成,它在提供地震波场时间-空间局部化性质的同时可以保留波场的运动学特性.通过对二维SEG/EAGE叠前、叠后数据的算例说明了Dreamlet变换用于地震数据压缩的有效性.利用压缩后的数据进行成像的结果更表明,与Curvelet变换方法相比,Dreamlet与多尺度Dreamlet方法可以提供更高的压缩比;在相同压缩比的条件下,使用Dreamlet与多尺度Dreamlet方法压缩重建后的数据进行成像能更好地保留成像结果中的重要结构.