电火花震源系统充电技术研究

电火花震源是海洋工程地震勘探作业中一种较好的震源,但"十五"之前我国的电火花震源主要依靠进口,仪器维护费用高、周期长。"十五"期间国家海洋局第一海洋研究所研制了海洋电火花系统PSS500J,并取得了很好的应用效果。结合PSS500J研究成果,介绍了电火花震源的构架,分析了电火花震源充电技术的两种设计方案,对震源工作中安全保障提出了技术建议。     

主能量脉冲反褶积

本文提出的主能量脉冲反褶积方法,旨在保持原反褶积相位特性的同时将有效频带拓宽为期望输出的主能量谱,以期获得更为可靠的地震分辨率。其首先根据地震波的有效频带范围设计出期望输出的主能量谱,然后通过主能量谱滤波获得主能量信号,并以其为输入求取初始的反褶积因子,再在保持初始反褶积因子相位特性的前提下以主能量谱为期望的振幅谱,求取1个优化的反褶积因子,最后即可通过褶积运算获得期望输出的地震信号。     

Multi-frequency seismic study of gas hydrate-bearing sediments in Lake Baikal,Siberia

In this paper we present and discuss the frequency-dependent behaviour of the acoustic characteristics of methane hydrate-bearing sediments in Lake Baikal, Siberia. Five different types of seismic sources (airgun-array, two types of single airguns, watergun and sparker) are used, encompassing a frequency bandwidth from 10 up to 1000 Hz. On low-frequency airgun-array data, the base of the hydrate stability zone (HSZ) is observed as a high-amplitude bottom-simulating reflection (BSR) with reversed polarity. The amplitude and continuity of the BSR decrease or even disappear on medium- to high-frequency data, a feature explained in terms of vertical and horizontal resolution. The increasing reflection amplitude of the BSR with increasing offset, the calculated reflection coefficient of the BSR and the occurrence of enhanced reflections below the BSR suggest the presence of free gas below the HSZ. The observation of some enhanced reflections extending above the BSR may be interpreted as an indication for free gas co-existing with hydrates within the HSZ. Amplitude blanking above the BSR is highly variable while the BSR itself appears to act as a low-pass frequency filter for medium- to high-frequency data.New single-channel airgun profiles provide the first seismic information across the Baikal Drilling Project (BDP-97) deep drilling site, at which hydrate-bearing sediments were retrieved at about 200 m above the base of the local HSZ. At the drilling site there are no seismic characteristics indicative of the presence of hydrates. Combination of the drilling and seismic information has allowed us to make a rough estimation of the volume of hydrates and carbon stored in the sediments of Lake Baikal, which lead us to conclude that the Lake Baikal gas hydrate reservoirs do not form a prospective energy resource.     

东海深部地层时深转换关系的分段优化拟合

在将钻井垂直地震剖面(VSP)数据的时深拟合公式应用于深部地层的时深转换时,拟合深度与计算的层速度常常不够准确。本文首先利用多项式和幂函数给出了东海陆架盆地中部某凹陷41口钻井VSP数据的时深拟合公式,并用双程旅行时(TWT)最深达8 s的三维地震速度体数据与多道地震剖面对拟合公式在深部地层的适用性进行分析。在TWT为8 s时,速度体数据表明41口钻井位置的平均深度为18 140 m,平均层速度为6 208 m/s,二次多项式的平均拟合深度较之偏高9.2%,计算的层速度偏高36.2%,幂函数则分别偏低28.9%与35.6%,拟合效果都不理想。对此,本文采用通过识别VSP数据的增速拐点并对增速拐点前的VSP数据进行二次多项式拟合,对增速拐点后的VSP数据进行幂函数拟合的分段拟合模型,将TWT为8 s时的平均拟合深度和层速度的误差降到3.3%与4.7%。地震剖面显示研究区莫霍面深度约为TWT=11 s,分段拟合模型在TWT=11 s的平均拟合深度为27 516 m,层速度为7 334 m/s,更接近前人研究成果,表明该模型能显著提高深部地层时深转换的精度。     

海洋地震资料相位控制处理——东海研究实例

为了满足定量解释的需要,地震资料处理必须是高分辨率、振幅保持和相位控制的.面向海洋地震资料的相位控制处理技术包括利用远场信号进行确定性零相位化、Q补偿、预测反褶积及剩余相位校正等.利用模拟远场信号求取确定性去气泡及零相位化算子.此文基于相位控制处理,讨论了海洋地震资料的相位和极性判别方法,探讨了大于24 ms算子长度的预测反褶积对于保持信号形状的必要性,介绍了东海三维地震资料相位控制流程.     

