Perturbation‐based moveout approximations in anisotropic media

The moveout approximations play an important role in seismic data processing. The standard hyperbolic moveout approximation is based on an elliptical background model with two velocities: vertical and normal moveout. We propose a new set of moveout approximations based on a perturbation series in terms of anellipticity parameters using the alternative elliptical background model defined by vertical and horizontal velocities. We start with a transversely isotropic medium with a vertical symmetry axis. Then, we extend this approach to a homogeneous orthorhombic medium. To define the perturbation coefficients for a new background, we solve the eikonal equation with horizontal velocities in transversely isotropic medium with a vertical symmetry axis and orthorhombic media. To stabilise the perturbation series and improve the accuracy, the Shanks transform is applied for all the cases. We select different parameterisations for both velocities and anellipticity parameters for an orthorhombic model. From the comparison in traveltime error, the new moveout approximations result in better accuracy comparing with the standard perturbation‐based methods and other approximations.     

Kinematic parameters of pure‐ and converted‐mode waves in elastic orthorhombic media

I derive the kinematic properties of single‐mode P, S1, and S2 waves as well as converted PS1, PS2, and S1S2 waves in elastic orthorhombic media including vertical velocity, two normal moveout velocities defined in vertical symmetry planes, and three anelliptic parameters (two of them are defined in vertical symmetry plane and one parameter is the cross‐term one). I show that the azimuthal dependence of normal moveout velocity and anellipticity is different in phase and group domains. The effects on‐vertical‐axis singularity and on‐vertical‐axis triplication are considered for pure‐mode S1 and S2 waves and converted‐mode S1S2 waves. The conditions and properties of on‐vertical‐axis triplication are defined in terms of kinematic parameters. The results are illustrated in four homogeneous orthorhombic models and one multilayered orthorhombic model with no variation in azimuthal orientation for all the layers.     

Post‐stack diffraction imaging in vertical transverse isotropy media using non‐hyperbolic moveout approximations

Diffractions carry valuable information about local discontinuities and small‐scale objects in the subsurface. They are still not commonly used in the process of geological interpretation. Many diffraction imaging techniques have been developed and applied for isotropic media, whereas relatively few techniques have been developed for anisotropic media. Ignoring anisotropy can result in low‐resolution images with wrongly positioned or spurious diffractors. In this article, we suggest taking anisotropy into account in two‐dimensional post‐stack domain by considering P‐wave non‐hyperbolic diffraction traveltime approximations for vertical transverse isotropy media, previously developed for reflection seismology. The accuracy of the final images is directly connected to the accuracy of the diffraction traveltime approximations. We quantified the accuracy of six different approximations, including hyperbolic moveout approximation, by the application of a post‐stack diffraction imaging technique on two‐dimensional synthetic data examples.     

水平界面上P-SV转换波转换点的精确解

转换波共转换点的叠加和道集选取都需要准确地计算转换点的位置.Tessmer和Behle、Taylor分别给出了水平单层介质中P-SV转换波转换点坐标的解析解,由于其表达式的复杂性,在应用中几乎不被采用.在本文中,运用Snell定律重新建立了在水平反射界面上反射的P-SV转换波的转换点坐标的四次方程,并严格地推导出与纵波速度、横波速度、炮检距和反射深度有关的转换点坐标的解析解,确定了惟一的解析表达式.将这一结果应用于P-SV转换波的速度分析和叠加处理中.简化的公式有较好的应用价值.     

VTI介质转换波动态聚焦束偏移成像方法研究

采用多分量地震数据进行处理,可以得到高质量的转换波成像结果,进而为流体识别和裂缝预测等提供有效支持.本文在前人研究的基础上,通过求解各向异性介质qPSV转换波运动学和动力学射线追踪方程,获得转换波中心射线走时、路径及动力学信息;同时,通过修改动态聚焦型传播算子,构建出各向异性介质动态聚焦束,提出了一种各向异性介质转换波...     

横向各向同性地层井孔伪瑞利波、弯曲波、螺旋波的低频极限速度

在多极源声测井中, 低频弯曲波或螺旋波被广泛应用于测量地层的横波速度. 前人的研究已证明在各向同性地层中井孔伪瑞利波、弯曲波及螺旋波的低频极限速度都等于地层横波速度. 大量的数值计算结果似乎表明此结论在横向各向同性（TI）地层情况下也能成立, 但缺乏理论证明. 本文在井孔平行于TI弹性地层对称轴的模型下, 考察了井内声源激发的波在流-固边界上的反射和透射情况, 阐述了非泄漏模式导波产生的必要条件并讨论了其速度上限值. 我们发现在各向异性参数满足一定条件的TI地层中,导波速度的低频极限值小于沿井孔方向上的横波速度. 通过对井孔导波的频散分析以及对时域全波列的数值模拟, 论证了在这类地层中进行多极源声测井时,不可能根据弯曲波或螺旋波的到达时间准确获取地层的横波速度值.     

