利用矩阵分解理论的双程波方程叠前深度偏移方法

叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践.     

利用矩阵分解理论的双程波方程叠前深度偏移方法

叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践.     

构造单程波算子的优化方法及其在偏移中的应用

波动方程叠前深度偏移具有处理复杂构造的特点,高精度偏移方法是勘探地球物理研究方向之一,其单程波算子的构造是波动方程叠前深度偏移的重要内容.本文着重对混合域中单程波算子构造的两种优化方法进行分析,讨论他们与其他单程波算子构造方法之间相互关系,并用上述方法计算变速介质的脉冲响应及偏移Marmousi模型数据,以此说明该类方法的偏移精度.     

最优可分表示法三维波场延拓

适于复杂介质的高精度波场延拓算子是叠前深度偏移研究的重要内容。本文采用最优可分表示方法,运用正反傅立叶变换构造了三维单程波场延拓算子,算子实现了波数域变量与空间(速度)域变量分离。波数域内进行相移计算,在空间域对因介质横向变速引起的时移作修正。脉冲响应显示在区域内各速度的脉冲计算值与理论值基本一致,说明最优可分表示法叠前深度偏移可适用于强变速条件下复杂介质的成像需求。SEG／EAGE模型和实测数据的成像结果验证了本文方法对复杂构造的成像能力。     

角度域叠前深度偏移叠加算子优化

叠前深度偏移是解决精细复杂构造成像的一种最佳方法,尤其是基于角度域叠前深度偏移成像方法的研究,由于射线路径的唯一性,被认为是一种能够消除假象,解决射线多路径问题的最优化方法．文章从单程波波动方程出发,推导了中点 偏移距域双平方根波动方程波场延拓算子;利用局部偏移距与反射角的关系,推导了角度域共成像点道集的构建方法．同时,为有效解决叠前深度偏移计算量大、耗时长等计算瓶颈,引入OpenMpi并行策略,利用共享存储环境下的多线程协同作业,提高了叠前深度偏移算子的计算效率;在角度域,通过调整改进叠加算子,一定程度上减少了计算量,提高了计算效率．最后,通过理论模型验证了角度域叠前偏移成像算子的正确性,评价了并行算法的计算效率及有效性;并将其应用于SEG/EAGE岩丘模型,利用“沉降观测”理论得到了不同CMP的中点 偏移距域共成像点道集和角度域共成像点道集,对比分析了两种共成像点道集的优缺点;并通过限制成像叠加角度的方法,改善了复杂模型的成像结果,验证了角度域叠前深度偏移成像算法对复杂模型的适应性．     

衰减干涉假象的海洋一次波与多次波同时成像方法

海水与空气间的强波阻抗差使得海洋地震资料普遍发育自由表面相关多次波,如何利用好多次波所携带的有效信息已成为提高海洋地震资料成像品质的新突破点.基于面炮偏移的一次波与多次波同时成像方法能够避免多次波预测精度的影响,但是,正向传播的震源子波与反向延拓的自由表面相关多次波所产生的干涉假象严重制约了该技术的应用,本文提出了一种基于单程波偏移算子,可在成像域压制干涉假象的一次波与多次波同时成像方法.其中包含了三个步骤:第一,传统单程波偏移成像方法中的震源子波替换为一次波、多次波与震源子波,初始上行延拓波场为一次波与多次波,基于单程波算子的波场延拓与互相关成像条件的应用得到包含干涉假象的一次波与多次波同时成像;第二,以子波为震源,自由表面相关多次波为记录,按照传统单程波偏移成像方法得到干涉假象;第三,基于最小二乘匹配滤波算法,将第一步的成像结果与第二步的干涉假象进行匹配相减,得到干涉假象衰减后的一次波与多次波同时成像,避开了由于实际资料子波无法准确提取而造成一次波与多次波对成像能量级的不一致性.Sigsbee2B模型测试验证了本方法的有效性,并在我国某探区深海实际资料处理中得到了成功应用,深层基底得到了清晰刻画,并且照明均衡度明显改善.     

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

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

三维各向异性介质中的波动方程叠前深度偏移方法

基于三维VTI各向异性介质的频散关系,构建波数项和空间项分离的单程波算子表达式,以优化算法,确定算子的待定系数,实现广角逼近三维VTI介质的广义相移算子,发展了可灵活处理强或弱各向异性介质的波动方程叠前深度偏移方法.文中同时也针对其工业应用建议了三维VTI各向异性介质中可提高计算效率的频率相关变步长波场深度延拓算法及稀...     

基于时空局域化dreamlet单程波算子的观测系统沉降法偏移

