首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Most amplitude versus offset (AVO) analysis and inversion techniques are based on the Zoeppritz equations for plane‐wave reflection coefficients or their approximations. Real seismic surveys use localized sources that produce spherical waves, rather than plane waves. In the far‐field, the AVO response for a spherical wave reflected from a plane interface can be well approximated by a plane‐wave response. However this approximation breaks down in the vicinity of the critical angle. Conventional AVO analysis ignores this problem and always utilizes the plane‐wave response. This approach is sufficiently accurate as long as the angles of incidence are much smaller than the critical angle. Such moderate angles are more than sufficient for the standard estimation of the AVO intercept and gradient. However, when independent estimation of the formation density is required, it may be important to use large incidence angles close to the critical angle, where spherical wave effects become important. For the amplitude of a spherical wave reflected from a plane fluid‐fluid interface, an analytical approximation is known, which provides a correction to the plane‐wave reflection coefficients for all angles. For the amplitude of a spherical wave reflected from a solid/solid interface, we propose a formula that combines this analytical approximation with the linearized plane‐wave AVO equation. The proposed approximation shows reasonable agreement with numerical simulations for a range of frequencies. Using this solution, we constructed a two‐layer three‐parameter least‐squares inversion algorithm. Application of this algorithm to synthetic data for a single plane interface shows an improvement compared to the use of plane‐wave reflection coefficients.  相似文献   

2.
The growing awareness of the significance of three-dimensional interpretation of seismic reflection data carries with it a need for better understanding of the role of the velocity configuration in the production of reflection patterns. Modeling the reflection responses of a simple dipping plane reflector through velocity models of various degrees of complexity demonstrates the importance of overlying velocities in determining the reflection pattern. Modeling is accomplished using a raytracing technique which determines total travel time of the normal incidence raypath through an arbitrary iso-velocity layer model.  相似文献   

3.
转换波和非转换波与射线参数p的关系   总被引:1,自引:1,他引:1  
本文对Aki固体-固体分界面的反射和透射公式进行适当变换,使之成为一组适合不同坐标系下,固体、液体、气体和真空之间所有可能形成的分界面上,入射P、SV和SH波所形成的反射和透射公式,从而把不同边界条件的反射和透射方程统一起来,并且阐明了转换波和非转换波动力学特征的差异.指出:在不同弹性介质分界面上入射平面弹性波,转换波的反射系数和透射系数是射线参数p的奇函数;非转换波的反射系数和透射系数是射线参数p的偶函数.根据奇函数和偶函数的性质可见:垂直入射时不存在转换波.同时也为简化Zoeppricz方程,开展AVO分析打开新思路.  相似文献   

4.
直接找油气的地球物理新原理的进展现状   总被引:2,自引:0,他引:2       下载免费PDF全文
近年来,国际上利用地震反射波勘探技术直接找油气的探索研究受到普遍重视.本文试图针对这一发展历史作一综述. 首先,介绍Knott(1899)和Zoeppritz(1919)分别导出的反射波振幅公式进行的数值计算,从而发现振幅与入射角的变化关系及其受反射面两侧岩层泊松比值的影响.然后,列举了几处油气区地震波勘探直接找油气的实例,以试验此项新原理在实践中所遇到的问题,说明直接找油气的地球物理新方法是可取的,是有前途的,但目前还不够有效,要达到实用还有一段较长的距离,因此,有进一步加强基础性实验研究的必要.  相似文献   

5.
地震属性分析技术在地球物理勘探领域的广泛应用,启发研究人员将其应用于人工源宽角反射/折射深地震测深剖面的资料预处理和震相识别。采用札达-泉水沟深地震测深资料,提取振幅、信噪比、主频、瞬时带宽、瞬时高频能量等地震属性参数,分析不同参数的物理含义,挑选其中对界面变化敏感的参数,对深地震测深资料进行预处理,并利用P波和S波的联合扫描,提高震相识别的准确性。走时互换结果显示,采用地震属性参数可有效提高震相拾取的准确性,进而提高后续地壳速度结构反演结果的精度。  相似文献   