转换波AVO反演速度比各横波反射系数

Based on the empirical Gardner equation describing the relationship between density and compressional wave velocity, the converted wave reflection coefficient extrema attributes for AVO analysis are proposed and the relations between the extrema position and amplitude, average velocity ratio across the interface, and shear wave reflection coefficient are derived. The extrema position is a monotonically decreasing function of average velocity ratio, and the extrema amplitude is a function of average velocity ratio and shear wave reflection coefficient. For theoretical models, the average velocity ratio and shear wave reflection coefficient are inverted from the extrema position and amplitude obtained from fitting a power function to converted wave AVO curves. Shear wave reflection coefficient sections have clearer physical meaning than conventional converted wave stacked sections and establish the theoretical foundation for geological structural interpretation and event correlation. ＂The method of inverting average velocity ratio and shear wave reflection coefficient from the extrema position and amplitude obtained from fitting a power function is applied to real CCP gathers. The inverted average velocity ratios are consistent with those computed from compressional and shear wave well logs.     

Observation of shear‐wave splitting in the multicomponent node data from Atlantis field,Gulf of Mexico

In 2005, a multicomponent ocean bottom node data set was collected by BP and BHP Billiton in the Atlantis field in the Gulf of Mexico. Our results are based on data from a few sparse nodes with millions of shots that were analysed as common receiver azimuthal gathers. A first‐order look at P‐wave arrivals on a common receiver gather at a constant offset reveals variation of P‐wave arrival time as a function of azimuth indicating the presence of azimuthal anisotropy at the top few layers. This prompted us to investigate shear arrivals on the horizontal component data. After preliminary processing, including a static correction, the data were optimally rotated to radial (R) and transverse (T) components. The R component shows azimuthal variation of traveltime indicating variation of velocity with azimuth; the corresponding T component shows azimuthal variation of amplitude and phase (polarity reversal). The observed shear‐wave (S‐wave) splitting, previously observed azimuthal P‐wave velocity variation and azimuthal P‐wave amplitude variation, all indicate the occurrence of anisotropy in the shallow (just below the seafloor) subsea sediment in the area. From the radial component azimuthal gather, we analysed the PP‐ and PS‐wave amplitude variation for the first few layers and determined corresponding anisotropy parameter and V_P/V_S values. Since fracture at this depth is not likely to occur, we attribute the observed azimuthal anisotropy to the presence of microcracks and grain boundary orientation due to stress. The evidence of anisotropy is ubiquitous in this data set and thus it argues strongly in favour of considering anisotropy in depth imaging for obtaining realistic subsurface images, at the least.     

多层倾斜层状介质P-SV转换波转换点计算的最佳角度搜索法

转换波转换点的位置对转换波道集的抽取和叠加都非常重要.目前对于单一倾斜反射界面或水平层状介质模型,已经能较好地计算转换波转换点的位置.本文针对多层倾斜层状介质,提出最佳角度搜索法求取P-SV转换波转换点的位置,此方法是通过搜索最符合Snell定律的入射角和反射角来确定最佳的P-SV转换波转换点的位置.在搜索最佳转换点位置的同时,该方法也能给出最佳的P-SV波的传播路径.本文在单层和双层倾斜层状介质模型中测试了该方法,获得了精确的P-SV转换波转换点的位置和P-SV波的传播路径.     

南海东北部及其邻近地区的Pn波速度结构与各向异性

