Multidimensional GPR array processing using Kirchhoff migration

We compare the ability of several practical ground-penetrating radar (GPR) array processing methods to improve signal-to-noise ratio (SNR), increase depth of signal penetration, and suppress out-of-plane arrivals for data with SNR of roughly 1. The methods include two-dimensional (2-D) monostatic, three-dimensional (3-D) monostatic, and 3-D bistatic Kirchhoff migration. The migration algorithm is modified to include the radiation pattern for interfacial dipoles. Results are discussed for synthetic and field data. The synthetic data model includes spatially coherent noise sources that yield nonstationary signal statistics like those observed in high noise GPR settings. Array results from the model data clearly indicate that resolution and noise suppression performance increases as array dimensionality increases. Using 50-MHz array data collected on a temperate glacier (Gulkana Glacier, AK), we compare 2-D and 3-D monostatic migration results. The data have low SNR and contain reflections from a complex, steeply dipping bed. We demonstrate that the glacier bed can only be accurately localized with the 3-D array. In addition, we show that the 3-D array increases SNR (relative to a 2-D array) by a factor of three.     

探地雷达小波域三维波动方程偏移

阐述了矩阵多分辨分析理论中的标准形式与非标准形式,并以Hilbert算子为例,说明了算子多分辨表示的压缩效果,为小波域偏移算法奠定了理论基础.从三维雷达波动方程出发,利用爆炸反射原理和浮动坐标变换,推导出三维探地雷达波动方程差分格式,并通过方程分裂算法及小波多分辨算法,在小波域求解波场外推矩阵,进而得到探地雷达小波域三维波动方程偏移算法,在此基础上,开发了探地雷达小波域偏移处理程序,并把该程序应用于三个球体空洞的3-D正演结果及实际的雷达数据中,通过对比偏移处理前后的雷达资料,得知该三维偏移算法能使3-D正演剖面中的反射波归位、绕射波收敛,极大地提高了雷达剖面的分辨率,有利于探地雷达资料的地质解释.     

Reverse time migration of ground-penetrating radar with full wavefield decomposition based on the Hilbert transform

The conventional reverse time migration of ground-penetrating radar data is implemented with the two-way wave equation. The cross-correlation result contains low-frequency noise and false images caused by improper wave paths. To eliminate low-frequency noise and improve the quality of the migration image, we propose to separate the left-up-going, left-down-going, right-up-going and right-down-going wavefield components in the forward- and backward-propagated wavefields based on the Hilbert transform. By applying the reverse time migration of ground-penetrating radar data with full wavefield decomposition based on the Hilbert transform, we obtain the reverse time migration images of different wavefield components and combine correct imaging conditions to generate complete migration images. The proposed method is tested on the synthetic ground-penetrating radar data of a tilt-interface model and a complex model. The migration results show that the imaging condition of different wavefield components can highlight the desired structures. We further discuss the reasons for incomplete images by reverse time migration with partial wavefields. Compared with the conventional reverse time migration methods for ground-penetrating radar data, low-frequency noise can be eliminated in images generated by the reverse time migration method with full wavefield decomposition based on the Hilbert transform.     

再论地震数据偏移成像

利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高.     

Integration of ground-penetrating radar and microgravimetric methods to map shallow caves

Ground-penetrating radar (GPR) and microgravimetric surveys have been conducted in the southern Jura mountains of western Switzerland in order to map subsurface karstic features. The study site, La Grande Rolaz cave, is an extensive system in which many portions have been mapped. By using small station spacing and careful processing for the geophysical data, and by modeling these data with topographic information from within the cave, accurate interpretations have been achieved.The constraints on the interpreted geologic models are better when combining the geophysical methods than when using only one of the methods, despite the general limitations of two-dimensional (2D) profiling. For example, microgravimetry can complement GPR methods for accurately delineating a shallow cave section approximately 10×10 m² in size. Conversely, GPR methods can be complementary in determining cavity depths and in verifying the presence of off-line features and numerous areas of small cavities and fractures, which may be difficult to resolve in microgravimetric data.     

Converted‐wave beam migration with sparse sources or receivers

Gaussian beam depth migration overcomes the single‐wavefront limitation of most implementations of Kirchhoff migration and provides a cost‐effective alternative to full‐wavefield imaging methods such as reverse‐time migration. Common‐offset beam migration was originally derived to exploit symmetries available in marine towed‐streamer acquisition. However, sparse acquisition geometries, such as cross‐spread and ocean bottom, do not easily accommodate requirements for common‐offset, common‐azimuth (or common‐offset‐vector) migration. Seismic data interpolation or regularization can be used to mitigate this problem by forming well‐populated common‐offset‐vector volumes. This procedure is computationally intensive and can, in the case of converted‐wave imaging with sparse receivers, compromise the final image resolution. As an alternative, we introduce a common‐shot (or common‐receiver) beam migration implementation, which allows migration of datasets rich in azimuth, without any regularization pre‐processing required. Using analytic, synthetic, and field data examples, we demonstrate that converted‐wave imaging of ocean‐bottom‐node data benefits from this formulation, particularly in the shallow subsurface where regularization for common‐offset‐vector migration is both necessary and difficult.     

