地震频发区域背景噪声提取格林函数的反褶积方法研究

针对云南这一地震频发区域记录到的地震背景噪声数据,设计了一种"最大值均一化"算法,该算法能有效地消除背景噪声中的地震信号、畸变信号及其它干扰信号。应用此算法对记录到的背景噪声数据进行预处理,对处理后的数据利用多窗谱反褶积计算方法可获得其具有较高信噪比和保留了相对幅值信息的格林函数。     

Geophones on blocks: a prototype towable geophone system for shallow land seismic investigations

In recent years some authors have given a certain amount of attention to towed seismic reflection acquisition systems. Based on some of these works we sought to design and test a system making use of currently available geophones instead of specifically designed sensors as employed in some previous works. Thus, bearing in mind that the geophone's characteristics are achieved in the conditions that they are originally designed to be employed in, i.e., coupled with a spike driven into the ground, we devoted our attention to some of the variables involved in the geophone's performance, namely the total weight, the effect of a spikeless geophone and the surfaces on which the geophone is placed. Previously, we had experimentally verified some variations in the signal response due to coupling geophones in different surface materials, such as hard soil, asphalt and concrete pavement and we noticed that these surface materials were in fact an important factor in the overall response. Hence, these materials, or as we also called them coupling agents, could be employed as a base material in the construction of a mobile seismic acquisition device composed of blocks of a certain size, on which the geophone would be then inserted and thus making it into a spikeless surface towable system. Therefore, various materials were tested in order to select one that could maintain a similar fidelity to that of the spike coupled geophone and thus contribute towards building a more time efficient and towable geophone and block system. Pinging tests revealed variations in the coupling frequency and damping characteristics of each coupling agent and from all of these tested materials one was selected for field comparative tests with the normally planted geophones with spike coupling. Finally a seismic reflection profile was acquired simultaneously with both systems, i.e., spike coupling versus cement block coupled geophones. This field test showed similar results in terms of signal levels and frequency content and therefore it was possible to observe the presence of the same reflectors and other seismic events in either of the time sections. With this experiment we propose not only a system that allows a time efficient seismic field operation but we also aim to encourage more research into the response dependency of the coupling agent of which the towable base block is made of.     

海陆两栖地区地震勘探新型检波器的研制及应用效果

本文分析了海陆两栖地带地震勘探中同时使 用的速度检波器和加速度检波器存在的差异及其对 地震信号的影响;在此基础上,设计研制出了陆用 压电检波器,并对其性能特点进行了分析;通过试 验资料分析,消除了海陆两种不同机理的检波器资 料的相位差的问题,陆用压电检波器和水中压电检 波器记录信号的频带和能量达到一致,提高了地震 资料的分辨率;实现了海陆地区可以同时采用相同 机理的检波器进行地震信号的接收,解决了滩海地 区速度检波器和加速度检波器长期混用的问题。     

Application of a New Qeneration of Geophonesto Improve Seismic Acquisitionin Onshore-offshore Transition Areas

In this paper, we review the differences between velocity geophones （VG） and acceleration geophones （AG） and their effect on seismic signals acquired in onshore-offshore transition areas. We present a new generation of Land Piezoelectric Geophone （LPG） and analyze its performance. Our field experiments demonstrate that our new LPG can be used to substitute for VGs in order to eliminate phase, frequency and energy differences between different geophone systems commonlv used in transition areas.     

Tomographic inversion of near-surface <Emphasis Type="Italic">Q</Emphasis> factor by combining surface and cross-hole seismic surveys

The estimation of the quality factor Q plays a fundamental role in enhancing seismic resolution via absorption compensation in the near-surface layer. We present a new geometry that can be used to acquire field data by combining surface and cross-hole surveys to decrease the effect of geophone coupling on Q estimation. In this study, we drilled number of receiver holes around the source hole, each hole has different depth and each geophone is placed geophones into the bottom of each receiver hole to avoid the effect of geophone coupling with the borehole wall on Q estimation in conventional cross-hole seismic surveys. We also propose a novel tomographic inversion of the Q factor without the effect of the source signature, and examine its stability and reliability using synthetic data. We estimate the Q factors of the near-surface layer in two different frequency bands using field data acquired in the Dagang Oilfield. The results show that seismic absorption in the near-surface layer is much greater than that in the subsurface strata. Thus, it is of critical practical importance to enhance the seismic solution by compensating for near-surface absorption. In addition, we derive different Q factors from two frequency bands, which can be treated, to some extent, as evidence of a frequency-dependent Q.     

