利用深反射地震数据构建的多阶面波频散曲线反演近地表横波速度结构——以跨班公湖—怒江缝合带深反射地震资料为例

深反射地震剖面法为了获取深部结构特征常常采取大的偏移距采集数据.目前公开发表的相关资料中,鲜有利用深反射地震炮集数据获取近地表的结构特征.为此,本文通过正演测试了相关数据处理流程,即利用有限差分正演了起伏地表模型的大偏移距地震单炮弹性波场特征,通过共检波点域面波信号F-K频谱叠加构建新方法,从深反射地震数据集中提取了高品质的多阶面波频散曲线,再利用多阶面波联合反演获得了近地表的结构特征.在前述正演流程基础上,利用跨越班公湖—怒江缝合带的SinoProbe深反射地震剖面中的实际炮集数据,求取了基阶和一阶瑞利波频散曲线,联合反演后得到近地表横波速度结构.该结果与初至波走时反演获取的纵波速度结构具有较好的一致性,且在近地表的浅层分辨率较纵波速度结构特征更高,而更与已有地质认识相吻合.本文提供的相关数据处理流程表明利用深反射地震炮集数据,也能够获取近地表浅层的横波速度结构.     

P and S reflection and P refraction: An integration for characterising shallow subsurface

Characterization of shallow structures was performed by using different approaches analysing both P- and S-wave seismic data with different resolution. The refraction tomography provided P and S velocity models of the first 80 m, while the reflection seismic processing gives a reasonable stacking velocity field until 300 m depth for both P- and S-wave data. So, we estimated the V_p/V_s ratio and an empirical relationship between the two velocities. We characterised the shallow layers using tomographic velocity models and the deeper layers using seismic images with different resolution. The seismic images were obtained by conventional CMP reflection seismic processing and by a novel multi-refractor imaging technique.     

The crustal structure under Sanjiang and its dynamic implications: Revealed by seismic reflection/refraction profile between Zhefang and Binchuan,Yunnan

The Sanjiang area in southwest China is considered as a tectonic intersection belt between the Tethys-Alps and the western Pacific, and has endured three-phase evolution processes: Proto-Tethys,Paleo-Tethys and Meso-Tethys[1―4]. In this area, its tectonics and struc- ture are extremely complicated, and intensively extru-sive deformation and faults are widely developed[1―3]. For that, the area is considered as the ideal na- ture-laboratory to study the evolution of Paleo-Tethys and also …     

Surface-wave analysis for static corrections in mineral exploration: A case study from central Sweden

In mineral exploration, increased interest towards deeper mineralizations makes seismic methods attractive. One of the critical steps in seismic processing workflows is the static correction, which is applied to correct the effect of the shallow, highly heterogeneous subsurface layers, and improve the imaging of deeper targets. We showed an effective approach to estimate the statics, based on the analysis of surface waves (groundroll) contained in the seismic reflection data, and we applied it to a legacy seismic line acquired at the iron-oxide mining site of Ludvika in Sweden. We applied surface-wave methods that were originally developed for hydrocarbon exploration, modified as a step-by-step workflow to suit the different geologic context of hard-rock sites. The workflow starts with the detection of sharp lateral variations in the subsurface, the existence of which is common at hard-rock sites. Their location is subsequently used, to ensure that the dispersion curves extracted from the data are not affected by strong lateral variations of the subsurface properties. The dispersion curves are picked automatically, windowing the data and applying a wavefield transform. A pseudo-2D time-average S-wave velocity and time-average P-wave velocity profile are obtained directly from the dispersion curves, after inverting only a reference curve. The time-average P-wave velocity profile is then used for the direct estimation of the one-way traveltime, which provides the static corrections. The resulting P-wave statics from the field data were compared with statics computed through conventional P-wave tomography. Their difference was mostly negligible with more than 91% of the estimations being in agreement with the conventional statics, proving the effectiveness of the proposed workflow. The application of the statics obtained from surface waves provided a stacked section comparable with that obtained by applying tomostatics.     

A METHOD TO DETERMINE THE VELOCITIES OF THE SEAFLOOR AND NEAR-SURFACE SEDIMENTS*