近海底高分辨率地震探测系统设计与实现

受大深度水体的影响,传统的海面拖曳式多道地震技术在进行深水地层探测时,目标地层深度处的菲涅耳半径非常大,水平分辨率低,难以满足海域天然气水合物高精度探测需求。针对海面拖曳式地震探测技术存在的上述问题,设计了一套可以在2 000 m水深近海底作业的地震探测系统。应用耐压透声发射阵技术,克服了20 MPa外压环境和瞬时内压冲击对等离子体震源子波幅频特性的不利影响,研制的深拖等离子体震源的声源级达到214 dB,主频低于1 000 Hz;水下控制中心采用集成SoC片上系统设计,可以对震源进行定距激发控制,进行近海底多道地震数据的连续采集。系统在2019年深海试验拖曳最大深度达2 025 m,测试剖面数据显示最大地层穿透深度达380 m,纵向分辨率<2 m,横向分辨率<10 m,为深水海域沉积地层的深拖高分辨率地震探测提供了技术支撑。     

Multiple seismic attribute analyses for determination of bottom simulating reflector of gas hydrate seismic data in the Ulleung Basin of Korea

The bottom simulating reflector (BSR), the boundary between the gas hydrate and the free gas zone, is considered to be the most common evidence in seismic data analysis for gas hydrate exploration. Multiple seismic attribute analyses of reflectivity and acoustic impedance from the post-stack deconvolution and complex analysis of instantaneous attribute properties including the amplitude envelope, instantaneous frequency, phase, and first derivative of the amplitude of seismic data have been used to effectively confirm the existence of a BSR as the base of gas hydrate stability zone. In this paper, we consider individual seismic attribute analysis and integrate the results of those attributes to locate the position of the BSR. The outputs from conventional seismic data processing of the gas hydrate data set in the Ulleung Basin were used as inputs for multiple analyses. Applying multiple attribute analyses to the individual seismic traces showed that the identical anomalies found in two-way travel time (TWT) between 3.1 and 3.2 s from the results of complex analyses and l ₁ norm deconvolution indicated the location of the BSR.     

潮汕坳陷多次波组合压制技术研究

潮汕坳陷中生界发育,是我国海上重要的油气勘探后备区。但该区多次波非常发育,多次波识别和压制一直是该区地震资料处理的关键和难点。合理地选取有效的多次波压制方法是获得高品质地震剖面的关键所在。针对原始资料中多次波异常发育的情况,文章分析了该区多种类型多次波的特点,尝试了零偏移距截距时间-射线参数(τ-p)域预测反褶积、自由界面多次波衰减法(Surface-Related Multiple Elimination,SRME)和双曲线Radon变换,对多次波进行组合压制,克服了仅依靠一种手段难以压制多种类型多次波的局限。从资料处理效果看,τ-p域预测反褶积方法对短周期的鸣震和微曲多次波的压制效果较好,SRME有效地压制自由表面多次波,且在近炮检距范围内压制多次波效果很好,而双曲线Radon变换是衰减长周期的层间多次波和未被充分压制的自由表面多次波很好的选择。组合压制多次波后的剖面质量明显得到改善。     

A fiber optic recording oscilloscope for the display of side-scan sonar and seismic data

Fibre optic recording oscilloscopes may be used to improve the display and reproduction of marine seismic data. The potential of the system includes the reproduction of true plan acoustic mosaics with compensation for skew or yawing and correction of slant range distortion, and sparker record improvement by printing, for example, on transparency for optical Fourier filtering to remove multiples.     

Blind marine seismic deconvolution using statistical MCMC methods

In order to improve the resolution of seismic images, a blind deconvolution of seismic traces is necessary, since the source wavelet is not known and cannot be considered as a stationary signal. The reflectivity sequence is modeled as a Gaussian mixture, depending on three parameters (high and low reflector variances and reflector density), on the wavelet impulse response, and on the observation noise variance. These parameters are unknown and must be estimated from the recorded trace, which is the reflectivity convolved with the wavelet, plus noise. Two methods are compared in this paper for the parameter estimation. Since we are considering an incomplete data problem, we first consider maximum likelihood estimation by means of a stochastic expectation maximization (SEM) method. Alternatively, proper prior distributions can be specified for all unknown quantities. Then, a Bayesian strategy is applied, based on a Monte Carlo Markov Chain (MCMC) method. Having estimated the parameters, one can proceed to the deconvolution. A maximum posterior mode (MPM) criterion is optimized by means of an MCMC method. The deconvolution capability of these procedures is checked first on synthetic signals and then on the seismic data of the IFREMER ESSR4 campaign, where the wavelet duration blurs the reflectivity, and on the SMAVH high-resolution marine seismic data.     