6.
Seismic amplitude variations with offset contain information about the elastic parameters. Prestack amplitude analysis seeks to extract this information by using the variations of the reflection coefficients as functions of angle of incidence. Normally, an approximate formula is used for the reflection coefficients, and variations with offset of the geometrical spreading and the anelastic attenuation are often ignored. Using angle of incidence as the dependent variable is also computationally inefficient since the data are recorded as a function of offset. Improved approximations have been derived for the elastic reflection and transmission coefficients, the geometrical spreading and the complex travel-time (including anelastic attenuation). For a 1 D medium, these approximations are combined to produce seismic reflection amplitudes (P-wave, S-wave or converted wave) as a Taylor series in the offset coordinate. The coefficients of the Taylor series are computed directly from the parameters of the medium, without using the ray parameter. For primary reflected P-waves, dynamic ray tracing has been used to compute the offset variations of the transmission coefficients, the reflection coefficient, the geometrical spreading and the anelastic attenuation. The offset variation of the transmission factor is small, while the variations in the geometrical spreading, absorption and reflection coefficient are all significant. The new approximations have been used for seismic modelling without ray tracing. The amplitude was approximated by a fourth-order polynomial in offset, the traveltime by the normal square-root approximation and the absorption factor by a similar expression. This approximate modelling was compared to dynamic ray tracing, and the results are the same for zero offset and very close for offsets less than the reflector depth.  相似文献   

7.
Utilising ambient seismic energy naturally propagating in the Earth as an alternative approach to active body-wave seismic investigations has been a topic of interest for a number of decades. However, because ambient surface-wave arrivals typically are of much greater amplitude than ambient body-wave energy, significant data signal processing and long recording times are required to mitigate this and other coherent noise sources, and to correlate sufficient reflected body-wave energy to converge to a stable image. Even for these scenarios, identifying and validating imaged body-wave reflection events remain challenging. In active-source investigations, extended imaging condition gathers are used to examine velocity (in)accuracy. Herein, we develop an ambient direct migration approach that uses a novel ambient (deconvolution) extended imaging condition. We simulate synthetic ambient-wavefield seismic data for two different models and use a field data set from Lalor Lake in Manitoba, Canada, to conduct a series of numerical experiments to demonstrate the velocity sensitivity and long-term stationarity of ambient-wavefield seismic data in the migration image domain. Tests with varying global velocity perturbations show a characteristic reflector moveout in deconvolution extended imaging condition gathers that can serve as a diagnostic of reflected ambient body-wave energy. We illustrate that this imaging formalism, under idealised circumstances, gives comparable results to conventional seismic methods, which extends the use of extended imaging condition gather-based image validation to ambient-wavefield seismic data scenarios. We assert that this may be a valuable tool for the validation of ambient migration techniques that to date have yielded largely inconclusive results.  相似文献   

8.
洱源-江川宽角地震剖面的地壳反射特征   总被引:6,自引:0,他引:6       下载免费PDF全文
通过对洱源-江川宽角反射地震资料的叠前偏移处理成像,得到了一个类似于近垂直入射多道地震数据的记录剖面.反射剖面图像与地壳速度结构共同揭示出地壳厚度由剖面西北端(约45 km)向东南端(约40 km)减薄.在洱源-楚雄西北附近深度约10 km处存在一组向东南倾斜的强反射震相,其东南约50 km处存在向西北倾斜的强射震相.易门断裂两侧地壳反射属性具有明显的差异,易门断裂之西北方向深度25~40 km处,中下地壳内存在两组强振幅,向东南方向上倾的地震反射同相轴,并被楚雄-建水断裂后期所错断,易门断裂之东南方向上,地壳内反射较为均一,未见特别明显的强间断面反射信息,这个反射结构被解释为印度板块东向俯冲与藏东缘地壳物质东向逃逸综合作用导致下地壳增厚和厚地壳变形的结果.  相似文献   

9.
When a seismic source is placed in the water at a height less than a wavelength from the water–solid interface, a prominent S‐wave arrival can be observed. It travels kinematically as if it was excited at the projection point of the source on the interface. This non‐geometric S‐wave has been investigated before, mainly for a free‐surface configuration. However, as was shown in a field experiment, the non‐geometric S‐wave can also be excited at a fluid–solid configuration if the S‐wave speed in the solid is less than the sound speed in the water. The amplitude of this wave exponentially decreases when the source is moved away from the interface revealing its evanescent character in the fluid. In the solid, this particular converted mode is propagating as an ordinary S‐wave and can be transmitted and reflected as such. There is a specific region of horizontal slownesses where this non‐geometric wave exists, depending on the ratio of the S‐wave velocity and the sound speed of water. Only for ratios smaller than 1, this wave appears. Lower ratios result in a wider region of appearance. Due to this property, this particular P‐S converted mode can be identified and filtered from other events in the Radon domain.  相似文献   