利用中国地震台网和ISC台站1980～2004年的地震数据,反演了南海东北部及其邻近地区的Pn波速度结构和各向异性.上地幔顶部的速度变化揭示出区域地质构造的深部特征：华南地区速度较高并且变化平缓,具有构造稳定地区的岩石层地幔特征;华南沿海尤其是滨海断裂带附近出现低速异常,表明该断裂可能穿过壳幔边界深达上地幔顶部.南海北部至台湾海峡较高的速度与华南地区类似,反映出大陆边缘和陆架地区的岩石层地幔性质;西沙海槽附近较高的速度不仅反映了华南大陆向南的延伸,而且与海槽裂谷拉张引起的地幔上拱有关,整个南海北部没有发现大规模地幔热流的活动痕迹.相比之下,南海东部次海盆的上地幔顶部存在明显的低速异常,对应于海底扩张中心的地幔上涌区,表明岩石层地幔强烈减薄甚至缺失;台湾东部－吕宋－菲律宾北部的低速异常与地震、火山活动以及岩浆作用紧密相关,揭示了西太平洋岛弧俯冲带的活动特征;南海东北部的洋－陆边界清晰,南海东部和菲律宾海西部较高的速度代表了海洋岩石层地幔的性质.Pn波各向异性反映出区域性构造应力状态及岩石层地幔的变形痕迹：华南地区的各向异性较小,说明这一构造稳定地区的岩石层地幔变形程度较弱;南海北部的快波方向与地壳浅表层构造的伸展方向一致,主要反映了中、新生代以来的大陆边缘张裂和剪切作用对岩石层地幔结构的影响;琉球－台湾－吕宋岛弧两侧各向异性十分强烈,平行于海沟的快波方向表明菲律宾海板块和欧亚大陆的相互作用导致俯冲板块前缘的岩石层地幔强烈变形;台湾东南海域快波方向的变化可能与欧亚大陆和菲律宾海板块俯冲机制的转换以及岩石层被撕裂有关.     

Frequency-dependence of velocity and attenuation of elastic waves in granite under uniaxial compression

Velocity as well as attenuation factorQ ^–1 ofP-wave in a dry granitic rock sample under uniaxial compressions were measured in the range of frequency between 100 kHz and 710 kHz by using the pulse transmission technique. Above the stress of 0.5 _f , where _f is the fracture stress, theP-wave velocity decreases with increasing axial stress, whereasQ ^–1 increases. Particularly, the change ofQ ^–1 is greater for high frequency than for low frequency. At a given stress level, the higher the frequency, the higher theP-wave velocity and the largerQ ^–1. This result means that the velocity decrease with increasing stress is smaller for higher frequency. Because of this frequency-dependence of velocity decrease, theP-wave in the rock under dilatant state shows dispersion. The body wave dispersion is more remarkable at higher stress, and is not found in a homogeneous material with no cracks. Thus the disperison is attributed to the generation of cracks. When the frequency-dependence ofQ ^–1 is approximated asf ⁿ in the present frequency range, the exponentn takes a value from 0.63 to 0.77.     

Group velocity distribution of Rayleigh waves and crustal and upper mantle velocity structure of the Chinese mainland and its vicinity

Introduction A lot of results have been achieved on the study of crust and upper mantle structure in Chinese mainland by the surface wave dispersion. The seismometer can be classified into three types based on the development from the analog to digital instrument. FENG, et al (1981) measured the dispersion curves of group and phase velocity along 28 paths, with 19 earthquakes recorded by 15 base analog seismic stations, subdivided the Chinese mainland into 5 blocks and obtained the average c…     

Surface‐wave inversion for a P‐velocity profile with a constant depth gradient of the squared slowness

Surface waves are often used to estimate a near‐surface shear‐velocity profile. The inverse problem is solved for the locally one‐dimensional problem of a set of homogeneous horizontal elastic layers. The result is a set of shear velocities, one for each layer. To obtain a P‐wave velocity profile, the P‐guided waves should be included in the inversion scheme. As an alternative to a multi‐layered model, we consider a simple smooth acoustic constant‐density velocity model, which has a negative constant vertical depth gradient of the squared P‐wave slowness and is bounded by a free surface at the top and a homogeneous half‐space at the bottom. The exact solution involves Airy functions and provides an analytical expression for the dispersion equation. If the ratio is sufficiently small, the dispersion curves can be picked from the seismic data and inverted for the continuous P‐wave velocity profile. The potential advantages of our model are its low computational cost and the fact that the result can serve as a smooth starting model for full‐waveform inversion. For the latter, a smooth initial model is often preferred over a rough one. We test the inversion approach on synthetic elastic data computed for a single‐layer P‐wave model and on field data, both with a small ratio. We find that a single‐layer model can recover either the shallow or deeper part of the profile but not both, when compared with the result of a multi‐layer inversion that we use as a reference. An extension of our analytic model to two layers above a homogeneous half‐space, each with a constant vertical gradient of the squared P‐wave slowness and connected in a continuous manner, improves the fit of the picked dispersion curves. The resulting profile resembles a smooth approximation of the multi‐layered one but contains, of course, less detail. As it turns out, our method does not degrade as gracefully as, for instance, diving‐wave tomography, and we can only hope to fit a subset of the dispersion curves. Therefore, the applicability of the method is limited to cases where the ratio is small and the profile is sufficiently simple. A further extension of the two‐layer model to more layers, each with a constant depth gradient of the squared slowness, might improve the fit of the modal structure but at an increased cost.     

