共查询到20条相似文献，搜索用时 31 毫秒
1.
Since the wave equation of magnetotelluric(MT)field is similar to the one of seis-mic,the migration techniques used in seismic can be applied to MT data.In thispaper we make use of the principle of reflector mapping(i.e.U/D imaging principle)to image MT data.That is,the MT wavefield observed on the surface of the earth canbe resolved into upgoing and downgoing waves,the waves are extrapolated downward bythe phase-shift method or the phase-shift plus interpolation(PSPI)method.Conductivi-ty interfaces of the medium could be found by using the time coincidence of the upgoingand downgoing waves.Theoretical calculations show that the migration technique of MTdata presented here is very effective.It can not only enhance the lateral resolution of MTdata,but also obtain the visual image of subsurface interfaces.As compared with theconventional 2-D inversion,this procedure is more simple in calculation and can be easi-ly put into practice on a personal computer and is able to obtain the MT depth section,which is 相似文献
2.
利用高阶交错网格有限差分模拟Kelvin-Voigt黏弹性介质中传播的地震波,同时将完全匹配层吸收边界条件引入到其边界处理中。数值模拟结果表明,完全匹配层吸收边界效果好,高阶有限差分能模拟得到的黏弹性介质波场精度较高。对模拟的黏弹性波场进行分析,表明介质的粘滞性使地震反射波的能量变弱,高频衰减明显,并比低频衰减得快,主频向低频方向移动,有效频带变窄,即降低了地震波的分辨率;并且反射转换波比反射纵波要衰减得快;而且还随着传播距离的增加,其峰值频率也逐渐降低。通过数值模拟分析具有不同的粘滞系数介质对地震波的吸收和衰减,结果表明随着粘滞系数的增大,地下介质对地震波的吸收衰减更明显。 相似文献
3.
4.
Asma Hadjadj Zahia Benaissa Abdelkader Benaissa Amar Boudella 《Arabian Journal of Geosciences》2017,10(9):212
The principal objective in Vertical Seismic Profile (VSP) processing is the separation of the downgoing and upgoing wavefields. Several methods have been suggested in this area. This paper presents a new approach based on the gray-scale Hough transform (GSHT) which is an extension of the conventional Hough transform used to detect straight lines and other curves. The technique, we suggest here, directly maps the gray-scale VSP image, including the downgoing and upgoing linear events, in image coordinate space (x,t,g) to the gray Hough parameter counting space (θ,ρ,g), where θ and ρ are the polar parameters and g is the gray-scale value. In this new space, the downgoing events appear in the negative angles θ quadrant and the upgoing in the positive quadrant, owning to their opposite apparent velocities. The inverse GSHT algorithm, we developed in this study, is performed for extracting separately these two wavefields by considering the straight lines that satisfy the corresponding filtering conditions. The experimental results on synthetic VSP datasets are convincing. The wave separation is well performed, even in the presence of loud noise levels, with signal to noise ratio improvement and amplitude preservation, in contrast to median filtering. 相似文献
5.
角度域照明可以更加细致地对地下复杂地质体进行照明分析。提出了一种适用于起伏地表条件下的角度域照明分析方法,其主要过程是利用局部倾斜叠加技术分解波场,对起伏地表处采用加权函数判定局部时窗内的振幅包络位置,以区分其位于地表之上或之下;并引入相似系数函数计算反假频局部平面波场,根据局部斜率与平面波场入射角关系,得到局部角度域地震波场;根据帕斯瓦尔等式下的时间-角度域照明关系式,计算得出角度域照明强度。通过起伏地表高速侵入矿体模型试算,得到给定观测系统下角度域照明强度与地下分析点照明随炮点、波场入射角的变化关系。此方法可用于起伏地表条件的观测系统优化、叠前AVA分析等方面。 相似文献
6.
井间地震串联滤波波场分离方法 总被引：2，自引：0，他引：2
严又生 《物探化探计算技术》1997,19(3):210-215
本文根据井间地震波场响应以及多种波场分离方法的性能，采用基于Ｋ－Ｌ变换的特征向量滤波与中值滤波相结合的串联滤波方法分离井间地震观测到的上、下行反射波。数值计算结果表明：在波场分离过程中，由于同时综合考虑了地震波能量主分量分布以及反射波同相轴时差两个参数，使得分离效果明显、处理更为简便，具有滤波特征矢量的选取条件更为宽松、无需精确确定拉平轴的时间等特点。 相似文献
7.
基于波动方程理论,以陇东地层结构为基础建立地质模型,针对页岩层这一重点研究对象,采用Tesseral软件进行有限差分法全波场正演技术,模拟计算陇东地区较为发育的反射波组,以便分析研究包括页岩层在内的各个反射波组特征。同时结合人工合成地震记录,并与实际地震时间剖面匹配对比,从而标定各个反射波组的地质属性,最终对页岩层反射波做出准确、合理的地质解释。 相似文献
8.
地震照明分析能够定量分析影响逆时偏移成像质量的因素。文章基于Poynting矢量方法进行地震定向照明分析及成像补偿方法研究,利用Poynting矢量进行波场传播角度分解,通过计算不同入射及散射波方向的部分照明建立局部照明矩阵。该矩阵包含了照明分析中波与地质结构相互作用所涉及的所有角度信息,可以方便地分析目标附近所有可能的入射波和散射波对局部照明的贡献。根据局部照明矩阵可以获得对地下结构的照明描述、采集倾角响应,并依此进行成像补偿计算。以SEG/EAGE模型作为数值计算实例,进行基于Poynting矢量的照明分析及成像补偿,有效地提高了成像质量。该方法为角度域逆时偏移成像校正及分辨率分析等提供了高效的处理手段。 相似文献
9.
10.
《Physics and Chemistry of the Earth, Part A: Solid Earth and Geodesy》1999,24(3):231-234
We study some effects that produce 3D topographical irregularities under incoming elastic waves from point sources using the indirect boundary element method. This technique is based on the representation of elastic waves in terms of single-layer boundary sources. In this way reflected and diffracted waves are constructed at the boundaries from where they are radiated by means of boundary sources. The field emited from the point source is computed with the moment tensor for a shear dislocation and analytical expressions of the elastodynamic 3D Green's functions. In this way we can construct a double couple with variable orientation. We compare the technique with that of Bouchon, who used the discrete wave number method for a shear dislocation in a halfspace with a triangular source function. We subsequently apply our method to simulate the seismic response of a mountain of ellipsoidal geometry. The source function that we have considered is a triangular pulse and we show results corresponding to the velocity registered over the surface of the irregularity and that of the halfspace. This motion is presented by means of snapshots showing the evolution of the wavefields that are present in the problem. 相似文献
11.
基于多相孔隙介质弹性理论,给出了非饱和土中不同弹性波的传播方程。根据分界面上的边界条件,建立了各势函数波幅值之间的关系式,讨论了入射剪切波在不同饱和度土层分界面上的反射与透射问题。在无限空间非饱和土体中存在3种压缩波和1种剪切波,因此,当剪切波传播到不同饱和度的非饱和土层分界面上将分别在上、下土层激发产生4种反射波和4种透射波。推导出不同反射波和透射波的振幅比例系数和能量比例系数的理论表达式,并且在此基础上进行数值分析。在数值算例中分别研究了各反射波与透射波的能量比例系数（即能量反射率和能量透射率）受入射频率、入射角度以及上、下土层土体饱和度变化的影响情况。计算结果表明：各能量反射率和能量透射率不仅与入射角和入射频率有关,而且其受上、下土层饱和度变化的影响也同样不能忽视。 相似文献
12.
Xuhui Zheng Jinzhong Sun Leihua Yao Lei Liu Anyang Shao 《Environmental Earth Sciences》2018,77(14):528
This paper presents a three-dimensional oblique incidence method based on stress time-history functions of seismic body waves. This method realizes the three-dimensional oblique incidence of a single body wave phase or superposed multiple phases, which considers the earthquake action in the analysis of the slope seismic response. To implement this method, first, the spatial relationship between slope and incident seismic wave rays is described, and the direction of incident seismic waves is defined with two parameters: the azimuth angle \(\alpha\), which reflects the azimuth of epicenter; and the incident angle \(\theta\), which is a concentrated expression of the hypocentral depth and epicentral distance. Second, at the bottom boundary of the slope numerical model, the vibration starting time of each node excited by an oblique incident seismic wave is calculated. By considering vibration time sequence, the successive disturbance of rock and soil medium at the bottom of the slope is achieved. Then the stress time-history functions of six body wave phases are derived. Numerical simulation results of the x demonstrate the effectiveness, applicability and necessity of the proposed method. 相似文献
13.
纵横波波场分离是弹性波逆时偏移成像中的关键步骤,可以有效消除纵横波串扰和成像畸变,提高成像精度。在各向异性介质中,常用的波场分离法有非平稳滤波器或者低秩近似的方法,但由于使用多次傅里叶变换,导致波场分离计算成本较高。借鉴基于Helmhotz算子构建各向同性介质纵横波解耦方程的思路,提出消除振幅畸变的修正伪Helmholtz分解算子,通过待定系数方法求解横向各向同性介质中P波和S波2种模式的表达式,并将其变换为一阶解耦的伪弹性波方程,实现时空域纵横波波场解耦。通过简单模型的波场分离测试,得到分离后的P波与S波波场,验证了该方法的波场分离的有效性。此外,将解耦得到的矢量P波与矢量S波应用到弹性逆时偏移之中,利用矢量点积互相关成像条件得到清晰的弹性波逆时偏移成像结果,说明本方法在复杂介质中有着较好的适用性,同时可以有效地应用到VTI介质弹性逆时偏移之中。 相似文献
14.
为了实现适应崎岖海底、大陡坡、海底火山等复杂海底地质条件且灵活、稳定、精度高的射线追踪方法,并基于该方法详细分析复杂海底条件下各种地震波型的运动学特征,综合多种算法的优势,实现了一种快速推进迎风插值射线追踪方法。首先,采用混合网格法剖分复杂海底地质模型;其次,通过融入迎风差分思想的线性插值策略来构建精度高、无条件稳定且灵活的局部走时和射线路径计算公式;然后,综合应用这些公式和多级次快速推进法,灵活计算整个模型各种类型地震波的走时,并基于逆向追踪方法计算射线路径;最后,对该方法进行了精度分析,发现其能够获得相对高的走时和射线路径计算精度,且反射波的计算精度远高于入射波。此外,计算实例分析显示,初至波中富含折射波和陡倾构造的反射波在很大偏移距还能被接收,崎岖海底各种波型传播路径复杂;基于此提出加大采集排列长度和采用直达波走时可提高复杂构造成像质量等复杂海底地震数据采集与处理方面的思考与建议。 相似文献
15.
单程波法地震波衰减特征数值模拟分析 总被引：2，自引：1，他引：1
在非均匀粘弹性介质中,地震波衰减研究能够有效揭示地层吸收特性对地震波的影响。这里从单程波方程出发,结合粘弹性介质中的复波数理论,考虑介质的横向变化,采用分步傅立叶波场延拓技术,推导出适用于横向变速介质的带Q地震记录计算公式,最后将此方法应用于四个不同含Q地层模型。试验表明：单程波法计算速度快,波场信息简单清晰,易于深入分析和反演算法的建立。并且,Q值对地震波波形、能量、频率等都产生了较为明显的影响,对实际资料能量补偿、地震资料分辨率的提高,以及利用衰减特征进行储层识别有一定指导价值。 相似文献
16.
利用时间分裂的错格伪谱法模拟地震波在基于改进BISQ模型的双相介质中的传播 总被引：1，自引：0，他引：1
改进的BISQ(Biot-Squirt)模型中各参数具有明确的物理意义和可实现性,在不引入特征喷流长度的情况下可将Biot流动和喷射流动两种力学机制有机地结合起来;而高精度的地震波场数值模拟技术是研究双相介质地震波传播规律的重要手段。本文从本构方程、动力学方程和动力学达西定律出发,推导了基于改进BISQ模型的双相各向同性介质的一阶速度--应力方程组;采用时间分裂错格伪谱法求该方程组的数值解,模拟半空间及层状双相介质中的地震波场。数值模拟结果表明:①与传统方法相比,时间分裂错格伪谱法波场数值模拟的精度更高,压制网格频散效果更好;②在非黏滞相界情况下,慢纵波呈传播性,而在黏滞相界情况下,慢纵波呈扩散性,以静态模式出现在震源位置;③双相介质分界面处,各类波型复杂的反射透射规律可由数值模拟结果清晰展现。 相似文献
17.
Stochastic seismic wave interaction with a slippery rock fault is studied, based on the principle of conservation of momentum
at the wave fronts along the fault. By using the displacement discontinuity method, the wave propagation equations are derived
for incident longitudinal-(P-) and shear-(S-) waves, respectively. This is an extension of the study by Li and Ma (2010) for
blast-induced wave propagation across a linear rock joint. Stochastic seismic waves are generated from a frequency spectrum
and used to analyze the seismic wave interaction with a rock fault having a Coulomb-slip behavior. Parametric studies are
carried out to investigate the effect of the intensity and impinging angle of the incident seismic waves on wave propagation
across a slippery rock fault. Results show that the transmission of the incident P-wave is almost not affected by the fault,
on the contrary, this is not the case for an incident S-wave, due to the occurrence of a relative slip which is related to
the impinging angle of the incident S-wave. A quantitative study is presented which is of help in understanding the propagation
and attenuation laws of seismic waves in discontinuous rock masses. 相似文献
18.
Jun Wang Jie Liu Qiao Liang Xiangsen Ouyang 《Geotechnical and Geological Engineering》2018,36(4):1991-2001
In the soil slope supported by bolts, longitude waves instead of transverse waves, generated by earthquakes, first reach the slope surface. With the dynamic response of the P (pressure) wave along the anchorage structure, first, a theoretical study was conducted to investigate the propagation characteristics of the interference superposition, generated by the SV (shear-vertical) and the P waves. The SV wave was formed by the wave, originating from the bottom and reflected from the free surface of the slope, whereas the latter was the incident P wave, propagating in the slope. In addition, the structural measures, restraining the seismic wave, and the characteristics of the restraint effect at the free segment of the bolt were investigated. According to the wave-way difference between the incident P wave and the reflected SV wave, the minimum critical slope angle, influenced by the interference at the shallow slope, and the maximum influencing depth of the dynamic response, acting vertically to the slope surface, were obtained. The results indicate that the maximum influencing depth linearly correlated with the slope angle. Furthermore, based on the propagation characteristics of the P wave along the bolt, and the coupled relation between the wave length and the anchorage design parameters, the axial acceleration of the wave propagating along the bolt axis was obtained. Then, the theoretical length of the anti-seismic bolt, subjected to seismic waves, and the compensation force of the anchorage structure were obtained. Finally, a numerical study, based on FlAC3D, properly verified the theoretical conclusions. 相似文献
19.
20.
The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations. 相似文献