三维复杂地壳结构非线性走时反演

中国大陆中西部乃至全球造山带普遍具有复杂地壳结构.随着矿产资源勘探和深部探测研究的深入,探测造山带及盆山耦合区下方地壳精细结构正逐渐成为当前面临的巨大挑战.人工源深地震测深方法正越来越清晰地揭示出不同构造域地壳速度结构的基本特征,然而传统的层状结构模型参数化方法难以准确描述复杂地质模型,通常情况下多忽略速度结构的精细间断面且采用层边界平滑处理,难以满足地壳精细结构成像的发展要求.针对上述困难,本文采用最近发展的块状结构建模方案构建三维复杂地壳模型,基于逐段迭代射线追踪正演走时计算方法,推导了走时对三角形界面深度以及网格速度的偏导数,开展了非线性共轭梯度走时反演方法研究.发展了利用直达波和反射波等多震相走时数据对界面深度和网格速度的多参数联合反演方法,并引人不同种类震相数据的权系数和不同类型参数偏导数归一化的方法.数值算例表明,基于块状结构的非线性共轭梯度走时反演方法适用于复杂地壳结构模型,在利用人工源走时数据反演复杂地壳精细结构领域具有良好的应用前景.     

三维地质模型中地震波共轭梯度非线性走时反演

地震体波走时层析成像是探测地球内部速度结构的重要方法之一。基于三维块状建模以及三角形拼接的界面描述方式,结合快速高效的逐段迭代射线追踪方法,获得三维复杂地质模型中的地震射线路径与走时信息,采用共轭梯度非线性反演算法,进行地震波走时反演。实验结果表明共轭梯度反演算法在三维层状模型中具有较高的有效性。     

Two-point ray tracing in 3-D

Algorithm for determination of all two-point rays of a given elementary wave by means of the shooting method is presented. The algorithm is designed for general 3-D models composed of inhomogeneous geological blocks separated by curved interfaces. It is independent of the initial conditions for rays and of the initial-value ray tracer. The algorithm described has been coded in Fortran 77, using subroutine packages MODEL and CRT for model specification and for initial-value ray tracing.     

弯曲射线追踪中Dijkstra算法的改进与实现

文章针对图论中寻找最短路径的Dijkstra算法内存占用量大,效率低的缺点,对该算法进行了改进,修改后的算法计算效率是原来的四倍,内存使用量和图中节点数呈线性关系.在此基础上,用新算法求出了激发点和接收点的最短走时路径,并由激发接收点的旅行时结合联合迭代法对理论模型和实际场地进行了反演.结果表明：和直射线追踪相比,弯曲射线路径能更好地反演出地质体内部的速度场分布.理论模型和实际探测结果证实改进后的算法是有效的.     

Common-ray tracing and dynamic ray tracing for S waves in a smooth elastic anisotropic medium

Anisotropic common S-wave rays are traced using the averaged Hamiltonian of both S-wave polarizations. They represent very practical reference rays for calculating S waves by means of the coupling ray theory. They eliminate problems with anisotropic-ray-theory ray tracing through some S-wave slowness-surface singularities and also considerably simplify the numerical algorithm of the coupling ray theory for S waves. The equations required for anisotropic-common-ray tracing for S waves in a smooth elastic anisotropic medium, and for corresponding dynamic ray tracing in Cartesian or ray-centred coordinates, are presented. The equations, for the most part generally known, are summarized in a form which represents a complete algorithm suitable for coding and numerical applications.     

Nonlinear ray perturbation theory with its applications to ray tracing and inversion in anisotropic media

A comprehensive approach, based on the general nonlinear ray perturbation theory (Druzhinin, 1991), is proposed for both a fast and accurate uniform asymptotic solution of forward and inverse kinematic problems in anisotropic media. It has been developed to modify the standard ray linearization procedures when they become inconsistent, by providing a predictable truncation error of ray perturbation series. The theoretical background consists in a set of recurrent expressions for the perturbations of all orders for calculating approximately the body wave phase and group velocities, polarization, travel times, ray trajectories, paraxial rays and also the slowness vectors or reflected/transmitted waves in terms of elastic tensor perturbations. We assume that any elastic medium can be used as an unperturbed medium. A total 2-D numerical testing of these expressions has been established within the transverse isotropy to verify the accuracy and convergence of perturbation series when the elastic constants are perturbed. Seismological applications to determine crack-induced anisotropy parameters on VSP travel times for the different wave types in homogeneous and horizontally layered, transversally isotropic and orthorhombic structures are also presented. A number of numerical tests shows that this method is in general stable with respect to the choice of the reference model and the errors in the input data. A proof of uniqueness is provided by an interactive analysis of the sensitivity functions, which are also used for choosing optimum source/receiver locations. Finally, software has been developed for a desktop computer and applied to interpreting specific real VSP observations as well as explaining the results of physical modelling for a 3-D crack model with the estimation of crack parameters.     

基于波前面三角形网格剖分的波前重建法三维射线追踪

在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用.     

