克希霍夫法VSP多波联合成像

VSP 转换波跟VSP 纵波或常规地面转换波相比,具有较高的分辨率和信噪比,但传统的VSP成像方法只利用了反射P波信息,而把转换波（反射S波、透射S波）以及透射P波当作影响成像质量的噪音.本文给出了一种VSP共炮点道集多分量地震资料克希霍夫法偏移成像的方法.本方法充分利用了多波（反射P波、反射S波、透射P波、透射S波）信息,根据转换点处四种波同时起跳,能量叠加最大的原理,从接收点分别用向绕射点延拓它们的能量,并将其叠加起来,求得的和最大的一点即反射点.通过模型试算和实际资料处理表明,此法成像精度高,信噪比高,且有利于改善剖面的频率特性.     

THE USE OF FORWARD- AND BACK-SCATTERED P-, S- AND CONVERTED WAVES IN CROSS-BOREHOLE IMAGING1

The principles of imaging, for example that of prestack migration, can be applied to cross-borehole seismic geometry just as they can to surface seismic configurations. However, when using actual cross-borehole data, a number of difficulties arise that are rarely or never encountered in imaging surface seismic data: discontinuities may reflect or diffract incident seismic waves in any direction. If a discontinuity lies between the lines of sources and receivers, forward-scattered, or interwell, events may be recorded. If a discontinuity lies outside the interwell region, back-scattered, or extra-well, events may be recorded. Many angles of incidence are possible, and all possible reflected modes (P–P, P–S, S–P and S–S) are present, frequently in nearly equal proportions. The planes of the reflectors dip from 0 to ±90°. In order to deal with these complexities we first separate propagation modes at the receiver borehole using both polarization and velocity. Next we compensate for phase distortion due to dispersion. Finally, and most importantly, we migrate or image the data in cross-borehole common-source gathers. To do this, a finite-difference solution to the 2D scalar wave equation, using reverse time, for an arbitrary distribution of velocities, is used to project the separated, reflected-diffracted wavefield back into the medium. There are four reflection modes (P–P, P–S, S–P and S–S), so we can apply four different imaging conditions. The zones outside the boreholes as well as inside the boreholes can be imaged with these conditions. These operations are repeated for each common-source gather: each common-source gather generates four partial images in each image space. This multiplicity of partial images can be stacked in various combinations to yield a final image of the subsurface. Our experiments using solid (not fluid) physical models indicate that when these procedures are correctly applied, high quality cross-borehole images can be obtained. These images appear with great clarity even though some of the weak diffractions causing diffraction images may be almost totally obscured by other high-amplitude events on the raw data.     

Wide-angle incidence and P-wave transmission

Polarity reversals may occur to transmitted P waves if the incidence angle is greater than the critical incidence angle. We analyze the characteristics of reflection and transmission coefficients under the condition of wide incidence angle based on Zoeppritz equations. We find that for specific conditions, as the incidence angle increases, the characteristic curve of the transmitted P-wave coefficient enters the third quadrant from the first quadrant through the origin, which produces a transition in the transmitted P wave and the corresponding coefficient experiences polarity reversal. We derive the incidence angle when the transmitted P-wave coefficient is zero and verify that it equals zero by using finite-difference forward modeling for a single-interface model. We replace the water in the model reservoir by gas and see that the reservoir P-wave velocity and density decrease dramatically. By analyzing the synthetic seismogram of the transmitted P wave in the single-interface model, we show that the gas-saturated reservoir is responsible for polarity reversal.     

共接收点倾斜叠加波动方程偏移

共接收点倾斜叠加波动方程偏移,本质上是一种叠前偏移方法.每给定一个斜率P,对经过叠前（动校正前）常规处理的地震记录中的各共接收点道集,沿直线t=τ+px进行倾斜叠加,就形成一个共接收点倾斜叠加剖面.对之进行波动方程偏移,该偏移剖面将代表地下真实构造.对一系列的p,我们可以得到一系列这样的偏移剖面.对它们作共接收点叠加,偏移叠加剖面的信噪比将超过水平叠加剖面.本文导出了在均匀、水平层状及非均匀介质条件下的共接收点倾斜叠加波动方程偏移算法.     