Seismic interferometry experiment in a shallow cased borehole using a seismic vibrator source‡

We present the results of a seismic interferometry experiment in a shallow cased borehole. The experiment is an initial study for subsequent borehole seismic surveys in an instrumented well site, where we plan to test other surface/borehole seismic techniques. The purpose of this application is to improve the knowledge of the reflectivity sequence and to verify the potential of the seismic interferometry approach to retrieve high‐frequency signals in the single well geometry, overcoming the loss and attenuation effects introduced by the overburden. We used a walkaway vertical seismic profile (VSP) geometry with a seismic vibrator to generate polarized vertical and horizontal components along a surface seismic line and an array of 3C geophones cemented outside the casing. The recorded traces are processed to obtain virtual sources in the borehole and to simulate single‐well gathers with a variable source‐receiver offset in the vertical array. We compare the results obtained by processing the field data with synthetic signals calculated by numerical simulation and analyse the signal bandwidth and amplitude versus offset to evaluate near‐field effects in the virtual signals. The application provides direct and reflected signals with improved bandwidth after vibrator signal deconvolution. Clear reflections are detected in the virtual seismic sections in agreement with the geology and other surface and borehole seismic data recorded with conventional seismic exploration techniques.     

SIMULTANEOUS VIBROSEIS RECORDING1

An experiment was undertaken at BP's Fulbeck Geophysical test site to compare the viability of various simultaneous vibroseis recording techniques, which are often recommended as a means of improving data acquisition production rates for 3D seismic surveys. Of particular interest were: (a) the ability to separate the signals from each source during processing, (b) the generation and suppression of harmonics and (c) the effects of any source interaction. Two vibrators were deployed with a baseplate separation of 10 m, about a borehole containing a vertical array of geophones. Our analysis concentrated on the groundforce signals measured at each vibrator and the far-field signatures measured using a vertical geo-phone at a depth of 204 m. By comparing single vibrator records with similar but separated records from a simultaneous recording sequence, signal separability, harmonic suppression and vibrator interaction could be fully studied. Separated far-field signatures from simultaneous vibroseis methods using combinations of up and downsweeps exhibited unsuppressed harmonics and substantial energy from the undesired source which leaked through the correlation process. The ‘up/down’ method was capable of separating the signal from each source by only 12.7 dB, and is therefore unsuitable as a field technique. The variphase simultaneous vibroseis methods studied afforded some harmonic suppression and gave signal separations of about 30.0 dB. Use of variphase simultaneous vibroseis methods will compromise the quality of the data recorded, when compared with single-source acquisition methods. None of the simultaneous vibroseis methods tested provided adequate signal separation and, therefore, cannot be recommended as data acquisition techniques. The ‘alternate sweeping’ method coupled with multispread recording will give the desired improvement in data acquisition rates, while preserving the necessary quality of our seismic data.     

Site characterization in the Vega Baja, SE Spain, using ambient-noise H/V analysis

New available geological, geotechnical, and geophysical data and new acquired ambient-noise surveys allowed for a detailed characterization of the soils in the Vega Baja region, located in the Lower Segura River basin, SE of the Iberian Peninsula. Single-station measurements were carried out at 90 sites, mainly over the central area of the basin, checking reproducibility by repeating some measurements under different conditions, on different dates, and through different acquisition systems, following well established standard procedures. The detailed data processing and analysis of the new acquired ambient-noise records provided the identification of multiple peaks in some of the horizontal-to-vertical spectral ratio curves, representing different impedance contrasts at depth that were 1D modelled to characterize the Vega Baja soils in terms of resonant frequencies, average shear-wave velocity, thickness of the different soil units, depth to the engineering bedrock and depth to the rock units when the model includes it. The depth to the engineering bedrock is around 30 m in the central part of the region, reaching maximum values of about 50 m in the north and to the east. Maximum values of the shear-wave velocity averaged over the soil units are close to 200 m/s while that averaged in the top 30 m reach values around 300 m/s south of the river bank. Contour mapping and cross-sections of the obtained model parameters illustrate the soil characteristics over the region, providing relevant information for the generation of representative sets of seismic hazard scenarios at different probability levels in future studies for their application in risk mitigation and emergency planning in the Vega Baja.     