莱州湾某区块浅部地层地震反射异常原因探讨

渤海莱州湾某区块采集的三维地震资料出现反射特征杂乱、同相轴连续性差等现象,对该区块浅部地层进行了电火花物探调查,并在地震资料异常区和非异常区分别选取1个位置进行钻孔,同时对海底面以下150 m范围内的土层进行连续的地质取样。通过对物探和地质资料的综合对比分析,判断造成研究区浅层地震资料反射异常与下述几个方面因素有关:(1)研究区处于莱州湾凹陷东部走滑断裂带内,浅部地层受深部断裂系统影响发育许多小断裂;(2)研究地层中富含大量的浅层气;(3)与邻近区块对比,所揭示浅部地层的土质性质十分异常,土层强度和剪切波速等均明显偏大。     

Measurement of sound-speed gradients in deep-ocean sediments usingl_{1}deconvolution techniques

A method is described for measuring the sound speed and the sound-speed gradient of surficial sea floor sediment from bottom-reflected signals recorded in marine seismic experiments. The technique makes use of the ocean-bottom impulse responses that are deconvolved from the data by means of a novel curve-fitting algorithm based on thel_{1}norm (least absolute value) criterion. The algorithm constructs the impulse response by extracting spikes one at a time in a manner that causes thel_{1}error to decrease by the maximum amount possible as each spike is chosen. Thel_{1}curve-fitting approach is a completely general strategy for deconvolution, and our algorithm can be used with data obtained from any type of marine seismic source. Since our experiments have been carried out with small explosive charges, we have also developed a method for estimating the bubble-pulse wavelet directly from the recorded bottom-reflected signal. In this paper, thel_{1}algorithm is used to deconvolve impulse responses for data obtained in an experiment in the Alaskan Abyssal Plain. The sediment-sound-speed gradient determined from these results is typical of other values reported for turbidite abyssal plains where the surficial sediments are composed of unconsolidated silty deposits.     

Reflection profiling in the deep ocean using a near bottom hydrophone

A 20 km long high resolution seismic reflection profile was carried out approximately 300 km southwest of Bermuda. The data were collected using a small airgun sound source and a single hydrophone receiver towed 100 m above the seafloor at a depth of 5400 m. Comparisons with nearby conventional seismic reflection profiles show the considerable improvement of resolving power provided by this method, particularly of the basement morphology beneath the 700 m thick sediment column. The data were recorded digitally and selected data examples show the enhancement provided by filtering, stacking, source deconvolution and corrections for hydrophone motion. The precise picture of basement topography that results from this data is compared with deep tow bathymetry profiles of the present day mid-Atlantic Ridge flanks, and is seen to be remarkably similar.     

Using at-sea experiments to study the effects of airguns on the foraging behavior of sperm whales in the Gulf of Mexico

Acoustic exposure and behavior of 8 sperm whales were recorded with acoustic and movement-recording tags before, during and after 5 separate 1–2 h controlled sound exposures of industry-provided airgun arrays. None of the 8 whales changed behavioral state (7 foraging, 1 resting) following the start of ramp-up at distances of 7–13 km, or full array exposures at 1–13 km. The most closely approached whale rested throughout exposure but started a foraging dive shortly after the airguns ceased, possibly indicating a delay in foraging during exposure. Using visual tracking and dead-reckoning of tagged animals, we found that their direction-of-movement was random with respect to the airguns, but correlated with their direction-of-movement just prior to the start of exposure, indicating that the tested whales did not show horizontal avoidance of the airguns. Oscillations in pitch generated by swimming movements during foraging dives were on average 6% lower during exposure than during the immediately following post-exposure period, with all 7 foraging whales exhibiting less pitching (p=0.014). Buzz rates, a proxy for attempts to capture prey, were 19% lower during exposure but given natural variation in buzz rates and the small numbers of whales, this effect was not statistically significant (P=0.141). Though additional studies are strongly recommended, these initial results indicate that sperm whales in the highly exposed Gulf of Mexico habitat do not exhibit avoidance reactions to airguns, but suggest they are affected at ranges well beyond those currently regulated due to more subtle effects on their foraging behavior.     

Blind Submarine Seismic Deconvolution for Long Source Wavelets