10.
基于照明补偿的单程波最小二乘偏移   总被引:7,自引:6,他引:1       下载免费PDF全文
最小二乘偏移是一种基于反射地震数据与地下反射率间线性关系而建立起来的地震数据线性反演方法,相比常规偏移成像具有更好的保幅性能.本文提出了一种基于照明补偿的单程波最小二乘偏移方法,首先利用单程波方程的稳定Born近似广义屏波场传播算子构建反射地震数据与地下反射率间的线性算子,然后再应用线性最优化方法求解最小二乘偏移所对应的线性反问题.在迭代求解最优化问题的过程中,以地震波场的地下照明强度作为迭代反演的预条件算子加快迭代的收敛速度.单程波传播过程中考虑了速度分界面产生的透射效应,并用单极震源代替常规偏移中的偶极震源.把本文提出的方法应用于层状理论模型和Marmosi模型地震数据的数值试验中均取得了理想的结果.  相似文献   

11.
汤阴地堑位于太行山脉与华北平原的过渡带,是太行山前重要的地质构造单元。为研究汤阴地堑上地壳结构和断裂特征,利用安阳市与新乡市活断层探测获得的深、浅地震反射剖面,结合研究区已有地震、地质资料,对汤阴地堑浅部结构、汤东断裂特征进行分析研究。结果表明,汤阴地堑北部与南部地壳结构差异明显。地堑北部是由汤东断裂控制的半地堑,地堑内反射震相丰富、反射波层组关系清晰,多组新生代、古生代沉积层强反射不整合地覆盖在自西向东倾伏的结晶基底反射Tg上,且随着深度增加,地层倾角增大,显示出明显的多期掀斜运动特征。地堑南部表现为由汤东、汤西断裂共同控制的断陷型地堑,地堑内新近系底界面反射波TN自西向东倾伏,其下为一些横向呈水平或东倾、延续性较短的反射震相。地堑南部与北部不同的反射震相特征表明,测线控制区域内的汤阴地堑沉积环境与运动特征可能有所差异。汤东断裂为走向NE,倾向NW的铲型正断层,汤西断裂为走向NE、倾向SE的正断层,北部中深层地震反射剖面上未发现该断层的存在。  相似文献   

12.
Modern seismic recording instruments allow precise measurements of the amplitude of reflected signals. Intuitively we would expect that this amplitude information could be used to increase our knowledge of the physical properties of the reflecting earth. The relevant factors defining the amplitude of a reflection signal are: spherical divergence, absorption, the reflection coefficient of the reflecting interface, the cumulative transmission loss at all interfaces above this, and the effect of multiple reflections. Of these factors, three—spherical divergence, the reflection coefficient and the transmission loss—are reasonably clear concepts (though the estimation of transmission loss from acoustic logs caused some difficulties in the hey-day of synthetic seismograms). Absorption still presents considerable problems of detail, but our understanding has increased significantly in recent years. The factor least well understood is undoubtedly the effect of multiple reflections. Multiple paths having an even number of bounces can have the effect of delaying, shaping and magnifying the pulse transmitted through a layered sequence. Simple demonstations of this phenomenon can be made using elementary thin plates, and these can be presented for various synthetic and real sequences of layers. Such demonstrations lead one to explore the relation between the spectrum of the transmitted pulse and the spectrum of the reflection coefficient series. If it were possible to isolate the amplitude and shape variations imposed by absorption within a layer, there would be a chance that this measure of absorption would be useful as a correlatable or diagnostic indication of rock properties. If it were possible to isolate the amplitude and shape variations imposed by multiple reflections, there would be a chance that this measure would be useful as an indication of cyclic sedimentation and of the dominant durations of the sedimentary cycles. However, the separation of these two effects constitutes a formidable challenge. The very difficulty of this separation suggests that it may be opportune to review the quantitative estimates of absorption made by field experiments.  相似文献   

