地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.     

三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.     

模拟退火方法在三维速度模型地震波走时反演中的应用

采用块状建模以及三角形拼接的界面描述方式,并通过立方体速度网格线性插值获得块体内部的速度分布。正演过程中采用逐段迭代射线追踪方法计算三维复杂地质模型中的射线走时,并采用模拟退火方法进行了三维模型中的地震波走时反演研究。模型测试结果表明,使用的射线追踪和走时反演算法有效。     

南北地震带南段远震P波走时层析成像研究

南北地震带南段位于青藏高原东南缘,是青藏高原与扬子克拉通的过渡地带.本文收集了该区域内90个固定台站和356个流动台站的远震波形数据,采用波形互相关方法拾取了88691个P波走时残差数据,应用FMTT(Fast Marching Teleseismic Tomography)层析成像方法获取了南北地震带南段深部的三维P波速度结构.结果显示了研究区深部的结构具有显著的不均匀性:腾冲火山地区深部400km以浅的深度内分布着明显的低速异常;四川盆地西南部下方300km内具有较强的高速异常;在上地幔顶部,沿川滇菱形块体周边的大型断裂带及川滇菱形块体南端分布着显著的低速异常,这些低速异常为青藏高原物质向东南方向挤出提供了必要的通道;保山地块下方存在一东倾的高速异常带,该高速异常带为印度板块岩石圈向东俯冲的体现.     

Seismic traveltime inversion of 3D velocity model with triangulated interfaces

Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks, which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block, and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization, which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models. The damped least squares method is employed in seismic traveltime inversion, which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.     

三维复杂介质中转换波走时快速计算

复杂介质中转换波走时计算是多波勘探地震学中重要内容之一.本项研究利用惠更斯原理和费玛原理,获得了三维复杂介质中转换波快速计算的改进型最小走时树方法.其中,在保证精度的条件下,为了提高三维转换波走时计算效率,首先对初至波最小走时树基本算法进行了改进.本方法通过将转换波分为上、下行波分别进行射线追踪以实现三维转换波走时的快速计算.模型计算表明,方法的计算速度快,而且稳定性强,对多波地震勘探具有较大的应用价值.     

Influence of background heterogeneity on traveltime shifts for compacting reservoirs

Compaction induced by pore‐pressure decrease inside a reservoir can be monitored by measuring traveltime shifts of reflection events on time‐lapse seismic data. Recently we introduced a perturbation‐based formalism to describe traveltime shifts caused by the 3D stress‐induced velocity field around a compacting reservoir. Application of this method to homogeneous background models showed that the offset variation of traveltime shifts is controlled primarily by the anisotropic velocity perturbations and can provide valuable information about the shear and deviatoric stresses. Here, we model and analyse traveltime shifts for compacting reservoirs whose elastic properties are different from those of the surrounding medium. For such models, the excess stress is influenced primarily by the contrast in the rigidity modulus μ across the reservoir boundaries. Synthetic examples demonstrate that a significant (25% or more) contrast in μ enhances the isotropic velocity perturbations outside the reservoir. Nevertheless, the influence of background heterogeneity is mostly confined to the reservoir and its immediate vicinity and the anisotropic velocity changes are still largely responsible for the offset dependence of traveltime shifts. If the reservoir is stiffer than the host rock, the background heterogeneity reduces anisotropic velocity perturbations inside the reservoir but increases them in the overburden. As a result, in this case, the magnitude of the offset variation of traveltime shifts is generally higher for reflections from interfaces above the reservoir. We also study compaction‐induced stress/strain and traveltime shifts for a stiff reservoir embedded in a softer layered model based on velocity profiles from the Valhall Field in the North Sea. Despite producing discontinuities in strain across medium interfaces, horizontal layering does not substantially alter the overall behaviour of traveltime shifts. The most pronounced offset variation of traveltime shifts is observed for overburden events recorded at common midpoints close to the reservoir edges. On the whole, prestack analysis of traveltime shifts should help better constrain compaction‐induced velocity perturbations in the presence of realistic background heterogeneity.     

Transmission + reflection anisotropic wave-equation traveltime and waveform inversion

A transmission + reflection wave-equation traveltime and waveform inversion method is presented that inverts the seismic data for the anisotropic parameters in a vertical transverse isotropic medium. The simultaneous inversion of anisotropic parameters and ε is initially performed using transmission wave-equation traveltime inversion method. Transmission wave-equation traveltime only provides the low-intermediate wavenumbers for the shallow part of the anisotropic model; in contrast, reflection wave-equation traveltime estimates the anisotropic parameters in the deeper section of the model. By incorporating a layer-stripping method with reflection wave-equation traveltime, the ambiguity between the background-velocity model and the depths of reflectors can be greatly mitigated. In the final step, multi-scale full-waveform inversion is performed to recover the high-wavenumber component of the model.  We use a synthetic model to illustrate the local minima problem of full-waveform inversion and how transmission and reflection wave-equation traveltime can mitigate this problem. We demonstrate the efficacy of our new method using field data from the Gulf of Mexico.     

起伏地形下的高精度反射波走时层析成像方法

全球造山带及中国大陆中西部普遍具有强烈起伏的地形条件.复杂地形条件下的地壳结构成像问题像一面旗帜引领了当前矿产资源勘探和地球动力学研究的一个重要方向.深地震测深记录中反射波的有效探测深度可达全地壳乃至上地幔顶部,而初至波通常仅能探测上地壳浅部.为克服和弥补初至波探测深度的不足,本文基于前人对复杂地形条件下初至波成像的已有研究成果,采用数学变换手段将笛卡尔坐标系的不规则模型映射到曲线坐标系的规则模型,并将快速扫描方法与分区多步技术相结合,发展了反射波走时计算和射线追踪的方法.进而利用反射波走时反演,实现起伏地形下高精度的速度结构成像,从而为起伏地形下利用反射波数据高精度重建全地壳速度结构提供了一种全新方案.数值算例从正演计算精度、反演中初始模型依赖性、反演精度、纵横向分辨率以及抗噪性等方面验证了算法的正确性和可靠性.     

Efficient traveltime compression for 3D prestack Kirchhoff migration

Kirchhoff 3D prestack migration, as part of its execution, usually requires repeated access to a large traveltime table data base. Access to this data base implies either a memory intensive or I/O bounded solution to the storage problem. Proper compression of the traveltime table allows efficient 3D prestack migration without relying on the usually slow access to the computer hard drive. Such compression also allows for faster access to desirable parts of the traveltime table. Compression is applied to the traveltime field for each source location on the surface on a regular grid using 3D Chebyshev polynomial or cosine transforms of the traveltime field represented in the spherical coordinates or the Celerity domain. We obtain practical compression levels up to and exceeding 20 to 1. In fact, because of the smaller size traveltime table, we obtain exceptional traveltime extraction speed during migration that exceeds conventional methods. Additional features of the compression include better interpolation of traveltime tables and more stable estimates of amplitudes from traveltime curvatures. Further compression is achieved using bit encoding, by representing compression parameters values with fewer bits.