Extraction of amplitude-preserving angle gathers based on vector wavefield reverse-time migration

Angle-domain common-image gathers (ADCIGs) transformed from the shotdomain common-offset gathers are input to migration velocity analysis (MVA) and prestack inversion. ADCIGs are non-illusion prestack inversion gathers, and thus, accurate. We studied the extraction of elastic-wave ADCIGs based on amplitude-preserving elastic-wave reversetime migration for calculating the incidence angle of P-and S-waves at each image point and for different source locations. The P-and S-waves share the same incident angle, namely the incident angle of the source P-waves. The angle of incidence of the source P-wavefield was the difference between the source P-wave propagation angle and the reflector dips. The propagation angle of the source P-waves was obtained from the polarization vector of the decomposed P-waves. The reflectors’ normal direction angle was obtained using the complex wavenumber of the stacked reverse-time migration (RTM) images. The ADCIGs of P-and S-waves were obtained by rearranging the common-shot migration gathers based on the incident angle. We used a horizontally layered model, the graben medium model, and part of the Marmousi-II elastic model and field data to test the proposed algorithm. The results suggested that the proposed method can efficiently extract the P-and S-wave ADCIGs of the elastic-wave reverse-time migration, the P-and S-wave incident angle, and the angle-gather amplitude fidelity, and improve the MVA and prestack inversion.     

接收函数方法研究进展

远震P波波形数据中包含了大量的地震台站下方地壳和上地慢速度间断面所产生的PS转换波及其多次反射波的信息，由此提取的接收函数是了解地壳上地慢速度精细结构的重要步骤之一。本介绍了目前在提高接收函数结果的稳定性和精度方面的研究进展，包括石油勘探中一些成熟的地震反射处理方法逐步渗透到接收函数的研究领域。这些处理方法以前所未有的分辨率展示了地壳上地幔结构的横向非均匀性。     

Teleseismic body waves extracted from ambient noise cross correlation between F-net and ChinArray phase II

The vertical-vertical noise cross-correlation functions (NCFs) between two seismic arrays, the Japan F-net and ChinArray phase II, are calculated using continuous recordings during 2013–2016. After array interferometry to obtain bin stacked NCFs, clear body waves are retrieved at different period bands. Teleseismic direct P waves for distance 15–40 degrees are observed between short period 3–10 ?s while core reflected PcP/ScS waves are more obvious for longer period 30–60s. The signal-to-noise-ratio (SNR) of the short period P waves reaches its highest point with bin widths around 20 ?km while SNRs of PcP and ScS increase slowly with bin width. All those body waves demonstrate clear directivity with strong signals traveling from the east. The time-lapse SNR variations for the PcP and ScS show correlation with the occurrence of major earthquakes, while the P-wave SNR demonstrates seasonal variations with additional contribution from major earthquakes. The present results suggest teleseismic body waves can be retrieved through bin stacking, though further processing is still necessary to obtain finer waveforms such as P wave triplications.     

2D and 3D amplitude-preserving elastic reverse time migration based on the vector-decomposed P- and S-wave records

The elastic reverse time migration approach based on the vector-wavefield decomposition generally uses the scalar product imaging condition to image the multicomponent seismic data. However, the resulting images contain the crosstalk artefacts and the polarity reversal problems, which are caused by the nonphysical wave modes and the angle-dependent reduction of image amplitudes, respectively. To overcome these two problems, we develop an amplitude-preserving elastic reverse time migration approach based on the vector-decomposed P- and S-wave seismic records. This approach includes two key points. The first is that we employ the vector-decomposed P- and S-wave multicomponent records to independently reconstruct the PP and PS reflection images to mitigate the crosstalk artefacts. The second is that we propose two schemes in addressing the issue of polarity reversal problem in the conventional PP image. One solution is to adopt the angle-dependent equation. Another one is to reconstruct an amplitude-preserving PP image with a separated scalar P-wave particle velocity, which has a clear physical meaning. Numerical examples using two-dimensional and three-dimensional models demonstrate that the proposed elastic reverse time migration approach can provide the images with better amplitude-preserving performance and fewer crosstalk artefacts, compared with the conventional elastic reverse time migration approach based on the scalar product imaging condition.     