傅里叶有限差分法保幅叠前深度偏移方法

地震数据中饱含有丰富的走时信息和振幅信息. 为解决传统偏移方法中几何扩散和入射角变化引起的振幅误差问题,本文提出了一种实用的波动方程保幅地震偏移方法. 该方法从全声波方程出发进行单程波保幅分解,得到直观、高效率的直接面对地震波传播波场的压力分量进行延拓的保幅偏移单程波方程,进而推导出一个含有6项的傅里叶有限差分法保幅偏移的算子方程;修改边界条件和成像条件,使修改后的边界条件和成像方程中考虑振幅补偿,从而从三方面补偿几何扩散损失和入射角变化对振幅的影响. 脉冲响应测试、单炮记录的数值试验以及SEG/EAGE盐丘模型的叠前偏移结果表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析(如AVO/AVA)提供更真实的地震信息.     

Modelling sand wave migration in shallow shelf seas

Sand waves form a prominent regular pattern in the offshore seabed of sandy shallow seas. The positions of sand-wave crests and troughs slowly change in time. Sand waves are usually assumed to migrate in the direction of the residual current. This paper considers the physical mechanisms that may cause sand waves to migrate and methods to quantify the associated migration rates. We carried out a theoretical study based on the assumption that sand waves evolve as free instabilities of the system. A linear stability analysis was then performed on a 2DV morphological model describing the interaction between the vertically varying water motion and an erodible bed in a shallow sea. Here, we disrupted the basic tidal symmetry by choosing a combination of a steady current (M₀) and a sinusoidal tidal motion (M₂) as the basic flow. We allowed for two different physical mechanisms to generate the steady current: a sea surface wind stress and a pressure gradient. The results show that similar sand waves develop for both flow conditions and that these sand waves migrate slowly in the direction of the residual flow. The rates of migration and wavelengths found in this work agree with theoretical and empirical values reported in the literature.     

横向各向同性介质拟声波方程及其在逆时偏移中的应用

地下岩石的速度各向异性影响地震波的传播与成像.横向各向同性(TI)介质为最普遍的等效各向异性模型.引入TI介质拟声波方程可以避免复杂的弹性波方程求解以及各向异性介质波场分离,以满足对纵波成像的实际需要.本文从垂直横向各向同性(VTI)介质弹性波方程出发,推导出正应力表达的拟声波方程以及相应的纵波分量的表达式,进而分析从频散关系得到的拟声波方程的物理意义,而后将拟声波方程扩展到更一般的倾斜横向各向同性(TTI)介质中.波前快照与群速度平面的对比验证了拟声波方程可以很好地近似描述qP波的运动学特征.在此基础上,将拟声波方程应用在逆时偏移中并与其特例声波近似方程进行对比,讨论了计算效率、稳定性等实际问题.数值试验表明VTI介质情况下采用声波近似方程可以提高计算效率,而TTI介质qP-qSV波方程则在效率相当的情况下可以保证稳定性.SEG/HESS模型和逆冲模型逆时偏移试验验证了本文TI介质拟声波方程的实用性.     

雷达波偏移方法研究

为了解决由于点状体的绕射和倾斜界面造成的失真现象,本文从探地雷达偏移处理的基本原理出发,分别介绍了以射线理论为基础的绕射扫描偏移和以波动方程为出发点的克希霍夫积分偏移、有限差分偏移以及有限元偏移．随着计算机的高速发展,运算速度得到了明显提高,目前,以波动方程为基础的偏移方法逐渐占据了主导地位,由于它的高精度等很多优点受到了人们的欢迎．总结了各种偏移方法的优缺点和适用范围,为探地雷达进一步的偏移研究奠定了理论基础．     

波动方程偏移成像阴影的照明补偿

受地下复杂构造和地震数据采集系统的影响,使地震波对地下目标的照明出现不均匀性,地震采集系统难以有效地获取地下某些目标的反射信息,进而使数据偏移成像在这些目标体上出现成像阴影. 根据波场和Green函数的窗口Fourier框架展开,利用角度域波动方程偏移成像和波动方程照明分析,并结合波动方程反演理论,提出一种角度域波动方程偏移成像阴影照明补偿方法. 这种补偿方法能同时考虑地震数据采集系统和波场传播路径对偏移成像的影响,消除复杂构造区的偏移成像阴影,改进波动方程叠前深度偏移成像在复杂构造区的成像效果.     

小尺度散射体稀疏最小二乘逆时偏移方法研究

