共查询到20条相似文献，搜索用时 25 毫秒
1.
2.
共反射面元（Common Reflection Surface）叠加是一种不依赖于宏观速度模型的零炮检距剖面成像方法,实现共反射面元叠加依赖于3个波场属性参数的确定,它们分别是零偏移距射线的出射角α、Normal波和Normal Incident Point波出射到地表的波前曲率半径RN和RNIP. 在CRS叠加的理论基础上,本文阐述如何在实际数据上实现CRS叠加. 首先,通过简洁的一维相关性分析在常规叠加剖面上找到对应该共反射面元的一组初始波场属性参数(α,RN,RNIP),然后在对应的叠前数据上应用最优化算法对这组参数进行优化处理,相比初始属性参数,优化后的属性参数能够更好地聚集来自地下反射层的能量,最后应用优化后的属性参数实现最优CRS叠加. 相似文献
3.
The method of common reflection surface (CRS) extends conventional stacking of seismic traces over offset to multidimensional stacking over offset‐midpoint surfaces. We propose a new form of the stacking surface, derived from the analytical solution for reflection traveltime from a hyperbolic reflector. Both analytical comparisons and numerical tests show that the new approximation can be significantly more accurate than the conventional CRS approximation at large offsets or at large midpoint separations while using essentially the same parameters. 相似文献
4.
共反射面元（Common Reflection Surface）叠加是一种独立于宏观速度模型的零偏移距剖面成像方法,传统的CRS叠加实现是以数据驱动的方式对属性参数进行自动搜索并对其进行优化合成相应的CRS叠加算子,通过该算子进行叠加能够得到信噪比和连续性更高的零偏移距剖面.但是数据驱动的实现方式带来了不可避免的“倾角歧视现象”,它造成了弱有效反射信号损失和运动学特征失真的问题.本文提出的倾角分解CRS叠加方法成功解决了上述问题,使CRS叠加方法更具实用价值. 相似文献
5.
共反射面元（Common Reflection Surface=CRS）叠加是一种特殊的零偏移距成像方法,实践中它具有独立于宏观速度模型和完全数据驱动实现的鲜明特色,CRS叠加理论认为在得到高质量的零偏移距剖面的同时,还可以得到三个有用的波场属性参数剖面反演宏观速度模型,CRS叠加剖面之后的叠后深度偏移质量将超过叠前深度偏移.虽然CRS叠加倡导的成像方式和承诺的上述理想境界带来了全新的启示,但是实践中这些特色同样带来了令人困扰的问题,为此我们提出了倾角分解CRS叠加方法解决这些问题.本文即是作者通过上述实践之后对CRS叠加方法形成的一些思考和总结. 相似文献
6.
在地表地形复杂的情况下,静校正不易做好,这是制约山地资料处理质量的一个很重要的因素.复杂地表共反射面元（CRS）叠加不需对叠前数据做静校正,而且在得到叠加剖面后可以利用叠加得到的波场参数剖面实现基准面重建.地震数据的试算表明,复杂地表CRS叠加得出的剖面与常规处理剖面相比有着较高的信噪比和同相轴连续性.与水平地表CRS叠加不同的是,在复杂地表CRS叠加的时距公式中,波场三参数耦合,难以通过简化CRS道集的方法将它们全部分离并逐个优化.引入模拟退火算法后,有效地解决了这一组合优化的难题. 相似文献
7.
Tiago A. Coimbra Jorge H. Faccipieri Dany S. Rueda Martin Tygel 《Studia Geophysica et Geodaetica》2016,60(3):500-530
Since the early days of seismic processing, time migration has proven to be a valuable tool for a number of imaging purposes. Main motivations for its widespread use include robustness with respect to velocity errors, as well as fast turnaround and low computation costs. In areas of complex geology, in which it has well-known limitations, time migration can still be of value by providing first images and also attributes, which can be of much help in further, more comprehensive depth migration. Time migration is a very close process to common-midpoint (CMP) stacking and, more recently, to zero-offset commonreflection- surface (CRS) stacking. In fact, Kirchhoff time migration operators can be readily formulated in terms of CRS parameters. In the nineties, several studies have shown advantages in the use of common-reflection-point (CRP) traveltimes to replace conventional CMP traveltimes for a number of stacking and migration purposes. In this paper, we follow that trend and introduce a Kirchhoff-type prestack time migration and velocity analysis algorithm, referred to as CRP time migration. The algorithm is based on a CRP operator together with optimal apertures, both computed with the help of CRS parameters. A field-data example indicates the potential of the proposed technique. 相似文献
8.
大量研究证明CRS叠加能提高地震勘探的信噪比,但是能否提高分辨率的关键在于如何确定CRS叠加孔径.本文详细探讨了地震波反射过程中菲涅耳带的特征,认为起伏地形下菲涅耳带可以采用椭圆予以近似,在此基础上提出了一种通过菲涅耳带来确定CRS叠加孔径的方法,并应用于泌阳凹陷陡坡带的地震剖面.结果表明,由于菲涅耳带确定的叠加范围使地震信号的能量达到最佳,相对于CMP叠加,菲涅耳带CRS叠加同时提高了地震资料的信噪比和分辨率,特别显示了中深部较弱的地震信息,而常规的CRS叠加则只在于提高资料的信噪比和改善浅部的地震信息. 相似文献
9.
叠加速度分析技术是常规地震资料处理中的重要环节,也是经典的时间域速度建模方法.叠加速度分析技术主要包括速度谱计算和拾取两个步骤.至今为止,多数研究工作通过提高速度谱的分辨率以及抗噪声能力,获得高质量的速度谱从而有利于拾取.本文的目标是将叠加速度分析技术转为一个全自动化的处理流程.从参数估计的角度出发,将叠加速度估计转化为稀疏反演框架下的模型参数估计问题,并通过稀疏反演算法自动反演叠加速度,进而提高叠加速度建模的效率.为实现这一目标,首先给出了正问题的定义,即层状介质中CMP道集的预测模型,利用叠加速度、垂向双程走时(t_0)以及反射子波以及CMP道集时距关系(如双曲时距关系)可以预测CMP道集.接着,速度分析反问题可以描述为已知观测的CMP道集,估计模型参数(叠加速度及t_0时间等).利用模型参数的稀疏性作为约束条件并用L_0范数作为模型稀疏性的度量准则,叠加速度分析可以转化为L_0范数约束下的稀疏反演问题.本文提出了一种基于预测校正思想的匹配追踪算法求解上述反问题,实现了自动叠加速度建模并为后续的高精度速度反演方法提供较好的初始模型.理论和实际资料的测试结果证明了本文方法的有效性. 相似文献
10.
基于虚拟偏移距方法的各向异性转换波保幅叠前时间偏移 总被引：2，自引：1，他引：1
In this paper, we use the method of pseudo-offset migration （POM） to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into common conversion scatter point （CCSP） gathers directly by POM, which simplifies the conventional processing procedure for converted waves. The POM gather fold and SNR are high, which is favorable for velocity analysis and especially suitable for seismic data with low SNR. We used equivalent anisotropic theory to compute anisotropic parameters. Based on the scattering wave traveltime equation in a VTI medium, the POM pseudo-offset migration in anisotropic media was deduced. By amplitude-preserving POM gather mapping, velocity analysis, stack processing, and so on, the anisotropic migration results were acquired. The forward modeling computation and actual data processing demonstrate the validity of converted wave pre-stack time migration with amplitude-preservation using the anisotropic POM method. 相似文献
11.
Common reflection surface stack using dip decomposition for rugged surface topography 总被引：1，自引：0，他引：1
We present an extension of the Common Reflection Surface （CRS） stack that provides support for an arbitrary top surface topography. CRS stacking can be applied to the original prestack data without the need for any elevation statics. The CRS-stacked zero- offset section can be corrected （redatumed） to a given planar level by kinematic wave field attributes. The seismic processing results indicate that the CRS stacked section for rugged surface topography is better than the conventional stacked section for S/N ratio and better continuity of reflection events. Considering the multiple paths of zero-offset rays, the method deals with reflection information coming from different dips and performs the stack using the method of dip decomposition, which improves the kinematic and dynamic character of CRS stacked sections. 相似文献
12.
Static shifts from near‐surface inhomogeneities very often represent the key problem in the processing of seismic data from arid regions. In this case study, the deep bottom fill of a wadi strongly degrades the image quality of a 2D seismic data set. The resulting static and dynamic problems are solved by both conventional and common‐reflection‐surface (CRS) processing. A straightforward approach derives conventional refraction statics from picked first breaks and then goes through several iterations of manual velocity picking and residual statics calculation. The surface‐induced static and dynamic inhomogeneities, however, are not completely solved by these conventional methods. In CRS processing, the local adaptation of the CRS stacking parameters results in very detailed dynamic corrections. They resolve the local inhomogeneities that were not detected by manual picking of stacking velocities and largely compensate for the surface‐induced deterioration in the stack. The subsequent CRS residual statics calculations benefit greatly from the large CRS stacking fold which increases the numbers of estimates for single static shifts. This improves the surface‐consistent averaging of static shifts and the convergence of the static solution which removes the remaining static shifts in the 2D seismic data. The large CRS stacking fold also increases the signal‐to‐noise ratio in the final CRS stack. 相似文献
13.
In conventional seismic processing, the classical algorithm of Hubral and Krey is routinely applied to extract an initial macrovelocity model that consists of a stack of homogeneous layers bounded by curved interfaces. Input for the algorithm are identified primary reflections together with normal moveout (NMO) velocities, as derived from a previous velocity analysis conducted on common midpoint (CMP) data. This work presents a modified version of the Hubral and Krey algorithm that is designed to extend the original version in two ways, namely (a) it makes an advantageous use of previously obtained common-reflection-surface (CRS) attributes as its input and (b) it also allows for gradient layer velocities in depth. A new strategy to recover interfaces as optimized cubic splines is also proposed. Some synthetic examples are provided to illustrate and explain the implementation of the method. 相似文献
14.
We review the multifocusing method for traveltime moveout approximation of multicoverage seismic data. Multifocusing constructs the moveout based on two notional spherical waves at each source and receiver point, respectively. These two waves are mutually related by a focusing quantity. We clarify the role of this focusing quantity and emphasize that it is a function of the source and receiver location, rather than a fixed parameter for a given multicoverage gather. The focusing function can be designed to make the traveltime moveout exact in certain generic cases that have practical importance in seismic processing and interpretation. The case of a plane dipping reflector (planar multifocusing) has been the subject of all publications so far. We show that the focusing function can be generalized to other surfaces, most importantly to the spherical reflector (spherical multifocusing). At the same time, the generalization implies a simplification of the multifocusing method. The exact traveltime moveout on spherical surfaces is a very versatile and robust formula, which is valid for a wide range of offsets and locations of source and receiver, even on rugged topography. In two‐dimensional surveys, it depends on the same three parameters that are commonly used in planar multifocusing and the common‐reflection surface (CRS) stack method: the radii of curvature of the normal and normal‐incidence‐point waves and the emergence angle. In three dimensions the exact traveltime moveout on spherical surfaces depends on only one additional parameter, the inclination of the plane containing the source, receiver and reflection point. Comparison of the planar and spherical multifocusing with the CRS moveout expression for a range of reflectors with increasing curvature shows that the planar multifocusing can be remarkably accurate but the CRS becomes increasingly inaccurate. This can be attributed to the fact that the CRS formula is based on a Taylor expansion, whereas the multifocusing formulae are double‐square root formulae. As a result, planar and spherical multifocusing are better suited to model the moveout of diffracted waves. 相似文献
15.
是否能够正确地建立深度域三维速度模型是三维叠前深度偏移成败的关键 .本文根据Deregowski循环 ,利用叠前深度域地震成像对速度模型变化的敏感性 ,采用偏移迭代逐次逼近最佳成像速度 ,研究开发了一套快捷有效的三维叠前深度偏移深度域速度模型建立技术 .借鉴时间域CDP(共深度点 )道集上常规叠加速度分析的策略 ,在深度域CRP(共反射点 )道集上 ,提出剩余慢度平方谱的概念并建立相应的实现技术 .导出深度域中均方根速度与层速度之间的关系 ;按照串级偏移原理确定偏移循环过程中初始速度、剩余速度及修改后速度之间的关系 ;采用蒙特卡洛非线性优化算法实现从剩余慢度平方谱中自动拾取层速度 ,讨论了其地质速度约束条件和蒙特卡洛非线性优化的收敛准则 ,使得所拾取的层速度模型具有合理的地质意义并获得最佳偏移成像效果 .SEG EAGE理论模型数值试算验证了方法的有效性 ,在海拉尔盆地霍多莫尔工区 ,5 8km2 三维资料的速度模型建立并获得满意的三维叠前深度偏移成像 . 相似文献
16.
立体层析成像是一种新的地震反射波层析成像方法,能为叠前深度偏移提供较为精确的宏观速度模型。本文研究了立体层析成像的实现方法,包括斜率与走时数据的拾取、离散速度模型构建和初始化、射线参数的确定、斜率和走时及射线计算以及反演问题解法等,建立立体层析成像的算法流程。并通过对Marmousi模型试验,对立体层析成像运行所需的主要参数,如初始速度模型、拾取数据量、离散网格尺寸、速度平滑权重等进行测试和分析,总结这些不同参数对立体层析反演结果的影响规律,用以指导生产实践。 相似文献
17.
Velocity model updating in prestack Kirchhoff time migration for PS converted waves: Part I – Theory
A velocity model updating approach is developed based on moveout analysis of the diffraction curve of PS converted waves in prestack Kirchhoff time migration. The diffraction curve can be expressed as a product of two factors: one factor depending on the PS converted‐wave velocity only, and the other factor depending on all parameters. The velocity‐dependent factor represents the hyperbolic behaviour of the moveout and the other is a scale factor that represents the non‐hyperbolic behaviour of the moveout. This non‐hyperbolic behaviour of the moveout can be corrected in prestack Kirchhoff time migration to form an inverse normal‐moveout common‐image‐point gather in which only the hyperbolic moveout is retained. This hyperbolic moveout is the moveout that would be obtained in an isotropic equivalent medium. A hyperbolic velocity is then estimated from this gather by applying hyperbolic moveout analysis. Theoretical analysis shows that for any given initial velocity, the estimated hyperbolic velocity converges by an iterative procedure to the optimal velocity if the velocity ratio is optimal or to a value closer to the optimal velocity if the velocity ratio is not optimal. The velocity ratio (VP/VS) has little effect on the estimation of the velocity. Applying this technique to a synthetic seismic data set confirms the theoretical findings. This work provides a practical method to obtain the velocity model for prestack Kirchhoff time migration. 相似文献
18.
19.
R.-G. Ferber 《Geophysical Prospecting》1994,42(2):99-112
The stacking velocity best characterizes the normal moveout curves in a common-mid-point gather, while the migration velocity characterizes the diffraction curves in a zero-offset section as well as in a common-midpoint gather. For horizontally layered media, the two velocity types coincide due to the conformance of the normal and the image ray. In the case of dipping subsurface structures, stacking velocities depend on the dip of the reflector and relate to normal rays, but with a dip-dependent lateral smear of the reflection point. After dip-moveout correction, the stacking velocities are reduced while the reflection-point smear vanishes, focusing the rays on the common reflection points. For homogeneous media the dip-moveout correction is independent of the actual velocity and can be applied as a dip-moveout correction to multiple offset before velocity analysis. Migration to multiple offset is a prestack, time-migration technique, which presents data sets which mimic high-fold, bin-centre adjusted, common-midpoint gathers. This method is independent of velocity and can migrate any 2D or 3D data set with arbitrary acquisition geometry. The gathers generated can be analysed for normal-moveout velocities using traditional methods such as the interpretation of multivelocity-function stacks. These stacks, however, are equivalent to multi-velocity-function time migrations and the derived velocities are migration velocities. 相似文献
20.
Local seismic event slopes contain subsurface velocity information and can be used to estimate seismic stacking velocity. In this paper, we propose a novel approach to estimate the stacking velocity automatically from seismic reflection data using similarity‐weighted k‐means clustering, in which the weights are local similarity between each trace in common midpoint gather and a reference trace. Local similarity reflects the local signal‐to‐noise ratio in common midpoint gather. We select the data points with high signal‐to‐noise ratio to be used in the velocity estimation with large weights in mapped traveltime and velocity domain by similarity‐weighted k‐means clustering with thresholding. By using weighted k‐means clustering, we make clustering centroids closer to those data points with large weights, which are more reliable and have higher signal‐to‐noise ratio. The interpolation is used to obtain the whole velocity volume after we have got velocity points calculated by weighted k‐means clustering. Using the proposed method, one obtains a more accurate estimate of the stacking velocity because the similarity‐based weighting in clustering takes into account the signal‐to‐noise ratio and reliability of different data points in mapped traveltime and velocity domain. In order to demonstrate that, we apply the proposed method to synthetic and field data examples, and the resulting images are of higher quality when compared with the ones obtained using existing methods. 相似文献