Construction of regional and local seismic anisotropic structures from wide-angle seismic data: crustal deformation in the southeast of China

We present a method to obtain spatial distributions of seismic anisotropy associated with regional stress and local faulting in the crust from wide-angle seismic data. The method contains three steps. The first step consists of obtaining radial- and transverse-component seismic sections using a pre-stack depth migration algorithm from the S-wave velocity model determined by conventional interpretation of picked intra-crustal seismic events. In the second step, we compute time delays between split shear-waves and polarizations of fast split shear-waves by minimizing the transverse-component seismic energy. The time delay and polarization in each layer are derived using a layer-stripping method. The final step is to estimate the average splitting parameters along the whole profile. Thus, the average time delay and polarization can be regarded as caused by the effects owing to regional structure and stress fields, whereas the residual values of the splitting parameters are considered to be related to local structures and local faulting. Our method allows us to construct multi-layer anisotropic images, which may later be interpreted in terms of intra-layer coupling/decoupling or deformation. We present results from a set of three-component seismic data acquired by a controlled source experiment in the southeast region of China. The results demonstrate that the average polarizations and time delays are consistent with the direction and strength of the stress field, and their lateral variations related to local anisotropy match the spatial distribution of surface faulting crossing the acquisition seismic profile.     

基于稀疏约束和多源激发的地震数据高效采集方法

地震勘探目标日趋复杂化和精细化,"两宽一高"等采集技术获得了广泛应用,从而导致当前地震数据采集周期越来越长、成本越来越高,如何解决日益增长的勘探成本问题成为当前地震采集领域的研究热点之一.针对上述问题,本文首先开展了基于稀疏性的地震数据高效采集方法理论研究,对地震数据稀疏性基本理论、稀疏约束下随机采样及其数据重建方法进行了深入探讨,提出使用改进的分段随机采样方法灵活地进行实际地震采集测网设计;详细阐述了多源地震激发方法,对多源地震数据分离方法开展了深入研究,提出了基于小窗口中值滤波与稀疏约束联合随机去噪的多源数据分离方法,并在数据分离处理中取得了较好的效果;将上述两种地震数据采集方案有机结合,提出了1）规则多源、随机检波点（DmsRg）、2）随机多源、规则检波点（RmsDg）和3）随机多源、随机检波点（RmsRg）等三种高效采集方案及相应的数据重建方案,满足了后续常规化数据处理的要求,并讨论了多源激发对数据成像的影响.基于Marmousi模型数据的数值试验表明,本文构建的基于稀疏约束和多源激发的高效采集方法理论对于提高地震数据采集效率、降低勘探成本具有重要的应用价值,建立的数据重建方法流程可以取得和常规数据接近的成像结果.本文方法虽然在数值试验中取得了较为理想的效果,但还需要得到野外实际数据采集的进一步检验.     

Geophone orientation and coupling in three-component sea-floor data: a case study

We analyse the geophone orientation and coupling in a data set from the North Sea. Based on the polarization of the water-break on the sea-floor, we have derived processing algorithms for determining the receiver orientation for gimballed and non-gimballed geophone systems. For a gimballed system, the problem reduces to a simple horizontal rotation. However, for a non-gimballed system, where all three geophone axes may vary due to varying acquisition conditions such as dipping sea-floor, twisting of recording cable, etc., the three orientation angles cannot be found directly from the recorded displacement vectors. Using the data redundancy within a common-receiver gather, a robust two-stage method is derived for the non-gimballed system in which all three orientations can initially be unknown. Testing on the North Sea data set acquired with a gimballed system shows that the three-component geophones in the data set are orientated satisfactorily within an error of 5°. However, there are some undesirable cross-couplings between the vertical and horizontal geophones, which results in leakage of shear-wave energy from the horizontal components to the vertical components.     

Orthogonal vibroseis sweeps

Vibroseis is a method that imparts coded seismic energy into the ground. The energy is recorded with geophones and then processed using the known (coded) input signal. The resulting time‐domain representation of vibroseis data is an impulsive wavetrain with wavelet properties consistent with the coded input signal convolved with the earth's reflectivity series. Historically, vibratory seismic surveys collect data from one source location at a time, summing one or more sources at each location. We present a method of designing orthogonal sweeps using the concept of combisweeps. The orthogonal sweeps allow simultaneous recording and later separation of two or more unique source locations. Orthogonality of sweeps permits separation of the data into unique source‐location field records by a conventional correlation procedure. The separation power of the orthogonal sweeps is demonstrated by a comparison between separated data and data acquired with one vibrator. Separation noise was at a negligible level for our demonstration data sets when two vibrators were located 50 m to 200 m apart. Coincident generation and recording of two vibroseis sweeps at different locations would allow almost double the amount of data to be recorded for a given occupation time and requires only half the storage medium.     

