共查询到20条相似文献，搜索用时 187 毫秒
1.
Reconstruction of the near‐offset gap in marine seismic data using seismic interferometric interpolation 
下载免费PDF全文
下载免费PDF全文 In conventional seismic exploration, especially in marine seismic exploration, shot gathers with missing near‐offset traces are common. Interferometric interpolation methods are one of a range of different methods that have been developed to solve this problem. Interferometric interpolation methods differ from conventional interpolation methods as they utilise information from multiples in the interpolation process. In this study, we apply both conventional interferometric interpolation (shot domain) and multi‐domain interferometric interpolation (shot and receiver domain) to a synthetic and a real‐towed marine dataset from the Baltic Sea with the primary aim of improving the image of the seabed by extrapolation of a near‐offset gap. We utilise a matching filter after interferometric interpolation to partially mitigate artefacts and coherent noise associated with the far‐field approximation and a limited recording aperture size. The results show that an improved image of the seabed is obtained after performing interferometric interpolation. In most cases, the results from multi‐domain interferometric interpolation are similar to those from conventional interferometric interpolation. However, when the source–receiver aperture is limited, the multi‐domain method performs better. A quantitative analysis for assessing the performance of interferometric interpolation shows that multi‐domain interferometric interpolation typically performs better than conventional interferometric interpolation. We also benchmark the interpolated results generated by interferometric interpolation against those obtained using sparse recovery interpolation. 相似文献
2.
Seismic interferometry is a relatively new technique to estimate the Green's function between receivers. Spurious energy, not part of the true Green's function, is produced because assumptions are commonly violated when applying seismic interferometry to field data. Instead of attempting to suppress all spurious energy, we show how spurious energy associated with refractions contains information about the subsurface in field data collected at the Boise Hydrogeophysical Research Site. By forming a virtual shot record we suppress uncorrelated noise and produce a virtual refraction that intercepts zero offset at zero time. These two features make the virtual refraction easy to pick, providing an estimate of refractor velocity. To obtain the physical parameters of the layer above the refractor we analyse the cross‐correlation of wavefields recorded at two receivers for all sources. A stationary‐phase point associated with the correlation between the reflected wave and refracted wave from the interface identifies the critical offset. By combining information from the virtual shot record, the correlation gather and the real shot record we determine the seismic velocities of the unsaturated and saturated sands, as well as the variable relative depth to the water‐table. Finally, we discuss how this method can be extended to more complex geologic models. 相似文献
3.
Fabio Ciabarri Alfredo Mazzotti Eusebio Stucchi Nicola Bienati 《Geophysical Prospecting》2015,63(2):296-314
Least squares Fourier reconstruction is basically a solution to a discrete linear inverse problem that attempts to recover the Fourier spectrum of the seismic wavefield from irregularly sampled data along the spatial coordinates. The estimated Fourier coefficients are then used to reconstruct the data in a regular grid via a standard inverse Fourier transform (inverse discrete Fourier transform or inverse fast Fourier transform). Unfortunately, this kind of inverse problem is usually under‐determined and ill‐conditioned. For this reason, the least squares Fourier reconstruction with minimum norm adopts a damped least squares inversion to retrieve a unique and stable solution. In this work, we show how the damping can introduce artefacts on the reconstructed 3D data. To quantitatively describe this issue, we introduce the concept of “extended” model resolution matrix, and we formulate the reconstruction problem as an appraisal problem. Through the simultaneous analysis of the extended model resolution matrix and of the noise term, we discuss the limits of the Fourier reconstruction with minimum norm reconstruction and assess the validity of the reconstructed data and the possible bias introduced by the inversion process. Also, we can guide the parameterization of the forward problem to minimize the occurrence of unwanted artefacts. A simple synthetic example and real data from a 3D marine common shot gather are used to discuss our approach and to show the results of Fourier reconstruction with minimum norm reconstruction. 相似文献
4.
Joint up/down decomposition and reconstruction using three‐component streamers with or without ghost model: the sampling theory 下载免费PDF全文
下载免费PDF全文 This paper addresses two artefacts inherent to marine towed‐streamer surveys: 1) ghost reflections and 2) too sparse a sampling in the crossline direction. A ghost reflection is generated when an upcoming reflection bounces off the sea surface back into the sensors and can, in principle, be removed by decomposing the measured wavefield into its up‐ and downgoing constituents. This process requires a dense sampling of the wavefield in both directions along and perpendicular to the streamers. A dense sampling in the latter direction is, however, often impossible due to economical and operational constraints. Recent multi‐component streamers have been designed to record the spatial gradients on top of the pressure, which not only benefits the wavefield decomposition but also facilitates a lower‐than‐Nyquist sampling rate of the pressure. In this paper, wavefield reconstruction and deghosting are posed as a joint inverse problem. We present two approaches to establish a system matrix that embeds both a deghosting and an interpolation operator. The first approach is derived with a ghost model, whereas the second approach is derived without a ghost model. The embodiment of a ghost model leads to an even lower sampling rate but relies on a more restrictive assumption on the sea surface. 相似文献
5.
三维表面多次波压制是海洋地震资料预处理中的重要研究课题,基于波动理论的三维表面多次波压制方法(3DSRME)是数据驱动的方法,理论上来说,可有效压制复杂构造地震数据表面多次波.但该方法因对原始地震数据采集要求高而很难在实际资料处理中广泛应用.本文基于贡献道集的概念,将稀疏反演方法引入到表面多次波压制中,应用稀疏反演代替横测线积分求和,无需对横测线进行大规模重建,进而完成三维表面多次波预测,这样可有效解决实际三维地震数据横测线方向稀疏的问题.基于纵测线多次波积分道集为抛物线的假设,为保证预测后三维表面多次波和全三维数据预测的多次波在运动学和动力学特征上基本一致,文中对预测数据实施基于稳相原理的相位校正.理论模型和实际数据的测试结果表明,本文基于稀疏反演三维表面多次波压制方法可在横测线稀疏的情况下,有效压制三维复杂介质地震资料中的表面多次波,从而更好地提高海洋地震资料的信噪比,为高分辨率地震成像提供可靠的预处理数据保障. 相似文献
6.
Passive seismic has recently attracted a great deal of attention because non‐artificial source is used in subsurface imaging. The utilization of passive source is low cost compared with artificial‐source exploration. In general, constructing virtual shot gathers by using cross‐correlation is a preliminary step in passive seismic data processing, which provides the basis for applying conventional seismic processing methods. However, the subsurface structure is not uniformly illuminated by passive sources, which leads to that the ray path of passive seismic does not fit the hyperbolic hypothesis. Thereby, travel time is incorrect in the virtual shot gathers. Besides, the cross‐correlation results are contaminated by incoherent noise since the passive sources are always natural. Such noise is kinematically similar to seismic events and challenging to be attenuated, which will inevitably reduce the accuracy in the subsequent process. Although primary estimation for transient‐source seismic data has already been proposed, it is not feasible to noise‐source seismic data due to the incoherent noise. To overcome the above problems, we proposed to combine focal transform and local similarity into a highly integrated operator and then added it into the closed‐loop surface‐related multiple elimination based on the 3D L1‐norm sparse inversion framework. Results proved that the method was capable of reliably estimating noise‐free primaries and correcting travel time at far offsets for a foresaid virtual shot gathers in a simultaneous closed‐loop inversion manner. 相似文献
7.
Least-squares reverse time migration is often formulated as an iterative updating process, where estimating the gradient of the misfit function is necessary. Traditional time domain shot-profile least-squares reverse time migration is computationally expensive because computing the gradient involves solving the two-way wave equation several times in every iteration. To reduce the computational cost of least-squares reverse time migration, we propose a double-plane-wave least-squares reverse time migration method based on a misfit function for frequency-domain double-plane-wave data. In double-plane-wave least-squares reverse time migration, the gradient is computed by multiplying frequency-domain plane-wave Green's functions with the corresponding double-plane-wave data residual. Because the number of plane-wave Green's functions used for migration is relatively small, they can be pre-computed and stored in a computer's discs or memory. We can use the pre-computed plane-wave Green's functions to obtain the gradient without solving the two-way wave equation in each iteration. Therefore, the migration efficiency is significantly improved. In addition, we study the effects of using sparse frequency sampling and sparse plane-wave sampling on the proposed method. We can achieve images with correct reflector amplitudes and reasonable resolution using relatively sparse frequency sampling and plane-wave sampling, which are larger than that determined by the Nyquist theorem. The well-known wrap-around artefacts and linear artefacts generated due to under-sampling frequency and plane wave can be suppressed during iterations in cases where the sampling rates are not excessively large. Moreover, implementing the proposed method with sparse frequency sampling and sparse plane-wave sampling further improves the computational efficiency. We test the proposed double-plane-wave least-squares reverse time migration on synthetic models to show the practicality of the method. 相似文献
8.
We propose a new method for removing sea-surface multiples from marine seismic reflection data in which, in essence, the reflection response of the earth, referred to a plane just above the sea-floor, is computed as the ratio of the plane-wave components of the upgoing wave and the downgoing wave. Using source measurements of the wavefield made during data acquisition, three problems associated with earlier work are solved: (i) the method accommodates source arrays, rather than point sources; (ii) the incident field is removed without simultaneously removing part of the scattered field; and (iii) the minimum-energy criterion to find a wavelet is eliminated. Pressure measurements are made in a horizontal plane in the water. The source can be a conventional array of airguns, but must have both in-line and cross-line symmetry, and its wavefield must be measured and be repeatable from shot to shot. The problem is formulated for multiple shots in a two-dimensional configuration for each receiver, and for multiple receivers in a two-dimensional configuration for each shot. The scattered field is obtained from the measurements by subtracting the incident field, known from measurements at the source. The scattered field response to a single incident plane wave at a single receiver is obtained by transforming the common-receiver gather to the frequency–wavenumber domain, and a single component of this response is obtained by Fourier transforming over all receiver coordinates. Each scattered field component is separated into an upgoing wave and a downgoing wave using the zero-pressure condition at the water-surface. The upgoing wave may then be expressed as a reflection coefficient multiplied by the incident downgoing wave plus a sum of scattered downgoing plane waves, each multiplied by the corresponding reflection coefficient. Keeping the upgoing scattered wave fixed, and using all possible incident plane waves for a given frequency, yields a set of linear simultaneous equations for the reflection coefficients which are solved for each plane wave and for each frequency. To create the shot records that would have been measured if the sea-surface had been absent, each reflection coefficient is multiplied by complex amplitude and phase factors, for source and receiver terms, before the five-dimensional Fourier transformation back to the space–time domain. 相似文献
9.
利用分形插值方法对稀疏采样的地震记录进行了插值加密处理,并对原始记录和分形插值重建记录进行了叠前逆时偏移处理．数值实验表明,分形插值重建的地震图是原始地震图的良好近似,而分形插值叠前偏移剖面是原始记录偏移剖面的良好近似,二者的分辨率几乎一样．利用分形插值重建技术可以在节约经济成本的情况下,得到地下结构的高精度成像结果． 相似文献
10.
In order to account for the effects of elastic wave propagation in marine seismic data, we develop a waveform inversion algorithm for acoustic‐elastic media based on a frequency‐domain finite‐element modelling technique. In our algorithm we minimize residuals using the conjugate gradient method, which back‐propagates the errors using reverse time migration without directly computing the partial derivative wavefields. Unlike a purely acoustic or purely elastic inversion algorithm, the Green's function matrix for our acoustic‐elastic algorithm is asymmetric. We are nonetheless able to achieve computational efficiency using modern numerical methods. Numerical examples show that our coupled inversion algorithm produces better velocity models than a purely acoustic inversion algorithm in a wide variety of cases, including both single‐ and multi‐component data and low‐cut filtered data. We also show that our algorithm performs at least equally well on real field data gathered in the Korean continental shelf. 相似文献
11.
Jing Sun Sigmund Slang Thomas Elboth Thomas Larsen Greiner Steven McDonald Leiv-J. Gelius 《Geophysical Prospecting》2020,68(3):845-871
Marine seismic interference noise occurs when energy from nearby marine seismic source vessels is recorded during a seismic survey. Such noise tends to be well preserved over large distances and causes coherent artefacts in the recorded data. Over the years, the industry has developed various denoising techniques for seismic interference removal, but although well performing, they are still time-consuming in use. Machine-learning-based processing represents an alternative approach, which may significantly improve the computational efficiency. In the case of conventional images, autoencoders are frequently employed for denoising purposes. However, due to the special characteristics of seismic data as well as the noise, autoencoders failed in the case of marine seismic interference noise. We, therefore, propose the use of a customized U-Net design with element-wise summation as part of the skip-connection blocks to handle the vanishing gradient problem and to ensure information fusion between high- and low-level features. To secure a realistic study, only seismic field data were employed, including 25,000 training examples. The customized U-Net was found to perform well, leaving only minor residuals, except for the case when seismic interference noise comes from the side. We further demonstrate that such noise can be treated by slightly increasing the depth of our network. Although our customized U-Net does not outperform a standard commercial algorithm in quality, it can (after proper training) read and process one single shot gather in approximately 0.02 s. This is significantly faster than any existing industry denoising algorithm. In addition, the proposed network processes shot gathers in a sequential order, which is an advantage compared with industry algorithms that typically require a multi-shot input to break the coherency of the noise. 相似文献
12.
We present a new depth‐imaging method for seismic data in heterogeneous anisotropic media. This recursive explicit method uses a non‐stationary extrapolation operator to allow lateral velocity variation, and it uses the relationship between phase angle and the spectral coordinates of seismic data to allow velocity variation with phase angle. A qualitative comparison of migration impulse responses suggests that, for an equivalent cost, the symmetric non‐stationary phase‐shift (SNPS) operator is superior to the phase‐shift plus interpolation (PSPI) operator, for very large depth intervals. To demonstrate the potential of the new method, seismic data from a physical model acquired over a transversely isotropic medium are imaged using a shot‐record migration based on the SNPS operator. 相似文献
13.
海水面的虚反射(鬼波)引起海上拖缆采集数据陷波,导致地震记录频带变窄,而近年发展的变深度缆采集技术,具有多样的陷波特征,通过专门的去虚反射处理方法可获得宽频数据.本文基于已有研究成果,将最小二乘反演迭代压制虚反射算法应用于某海上变深度缆宽频处理.基于频率波数域镜像记录生成方法获得镜像炮集记录,并采用最小二乘解从变深度缆原始和镜像炮集记录中提取上行波.针对镜像炮集记录生成受初始速度模型精度的影响,使得某深度缆接收的上行波和下行波之间的实际延迟时间存在误差,采用最小二乘反演迭代算法最优化计算下行波与上行波之间的平均延迟时间和上行波记录,并采用时空数据窗口滑动克服延迟时间随炮检距和目的层深度变化问题.合成数据及某海上实际变深度缆数据处理测试结果表明,该方法能较好地压制变深度缆由海水面产生的虚反射,能达到拓宽地震记录频带目的. 相似文献
14.
A calendar time interpolation method for 2D seismic amplitude maps, done in two steps, is presented. The contour interpolation part is formulated as a quadratic programming problem, whereas the amplitude value interpolation is based on a conditional probability formulation. The method is applied on field data from the Sleipner CO2 storage project. The output is a continuous image (movie) of the CO2 plume. Besides visualization, the output can be used to better couple 4D seismic to other types of data acquired. The interpolation uncertainty increases with the time gap between consecutive seismic surveys and is estimated by leaving a survey out (blind test). Errors from such tests can be used to identify problems in understanding the flow and possibly improve the interpolation scheme for a given case. Field‐life cost of various acquisition systems and repeat frequencies are linked to the time‐lapse interpolation errors. The error in interpolated amplitudes increased by 3%‐4% per year of interpolation gap for the Sleipner case. Interpolation can never fully replace measurements. 相似文献
15.
Kurt Eggenberger Philip Christie Massimiliano Vassallo Ali Özbek Everhard Muyzert Dirk‐Jan van Manen Ed Kragh 《Geophysical Prospecting》2014,62(5):994-1008
Wave field reconstruction – the estimation of a three‐dimensional (3D) wave field representing upgoing, downgoing or the combined total pressure at an arbitrary point within a marine streamer array – is enabled by simultaneous measurements of the crossline and vertical components of particle acceleration in addition to pressure in a multicomponent marine streamer. We examine a repeated sail line of North Sea data acquired by a prototype multicomponent towed‐streamer array for both wave field reconstruction fidelity (or accuracy) and reconstruction repeatability. Data from six cables, finely sampled in‐line but spaced at 75 m crossline, are reconstructed and placed on a rectangular data grid uniformly spaced at 6.25 m in‐line and crossline. Benchmarks are generated using recorded pressure data and compared with wave fields reconstructed from pressure alone, and from combinations of pressure, crossline acceleration and vertical acceleration. We find that reconstruction using pressure and both crossline and vertical acceleration has excellent fidelity, recapturing highly aliased diffractions that are lost by interpolation of pressure‐only data. We model wave field reconstruction error as a linear function of distance from the nearest physical sensor and find, for this data set with some mismatched shot positions, that the reconstructed wave field error sensitivity to sensor mispositioning is one‐third that of the recorded wave field sensitivity. Multicomponent reconstruction is also more repeatable, outperforming single‐component reconstruction in which wave field mismatch correlates with geometry mismatch. We find that adequate repeatability may mask poor reconstruction fidelity and that aliased reconstructions will repeat if the survey geometry repeats. Although the multicomponent 3D data have only 500 m in‐line aperture, limiting the attenuation of non‐repeating multiples, the level of repeatability achieved is extremely encouraging compared to full‐aperture, pressure‐only, time‐lapse data sets at an equivalent stage of processing. 相似文献
16.
Wave‐equation based shot‐record migration provides accurate images but is computationally expensive because every shot must be migrated separately. Shot‐encoding migration, such as random shot‐encoding or plane‐wave migration, aims to reduce the computational cost of the imaging process by combining the original data into synthesized common‐source gathers. Random shot‐encoding migration and plane‐wave migration have different and complementary features: the first recovers the full spatial bandwidth of the image but introduces strong artefacts, which are due to the interference between the different shot wavefields; the second provides an image with limited spatial detail but is free of crosstalk noise. We design a hybrid scheme that combines linear and random shot‐encoding in order to limit the drawbacks and merge the advantages of these two techniques. We advocate mixed shot‐encoding migration through dithering of plane waves. This approach reduces the crosstalk noise relative to random shot‐encoding migration and increases the spatial bandwidth relative to conventional plane‐wave migration when the take‐off angle is limited to reduce the duration of the plane‐wave gather. In turn, this decreases the migration cost. Migration with dithered plane waves operates as a hybrid encoding scheme in‐between the end members represented by plane‐wave migration and random shot‐encoding. Migration with dithered plane waves has several advantages: every synthesized common‐source gather images in a larger aperture, the crosstalk noise is limited and higher spatial resolution is achievable compared to shot‐record migration, random shot‐encoding and linear shot‐encoding, respectively. Computational cost is also reduced relative to both random and linear shot‐encoding migration since fewer synthesized common‐source gathers are necessary to obtain a high signal‐to‐noise ratio and high spatial resolution in the final image. 相似文献
17.
The utilization of the double focal transformation for sparse data representation and data reconstruction 下载免费PDF全文
下载免费PDF全文 In many cases, seismic measurements are coarsely sampled in at least one dimension. This leads to aliasing artefacts and therefore to problems in the subsequent processing steps. To avoid this, seismic data reconstruction can be applied in advance. The success and reliability of reconstruction methods are dependent on the assumptions they make on the data. In many cases, wavefields are assumed to (locally) have a linear space–time behaviour. However, field data are usually complex, with strongly curved events. Therefore, in this paper, we propose the double focal transformation as an efficient way for complex data reconstruction. Hereby, wavefield propagation is formulated as a transformation, where one‐way propagation operators are used as its basis functions. These wavefield operators can be based on a macro velocity model, which allows our method to use prior information in order to make the data decomposition more effective. The basic principle of the double focal transformation is to focus seismic energy along source and receiver coordinates simultaneously. The seismic data are represented by a number of localized events in the focal domain, whereas aliasing noise spreads out. By imposing a sparse solution in the focal domain, aliasing noise is suppressed, and data reconstruction beyond aliasing is achieved. To facilitate the process, only a few effective depth levels need to be included, preferably along the major boundaries in the data, from which the propagation operators can be calculated. Results on 2D and 3D synthetic data illustrate the method's virtues. Furthermore, seismic data reconstruction on a 2D field dataset with gaps and aliased source spacing demonstrates the strength of the double focal transformation, particularly for near‐offset reflections with strong curvature and for diffractions. 相似文献
18.
人工地震方法由于受到野外观测系统和经济因素等的限制,采集的数据在空间方向总是不规则分布.但是,许多地震数据处理技术的应用（如：多次波衰减,偏移和时移地震）都基于空间规则分布条件下的地震数据体.因此,数据插值技术是地震数据处理流程中关键环节之一.失败的插值方法往往会引入虚假信息,给后续处理环节带来严重的影响.迭代插值方法是目前广泛应用的地震数据重建思路,但是常规的迭代插值方法往往很难保证插值精度,并且迭代收敛速度较慢,尤其存在随机噪声的情况下,插值地震道与原始地震道之间存在较大的信噪比差异.因此开发快速的、有效的迭代数据插值方法具有重要的工业价值.本文将地震数据插值归纳为数学基追踪问题,在压缩感知理论框架下,提出新的非线性Bregman整形迭代算法来求解约束最小化问题,同时在迭代过程中提出两种匹配的迭代控制准则,通过有效的稀疏变换对缺失数据进行重建.通过理论模型和实际数据测试本文方法,并且与常规迭代插值算法进行比较,结果表明Bregman整形迭代插值方法能够更加有效地恢复含有随机噪声的缺失地震信息. 相似文献
19.
A local plane-wave approach of generalized diffraction tomography in heterogeneous backgrounds, equivalent to Kirchhoff summation techniques when applied in seismic reflection, is re-programmed to act as repeated synthetic aperture radar (SAR) imaging for seismic prestack depth migration. Spotlight-mode SAR imaging quickly provides good images of the electromagnetic reflectivity of the ground via fast Fourier transform (FFT)-based signal processing. By calculating only the Green's functions connecting the aircraft to the centre of the illuminated patch, scattering structures around that centre are also recovered. SAR technology requires us to examine seismic imaging from the local point of view, where the quantity and quality of the available information at each image point are what are important, regardless of the survey geometry. When adapted to seismics, a local image of arbitrary size and sampling is obtained by FFT of seismic energy maps in the scattering wavenumber domain around each node of a pre-calculated grid of Green's functions. These local images can be used to generate a classic prestack depth-migrated section by collecting only their centres. However, the local images also provide valuable information around the centre, as in SAR. They can therefore help to pre-analyse prestack depth migration efficiently, and to perform velocity analysis at a very low cost. The FFT-based signal-processing approach allows local, efficient and automatic control of anti-aliasing, noise and resolution, including optimized Jacobian weights. Repeated local imaging could also be used to speed up migration, with interpolation between local images associated with a coarse grid of Green's functions, as an alternative to interpolation of Green's functions. The local images may, however, show distortions due to the local plane-wave approximation, and the velocity variations across their frame. Such effects, which are not necessarily a problem in SAR, should be controlled and corrected to further enhance seismic imaging. Applications to realistic models and to real data show that, despite the distortion effects, the local images can yield similar information to prestack depth migration, including common-image-point gathers for velocity analyses and AVO/AVA effects, at a much lower cost when a small target is considered. 相似文献
20.
常规虚源点Marchenko自聚焦多次波预测方法只适用于预测不含自由表面的多次波模型,局限于压制层间多次波,该方法在构建上下行格林函数场前,必须从反射响应中去除所有与表面相关的多次波.本文对构建上下行Marchenko格林函数方程进行改进,得到了包含一次波、层间多次波和自由表面多次波的格林函数,利用改进的Marchenko自聚焦预测方法预测自由表面多次波.本文利用水平层状模型数据及SMARRT模型数据证明,改进后的Marchenko法预测海底相关的自由表面多次波效果较为理想,该方法避免了常规SRME自由表面多次波预测方法需要近道重构的缺陷,能够有效提高地震资料的信噪比和分辨率. 相似文献