Dreamlet偏移成像目的是探索一类能够对地震波场和单程波传播算子同时分解和压缩的理论和方法,也即实现在压缩域的传播与成像、地震数据在时间和空间的非平稳性质,决定了要实现地震数据的有效稀疏表示,分解方法必须在时间和空间上同时具有局域化性质.Dreamlet由时间和空间局部分解原子的张量积构成,可以看作一种脉冲-小波束形式的波场分解原子.时空局域化的dreamlet单程波传播算子在对波场沿深度方向延拓时,地震数据在时间轴上总是向同一方向流动.随着深度的增加,部分用于成像浅层结构的数据归位至其空间位置后被dreamlet算子丢弃,波场的有效记录时间变短,每一步用于波场延拓的计算量也相应下降.为了充分发挥这一优势,本文介绍dreamlet传播算子与观测系统沉降法偏移相结合的理论与方法.观测系统沉降法偏移每一步都将记录到的所有数据向下延拓,沉降后的波场等效于把源和检波器都放置在目标深度所能接收的反射数据.Dreamlet观测系统沉降过程只保留用于成像观测系统下部地质结构的有效数据,自动丢弃已经用于成像观测系统上部而对下部成像没有贡献的信号.本文通过二维SEG/EAGE叠后和Marmousi叠前数据算例展示了dreamlet传播算子应用于观测系统沉降法偏移的这一特点.数值算例结果显示,在不影响成像质量的前提下,该偏移方法能够有效减少传播数据量,为发展一种快速高效的偏移方法提供了新的思路.     

起伏地表条件下基于复Pade逼近的叠前深度偏移

叠前深度偏移是解决复杂地表和复杂构造地震成像的有效技术,而波场“直接下延”法实现了复杂地表条件下的地震成像.基于上述成果,结合高精度的波场延拓算子,本文提出了一种新的叠前深度偏移方法,这种方法是在波场延拓时,对声波方程中的平方根项进行复Pade逼近,通过推导得到基于复Pade逼近的傅里叶有限差分算子,结合波场“直接下延”法,实现了起伏地表条件下的叠前深度偏移,该算法减少了偏移噪音,从而得到准确、稳定的偏移成像结果.通过理论模型试算和实际资料试处理,验证了该方法的有效性.     

Mars in depth

P Lognonné, T Spohn and D Giardini describe the state of knowledge of Mars' interior and make the case for deep seismology to address the unanswered questions.     

深度震相sSmS特征及其在震源深度确定中的应用

通过一系列理论地震图模拟,研究了莫霍面超临界反射深度震相sSmS的特征,分析了影响该震相的各种因素。结果表明,SmS和sSmS属于高频波,一般情况下在高频段(1Hz左右)可被清晰地观测到;而在更长周期的地震图上,SmS和sSmS的强度比S波或者S多次波弱,不易辨认;地壳结构复杂地区且震源深度较浅时,sSmS震相也不容易被观测到。本文以2011年6月20日腾冲MS5.2地震为研究实例,利用sSmS深度震相确定其震源深度为6km,与其它方法所得结果一致。在利用深度震相测定震源深度的研究中,sSmS震相可以作为震源深度精确测定的手段之一。     

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the contrary, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending on parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). This effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance. Published 2020. This article is a U.S. Government work and is in the public domain in the USA     

Diffraction imaging in depth

High resolution imaging is of great value to an interpreter, for instance to enable identification of small scale faults, and to locate formation pinch-out positions. Standard approaches to obtain high-resolution information, such as coherency analysis and structure-oriented filters, derive attributes from stacked, migrated images. Since they are image-driven, these techniques are sensitive to artifacts due to an inadequate migration velocity; in fact the attribute derivation is not based on the physics of wave propagation. Diffracted waves on the other hand have been recognized as physically reliable carriers of high- or even super-resolution structural information. However, high-resolution information, encoded in diffractions, is generally lost during the conventional processing sequence, indeed migration kernels in current migration algorithms are biased against diffractions. We propose here methods for a diffraction-based, data-oriented approach to image resolution. We also demonstrate the different behaviour of diffractions compared to specular reflections and how this can be leveraged to assess characteristics of subsurface features. In this way a rough surface such as a fault plane or unconformity may be distinguishable on a diffraction image and not on a traditional reflection image.
We outline some characteristic properties of diffractions and diffraction imaging, and present two novel approaches to diffraction imaging in the depth domain. The first technique is based on reflection focusing in the depth domain and subsequent filtering of reflections from prestack data. The second technique modifies the migration kernel and consists of a reverse application of stationary-phase migration to suppress contributions from specular reflections to the diffraction image. Both techniques are proposed as a complement to conventional full-wave pre-stack depth migration, and both assume the existence of an accurate migration velocity.     

Velocity–depth and time–depth relationships for a decompacted uplifted unit

A linear instantaneous velocity model is used to describe the velocity variations in an uplifted unit that has been partly decompacted as a result of the reduction in overburden that often accompanies uplift. The model results in a series of equations for deriving values for the function parameters in the velocity–depth and the time– depth domains and for carrying out time‐to‐depth conversions. The formulation uses the base of the unit as a reference level to generate the reference datum from a combination of the depth of the base of the unit and a parameter that represents the decompaction factor.     