Shotpoint gathers from conventional reflection seismic surveys contain both reflected and refracted waves. In this study shot records were processed and analyzed, and the data were modeled with reflected, refracted, and reflected-refracted waves to fit the recorded data. The result is a detailed velocity model. The inverse problem for refracted waves was solved by using the Wiechert-Herglotz inversion. A 500-km-long 26-fold reflection seismic line from the Barents Sea, north of Norway, has been investigated. The data show high velocities, multiple reflections, and various types of noise. To test the method a total of 34 shot gathers were analyzed along this line. The aim of the interpretation was to determine the velocity in the seafloor and the near-surface sediments. It is possible to map the vertical as well as the lateral velocity distribution in detail. Depending on the length of the streamer and the velocity gradient in the sediments, the calculated depth varies between 300 and 500 m below the seafloor. These velocities were also compared to the stacking velocities obtained from the reflection seismic data to see how the velocities determined by different methods were related. The velocity distribution in the sediments is one of the key factors in seismic interpretation. The technique discussed in this paper can contribute to velocity information both in the processing and interpretation of seismic data.     

Experimental measurements of seismic velocities on core samples and their dependence on mineralogy and stress; Witwatersrand Basin (South Africa)

Physical property measurements were integrated with mineralogical analyses to better understand the nature of the seismic reflectivity of the deepest (>3.5 km depth) gold ore body (Carbon Leader Reef). The CLR lies at depths between 3.5 km and 4.5 km below the surface. Over 50 drill-core samples were selected for geochemical analyses, density and seismic velocity measurements. Ultrasonic measurements were conducted at ambient and elevated stresses, using transducers operating at 0.5 MHz. The study reveals that P-wave velocities generally increase with increasing bulk density. The CLR conglomerate, the gold-bearing reef, has slightly higher P-wave velocity (～5070–5468 m/s) and density values (～2.78 g/cm³) amongst the quartzitic units, possibly due to its massive pyrite content. The quartzite hangingwall and footwall rocks to the CLR exhibit similar P-wave velocity (～5028–5480 and ～4777–5211 m/s, respectively) and density values (～2.68 and 2.66 g/cm³, respectively). The reflection coefficients calculated at the interface between the CLR conglomerate and its hangingwall and footwall units range between ～0.02 and 0.05 which is below the required minimum reflection coefficient value of 0.06 to produce a strong reflection between two lithological boundaries. This suggests that seismic reflection methods might not be able to directly image the CLR, as observed from its poor reflectivity in the 3D seismic data. Samples were also subjected to stresses of up to 65 MPa to simulate in situ-like conditions and to investigate the dependence of seismic velocities on applied stresses. P-wave velocities increase with progressive loading, but at different rates in shale and quartzite rocks as a result of the presence of micro-defects.     

PREDICTING ABNORMALLY PRESSURED SEDIMENTARY ROCKS*

Modern seismic processing techniques developed in recent years have provided the explorationist with more meaningful data than would have been predicted even by optimists. Correct migration of seismic data, relative amplitude preservation of reflections, and seismic trace inversion represent the necessary efforts to ensure that the best possible picture of in situ physical properties of the subsurface section is revealed. Moreover, compacted and over-pressured zones can be predicted from surface data prior to drilling a well through them. The basic tool for predicting overpressured zones from the surface is still the velocity analysis derived from good reflection data with few erratic multiples. The extraction of regional normal compaction trends from the seismic velocities allows one—where velocities deviate from the trend—to locate the top of overpressure. Moreover, the statistical behavior of the ratios of the sonic log vs pore pressure in existing boreholes enables one to convert the deviation from the trend of the seismic velocities into overpressure rates expected at the seismic reflection horizon. This paper presents a field case study showing how the knowledge of well site lithology together with the more detailed information extracted from inverted seismic data enables the prediction to match well conditions with high reliability.     

Velocity analysis from large offset seismic records

One of the most important problems in applied geophysics is to extract velocities of compressional and shear waves, using the observed data collected at the Earth's surface or in boreholes. Unfortunately, in a typical seismic experiment, we do not have enough information to uniquely recover seismic velocities as functions ofx, y, andz. Thus, in the paper, a simplified model of the Earth (a stack of horizontal homogeneous layers) is considered and a critical discussion of modern techniques for processing reflection arrivals is presented.     

BOREHOLE SEISMIC PROFILING AND TUBE WAVE APPLICATIONS IN A DAM SITE INVESTIGATION*

Continuous, single-channel reflection profiling has been carried out in PVC-lined boreholes, primarily with the aim of ascertaining the position of an old subsurface gas storage tunnel on a proposed dam site. Tube wave reflection patterns thus generated have been interpreted in terms of sediment rigidity and shear wave velocity, and these results could be compared with some independent data. It is interesting to note that, within the well section penetrating Tertiary clays, the velocity of the hydraulic transients apparently was not affected by the PVC casing, which might be explained by a tight coupling between casing and clay wall. In such situations, tube waves turn out a straightforward tool for the determination of shear wave velocity and the derivation of dynamic elastic moduli of unconsolidated sediments. Further applications of the study of the distribution of seismic velocities on the dam site dealt with the consolidation history of the clays. A level of abnormally low P-wave velocities has been detected and interpreted as a gas-charged horizon which, by its coincidence with the base level of clay diapirs, might be considered to have contributed to clay flowage in past geological times. Data about maximum past burial depth, derived from shear wave velocities, turned out to be in agreement with results from consolidation testing.     