青藏高原东南部地区瑞雷波相速度层析成像

本研究收集了"中国地震科学探测台阵-南北地震带南段"项目325个流动宽频带台站于2011年8月至2012年9月记录的远震垂直向资料,利用双台法测得了3594条独立路径上的瑞雷波相速度频散曲线,反演得到了青藏高原东南部地区周期10~60 s瑞雷波的相速度分布图像.空间分辨尺度图表明,在台站覆盖范围内的绝大部分地区横向分辨率达到50 km.2D相速度分布图显示,青藏高原东南部地区地壳上地幔S波速度结构存在较明显的横向非均匀性.短周期(如10 s)的相速度分布主要受地表沉积层厚度的影响.绝大多数地震发生在周期15 s相速度图上的低速区或高低速的陡变梯度带附近,充分说明该区的强震活动与中上地壳速度结构的变化有直接关系.中等周期(如20~30 s)的相速度分布主要与中下地壳速度结构、地壳厚度密切相关,小江断裂、松潘-甘孜块体呈现最显著的低速,可能暗示这两处的中、下地壳存在低速层.较长周期(如40~60 s)的相速度分布与上地幔顶部热状态和构造活动(如岩浆作用)有关.滇西南地区表现为大范围的显著低速,可能暗示滇西南地区上地幔顶部物质存在部分熔融.不同构造块体下方的频散曲线,具有不同的相速度特征.腾冲火山下方的频散曲线在10~60 s一直为较低的速度,尤其是到40 s以后,相速度随周期的变大增速明显放缓,至60 s比其他任何块体速度都低,暗示腾冲火山区下方的低速至少来自上地幔顶部(约100 km).     

郯庐断裂带及其邻区上地幔顶部Pn波速度与各向异性层析成像

本研究拾取了中国数字测震台网固定台站记录的2008-2016年2级以上地震事件中的27233条高质量Pn到时资料,反演得到了郯庐断裂带及其邻区上地幔顶部Pn波速度和各向异性结构模型.结果显示,研究区上地幔顶部Pn波速度结构存在强烈的横向不均匀性,速度异常形态与区域地质构造较为吻合.太行山造山带、鲁西隆起、大别造山带、苏鲁褶皱带、胶辽隆起和华北盆地南端等隆起区表现为低波速异常,而黄海北、南部盆地、渤海湾和华北盆地北部等凹陷区均为高波速异常.壳内强震主要发生在Pn低波速异常和高低波速异常的横向过渡地带,说明强震的发生与上地幔结构的横向变化之间存在有一定关联.郯庐断裂带两侧Pn波速度以郯城地震为界其东北侧和西南侧分别分布有与断裂带近平行的低波速异常条带,而西北侧和东南侧分别分布有高波速异常条带,各向异性快波方向近乎沿断裂带走向,可能由于上地幔热物质沿郯庐断裂带上涌形成低速异常后断裂带发生左旋平移运动所致.华北盆地内上地幔顶部Pn波速度结构和各向异性的明显变化,反映华北克拉通破坏过程中经历了地幔热物质上涌、莫霍面隆升以及岩石圈拆沉等复杂构造变形.

     

Multi‐block finite‐difference method for 3D elastodynamic simulations in anisotropic subhorizontally layered media

Prediction of elastic full wavefields is required for reverse time migration, full waveform inversion, borehole seismology, seismic modelling, etc. We propose a novel algorithm to solve the Navier wave equation, which is based on multi‐block methodology for high‐order finite‐difference schemes on curvilinear grids. In the current implementation, the blocks are subhorizontal layers. Smooth anisotropic heterogeneous media in each layer can have strong discontinuities at the interfaces. A curvilinear adaptive hexahedral grid in blocks is generated by mapping the original 3D physical domain onto a parametric cube with horizontal layers and interfaces. These interfaces correspond to the main curvilinear physical contrast interfaces of a subhorizontally layered formation. The top boundary of the parametric cube handles the land surface with smooth topography. Free‐surface and solid–solid transmission boundary conditions at interfaces are approximated with the second‐order accuracy. Smooth media in the layers are approximated up to sixth‐order spatial schemes. All expected properties of the developed algorithm are demonstrated in numerical tests using corresponding parallel message passing interface code.     

郯庐断裂带中南段及邻区Pn波速度结构与各向异性