Directional‐oriented wavefield imaging: a new wave‐based subsurface illumination imaging condition for reverse time migration

The key objective of an imaging algorithm is to produce accurate and high‐resolution images of the subsurface geology. However, significant wavefield distortions occur due to wave propagation through complex structures and irregular acquisition geometries causing uneven wavefield illumination at the target. Therefore, conventional imaging conditions are unable to correctly compensate for variable illumination effects. We propose a generalised wave‐based imaging condition, which incorporates a weighting function based on energy illumination at each subsurface reflection and azimuth angles. Our proposed imaging kernel, named as the directional‐oriented wavefield imaging, compensates for illumination effects produced by possible surface obstructions during acquisition, sparse geometries employed in the field, and complex velocity models. An integral part of the directional‐oriented wavefield imaging condition is a methodology for applying down‐going/up‐going wavefield decomposition to both source and receiver extrapolated wavefields. This type of wavefield decomposition eliminates low‐frequency artefacts and scattering noise caused by the two‐way wave equation and can facilitate the robust estimation for energy fluxes of wavefields required for the seismic illumination analysis. Then, based on the estimation of the respective wavefield propagation vectors and associated directions, we evaluate the illumination energy for each subsurface location as a function of image depth point and subsurface azimuth and reflection angles. Thus, the final directional‐oriented wavefield imaging kernel is a cross‐correlation of the decomposed source and receiver wavefields weighted by the illuminated energy estimated at each depth location. The application of the directional‐oriented wavefield imaging condition can be employed during the generation of both depth‐stacked images and azimuth–reflection angle‐domain common image gathers. Numerical examples using synthetic and real data demonstrate that the new imaging condition can properly image complex wave paths and produce high‐fidelity depth sections.     

Features of different types of seismic events in China’s Capital Region

Seismic records produced by different seismic sources vary. In this study, we compared the waveform records and time-frequency characteristics of tectonic earthquakes, artificial explosions, and mine collapses in China’s Capital Region. The results show that tectonic earthquakes are characterized by stronger S-wave energy than P-wave energy, obvious high-frequency components, and wide frequency bands of P and S waves. Artificial explosions are characterized by greater P-wave amplitude than S-wave amplitude and near-station surface wave development. Mine collapses are characterized by lower overall frequency, more obvious surface waves, and longer duration. We extracted quantitative discriminants based on the analysis of different event records, with 31 feature values in 7 categories (P/S maximum amplitude ratio, high/low frequency energy ratio, P/S spectral ratio, corner frequency, duration, the second-order moment of spectrum, and energy strongest point). A comparison of the ability of these feature values to recognize distinct events showed that the 6–17 Hz P/S spectral ratio was able to completely distinguish artificial explosions from the other two types of events. The S-wave corner frequency performed relatively well in identifying all three types of events, with an accuracy of over 90%. Additionally, a support vector machine was used to comprehensively distinguish multiple features, with an accuracy for all three types of events reaching up to 100%.     

Reflection point sideslip and smear in imaging below dipping anisotropic media