复杂地质层面的网格化模型及两点射线追踪方法

本文针对三维复杂介质的两点射线追踪的适用性和有效性,提出了网格模型和层状模型相结合的建模方式,文中利用规则网格化的层状模型对复杂地质单元进行数值描述,一方面避免了用纯粹复杂地质网格体描述时需要的大量内存,另一方面结合运用规则化的层状模型,必有利于提高计算效率,在计算单元自动生成的射线追踪方法基础上,本文结合波前法提出一种快速、有效的两点射线追踪方法,利用反射波波前在地表形成的射线方位角和倾角坐标体系,能快速、准确地收敛到所有从源点到接收点的射线路径,亦可有效地解决射线路径的多值解难题,从而可适用于复杂地质模型的特点.     

Theory of anisotropic dynamic ray tracing in ray-centred coordinates

A 4×4-propagator matrix formalism is presented for anisotropic dynamic ray tracing, including the propagation across curved interfaces. The computations are organised in the same way as in ervený's well-known isotropic propagator matrix formalism. Attention is paid to cases where double eigenvalues of the Christoffel matrix result in unstable expressions in the dynamic ray tracing system, but where geometrical spreading is well-behaved.     

三维非均匀地质模型中的逐段迭代射线追踪

地震射线追踪是地震定位、层析成像、偏移等领域的重要正演环节.随着这些领域研究的深入,针对传统的网格结构和层状结构在描述复杂地质模型遇到的很大困难,我们采用大小不等、形状各异的地质块组成的集合体来描述三维复杂地质模型,并用三角形面片来描述地质块之间的物性间断面,理论上可以描述任意复杂的地质模型.为适应任意非均匀速度分布的地质模型,基于费马原理,本文发展了与之相适应的逐段迭代射线追踪方法.该方法属于弯曲法范畴,对路径点采用一阶显式增量修正,相对于传统的迭代法,高效省时.数值试验表明,联合逐段迭代法和伪弯曲法的射线追踪扰动修正方案在三维复杂非均匀块状模型中有适用性和高效性.     

Grid travel-time tracing: Second-order method for the first arrivals in smooth media

A new computational scheme for calculating the first-arrival travel times on a rectangular grid of points is proposed. The new proposed method is of second-order accuracy. This means that the error of the calculated travel time is proportional to the second power of the grid spacing. The method should be sufficiently accurate for all applications in smooth seismic models. On the other hand, the method is not, in its present form, proposed for models with structural interfaces which make the method unstable and generate travel-time errors of the first order. Equations are also presented for the appropriate evaluation of the errors of calculated travel times to check their accuracy, and the proposed method is compared with other numerical methods. The method is developed, described and demonstrated in 2-D, but may also be extended to 3-D models and to general models with structural interfaces.     

椭球坐标系下多震相地震射线追踪

地震射线追踪方法技术在地震学领域有着较为广泛的应用,然而大多数算法建立在直角坐标系或球坐标系下,实际地球并非完美的球体,而是两极略扁的椭球体,因此,球坐标系下计算结果与真实情况存在一定误差.传统的做法一般是在球坐标系下进行计算,而后进行椭球校正.本文提出了一种直接在椭球体模型中采用分区多步最短路径算法进行多震相地震射线追踪的方法技术,实现了椭球坐标系下多震相地震波射线路径追踪和走时计算.与解析解的对比表明：该算法具有较高的计算精度,适用于任意形状的椭球体,且不需要进行额外的走时校正.数值模拟结果表明,计算所得P波和PcP反射波的走时与AK135走时表的误差小于0.1 s.当震中距较大时,使用球对称模型和椭球体模型计算所得的走时差异显著,说明采用椭球坐标系的必要性.

     

A rapid and accurate two-point ray tracing method in horizontally layered velocity model

Introduction The determination of the ray path between a source and a receiver is fundamental to many seismic problems such as earthquake location, travel time inversion, computation of synthetic seismograms, depth migration, and seismic tomography. The calculation of travel times and their derivatives depend on the determination of ray paths, which is the most time-consuming work in seismic tomography. Therefore, a fast and accurate ray tracing method is particularly important. Many ray-traci…     

基于三维层状介质模型地震层析成像正演网格研究

地震层析成像技术是研究地层速度结构的一个重要手段,该技术是利用震源到接收点的地震波走时或波形的观测数据,结合建立的数学物理模型来重建地层速度结构。而层析成像正演模拟的精度将直接影响反演出的速度与地层真实速度的拟合度,因此找出一种精度较高的层析成像正演模拟算法是非常重要的。本文针对三维层状介质,通过网格化模型,采用最小走时射线追踪方法展开层析成像正演研究,分别采用三角形网格、矩形网格和六边形网格进行模型参数化,并对比分析各自正演结果的精度,总结出针对三维层状介质模型的最佳网格划分方式。     

高分辨率非线性三维整体反演研究(英文)