A comparison of methods for the analysis of compressional,shear, and surface wave seismic data,and determination of the shear modulus

The seismic refraction method is commonly used to determine the lithology and stratigraphic geometry of geological sites. Beyond this application there is also the potential to extract additional velocity-related information such as mechanical properties of soils and rocks. However, this requires a reliable model of the subsurface velocity variations. Refraction data, P- and SH-wave first arrivals, and surface waves were analyzed using three different techniques: delay-time in combination with ray-tracing, tomography and multichannel analysis of surface waves (MASW). Results from the first two techniques were compared, which showed that sharp high-contrasting layering is best imaged by the traditional method, delay-time followed by ray-tracing. The tomographic method was unable to detect the water table in the P-wave survey but resolved near-surface gradational velocity changes. On the other hand, in the SH-wave survey the traditional method was not useful because of gradually increasing velocities, which were better suited to the tomographic method. Furthermore, to produce spatially detailed velocity-variation models the tomographic or the MASW methods are applicable. The MASW model showed somewhat lower velocities compared to the SH-refraction tomographic model and, in contrast, showed inverted velocity gradients. This study also presents a comparison between the shear moduli measured in situ, i.e. calculated from shear wave velocities, and determined using empirical relationships. The empirical relationship for sand gives higher values for shear moduli than those measured in situ.     

Some SH-wave seismic reflections from depths of less than three metres

Shallow SH-wave reflections are far from routine, although their study can provide insights into important properties of near-surface materials that cannot be inferred from P-wave data alone. Difficulties in separating SH-wave reflections from Love waves are generally considered the major obstacle to progress in shallow SH-wave seismic reflection. This may be the case in surveys undertaken at great depths, but it is not necessarily true for reflection data gathered at shallow and ultra-shallow depths. This paper shows that when SH-wave data possess wavelengths greater than the thickness of the superficial layer, Love waves are not greatly dispersed. In this case, misinterpretation between parts of reflection hyperbolae and waveguide arrivals is sufficiently limited. In a one-layer model earth, which well approximates typical situations of the near-surface underground, the most energetic modes (the lowermost modes) of the dispersed surface waves have a dominant frequency band that falls below the wavelet spectrum of the shallow reflections; therefore, they can be filtered out in the frequency domain. Higher modes, although their spectral content overlaps that of the reflections, exhibit small amplitudes on seismograms and leave strong reflections unaffected.We present field examples from three different sites where we were able to obtain ultra-shallow reflections (< 3 m) in unconsolidated sediments. The high level of resolution (vertical resolution up to 15 cm) suggests that SH-wave reflection imaging has the potential to complement other high-resolution techniques, such as P-wave reflection and ground-penetrating radar (GPR) imaging, allowing a better and more complete characterization of the near-surface environments.     

华北克拉通北缘(怀来-苏尼特右旗)地壳结构

2009年,中国地质科学院地质研究所与美国俄克拉荷马大学合作实施了一条长453 km的深地震反射、宽角反射与折射、三分量反射地震联合探测剖面. 剖面南起怀来盆地,向北依次穿过燕山造山带西缘、内蒙地轴、白乃庙弧带、温都尔庙杂岩带,到达索伦缝合带. 其中,宽角反射与折射剖面采用8个0.5~1.5 t炸药震源激发,使用300套Texan单分量数字检波器接收,获得了高质量的地震资料. 通过资料分析和处理,识别出沉积层及结晶基底的折射波（Pg）、来自上地壳底界面的反射波（Pcp）,中地壳底界面的反射波（Plp）,莫霍界面的反射波（Pmp）及上地幔顶部的折射波（Pn）等5个震相. 分别采用Hole有限差分层析成像和Rayinvr算法对华北克拉通北缘及中亚造山带南部进行了上地壳P波速度结构成像和全地壳二维射线追踪反演成像. 结果显示：（1）中亚造山带地壳厚度~40 km,变化平缓,低于全球平均造山带地壳平均厚度,可能为造山后区域伸展的结果. 阴山-燕山带附近莫霍明显加深,推测其为燕山期造山过程形成的山根,但该山根很可能在后期被改造. （2）测线中部地壳上部速度较高,对应地表大面积花岗岩出露,而下地壳速度较低,速度梯度低,呈通道状,推测其可能曾为古亚洲洋向南俯冲消亡的主动陆缘,并在碰撞后演变为伸展环境下岩浆侵入的通道. （3）华北克拉通北缘与中亚造山带显示出不同速度变化特征,前者变化相对缓而后者则变化剧烈,二者的分界出现在赤峰-白云鄂博断裂附近.     

