An Attempt to Integrate Reflection Seismics and Balanced Profiles

—Different techniques in Geophysics and Geology are used to derive the structure of the subsurface. They are based on different data sets, i.e., seismic and geological data, and a combination of these techniques should produce better earth models. The case study presented in this paper is based on data of the German Continental Reflection program (DEKORP) collected in the Münsterland basin and the Rhenish Massif located at the northern border of the Rhenohercynian fold and thrust belt of the Mid-European Variscides. In this study we present an attempt to integrate balanced profiles, i.e., structural geology, and reflection seismics. The integration is performed by synthetically modelling seismic waves according to the acquisition of the field data, where the velocity model is based on the balanced profile. The synthetic data are compared with the field observations. Differences between observed data and field data are either caused by velocity errors in the model or by errors in the balanced profile. Criteria are developed to interpret these differences in order to improve the joint model of geologists and geophysicists. The case study presented in this paper shows that the combination of balanced profiles and reflection seismics may lead to shortcuts in the determination of seismic velocities of the subsurface. These shortcuts can reduce processing times and processing costs of reflection seismic data.     

A THEORY OF ACOUSTIC DIFFRACTORS APPLIED TO 2-D MODELS*

Finite-offset seismic reflection modeling of acoustic waves, propagating in a two-dimensional depth section of arbitrary complexity, is discussed. The procedure developed employs the principles of simplified (far-field) diffractor theory and ray tracing. Each reflector is represented by a set of discrete secondary sources or diffractors and the wavefield associated with each diffractor is calculated directly in the time domain by ray tracing. Reflections and diffractions are subsequently built up by the numerical superposition of these wavefields. This superposition is nondispersive for all frequencies for which the Fresnel zones are large compared with the diffractor separation. All primary travel paths connecting the shot to diffractor and diffractor to geophone are accounted for together with phase changes induced by focal events. The method allows the modeling of arbitrary trace gathers for energy originating from selected reflectors. The nonsequential nature of the algorithm makes it suited to machines capable of carrying out many similar operations in parallel or concurrently. Diffractor theory also provides physical insight into wave scattering and focusing. In particular, the half-differential waveform associated with a line diffractor leads to an explanation of the 90° phase lead induced by a cylindrical focus and, similarly, the full differential waveform of a point diffractor can be used to explain the 180° phase shift induced by a point focus.     

Diffraction imaging by prestack reverse‐time migration in the dip‐angle domain

Prestack image volumes may be decomposed into specular and non‐specular parts by filters defined in the dip‐angle domain. For space‐shift extended image volumes, the dip‐angle decomposition is derived via local Radon transform in depth and midpoint coordinates, followed by an averaging over space‐shifts. We propose to employ prestack space‐shift extended reverse‐time migration and dip‐angle decomposition for imaging small‐scale structural elements, considered as seismic diffractors, in models with arbitrary complexity. A suitable design of a specularity filter in the dip‐angle domain rejects the dominant reflectors and enhances diffractors and other non‐specular image content. The filter exploits a clear discrimination in dip between specular reflections and diffractions. The former are stationary at the specular dip, whereas the latter are non‐stationary without a preferred dip direction. While the filtered image volume features other than the diffractor images (for example, noise and truncation artefacts are also present), synthetic and field data examples suggest that diffractors tend to dominate and are readily recognisable. Averaging over space‐shifts in the filter construction makes the reflectors? rejection robust against migration velocity errors. Another consequence of the space‐shift extension and its angle‐domain transforms is the possibility of exploring the image in a multiple set of common‐image gathers. The filtered diffractions may be analysed simultaneously in space‐shift, scattering‐angle, and dip‐angle image gathers by means of a single migration job. The deliverables of our method obviously enrich the processed material on the interpreter's desk. We expect them to further supplement our understanding of the Earth's interior.     

Wave refocusing for linear diffractor using the time-reversal principle

