共查询到20条相似文献，搜索用时 375 毫秒
1.
地震资料含有各种类型多次波，而传统成像方法仅利用地震一次反射波成像，在地震成像前需将多次波去除.然而，多次波携带了丰富的地下结构信息，多次波偏移能够提供除反射波外的额外地下照明.修改传统逆时偏移方法，用包含一次反射波和多次波的原始记录代替震源子波，将SRME方法预测的表面多次波代替一次反射波作为输入数据，可将表面多次波成像.多次波成像的挑战和困难在于大量串扰噪声的产生，针对表面多次波成像中的成像噪声问题，将最小二乘逆时偏移方法与多次波分阶思想结合起来，发展可控阶数的表面多次波反演成像方法，有望初步实现高精度的表面多次波成像.在消除原始记录中的表面多次波后，通过逆散射级数方法预测得到层间多次波，将层间多次波作为逆时偏移方法的输入数据可将其准确归位到地下反射位置.数值实验表明，多次波成像能够有效地为地下提供额外照明，而可控阶表面多次波最小二乘逆时偏移成像方法几乎完全避免成像噪声. 相似文献
2.
3.
在常规地震资料处理中,多次反射波被视为噪声并从地震数据中去除,以免在之后的地震资料解释中造成误解.而事实上,多次波也是地震信号,是照明波场的一部分,能够对地下构造成像的精度做出贡献.本文分析了多次波在传统单程波叠前深度偏移中产生构造假象的机制和表现,为实现基于单程波偏移算子的多次波成像,修改了单程波叠前深度偏移的边界条件,即将输入的震源波场用包含多次波的记录来替代,输入的记录波场用预测出的表层相关多次波来替代,实现了基于单程波偏移算子的地表相关多次波成像,并从理论上给出了其成像依据.通过基于二范式最小能量差原则求取的匹配因子,将多次波成像结果与一次波成像结果进行匹配叠加,应用多次波成像来弥补一次波成像的不足.简单模型验证了基于单程波偏移算子的多次波成像方法的有效性,最后对Sigsbee2B模型进行了一次波与多次波联合成像试算,盐边界高陡构造成像质量得到了明显改善. 相似文献
4.
在常规偏移方法中一般都需要压制地震数据中的多次波,仅利用一次波信息成像,把自由表面反射的多次波视为噪声,但是在多次波中也包含着地下结构信息,应该将其充分利用到成像中来.事实上,已经有不少成像方法试图利用多次波信息,但是大部分方法都需要对多次波进行预测.本文提出了基于傅里叶有限差分偏移算子的数据自相关偏移方法.在这种偏移方法中,对含有一次波和多次波的地震数据,分别进行下行和上行延拓,然后直接利用常规的互相关成像条件成像.由于波场延拓采用了傅里叶有限差分算子,其计算效率高,能够很好地对复杂介质中的地震数据进行延拓.在数值试验中,使用了一个含散射点的三层模型和Marmousi模型.合成数据测试结果表明,这种方法可以对更大范围的地下构造成像,比常规的只利用一次波的傅里叶有限差分法照明度更好,并且在浅层可以提供更高的分辨率.我们提出的数据自相关策略易于实现且避免了繁杂的多次波预测,这对于复杂地下构造成像可能有着重大意义. 相似文献
5.
A.J. Berkhout 《Geophysical Prospecting》2017,65(1):106-145
Surface removal and internal multiple removal are explained by recursively separating the primary and multiple responses at each depth level with the aid of wavefield prediction error filtering. This causal removal process is referred to as “data linearization.” The linearized output (primaries only) is suitable for linear migration algorithms. Next, a summary is given on the migration of full wavefields (primaries + multiples) by using the concept of secondary sources in each subsurface gridpoint. These secondary sources are two‐way and contain the gridpoint reflection and the gridpoint transmission properties. In full wavefield migration, a local inversion process replaces the traditional linear imaging conditions. Finally, Marchenko redatuming is explained by iteratively separating the full wavefield response from above a new datum and the full wavefield response from below a new datum. The redatuming output is available for linear migration (Marchenko imaging) or, even better, for full wavefield migration. Linear migration, full wavefield migration, and Marchenko imaging are compared with each other. The principal conclusion of this essay is that multiples should not be removed, but they should be utilized, yielding two major advantages: (i) illumination is enhanced, particularly in the situation of low signal‐to‐noise primaries; and (ii) both the upper side and the lower side of reflectors are imaged. It is also concluded that multiple scattering algorithms are more transparent if they are formulated in a recursive depth manner. In addition to transparency, a recursive depth algorithm has the flexibility to enrich the imaging process by inserting prior geological knowledge or by removing numerical artefacts at each depth level. Finally, it is concluded that nonlinear migration algorithms must have a closed‐loop architecture to allow successful imaging of incomplete seismic data volumes (reality of field data). 相似文献
6.
The key processes in marine seismic imaging include (i) removing from seismic data all seismic events (free-surface multiples and ghosts) which contain at least one reflection at the sea surface in their wave-propagation path, and leaving those with no reflection at the free surface (internal multiples and primaries), (ii) removing events with at least two reflections in the subsurface (internal multiples), and leaving events with only one reflection in the subsurface (primaries), and then (iii) locating the scattering points and reflectors inside the subsurface which are the sources of primaries and internal multiple events. All these processes are here explained, derived, and optimized via scattering diagrams (diagrammatica) in a way similar to the way the quantum field theory is often explained via Feynman diagrams. Our discussion of the removal of events with free-surface reflections from the data will be brief, as the diagrammatica of these events are now well understood.The main focus of this paper is the diagrammatica of internal multiples and primaries. Although these events do not contain any reflection at the sea surface, it is important to reconstruct them with scattering points near the sea surface, where seismic data are recorded. So our diagrammatica of primaries and internal multiples include events which are not directly recorded in seismic data but which can be constructed from seismic data. These events have allowed us to construct scattering diagrams of primaries and internal multiples with scattering points near the sea surface. Furthermore, these new diagrammatica of internal multiples and primaries can be used to remove internal multiples from the data. 相似文献
7.
本文深入分析了多次波逆时偏移的成像原理和串扰假象产生机制,并提出了多次波分阶逆时偏移策略,即首先应用自由界面多次波衰减方法对原始炮集记录中的多次波进行剔除得到一次波记录,然后应用一次波记录预测获得各阶多次波记录,最后将各阶多次波记录分别进行逆时偏移成像.模型实验结果表明,多次波分阶逆时偏移其各阶多次波成像剖面均可对真实界面正确成像且其能够有效压制串扰假象,其中一阶多次波逆时偏移剖面的成像精度最高,其深层成像质量明显优于常规多次波逆时偏移. 相似文献
8.
Reverse time migration of multiples can be used to construct subsurface structures where primaries cannot illuminate well. However, the images generated using multiples suffer from severe artefacts due to the cross-talks created by interference among unrelated multiples. We developed a migration approach using water-bottom-related multiples to reduce these cross-talk artefacts. This approach first isolates primaries from the original data and predicts water-column primaries. The nth-order water-column multiples can be obtained by auto-convolving the water-column primaries n times, followed by convolving the nth-order water-column multiples with the primaries to extract the (n+1)th-order water-bottom-related multiples. The approach takes the nth-order water-column multiples as the secondary source and regards the (n+1)th-order water-bottom-related multiples as the receiver wavefield, followed by a cross-correlation imaging condition. Numerical examples from synthetic and field data sets reveal that our approach can provide images with substantially fewer cross-talk artefacts than conventional reverse time migration using multiples, as well as greatly improving shallow imaging compared with reverse time migration of primaries. 相似文献
9.
Migration methods for imaging different-order multiples 总被引：2，自引：0，他引：2
Zhiyong Jiang Jianming Sheng Jianhua Yu Gerard T. Schuster Brian E. Hornby 《Geophysical Prospecting》2007,55(1):1-19
Multiples contain valuable information about the subsurface, and if properly migrated can provide a wider illumination of the subsurface compared to imaging with VSP primary reflections. In this paper we review three different methods for migrating multiples. The first method is model-based, and it is more sensitive to velocity errors than primary migration; the second method uses a semi-natural Green's function for migrating multiples, where part of the traveltimes are computed from the velocity model, and part of the traveltimes (i.e., natural traveltimes) are picked from the data to construct the imaging condition for multiples; the third method uses cross-correlation of traces. The last two methods are preferred in the sense that they are significantly less sensitive to velocity errors and statics because they use “natural data” to construct part of the migration imaging conditions. Compared with the interferometric (i.e., crosscorrelation) imaging method the semi-natural Green's function method is more computationally efficient and is sometimes less prone to migration artifacts. Numerical tests with 2-D and 3-D VSP data show that a wider subsurface coverage, higher-fold and more balanced illumination of the subsurface can be achieved with multiple migration compared with migration of primary reflections only. However, there can be strong interference from multiples with different orders or primaries when multiples of high order are migrated. One possible solution is to filter primaries and different orders of multiples before migration, and another possible solution is least squares migration of all events. A limitation of multiple migration is encountered for subsalt imaging. Here, the multiples must pass through the salt body more than twice, which amplifies the distortion of the image. 相似文献
10.
使用相同的炮记录,多次波偏移能提供比反射波偏移更广的地下照明和更多的地下覆盖但是同时产生很多的串声噪声.相比传统逆时偏移,最小二乘逆时偏移反演的反射波成像结果具有更高的分辨率和更均衡的振幅.我们主要利用最小二乘逆时偏移压制多次波偏移产生的串声噪声.多次波最小二乘逆时偏移通常需要一定的迭代次数以较好地消除串声噪声.若提前将一阶多次波从所有阶数的多次波中过滤出来,使用相同的迭代次数,一阶多次波的最小二乘逆时偏移能够得到具有更高信噪比的成像剖面,而且能够提供与多次波最小二乘逆时偏移相似的有效地下结构成像. 相似文献
11.
Multiple scattering is usually ignored in migration algorithms, although it is a genuine part of the physical reflection response. When properly included, multiples can add to the illumination of the subsurface, although their crosstalk effects are removed. Therefore, we introduce full‐wavefield migration. It includes all multiples and transmission effects in deriving an image via an inversion approach. Since it tries to minimize the misfit between modeled and observed data, it may be considered a full waveform inversion process. However, full‐wavefield migration involves a forward modelling process that uses the estimated seismic image (i.e., the reflectivities) to generate the modelled full wavefield response, whereas a smooth migration velocity model can be used to describe the propagation effects. This separation of modelling in terms of scattering and propagation is not easily achievable when finite‐difference or finite‐element modelling is used. By this separation, a more linear inversion problem is obtained. Moreover, during the forward modelling, the wavefields are computed separately in the incident and scattered directions, which allows the implementation of various imaging conditions, such as imaging reflectors from below, and avoids low‐frequency image artefacts, such as typically observed during reverse‐time migration. The full wavefield modelling process also has the flexibility to image directly the total data (i.e., primaries and multiples together) or the primaries and the multiples separately. Based on various numerical data examples for the 2D and 3D cases, the advantages of this methodology are demonstrated. 相似文献
12.
海水与空气间的强波阻抗差使得海洋地震资料普遍发育自由表面相关多次波,如何利用好多次波所携带的有效信息已成为提高海洋地震资料成像品质的新突破点.基于面炮偏移的一次波与多次波同时成像方法能够避免多次波预测精度的影响,但是,正向传播的震源子波与反向延拓的自由表面相关多次波所产生的干涉假象严重制约了该技术的应用,本文提出了一种基于单程波偏移算子,可在成像域压制干涉假象的一次波与多次波同时成像方法.其中包含了三个步骤:第一,传统单程波偏移成像方法中的震源子波替换为一次波、多次波与震源子波,初始上行延拓波场为一次波与多次波,基于单程波算子的波场延拓与互相关成像条件的应用得到包含干涉假象的一次波与多次波同时成像;第二,以子波为震源,自由表面相关多次波为记录,按照传统单程波偏移成像方法得到干涉假象;第三,基于最小二乘匹配滤波算法,将第一步的成像结果与第二步的干涉假象进行匹配相减,得到干涉假象衰减后的一次波与多次波同时成像,避开了由于实际资料子波无法准确提取而造成一次波与多次波对成像能量级的不一致性.Sigsbee2B模型测试验证了本方法的有效性,并在我国某探区深海实际资料处理中得到了成功应用,深层基底得到了清晰刻画,并且照明均衡度明显改善. 相似文献
13.
Multiples have longer propagation paths and smaller reflection angles than primaries for the same source–receiver combination, so they cover a larger illumination area. Therefore, multiples can be used to image shadow zones of primaries. Least-squares reverse-time migration of multiples can produce high-quality images with fewer artefacts, high resolution and balanced amplitudes. However, viscoelasticity exists widely in the earth, especially in the deep-sea environment, and the influence of Q attenuation on multiples is much more serious than primaries due to multiples have longer paths. To compensate for Q attenuation of multiples, Q-compensated least-squares reverse-time migration of different-order multiples is proposed by deriving viscoacoustic Born modelling operators, adjoint operators and demigration operators for different-order multiples. Based on inversion theory, this method compensates for Q attenuation along all the propagation paths of multiples. Examples of a simple four-layer model, a modified attenuating Sigsbee2B model and a field data set suggest that the proposed method can produce better imaging results than Q-compensated least-squares reverse-time migration of primaries and regular least-squares reverse-time migration of multiples. 相似文献
14.
将多次波转换成反射波并按传统反射波偏移算法成像,是多次波成像的一种方法.聚焦变换能准确的将多次波转换为纵向分辨率更高的新波场记录,其中一阶多次波转换为反射波.本文对聚焦变换提出了两点改进:1)提出局部聚焦变换,以减小存储量和计算量,增强该方法对检波点随炮点移动的采集数据的适应性;2)引入加权矩阵,理论上证明原始记录的炮点比检波点稀疏时,共检波点道集域的局部聚焦变换可以将多次波准确转换成炮点与检波点有相同采样频率的新波场记录.本文在第一个数值实验中对比了对包含反射波与多次波的原始记录做局部聚焦变换和直接对预测的多次波做局部聚焦变换两种方案,验证了第二种方案转换得到的波场记录信噪比更高且避免了第一个方案中切聚焦点这项比较繁杂的工作.第二个数值实验表明:在炮点采样较为稀疏时,该方法能有效的将一阶多次波转换成反射波;转换的反射波能提供更丰富的波场信息,成像结果更均衡、在局部有更高的信噪比,以及较高的纵向分辨率. 相似文献
15.
Giovanni Angelo Meles Lele Zhang Jan Thorbecke Kees Wapenaar Evert Slob 《Geophysical Prospecting》2020,68(6):1834-1846
Seismic images provided by reverse time migration can be contaminated by artefacts associated with the migration of multiples. Multiples can corrupt seismic images, producing both false positives, that is by focusing energy at unphysical interfaces, and false negatives, that is by destructively interfering with primaries. Multiple prediction/primary synthesis methods are usually designed to operate on point source gathers and can therefore be computationally demanding when large problems are considered. A computationally attractive scheme that operates on plane-wave datasets is derived by adapting a data-driven point source gathers method, based on convolutions and cross-correlations of the reflection response with itself, to include plane-wave concepts. As a result, the presented algorithm allows fully data-driven synthesis of primary reflections associated with plane-wave source responses. Once primary plane-wave responses are estimated, they are used for multiple-free imaging via plane-wave reverse time migration. Numerical tests of increasing complexity demonstrate the potential of the proposed algorithm to produce multiple-free images from only a small number of plane-wave datasets. 相似文献
16.
海水与空气间的强波阻抗界面使得海洋地震数据普遍发育自由表面相关多次波,多次波信息的利用是提高海洋地震资料成像品质的新突破点.近年来发展了一系列多次波成像方法,干涉假象是制约其应用推广的关键问题之一.为了避免假象影响,本文提出了不同阶次自由表面相关多次波预测与成像方法,首先,修改了传统SRME(表面相关多次波衰减)方法中的边界条件,通过多次波升阶次与匹配相减的方法预测出不同阶次自由表面相关多次波;其次,基于单程波偏移算子和"面炮"偏移策略,以一次反射波或第(N-1)阶自由表面相关多次波为下行波场正向延拓,以第1阶多次或第N阶多次波为上行波场逆向延拓,并在每一层互相关成像得到第1阶或N阶多次波单独成像.本方法避免了低阶多次波和高阶多次波产生的相关假象,且相对于全波算子的偏移方法具有较高的计算效率,增强了多次波成像方法的实用性.单层模型和三层模型测试验证了本方法的正确性,并在我国某深海探区实际资料处理中得到了成功应用.相对于传统一次波成像,分阶次多次波成像具有更高的照明均衡度、垂向分辨率和信噪比.本研究表明,海洋多次波成像是一次波成像的有力补充,对于稳定海底沉积的深海地区,具有一定的应用前景. 相似文献
17.
Extending illumination using all multiples: application to 3D acquisition geometry analysis 
下载免费PDF全文
下载免费PDF全文 Recent advances in survey design have led to conventional common‐midpoint‐based analysis being replaced by subsurface‐based seismic acquisition analysis, with emphasis on advanced techniques of illumination analysis. Among them is the so‐called focal beam method, which is a wave‐equation‐based seismic illumination analysis method. The objective of the focal beam method is to provide a quantitative insight into the combined influence of acquisition geometry, overburden structure, and migration operators on the resolution and angle‐dependent amplitude fidelity of the image. The method distinguishes between illumination and sensing capability of a particular acquisition geometry by computing the focal source beam and the focal detector beam, respectively. Sensing is related to the detection properties of a detector configuration, whereas illumination is related to the emission properties of a source configuration. The focal source beam analyses the incident wavefield at a specific subsurface grid point from all available sources, whereas the focal detector beam analyses the sensing wavefield reaching at the detector locations from the same subsurface grid point. In the past, this method could only address illumination by primary reflections. In this paper, we will extend the concept of the focal beam method to incorporate the illumination due to the surface and internal multiples. This in fact complies with the trend of including multiples in the imaging process. Multiple reflections can illuminate a target location from other angles compared with primary reflections, resulting in a higher resolution and an improved illumination. We demonstrate how an acquisition‐related footprint can be corrected using both the surface and the internal multiples. 相似文献
18.
VSP资料上下行波场发育丰富.本文在分析VSP直达波、上行反射波、下行反射波传播路径及其照明范围的基础上,指出了常规VSP波动方程偏移方法缺陷,进而通过修改波场延拓方式,提出了上下行反射波联合成像方法,并在高频近似下分析了该方法的成像原理.该方法不需要进行VSP上下行反射波场分离,能够同时对VSP资料中的一次反射波、自由表面多次波、层间多次波进行成像,比常规成像剖面具有更宽的成像范围和更好的成像效果.该方法能够对下行一次反射波进行成像,从而可以实现常规偏移方法难以处理的高陡倾角构造成像.模拟资料和实际资料处理证明了本文方法的正确性. 相似文献
19.
Antoine Guitton 《Geophysical Prospecting》2006,54(2):135-152
Surface‐related multiples are attenuated for one sail line and one streamer of a 3D data set (courtesy of Compagnie Générale de Géophysique). The survey was carried out in the Gulf of Mexico in the Green Canyon area where salt intrusions close to the water‐bottom are present. Because of the complexity of the subsurface, a wavefield method incorporating the full 3D volume of the data for multiple removal is necessary. This method comprises modelling of the multiples, where the data are used as a prediction operator, and a subtraction step, where the model of the multiples is adaptively removed from the data with matching filters. The accuracy of the multiple model depends on the source/receiver coverage at the surface. When this coverage is not dense enough, the multiple model contains errors that make successful subtraction more difficult. In these circumstances, one can either (1) improve the modelling step by interpolating the missing traces, (2) improve the subtraction step by designing methods that are less sensitive to modelling errors, or (3) both. For this data set, the second option is investigated by predicting the multiples in a 2D sense (as opposed to 3D) and performing the subtraction with a pattern‐based approach. Because some traces and shots are missing for the 2D prediction, the data are interpolated in the in‐line direction using a hyperbolic Radon transform with and without sparseness constraints. The interpolation with a sparseness constraint yields the best multiple model. For the subtraction, the pattern‐based technique is compared with a more standard, adaptive‐subtraction scheme. The pattern‐based approach is based on the estimation of 3D prediction‐error filters for the primaries and the multiples, followed by a least‐squares estimation of the primaries. Both methods are compared before and after prestack depth migration. These results suggest that, when the multiple model is not accurate, the pattern‐based method is more effective than adaptive subtraction at removing surface‐related multiples while preserving the primaries. 相似文献