DECOMPOSITION AND INVERSION OF SEISMIC DATA–AN INSTANTANEOUS SLOWNESS APPROACH*

Interpretation techniques are presented that aim at the estimation of seismic velocities. The application of localized slant stacks, weighted by coherency, produces a decomposition of multichannel seismic data into single trace instantaneous slowness p(x, t) components. Colour displays support the interpretation of seismic data relevant to the near surface velocity structure. Since p(x, t) is directly related to stacking velocities and the depth of reflection, or bottoming points, in the subsurface, this data transformation provides a powerful tool for the inversion of reflection and refraction data.     

起伏地形下的高精度反射波走时层析成像方法

全球造山带及中国大陆中西部普遍具有强烈起伏的地形条件.复杂地形条件下的地壳结构成像问题像一面旗帜引领了当前矿产资源勘探和地球动力学研究的一个重要方向.深地震测深记录中反射波的有效探测深度可达全地壳乃至上地幔顶部,而初至波通常仅能探测上地壳浅部.为克服和弥补初至波探测深度的不足,本文基于前人对复杂地形条件下初至波成像的已有研究成果,采用数学变换手段将笛卡尔坐标系的不规则模型映射到曲线坐标系的规则模型,并将快速扫描方法与分区多步技术相结合,发展了反射波走时计算和射线追踪的方法.进而利用反射波走时反演,实现起伏地形下高精度的速度结构成像,从而为起伏地形下利用反射波数据高精度重建全地壳速度结构提供了一种全新方案.数值算例从正演计算精度、反演中初始模型依赖性、反演精度、纵横向分辨率以及抗噪性等方面验证了算法的正确性和可靠性.     

The seismic velocity distribution in the vicinity of a mine tunnel at Thabazimbi, South Africa

Analysis of the refracted arrivals on a seismic reflection profile recorded along the wall of a tunnel at an iron mine near Thabazimbi, South Africa, shows variations in P-wave velocity in dolomite away from the de-stressed zone that vary between 4.4 and 7.2 km/s, though values greater than 5.8 km/s predominate along most of the profile. The seismic velocities at the tunnel wall, however, vary between 4.2 and 5.2 km/s. Time–depth terms are in the range from 0.1 to 0.9 ms, and yield thicknesses of the zone disturbed by the tunnel excavations of between 2 and 9 m. The very low seismic velocities away from the tunnel wall in two regions are associated with alcoves or ‘cubbies’ involving offsets in the wall of up to 10 m. The large variations in seismic velocity resolved over distances less than 15 m with signals of wavelength around 6–9 m are attributed to variations in the sizes and concentrations of fracture systems and cracks, and in the degree of groundwater saturation of the fracture systems. The results suggest that seismic velocity variations from reflection surveys may also assist modelling studies of the stress regime in deep mines, particularly if both P and S wave velocity variations can be determined. The seismic velocity variations inferred also show that application of refraction static corrections in the processing of ‘in-mine’ seismic reflection profiles is as important as in surface surveys, because of the higher frequencies of the seismic energy recorded in the deep mine environment.     

Body-wave passive seismic interferometry revisited: mining exploration using the body waves of local microearthquakes

As the global need for mineral resources is constantly rising and the exploitable concentrations of these resources tend to become increasingly complex to explore and exploit, the mining industry is in a constant quest for innovative and cost-effective exploration solutions. In this context, and in the framework of the Smart Exploration action, an integrated passive seismic survey was launched in the Gerolekas bauxite mining site in Central Greece. A passive seismic network, consisting of 129 three-component short-period stations was installed and operated continuously for 4 months. The acquired data permitted detection of approximately 1000 microearthquakes of very small magnitude (duration magnitude ranging between –1.5 and 2.0), located within or at a very close distance from the study area. We use this microseismicity as input for the application of passive seismic interferometry for reflection retrieval, using the body waves (P- and S-wave coda) of the located microearthquakes. We retrieve by autocorrelation zero-offset virtual reflection responses, per component, below each of the recording stations. We process the acquired results using reflection processing techniques to obtain zero-offset time and depth sections, both for P- and for S-waves. In the context of the present work, we evaluate one of the acquired depth sections, using an existing seismic line passing through the Gerolekas passive seismic network, and we perform forward modelling to assess the quality and value of the acquired results. We confirm that passive seismic reflected-wave interferometry could constitute a cost-effective and environmentally friendly innovative exploration alternative, especially in cases of difficult exploration settings.     