One of the problems encountered in a variety of near-surface investigations is detecting and mapping localized heterogeneities. The heterogeneities may be classified under two kinds of objects: (1) a point diffractor that can be considered as an approximation of a small quasi-isometric, such as small karstic cavities and caves; (2) a linear diffractor roughly approximating an elongated object, such as a tube or fault plane. The point and linear diffractors generate two types of seismic diffraction: tip and edge waves, respectively. During the last few decades, different methods were proposed by many researchers for detecting these heterogeneities utilizing seismic waves diffracted by them. An alternative method for detecting point diffractors using a time-reversal principle combined with focusing analysis is proposed in this study: we present an extension of the time-reversal method for linear diffractors. It consists of a coherent summation of seismic energy along edge-diffraction traveltimes. Real data examples show the feasibility and efficiency of the proposed method.     

Geological interpretation of gravity profiles through the Karlovy Vary granite massif (Czech Republic)

We examined the shape of the Late Variscan Karlovy Vary granite massif located south of the Ohre/Eger graben in Northern Bohemia by reinterpretation of existing gravity data on two perpendicular profiles. The granite body of about 360 km² total outcrop size has the elongation ratio 0.35 with the major axis trending NE-SW. The SW part of the body was crossed in the nineties by the seismic profile 9HR which localized the bottom of granites in a depth of about 10 km. We used this value as a reference datum in our gravity profiles. We positioned one of our profiles along the seismic profile 9HR and the other one perpendicularly, i.e. parallel with the elongation of the outcrop surface. We interpret the shape of the main granite body in the vicinity of Karlovy Vary as a continuous desk whose floor is horizontal (or subhorizontal) and varies along its whole extension about a depth of 10 km. This thickness is approximately identical with that of the Saxothuringian nappes imaged by seismic reflection. The near surface upper contact of the granite body is mildly inclined, and outward dipping. It changes to steep sides or inward inclined contacts in deeper levels. The Lesny-Lysina (Kynžvart) massif is a separate granite body about 324 km thick, not continuously connected with the main Karlovy Vary massif. The gravity curve suggests that granites often enclose in their endocontact large blocks of country metasediments or metabasites the existence of which is partly evidenced by their outcrops outside the line of the profile. The granite body is found density-homogenous. Minor density differences between granite varieties are caused mainly by more intense hydrothermal alterations in younger suite granites. We interpret vertical conduits for the ascent of granitic magmas to be parallel to the Jáchymov-Gera and Ohře (Eger) lineaments or the Mariánské Lázně fault zone as indicated by the elongation of some outcrops. However, they are not clearly imaged from the gravity data. The effect of the depression of the Sokolov basin along the faults parallel with the Ohře (Eger) lineament is shallow and it is not indicated by any change in the floor depth of the granite body. Comparison of the seismicity distribution suggests that the hypocenters occur mostly outside of the granite bodies or near their contact with the country rock.     

线形区域的主辅线反射地震调查方法

本文提出主辅线反射地震法, 它由一条主测线和1～2条辅测线共同组成一个反射地震观测系统,其中辅线位于主线的一侧或两侧,大致与主线平行,只在辅线上观测而不放炮.这种系统的成本仅比原常规单线增加20%左右,但是,它却可以同时取得三方面的调查成果：(1)主线常规反射地震剖面,(2)主线与辅线之间区域的表层地震波速层析图,(3)主线与辅线之中线上的反射地震剖面.用比地震剖面略高的费用,主辅线反射地震实现了了解线形区域下方侧向可对比的地质构造变化的目标, 对线形区域的调查尤其适用,而且在中国大陆科学深钻孔区完成的主辅线反射地震调查证实了这一论断.这次主辅线反射地震法调查说明科学深钻孔区地下3000～3600m 深度还可能有规模较大的榴辉岩体存在, 现有金红石矿体储量有可能扩大若干倍.同时,主辅线反射地震法调查提供了一个后备孔位,它是将来扩大矿区储量的首选勘探孔位.     

3-D Prestack-migration of Wide-angle Data from a Variscan Transition Zone

—In addition to the near normal-incidence observations within the German DEKORP 2 project in 1984, wide-angle observations have been carried out on a parallel profile across the boundary between the Saxothuringian and Moldanubian crust, approximately 50 km NE of the main transect to control three-dimensional variations. Explosion sources have been used for the entire survey, providing excellent conditions for wide-angle registrations. A velocity model has been derived on the basis of in- and off-line refraction measurements using a kinematic raytracer which was extended to three dimensions by interpolation of 2-D velocity fields between parallel sections. Although prestack-migration of the data led to aliasing effects due to large shot and geophone spacing, stable results were obtained by forming envelopes after single-shot migration. The migrated sections reveal a strongly reflective Moho at about 31 km depth and a steeply (50°) dipping intracrustal reflector, which seems to be related to the border between the two Variscan units.     

