共查询到20条相似文献，搜索用时 375 毫秒
1.
Multicomponent seismic data are acquired by orthogonal geophones that record a vectorial wavefield. Since the single components are not independent, the processing should be performed jointly for all the components. In this contribution, we use hypercomplex numbers, specifically quaternions, to implement the Wiener deconvolution for multicomponent seismic data. This new approach directly derives from the complex Wiener filter theory, but special care must be taken in the algorithm implementation due to the peculiar properties of quaternion algebra. Synthetic and real data examples show that quaternion deconvolution, either spiking or predictive, generally performs superiorly to the standard (scalar) deconvolution because it properly takes into account the vectorial nature of the wavefields. This provides a better wavelet estimation and thus an improved deconvolution performance, especially when noise affects differently the various components. 相似文献
2.
3.
应用多分量地震资料进行成像时通常需要先做波场分离,然后再对分离的波型进行成像。其中,波场分离可以在空间域或波数域实现。然而,由于用交错网格有限差分进行弹性波场数值模拟时,用来进行波数域波场分离的质点振动速度分量定义在不同网格节点上,本文提出了利用波数域插值方法来估算同一网格节点所需质点振动速度值;进而给出了先进行波数域插值后进行波场分离的波数域保幅波场分离方案。数值实验结果表明波数域插值方法具有较高的插值精度且保幅波场分离方法具有较好的保幅性,将本文方法进一步应用于弹性波逆时偏移可以获得保幅性较好的成像结果且对存在一定程度速度误差情况具有较好的适应性。 相似文献
4.
The most difficult part of multicomponent processing is the estimation of the shear-wave velocity map for migration. We used refracted shear waves and a simple iterative method called wavefield continuation (WFC) to evaluate the shallow shear-wave velocity profile on a real data example. The WFC was developed in 1981 by Clayton and McMechan to determine compressional-wave velocity profiles from refracted compressional waves. The application to refracted shear waves is straightforward. The real data example shows that shear structure can be easily determined independently of the compressional structure. 相似文献
5.
P. B. DILLON 《Geophysical Prospecting》1990,38(7):757-777
A modern approach to migration is to perform wavefield extrapolation, subject to an imaging condition. Correct wavefield extrapolation requires that the boundary conditions at the array of geophones satisfy the wave equation. A sufficient condition is to perform the survey with a single stationary source. Contrary to this condition, many VSPs are conducted in deviated wells, where the source is maintained vertically above the down-hole geophone at each well station. Such a survey fails to provide the boundary conditions theoretically necessary for wave-equation migration. A recently published inversion scheme, referred to as acoustic generalized Radon transform migration (GRT migration), was developed to handle any configuration of sources and geophones, including moving-source deviated-well VSP surveys. GRT migration may be viewed as a weighted version of the generalized Kirchhoff migration, derived in this paper from the exploding-reflector model. When a VSP-survey geometry has been specified, GRT migration can be expressed in terms of array parameters, and compared with the equivalent expression for Kirchhoff (wave-equation) migration. The differences between the two integrals are significant and their effect is demonstrated on VSP data. 相似文献
6.
Reflection and refraction data are normally processed with tools designed to deal specifically with either near- or far-offset data. Furthermore, the refraction data normally require the picking of traveltimes prior to analysis. Here, an automatic processing algorithm has been developed to analyse wide-angle multichannel streamer data without resorting to manual picking or traveltime tomography. Time–offset gathers are transformed to the tau–p domain and the resulting wavefield is downward continued to the depth–p domain from which a velocity model and stacked section are obtained. The algorithm inputs common-depth-point (CDP) gathers and produces a depth-converted stacked section that includes velocity information. The inclusion of long-offset multichannel streamer data within the tau–p transformation enhances the signal from high-velocity refracted basalt arrivals. Downward continuation of the tau–p transformed wavefield to the depth–p domain allows the reflection and refraction components of the wavefield to be treated simultaneously. The high-slowness depth–p wavefield provides the velocity model and the low-slowness depth–p wavefield may be stacked to give structural information. The method is applied to data from the Faeroe Basin from which sub-basalt velocity images are obtained that correlate with an independently derived P-wave model from the line. 相似文献
7.
S. Parolai D. Bindi S. Ullah S. Orunbaev S. Usupaev B. Moldobekov H. Echtler 《Journal of Seismology》2013,17(2):707-719
We present results from a vertical array of accelerometers that was recently installed in Bishkek (Kyrgyzstan) with the long-term aim of recording strong motion data. Taking advantage of recordings of a Mb 4.7 earthquake that occurred 40 km from the array site during the installation phase, we provide results of some preliminary data analysis. First, estimates of the S-wave velocity and Qs structure are deduced by the inversion of the deconvolved wavefield between the sensors in the borehole. Furthermore, the application of the nonstationary ray decomposition Kinoshita (Earth Planets Space 61:1297-1312, 2009) allowed at least three reflectors in the shallow velocity structure below the array to be identified. The complex nature of the wavefield (with up-going, down-going waves, and converted phases) due to the coarse, unconsolidated subsoil structure is highlighted by means of numerical simulations of ground motion. 相似文献
8.
双程波方程逆时深度偏移是复杂介质高精度成像的有效技术,但其结果中通常包含成像方法引起的噪音和假象,一般的滤波方法会破坏成像剖面上的振幅,其中的假象也会给后续地质解释带来困扰.将波场进行方向分解然后实现入射波与反射波的相关成像能够有效地消除这类成像噪音,并提高逆时偏移成像质量.波传播方向的分解通常在频率波数域实现,它会占用大量的存储和计算资源,不便于在沿时间外推的逆时深度偏移中应用.本文提出解析时间波场外推方法,可以在时间外推的每个时间片上实现波传播方向的显式分解,逆时深度偏移中利用分解后的炮检波场进行对应的相关运算,实现成像噪音和成像信号的分离.在模型和实际数据上的测试表明,相比于常规互相关逆时偏移成像结果,本文方法能够有效地消除低频成像噪音和特殊地质构造导致的成像假象. 相似文献
9.
A method for interval velocity analysis is formulated on the basis of wavefield extrapolation, i.e., on the basis of wave-equation migration. When this scheme is applied to multioffset seismic sections or to an ensemble of CMP gathers, it allows for the proper treatment of dipping events. The underlying assumptions are that local velocities should be derived from data associated with events within the interval under consideration. To minimize the effect of the region above the layer of interest, the data are first extrapolated to the top of the analysis interval. Subsequent analysis of these data then pertains to the events within this interval. Velocity estimation consists of repeated wavefield extrapolations through the analysis interval using a set of trial velocities. The optimal velocity is chosen on the basis of coherency measures designed to express the collective phase agreement among a set of offset Fourier modes. The reliability of this approach to interval velocity estimation is demonstrated on synthetic multi-offset data. 相似文献
10.
Up–down wavefield decomposition is effectuated by a scaled addition or subtraction of the pressure and vertical particle velocity, generally on horizontal or vertical surfaces, and works well for data given on such surfaces. The method, however, is not applicable to decomposing a wavefield when it is given at one instance in time, i.e. on snapshots. Such situations occur when a wavefield is modelled with methods like finite-difference techniques, for the purpose of, for example, reverse time migration, where the entire wavefield is determined per time instance. We present an alternative decomposition method that is exact when working on snapshots of an acoustic wavefield in a homogeneous medium, but can easily be approximated to heterogeneous media, and allows the wavefield to be decomposed in arbitrary directions. Such a directional snapshot wavefield decomposition is achieved by recasting the acoustic system in terms of the time derivative of the pressure and the vertical particle velocity, as opposed to the vertical derivative in up–down decomposition for data given on a horizontal surface. As in up–down decomposition of data given at a horizontal surface, the system can be eigenvalue decomposed and the inverse of the eigenvector matrix decomposes the wavefield snapshot into fields of opposite directions, including up–down decomposition. As the vertical particle velocity can be rotated at will, this allows for decomposition of the wavefield into any spatial direction; even spatially varying directions are possible. We show the power and effectiveness of the method by synthetic examples and models of increasing complexity. 相似文献
11.
Xiaoxue Jiang Fan Zheng Haiqing Jia Jun Lin Hongyuan Yang 《Studia Geophysica et Geodaetica》2016,60(1):91-111
A time-domain hyperbolic Radon transform based method for separating multicomponent seismic data into P-P and P-SV wavefields is presented. This wavefield separation method isolates P-P and P-SV wavefields in the Radon panel due to their differences in slowness, and an inverse transform of only part of the data leads to separated wavefields. A problem of hyperbolic Radon transform is that it works in the time domain entailing the inversion of large operators which is prohibitively time-consuming. By applying the conjugate gradient algorithm during the inversion of hyperbolic Radon transform, the computational cost can be kept reasonably low for practical application. Synthetic data examples prove that P-P and P-SV wavefield separation by hyperbolic Radon transform produces more accurate separated wavefields compared with separation by high-resolution parabolic Radon transform, and the feasibility of the proposed separation scheme is also verified by a real field data example. 相似文献
12.
In this paper, regionalP-wave upper mantle structure is investigated using slant-stack velocity analysis of short-period earthquake data recorded at station MAJO (Matsushiro, Japan). Shallow earthquakes from 1980–1986 within 35° of MAJO are used to construct a common receiver gather. Processing of the wavefield data includes focal depth and static time corrections, as well as deterministic deconvolution, in order to equalize pulse shapes and align wavelets on the first arrivals. The processed wavefield data are slant stacked and interatively downward continued to obtain a regional upper mantle velocity model. The model includes a low velocity zone between 107 and 220 km. Beneath the LVZ, the velocity increases smoothly down to the discontinuity at 401 km. In the transition zone, the velocity model again increases linearly, although there is some suggestion of further complexity in the downward continued wavefield data. At the base of the transition zone, a second velocity discontinuity occurs at 660 km, with a linear velocity gradient below. In addition to slant-stack analysis, travel times and synthetic seismograms are computed and compared with the processed and unprocessed wavefield data. 相似文献
13.
14.
Multichannel analysis of surface wave method with the autojuggie 总被引：2，自引：0，他引：2
Gang Tian Don W. Steeples Jianghai Xia Richard D. Miller Kyle T. Spikes Matthew D. Ralston 《Soil Dynamics and Earthquake Engineering》2003,23(3):61-65
The shear (S)-wave velocity of near-surface materials and its effect on seismic-wave propagation are of fundamental interest in many engineering, environmental, and groundwater studies. The multichannel analysis of surface wave (MASW) method provides a robust, efficient, and accurate tool to observe near-surface S-wave velocity. A recently developed device used to place large numbers of closely spaced geophones simultaneously and automatically (the ‘autojuggie’) is shown here to be applicable to the collection of MASW data. In order to demonstrate the use of the autojuggie in the MASW method, we compared high-frequency surface-wave data acquired from conventionally planted geophones (control line) to data collected in parallel with the automatically planted geophones attached to steel bars (test line). The results demonstrate that the autojuggie can be applied in the MASW method. Implementation of the autojuggie in very shallow MASW surveys could drastically reduce the time required and costs incurred in such surveys. 相似文献
15.
基于弹性波解耦延拓方程的波场分离方法不仅可以得到解耦的纵、横波质点振动速度场,还可以获得纵、横波应力场.针对横波纯应力场在利用单一分量成像时不具有明确物理意义的问题,本文将应力偏张量引入到横波应力场中,基于应力偏量第二不变量构建得到横波应力不变量并将其用于应力场的逆时偏移成像中,获得了可完整表征横波应力场的成像结果.模型试算表明,本文构建的横波应力不变量可以有效利用横波应力张量中的波场信息,并得到准确的弹性逆时偏移成像结果. 相似文献
16.
Xu Zhaotao Wang Huazhong 《应用地球物理》2006,3(4):226-233
Based on perturbation theory, the wave equation extrapolation operator with mixed domains has the ability to deal with lateral velocity variations. It is the image method that has undergone much research in seismology. All extrapolation operators face the problem of choosing the reference velocity due to continuation in depth. The wavefield extrapolation operator with a single reference velocity is suitable for media with weak lateral variation. The multi-reference velocity extrapolation operator can cope with severe lateral velocity variations and improve image accuracy. However, the calculation cost is large. We present a self-adaptive approach to automatically determine the number of selected reference velocities according to the complexity of structure and the given velocity threshold value. The approach can be used to construct the SSF, FFD, WXFD, and GSP multi-reference velocity wavefield extrapolation image algorithms. The result of a salt-dome model data test demonstrates that the self-adoptive multi-reference wavefield extrapolation algorithm has the ability to deal with severe lateral velocity variations and can also be used for structure edge detection. The method is flexible and computationally cost-effective. 相似文献
17.
Field testing of modular borehole monitoring with simultaneous distributed acoustic sensing and geophone vertical seismic profiles at Citronelle,Alabama 
下载免费PDF全文
下载免费PDF全文 A modular borehole monitoring concept has been implemented to provide a suite of well‐based monitoring tools that can be deployed cost effectively in a flexible and robust package. The initial modular borehole monitoring system was deployed as part of a CO2 injection test operated by the Southeast Regional Carbon Sequestration Partnership near Citronelle, Alabama. The Citronelle modular monitoring system transmits electrical power and signals, fibre‐optic light pulses, and fluids between the surface and a reservoir. Additionally, a separate multi‐conductor tubing‐encapsulated line was used for borehole geophones, including a specialized clamp for casing clamping with tubing deployment. The deployment of geophones and fibre‐optic cables allowed comparison testing of distributed acoustic sensing. We designed a large source effort (>64 sweeps per source point) to test fibre‐optic vertical seismic profile and acquired data in 2013. The native measurement in the specific distributed acoustic sensing unit used (an iDAS from Silixa Ltd) is described as a localized strain rate. Following a processing flow of adaptive noise reduction and rebalancing the signal to dimensionless strain, improvement from repeated stacking of the source was observed. Conversion of the rebalanced strain signal to equivalent velocity units, via a scaling by local apparent velocity, allows quantitative comparison of distributed acoustic sensing and geophone data in units of velocity. We see a very good match of uncorrelated time series in both amplitude and phase, demonstrating that velocity‐converted distributed acoustic sensing data can be analyzed equivalent to vertical geophones. We show that distributed acoustic sensing data, when averaged over an interval comparable to typical geophone spacing, can obtain signal‐to‐noise ratios of 18 dB to 24 dB below clamped geophones, a result that is variable with noise spectral amplitude because the noise characteristics are not identical. With vertical seismic profile processing, we demonstrate the effectiveness of downgoing deconvolution from the large spatial sampling of distributed acoustic sensing data, along with improved upgoing reflection quality. We conclude that the extra source effort currently needed for tubing‐deployed distributed acoustic sensing vertical seismic profile, as part of a modular monitoring system, is well compensated by the extra spatial sampling and lower deployment cost as compared with conventional borehole geophones. 相似文献
18.
基于扰动理论的混合域波动方程波场外推算子,具有一定介质横向速度变化适应能力,是反射地震学中研究较多的成像方法。此类波场外推算子沿深度层进行波场外推,都存在参考速度选择问题。单参考速度波场外推算子,适应地下介质弱横向变速,而多参考速度波场外推算子可以提高横向变速的适应能力和成像精度,但要以大量计算为代价。本文提出的自适应多参考速度选择策略,根据外推层地质构造的复杂度和给定的速度门槛值自动选择参考速度个数,利用该策略构造混合域SSF、FFD、WXFD和GSP等多参考速度波场外推成像算法。盐丘模型理论数据测试结果表明,自适应多参考速度波场外推算法具有强横向变速适应能力和较高成像精度。 相似文献
19.
Everhard Muyzert Nihed Allouche Pascal Edme Nicolas Goujon 《Geophysical Prospecting》2019,67(1):97-113
We built a five-component (5C) land seismic sensor that measures both the three-component (3C) particle acceleration and two vertical gradients of the horizontal wavefield through a pair of 3C microelectromechanical accelerometers. The sensor is a small cylindrical device planted vertically just below the earth's surface. We show that seismic acquisition and processing 5C sensor data has the potential to replace conventional seismic acquisition with analogue geophone groups by single 5C sensors placed at the same station interval when combined with a suitable aliased ground roll attenuation algorithm. The 5C sensor, therefore, allows for sparser, more efficient, data acquisition. The accuracy of the 5C sensor wavefield gradients depends on the 3C accelerometers, their sensitivity, self-noise and their separation. These sensor component specifications are derived from various modelling studies. The design principles of the 5C sensor are validated using test data from purpose-built prototypes. The final prototype was constructed with a pair of 3C accelerometers separated by 20 cm and with a self-noise of 35 ng Hz−1/2. Results from a two-dimensional seismic line show that the seismic image of 5C sensor data with ground roll attenuated using 5C sensor gradient data was comparable to simulated analogue group data as is the standard in the industry. This field example shows that up to three times aliased ground roll was attenuated. The 5C sensor also allows for correcting vertical component accelerometer data for sensor tilt. It is shown that a vertical component sensor that is misaligned with the vertical direction by 10° introduces an error in the seismic data of around –20 dB with respect to the seismic signal, which can be fully corrected. Advances in sensor specifications and processing algorithms are expected to lead to even more effective ground roll attenuation, enabling a reduction in the receiver density resulting in a smaller number of sensors that must be deployed and, therefore, improving the operational efficiency while maintaining image quality. 相似文献