高分辨率非线性三维整体反演方法是基于非线性理论,在层位控制下,将工区多井(或全部井)的测井数据与井旁地震道数据输入具有多输入多输出的网络,同时进行整体训练,可获得整个工区的自适应权函数,并建立综合非线性映射关系,并根据储层在纵横方向上的地质变化特征更新这种非线性映射关系,这样,就能对反演过程及其反演结果起到约束和控制的作用,从而获得稳定且分辨率高的地震反演剖面(速度反演剖面/波阻抗反演剖面/密度反演剖面),实现整体反演,该方法通过模型试算和实际资料处理,获得较好的地质效果,证明该方法精度高、实用性强,可用于储层的定量分析。     

Detailed reservoir inversion addressing geological problems in reservoir development

垦71开发区块的主要产油层系为新近系的馆陶组河流相和古近系的东营组三角洲相地层,具有被断层复杂化的多油水系统。该区井网密集,储层厚度薄,横向关系复杂,储层与非储层速度差异小,后续开发方案调整所要求的精细储层描述难度较大。本文展示了以精细储层描述为目标的地震约束反演实例,利用多井地震、地质标定,结合三维构造解释成果,建立三维多断层精细约束地质模型。在此基础上,采用随机地震反演方法,进行了波阻抗、岩性等多种地层参数反演,较好地解决了该区薄互层储层精细描述问题,证实了精细储层反演技术在油田开发中的重要作用。     

面向油藏开发地质问题的精细储层反演研究

垦71开发区块的主要产油层系为新近系的馆陶组河流相和古近系的东营组三角洲相地层,具有被断层复杂化的多油水系统。该区井网密集,储层厚度薄,横向关系复杂,储层与非储层速度差异小,后续开发方案调整所要求的精细储层描述难度较大。本文展示了以精细储层描述为目标的地震约束反演实例,利用多井地震、地质标定,结合三维构造解释成果,建立三维多断层精细约束地质模型。在此基础上,采用随机地震反演方法,进行了波阻抗、岩性等多种地层参数反演,较好地解决了该区薄互层储层精细描述问题,证实了精细储层反演技术在油田开发中的重要作用。     

Quasi-shear wave ray tracing by wavefront construction in 3-D, anisotropic media

Wavefront construction (WFC) methods provide robust tools for computing ray theoretical traveltimes and amplitudes for multivalued wavefields. They simulate a wavefront propagating through a model using a mesh that is refined adaptively to ensure accuracy as rays diverge during propagation. However, an implementation for quasi-shear (qS) waves in anisotropic media can be very difficult, since the two qS slowness surfaces and wavefronts often intersect at shear-wave singularities. This complicates the task of creating the initial wavefront meshes, as a particular wavefront will be the faster qS-wave in some directions, but slower in others. Analogous problems arise during interpolation as the wavefront propagates, when an existing mesh cell that crosses a singularity on the wavefront is subdivided. Particle motion vectors provide the key information for correctly generating and interpolating wavefront meshes, as they will normally change slowly along a wavefront. Our implementation tests particle motion vectors to ensure correct initialization and propagation of the mesh for the chosen wave type and to confirm that the vectors change gradually along the wavefront. With this approach, the method provides a robust and efficient algorithm for modeling shear-wave propagation in a 3-D, anisotropic medium. We have successfully tested the qS-wave WFC in transversely isotropic models that include line singularities and kiss singularities. Results from a VTI model with a strong vertical gradient in velocity also show the accuracy of the implementation. In addition, we demonstrate that the WFC method can model a wavefront with a triplication caused by intrinsic anisotropy and that its multivalued traveltimes are mapped accurately. Finally, qS-wave synthetic seismograms are validated against an independent, full-waveform solution.     

A note on dynamic ray tracing in ray-centered coordinates in anisotropic inhomogeneous media

Dynamic ray tracing plays an important role in paraxial ray methods. In this paper, dynamic ray tracing systems for inhomogeneous anisotropic media, consisting of four linear ordinary differential equations of the first order along the reference ray, are studied. The main attention is devoted to systems expressed in a particularly simple choice of ray-centered coordinates, here referred to as the standard ray-centered coordinates, and in wavefront orthonormal coordinates. These two systems, known from the literature, were derived independently and were given in different forms. In this paper it is proved that both systems are fully equivalent. Consequently, the dynamic ray tracing system, consisting of four equations in wavefront orthonormal coordinates, can also be used if we work in ray-centered coordinates, and vice versa. vcerveny@seis.karlov.mff.cuni.cz     

射线法模拟分析井间地震观测的波场特征

按照井间地震的观测系统,用改进的突变点加插值射线追踪方法,追踪每炮每道的射线路径,计算几种主要类型的波沿射线路径的波至时间和射线振幅,制作井间地震多炮多道水平分量和垂直分量的合成记录.并将合成记录选排为井间共炮点道集、共接收点道集、共偏移距道集和共中心深度点道集,系统地分析了不同道集内几种主要类型的地震波的传播特征.对野外观测的实际井间地震记录进行了模拟,从复杂的井间地震记录中,识别出井间地震实际观测到的不同类型的波场,为随后的井间地震资料处理和应用提供了依据.