陆上高分辨率地震勘探检波器性能及应用效果分析(英文)

地震检波器的性能是高分辨率地震数据采集中的一个重要因素,对资料的品质影响很大。为此,作者在不同类型地区对目前在高分辨率地震勘探中常用的动圈式检波器、涡流检波器和数字检波器的性能进行了对比试验。通过野外工作的实际,总结了动圈式检波器、涡流检波器和数字检波器在不同表层地质条件地区的应用效果,提出了适合该类地区高分辨率地震检波器的性能指标,指出合理选择检波器类型和检波器的联接方式,能提高地震资料的信噪比和分辨率。     

Vertical seismic profiling using a daisy‐chained deployment of fibre‐optic cables in four wells simultaneously – Case study at the Ketzin carbon dioxide storage site

The geological storage of carbon dioxide is considered as one of the measures to reduce greenhouse gas emissions and to mitigate global warming. Operators of storage sites are required to demonstrate safe containment and stable behaviour of the storage complex that is achieved by geophysical and geochemical monitoring, combined with reservoir simulations. For site characterization, as well as for imaging the carbon dioxide plume in the reservoir complex and detecting potential leakage, surface and surface‐borehole time‐lapse seismic monitoring surveys are the most widespread and established tools. At the Ketzin pilot site for carbon dioxide storage, permanently installed fibre‐optic cables, initially deployed for distributed temperature sensing, were used as seismic receiver arrays, demonstrating their ability to provide high‐resolution images of the storage formation. A vertical seismic profiling experiment was acquired using 23 source point locations and the daisy‐chained deployment of a fibre‐optic cable in four wells as a receiver array. The data were used to generate a 3D vertical seismic profiling cube, complementing the large‐scale 3D surface seismic measurements by a high resolution image of the reservoir close to the injection well. Stacking long vibro‐sweeps at each source location resulted in vertical seismic profiling shot gathers characterized by a signal‐to‐noise ratio similar to gathers acquired using geophones. A detailed data analysis shows strong dependency of data quality on borehole conditions with significantly better signal‐to‐noise ratio in regions with good coupling conditions.     

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

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

VSP地震干涉成像及应用研究

受井中检波器串级数局限,垂直地震剖面（VSP）反射波成像范围窄,且不能对井中最浅接收点上方区域有效成像.虽然多次波成像能扩大成像范围,但在实际应用中尚有诸多困难和挑战.本文根据Wapenaar的地震干涉理论,基于上下行波场分离技术,研发了VSP地震干涉成像方法.该方法将VSP自由表面多次波重建为在地表震源位置激发（虚震源）接收的拟地面地震反射波,然后偏移成像,以达到对多次波间接成像的目的.通过数值模型实验,测试了VSP干涉成像的极限分辨率,并讨论了主要采集参数的影响,结果表明：该方法的垂向和水平极限分辨率分别达约10 m和20 m,且能分辨深度达6500 m处的50 m×100 m溶洞;采用12至24道井中检波器串采集的VSP资料,其干涉成像结果显著优于VSP反射波成像,与相应的地面地震成像效果相当.将本文方法应用于新疆地区采集的VSP资料,结果表明：与VSP反射波成像相比,成像同相轴更加连续,成像范围显著扩大;与地面地震成像相比,成像结果相当,尤其在浅中部甚至更好.新方法不仅无需进行井中接收点静校正,且能显著增大成像范围,有利于成像同相轴的追踪对比、地震属性提取和地质解释,尤其对中国新疆地区深部缝洞型储层的成像,具有广泛的实际应用前景.

     

沙漠地区地震勘探随机噪声建模及其在噪声压制中的应用