In seismic deconvolution, blind approaches must be considered in situations where reflectivity sequence, source wavelet signal, and noise power level are unknown. In the presence of long source wavelets, strong interference among the reflectors contributions makes the wavelet estimation and deconvolution more complicated. In this paper, we solve this problem in a two-step approach. First, we estimate a moving average (MA) truncated version of the wavelet by means of a stochastic expectation-maximization (SEM) algorithm. Then, we use Prony's method to improve the wavelet estimation accuracy by fitting an autoregressive moving average (ARMA) model with the initial truncated wavelet. Moreover, a solution to the wavelet initialization problem in the SEM algorithm is also proposed. Simulation and real-data experiment results show the significant improvement brought by this approach.     

Assessment of Offshore New Jersey Sources of Beach Replenishment Sand by Diversified Application of Geologic and Geophysical Methods

Beach replenishment serves the dual purpose of maintaining a source of tourism and recreation while protecting life and property. For New Jersey, sources for beach sand supply are increasingly found offshore. To meet present and future needs, geologic and geophysical techniques can be used to improve the identification, volume estimation, and determination of suitability, thereby making the mining and manag ing of this resource more effective. Current research has improved both data collection and interpretation of seismic surveys and vibracore analysis for projects investigating sand ridges offshore of New Jersey. The New Jersey Geological Survey in cooperation with Rutgers University is evaluating the capabilities of digital seismic data (in addition to analog data) to analyze sand ridges. The printing density of analog systems limits the dynamic range to about 24 dB. Digital acquisition systems with dynamic ranges above 100 dB can permit enhanced seismic profiles by trace static correction, deconvolution, automatic gain scaling, horizontal stacking and digital filtering. Problems common to analog data, such as wave-motion effects of surface sources, water-bottom reverberation, and bubble-pulse-width can be addressed by processing. More than 160 line miles of digital high-resolution continuous profiling seismic data have been collected at sand ridges off Avalon, Beach Haven, and Barnegat Inlet. Digital multichannel data collection has recently been employed to map sand resources within the Port of New York New Jersey expanded dredge-spoil site located 3 mi offshore of Sandy Hook, New Jersey. Multichannel data processing can reduce multiples, improve signal-to-noise calculations, enable source deconvolution, and generate sediment acoustic velocities and acoustic impedance analysis. Synthetic seismograms based on empirical relationships among grain size distribution, density, and velocity from vibracores are used to calculate proxy values for density and velocity. The seismograms are then correlated to the digital seismic profile to confirm reflected events. They are particularly useful where individual reflection events cannot be detected but a waveform generated by several thin lithologic units can be recognized. Progress in application of geologic and geophysical methods provides advan tages in detailed sediment analysis and volumetric estimation of offshore sand ridges. New techniques for current and ongoing beach replenishment projects not only expand our knowledge of the geologic processes involved in sand ridge origin and development, but also improve our assessment of these valuable resources. These reconnaissance studies provide extensive data to the engineer regarding the suitability and quantity of sand and can optimize placement and analysis of vibracore samples.     

Chalk fields along the Lindesnes Ridge, Eldfisk

Eldfisk field lies on the Lindesnes Ridge northwest of the Valhall and Hod fields and to the southeast of Edda in Block 2/7 of the Norwegian North Sea. The field has a complex geological history and reservoir geometry. Seismic reflection data is strongly gas affected and degraded over the majority of the field. The geological model of the field was based exclusively on well data and the model was inadequate to explain some of the production characteristics or to identify additional prospective well locations. In 1982 a programme of Vertical Seismic Profiling (VSP) was started in the northern part of the field. The seismic data provided by VSP have improved and confirmed the geological model, provided some explanations for production variances and identified further areas of interest.     

Assessment of Offshore New Jersey Sources of Beach Replenishment Sand by Diversified Application of Geologic and Geophysical Methods

