共查询到20条相似文献，搜索用时 31 毫秒
1.
最短路径射线追踪方法及其改进 总被引：35，自引：9，他引：35
综述了用网络最短路径算法求解地震射线追踪问题的原理、方法技术以及存在问题和改进措施。特别介绍了作者在最短路径算法基础上，提出的动态网络最短路径地震射线追踪方法。该方法先采集从炮点到整个模型所有节点上的初至旅行时，其中，在一个单元内，对相邻每对已计算出最小旅行时的节点进行线性插值，并利用Fermat原理计算未知节点的最小旅行时；然后，利用同样的方法，从接收点开始，反向追踪炮点到接收点的射线路径、该方法能适于各种复杂的非均匀介质，极大地提高了射线追踪的精度。 相似文献
2.
3.
地震波传播的微分几何描述 总被引：6，自引：2，他引：6
本文以复杂介质中的走向函数为基础,引入走时场的微分流形,从微分几何观点看,复杂介质中的射线方程即为流形上的测地线方程,复杂介质中的标量波动方程即为流形上的协变标量场方程,从而建立起复杂介质中地震波传播的微分几何描述,文中还讨论该方法在射线追踪,波场变换等方面的应用。 相似文献
4.
以Moser方法为代表的最短路径射线追踪算法可以快速稳定地获得整个追踪区域的全局最小走时和路径，但它存在两个缺陷：一是射线大多由折线呈锯齿状相连，长度和位置偏离真实射线路径；二是在低变速区容易出现射线路径多值现象．本文提出的界面二次源波前扩展法全局最小走时射线追踪技术（以下简称界面源法）旨在解决上述两个问题．不同于Moser方法，界面源法只在物性分界面上设置子波源点，子波出射射线可以到达任何不穿越物性界面而直接到达的空间点和界面离散点，在均匀块体内或层内地震波以精确的射线路径传播．显然，界面源法的子波出射方向数远远大于传统方法，算法的追踪误差主要由界面离散引起的，因此，界面源法很好地解决了Moser法存在的问题，大大提高了追踪的精度．同时，由于界面源法的子波源点数远远小于Moser法，因而效率也很高．模型实算证实了该算法的高效性． 相似文献
5.
Simulating the propagation of seismic energy in an inhomogeneous, possibly anisotropic earth model is a challenging problem, both theoretically and from the computer implementation side. While most publications in the geophysical literature focus on the theory, this paper illustrates how recent advances in computer science can be of great value to researchers who try to create numeric realizations of complex geophysical processes such as seismic raytracing. Modern concepts including object- oriented programming, combined with new and often free-of-charge development tools can significantly improve productivity and flexibility in implementing and exploring new ideas.A variety of issues important in wavefront-construction based raytracing are ad- dressed and it is demonstrated how modern program design may help researchers to more quickly produce better software implementations and free up time to focus on the fundamental theoretical problems. 相似文献
6.
Abstract In dealing with the transient sediment transport problem, the commonly used uncoupled model may not be suitable. The uncoupling technique is intended to separate the physical coupling phenomenon of water flow and sediment transport into two independent processes. Very often, as a result, severe numerical oscillation and solution instability problems appear in the simulation of transient sediment transport in alluvial channels. The coupled model, which simultaneously solves water flow continuity, momentum and sediment continuity equations, gives fewer numerical oscillation and solution instability problems. In this article, a coupled model using a matrix double-sweep method to solve the system of nonlinear algebraic equations has been developed. Several test runs designed on the basis of a schematic model have been performed. The numerical oscillation and solution instability problems have been investigated through a comparison with those obtained from an uncoupled model. Based on the proposed case studies, it can be concluded that, for transient bed evolution, the performance of the coupled model is much better than that of the uncoupled model. The numerical oscillation is reduced and the solution is more stable. This newly developed coupled model was also applied to the Cho-Shui River in Taiwan. This application study implied that the effect of the peaky flood wave propagation on the bed evolution could be simulated better by the coupled model than by the uncoupled model. 相似文献
7.
《中国地震研究》2021,(1)
Accurate calculations of travel times and raypaths of reflection waves are important for reflection travel time tomography. The multistage shortest path method(MSPM) and multistage fast marching method(MFMM) have been widely used in reflection wave raytracing, and both of them are characterized by high efficiency and accuracy. However, the MSPM does not strictly follow Snell's law at the interface because it treats the interface point as a sub-source, resulting in a decrease in accuracy. The MFMM achieves high accuracy by solving the Eikonal equation in local triangular mesh. However, the implementation process is complex. Here we propose a new method which uses linear interpolation to compute the incident travel time of interface points and then using Snell's law to compute the reflection travel time of grid points just above the interface. Our new method is much simpler than the MFMM; furthermore, numerical simulations show that the accuracy of the MFMM and our new method are basically the same, thus the reflection tomography algorithms which use our new method are easier to implement without decreasing accuracy. Besides, our new method can be extended easily to other grid-based raytracing methods. 相似文献
8.
In this paper, the dynamic stress concentration and scattering of SH-waves by bi-naterial structures that possess an interface elliptic cavity are investigated. First, by using the complex function method, the Green‘s function is constructed. This yields the solution of the displacement field for an elastic half space with a semi-elliptic canyon impacted by an anti-plane harmonic line source loading on the horizontal surface. Then, the problem is divided into an upper and lower half space along the horizontal interface, regarded as a harmony model. In order to satisfythe integral continuity condition,the unknown anti-plane forces are applied to the interface. The integral equations with unknown forces can be established through the continuity condition, and after transformation, the algebraic equations are solved numerically. Finally, the distribution of the dynamic stress concentration factor (DSCF) around the elliptic cavity is given and the effect of different parameters on DSCF is discussed. 相似文献
9.
Dipping anisotropic clastic strata are ubiquitous in fold and thrust belts. Geological structures below these strata will be mispositioned laterally and vertically on seismic images if we do not properly correct for seismic anisotropy during migration. The magnitude of this lateral mispositioning of a target structure varies with source‐receiver offset, so reflection points will be smeared in the final stacked image. Raytracing demonstrates the lateral‐position and smear phenomena when imaging structures below tilted transversely isotropic media. Analysis of the raytracing results predicts the quantity of lateral‐position error and reflection‐point smear on a seismic image. We created numerical‐model seismic data to show reflection‐point smear on synthetic seismic images and to evaluate the accuracy of the predictions from raytracing. 相似文献
10.
11.
SH波入射时多个半圆形沉积谷地附近浅埋圆形孔洞的动力分析 总被引：3，自引：0，他引：3
本文采用辅助函数的思想,利用复变函数和多级坐标的方法给出了SH波入射条件下多个半圆形沉积谷地附近浅埋圆形孔洞动力分析问题的解答。将整个求解区域分割成两部分来处理,区域I为多个半圆形沉积谷地,区域II为浅埋圆形孔洞附近带半圆形凹陷的半无限弹性空间。在区域I和II中分别构造位移解,并在二个区域的“公共边界”上实施位移应力的连续条件,建立求解该问题的无穷代数方程组。最后,本文给出了算例和数值结果,并对其进行了讨论。 相似文献
12.
J.H. Leven 《Physics of the Earth and Planetary Interiors》1985,38(1):9-27
During the period 1974 to 1977, a long range seismic refraction project was conducted in Central Australia, along a profile extending south from Darwin. Earthquakes from the Banda Sea region were used as seismic sources for this experiment. An analysis by Hales and co-workers of the resulting data based on travel times, and using geometric raytracing techniques, has resulted in the construction of an upper mantle velocity model. Using synthetic seismograms to model amplitudes, it is shown that additional constraints can be placed on the derived velocity profile. The low velocity zone beneath the “200 km” discontinuity is found to have a more abrupt onset than was previously suggested. A smaller discontinuity at 325 km depth is now implied. The analysis suggests that the “400 km” discontinuity is a first order velocity increase, whereas all other observed upper mantle discontinuities are more satisfactorily modelled as second order type structures. 相似文献
13.
最短路径射线追踪算法,用预先设置的网络节点的连线表示地震波传播路径,当网络节点稀疏时,获得的射线路径呈之字形,计算的走时比实际走时系统偏大. 本文在波前扩展和反向确定射线路径的过程中,在每个矩形单元内,通过对某边界上的已知走时节点的走时进行线性插值,并利用Fermat原理即时求出从该边界到达其他边界节点的最小走时及其子震源位置和射线路径,发展了相应的动态网络算法. 从而克服了最短路径射线追踪算法的缺陷,大大提高了最小走时和射线路径的计算精度. 相似文献
14.
Fractal velocity models in refraction seismology 总被引：4，自引：0，他引：4
The introduction of noise components to a simple crustal velocity model is shown to markedly affect the appearance of synthetic seismograms calculated according to ray theory applied to refraction experiments. Here we simulate noise by a self-similar process with a power spectral density which falls off as inverse wavenumber to a simple power (0–2). The major effect is to destroy the coherency of the arrival branches normally expected from deterministic velocity models; the arrival amplitudes also show large trace-to-trace variations and considerable sensitivity to shot position. Some of these differences can be ascribed to the variety of noise model chosen (i.e., white noise, flicker noise and brown noise). It is argued that there is no clear distinction between coherent noise and geological structure. 相似文献
15.
In this study, the null-field boundary integral equation method (BIEM) and the image method are used to solve the SH wave scattering problem containing semi-circular canyons and circular tunnels. To fully utilize the analytical property of circular geometry, the polar coordinates are used to expand the closed-form fundamental solution to the degenerate kernel, and the Fourier series is also introduced to represent the boundary density. By collocating boundary points to match boundary condition on the boundary, a linear algebraic system is constructed. The unknown coefficients in the algebraic system can be easily determined. In this way, a semi-analytical approach is developed. Following the experience of near-trapped modes in water wave problems of the full plane, the focusing phenomenon and near-trapped modes for the SH wave problem of the half-plane are solved, since the two problems obey the same mathematical model. In this study, it is found that the SH wave problem containing two semi-circular canyons and a circular tunnel has the near-trapped mode and the focusing phenomenon for a special incident angle and wavenumber. In this situation, the amplification factor for the amplitude of displacement is over 300. 相似文献
16.
浅埋圆柱形弹性夹杂附近等腰三角形凸起地形引起的SH波的散射 总被引：4，自引：0，他引：4
文章利用“契合”思想,给出了地下弹性夹杂与地面上的等腰三角形凸起地形引起的SH波散射问题的解析解答。利用复平面下坐标移动,通过区域Ⅰ和区域Ⅱ以及区域Ⅱ和区域Ⅲ的“公共边界”位移应力连续条件,建立起求解该问题的无穷代数方程组并截断有限项进行求解,最后通过具体算例及结果分析得出相应结论。 相似文献
17.
18.
This paper deals with an analysis method for the response and motion of soil-like rigid-plastic bodies under seismic loading conditions. A continuity condition to determine the acceleration distribution within the rigid-plastic body when the failure occurs during seismic motions is proposed. Combining this continuity condition of acceleration and the ‘Generalized Limit Equilibrium Method (GLEM)’, the responses of the earth structure during seismic motions as well as the permanent displacements can be obtained, where GLEM is one of the limit equilibrium methods proposed by the authors for static problems and providing the approximate solution for Kötter's equation. The theoretical formulation of the method, the illustrative examples, and some comparisons between the analytical and experimental results are demonstrated. 相似文献
19.
Iterative solvers preconditioned with algebraic multigrid have been devised as an optimal technology to speed up the response of large sparse linear systems. In this work, this technique was implemented in the framework of the dual delineation approach. This involves a single groundwater flow linear solution and a pure advective transport solution with different right-hand sides. The new solver was compared with other preconditioned iterative methods, the MODFLOW's GMG solver, and direct sparse solvers. Test problems include two- and three-dimensional benchmarks spanning homogeneous and highly heterogeneous and anisotropic formations. For the groundwater flow problems, using the algebraic multigrid preconditioning speeds up the numerical solution by one to two orders of magnitude. The algebraic multigrid preconditioner efficiency was preserved for the three dimensional heterogeneous and anisotropic problem unlike for the MODFLOW's GMG solver. Contrarily, a sparse direct solver was the most efficient for the pure advective transport processes such as the forward travel time simulations. Hence, the best sparse solver for the more general advection-dispersion transport equation is likely to be Péclet number dependent. When equipped with the best solvers, processing multimillion grid blocks by the dual delineation approach is a matter of seconds. This paves the way for its routine application to large geological models. The paper gives practical hints on the strategies and conditions under which algebraic multigrid preconditioning would remain competitive for the class of nonlinear and/or transient problems. 相似文献
20.
一个以遗传算法为基础的结构可靠性分析方法 总被引：5，自引：0，他引：5
本文首先总结了一次二阶矩法在结构可靠性分析时的五个弱点。针对其在寻求验算点时需要求极限状态函数的导数以及在处理多峰性极限状态函数时存在的不收敛或收敛于局部验算点等弱点，探讨了将遗传算法应用于结构可靠性分析的可能性，并在分析其应用时存在的具体问题的基础上提出了一个以智能生物为基础的遗传算法。本文的分析计算表明：遗传算法在结构可靠性分析中是适用的，它可以克服一次二阶矩法在求解验算点时存在的几个弱点。本 相似文献