共查询到20条相似文献,搜索用时 31 毫秒
1.
2.
3.
从含Thomsen各向异性参数的qP波相速度表示式出发,建立并求解三维VTI介质中的频散方程,得到三维VTI介质中的相移算子,进而将以相移算子为基础的最优分裂Fourier方法推广到三维VTI介质,发展了一个三维VTI介质的深度偏移方法.文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中.文中提出的偏移算法对相移法引入了高阶校正项来补偿介质横向变化的影响,使该方法可应用于横向非均匀VTI介质的陡角度成像,文中给出的偏移脉冲响应很好地证明了这一点. 相似文献
4.
Global SH-wavefield calculation for a two-dimensional whole-Earth model with the parallel hybrid PSM/FDM algorithm 下载免费PDF全文
We present a parallel hybrid algorithm based on pseudospectral method (PSM) and finite difference method (FDM) for two-dimensional (2-D) global SH-wavefield simulation. The whole-Earth model is taken as a cross section of spherical Earth, and corresponding wave equations are defined in 2-D cylindrical coordinates. Spatial derivatives in the wave equations are approximated with efficient and high accuracy PSM in the lateral and high-order FDM in the radial direction on staggered grids. This algorithm allows us to divide the whole-Earth into sub-domains in radial direction and implement efficient parallel computing on PC cluster, while retains high accuracy and efficiency of PSM in lateral direction. A transformation of moment tensor between 3-D spherical Earth and our 2-D model was proposed to give corresponding moment tensor components used in 2-D modeling. Comparison of modeling results with those obtained by direct solution method shows very good accuracy of our algorithm. We also demonstrate its feasibility with a lateral heterogeneous whole-Earth model with localized velocity perturbation. 相似文献
5.
L.A. Tabarovsky M.M. Goldman M.B. Rabinovich K.-M. Strack 《Journal of Applied Geophysics》1996,35(4):261-284
Understanding, using, or eliminating three-dimensional (3-D) effects in electromagnetic methods of geophysics are critical requirements. Numerous achievements in 3-D modeling sometimes give the impression that they are widely available today in geophysical practice. This is not necessarily true. Existing 3-D modeling packages prove that we know how to perform 3-D modeling. However, the computer resources and costs involved make the practical application of 3-D EM modeling in geophysical applications very limited.A practical compromise, or even alternative, is represented by 2.5-D modeling characterized by the use of a 3-D source in a 2-D medium. This combination allows one to mathematically describe the related boundary value problem as a sequence of independent two-dimensional problems. The typical technique leading to such a split formulation is Fourier analysis. That is why the individual terms of a split solution are often referred to as harmonics.Although each independent problem is two-dimensional, the algorithmic implementation of finite differences or integral equations for the higher harmonics has some specific features not present in the classical 2-D cases. In this paper, a hybrid scheme consisting of a combination of the finite difference technique with the integral equation approach for transient fields is described. Evaluation of algorithm accuracy is presented and a transient logging technique application is considered. The algorithm is fast and easily implemented on a personal computer 相似文献
6.
对三维电阻率反演问题进行了深入研究,提供了一种利用地表观测数据实现三维反演的实用算法.该方法应用有限差分求正演解,并通过对粗糙度矩阵元素进行适当改进,使之适用于各种情况下粗糙度矩阵的求取,进而建立在模型的总粗糙度极小条件下的反演方程.对反演方程采用收敛速度快且稳定的最小二乘正交分解(LSQR)法进行迭代求解,在迭代求解过程中只需利用偏导数矩阵和其转置矩阵乘以一个向量的结果,回避了直接求偏导数矩阵的繁琐计算,节省了内存,加快了反演的计算速度.不同的计算实例表明上述方法是求解大规模三维电阻率反演问题的有效方法. 相似文献
7.
三维地震数据可视化原理及方法 总被引:5,自引:0,他引:5
可视化是实现三维地震数据真三维解释的有效方法,本文介绍了三维地震数据可视化的基本原理,论述了用光线投射算法实现可视化的主要步骤,详细推导了合成可视化图像的计算公式,分析了阻光度曲线的物理意义及对可视化图像的调节作用,并用两个实际三维地震数据进行了试验.结果表明该方法能够由原始三维地震数据直接对地质体进行成像,能将分散的、孤立的信息互相联系起来,揭示隐藏在数据中的地质现象和规律. 相似文献
8.
基于稳定电流场的基本方程、三维区域满足的边值问题以及相应的变分问题,研究了三维起伏地形条件下电阻率的有限单元数值模拟算法. 离散积分区域时,以三棱柱为最小研究单元,推导了含有地形特征信息的三线性插值型函数以及单元刚度矩阵. 采用变带宽、一维数组方式只存储稀疏刚度矩阵中非零元素,能够节约内存. 利用Cholesky分解法只分解一次大型稀疏矩阵,通过回代可以求出方程组的全部解,当求解有多个供电点的测深问题时可以缩短计算时间. 模型计算表明,在水平层状介质模型上,三维计算结果与解析解或二维数值解十分吻合,计算精度满足误差要求. 在二维山脊上的二极剖面或三维山谷上的中间梯度剖面上,其三维计算结果与相应模型的土槽实验结果或边界元法计算结果也非常接近. 相似文献
9.
10.
为提高水平层状介质中三维电磁波散射和逆散射数值模拟的效率,在对角张量近似(DTA)的基础上根据不同回代方式得到了求解积分方程的DTA1和DTA2两种近似. 这两种近似可以作为计算积分方程稳定型双共轭梯度快速Fourier变换(BCGS-FFT)算法的初始猜测值和预条件因子,从而形成效率更高的混合DTA-BCGS算法. 散射实例说明了DTA2的高精度和混合DTA-BCGS算法尤其是混合DTA2-BCGS算法的高效率. 由于DTA2近似程度更高,将DTA2与变型Born迭代反演方法(DBIM)相结合形成了一种对三维异常体进行重构的快速电磁波逆散射技术. 文中的逆散射实例说明所开发的逆散射技术对重构水平层状介质中的任意三维异常体是非常有效的. 相似文献
11.
A three-dimensional (3-D) kinematic migration algorithm for media in which migration velocity varies linearly with depth is developed, implemented and tested. The algorithm is based on the concept that a single reflection or diffraction in a (zero- or finite-offset) trace may have originated at any point on a constant traveltime surface within the Earth defined by the observed two-way traveltime. The envelope of all such constant time surfaces, for all observed reflections and diffractions produced by one reflector, is the desired migrated 3-D image. The optimal envelope position in depth is determined, beneath each point on a regular grid, by a statistical imaging condition; an incremental function of depth containing the number of constant time surfaces passing through that depth increment is cross-correlated with a Gaussian function whose width is chosen to correspond to the vertical scale of the features of interest. The numerical procedures are based on the observation that, in a medium in which velocity varies linearly with depth, ray segments are circular so traveltimes can be computed analytically. Also, traveltimes are independent of azimuth so the 3-D problem can be collapsed into an equivalent 2-D problem. The algorithm is illustrated and tested by application to synthetic data and to scale-model data from the Seismic Acoustics Laboratory at the University of Houston. 相似文献
12.
13.
14.
15.
地壳三维构造反演和速度层析成像 总被引:5,自引:1,他引:5
本文给出利用空间深地震测深反射资料重建地壳三维构造和速度分布的方法,适用于任意分层并且每层可由若干断块组成的地壳模型。通过采用适当的反演方法获得地壳界面的三维构造,并在此结果基础之上,进一步将剩余走时残差归因于地壳速度的不均匀性,采用东分块的反演方法,重建地壳三维速度分布,数值模拟的结果表明算法是有效的。 相似文献
16.
17.
Multidimensional GPR array processing using Kirchhoff migration 总被引:1,自引:0,他引:1
Mark L. Moran Roy J. Greenfield Steven A. Arcone Allan J. Delaney 《Journal of Applied Geophysics》2000,43(2-4)
We compare the ability of several practical ground-penetrating radar (GPR) array processing methods to improve signal-to-noise ratio (SNR), increase depth of signal penetration, and suppress out-of-plane arrivals for data with SNR of roughly 1. The methods include two-dimensional (2-D) monostatic, three-dimensional (3-D) monostatic, and 3-D bistatic Kirchhoff migration. The migration algorithm is modified to include the radiation pattern for interfacial dipoles. Results are discussed for synthetic and field data. The synthetic data model includes spatially coherent noise sources that yield nonstationary signal statistics like those observed in high noise GPR settings. Array results from the model data clearly indicate that resolution and noise suppression performance increases as array dimensionality increases. Using 50-MHz array data collected on a temperate glacier (Gulkana Glacier, AK), we compare 2-D and 3-D monostatic migration results. The data have low SNR and contain reflections from a complex, steeply dipping bed. We demonstrate that the glacier bed can only be accurately localized with the 3-D array. In addition, we show that the 3-D array increases SNR (relative to a 2-D array) by a factor of three. 相似文献
18.
Leiv -J. Gelius 《Studia Geophysica et Geodaetica》2012,56(2):433-456
An accurate, fast, and simple algorithm for 3-D acoustic modelling of seismic edge diffractions, originally developed in the 1980s, is revisited in this paper. The main objective is to reintroduce this simple approach to edge-diffraction modelling and for the first time give the details of the theory in the open literature. The method is based on a combination of Kirchhoff theory and uniform asymptotic techniques developed within a high-frequency assumption. The diffraction contributions are then computed at stationary edge points only, by analogy with the geometrical ray contributions associated with internal stationary points or specular points. To be able to handle sampling inaccuracies of the critical edge points, a modified algorithm is proposed. Its robustness is verified in case of scattering from a circular edge. Also the extension from rigid or free boundary conditions to the case of edges defined by two penetrable surfaces is discussed in this paper. Both experimental and synthetic 3-D data are presented to demonstrate the potential of this edge-diffraction modelling technique. Since all parameters needed in the computations are obtained from dynamic ray tracing, the algorithm can readily be incorporated in existing software packages for 3-D seismic ray modelling. 相似文献
19.
用可视化技术分析虚拟现实环境下三维地震数据场的研究 总被引:2,自引:0,他引:2
目的 用可视化方法分析三维地震数据,直接观察数据中隐藏的地质现象。方法:将第每个反射子波用一个粒子模型表示,来自同一反射面的子波粒子相互叠合成一个同相面。为加快显示速度,将一个三维地震数据体转换为三维二值图像,用有序扩展算法追踪同相面所在的三维连通域,细化后得到反射面的坐标数据。再在虚拟现实系统下的对反射面进行立体显示,直观地分析地下地层的结构。结果 实际数据试验得到了反射面的三维形态及断层的分布规律,立体显示更深刻地揭示了复杂地质现象及其规律。结论 只有在虚拟现实系统中进行三维地震数据的可视化,才能直观地看到地下地层的真实场景,真正实现三维数据的全三维解释,所提出的方法是有效。 相似文献
20.
提出了针对轴对称二维井间模型的一维、二维联合反演方法.该方法基于纵向成层背景地层的Green函数,以一维反演结果作为联合反演的迭代初始值,同时反演一维背景地层和二维异常剖面的电导率.采用递推技术计算Green函数的系数,可以很方便地同时得到Green函数对空间坐标的积分和对背景层电导率的微分,大大加快了计算雅可比矩阵〖WTHX〗M〖WTBZ〗的速度并使计算方便、准确.反演实例说明了该方法的有效性.反演实例还显示,背景层电导率较二维成像剖面的电导率收敛速度要快,测量数据的精度对二维剖面电导率成像质量的影响大于对背景层电导率分布的影响. 相似文献