13.
The quantitative estimate of the anomalies in the seismic signal associated with a realistic subsurface hydrocarbon bearing anticline is discussed. This estimate is derived by the calculation of synthetic seismograms for a dissipative media. The reflection amplitude anomaly associated with gas deposit are as much as three times the background level. However, oil deposits give rise to much weaker anomalies. A decrease of the reflection amplitude by about fifty percent occurs for horizons lying below the gas bearing zone, causing a shadow zone. Time sag anomaly is dependent chiefly on the thickness of the gas bearing zone. The appearance of the anomaly on the time section is modified by the structure, and flat reflectors may be seen over the structure owing to the mutual cancellation of the structural rise and the time sag. Flat spots corresponding to the horizontal fluid interfaces may be expected at the edges of the pool only. The presence of hydrocarbon bearing layers causes a diagnostic character in the composite reflection signal which is very clearly seen in the spectra, the cross correlation functions, and the cross-spectra in profiles crossing the hydrocarbon bearing structures. These are potential diagnostic tools for locating hydrocarbon deposits, particularly when entrapped in thin layers.  相似文献   

14.
地震波场的复射线展开   总被引:2,自引:0,他引:2  
利用复射线理论将惠更斯原理推广到复空间,发展出一种简便的波场计算方法.本文以均匀柱面波在平界面上的反射场为例,用复射线展开法计算了时谐地震波及瞬变地震波在部分反射、全反射和临界角过渡区的场分布,并与平面波谱积分结果进行了比较.结果表明,用复射线展开法可以自动平滑过渡区反射场的奇异性,因而该方法可直接处理焦点、焦散线和全反射过渡区等复杂的波场分析与综合问题,并可简化数值计算过程.  相似文献   

15.
The main problem in seismic prospecting is to infer from the observed reflection response the distribution of density and seismic velocity with depth. This process is generally called the inversion of the reflection data. For plane waves propagating through plane parallel stratification, it can be shown that at any depth the ratio between the amplitude of the transmitted and reflected wave satisfies the Riccati equation. Based on this equation we have formulated an iterative inversion method, which is found to be suitable for numerical computations. We have applied this method on synthetic reflection data, and found that it provides a very fast and accurate inversion.  相似文献   

16.
Cachar area, a part of the Assam–Arakan frontal fold belt in the eastern part of India is characterized by a series of narrow anticlines which are dissected by a number of faults including thrusts, and broad synclines. Out of several exposed structures, Badarpur anticline is most important because of its major share of hydrocarbon strikes in Cachar. However, the seismic image of this area so far has been extremely poor for deciding locations for exploratory drilling. The problems in obtaining a good seismic image in this area include proper spread configuration including far and near offsets to capture reflected signals from the complex subsurface, variation of charge-sizes and depths along a line, large variation in statics due to rapid changes in elevation and near surface properties. The presence of coherent noise due to scattering and proper migration are major processing problems. This paper deals with the problem of spread configuration. For areas of complex subsurface, spread configuration obtained through field experiment at one spot may be quite different from other spots, requiring a large number of costly experiments. In this work, therefore, we used numerical modeling to derive an optimum spread configuration. For this, we chose a line passing through four wells on Badarpur anticline, and prepared a depth model using available seismic, surface geological and well data. Continuous synthetic common shot records were generated through oblique incidence ray tracing and were processed to obtain stacks with various spread configurations. Comparison of these stacks of synthetic data shows that end-on shooting configuration is inadequate for acquiring a good seismic image of the complex subsurface. From the present study it is evident that a split–spread configuration with near and far offsets of 50 m and 2400 m, respectively, can bring out optimum seismic images in Badarpur and other anticlines in this area.  相似文献   

17.
平面声波在粗糙界面上的反射特征研究   总被引:4,自引:1,他引:3       下载免费PDF全文
基于有关粗糙界面的Rayleigh假设,讨论了平面声波按余弦规律快速变化的小尺度粗糙界面上的反射特征.研究表明:这类界面与位于该位置的一个过渡地层的作用相当.该过渡层的厚度为粗糙界面的起伏幅度,速度和密度为上下两层介质相应量的平均值.研究了埋藏很深的微粗糙界面所引起的地震反射(绕射)波的频散特性和走时的构成,即包含零炮检距反射时间、正常时差和界面粗糙时差三部分内容.该粗糙时差与空间坐标和时间坐标无关,与绕射波的阶次有关,绕射波尾随在反射波之后以某一固定的时差出现.且只有当界面的粗糙波长与地震波的波长相当时,才能观测到这类绕射波.该结论为粗糙界面地震反射资料的处理方法提供了理论依据.  相似文献   

