共查询到20条相似文献，搜索用时 31 毫秒
1.
We investigated the seismic attenuation of compressional (P‐) and converted shear (S‐) waves through stacked basalt flows using short‐offset vertical seismic profile (VSP) recordings from the Brugdan (6104/21–1) and William (6005/13–1A) wells in the Faroe‐Shetland Trough. The seismic quality factors (Q) were evaluated with the classical spectral ratio method and a root‐mean‐square time‐domain amplitude technique. We found the latter method showed more robust results when analysing signals within the basalt sequence. For the Brugdan well we calculated effective Q estimates of 22–26 and 13–17 for P‐ and S‐waves, respectively, and 25–33 for P‐waves in the William well. An effective QS/QP ratio of 0.50–0.77 was found from a depth interval in the basalt flow sequence where we expect fully saturated rocks. P‐wave quality factor estimates are consistent with results from other VSP experiments in the North Atlantic Margin, while the S‐wave quality factor is one of the first estimates from a stacked basalt formation using VSP data. Synthetic modelling demonstrates that seismic attenuation for P‐ and S‐waves in the stacked basalt flow sequence is mainly caused by one‐dimensional scattering, while intrinsic absorption is small. 相似文献
2.
Alexander Droujinine 《Geophysical Prospecting》2005,53(1):37-63
Walkaway vertical seismic profile (VSP) data acquired in basalt‐covered areas can be used to improve knowledge of the sub‐basalt structure. A synthetic example and a case study from the North Atlantic (UK) show that elastic two‐way downward‐continuation migration combined with the stationary‐phase principle is well suited to the processing of VSP data. Vector data are processed using decoupled elastic migration algorithms in both isotropic and anisotropic media. To illustrate the value of decoupled imaging equations, conventional PP imaging is carried out on the enhanced VSP data and compared with the decoupled scheme. Decoupled vector migration operates directly on the displacement vector, and uses various wave modes. Downgoing waves are migrated to image basalt lava flows and measure their anisotropy. Upgoing waves are used for high‐resolution sub‐basalt imaging. 相似文献
3.
Due to the complicated geophysical character of tight gas sands in the Sulige gasfield of China, conventional surface seismic has faced great challenges in reservoir delineation. In order to improve this situation, a large‐scale 3D‐3C vertical seismic profiling (VSP) survey (more than 15 000 shots) was conducted simultaneously with 3D‐3C surface seismic data acquisition in this area in 2005. This paper presents a case study on the delineation of tight gas sands by use of multi‐component 3D VSP technology. Two imaging volumes (PP compressional wave; PSv converted wave) were generated with 3D‐3C VSP data processing. By comparison, the dominant frequencies of the 3D VSP images were 10–15 Hz higher than that of surface seismic images. Delineation of the tight gas sands is achieved by using the multi‐component information in the VSP data leading to reduce uncertainties in data interpretation. We performed a routine data interpretation on these images and developed a new attribute titled ‘Centroid Frequency Ratio of PSv and PP Waves’ for indication of the tight gas sands. The results demonstrated that the new attribute was sensitive to this type of reservoir. By combining geologic, drilling and log data, a comprehensive evaluation based on the 3D VSP data was conducted and a new well location for drilling was proposed. The major results in this paper tell us that successful application of 3D‐3C VSP technologies are only accomplished through a synthesis of many disciplines. We need detailed analysis to evaluate each step in planning, acquisition, processing and interpretation to achieve our objectives. High resolution, successful processing of multi‐component information, combination of PP and PSv volumes to extract useful attributes, receiver depth information and offset/ azimuth‐dependent anisotropy in the 3D VSP data are the major accomplishments derived from our attention to detail in the above steps. 相似文献
4.
The main objective of this work is to establish the applicability of shallow surface‐seismic traveltime tomography in basalt‐covered areas. A densely sampled ～1300‐m long surface seismic profile, acquired as part of the SeiFaBa project in 2003 ( Japsen et al. 2006 ) at Glyvursnes in the Faroe Islands, served as the basis to evaluate the performance of the tomographic method in basalt‐covered areas. The profile is centred at a ～700‐m deep well. VP, VS and density logs, a zero‐offset VSP, downhole‐geophone recordings and geological mapping in the area provided good means of control. The inversion was performed with facilities of the Wide Angle Reflection/Refraction Profiling program package ( Ditmar et al. 1999 ). We tested many inversion sequences while varying the inversion parameters. Modelled traveltimes were verified by full‐waveform modelling. Typically an inversion sequence consists in several iterations that proceed until a satisfactory solution is reached. However, in the present case with high velocity contrasts in the subsurface we obtained the best result with two iterations: first obtaining a smooth starting model with small traveltime residuals by inverting with a high smoothing constraint and then inverting with the lowest possible smoothing constraint to allow the inversion to have the full benefit of the traveltime residuals. The tomogram gives usable velocity information for the near‐surface geology in the area but fails to reproduce the expected velocity distribution of the layered basalt flows. Based on the analysis of the tomogram and geological mapping in the area, a model was defined that correctly models first arrivals from both surface seismic data and downhole‐geophone data. 相似文献
5.
Vertical geophone arrays in boreholes have been used for many years to study seismic velocities by investigating the first arrivals of records. The development of the vertical seismic profiling (VSP) technique shows possibilities of using the reflected events to close the gap between interpretation of conventional seismic data and physical observations made in the well. Reflected events recorded by vertical arrays (as in VSP) generally have higher signal-to-noise ratio, larger bandwidth and can easily be separated from multiples. The new Continuous Vertical Array (CVA) technique combines vertical arrays in several boreholes with a line of source points near the surface. The result is a multi-covered seismic line similar to that of a conventional seismic survey, but it retains the benefits of observations with vertical arrays. The possibilities of the new technique are discussed with the aid of theoretical considerations, model studies, and a first field case using nine boreholes 500 m apart with depths of 400 m. New data acquisition and processing techniques (mainly migration before stack) have been developed. The CVA-seismic method is still in the development stage but promises new possibilities for detailed surveys in difficult areas. 相似文献
6.
F. AMINZADEH 《Geophysical Prospecting》1989,37(8):893-906
Data from offshore Norway is used to study applications of elastic VSP modelling in detecting shear waves and observing the effects of successive mode conversion in field-recorded VSP data. The shear-wave velocities and densities from log data are used in conjunction with compressional wave velocities determined from surface seismic and log data in the VSP modelling. The time domain non-normal incidence elastic VSP modelling technique of Aminzadeh and Mendel is used as the modelling algorithm. Two surface seismograms are computed first. One is the vertical component and the other is the horizontal component for plane waves that have specified incident angles. A downward continuation method is then applied to generate seismograms at different depth points. The collection of these seismograms constitutes non-normal incidence VSPs. Both vertical and horizontal components of VSP data can be obtained by this procedure. In this paper non-normal incidence VSPs are generated for a 12.5° incident plane wave. The modelling results of layered earth systems of thin layers and thick layers are both compared with field data, and the effect of mode conversions in thin layers is observed. Several events in the field data can be explained by this elastic VSP modelling. Comparison of the model data and field data enabled a probable tube wave or out-of-plane event to be identified, the removal of which significantly improved the final VSP section. This study also shows how the VSP data helped the interpretation of the surface 3D data. 相似文献
7.
A Method of estimating attenauation from the first arrivals of VSP data is presented. The motivation is the desire to investigate the effects of scattering on wave propagation, and particularly the apparent attenuation and associated phase delay due to fine layering (the O'Doherty-Anstey effect). In order to take account of the frequency dependence of the predicted scattering attenuation, and to provide robust statistics for the estimates, a beam-forming method is used to measure the attenuation. This simularaneously estimates the slowness and polarization angle of the different wave modes, and results in attenuation measurements which are largely free of interference from reflected and mode-converted energy. By working in the frequency domain and measuring amplitude decay with depth, the frequency dependence of the attenuation is also accounted for. The beam-forming algorithm works in two passes, the first of which estimates slownesses and polarization angles over a small depth range, while the second uses the information from the first pass over a larger depth range to estimate attenuation. An approximate error analysis of the method shows that the standard variance of the estimated Q values is proportional to Q2 and the data quality (measured by its spectral coherence), and inversely proportional to the square of the analysis depth range and the square of the frequency. Hence the depth resolution is traded against the stability of the results. The method is applied to a zero-offset three-component VSP. The data are of good quality, with a bandwidth ranging from 180 Hz in the shallow part to 100 Hz in the deepest part. Stable results were obtained using a 450 m depth range. Above about 50 Hz, there is little evidence of frequency dependence in the attenuation. There is a clear division in depth into layers of higher and lower attenuation, with values of Q typically between 50 and 200. Below 50 Hz, however, attenuation increases rapidly with decreasing frequency throughout the depth range, with values of Q of less than 10 at 15 Hz. This behaviour appears anomalous since on physical grounds we expect very high values of Q at low frequency, and we have no explanation for these observations. 相似文献
8.
南黄海海相地层的地震波场特征和层位标定一直是困扰地震勘探的重要问题.为了近距离、高精度和高分辨率地观测井周围构造特征和岩石性质引起的波场变化,为地震资料的采集、处理与解释提供地震波衰减规律、速度与层位标定等信息,对大陆架科学钻探CSDP-2井实施了近零偏移距垂直地震剖面（VSP）观测.针对海相地层顶部强反射界面地震波穿透难的问题,采用了大容量气枪震源并设计了气枪阵列组合方式,提高了激发地震波的能量,获得了强反射界面之下清晰的PP、PS下行波和上行波信号.采用了三分量偏振合成、组合滤波和波场分离等处理方法,对VSP观测数据进行处理,获得了海相三叠系—志留系的精细的纵波、横波速度结构和地层吸收因子等物性数据,建立了钻井地层、测井、VSP上行波和多道地震剖面对应关系,实现了不同尺度的地质和地球物理属性资料的有效衔接,标定了钻井地质剖面上各深度地质体的地震反射特性,厘定了过井地震剖面上反射同相轴的地质属性.此次观测取得的纵波、横波速度信息,成为建立南黄海海相地层速度模型主要的资料来源,也是地震资料的岩性反演处理不可缺少的信息. 相似文献
9.
Three‐dimensional seismic survey design should provide an acquisition geometry that enables imaging and amplitude‐versus‐offset applications of target reflectors with sufficient data quality under given economical and operational constraints. However, in land or shallow‐water environments, surface waves are often dominant in the seismic data. The effectiveness of surface‐wave separation or attenuation significantly affects the quality of the final result. Therefore, the need for surface‐wave attenuation imposes additional constraints on the acquisition geometry. Recently, we have proposed a method for surface‐wave attenuation that can better deal with aliased seismic data than classic methods such as slowness/velocity‐based filtering. Here, we investigate how surface‐wave attenuation affects the selection of survey parameters and the resulting data quality. To quantify the latter, we introduce a measure that represents the estimated signal‐to‐noise ratio between the desired subsurface signal and the surface waves that are deemed to be noise. In a case study, we applied surface‐wave attenuation and signal‐to‐noise ratio estimation to several data sets with different survey parameters. The spatial sampling intervals of the basic subset are the survey parameters that affect the performance of surface‐wave attenuation methods the most. Finer spatial sampling will reduce aliasing and make surface‐wave attenuation easier, resulting in better data quality until no further improvement is obtained. We observed this behaviour as a main trend that levels off at increasingly denser sampling. With our method, this trend curve lies at a considerably higher signal‐to‐noise ratio than with a classic filtering method. This means that we can obtain a much better data quality for given survey effort or the same data quality as with a conventional method at a lower cost. 相似文献
10.
A. Pirogova R. Pevzner B. Gurevich S. Glubokovskikh K. Tertyshnikov 《Geophysical Prospecting》2019,67(7):1778-1797
At the CO2CRC Otway geosequestration site, the abundance of borehole seismic and logging data provides a unique opportunity to compare techniques of Q (measure of attenuation) estimation and validate their reliability. Specifically, we test conventional time-domain amplitude decay and spectral-domain centroid frequency shift methods versus the 1D waveform inversion constrained by well logs on a set of zero-offset vertical seismic profiles. The amplitude decay and centroid frequency shift methods of Q estimation assume that a seismic pulse propagates in a homogeneous medium and ignore the interference of the propagating wave with short-period multiples. The waveform inversion explicitly models multiple scattering and interference on a stack of thin layers using high-resolution data from sonic and density logs. This allows for stable Q estimation in small depth windows (in this study, 150 m), and separation of the frequency-dependent layer-induced scattering from intrinsic absorption. Besides, the inversion takes into account band-limited nature of seismic data, and thus, it is less dependent on the operating frequency bandwidth than on the other methods. However, all considered methods of Q estimation are unreliable in the intervals where subsurface significantly deviates from 1D geometry. At the Otway site, the attenuation estimates are distorted by sub-vertical faults close to the boreholes. Analysis of repeated vertical seismic profiles reveals that 15 kt injection of the CO2-rich fluid into a thin saline aquifer at 1.5 km depth does not induce detectable absorption of P-waves at generated frequencies 5–150 Hz, most likely because the CO2 plume in the monitoring well is thin, <15 m. At the Otway research site, strong attenuation Q ≈ 30–50 is observed only in shaly formations (Skull Creek Mudstone, Belfast Mudstone). Layer-induced scattering attenuation is negligible except for a few intervals, namely 500–650 m from the surface, and near the injection interval, at around 1400–1550 m, where Qscat ≈ 50–65. 相似文献
11.
常规陆上VSP(Vertical Seismic Profiling)勘探普遍采用纵波震源激发,三分量检波器接收,主要利用的是纵波和转换横波信息。已有的研究表明,炸药震源在井下激发、可控震源在地面垂向振动,均会产生较强的纯纵波和一定强度的纯横波;泊松比差别较大的分界面有利于形成较强的透射转换横波。本文通过对激发形成的纯横波和下行转换形成的横波进行对比分析,认为纯横波的主频往往低于纯纵波的主频,而下行转换横波的主频通常接近纵波的主频。本文分别对两个陆上纵波源零偏和非零偏VSP资料进行分析,结果表明这些资料中普遍存在纯横波,只是横波的强弱存在不同程度的变化。利用纵波源零偏VSP资料,可以获得横波速度信。最后对VSP纵波和横波联合应用前景进行了分析,应该充分利用纵波源VSP资料中的横波信息。 相似文献
12.
Fluid flow in many hydrocarbon reservoirs is controlled by aligned fractures which make the medium anisotropic on the scale of seismic wavelength. Applying the linear‐slip theory, we investigate seismic signatures of the effective medium produced by a single set of ‘general’ vertical fractures embedded in a purely isotropic host rock. The generality of our fracture model means the allowance for coupling between the normal (to the fracture plane) stress and the tangential jump in displacement (and vice versa). Despite its low (triclinic) symmetry, the medium is described by just nine independent effective parameters and possesses several distinct features which help to identify the physical model and estimate the fracture compliances and background velocities. For example, the polarization vector of the vertically propagating fast shear wave S1 and the semi‐major axis of the S1‐wave normal‐moveout (NMO) ellipse from a horizontal reflector always point in the direction of the fracture strike. Moreover, for the S1‐wave both the vertical velocity and the NMO velocity along the fractures are equal to the shear‐wave velocity in the host rock. Analysis of seismic signatures in the limit of small fracture weaknesses allows us to select the input data needed for unambiguous fracture characterization. The fracture and background parameters can be estimated using the NMO ellipses from horizontal reflectors and vertical velocities of P‐waves and two split S‐waves, combined with a portion of the P‐wave slowness surface reconstructed from multi‐azimuth walkaway vertical seismic profiling (VSP) data. The stability of the parameter‐estimation procedure is verified by performing non‐linear inversion based on the exact equations. 相似文献
13.
14.
Understanding fracture orientations is important for optimal field development of fractured reservoirs because fractures can act as conduits for fluid flow. This is especially true for unconventional reservoirs (e.g., tight gas sands and shale gas). Using walkaround Vertical Seismic Profiling (VSP) technology presents a unique opportunity to identify seismic azimuthal anisotropy for use in mapping potential fracture zones and their orientation around a borehole. Saudi Aramco recently completed the acquisition, processing and analysis of a walkaround VSP survey through an unconventional tight gas sand reservoir to help characterize fractures. In this paper, we present the results of the seismic azimuthal anisotropy analysis using seismic traveltime, shear‐wave splitting and amplitude attenuation. The azimuthal anisotropy results are compared to the fracture orientations derived from dipole sonic and image logs. The image log interpretation suggests that an orthorhombic fracture system is present. VSP data show that the P‐wave traveltime anisotropy direction is NE to SW. This is consistent with the cemented fractures from the image log interpretation. The seismic amplitude attenuation anisotropy direction is NW to SE. This is consistent with one of the two orientations obtained using transverse to radial amplitude ratio analysis, with the dipole sonic and with open fracture directions interpreted from image log data. 相似文献
15.
本文首先解决了声波方程的非Born近似的正演计算问题,从而获得理论上不带近似的正演数据;然后,推导了井间（CBP）、垂直地震剖面（VSP）和地面反射（SRP）三种不同的数据采集方式下的衍射CT的重建公式;利用这些重建算法和正演数据,系统地研究了影响到地球物理CT成象质量的三种因素,即:（1）数据采集方式,（2）异常程度和（3）成象区域的尺寸,对重建图象的影响;并比较了衍射地震CT和射线地震CT的成象质量。 相似文献
16.
We describe a method to invert a walkaway vertical seismic profile (VSP) and predict elastic properties (P‐wave velocity, S‐wave velocity and density) in a layered model looking ahead of the deepest receiver. Starting from Bayes's rule, we define a posterior distribution of layered models that combines prior information (on the overall variability of and correlations among the elastic properties observed in well logs) with information provided by the VSP data. This posterior distribution of layered models is sampled by a Monte‐Carlo method. The sampled layered models agree with prior information and fit the VSP data, and their overall variability defines the uncertainty in the predicted elastic properties. We apply this technique first to a zero‐offset VSP data set, and show that uncertainty in the long‐wavelength P‐wave velocity structure results in a sizable uncertainty in the predicted elastic properties. We then use walkaway VSP data, which contain information on the long‐wavelength P‐wave velocity (in the reflection moveout) and on S‐wave velocity and density contrasts (in the change of reflectivity with offset). The uncertainty of the look‐ahead prediction is considerably decreased compared with the zero‐offset VSP, and the predicted elastic properties are in good agreement with well‐log measurements. 相似文献
17.
本文采用2008—2017年宁夏及邻区内15个数字地震台站记录的131次地震的波形资料,基于地震波衰减与场地响应联合反演方法,分析频率介于1—7 Hz之间的Lg波衰减及场地响应。采用2.2—3.6 km/s的速度窗截取了1 069条Lg波,计算得到Lg波各分量的地震波衰减品质因子Q（ f ）与频率f之间的关系。结果显示:Lg波垂直、东西及南北向各分量的Q0分别为237.1,201.8及245.9;Q（ f ）对f的依赖性指数分别为0.44,0.52及0.44,该结果与全球其它地震活跃地区的相应研究结果相一致。对Q（ f ）在时间域的分析结果表明,研究时段内该地区地壳介质的衰减属性并未改变,15个台站的最大场地响应幅值不大于6。 相似文献
18.
Claudio Strobbia Alexander Zarkhidze Roger May Fatma Ibrahim 《Geophysical Prospecting》2014,62(5):1143-1161
Coherent noise in land seismic data primarily consists in source‐generated surface‐wave modes. The component that is traditionally considered most relevant is the so‐called ground roll, consisting in surface‐wave modes propagating directly from sources to receivers. In many geological situations, near?surface heterogeneities and discontinuities, as well as topography irregularities, diffract the surface waves and generate secondary events, which can heavily contaminate records. The diffracted and converted surface waves are often called scattered noise and can be a severe problem particularly in areas with shallow or outcropping hard lithological formations. Conventional noise attenuation techniques are not effective with scattering: they can usually address the tails but not the apices of the scattered events. Large source and receiver arrays can attenuate scattering but only in exchange for a compromise to signal fidelity and resolution. We present a model?based technique for the scattering attenuation, based on the estimation of surface‐wave properties and on the prediction of surface waves with a complex path involving diffractions. The properties are estimated first, to produce surface?consistent volumes of the propagation properties. Then, for all gathers to filter, we integrate the contributions of all possible diffractors, building a scattering model. The estimated scattered wavefield is then subtracted from the data. The method can work in different domains and copes with aliased surface waves. The benefits of the method are demonstrated with synthetic and real data. 相似文献
19.
Migration for velocity and attenuation perturbations 总被引：1，自引：0，他引：1
Migration maps seismic data to reflectors in the Earth. Reflections are not only caused by small‐scale variations of the velocity and density but also of the quality factor that describes attenuation. We investigated scattering due to velocity and attenuation perturbations by computing the resolution function or point‐spread function in a homogeneous background model. The resolution function is the migration image of seismic reflection data generated by a point scatterer. We found that the resolution function mixes velocity and attenuation parameter perturbations to the extent that they cannot be reconstructed independently. This is true for a typical seismic setting with sources and receivers at the surface and a buried scatterer. As a result, it will be impossible to simultaneously invert for velocity and attenuation perturbations in the scattering approach, also known as the Born approximation. We proceeded to investigate other acquisition geometries that may resolve the ambiguity between velocity and attenuation perturbations. With sources and receivers on a circle around the scatterer, in 2D, the ambiguity disappears. It still shows up in a cross‐well setting, although the mixing of velocity and attenuation parameters is less severe than in the surface‐to‐surface case. We also consider illumination of the target by diving waves in a background model that has velocity increasing linearly with depth. The improvement in illumination is, however, still insufficient to remove the ambiguity. 相似文献
20.
Sensitivity of time-lapse seismic to reservoir stress path 总被引：1，自引：1，他引：1
Colin M. Sayers 《Geophysical Prospecting》2006,54(3):369-380
The change in reservoir pore pressure due to the production of hydrocarbons leads to anisotropic changes in the stress field acting on the reservoir. Reservoir stress path is defined as the ratio of the change in effective horizontal stress to the change in effective vertical stress from the initial reservoir conditions, and strongly influences the depletion‐induced compaction behaviour of the reservoir. Seismic velocities in sandstones vary with stress due to the presence of stress‐sensitive regions within the rock, such as grain boundaries, microcracks, fractures, etc. Since the response of any microcracks and grain boundaries to a change in stress depends on their orientation relative to the principal stress axes, elastic‐wave velocities are sensitive to reservoir stress path. The vertical P‐ and S‐wave velocities, the small‐offset P‐ and SV‐wave normal‐moveout (NMO) velocities, and the P‐wave amplitude‐versus‐offset (AVO) are sensitive to different combinations of vertical and horizontal stress. The relationships between these quantities and the change in stress can be calibrated using a repeat seismic, sonic log, checkshot or vertical seismic profile (VSP) at the location of a well at which the change in reservoir pressure has been measured. Alternatively, the variation of velocity with azimuth and distance from the borehole, obtained by dipole radial profiling, can be used. Having calibrated these relationships, the theory allows the reservoir stress path to be monitored using time‐lapse seismic by combining changes in the vertical P‐wave impedance, changes in the P‐wave NMO and AVO behaviour, and changes in the S‐wave impedance. 相似文献