Diffraction separation by plane-wave prediction filtering

Seismic data processing typically deals with seismic wave reflections and neglects wave diffraction that affect the resolution. As a general rule, wave diffractions are treated as noise in seismic data processing. However, wave diffractions generally originate from geological structures, such as fractures, karst caves, and faults. The wave diffraction energy is much weaker than that of the reflections. Therefore, even if wave diffractions can be traced back to their origin, their energy is masked by that of the reflections. Separating and imaging diffractions and reflections can improve the imaging accuracy of diffractive targets. Based on the geometrical differences between reflections and diffractions on the plane-wave record; that is, reflections are quasi-linear and diffractions are quasi-hyperbolic, we use plane-wave prediction filtering to separate the wave diffractions. First, we estimate the local slope of the seismic event using planewave destruction filtering and, then, we predict and extract the wave reflections based on the local slope. Thus, we obtain the diffracted wavefield by directly subtracting the reflected wavefield from the entire wavefield. Finally, we image the diffracted wavefield and obtain high-resolution diffractive target results. 2D SEG salt model data suggest that the plane-wave prediction filtering eliminates the phase reversal in the plane-wave destruction filtering and maintains the original wavefield phase, improving the accuracy of imaging heterogeneous objects.     

Focusing in Prestack Isochrone Migration Using Instantaneous Slowness Information

—Prestack migration finds increasing application in processing crustal seismic data. However, less effort has been made to incorporate slowness information in the imaging process. The combination of slowness information with migration leads to an improved image in the depth domain, especially by reducing migration artefacts and noise. A slowness-driven isochrone migration scheme is introduced for migration of 2-D seismic data. Instantaneous slowness information p(x, t) is extracted from the data using correlation analysis in moving time and space windows. Slowness values resulting from spatial coherent energy (signal) and incoherent background noise are distinguished by the simultaneous evaluation of an instantaneous coherence criterion g(x, t). In slowness-driven isochrone migration this information is used for locally weighting the amplitude A(x, t) smearing on the isochrone surface. In particular, slowness p and coherence criterion g determine position and sharpness of a Gaussian weighting function. The method is demonstrated using two synthetic data examples and is subsequently applied to two deep crustal data sets, one wide-angle (along DEKORP4) and one steep-angle reflection seismic observation (KTB8506). Both data sets were collected in the surroundings of the KTB drill site, Oberpfalz, as part of the German DEKORP project.     

地震绕射波弱信号U-net网络提取方法

绕射波携带大量小尺度非均匀地质体信息,对于提高地震勘探分辨率具有重要意义.绕射波能量远小于反射波,在地震记录中常被强反射波掩盖,因此分离并单独成像绕射波,为探测小尺度非均匀地质体的关键问题.传统绕射波分离方法受限于理论模型假设,对陡倾角反射波去除效果不佳,且易对绕射波造成损伤.基于经典编码-解码框架下的U-net网络和注意力机制,本文提出了一种绕射波智能分离方法,通过编码器自动提取地震数据中的绕射波特征,再由解码器恢复绕射波,从而隐性去除反射波.该方法作为端到端的机器学习,训练后的U-net网络可自适应地分离绕射波.本文通过数值模拟数据与实际数据构建训练数据集,利用训练后的U-net网络分离绕射波,并将结果偏移成像.数值模型测试和实际资料应用表明,融合了注意力机制的U-net网络能够有效压制反射波能量,保留绕射波动力学特征,克服了传统绕射波分离方法难以去除陡倾角反射的局限性,其提取的绕射波弱信号特征较为完整,能够进一步提高地震成像分辨率,在小尺度断裂刻画上具有优势.

     

局部倾角滤波和预测反演联合分离绕射波