18.
Currently, the study on the active fault in the land areas is relatively mature, while there is still lack of detection and research on active faults in the sea areas. Marine exploration, which is different from land areas, has a prominent problem due to the existence of strong reflecting interfaces such as water surface and seafloor in the sea, thus the recording is often accompanied by interference of multiples on seafloor reflections. In addition, because of the characteristics of marine seismic exploration, the source exciting in the water and the geophone receiving in the water, ghost wave usually can be recorded simultaneously during the reflected wave propagation. This phenomenon makes it difficult to distinguish the effective waves and the noise, and has always plagued the data and seriously affects the quality of records. In the offshore and other regions of complex structures, such as inclined interfaces, it is difficult to eliminate the interference of multiples accurately by traditional multiples suppression methods, which are based on the horizontal interface assumption. This paper combines the sea area seismic data and its acquisition method, uses simplified model to simulate the multiples based on the time-distance analysis of multiples and their ghost wave in inclined interface. The time-distance characteristics of the multiples and their ghost waves from different interfaces(including the inclined interface)are obtained, and they are consistent with the actual records. The multiples time-distance simulation can help to distinguish the causes of reflected waves, summarize the multiple-wave time-distance characteristics from different interfaces(including inclined interfaces), and analyze the relationship between the characteristics of multiple waves and primary waves. In particular, this simulation has a significant effect on characterizing the internal multiples that are difficult to identify due to inconspicuous periodicity and the multiples of the inclined interface which present the phenomenon that the vertex of the time-distance curve is shifted. On this basis, relying on the time-distance analysis of ghost wave, we analyze the travel time difference characteristics between reflected waves and their accompanying ghost waves. The differences of the travel time characteristics of different orders ghost wave and reflected wave are summarized and the symmetry of the travel time difference between inclined interface and horizontal interface of ghost waves and reflected waves is analyzed. We simulate the distraction of the ghost wave event with the event of the reflected wave and analyze the influence of the ghost wave on the sea area seismic records. These results can improve the practical interpretation of seismic data. At last, the time-distance information is used to synthesize sea area seismic records, which can help us carry out the effective data processing and understand the characteristics of the time-distance and velocity of multiples in different interfaces and the layer artifact caused by multiples. This study combines the time-distance simulation of multiples and their ghost wave with conventional seismic data processing to analyze the pre-stack and post-stack features of multiple waves and their ghost waves in the seismic records of the sea area. The results of this study are conducive to the effective identification of multiples in seismic records in the sea, provide a theoretical basis for multi-wave suppression and prediction, and may facilitate the future study of sea-area seismic activity detection.  相似文献   

19.
地震反射法中薄煤层分辨能力的研究   总被引:16,自引:2,他引:16       下载免费PDF全文
本文以非线性地震层状构造理论为基础,对薄层分辨率进行了研究,研究中采用理论计算和合成地震记录,并用超声地震模型试验加以验证。结果表明,薄层可检测性准则,应是它的复合反射波的振幅不小于围岩中较强反射波的振幅。由Farr提出的可检测薄层的极限厚度为λ/12(λ为薄层中地震波波长),是在薄层顶、底界面反射系数很小时的近似值。指出,当薄层顶、底界面反射系数较大时,如煤层,其可检测性高于Farr定义,且随界面反射系数的增大或频率的提高而提高;适当地提高地震反射波的频率,就能得到厚2-3m薄煤层的大振幅、高信噪比的反射波。上述结论均用煤田地震勘探实例证实。  相似文献   

20.
南阳盆地特殊震相研究   总被引:1,自引:0,他引:1  
南阳盆地内发生的地震在仪器记录图上与其他地震波形有所不同,在直达横波正常出现后立即开始衰减,紧接着二次出现大振幅、大周期波形。对此进行了分析,证明在南阳地区地壳中存在一个硬度较大的界面,地震波到达该界面时发生反射,反射波与直达波叠加,即出现了第二个较大振幅的波形。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号