共查询到20条相似文献，搜索用时 15 毫秒
1.
A fast and global two point low storage optimization technique for tracing rays in 2D and 3D isotropic media 总被引：1，自引：0，他引：1
We present the problem of tracing rays in 2D and 3D heterogeneous isotropic media as a set of optimization problems. Each optimization problem is obtained by applying Fermat's principle to an approximation of the travel time equation from a fixed source to a fixed receiver. We assume a piecewise linear ray path that simplifies the computations of the problem, in the same way Mao and Stuart suggested in a very recent paper. Here, instead, the reflector geometry and the velocity function are computed by using nonuniformly biharmonic splines. On the other hand, to solve the optimization problem we use the Global Spectral Gradient method. This recent developed optimization scheme is a low storage optimization technique that requires very few floating point operations. It only requires the gradient of the travel time function, and it is global because it converges independently of the initial guess, that is, it does not require a close initial ray path. These three properties of the optimization method and the assumption of piecewise linear rays make this ray tracing scheme a very fast, global and effective method when estimating velocities via tomography. Moreover, in a homogeneous stratified or dipped media, any solution of the optimization problem is the best solution, i.e., it is the global minimum, no matter what numerical approach is used. We present some numerical results that show the computational advantages and the performance of this ray tracing in homogeneous and heterogeneous media. 相似文献
2.
A new scheme for fast calculation of seismic traveltimes and ray paths in heterogeneous media 总被引：2，自引：0，他引：2
IntroductionBoth traveltimes and ray paths are vital information for seismic theoretical research and practice such as stack migration, traveltime inverse, calculation of covering time and so on (Xu, et al,1992). At present, they are obtained usually by either ray shooting (Liu, el al, 1986) or finite difference solution of the eikonal equation (Vidale, 1988; Zhang, et al, 1996). The ray shootingmethod can be understood and programmed easily. However, it is difficult for the method to treatc… 相似文献
3.
In order to trace a ray between known source and receiver locations in a perfectly elastic medium, the take-off angle must be determined, or equialently, the ray parameter. In a viscoelastic medium, the initial value of a second angle, the attenuation angle (the angle between the normal to the plane wavefront and the direction of maximum attenuation), must also be determined. There seems to be no agreement in the literature as to how this should be done. In computing anelastic synthetic seismograms, some authors have simply chosen arbitrary numerical values for the initial attenuation angle, resulting in different raypaths for different choices. There exists, however, a procedure in which the arbitrariness is not present, i.e., in which the raypath is uniquely determined. It consists of computing the value of the anelastic ray parameter for which the phase function is stationary (Fermat's principle). This unique value of the ray parameter gives unique values for the take-off and attenuation angles. The coordinates of points on these stationary raypaths are complex numbers. Such rays are known as complex rays. They have been used to study electromagnetic wave propagation in lossy media. However, ray-synthetic seismograms can be computed by this procedure without concern for the details of complex raypath coordinates. To clarify the nature of complex rays, we study two examples involving a ray passing through a vertically inhomogeneous medium. In the first example, the medium consists of a sequence of discrete homogeneous layers. We find that the coordinates of points on the ray are generally complex (other than the source and receiver points which are usually assumed to lie in real space), except for a ray which is symmetric about an axis down its center, in which case the center point of the ray lies in real space. In the second example, the velocity varies continuously and linearly with depth. We show that, in geneneral, the turning point of the ray lies in complex space (unlike the symmetric ray in the discrete layer case), except if the ratio of the velocity gradient to the complex frequency-dependent velocity at the surface is a real number. We also present a numerical example which demonstrates that the differences between parameters, such as arrival time and raypath angles, for the stationary ray and for rays computed by the above-mentioned arbitrary approaches can be substantial. 相似文献
4.
5.
地震射线追踪是地震定位、层析成像、偏移等领域的重要正演环节.随着这些领域研究的深入,针对传统的网格结构和层状结构在描述复杂地质模型遇到的很大困难,我们采用大小不等、形状各异的地质块组成的集合体来描述三维复杂地质模型,并用三角形面片来描述地质块之间的物性间断面,理论上可以描述任意复杂的地质模型.为适应任意非均匀速度分布的地质模型,基于费马原理,本文发展了与之相适应的逐段迭代射线追踪方法.该方法属于弯曲法范畴,对路径点采用一阶显式增量修正,相对于传统的迭代法,高效省时.数值试验表明,联合逐段迭代法和伪弯曲法的射线追踪扰动修正方案在三维复杂非均匀块状模型中有适用性和高效性. 相似文献
6.
三维场地波动传播的快速射线追踪法 总被引：7，自引：0，他引：7
本文基于Snell定律和Fermat原理对三维任意界面情况下的两点间射线追踪问题进行了研究，从同一条射线满足相同的射线参数出发，推得一个适用于任意界面情况下计算反（折）射点的一阶近似公式。结合迭代技术，给出了三维场地条件下射线追踪迭代算法的计算格式，并进行了三维场地射线追踪模拟计算。计算表明：计算速度相当快，且其计算精度可以根据需要满足要求。 相似文献
7.
地震波CT成果具有图像直观可靠、信息量丰富及适用性强等优点。实践成果表明,地震波CT技术在工程地质勘查、建筑物无损检测、大坝安全检测以及防渗墙质量检测等方面都有良好的探测效果。SIRT算法是地震波CT理论较为成熟的算法。本文提出一种基于Eikonal解的思想,费马原理,惠更斯原理及弯曲波前假设的有限元分析的FDM算法的Matlab计算机程序。该算法程序克服了SIRT算法的缺点,可以成功地处理任意变化的速度、高速对比、尖锐的边界和任意的测量布局。给出几个数值模拟实例进行实践,将两者结合在一起,前后对比,得出结论,并指导工程实践应用。 相似文献
8.
矩形网格三点Fermat射线追踪技术 总被引：4，自引：2，他引：2
矩形网格三点Fermat射线追踪法是基于矩形网格三点扰动法的一种提高计算速度的方法．取矩形网格三个点,在Fermat最小旅行时原则下求取扰动中间点的位置,而不象扰动法那样依次扰动．因此,计算速度比扰动法提高2倍多,同时不受扰动摄动量大小选择的困扰．该方法继承了矩形网格三点扰动法的优点,对任意离散的速度场,总能找到最短时间路径,避免了射线盲区和追踪路径并非时间最短路径问题． 相似文献
9.
地震波走时和射线的有限差分计算 总被引：5，自引：0，他引：5
以往都是采用射线追踪的方法计算地震波的走时和射线，但是当速度模型复杂时这种方法存在一些问题。本文提出另一种计算地震波走时和射线的方法。该方法从程函方程出发，利用互换原理和Ｆｅｒｍａｔ原理计算出各种波的到时和射线。解决了射线追踪方法存在的问题。为地震波走时和射线的计算以及地震波走时反演开辟了一条新途径。 相似文献
10.
分别在湖北省药姑山和九宫山的山头上开展GPS掩星观测实验，成功获取山基掩星观测数据，对掩星事件进行了分析和统计.给出利用山基掩星观测数据反演大气折射指数剖面和电波弯曲角的原理和算法.利用山基GPS掩星观测模拟数据，对反演方法进行试算和检验，结果表明反演方法准确可行.将该反演方法应用于观测数据的反演，获得了观测点高度以下的大气折射率剖面，以及电波弯曲角.实验结果和原理研究表明，山基掩星观测技术是一种潜在的低层大气环境监测新技术. 相似文献
11.
The behaviour of the actual polarization of an electromagnetic wave or elastic S–wave is described by the coupling ray theory, which represents the generalization of both the zero–order isotropic and anisotropic ray theories and provides continuous transition between them. The coupling ray theory is usually applied to anisotropic common reference rays, but it is more accurate if it is applied to reference rays which are closer to the actual wave paths. In a generally anisotropic or bianisotropic medium, the actual wave paths may be approximated by the anisotropic–ray–theory rays if these rays behave reasonably. In an approximately uniaxial (approximately transversely isotropic) anisotropic medium, we can define and trace the SH (ordinary) and SV (extraordinary) reference rays, and use them as reference rays for the prevailing–frequency approximation of the coupling ray theory. In both cases, i.e. for the anisotropic–ray–theory rays or the SH and SV reference rays, we have two sets of reference rays. We thus obtain two arrivals along each reference ray of the first set and have to select the correct one. Analogously, we obtain two arrivals along each reference ray of the second set and have to select the correct one. In this paper, we suggest the way of selecting the correct arrivals. We then demonstrate the accuracy of the resulting prevailing–frequency approximation of the coupling ray theory using elastic S waves along the SH and SV reference rays in four different approximately uniaxial (approximately transversely isotropic) velocity models. 相似文献
12.
Traveltime computation by wavefront-orientated ray tracing 总被引：1，自引：0，他引：1
For multivalued traveltime computation on dense grids, we propose a wavefront‐orientated ray‐tracing (WRT) technique. At the source, we start with a few rays which are propagated stepwise through a smooth two‐dimensional (2D) velocity model. The ray field is examined at wavefronts and a new ray might be inserted between two adjacent rays if one of the following criteria is satisfied: (1) the distance between the two rays is larger than a predefined threshold; (2) the difference in wavefront curvature between the rays is larger than a predefined threshold; (3) the adjacent rays intersect. The last two criteria may lead to oversampling by rays in caustic regions. To avoid this oversampling, we do not insert a ray if the distance between adjacent rays is smaller than a predefined threshold. We insert the new ray by tracing it from the source. This approach leads to an improved accuracy compared with the insertion of a new ray by interpolation, which is the method usually applied in wavefront construction. The traveltimes computed along the rays are used for the estimation of traveltimes on a rectangular grid. This estimation is carried out within a region bounded by adjacent wavefronts and rays. As for the insertion criterion, we consider the wavefront curvature and extrapolate the traveltimes, up to the second order, from the intersection points between rays and wavefronts to a gridpoint. The extrapolated values are weighted with respect to the distances to wavefronts and rays. Because dynamic ray tracing is not applied, we approximate the wavefront curvature at a given point using the slowness vector at this point and an adjacent point on the same wavefront. The efficiency of the WRT technique is strongly dependent on the input parameters which control the wavefront and ray densities. On the basis of traveltimes computed in a smoothed Marmousi model, we analyse these dependences and suggest some rules for a correct choice of input parameters. With suitable input parameters, the WRT technique allows an accurate traveltime computation using a small number of rays and wavefronts. 相似文献
13.
We describe the behaviour of the anisotropic–ray–theory S–wave rays in a velocity model with a split intersection singularity. The anisotropic–ray–theory S–wave rays crossing the split intersection singularity are smoothly but very sharply bent. While the initial–value rays can be safely traced by solving Hamilton’s equations of rays, it is often impossible to determine the coefficients of the equations of geodesic deviation (paraxial ray equations, dynamic ray tracing equations) and to solve them numerically. As a result, we often know neither the matrix of geometrical spreading, nor the phase shift due to caustics. We demonstrate the abrupt changes of the geometrical spreading and wavefront curvature of the fast anisotropic–ray–theory S wave. We also demonstrate the formation of caustics and wavefront triplication of the slow anisotropic–ray–theory S wave.Since the actual S waves propagate approximately along the SH and SV reference rays in this velocity model, we compare the anisotropic–ray–theory S–wave rays with the SH and SV reference rays. Since the coupling ray theory is usually calculated along the anisotropic common S–wave rays, we also compare the anisotropic common S–wave rays with the SH and SV reference rays. 相似文献
14.
The coupling ray theory is usually applied to anisotropic common reference rays, but it is more accurate if it is applied to reference rays which are closer to the actual wave paths. If we know that a medium is close to uniaxial (transversely isotropic), it may be advantageous to trace reference rays which resemble the SH–wave and SV–wave rays. This paper is devoted to defining and tracing these SH and SV reference rays of elastic S waves in a heterogeneous generally anisotropic medium which is approximately uniaxial (approximately transversely isotropic), and to the corresponding equations of geodesic deviation (dynamic ray tracing). All presented equations are simultaneously applicable to ordinary and extraordinary reference rays of electromagnetic waves in a generally bianisotropic medium which is approximately uniaxially anisotropic. The improvement of the coupling–ray–theory seismograms calculated along the proposed SH and SV reference rays, compared to the coupling–ray–theory seismograms calculated along the anisotropic common reference rays, has already been numerically demonstrated by the authors in four approximately uniaxial velocity models. 相似文献
15.
Galactic cosmic rays-clouds effect and bifurcation model of the Earth global climate. Part 1. Theory
Vitaliy D. Rusov Alexandr V. Glushkov Vladimir N. Vaschenko Оksana T. Myhalus Yuriy A. Bondartchuk Vladimir P. Smolyar Elena P. Linnik Strachimir Cht. Mavrodiev Boyko I. Vachev 《Journal of Atmospheric and Solar》2010,72(5-6):398-408
16.
Einar Maeland 《Geophysical Prospecting》1997,45(4):641-651
Migration using an erroneous velocity gives a curve along which the energy is smeared. Associated with this curve is the caustic enveloped by the normal rays. It is possible to compensate for an erroneous velocity by a simple modification of the imaging principle. Formulae are derived for the general case when the velocity changes laterally, and the position of the caustic suggests how to modify the imaging principle so as to obtain an estimate of the NMO velocity. A synthetic example is used to illustrate the results of the analysis. 相似文献
17.
用慢度分块均匀正方形模型将介质参数化，仅在正方形单元的边界上设置计算结点，这些结点构成界面网．根据Ｈｕｖｓｅｎｓ和Ｆｅｒｍａｔ原理，由不断扩张、收缩的波前点扫描代替波前面搜索，在波前点附近点的局部最小走时计算中对波前点之间的走时使用双曲线近似，通过比较确定最小走时和相应的次级源位置，记录在以界面网点位置为指针的３个一维数组中．借助这些数组通过向源搜索可计算任意点（包括界面网以外的点）上的全局最小走时和射线路径．这一方法不受介质慢度差异大小限制，占内存少，计算速度较快，适于走时反演和以Ｍａｓｌｏｖ射线理论为基础的波场计算． 相似文献
18.
Aldo L. Vesnaver 《Geophysical Prospecting》1996,44(5):741-760
A new method to trace rays in irregular grids based on Fermat's principle of minimum time is introduced. Besides the usual transmitted and reflected waves, refracted, diffracted and converted waves can also be simulated. The proposed algorithm is fast and stable, and respects the reciprocity principle between source and receiver better than procedures adopting the shooting method. It is particularly suited to form part of a traveltime inversion procedure. The use of irregular grids allows adaptation of the earth discretization to the available acquisition geometry and ray distribution, to obtain more stable and reliable tomographic images. 相似文献
19.
用于势场反演的特殊解法 总被引：1，自引：0，他引：1
本文以无旋场的积分路径为射线，以边界上已知的势值为输入、输出量、利用反拉东变换求一个梯度场的势函数，用反演的方法求解可以化为拉普拉斯方程的一类二阶偏微分方程。其思路新颖，原理独特，运算正确，边界条件处理方便，数值实现容易，在重力场和地电场的研究中可能得到应用。 相似文献
20.
Diffraction and anelasticity problems involving decaying, “evanescent” or “inhomogeneous” waves can be studied and modelled using the notion of “complex rays”. The wavefront or “eikonal” equation for such waves is in general complex and leads to rays in complex position-slowness space. Initial conditions must be specified in that domain: for example, even for a wave originating in a perfectly elastic region, the ray to a real receiver in a neighbouring anelastic region generally departs from a complex point on the initial-values surface. Complex ray theory is the formal extension of the usual Hamilton equations to complex domains. Liouville's phase-space-incompressibility theorem and Fermat's stationary-time principle are formally unchanged. However, an infinity of paths exists between two fixed points in complex space all of which give the same final slowness, travel time, amplitude, etc. This does not contradict the fact that for a given receiver position there is a unique point on the initial-values surface from which this infinite complex ray family emanates.In perfectly elastic media complex rays are associated with, for example, evanescent waves in the shadow of a caustic. More generally, caustics in anelastic media may lie just outside the real coordinate subspace and one must trace complex rays around the complex caustic in order to obtain accurate waveforms nearby or the turning waves at greater distances into the lit region. The complex extension of the Maslov method for computing such waveforms is described. It uses the complex extension of the Legendre transformation and the extra freedom of complex rays makes pseudocaustics avoidable. There is no need to introduce a Maslov/KMAH index to account for caustics in the geometrical ray approximation, the complex amplitude being generally continuous. Other singular ray problems, such as the strong coupling around acoustic axes in anisotropic media, may also be addressed using complex rays.Complex rays are insightful and practical for simple models (e.g. homogeneous layers). For more complicated numerical work, though, it would be desirable to confine attention to real position coordinates. Furthermore, anelasticity implies dispersion so that complex rays are generally frequency dependent. The concept of group velocity as the velocity of a spatial or temporal maximum of a narrow-band wave packet does lead to real ray/Hamilton equations. However, envelope-maximum tracking does not itself yield enough information to compute synthetic seismogramsFor anelasticity which is weak in certain precise senses, one can set up a theory of real, dispersive wave-packet tracking suitable for synthetic seismogram calculations in linearly visco-elastic media. The seismologically-accepiable constant-Q rheology of Liu et al. (1976), for example, satisfies the requirements of this wave-packet theory, which is adapted from electromagnetics and presented as a reasonable physical and mathematical basis for ray modelling in inhomogeneous, anisotropic, anelastic media. Dispersion means that one may need to do more work than for elastic media. However, one can envisage perturbation analyses based on the ray theory presented here, as well as extensions like Maslov's which are based on the Hamiltonian properties. 相似文献