地震绕射波源于介质非连续性,从地震记录中将绕射波分离出来并进行成像,其结果对研究诸如碳酸盐岩缝洞储层这类复杂非均质储层具有重要意义.对炮集记录进行平面波分解,在地层倾角不大的假设下,反射波和绕射波同相轴在平面波分解剖面上存在较大的倾角差异.基于此,我们提出分步进行绕射波分离的方法:(1)利用局部倾角滤波方法将绕射波的较大倾角信息成分分离出来,此时,余下的部分包含有反射波和残留的低倾角绕射波信息;(2)利用频率-空间域预测反演方法从上述含有反射波和残留的低倾角绕射波信息中分离出残留绕射波成分;(3)将两次分离的绕射波信息相加得到最终的绕射波估计.用该方法能够得到相对完整的绕射波信息,有效地克服了靠单一的倾角差异进行绕射波分离时明显损失低倾角信息,从而影响绕射波成像结果横向分辨率这一问题.理论与实际资料试算验证了该方法的有效性.     

二维层状介质中衍射波的计算

本文讨论了衍射波的计算,它是文献[5]的继续和发展,是为了应用射线法解决二维地震人机联作工作的一部分。根据波的物理性质,引入了等效炮点、等效检波点和等效界面等概念,将Berryhill导出的衍射波计算公式推广到倾斜界面和楔形结构等情形。这一工作与文献[1—3]中的工作相结合,在实用上解决了二维地震人机联作解释中的正演计算方法问题。     

二维层状介质中衍射波的计算

本文讨论了衍射波的计算,它是文献[5]的继续和发展,是为了应用射线法解决二维地震人机联作工作的一部分。根据波的物理性质,引入了等效炮点、等效检波点和等效界面等概念,将Berryhill导出的衍射波计算公式推广到倾斜界面和楔形结构等情形。这一工作与文献[1-3]中的工作相结合,在实用上解决了二维地震人机联作解释中的正演计算方法问题。     

腾冲火山区地壳结构的人工地震探测

介绍在腾冲地区完成的由一条近南北方向纵测线和两条近东西向的非纵测线组成的人工地震测深工程以及对纵测线资料解释的初步结果。上地壳的平均速度为５．９ｋｍ／ｓ,顶界面深度为２３ｋｍ左右;中地壳可分为２层,平均速度为６．１７ｋｍ／ｓ;纵线的平均地壳厚度为４０ｋｍ。在团田至腾冲之间的基底速度相当低（５．８０ｋｍ／ｓ）。在腾冲与固东之间的中地壳界面有局部上拱现象。局部地区的地壳速度偏低以及壳低界面上拱可能与岩     

基于倾角-偏移距域道集的绕射波成像

地震绕射波是地下非连续性地质体的地震响应,绕射波成像对地下断层、尖灭和小尺度绕射体的识别具有重要的意义.在倾角域共成像点道集中,反射波同相轴表现为一条下凸曲线,能量主要集中在菲涅耳带内,绕射波能量则比较发散.由于倾角域菲涅耳带随偏移距变化而存在差异,因此本文提出一种在倾角-偏移距域道集中精确估计菲涅耳带的方法,在各偏移距的倾角域共成像点道集中实现菲涅耳带的精确切除,从而压制反射波.在倾角-偏移距域道集中还可以分别实现绕射波增强,绕射波同相轴相位校正,因此能量弱的绕射波可以清晰地成像.在倾角域共成像点道集中,反射波同相轴的最低点对应于菲涅耳带估计所用的倾角,因此本文提出一种在倾角域共成像点道集中直接自动拾取倾角场的方法.理论与实际资料试算验证了本文绕射波成像方法的有效性.

     

平面声波在粗糙界面上的反射特征研究

基于有关粗糙界面的Rayleigh假设,讨论了平面声波按余弦规律快速变化的小尺度粗糙界面上的反射特征.研究表明：这类界面与位于该位置的一个过渡地层的作用相当.该过渡层的厚度为粗糙界面的起伏幅度,速度和密度为上下两层介质相应量的平均值.研究了埋藏很深的微粗糙界面所引起的地震反射（绕射）波的频散特性和走时的构成,即包含零炮检距反射时间、正常时差和界面粗糙时差三部分内容.该粗糙时差与空间坐标和时间坐标无关,与绕射波的阶次有关,绕射波尾随在反射波之后以某一固定的时差出现.且只有当界面的粗糙波长与地震波的波长相当时,才能观测到这类绕射波.该结论为粗糙界面地震反射资料的处理方法提供了理论依据.     

利用反射地震数据和XBT数据联合反演海水的温盐分布