Heat flow at standard depth

Secular and long-term periodic changes in surface temperature cause perturbations to the geothermal gradient which may be significant to depths of at least 1000 m, and major corrections are required to determine absolute values of heat flow from the Earth's interior. However, detailed climatic models remain contentious and estimates of error in geothermal gradients differ widely. Consequently, regions of anomalous heat flow which could contain geothermal resources may be more easily resolved by measuring relative values at a standard depth (e.g. 100 m) so that all data are subject to similar corrections.Regional heat flow data obtained in existing deep holes show reasonable correlation with values determined at shallow depth. Hence geothermal resources of low enthalpy can be characterised by extrapolating temperatures from relative heat flow data readily obtained from shallow boreholes. Regional control can be provided by casing deep boreholes drilled for other purposes.For routine geothermal exploration, borehole temperatures can be measured using gradient probes with fixed sensor separation (e.g. 5 m), allowing very accurate determinations of the geothermal gradient at a single depth. Values of relative heat flow can then be obtained after determining the thermal resistivity of the corresponding core interval. Sampling errors can be minimised by multiple determinations of thermal conductivity over the complete interval.     

Resolution of prestack depth migration

The resolution of a general 3-D common-shot elastic prestack depth migration in a heterogeneous anisotropic medium is studied approximately, using the ray theory. It is demonstrated that the migrated section can approximately be represented by the convolution of the reflectivity function with the corresponding local resolution function. Alternatively, it can also be approximately represented by the convolution of the spatial distribution of the weak-contrast displacement reflection-transmission coefficient with the corresponding local resolution function. The derived explicit approximate equations enable us to predict the migration resolution approximately without doing the whole and expensive migration. The equations are applicable to 3-D elastic migrations in 3-D isotropic or anisotropic, heterogeneous velocity models. Both the reflectivity function and the spatial distribution of the weak-contrast displacement reflection-transmission coefficient approximately determine the linear combination of the perturbations of elastic moduli and density to which the migrated section is sensitive. The imaged linear combination of the perturbations of elastic parameters depends on the selection of the polarizations (wave types) of the incident and back-propagated wavefields and on the directions of propagation. The resolution of the linear combination of the perturbations of elastic moduli and density in the migrated section is determined by the above mentioned local resolution functions. The local resolution functions depend on the aperture and on the imaging function. The imaging function is determined by the source time function and by the form of the imaging functional. The local resolution functions are considerably sensitive to heterogeneities. The local resolution functions in elastic media differ from their acoustic counterparts, especially by the existence of converted scattered waves in elastic media.     

Decoupled elastic prestack depth migration

This paper presents a new decoupled form of the formula for common-shot or common-receiver amplitude-preserving elastic prestack depth migration (PreSDM), which can be used for estimating angle-dependent elastic reflection coefficients in laterally inhomogeneous anisotropic media. The multi-shot or multi-receiver extension of this formula is suitable for automated prestack amplitude-versus-angle (AVA) elastic inversion of ocean-bottom cable (OBC), walkaway VSP (WVSP) or standard towed-cable data at any subsurface location. The essence of the theory is a systematic application of the stationary-phase principle and high-frequency approximations to the basic elastic Green's theorem. This leads to nonheuristic explicit wave mode decoupling and scalarization of vector PreSDM. Used in combination, ray-trace and finite-difference (FD) eikonal solvers create a useful tool to calculate accurate Green's function travel time and amplitude maps. Examples of synthetic OBC data and applications to field WVSP data show that the new imaging technique can produce a clear multi-mode elastic image.     

High-resolution depth imaging with sparseness-constrained inversion

An imaging technique is developed which exceeds the resolution limitation prescribed by conventional seismic imaging methods. The high‐resolution imaging is obtained by introducing a sparseness‐constrained least‐squares inversion into the imaging process of prestack depth migration. This is implemented by a proposed interference technique. In contrast to conventional depth migration, a decomposed signal or combined event, instead of the source wavelet, is needed in the proposed scheme. The proposed method aims to image a small local region with a higher resolution using the prestack data set. It should be applied following conventional depth imaging if a higher resolution is needed in a target zone rather than replacing the conventional depth imaging for the entire medium. Synthetic examples demonstrate the significant improvements in the resolution using the proposed scheme.     

Prestack depth imaging via model-independent stacking

Most seismic reflection imaging methods are confronted with the difficulty of accurately knowing input velocity information. To eliminate this, we develop a special prestack depth migration technique which avoids the necessity of constructing a macro-velocity model. It is based upon the weighted Kirchhoff-type migration formula expressed in terms of model-independent stacking velocity and arrival angle. This formula is applied to synthetic sub-basaltic data. Numerical results show that the method can be used to successfully image beneath basalts.