地震勘探方法在深部固体矿产资源勘探中发展潜力巨大,同时也面临挑战.由于固体矿产资源地下分布呈现陡峭构造、尺度小,物性差异小的特点,常规偏移方法对小尺度矿体成像的分辨率提高有限.本文研究了一种基于稀疏促进约束的最小二乘逆时偏移方法.首先,将非均匀分布的矿体等效为随机介质,建立小尺度扰动的矿体模型;其次,改进现有最小二乘偏移方法,以稀疏模型为先验信息约束成像结果,并通过Curvelet变换压缩成像空间,经过多次迭代计算,可以提高小尺度散射体的成像分辨率;再次,对炮域记录进行随机震源编码,减少成像所需的炮集个数,通过稀疏促进约束条件,降低串扰噪声引起的成像误差.通过庐枞金属矿模型数值计算,验证本文方法可以较好的成像包含小尺度散射体的金属矿地质模型.

     

矢量黏弹性衰减补偿高斯束偏移

基于弹性波动理论的多波多分量高斯束偏移具有计算效率高和成像准确等优点.但是目前此方法没有考虑实际地下介质的黏弹性对地震波传播的影响,从而无法补偿能量衰减和校正相位畸变,这使得该方法对一些含高黏弹性地层的成像效果不佳.针对衰减区域的成像问题,本文提出一种黏弹性衰减补偿高斯束偏移方法,该方法以多波多分量矢量波场弹性高斯束偏移方法为基础,在偏移过程中沿射线路径通过引入品质因子Q来考虑黏弹性影响并进行衰减补偿.该方法能够在偏移过程中实现PP波和PS波的自动分离及分别成像.同时,本文给出了在矢量波场偏移过程中提取角度域共成像点道集的方法,以便用于成像质量控制,并为后续速度和黏弹性参数反演提供所需的数据.本文利用2D层状模型和洼陷模型进行了方法测试,其成像结果验证了本文所提出的黏弹性衰减补偿高斯束偏移方法的可行性和有效性.

     

Double plane‐wave reverse‐time migration

We develop a new time‐domain reverse‐time migration method called double plane‐wave reverse‐time migration that uses plane‐wave transformed gathers. Original shot gathers with appropriate data acquisition geometry are double slant stacked into the double plane‐wave domain with minimal slant stacking artefacts. The range of plane‐wave components needed for migration can be determined by estimating the maximum time dips present in shot gathers. This reduces the total number of input traces for migration and increases migration efficiency. Unlike the pre‐stack shot‐profile reverse‐time migration where the number of forward propagations is proportional to the number of shots, the number of forward propagations needed for the proposed method remains constant and is relatively small even for large seismic datasets. Therefore, the proposed method can improve the efficiency of the migration and be suitable for migrating large datasets. Double plane‐wave reverse‐time migration can be performed for selected plane‐wave components to obtain subsurface interfaces with different dips, which makes the migration method target oriented. This feature also makes the method a useful tool for migration velocity analysis. For example, we are able to promptly obtain trial images with nearly horizontal interfaces and adjust velocity models according to common image gathers. Seismic signal coming from steeply dipping interfaces can be included into the migration to build images with more detailed structures and higher spatial resolution as better velocity models become available. Illumination compensation imaging conditions for the proposed method are also introduced to obtain images with balanced amplitudes.     

单程波角度域内压制多次波偏移假象

多次波偏移中的假象主要来自于不同地震事件之间的互相关,由于这种互相关满足成像条件,很难直接在偏移过程中去除.但是对于准确的速度模型,真实的成像结果在角度域内应该是平直的.根据这个判断准则,可以在角度域内移除多次波偏移中的假象.本文以数据自相关偏移为例,提出了在单程波多次波偏移中移除假象的主要流程：首先在在单程波偏移过程中高效地提取角度域共成像点道集,然后对角度域共成像点道集应用高分辨率的抛物线型Radon变换,用合适的切除函数处理后,反变换回到角度域,最后叠加各个角度成分,得到偏移结果.Marmousi模型的合成数据测试表明,这种方法可以很好地压制多次波偏移过程中产生的假象,有效地提高成像结果的信噪比.     

Coal top detection by conductively guided borehole radar wave imaging for open cut blast-hole drilling

