共查询到20条相似文献，搜索用时 30 毫秒
1.
2.
Modified image algorithm to simulate seismic channel waves in 3D tunnel model with rugged free surfaces 总被引：1，自引：0，他引：1 
下载免费PDF全文
下载免费PDF全文 The existence of rugged free‐surface three‐dimensional tunnel conditions in the coal seams, caused either by geological or mining processes, will inevitably influence wave propagation characteristics when the seismic waves go through the coal mines. Thus, a modified image algorithm has been developed to account for seismic channel waves propagating through this complicated topography with irregular free surfaces. Moreover, the seismic channel waves commonly exhibit damped and dispersive signatures, which is not only because of their own unique sandwich geometry of rock–coal–rock but also because of the viscoelastic behavior of coal. Considering the complexity of programming in three‐dimensional tunnel models with rugged free surfaces, an optimized vacuum grid search algorithm, enabling to model highly irregular topography and to compute efficiently, is also proposed when using high‐order staggered finite‐difference scheme to simulate seismic channel wave propagations in viscoelastic media. The numerical simulations are implemented to investigate the accuracy and stability of the method and the impact of coal's viscoelastic behavior on seismic channel wave propagation characteristics. The results indicate that the automatic vacuum grid search algorithm can be easily merged into high‐order staggered finite‐difference scheme, which can efficiently be applied to calculate three‐dimensional tunnel models with rugged free surfaces in the viscoelastic media. The simulation also suggests that the occurrence of a three‐dimensional tunnel with free surfaces has a remarkable influence on the seismic channel wave propagation characteristics and elastic energy distribution. 相似文献
3.
本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础，推导出了曲线坐标系下对应的方程，然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程，并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件，实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面，地表的网格线紧贴地形，避免了台阶近似造成的数值散射.数值模拟结果表明，在双相介质起伏自由地表和分界面处，各类波型复杂的反射透射规律可以清晰展现，曲线网格有限差分法可以精确地解决地震波在含起伏地表的双相各向同性介质中的传播问题. 相似文献
4.
We have pursued two-dimensional (2D) finite-difference (FD) modelling of seismic scattering from free-surface topography. Exact free-surface boundary conditions for the particle velocities have been derived for arbitrary 2D topographies. The boundary conditions are combined with a velocity–stress formulation of the full viscoelastic wave equations. A curved grid represents the physical medium and its upper boundary represents the free-surface topography. The wave equations are numerically discretized by an eighth-order FD method on a staggered grid in space, and a leap-frog technique and the Crank–Nicholson method in time.
In order to demonstrate the capabilities of the surface topography modelling technique, we simulate incident point sources with a sinusoidal topography in seismic media of increasing complexities. We present results using parameters typical of exploration surveys with topography and heterogeneous media. Topography on homogeneous media is shown to generate significant scattering. We show additional effects of layering in the medium, with and without randomization, using a von Kármán realization of apparent anisotropy. Synthetic snapshots and seismograms indicate that prominent surface topography can cause back-scattering, wave conversions and complex wave patterns which are usually discussed in terms of inter-crust heterogeneities. 相似文献
In order to demonstrate the capabilities of the surface topography modelling technique, we simulate incident point sources with a sinusoidal topography in seismic media of increasing complexities. We present results using parameters typical of exploration surveys with topography and heterogeneous media. Topography on homogeneous media is shown to generate significant scattering. We show additional effects of layering in the medium, with and without randomization, using a von Kármán realization of apparent anisotropy. Synthetic snapshots and seismograms indicate that prominent surface topography can cause back-scattering, wave conversions and complex wave patterns which are usually discussed in terms of inter-crust heterogeneities. 相似文献
5.
中国大陆中西部普遍具有强烈的地形起伏,起伏地形会对地震资料的处理分析产生严重干扰.精细处理起伏地形成为高精度地震成像的必然要求.传统方法通过填充低速介质将不规则模型扩展为规则模型来处理起伏地形.近年来,借助坐标变换将物理空间不规则模型转换为计算空间规则模型的地形平化方法,为解决起伏地形问题提供了新思路.本文基于经典的模型扩展和新发展的地形平化方法分别处理起伏地形,从走时正演、射线追踪和反演成像三个方面,全面细致地评判了两种地形处理方法在起伏地形层析成像中的适用性和有效性.结果表明,模型扩展中阶梯状近似和填充介质速度参与计算,会造成起伏地形走时计算精度损失,出现虚假射线路径和错误出射角,导致反演分辨率降低,成像结果模糊甚至失真;地形平化中采用贴体网格参数化,能够保证离散模型完全匹配起伏地形,并且保持起伏地形在物理空间和计算空间中均为自由表面.在此基础上发展的层析成像技术具有高度的保真性,有效地处理了地形起伏效应,为起伏地形区域精细速度成像提供了有力的技术保障. 相似文献
6.
7.
New alternative formulations of exact boundary conditions for arbitrary three-dimensional (3D) free-surface topographies on seismic media have been derived. They are shown to be equivalent to previously published formulations, thereby verifying the validity of each set of formulations. The top of a curved grid represents the free-surface topography while the interior of the grid represents the physical medium. We assume the velocity–stress version of the viscoelastic wave equations to be valid in this grid before transforming the equations to a rectangular grid. In order to perform the numerical discretization we apply the latter version of the equations for seismic wave propagation simulation in the medium. The numerical discretization of the free-surface topography boundary conditions by second-order finite differences (FDs) is shown, as well as the spatially unconditional stability of the resulting system of equations. The FD order is increased by two for each point away from the free surface up to eight, which is the order used in the interior. We use staggered grids in both space and time and the second-order leap-frog and Crank– Nicholson methods for wavefield time propagation. An application using parameters typical of teleseismic earthquakes and explosions is presented using a 200 × 100 km2 area of real topography from southwestern Norway over a homogeneous medium. A dipping plane wave simulates a teleseismic P-wave incident on the surface topography. Results show conversion from P- to Rg- (short period fundamental mode Rayleigh) waves in the steepest and/or roughest topography, as well as attenuated waves in valleys and fjords. The codes are parallelized for simulation on fast supercomputers and PC-clusters to model high frequencies and/or large areas. 相似文献
8.
9.
10.
We have implemented a 3D finite‐difference scheme to simulate wave propagation in arbitrary anisotropic media. The anisotropic media up to orthorhombic symmetry were modelled using a standard staggered grid scheme and beyond (monoclinic and triclinic) using a rotated staggered grid scheme. The rationale of not using rotated staggered grid for all types of anisotropic media is that the rotated staggered grid schemes are more expensive than standard staggered grid schemes. For a 1D azimuthally anistropic medium, we show a comparison between the seismic data generated by our finite‐difference code and by the reflectivity algorithm; they are in excellent agreement. We conducted a study on zero‐offset shear‐wave splitting using the finite‐difference modelling algorithm using the rotated staggered grid scheme. Our S‐wave splitting study is mainly focused on fractured media. On the scale of seismic wavelenghts, small aligned fractures behave as an equivalent anisotropic medium. We computed the equivalent elastic properties of the fractures and the background in which the fractures were embedded, using low‐frequency equivalent media theories. Wave propagation was simulated for both rotationally invariant and corrugated fractures embedded in an isotropic background for one, or more than one, set of fluid‐filled and dry fractures. S‐wave splitting was studied for dipping fractures, two vertical non‐orthogonal fractures and corrugated fractures. Our modelling results confirm that S‐wave splitting can reveal the fracture infill in the case of dipping fractures. S‐wave splitting has the potential to reveal the angle between the two vertical fractures. We also notice that in the case of vertical corrugated fractures, S‐wave splitting is sensitive to the fracture infill. 相似文献
11.
A hybrid method combining finite element and 4th-order finite difference techniques is developed to model SH and P-SV seismic wave propagation in a 2D elastic medium with irregular surface topography. Both the classic staggered grid finite difference scheme and the partially staggered grid scheme are tested. The accuracy of the hybrid method is studied by comparison with a semi-analytical and another numerical method. Subsequently, to study the amplification, numerical simulations of seismic wave propagation in a series of hills are carried out and compared with the single-hill case. Depending on the position of the source in relation to the topography, the ratio between the heights and lengths of the hills or the ratio between the lengths of the hills and the wavelength, the presence of several hills as opposed to a single one can increase the amplification effect due to topography. This study highlights the fact that, when evaluating topographic site effects, surrounding topography must be taken into account in addition to local topography. 相似文献
12.
Modelling the wave phenomena in acoustic and elastic media with sharp variations of physical properties using the grid‐characteristic method 下载免费PDF全文
下载免费PDF全文 Alena V. Favorskaya Michael S. Zhdanov Nikolay I. Khokhlov Igor B. Petrov 《Geophysical Prospecting》2018,66(8):1485-1502
This paper introduces a novel method of modelling acoustic and elastic wave propagation in inhomogeneous media with sharp variations of physical properties based on the recently developed grid‐characteristic method which considers different types of waves generated in inhomogeneous linear‐elastic media (e.g., longitudinal, transverse, Stoneley, Rayleigh, scattered PP‐, SS‐waves, and converted PS‐ and SP‐waves). In the framework of this method, the problem of solving acoustic or elastic wave equations is reduced to the interpolation of the solutions, determined at earlier time, thus avoiding a direct solution of the large systems of linear equations required by the FD or FE methods. We apply the grid‐characteristic method to compare wave phenomena computed using the acoustic and elastic wave equations in geological medium containing a hydrocarbon reservoir or a fracture zone. The results of this study demonstrate that the developed algorithm can be used as an effective technique for modelling wave phenomena in the models containing hydrocarbon reservoir and/or the fracture zones, which are important targets of seismic exploration. 相似文献
13.
针对山地地震勘探数据低信噪比问题,近地表散射波分离意义显得尤为突出,地震干涉测量法为此提供了一种技术手段.本文将地震干涉测量理论和散射理论结合起来,导出了近地表散射波地震干涉测量表达式,分为互相关型和褶积型表达式,它们由实际波场和背景波场干涉测量构成.根据近地表散射波分离理论,结合陆上地震勘探实际观测系统,采用褶积和反褶积混合型地震干涉测量配置,用实际地震资料展示了近地表散射波分离技术的应用效果.经过理论分析和砾石区实际资料试验,表明地震干涉测量不仅能分离测线上散射源产生的散射波,而且能分离部分侧面散射波.该技术的优点在于它适应于起伏地形和不均匀近地表结构,并且不需要起伏地形和近地表速度信息.为了从实际资料中消除近地表散射波,本文采用多道匹配滤波自适应减法,在砾石区见到较好效果. 相似文献
14.
一般而言,由于地下非均质体的存在所产生的二次波源,由它再生成新的波场,叫散射波场。目前,将散射波作为有效波来成像,已开始在溶洞和裂缝等特殊地质体的识别中得到应用。但对于野外资料采集来说,地表复杂地区,如戈壁、砾石区和山前带,大量存在的散射波却是干扰波,它们的存在会严重影响资料的品质,而其研究与实际应用国内外还很少。因此,通过正演模拟,分析散射波的基本特征,在此基础上研究散射波的去噪方法显得十分必要。本文从地震波运动学时距关系出发,研究了反射波和散射波的几何特征;然后用有限差分正演,模拟了散射波场,用理论模型研究并测试了局域双曲线Radon变换散射波去噪新方法。对于实际炮集资料,分析了F-K滤波方法压制散射噪声的局限,采用局域双曲Radon变换有效地去除了炮集中存在的散射噪声,取得了较好的应用效果。 相似文献
15.
地球深部圈层及沉积盆地是一种分区非均匀介质系统,其中不规则地层边界(含起伏地表)对地震波的主要特征有显著影响,而地层的随机非均匀性则主要影响地震波的散射和衰减特征.为了精确刻画不规则地层边界对地震波的反射、透射效应以及非均质体散射引起的地震波衰减效应,全局广义R/T递推传播矩阵法(GGRTM)被提出并逐步发展成为继有限元和有限差分方法之后的另一种复杂介质高精度地震波传播半解析求解方法.在已有的此类方法中,不规则边界均匀地层GGRTM法的优势在于对不规则地层边界的反射和透射效应的准确模拟,而非均质地层薄板化GGRTM法则能准确描述非均质体散射对地震波衰减的影响.本文吸收这两种已有方法的优势,提出了一种考虑非均匀介质、不规则边界的全局广义R/T递推传播矩阵混合方法,并将其用于对边界不规则、层内非均质的复杂模型的二维SH波场模拟.随后在本文方法与边界元法对比研究的基础上讨论了方法的模拟精度.研究结果表明本文提出的混合法是一种解决复杂模型高精度地震模拟的有效方法. 相似文献
16.
地震诱导电磁现象是国内外地学领域十分关注的前沿问题,前人对地震波和电磁场耦合波场的认识主要是基于规则模型获得的.为研究含起伏地表和地下界面的地层中震电波场激发、传播特性,本文采用有限元软件COMSOL Multiphysics模拟点震源激发的电磁场.首先给出频率域二维SHTE模式震电耦合方程组,然后利用COMSOL软件建立计算模型,并求解出点力源激发震电波场的频率域响应,最后利用FFT变换得到地震波场和电磁场的时间域波形.模拟结果表明,震电波场中存在三种类型的电磁信号,第一种是震源直接激发的电磁波;第二种是地震波在分界面处激发的电磁波(包括自由表面、地下不同介质分界面);第三种是伴随地震波的同震信号,前两种电磁波比地震波更早到达远处观测台站,对地震预警有重要意义.此外,研究还发现:当地震波传播至地表并沿着地表传播时,在地表附近空气层中同样记录到了伴随地震波传播的电磁扰动信号,该信号与相同水平源距条件下、地下观测点接收到的电磁信号相同,这与前人的一些观测结果相符.本文研究结果为今后地震电磁信号的解释提供了理论证据. 相似文献
17.
The scattering of elastic energy by random clusters of fractures and/or cavities in a massive rock is studied. The interpretation of the scattered seismic response reveals crucial information about the clusters of inhomogeneities (fractures/cavities), which may correspond to reservoirs. The study is based on a new two‐dimensional numerical‐modelling method that relaxes the constraints on the location and orientation of the inhomogeneities, accounts for inhomogeneities that have almost no volume but a finite surface area (fractures) and improves the accuracy of the calculation when the size of the inhomogeneities is comparable to the mesh size. It is shown that the nature of the seismic response of zones of diffuse fracturing and/or cavities is associated with the non‐uniformity of micro‐inhomogeneities in such zones; accumulations of these micro‐inhomogeneities are known as clusters. The relationship between the non‐uniformity of micro‐inhomogeneities and the strength of the seismic response has been established and measured. Considerable differences in the structure of the seismic response of zones of diffuse fracturing and diffuse cavities have been identified. Converted PS‐waves dominate in the scattered wavefield associated with fractures. This is explained, as the modelling results show, by a greater transparency of fluid‐filled fractures, which reduces the reflected energy of compressional waves. The wavefield associated with cavities is characterized by the predominance (in terms of strength) of compressional PP‐waves. The strength of converted PS‐waves in the scattered wavefields for both media is approximately the same. On the whole, according to the results of the modelling, the energy of the scattered response of fractured reservoirs is considerably less (about two times) than that of cavernous reservoirs. 相似文献
18.
Boris Gurevich 《Geophysical Prospecting》2015,63(1):141-150
Naturally fractured reservoirs are becoming increasingly important for oil and gas exploration in many areas of the world. Because fractures may control the permeability of a reservoir, it is important to be able to find and characterize fractured zones. In fractured reservoirs, the wave‐induced fluid flow between pores and fractures can cause significant dispersion and attenuation of seismic waves. For waves propagating normal to the fractures, this effect has been quantified in earlier studies. Here we extend normal incidence results to oblique incidence using known expressions for the stiffness tensors in the low‐ and high‐frequency limits. This allows us to quantify frequency‐dependent anisotropy due to the wave‐induced flow between pores and fractures and gives a simple recipe for computing phase velocities and attenuation factors of quasi‐P and SV waves as functions of frequency and angle. These frequency and angle dependencies are concisely expressed through dimensionless velocity anisotropy and attenuation anisotropy parameters. It is found that, although at low frequencies, the medium is close to elliptical (which is to be expected as a dry medium containing a distribution of penny‐shaped cracks is known to be close to elliptical); at high frequencies, the coupling between P‐wave and SV‐wave results in anisotropy due to the non‐vanishing excess tangential compliance. 相似文献
19.
提出一种新的三维空间不规则网格有限差分方法,模拟具有地形构造的非均匀各向异性介质中弹性波传播过程. 该方法通过具有二阶时间精度和四阶空间精度的不规则交错网格差分算子来近似一阶弹性波动方程,与多重网格不同,无需在精细网格和粗糙网格间进行插值,所有网格点上的计算在同一次空间迭代中完成. 针对具有复杂物性参数和复杂几何特征的地层结构,使用精细不规则网格处理粗糙界面、断层和空间界面等复杂几何构造, 理论分析和数值算例表明,该方法不但节省了大量计算机内存和计算时间,而且具有令人满意的稳定性和精度. 相似文献
20.
An upwind fast sweeping scheme for calculating seismic wave first‐arrival travel times for models with an irregular free surface 下载免费PDF全文
下载免费PDF全文 The topography‐dependent eikonal equation formulated in a curvilinear coordinate system has recently been established and revealed as being effective in calculating first‐arrival travel times of seismic waves in an Earth model with an irregular free surface. The Lax–Friedrichs sweeping scheme, widely used in previous studies as for approximating the topography‐dependent eikonal equation viscosity solutions, is more dissipative and needs a much higher number of iterations to converge. Furthermore, the required number of iterations grows with the grid refinement and results in heavy computation in dense grids, which hampers the application of the Lax–Friedrichs sweeping scheme to seismic wave travel‐time calculation and high‐resolution imaging. In this paper, we introduce a new upwind fast sweeping solver by discretising the Legendre transform of the numerical Hamiltonian of the topography‐dependent eikonal equation using an explicit formula. The minimisation related to the Legendre transform in the sweeping scheme is solved analytically, which proved to be much more efficient than the Lax–Friedrichs algorithm in solving the topography‐dependent eikonal equation. Several numerical experiments demonstrate that the new upwind fast sweeping method converges and achieves much better accuracy after a finite number of iterations, independently of the mesh size, which makes it an efficient and robust tool for calculating travel times in the presence of a non‐flat free surface. 相似文献