随机噪声是影响地震勘探有效信号的主要因素,其存在大大降低了地震记录的信噪比.在噪声压制方法不断被改进的同时,对随机噪声特性进行研究,了解噪声的产生机制是对其进行压制的先决条件,目前对噪声的研究主要是特性研究以寻找规律性,对其进行定性定量的分析还比较少.本文根据塔里木沙漠地区实际采集环境,考虑到噪声的连续性给计算带来的不便,假设各类噪声源以点源的形式分布在检波器周围,依据相应理论确定各类噪声源的源函数,其激发的噪声经由波动方程传播,将随机噪声作为各类噪声源共同作用的综合波场,建立随机噪声的理论模型.通过分析不同种噪声对地震记录的影响,选取合适的滤波方法对其进行压制,实验结果表明,通过建立沙漠地区随机噪声的理论模型,为选择有效的滤波方法,提高地震记录信噪比起到理论指导作用.     

ON THE MINIMIZATION OF THE DISTORTION CAUSED BY THE GEOPHONE-GROUND COUPLING *

It is well known that the application of the “bright spot’ technique has been more successful in marine prospecting than in land prospecting. This is due partly to the problem of distortion of the seismic signal caused by the geophone-ground coupling, especially when carrying out high resolution, shallow seismic surveys in swampy terrain. The effect of geophone-ground coupling on the response of a single geophone to the incident compressional waves has been treated by several authors. However, they have always neglected the influence of mutual interaction between an array of geophones on the response of each geophone forming the array. We show that mutual interaction, which results from the re-radiation of the incident compressional waves by the geophones forming the array, can have considerable effect on the response of each geophone. The effect of the geophone-ground coupling on the response of a seismic channel is considered in the absence and presence of mutual interaction between a group of geophones for the case when the shear wave velocity of the soil varies by a factor of three.     

基于多用户峰度准则的海洋强噪声衰减方法(英文)

海洋地震勘探过程中,由于采集设备的老化或电源的不稳定而造成的漏电,在地震记录表现为强噪音干扰,利用常规噪音衰减方法处理此类强噪音效果并不理想。鉴于强噪音在统计学上具有相同的特性,本文在基于峰度的盲分离(blind source separation,BSS)算法研究基础上,推导出一种基于多用户峰度(multiuser kurtosis,MUK)准则的噪音衰减算法来估计地震记录中具有相同统计特性的强噪音,并将其从地震记录中分离,从而达到衰减强噪音的目的。模型试验与实际资料的处理表明:该方法能够在好的衰减海洋地震勘探记录中的强噪音,保留了更多的有效信息,提高海洋地震数据的信噪比,具有可行性和应用前景。     

Burying receivers for improved time‐lapse seismic repeatability: CO2CRC Otway field experiment

4D seismic is widely used to remotely monitor fluid movement in subsurface reservoirs. This technique is especially effective offshore where high survey repeatability can be achieved. It comes as no surprise that the first 4D seismic that successfully monitored the CO₂ sequestration process was recorded offshore in the Sleipner field, North Sea. In the case of land projects, poor repeatability of the land seismic data due to low S/N ratio often obscures the time‐lapse seismic signal. Hence for a successful on shore monitoring program improving seismic repeatability is essential. Stage 2 of the CO2CRC Otway project involves an injection of a small amount (around 15,000 tonnes) of CO₂/CH₄ gas mixture into a saline aquifer at a depth of approximately 1.5 km. Previous studies at this site showed that seismic repeatability is relatively low due to variations in weather conditions, near surface geology and farming activities. In order to improve time‐lapse seismic monitoring capabilities, a permanent receiver array can be utilised to improve signal to noise ratio and hence repeatability. A small‐scale trial of such an array was conducted at the Otway site in June 2012. A set of 25 geophones was installed in 3 m deep boreholes in parallel to the same number of surface geophones. In addition, four geophones were placed into boreholes of 1–12 m depth. In order to assess the gain in the signal‐to‐noise ratio and repeatability, both active and passive seismic surveys were carried out. The surveys were conducted in relatively poor weather conditions, with rain, strong wind and thunderstorms. With such an amplified background noise level, we found that the noise level for buried geophones is on average 20 dB lower compared to the surface geophones. The levels of repeatability for borehole geophones estimated around direct wave, reflected wave and ground roll are twice as high as for the surface geophones. Both borehole and surface geophones produce the best repeatability in the 30–90 Hz frequency range. The influence of burying depth on S/N ratio and repeatability shows that significant improvement in repeatability can be reached at a depth of 3 m. The level of repeatability remains relatively constant between 3 and 12 m depths.