长白山天池火山区三维地壳结构层析成像

利用长白山天池火山区三维空间深地震测观测系统所采集的反射P波走时资料,采用层析成像技术,重建了该区地壳界面构造形态和速度分析图像。界面成像结果表明：研究区地壳界面总体上显示了由北西向东南加深的趋势;马鞍山—三道白河断裂和富尔河—红旗河断裂是本区两条主要的深部构造,尤其是马鞍山—三道白河断裂,北北东向穿越天池火山口,其两侧的地壳界面存在明显的错断,预示了该处地壳厚度陡变或深大断裂带的存在,速度成像结果显示在10km深度,明显的P波低速异常分布在天池周围;15km深度上它表现为一个近南北向的P波低速异常条带,其延展尺度南北向为80~90km,东西向30~40km;随着深度的增加,P波低速异常分布在天池西侧,其尺度有明显的缩小,分布范围更加集中,而且低速扰动幅度更大。这种P波速度异常的变化图像在一定程度上反映了天池火山口下方壳内岩浆系统的空间分状况。     

Lunar noise correlation,imaging and monitoring

Passive seismic techniques have revolutionarised seismology,leading for example to increased resolution in surface wave tomography,to the possibility to monitor changes in the propagation medium,and to many new processing strategies in seismic exploration.Here we review applications of the new techniques to a very particular dataset,namely data from the Apollo 17 lunar network.The special conditions of the lunar noise environment are investigated,illustrating the interplay between the properties of the noise and the ability to reconstruct Green's functions.With a dispersion analysis of reconstructed Rayleigh waves new information about the shallow shear velocity structure of the Moon are obtained.Passive image interferometry is used to study the effect of temperature changes in the subsurface on the seismic velocities providing direct observation of a dynamic process in the lunar environment.These applications highlight the potential of passive techniques for terrestrial and planetary seismology.     

浅层地震资料解释陷阱(英文)

高分辨率浅层地震方法是在近地表调查中使用最为广泛的方法。然而,在许多情况下,地震资料的解释经常会出现错误。在本文中,我们介绍了三个例子,分析了造成P波,SH波,多道的面波(MASW)地震资料解释的错误原因,大都是由于在表面或地下条件约束不确当引起的。第一个例子是P波反射剖面上的一个波的特征被解释为浅层断裂带,但后来证实它是由高水平的背景噪音引起的,因为采集测线通过了一个公路交叉口。第二个例子是SH波反射地震剖面上一个波特征被解释为是逆倾向滑断层,但有针对性的钻探表明,它是一个侵入到基岩面的一个深层局部侵蚀。最后,第三个例子,MASW调查剖面上,一个陡倾特征一开始被解释为基岩谷。然而,后来的钻探表明这是一个非常软的湖泊沉积物,后者严重损坏了应用面波频段。虽然最初的解释是不正确的,但这刺激地球物理学家和地质学家之间的讨论,并强调地球物理数据采集的时候,采集之前以及采集之后需要科学家之间有意义的合作与讨论。     

Surface and downhole shear wave seismic methods for thick soil site investigations

Shear wave velocity–depth information is required for predicting the ground motion response to earthquakes in areas where significant soil cover exists over firm bedrock. Rather than estimating this critical parameter, it can be reliably measured using a suite of surface (non-invasive) and downhole (invasive) seismic methods. Shear wave velocities from surface measurements can be obtained using SH refraction techniques. Array lengths as large as 1000 m and depth of penetration to 250 m have been achieved in some areas. High resolution shear wave reflection techniques utilizing the common midpoint method can delineate the overburden-bedrock surface as well as reflecting boundaries within the overburden. Reflection data can also be used to obtain direct estimates of fundamental site periods from shear wave reflections without the requirement of measuring average shear wave velocity and total thickness of unconsolidated overburden above the bedrock surface. Accurate measurements of vertical shear wave velocities can be obtained using a seismic cone penetrometer in soft sediments, or with a well-locked geophone array in a borehole. Examples from thick soil sites in Canada demonstrate the type of shear wave velocity information that can be obtained with these geophysical techniques, and show how these data can be used to provide a first look at predicted ground motion response for thick soil sites.