Dipping anisotropic clastic strata are ubiquitous in fold and thrust belts. Geological structures below these strata will be mispositioned laterally and vertically on seismic images if we do not properly correct for seismic anisotropy during migration. The magnitude of this lateral mispositioning of a target structure varies with source‐receiver offset, so reflection points will be smeared in the final stacked image. Raytracing demonstrates the lateral‐position and smear phenomena when imaging structures below tilted transversely isotropic media. Analysis of the raytracing results predicts the quantity of lateral‐position error and reflection‐point smear on a seismic image. We created numerical‐model seismic data to show reflection‐point smear on synthetic seismic images and to evaluate the accuracy of the predictions from raytracing.     

MULTICOMPONENT COMMON-RECEIVER GATHER MIGRATION OF SINGLE-LEVEL WALK-AWAY SEISMIC PROFILES1

Seismic data are usually separated into P-waves and S-waves before being put through a scalar (acoustic) migration. The relationship between polarization and moveout is exploited to design filters that extract the desired wavetype. While these filters can always be applied to shot records, they can only be applied to a triaxial common-receiver gather in special cases since the moveout of scattered energy on the receiver gather relates to path differences between the surface shots and the scatterer while the polarization is determined by the path from scatterer to downhole geophone. Without the ability to separate wavefields before migration, a ‘vector scalar’ or an elastic migration becomes a necessity. Here the propagation of the elastic wavefield for a given mode (e.g. P-S) is approximated by two scalar (acoustic) propagation steps in a ‘vector scalar’ migration. ‘Vector’ in that multicomponent data is migrated and 'scalar’ in that each propagation step is based on a scalar wave equation for the appropriate mode. It is assumed that interaction between the wavefields occurs only once in the far-field of both the source and receiver. Extraction of the P, SV and SH wavefields can be achieved within the depth migration (if one assumes isotropy in the neighbourhood of the downhole receiver) by a projection onto the polarization for the desired mode. Since the polarization of scattered energy is only a function of scatterer position and receiver position (and not source position), the projection may be taken outside the migration integral in the special case of the depth migration of a common-receiver gather. The extraction of the desired mode is then performed for each depth migration bin after the separate scalar migration of each receiver gather component. This multicomponent migration of triaxial receiver gathers is conveniently implemented with a hybrid split-step Fourier-excitation-time imaging condition depth migration. The raytracing to get the excitation-time imaging condition also provides the expected polarization for the post-migration projection. The same downward extrapolated wavefield can be used for both the P-P and P-S migrations, providing a flexible and efficient route to the migration of multicomponent data. The technique is illustrated on a synthetic example and a single-level Walk-away Seismic Profile (WSP) from the southern North Sea. The field data produced images showing a P-P reflector below the geophone and localized P-P and P-S scatterers at the level of the geo-phone. These scatterers, which lie outside the zone of specular illumination, are interpreted as faults in the base Zechstein/top Rotliegendes interface.     

Quantitative imaging of complex structures from dense wide-aperture seismic data by multiscale traveltime and waveform inversions: a case study

An integrated multiscale seismic imaging flow is applied to dense onshore wide‐aperture seismic data recorded in a complex geological setting (thrust belt). An initial P‐wave velocity macromodel is first developed by first‐arrival traveltime tomography. This model is used as an initial guess for subsequent full‐waveform tomography, which leads to greatly improved spatial resolution of the P‐wave velocity model. However, the application of full‐waveform tomography to the high‐frequency part of the source bandwidth is difficult, due to the non‐linearity of this kind of method. Moreover, it is computationally expensive at high frequencies since a finite‐difference method is used to model the wave propagation. Hence, full‐waveform tomography was complemented by asymptotic prestack depth migration to process the full‐source bandwidth and develop a sharp image of the short wavelengths. The final traveltime tomography model and two smoothed versions of the final full‐waveform tomography model were used as a macromodel for the prestack depth migration. In this study, wide‐aperture multifold seismic data are used. After specific preprocessing of the data, 16 frequency components ranging from 5.4 Hz to 20 Hz were inverted in cascade by the full‐waveform tomography algorithm. The full‐waveform tomography successfully imaged SW‐dipping structures previously identified as high‐resistivity bodies. The relevance of the full‐waveform tomography models is demonstrated locally by comparison with a coincident vertical seismic profiling (VSP) log available on the profile. The prestack depth‐migrated images, inferred from the traveltime, and the smoothed full‐waveform tomography macromodels are shown to be, on the whole, consistent with the final full‐waveform tomography model. A more detailed analysis, based on common‐image gather computations, and local comparison with the VSP log revealed that the most accurate migrated sections are those obtained from the full‐waveform tomography macromodels. A resolution analysis suggests that the asymptotic prestack depth migration successfully migrated the wide‐aperture components of the data, allowing medium wavelengths in addition to the short wavelengths of the structure to be imaged. The processing flow that we applied to dense wide‐aperture seismic data is shown to provide a promising approach, complementary to more classical seismic reflection data processing, to quantitative imaging of complex geological structures.     

