Combination of common-midpoint-refraction seismics with the generalized reciprocal method

A useful method for increasing the signal/noise ratio of refracted waves is Common-Midpoint (CMP)-refraction seismics. With this technique the shallow underground can be described in detail using all information (amplitude, frequency, phase characteristics) of the wavetrain following the first break (first-break phase). Thus, the layering can be determined and faults, weak zones, and clefts can be identified. This paper deals with the optimization of CMP-refraction seismics used in combination with the Generalized Reciprocal Method (GRM). Theoretical studies show a close relationship of both methods to the kinematics of wave propagation. Velocities and optimum offsets determined by the GRM can be used directly in the partial Radon transformation in CMP-refraction seismics. The integration of refracted waves leads to an increase in the signal/noise ratio but simultaneously the integration boundaries must be restricted to deal only with selective parts of the investigated refractor. The result of this process is an intercept-time section which can be converted directly to a depth section using standard refraction seismic techniques. Another possibility of depth conversion is the transformation of this intercept-time section to a `pseudo-zero-offset section', known from reflection seismics. Thus, zero-offset sections can be migrated using wave-equation techniques such as Kirchhoff migration.     

Seismic stratigraphy of the Mesozoic and Cenozoic in northern Belgium: main results of a high-resolution reflection seismic survey along rivers and canals

This paper presents the results of high-resolution reflection seismic surveys carried out between 1989 and 1996 along rivers and canals in northern Belgium. The seismic data penetrate down to 900 m in the sedimentary cover or to the Paleozoic basement. The reflection response of the acoustic basement provides clear indications with regard to the top of the Paleozoic: crystalline basement and Lower Paleozoic metasediments and volcanics of the London-Brabant Massif and NE-dipping Devonian and Carboniferous strata. The subhorizontal Mesozoic and Cenozoic sedimentary cover comprises 20 unconformity-bound seismic units: 5 in the Cretaceous and 15 in the Cenozoic. Based on borehole information, these units are correlated with lithostratigraphically defined formations or groups. Some of the unit-bounding unconformities are of regional importance. They are attributed i) to eustatic sea-level changes causing regional flooding during the Late Cretaceous or incision of deep valleys during the Late Oligocene and Late Miocene, ii) to regional tectonic tilting between Late Eocene and Early Oligocene, or iii) to a combination of eustasy and tectonics causing valley incisions during the Lutetian. Faults of the Roer Valley Graben have offset different stratigraphic levels by sometimes considerable amounts (up to 230 m in the Oligocene to Quaternary succession). Although the main tectonic phase took place during the Miocene, the activity has varied considerably through time, and also from fault to fault. Most faults seem to have a 10 to 30-m displacement since the Late Pliocene.     

Acquisition and Inversion of Dispersive Seismic Waves in Shallow Marine Environments

Two types of dispersive seismic waves have been acquired in different geological settings to investigate the potential to reveal the elastic parameters of the shallow marine subsurface. Scholte waves as well as acoustic guided waves are excited by a near-surface towed airgun, and recorded using two acquisition methods: (1) the towed-acquisition system using a hydrophone streamer towed close to the sea floor, and (2) the stationary-receiver method using Ocean-Bottom Seismometers and/or Hydrophones (OBS/OBH). Our diverse data sets reveal that the spatial sampling of the wavefield required to avoid aliasing may vary significantly for different geological settings. Scholte waves are characterised by a few distinct modes observed at low frequencies and low phase velocities. Their dispersion is mainly controlled by the depth profile of the shear-wave velocity. Acoustic guided waves show profound amplitude variations of numerous higher modes over a broad frequency range. These are sensitive to shear-wave velocity, but more sensitive to compressional-wave velocity than Scholte waves are. To avoid the identification of distinct modes we infer 1-D models of elastic parameters of the subsurface from the inversion of the full wavefield spectra of acoustic guided waves. In the Siberian Laptev Sea we infer the presence of a soft sediment layer (8–10 m) with a well resolved strong S-velocity gradient (150–450 m/s). In the Baltic Sea a low P-velocity layer with a strong vertical gradient (1250–1440 m/s) corresponding to a post-glacial gassy mud layer could be resolved, which agrees well with the sediment stratigraphy derived from a gravity core.     

瑞雷面波勘探的过去、现在和未来

回顾了瑞雷面波勘探发展的历史,剖析了面波频散曲线正反演的进展及存在的问题,重点指出在石油和煤炭反射地震中,应从剔除面波转变为利用面波,并借助现有石油“大炮”反射地震资料中“有害的”瑞雷面波提取浅层地质结构信息,为反射地震静校正、近地表工程地质评价等提供物性参数,把反射地震的面波干扰变为对勘探十分有利的参数提取源,实现以强干扰面波为主要数据源的浅层结构探测,以促进瑞雷面波勘探技术的进一步发展。     

Converted wave seismic imaging in the Flin Flon mining camp, Canada

Seismic imaging for mineral exploration in hardrock environment has typically relied on P-waves recorded using vertical component geophones. Here we demonstrate the potential of using P-to-S converted waves for mineral exploration by applying converted wave processing to high-resolution 3-component data acquired in the Flin Flon mining camp, Canada. Inferences regarding the relative value of PS data are made through comparison with corresponding PP images, interpretation with geological constraints, and by comparison with modeling results. Modeling results show that the horizontal component data provide a robust means of separating the PS conversions from the PP reflections in converted wave processing whereas contamination by PP energy is more severe if the vertical component data are used. A customized processing workflow has been developed in which horizontal component data are sorted into common-reflection point gathers by the application of the converted wave dip-moveout operator. The PS images are characterized by lower signal-to-noise ratios than the PP images due to the smaller S-wave impedance contrasts associated with the main lithologic units and the angular dependence of P-to-S conversion. An exception to this occurs for PS reflections from the ore zones which are accentuated on the PS images relative to the background host-rock reflectivity. There are discrepancies in the reflection dips and positions observed on corresponding PP and PS sections that are not fully explained. They can be partly attributed to the non-coincidence of subsurface spatial sampling associated with common-depth points versus common-conversion points, with an additional complication associated with the known 3D geology. Vp/Vs ratios determined by correlation of reflections observed on PP and PS images can provide a limited means of lithological discrimination.     

The dynamic response of a fluid‐filled borehole under a normal point load on the free surface

The dynamic response of a semi‐infinite fluid‐filled borehole embedded in an elastic half‐space under a concentrated normal surface load is analysed in the long‐wavelength limit. The solution of the problem is obtained with integral transforms in the form of a double integral with respect to the slowness and frequency. The partial P‐ and SV‐wave responses are further transformed to path integrals along Cagniard paths in the complex slowness plane. Unlike the traditional Cagniard‐de Hoop technique based on the Laplace transform of time dependence, this paper is based on the Fourier transform. The tube‐wave response is presented as a causal integral over a slowness range. The resultant representation in the time‐domain is suitable for the numerical evaluation of the complete response in the fluid‐filled borehole, especially at large distances. Asymptotic analysis of seismic phases arising in the borehole is performed on the basis of the obtained solution. The complete asymptotic wavefield consists in P‐ and SV‐waves, the Rayleigh wave and the low‐frequency Stoneley (tube) wave. Pressure synthetics obtained by the use of the asymptotic formulas are shown to be in good agreement with straightforward calculations.     

浙江省地震台网台站地动噪声分析

介绍浙江省数字地震台网基本情况,阐述功率谱分析的计算方法,并对各台站地动噪声计算分析。结果表明,各台站的台址环境和观测仪器的性能基本优良,能够适应宽频带地震观测要求。部分台址观测条件由于受到人为因素的影响,可以有针对性地对部分台站进行环境改造,为地震监测预报创造条件。