郯庐断裂带是一条纵贯我国大陆东部NNE走向的巨型深断裂,其中南段及邻区(115°E-122°E,29°N-38°N)跨越了华北断块区、扬子断块区和华南褶皱系三大一级构造单元,由于其重要性和复杂性,长期以来一直是地学家们研究的热点.本文从国际地震中心(ISC)、中国地震台网及区域地震台网的地震观测报告中精心挑选出6381个Pn震相数据,用Pn波时间项层析成像法反演得到了郯庐断裂带中南段及邻区上地幔顶部Pn波速度结构和各向异性.结果显示,研究区上地幔顶部具有显著的横向非均匀性,相对于7.95 km·s-1的平均速度而言,Pn波速度值在7.68~8.24 km·s-1范围内变化.Pn波速度分布在郯庐断裂带中段和南段具有分段性:沿中段及周边存在一NE向低速异常带,低速可能是由于岩石圈的减薄和软流圈的高温物质沿郯庐带上涌导致;沿南段表现为一NNE向弱高波速异常带,作为高低速的边界带清晰地勾勒出了华北与扬子这两个不同块体,该边界在江苏域向华北地块NW方向凹进.Pn波速度各向异性的强弱与速度分布存在一定的相关性.总体上,如鲁西隆起及以南等低速区、茅山断裂附近的高低速过渡带,其速度各向异性较为强烈;而在具有高速异常的苏北盆地、合肥盆地等稳定区域下方其各向异性较弱.本文通过Pn波震相基本未能探测到郯庐断裂带中段的方位各向异性,推测是上地幔顶部被"冻结"下来的各向异性痕迹被软流圈热物质上涌这一强烈构造运动削弱所导致.南段具有与断裂伸展方向近乎平行的快波速方向.Pn波速度横向变化和强震活动存在一定关联.强震主要发生在Pn波低速异常区或高低速过渡带上.郯城8.5级地震震中位于中段和南段高低速过渡带,该区域也是速度横向变化最大的地方,最容易集中应力和产生应力差.     

正交各向异性介质中qP波入射的二阶近似反射系数与透射系数

地震波在各向异性介质中以一个准P波（qP）和两个准S波（qS1和qS2）的形式传播.研究三种波的相速度、群速度以及偏振方向等传播性质能够为各向异性介质中的正反演问题提供有效支撑.具有比横向各向同性（TI）介质更一般对称性的正交各向异性介质通常需要9个独立参数对其进行描述,这使得对传播特征的计算更为复杂.当两个准S波速度相近时具有耦合性,从而令慢度的计算产生奇异性.因此,奇异点（慢度面的鞍点和交叉点）附近的反射与透射（R/T）系数的求解不稳定,会导致波场振幅不准确.本文首次通过结合耦合S波射线理论和基于迭代的各向异性相速度与偏振矢量的高阶近似解,得到了适用于正交各向异性介质以qP波入射所产生的二阶R/T系数的计算方法.与基于一阶近似的结果相比,基于二阶近似的方法提高了qP波R/T系数的精度,能得到一阶耦合近似无法表达的准确的qP-qS转换波的R/T系数解,且方法适用于较强的各向异性介质.

     

Influence of crack distribution of rocks on P‐wave velocity anisotropy – a laboratory and field scale study‡

The purpose of this paper is the comparison of P‐wave velocity and velocity anisotropy, measured at different scales under laboratory and field conditions. A shallow seismic refraction survey with shot/receiver spacing of up to 10 m was carried out on a flat outcrop of lhertzolite in the southern part of the Balmuccia massif. Oriented rock samples were also obtained from the locality. The particular advantage of the laboratory method used is the possibility of measuring velocity in any direction under controlled conditions. Laboratory tests were made on spherical peridotite samples, 50 mm in diameter, by ultrasonic velocity measurements in 132 directions (meridian and parallel networks) under confining stress ranging from atmospheric to 400 MPa. The mean P‐wave velocity of the field and laboratory data differed by between 20–30%. In addition, P‐wave velocity anisotropy of 25% was detected in the field data. Whereas the anisotropy in the laboratory samples in the same orientation as the field surveys was less than 2%. This observed scaling factor is related to the different sampling sizes and the difference in frequencies of applied elastic waves. With an ultrasonic wavelength of 10 mm, laboratory samples represent a continuum. The field velocities and velocity anisotropy reflect the presence of cracks, which the laboratory rock samples do not contain. Three sub‐vertical fracture sets with differing strikes were observed in the field outcrop. Estimates of fracture stiffness from the velocity anisotropy data are consistent with other published values. These results highlight the difficulty of using laboratory velocity estimates to interpret field data.