The common focal point and common reflection surface methodologies: subsets or complements?

The common focal point (CFP) method and the common reflection surface (CRS) stack method are compared. The CRS method is a fast, highly automated procedure that provides high S/N ratio simulation of zero‐offset (ZO) images by combining, per image point, the reflection energy of an arc segment that is tangential to the reflector. It uses smooth parametrized two‐way stacking operators, based on a data‐driven triplet of attributes in 2D (eight parameters in 3D). As a spin‐off, the attributes can be used for several applications, such as the determination of the geometrical spreading factor, multiple prediction, and tomographic inversion into a smooth background velocity model. The CFP method aims at decomposing two‐way seismic reflection data into two full‐aperture one‐way propagation operators. By applying an iterative updating procedure in a half‐migrated domain, it provides non‐smooth focusing operators for prestack imaging using only the energy from one focal point at the reflector. The data‐driven operators inhibit all propagation effects of the overburden. The CFP method provides several spin‐offs, amongst which is the CFP matrix related to one focal point, which displays the reflection amplitudes as measured at the surface for each source–receiver pair. The CFP matrix can be used to determine the specular reflection source–receiver pairs and the Fresnel zone at the surface for reflection in one single focal point. Other spin‐offs are the prediction of internal multiples, the determination of reflectivity effects, velocity‐independent redatuming and tomographic inversion to obtain a velocity–depth model. The CFP method is less fast and less automated than the CRS method. From a pointwise comparison of features it is concluded that one method is not a subset of the other, but that both methods can be regarded as being to some extent complementary.     

Common-Reflection-Surface Stack for Converted Waves

The finite-offset (FO) common-reflection-surface (CRS) stack has been shown to be able to handle not only P-P or S-S but also arbitrarily converted reflections. It can provide different stack sections such as common-offset (CO), common-midpoint (CMP) and common-shot (CS) sections with significantly increased signal-to-noise ratio from the multi-coverage pre-stack seismic data in a data-driven way. It is our purpose in this paper to demonstrate the performance of the FO CRS stack on data involving converted waves in inhomogeneous layered media. In order to do this we apply the FO CRS stack for common-offset to a synthetic seismic data set involving P-P as well as P-S converted primary reflections. We show that the FO CRS stack yields convincing improvement of the image quality in the presence of noisy data and successfully extracts kinematic wavefield attributes useful for further analyses. The extracted emergence angle information is used to achieve a complete separation of the wavefield into its P-P and P-S wave components, given the FO CRS stacked horizontal and vertical component sections.     

华北地区中部地震精定位与构造应力场研究

选取华北地区中部2008年1月至2012年12月741次地震事件的波形文件,读取P波初动极性以及P、S波震相到时.利用波形互相关技术,精确计算地震对到时差,结合读取出的震相到时数据,采用双差定位法,对这些地震进行了重新定位,获得了468个地震的精确重定位结果.相比于初始结果,精定位结果在平面分布上更为集中,沿断裂带分布特征更明显;深度也更为合理.在新河断裂附近,存在明显的地震集中条带,整体走向约为北偏东35°,通过剖面分析,发现该断裂倾角很高.依据地震精定位获取P波射线参数,利用P波初动极性,采用改进的格点尝试法,计算了区内单个地震震源机制解及小震综合断层面解,并结合已有的应力数据,综合分析了区内构造应力环境.结果表明,华北地区中部现今构造应力场保持稳定,为最大主应力轴北东东一南西西向,最小主应力轴北北西-南南东向的走滑型应力状态.     