利用地震数据反演海水的温盐分布是地震海洋学研究的一个重要方向.本文分析了采集于伊比利亚半岛西南缘海域的GOLR12低频地震数据和同步观测的XBT(eXpendable BathyThermograph)、CTD(Conductivity-Temperature-Depth)资料,将叠后约束波阻抗反演方法应用到海水温度、盐度剖面的计算当中,计算结果显示这种反演方法可以提供水平分辨率为6.25 m的温度和盐度资料,温度的精度约为0.16 ℃,可以为物理海洋学研究提供高横向分辨率的基础数据.     

Diffraction problems of 3D seismic imaging

Migration is essential to seismic imaging. It is carried out by backward extrapolation of the wavefield registered on the observation surface. The quality of images depends on the accuracy of the wavefield reconstruction at interior subsurface points. From the theory based on the exact solution of the scalar wave equation it is known that, for accurate wave extrapolation, data must be obtained from an infinite observation surface. Limiting of migration apertures, which is inevitable in practice, leads to artefacts in extrapolated fields. The distortion they cause in 2D and 3D imaging is different. In 2D migration, the artefacts known as truncation effects are much weaker than the signals being extrapolated and for this reason attract no special attention. In 3D migration, diffractions caused by an aperture edge are stronger and may create serious problems. For a circular aperture, their amplitudes are comparable to the amplitudes of the signals themselves. The study of aperture diffractions is intended to help in the search for ways of either suppressing them efficiently or deliberately utilizing them in order to improve imaging.
In optics, diffractions by an aperture play a constructive role in image making. This research shows that the same may take place in seismic imaging.     

南海东北部973剖面BSR及其热流特征

对南海东北部973项目采集的地震测线进行了处理,阐明了恒春海脊处天然气水合物似海底反射（BSR）的特征.不同类型剖面的分析表明,单道地震剖面可以揭示南海东北部地区BSR一定的分布,但地震处理对揭示BSR分布全貌起重要作用.南海东北部的BSR具有世界大陆边缘BSR典型的特征,切穿了沉积层理反射,与海底起伏大体平行,为一强振幅的负极性反射.该BSR特征明显,意味着南海东北部地区存在含天然气水合物沉积.利用甲烷水合物与多组分天然气水合物相平衡曲线,计算了稳定带底界的埋深,并与BSR深度进行了对比分析.无论是甲烷水合物稳定带底界还是多组分天然气水合物稳定带的底界,单一地温梯度计算的结果不可能与BSR深度在整个剖面上对应.可知本区横向上地温梯度变化较大.利用BSR资料估算了地温梯度并求得热流值.计算表明,地温梯度与热流值由西向东,随着离海沟距离的增大、离岛弧距离的减小而减小.BSR计算得到的热流值为28～64 mW/m2,与台湾西南实测热流值的结果基本可以对比.     

Low-frequency microtremor anomalies at an oil and gas field in Voitsdorf, Austria

Results of a passive microtremor survey at an oil and gas field in Voitsdorf, Austria, are presented. The survey consists in six parallel profiles approximately 9 km long over two hydrocarbon reservoirs. For each profile the seismic wavefield was recorded synchronously at 11 in-line stations. The measurements were conducted with broadband seismometers and lasted, for each profile, at least 12 hours overnight. Data interpretation is based on a comprehensive data set and on the analysis of four different spectral attributes. These attributes quantify the characteristic features of the wavefield's Fourier spectra in the low-frequency range (<10Hz). One attribute quantifies the spectral energy in the vertical wavefield component, another attribute quantifies the maxima in vertical-to-horizontal spectral ratios and two attributes describe the frequency shifts of peaks within the spectra of vertical and horizontal wavefield components. Due to temporal variations of the signals we combine the long-term measurements (several hours of continuous records) of multiple profiles. This procedure considerably enhances the consistency of each spectral attribute and makes them suitable to quantify lateral variations of the wavefield. The results show that using a combination of several attributes significantly increases the reliability of detecting anomalies in the microtremor wavefield that are presumably caused by hydrocarbon reservoirs. A numerical study of two-dimensional seismic wave propagation is applied to investigate the peak frequency shift attributes. The results of the study indicate that the attributes may contain information on the depth of hydrocarbon reservoirs, assuming that the reservoir acts as a (secondary) source of low-frequency seismic waves.