Damage to the top of coal seams, caused by incorrect blast stand-off distances, results in coal losses of up to 10–15% to the Australian open cut coal mining operations. This is a serious issue to be addressed. Here we propose to use a new forward-looking imaging technique based on the borehole radar technology to predict the coal seam top in real time while drilling blast holes. This is achieved by coupling the conventional borehole radar waves on to a steel drill rod to induce a guided wave along the axial drill rod. The drill rod ahead of the borehole radar behaves as a forward-looking antenna for the guided waves. Both numerical modelling and field trials simulating a drill rod as an antenna are used to investigate the feasibility of the proposed technique for prediction of the coal top under typical open cut environments. Numerical modelling demonstrated that conductivity of the overburden is the most important factor affecting our ability to see coal seams ahead of the drill bit, the guided borehole radar waves could be used for top coal prediction and a theoretical prediction error less than 10 cm and a forward-looking capability of 4–6 m can be achieved. Field trials at Australian open cut coal mines also demonstrated that guided borehole radar waves can be observed and used for prediction of coal top ahead of drill bit during blast-hole drilling in resistive, open cut environments (the average resistivity should be higher than 75 Ωm).     

用地质雷达数据资料反演二维地下介质的方法

从二维麦克斯韦方程组出发推导出反演介电常数和电导率等二维介质物性参数的反演公式.反演的步骤是: 建立初始猜测模型,利用电磁波时间域有限差分法模拟正演数据,用正演数据与观测数据之间的数据残差建立目标函数,通过引入一个由麦克斯韦方程计算的伴随场,将目标函数对介质参数的导数表示成显式形式,应用最优化理论得出对初始猜测模型的修改,用共轭梯度法迭代,最终得到反演结果.用合成数据反演具有粗糙地表的非导电介质的介电常数,用实验数据同时反演介电常数和电导率,并比较了麦克斯韦方程反演结果与声波方程反演结果、波动方程偏移剖面的差异.     

基于波场传播方向的波动方程偏移有效去噪方法(英文)

本文研究分析了双程波波动方程偏移成像中广泛存在的三种主要噪声,特别是针对过去研究中没能很好解决的存在于高速盐丘悬垂边界附近的射线状噪声,提出了基于优化成像条件的有效去噪方法。射线状噪声主要来自于震源一侧波场的下行透射波分量和接收阵列一侧波场的上行散射波互相关成像。这一部分能量具有较强的互相关性,但并不携带真实的反射面信息。它广泛存在叠前偏移成像中,与信号在强度上同量级。多数情况下偏移成像中的相关噪声由方向性传播的波场能量产生。利用波场梯度得到的波场传播角度,可以分离出噪声对应的波场能量,并在成像条件中减去。采用这一方法可以有效地去除多种噪声,包括直达波噪声、散射波噪声和射线状噪声。该去噪方法不依赖波场外推算子,在需要时可以方便地运用到几乎所有的波动方程偏移中去。并且该去噪方法针对噪声的物理根源,对信号的损害很小。对去噪后的偏移成像结果额外地进行波数域滤波处理,可以进一步提高叠加图像的质量。这一去噪方法在超广角单程波偏移成像中取得良好效果,我们同时期待其在其他双程波波动方程偏移特别是逆时偏移(RTM)中的成功运用。     

弹性波逆时偏移中的稳定激发振幅成像条件

本文针对弹性波逆时偏移,提出稳定的激发振幅成像条件.在震源波场的正向传播过程中,计算每个网格点的能量,并保存最大能量密度的时刻和相应的波场值;在检波器波场的逆时传播过程中,在每个网格点提取最大能量密度时刻的检波器波场值,并利用保存的最大能量震源波场做归一化,获得角度依赖的反射系数成像剖面.相比于归一化互相关成像条件,该成像条件在震源波场的正向传播过程中无需存储波场快照,节省大量磁盘空间和I/O吞吐任务,提高了计算效率;相比于弹性波的激发时间成像条件,该成像条件自动校正了水平分量在震源两侧的极性反转,在多炮叠加时避免振幅损失.数值试验表明,与归一化成像条件相比,稳定激发振幅成像条件具有更小的计算量,偏移剖面的低频假象更弱,水平分量的成像能力更优,具有更高的空间分辨率.     

数据自相关多次波偏移成像

在常规偏移方法中一般都需要压制地震数据中的多次波,仅利用一次波信息成像,把自由表面反射的多次波视为噪声,但是在多次波中也包含着地下结构信息,应该将其充分利用到成像中来.事实上,已经有不少成像方法试图利用多次波信息,但是大部分方法都需要对多次波进行预测.本文提出了基于傅里叶有限差分偏移算子的数据自相关偏移方法.在这种偏移方法中,对含有一次波和多次波的地震数据,分别进行下行和上行延拓,然后直接利用常规的互相关成像条件成像.由于波场延拓采用了傅里叶有限差分算子,其计算效率高,能够很好地对复杂介质中的地震数据进行延拓.在数值试验中,使用了一个含散射点的三层模型和Marmousi模型.合成数据测试结果表明,这种方法可以对更大范围的地下构造成像,比常规的只利用一次波的傅里叶有限差分法照明度更好,并且在浅层可以提供更高的分辨率.我们提出的数据自相关策略易于实现且避免了繁杂的多次波预测,这对于复杂地下构造成像可能有着重大意义.