接收函数方法及研究进展

远震P波波形数据中包含了大量在台站下方地壳上地幔速度间断面所产生的P-S转换波及其多次反射波的信息，是研究台站下方局部区域S速度分布理想的震相，由此产生的接收函数方法是反演台站下方S波速度结构的有效手段。接收函数方法可以通过波形反演拟合接收函数的径向分量，对观测台站下方地球介质的S波速度结构进行估计，也可以通过偏移叠加获得的接收函数道集（地震剖面图）追踪速度间断面。这种方法避免了对天然地震震源及其附近结构混响效应等复杂因素的影响，对S波速度的垂向分布敏感，垂向分辨率高。由于宽频带流动地震台阵的发展，用此方法还可获得研究区域速度结构的横向变化，横向分辨能力主要取决于台站的间距。本文回顾20年来接收函数研究的进展，探讨了方法研究的发展趋势，介绍了对地壳-上地幔结构的部分研究结果。     

P波斜入射角对沉管隧道地震响应的影响

为了研究P波斜入射对沉管隧道地震响应的影响,以港珠澳大桥沉管隧道为工程背景,考虑上覆海水与海床、沉管隧道之间耦合作用,采用粘弹性边界和等效力的地震荷载输入方式,利用ADINA软件建立三维有限元模型进行地震响应分析。分析入射角为0°、20°、40°、50°、60°时P波对沉管隧道环向应力峰值（正应力峰值、剪应力峰值）和位移峰值的影响,结果表明:入射角为40°时,沉管隧道应力峰值最大;入射角为0°—40°时,隧道的应力峰值逐渐增大,入射角为40°—60°时,隧道的应力峰值逐渐减小;隧道截面4个转角处及隔墙与顶板、底板的连接处为隧道剪应力峰值最大处;隧道截面左侧剪应力峰值远大于右侧;隧道顶板正应力峰值最大,顶板的正应力峰值大约为底板的2倍;隧道截面左侧位移峰值远大于隧道截面右侧。     

利用接收函数研究程海断裂带莫霍面深度和地壳内S波速度结构特征

利用程海断裂带附近27个数字地震台站远震波形资料,提取每一个台站的接收函数,计算出各台站莫霍面深度同时利用时间域线性反演方法,获得了各个台站下方的横波速度。结果显示:程海断裂带莫霍面深度从南部42km增至北部的54km,南部和北部莫霍面深度有明显的不同。从程海断裂带下方不同深度S波速度剖面可以看出,宾川及其北东部地区中下地壳存在明显的低速层,此低速层可能与还没有固结的热物质有关。而永胜南部地区,地壳中S波速度垂直变化剧烈,低速异常高速异常交替丛生,这可能是此区地震频发的主要原因。同时,本文对宽频带地震仪和短周期地震仪得到的接收函数进行了初步的对比分析。     

一种同时反演纵波速度和泊松比的方法

本文提出一种角度部分叠加资料同时反演纵波速度和泊松比的方法.以角度部分叠加资料为基础,利用地震波振幅随入射角变化与弹性参数间的数学关系,基于非线性最优化理论,通过将最小平方问题转化为大型带状矩阵的求解问题,采用测井约束逐道外推技术依次求得该角度剖面每个点的纵波和泊松比值.通过由每个角度叠加剖面反演得到的相应角度范围的纵波速度和泊松比剖面的对比分析,进一步得到地震波垂直入射时的纵波速度和泊松比,为地震资料的岩性及含气性解释提供丰富的参数信息.