Beach replenishment serves the dual purpose of maintaining a source of tourism and recreation while protecting life and property. For New Jersey, sources for beach sand supply are increasingly found offshore. To meet present and future needs, geologic and geophysical techniques can be used to improve the identification, volume estimation, and determination of suitability, thereby making the mining and manag ing of this resource more effective. Current research has improved both data collection and interpretation of seismic surveys and vibracore analysis for projects investigating sand ridges offshore of New Jersey. The New Jersey Geological Survey in cooperation with Rutgers University is evaluating the capabilities of digital seismic data (in addition to analog data) to analyze sand ridges. The printing density of analog systems limits the dynamic range to about 24 dB. Digital acquisition systems with dynamic ranges above 100 dB can permit enhanced seismic profiles by trace static correction, deconvolution, automatic gain scaling, horizontal stacking and digital filtering. Problems common to analog data, such as wave-motion effects of surface sources, water-bottom reverberation, and bubble-pulse-width can be addressed by processing. More than 160 line miles of digital high-resolution continuous profiling seismic data have been collected at sand ridges off Avalon, Beach Haven, and Barnegat Inlet. Digital multichannel data collection has recently been employed to map sand resources within the Port of New York New Jersey expanded dredge-spoil site located 3 mi offshore of Sandy Hook, New Jersey. Multichannel data processing can reduce multiples, improve signal-to-noise calculations, enable source deconvolution, and generate sediment acoustic velocities and acoustic impedance analysis. Synthetic seismograms based on empirical relationships among grain size distribution, density, and velocity from vibracores are used to calculate proxy values for density and velocity. The seismograms are then correlated to the digital seismic profile to confirm reflected events. They are particularly useful where individual reflection events cannot be detected but a waveform generated by several thin lithologic units can be recognized. Progress in application of geologic and geophysical methods provides advan tages in detailed sediment analysis and volumetric estimation of offshore sand ridges. New techniques for current and ongoing beach replenishment projects not only expand our knowledge of the geologic processes involved in sand ridge origin and development, but also improve our assessment of these valuable resources. These reconnaissance studies provide extensive data to the engineer regarding the suitability and quantity of sand and can optimize placement and analysis of vibracore samples.     

Morphosedimentary features and recent depositional architectural model of the Cantabrian continental margin

Multibeam bathymetry, high (sleeve airguns) and very high resolution (parametric system-TOPAS-) seismic records were used to define the morphosedimentary features and investigate the depositional architecture of the Cantabrian continental margin. The outer shelf (down to 180–245 m water depth) displays an intensively eroded seafloor surface that truncates consolidated ancient folded and fractured deposits. Recent deposits are only locally present as lowstand shelf-margin deposits and a transparent drape with bedforms. The continental slope is affected by sedimentary processes that have combined to create the morphosedimentary features seen today. The upper (down to 2000 m water depth) and lower (down to 3700–4600 m water depth) slopes are mostly subject to different types of slope failures, such as slides, mass-transport deposits (a mix of slumping and mass-flows), and turbidity currents. The upper slope is also subject to the action of bottom currents (the Mediterranean Water — MW) that interact with the Le Danois Bank favouring the reworking of the sediment and the sculpting of a contourite system. The continental rise is a bypass region of debris flows and turbidity currents where a complex channel-lobe transition zone (CLTZ) of the Cap Ferret Fan develops.The recent architecture depositional model is complex and results from the remaining structural template and the great variability of interconnected sedimentary systems and processes. This margin can be considered as starved due to the great sediment evacuation over a relatively steep entire depositional profile. Sediment is eroded mostly from the Cantabrian and also the Pyrenees mountains (source) and transported by small stream/river mountains to the sea. It bypasses the continental shelf and when sediment arrives at the slope it is transported through a major submarine drainage system (large submarine valleys and mass-movement processes) down to the continental rise and adjacent Biscay Abyssal Plain (sink). Factors controlling this architecture are tectonism and sediment source/dispersal, which are closely interrelated, whereas sea-level changes and oceanography have played a minor role (on a long-term scale).     

The Golo submarine turbidite system (east Corsica margin): morphology and processes of terrace formation from high-resolution seismic reflection profiles

New high-resolution seismic reflection data collected along the eastern margin of Corsica have been analysed to describe the morphology of the turbidite systems located seaward of the Golo River mouth. The boomer data reveal that there is not only one turbidite system directly associated with the river, but four additional, non-coalescing systems which grew simultaneously. In the south, the system has the typical morphology of a turbidite deposit rich in mud and sand with a well-developed meandering canyon and channel morphology. In the north, they have the morphology of sand-rich turbidite systems with shorter straight channels. The southernmost deposits are interpreted to represent a more advanced stage of turbidite system development. Terraces, recognised by their particular seismic facies on boomer profiles which clearly differs from the surrounding levee facies, are observed in the channel meanders. They are interpreted as confined levees built by vertical accretion due to deposition from low-energy flows. Despite limited penetration, boomer seismics are demonstrated to be a useful complement to lower-resolution sparker data. The boomer data are superior (1) for the characterisation of fine-grained turbidite deposits by extending seawards the limits of the turbidite systems commonly defined by the acoustic response of sands, (2) in demonstrating the persistence of turbiditic processes farther towards the basin, and (3) for proposing conceptual models for the formation of terraces in fan valleys.