共查询到20条相似文献，搜索用时 31 毫秒
1.
Xing-Yan Zhang Dong-Ming Pan Wen-Ying Shi Zhong-Yu Fang Zhi-Wei Dan Li-Xia Zhang 《应用地球物理》2015,12(4):573-584
Currently, the deghosting of towed streamer seismic data assumes a flat sea level and a sea-surface reflection coefficient of ?1; this decreases the precision of deghosting. A new method that considers the rough sea surface is proposed to suppress ghost reflections. The proposed deghosting method obtains the rough sea surface reflection coefficient using Gaussian statistics, and calculates the optimized deghosting operator in the τ/p domain. The proposed method is closer to the actual sea conditions, offers an improved deghosting operator, removes the ghost reflections from marine towed seismic data, widens the bandwidth and restores the low-frequency information, and finally improves the signal-tonoise ratio and resolution of the seismic data. 相似文献
2.
Three‐dimensional receiver ghost attenuation (deghosting) of dual‐sensor towed‐streamer data is straightforward, in principle. In its simplest form, it requires applying a three‐dimensional frequency–wavenumber filter to the vertical component of the particle motion data to correct for the amplitude reduction on the vertical component of non‐normal incidence plane waves before combining with the pressure data. More elaborate techniques use three‐dimensional filters to both components before summation, for example, for ghost wavelet dephasing and mitigation of noise of different strengths on the individual components in optimum deghosting. The problem with all these techniques is, of course, that it is usually impossible to transform the data into the crossline wavenumber domain because of aliasing. Hence, usually, a two‐dimensional version of deghosting is applied to the data in the frequency–inline wavenumber domain. We investigate going down the “dimensionality ladder” one more step to a one‐dimensional weighted summation of the records of the collocated sensors to create an approximate deghosting procedure. We specifically consider amplitude‐balancing weights computed via a standard automatic gain control before summation, reminiscent of a diversity stack of the dual‐sensor recordings. This technique is independent of the actual streamer depth and insensitive to variations in the sea‐surface reflection coefficient. The automatic gain control weights serve two purposes: (i) to approximately correct for the geometric amplitude loss of the Z data and (ii) to mitigate noise strength variations on the two components. Here, Z denotes the vertical component of the velocity of particle motion scaled by the seismic impedance of the near‐sensor water volume. The weights are time‐varying and can also be made frequency‐band dependent, adapting better to frequency variations of the noise. The investigated process is a very robust, almost fully hands‐off, approximate three‐dimensional deghosting step for dual‐sensor data, requiring no spatial filtering and no explicit estimates of noise power. We argue that this technique performs well in terms of ghost attenuation (albeit, not exact ghost removal) and balancing the signal‐to‐noise ratio in the output data. For instances where full three‐dimensional receiver deghosting is the final product, the proposed technique is appropriate for efficient quality control of the data acquired and in aiding the parameterisation of the subsequent deghosting processing. 相似文献
3.
作为一种特殊的噪声,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率及保真度的重要因素.以格林公式为基础,详细论述了基于格林函数理论的鬼波压制方法,在不需要地下介质信息的条件下,进行地震数据驱动鬼波压制,并根据"Double Dirichlet"(双狄利克雷)边界条件,预测压力波场和垂直速度波场.建立了基于格林函数理论鬼波压制的处理流程,数值模拟和实际资料处理结果表明,基于格林函数理论鬼波压制方法在很好地去除鬼波的同时极大地拓宽了地震资料的频带,尤其提升了低频端能量,有利于后续资料的处理解释. 相似文献
4.
The rough sea surface causes perturbations in the seismic data that can be significant for time‐lapse studies. The perturbations arise because the reflection response of the non‐flat sea perturbs the seismic wavelet. In order to remove these perturbations from the received seismic data, special deconvolution methods can be used, but these methods require, as input, the time varying wave elevation above each hydrophone in the streamer. In addition, the vertical displacement of the streamer itself must also be known at the position of each hydrophone and at all times. This information is not available in conventional seismic acquisition. However, it can be obtained from the hydrophone measurements provided that the hydrophones are recorded individually (not grouped), that the recording bandwidth is extended down to 0.05 Hz and that data are recorded without gaps between the shot records. The sea surface elevation, and also the wave‐induced vertical displacement of the streamer, can be determined from the time‐varying pressure that the sea waves cause in the hydrophone measurements. When this was done experimentally, using a single sensor seismic streamer without a conventional low cut filter, the wave induced pressure variations were easily detected. The inversion of these experimental data gives results for the sea surface elevation that are consistent with the weather and sea state at the time of acquisition. A high tension approximation allows a simplified solution of the equations that does not demand a knowledge of the streamer tension. However, best results at the tail end of the streamer are obtained using the general equation. 相似文献
5.
Wasiu Makinde Nathalie Favretto-Cristini Eric de Bazelaire 《Geophysical Prospecting》2005,53(3):373-397
Seismic wavefield scattering from a statistically randomly rough interface in a multilayered piecewise homogeneous medium is studied in 3D. The influence of the surface roughness on the scattered wavefield is analysed numerically by using a finite‐difference operator in the acoustic domain. Since interface scattering in the real practical sense is a 3D physical phenomenon, we show in this work that the scattering response of a randomly rough interface is not the same in 3D situations as in the 2D cases described in some earlier works. For a given interface roughness height in 3D, an interface roughness height at least three times greater is required to produce an equivalent phase scattering effect in 2D situations, for a given correlation length of the interface roughness scale. Based on observations from spectral analysis, we show that scattering results principally in de‐phasing and frequency band‐limiting of the incident wavefront, the frequency band‐limiting properties being comparable to cases reported in the literature for absorption and thin‐layer filtering. The interface scattering phenomenon should be critically considered when using amplitude and phase information from seismic signal during inversion processes. 相似文献
6.
Comparing the performances of four stochastic optimisation methods using analytic objective functions, 1D elastic full‐waveform inversion,and residual static computation 
下载免费PDF全文
下载免费PDF全文 Angelo Sajeva Mattia Aleardi Bruno Galuzzi Eusebio Stucchi Emmanuel Spadavecchia Alfredo Mazzotti 《Geophysical Prospecting》2017,65(Z1):322-346
We compare the performances of four stochastic optimisation methods using four analytic objective functions and two highly non‐linear geophysical optimisation problems: one‐dimensional elastic full‐waveform inversion and residual static computation. The four methods we consider, namely, adaptive simulated annealing, genetic algorithm, neighbourhood algorithm, and particle swarm optimisation, are frequently employed for solving geophysical inverse problems. Because geophysical optimisations typically involve many unknown model parameters, we are particularly interested in comparing the performances of these stochastic methods as the number of unknown parameters increases. The four analytic functions we choose simulate common types of objective functions encountered in solving geophysical optimisations: a convex function, two multi‐minima functions that differ in the distribution of minima, and a nearly flat function. Similar to the analytic tests, the two seismic optimisation problems we analyse are characterised by very different objective functions. The first problem is a one‐dimensional elastic full‐waveform inversion, which is strongly ill‐conditioned and exhibits a nearly flat objective function, with a valley of minima extended along the density direction. The second problem is the residual static computation, which is characterised by a multi‐minima objective function produced by the so‐called cycle‐skipping phenomenon. According to the tests on the analytic functions and on the seismic data, genetic algorithm generally displays the best scaling with the number of parameters. It encounters problems only in the case of irregular distribution of minima, that is, when the global minimum is at the border of the search space and a number of important local minima are distant from the global minimum. The adaptive simulated annealing method is often the best‐performing method for low‐dimensional model spaces, but its performance worsens as the number of unknowns increases. The particle swarm optimisation is effective in finding the global minimum in the case of low‐dimensional model spaces with few local minima or in the case of a narrow flat valley. Finally, the neighbourhood algorithm method is competitive with the other methods only for low‐dimensional model spaces; its performance sensibly worsens in the case of multi‐minima objective functions. 相似文献
7.
Three‐dimensional receiver deghosting of seismic streamer data using L1 inversion and redundant extended radon dictionary 下载免费PDF全文
下载免费PDF全文 In this paper, we propose a novel three‐dimensional receiver deghosting algorithm that is capable of deghosting both horizontal and slanted streamer data in a theoretically consistent manner. Our algorithm honours wave propagation phenomena in a true three‐dimensional sense and frames the three‐dimensional receiver deghosting problem as a Lasso problem. The ultimate goal is to minimise the mismatch between the actual measurements and the simulated wavefield with an L1 constraint applied in the extended Radon space to handle the underdetermined nature of this problem. We successfully demonstrate our algorithm on a modified three‐dimensional EAGE/SEG Overthrust model and a Red Sea marine dataset. 相似文献
8.
Denis Mougenot 《Geophysical Prospecting》2018,66(5):922-933
The broadband capabilities of marine, seabed, and land seismic equipment are reviewed with respect to both the source and the receiver sides. In marine acquisition, the main issue at both ends of the spectrum relates to ghosts occurring at the sea surface. Broadband deghosting requires towing at variable depth to introduce notch diversity or using new equipment like multi‐component and/or low‐noise streamers. As a result, a doubling of the bandwidth from about three to six octaves (2.5–200 Hz) has been achieved. Such improvement is not yet observed for seabed surveys in spite of deghosting being a standard process on the receiver side. One issue may be related to the coupling of the particle motion sensor, particularly at high frequencies. For land acquisition, progress came from the vibrators. New shakers and control electronics using broadband sweeps made it possible to add two more octaves to the low‐frequency signal (from 8 to 2 Hz). Whereas conventional 10 Hz geophones are still able to record such low frequencies, 5 Hz high gain geophones or digital accelerometers enhance them to keep the signal above the noise floor. On the high end of the bandwidth, progress is not limited by equipment specifications. Here, the issue is related to a low signal‐to‐noise ratio due to the strong absorption that occurs during signal propagation. To succeed in enlarging the bandwidth, these improved equipment and sweeps must be complemented by a denser spatial sampling of the wavefield by point–source and point–receiver acquisition. 相似文献
9.
Single‐component towed‐streamer marine data acquisition records the pressure variations of the upgoing compressional waves followed by the polarity‐reversed pressure variations of downgoing waves, creating sea‐surface ghost events in the data. The sea‐surface ghost for constant‐depth towed‐streamer marine data acquisition is usually characterised by a ghost operator acting on the upgoing waves, which can be formulated as a filtering process in the frequency–wavenumber domain. The deghosting operation, usually via the application of the inverse Wiener filter related to the ghost operator, acts on the signal as well as the noise. The noise power transfer into the deghosted data is proportional to the power spectrum of the inverse Wiener filter and is amplifying the noise strongly at the notch wavenumbers and frequencies of the ghost operator. For variable‐depth streamer acquisition, the sea‐surface ghost cannot be described any longer as a wavenumber–frequency operator but as a linear relationship between the wavenumber–frequency representation of the upgoing waves at the sea surface and the data in the space–frequency domain. In this article, we investigate how the application of the inverse process acts on noise. It turns out that the noise magnification is less severe with variable‐depth streamer data, as opposed to constant depth, and is inversely proportional to the local slant of the streamer. We support this statement via application of the deghosting process to real and numerical random noise. We also propose a more general concept of a wavenumber–frequency ghost power transfer function, applicable for variable‐depth streamer acquisition, and demonstrate that the inverse of the proposed variable‐depth ghost power transfer function can be used to approximately quantify the action of the variable‐depth streamer deghosting process on noise. 相似文献
10.
11.
In marine acquisition, reflections of sound energy from the water–air interface result in ghosts in the seismic data, both in the source side and the receiver side. Ghosts limit the bandwidth of the useful signal and blur the final image. The process to separate the ghost and primary signals, called the deghosting process, can fill the ghost notch, broaden the frequency band, and help achieve high‐resolution images. Low‐signal‐to‐noise ratio near the notch frequencies and 3D effects are two challenges that the deghosting process has to face. In this paper, starting from an introduction to the deghosting process, we present and compare two strategies to solve the latter. The first is an adaptive mechanism that adjusts the deghosting operator to compensate for 3D effects or errors in source/receiver depth measurement. This method does not include explicitly the crossline slowness component and is not affected by the sparse sampling in the same direction. The second method is an inversion‐type approach that does include the crossline slowness component in the algorithm and handles the 3D effects explicitly. Both synthetic and field data examples in wide azimuth acquisition settings are shown to compare the two strategies. Both methods provide satisfactory results. 相似文献
12.
鬼波问题是影响海上地震资料分辨率和保真度提高的最重要因素之一.详细论述了逆散射级数理论和逆散射级数法鬼波压制原理,说明了逆散射级数方法进行鬼波压制理论的完善性和对鬼波描述的精确性.实现了基于逆散射级数理论的鬼波压制方法,方法以波动方程和Lippmann-Schwinger方程为基础,在频率-波数-波数域内构造与鬼波相关的压制算子,在不需要对地下介质作任何假设条件下实现地震数据驱动鬼波压制,并通过改善消除鬼波的压制算子,提高算法的稳定性.资料处理试验与处理结果分析表明,基于逆散射级数鬼波压制方法能在实现鬼波压制的同时较好保留有效反射波的信息,从而补偿地震资料低频损失和提高地震数据的保真度.数据处理试验还表明,研究方法能对低信噪比的地震资料进行有效的鬼波压制处理.建立了基于逆散射级数鬼波压制处理流程. 相似文献
13.
Time‐lapse seismic surveying has become an accepted tool for reservoir monitoring applications, thus placing a high premium on data repeatability. One factor affecting data repeatability is the influence of the rough sea‐surface on the ghost reflection and the resulting seismic wavelets of the sources and receivers. During data analysis, the sea‐surface is normally assumed to be stationary and, indeed, to be flat. The non‐flatness of the sea‐surface introduces amplitude and phase perturbations to the source and receiver responses and these can affect the time‐lapse image. We simulated the influence of rough sea‐surfaces on seismic data acquisition. For a typical seismic line with a 48‐fold stack, a 2‐m significant‐wave‐height sea introduces RMS errors of about 5–10% into the stacked data. This level of error is probably not important for structural imaging but could be significant for time‐lapse surveying when the expected difference anomaly is small. The errors are distributed differently for sources and receivers because of the different ways they are towed. Furthermore, the source wavelet is determined by the sea shape at the moment the shot is fired, whereas the receiver wavelet is time‐varying because the sea moves significantly during the seismic record. 相似文献
14.
高分辨率的宽带地震勘探技术是最近几年海上油气地震勘探的热点问题,鬼波压制是其中的核心议题.石油工业界,除了提出各种压制鬼波的采集方式外,资料处理过程中压制鬼波的方法也在不断地推陈出新.源、检鬼波的存在使得反射子波的有效频带变窄,成像分辨率降低,且干扰后续的自由表面多次波压制、FWI等.不同于常规的鬼波压制模型,本文基于编码与解码理论框架,用编码建立起鬼波预测模型,在Bayes反演框架下建立起解码估计一次波、从而压制鬼波的方法.基于此鬼波预测模型及相应反演理论的鬼波压制方法对鬼波的压制更为彻底.理论模型数据与实际资料测试结果验证了本文提出的理论框架和方法技术的有效性与优越性. 相似文献
15.
The effects of the topographic data source and resolution on the hydraulic modelling of floods were analysed. Seven digital terrain models (DTMs) were generated from three different altimetric sources: a global positioning system (GPS) survey and bathymetry; high‐resolution laser altimetry data LiDAR (light detection and ranging); and vectorial cartography (1:5000). Hydraulic results were obtained, using the HEC‐RAS one‐dimensional model, for all seven DTMs. The importance of the DTM's accuracy on the hydraulic modelling results was analysed within three different hydraulic contexts: (1) the discharge and water surface elevation results from the hydraulic model; (2) the delineation of the flooded area; and (3) the relative sensitivity of the hydraulic model to changes in the Manning's n roughness coefficient. The contour‐based DTM was the least accurate with a root mean square error (RMSE) of 4·5 m in the determination of the water level and a variation of up to 50 per cent in the estimation of the inundated area of the floodplain. The GPS‐based DTM produced more realistic water surface elevation results and variations of up to 8 per cent in terms of the flooded area. The laser‐based model's RMSE for water level was 0·3 m, with the flooded area varying by less than 1 per cent. The LiDAR data also showed the greatest sensitivity to changes in the Manning's roughness coefficient. An analysis of the effect of mesh resolution indicated an influence on the delineation of the flooded area with variations of up to 7·3 per cent. In addition to determining the accuracy of the hydraulic modelling results produced from each DTM, an analysis of the time–cost ratio of each topographic data source illustrates that airborne laser scanning is a cost‐effective means of developing a DTM of sufficient accuracy, especially over large areas. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
16.
An adaptive free‐surface expression for three‐dimensional finite‐difference frequency‐domain modelling of elastic wave 下载免费PDF全文
下载免费PDF全文 Finite‐difference frequency‐domain modelling of seismic wave propagation is attractive for its efficient solution of multisource problems, and this is crucial for full‐waveform inversion and seismic imaging, especially in the three‐dimensional seismic problem. However, implementing the free surface in the finite‐difference method is nontrivial. Based on an average medium method and the limit theorem, we present an adaptive free‐surface expression to describe the behaviour of wavefields at the free surface, and no extra work for the free‐surface boundary condition is needed. Essentially, the proposed free‐surface expression is a modification of density and constitutive relation at the free surface. In comparison with a direct difference approximate method of the free‐surface boundary condition, this adaptive free‐surface expression can produce more accurate and stable results for a broad range of Poisson's ratio. In addition, this expression has a good performance in handling the lateral variation of Poisson's ratio adaptively and without instability. 相似文献
17.
A marine source generates both a direct wavefield and a ghost wavefield. This is caused by the strong surface reflectivity, resulting in a blended source array, the blending process being natural. The two unblended response wavefields correspond to the real source at the actual location below the water level and to the ghost source at the mirrored location above the water level. As a consequence, deghosting becomes deblending (‘echo‐deblending’) and can be carried out with a deblending algorithm. In this paper we present source deghosting by an iterative deblending algorithm that properly includes the angle dependence of the ghost: It represents a closed‐loop, non‐causal solution. The proposed echo‐deblending algorithm is also applied to the detector deghosting problem. The detector cable may be slanted, and shot records may be generated by blended source arrays, the blending being created by simultaneous sources. Similar to surface‐related multiple elimination the method is independent of the complexity of the subsurface; only what happens at and near the surface is relevant. This means that the actual sea state may cause the reflection coefficient to become frequency dependent, and the water velocity may not be constant due to temporal and lateral variations in the pressure, temperature, and salinity. As a consequence, we propose that estimation of the actual ghost model should be part of the echo‐deblending algorithm. This is particularly true for source deghosting, where interaction of the source wavefield with the surface may be far from linear. The echo‐deblending theory also shows how multi‐level source acquisition and multi‐level streamer acquisition can be numerically simulated from standard acquisition data. The simulated multi‐level measurements increase the performance of the echo‐deblending process. The output of the echo‐deblending algorithm on the source side consists of two ghost‐free records: one generated by the real source at the actual location below the water level and one generated by the ghost source at the mirrored location above the water level. If we apply our algorithm at the detector side as well, we end up with four ghost‐free shot records. All these records are input to migration. Finally, we demonstrate that the proposed echo‐deblending algorithm is robust for background noise. 相似文献
18.
Nicolas Baghdadi Olivier Cerdan Mehrez Zribi Véronique Auzet Frédéric Darboux Mahmoud El Hajj Rania Bou Kheir 《水文研究》2008,22(1):9-20
Synthetic aperture radar (SAR) sensors are often used to characterize the surface of bare soils in agricultural environments. They enable the soil moisture and roughness to be estimated with constraints linked to the configurations of the sensors (polarization, incidence angle and radar wavelength). These key soil characteristics are necessary for different applications, such as hydrology and risk prediction. This article reviews the potential of currently operational SAR sensors and those planned for the near future to characterize soil surface as a function of users' needs. It details what it is possible to achieve in terms of mapping soil moisture and roughness by specifying optimal radar configurations and the precision associated with the estimation of soil surface characteristics. The summary carried out for the present article shows that mapping soil moisture is optimal with SAR sensors at low incidence angles (<35 ). This configuration, which enables an estimated moisture accuracy greater than 6% is possible several times a month taking into account all the current and future sensors. Concerning soil roughness, it is best mapped using three classes (smooth, moderately rough, and rough). Such mapping requires high‐incidence data, which is possible with certain current sensors (RADARSAT‐1 and ASAR both in band C). When L‐band sensors (ALOS) become available, this mapping accuracy should improve because the sensitivity of the radar signal to Soil Surface Characteristics (SSC) increases with wavelength. Finally, the polarimetric mode of certain imminent sensors (ALOS, RADARSAT‐2, TerraSAR‐X, etc.), and the possibility of acquiring data at very high spatial resolution (metre scale), offer great potential in terms of improving the quality of SSC mapping. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
19.
Hakan Karslı 《Geophysical Prospecting》2011,59(1):56-65
Interpreting a post‐stack seismic section is difficult due to the band‐limited nature of the seismic data even post deconvolution. Deconvolution is a process that is universally applied to extend the bandwidth of seismic data. However, deconvolution falls short of this task as low and high frequencies of the deconvolved data are either still missing or contaminated by noise. In this paper we use the autoregressive extrapolation technique to recover these missing frequencies, using the high signal‐to‐noise ratio (S/N) portions of the spectrum of deconvolved data. I introduce here an algorithm to extend the bandwidth of deconvolved data. This is achieved via an autoregressive extrapolation technique, which has been widely used to replace missing or corrupted samples of data in signal processing. This method is performed in the spectral domain. The spectral band to be extrapolated using autoregressive prediction filters is first selected from the part of the spectrum that has a high signal‐to‐noise ratio (S/N) and is then extended. As there can be more than one zone of good S/N in the spectrum, the results of prediction filter design and extrapolation from three different bands are averaged. When the spectrum of deconvolved data is extended in this way, the results show higher vertical resolution to a degree that the final seismic data closely resemble what is considered to be a reflectivity sequence of the layered medium. This helps to obtain acoustic impedance with inversion by stable integration. The results show that autoregressive spectral extrapolation highly increases vertical resolution and improves horizon tracking to determine continuities and faults. This increase in coherence ultimately yields a more interpretable seismic section. 相似文献
20.
The classical finite‐difference methods for seismic wave modelling are very accurate at low wavenumbers but suffer from inaccuracies at high wavenumbers, particularly at Nyquist wavenumber. In contrast, the optimisation finite‐difference methods reduce inaccuracies at high wavenumbers but suffer from inaccuracies at low wavenumbers, particularly at zero wavenumber when the operator length is not long and the whole range of wavenumbers is considered. Inaccuracy at zero wavenumber means that the optimisation methods only have a zeroth‐order accuracy of truncation and thus are not rigorously convergent. To guarantee the rigorous convergence of the optimisation methods, we have developed accuracy‐constrained optimisation methods. Different‐order accuracy‐constrained optimisation methods are presented. These methods not only guarantee the rigorous convergence but also reduce inaccuracies at low wavenumbers. Accuracy‐constrained optimisation methods are applied to staggered‐grid elastic wave modelling. 相似文献