共查询到16条相似文献，搜索用时 247 毫秒
1.
2.
给出了计算多层垂直对称轴横向各向同性（VTI）介质精确射线路径和走时的方法,所用的体波相速度公式、群速度公式和Snell定律都是严格的显式解析公式. 任意基本波的射线路径和走时计算问题都可以转化成一个等效的透射问题,再用文中的公式来计算,具体实现方法用一个多次波和一个首波的实例给出. 最后分别用精确公式和Thomsen近似公式计算了相同模型相同基本波的走时曲线. 比较两者计算结果可发现, 近似公式反复使用会使误差积累,同时揭示了近似公式适用范围的局限性,强调了使用近似公式需要注意其适用范围的重要性. 相似文献
3.
4.
5.
最短路径射线追踪算法,用预先设置的网络节点的连线表示地震波传播路径,当网络节点稀疏时,获得的射线路径呈之字形,计算的走时比实际走时系统偏大. 本文在波前扩展和反向确定射线路径的过程中,在每个矩形单元内,通过对某边界上的已知走时节点的走时进行线性插值,并利用Fermat原理即时求出从该边界到达其他边界节点的最小走时及其子震源位置和射线路径,发展了相应的动态网络算法. 从而克服了最短路径射线追踪算法的缺陷,大大提高了最小走时和射线路径的计算精度. 相似文献
6.
本文基于宾川气枪地震信号发射台激发的地震波信号,利用线性叠加和时频域相位加权叠加方法提高信噪比,通过反褶积、插值拟合和波形互相关等方法,获得2021年漾濞M;6.4地震前后气枪震源初至波信号走时变化特征.结果表明两种叠加方法得到的走时变化趋势基本一致,在沿发震断裂带附近的3个台站观测到的初至波走时延迟为7.3~14.4 ms,在射线路径上平均波速降低了0.08%~0.12%,距离震中位置较远的台站走时延迟为2~6 ms.分析认为,观测到的走时变化主要是漾濞M;6.4地震引起强地面运动造成浅层介质疏松、同震破裂导致震源区地下介质裂隙增加和地下流体侵入等共同作用造成,相关研究成果有助于加深和促进对此次地震震源物理过程的认识和了解. 相似文献
7.
8.
《地球物理学进展》2016,(1)
CT技术已经广泛应用到科学和工程领域,特别是在地球物理学和土木工程质量无损检测方面.本文用精密可控的地震波雷达在芦山震中区持续300 s向地壳主动发射低能量的线性调频弹性波,通过高灵敏度的地震仪器接收信号,实现对上地壳剖面波速结构进行无破坏性探测.用短时傅里叶变换分析记录到的线性调频信号的时间和频率特征,并用能量累积法提取信号的走时形成走时曲线图,识别震相后用射线追踪方法进行P波走时反演,得到了芦山地震震中区上地壳剖面二维波速结构.分析结果表明地震波雷达发出的信号符合控制要求,重复性高达99.7%;用能量累积法得到的地震波走时曲线易于识别震相.从芦山地震震中区上地壳剖面二维波速结构中可以看出,在芦山震中区附近有明显的波速差异,并且出现了向上凸起、向下急剧下降,再缓慢上升,呈"铲形"特征,存在低速层,分析出芦山地震是一次逆冲型地震,并且距震中区近的检测点波速比偏低,距震中区远的检测点波速比偏高.通过与其他方法研究结果对比,证明了用地震波雷达信号走时反演地壳波速结构是可行的. 相似文献
9.
基于图形结构的三维射线追踪方法 总被引：39，自引：16，他引：23
在地震层析成像研究中,为了克服最小走时射线路径追踪方法存在的问题,对该方法计算过程中的关键步骤进行了改进．在节点走时的计算中引入Ｂｒｅｓｅｎｈａｍ画线算法;在最小走时节点查寻中,结合使用快速排序算法与插入排序算法,替代以往方法中多采用的堆排序算法;所采用的节点设置方式,可以引入速度界面,还可以实现反射波射线追踪．模型计算证明,改进的最小走时射线路径方法具有精度高,速度快的特点,所提出的三维空间反射波射线追踪算法简便易行。 相似文献
10.
本文叙述无水注入情况下,中心受压的大理石板(60060020毫米)破裂前,P波走时变化的实验结果。 路径穿过未来裂缝、靠近受压中心时,走吋上升一回降;不靠近受压中心则走时只上升。路径近于平行且靠近未来裂缝,走时一般不变。路径不穿过未来裂缝和受压中心,发射点附近形变小时,走时不变;发射点附近形变大时,走时一般下降。路径穿过受压中心吋,有些走时先上升后恢复再下降。有些则下降后回升。 这说明在岩石破裂前,破裂孕育区附近同时存在P波速正负异常区和正常区。 同时还观测到,P波速度负异常区随中心压力增加逐步扩大。 由裂缝一侧波源所得的正常区正好是由另一侧波源所得的负异常区。 相似文献
11.
We present a new method of three-dimensional (3-D) seismic ray tracing, based on an improvement to the linear traveltime interpolation (LTI) ray tracing algorithm. This new technique involves two separate steps. The first involves a forward calculation based on the LTI method and the dynamic successive partitioning scheme, which is applied to calculate traveltimes on cell boundaries and assumes a wavefront that expands from the source to all grid nodes in the computational domain. We locate several dynamic successive partition points on a cell's surface, the traveltimes of which can be calculated by linear interpolation between the vertices of the cell's boundary. The second is a backward step that uses Fermat's principle and the fact that the ray path is always perpendicular to the wavefront and follows the negative traveltime gradient. In this process, the first-arriving ray path can be traced from the receiver to the source along the negative traveltime gradient, which can be calculated by reconstructing the continuous traveltime field with cubic B-spline interpolation. This new 3-D ray tracing method is compared with the LTI method and the shortest path method (SPM) through a number of numerical experiments. These comparisons show obvious improvements to computed traveltimes and ray paths, both in precision and computational efficiency. 相似文献
12.
13.
在地震定位中常常需要求解震源轨迹,但由于复杂介质中的震源轨迹较为复杂,难以给出其解析解,因此震源轨迹的计算通常仅限于简单介质模型.本文基于最小走时树射线追踪技术,提出了一种计算复杂介质中震源轨迹的方法.为回避发震时间问题,以观测到时差作为震源轨迹的约束条件.首先从模型节点中选出少量理论到时差与观测到时差之绝对差,即双重时差较小的点作为震源轨迹的代表点,然后以其中双重时差最小的点为初始点,在双重时差场中利用最小走时树射线追踪方法计算出初始点到其他震源轨迹代表点的射线路径作为震源轨迹.当选的震源轨迹代表点较多时,得到的震源轨迹较为粗略,此时可去掉射线经过次数较少的代表点的射线路径使震源轨迹更为精细.为减少计算量,对最小走时树射线追踪方法的终止条件做了修正.以一个复杂介质模型中的地震为例,计算了包括速度扰动、到时扰动等不同情况下的震源轨迹,结果表明所提出的震源轨迹计算方法切实可行. 相似文献
14.
15.
走时计算已广泛用于地震建模、成像及速度分析等诸方面．基于地震波场的Ｈａｍｉｌｔｏｎ力学性质,本文探索应用适于Ｈａｍｉｌｔｏｎ力学的计算方法──辛几何算法对地震射线 进行走时及路径的计算．利用辛几何算法和属于耗散算法的四阶 Ｒｕｎｇｅ－Ｋｕｔｔａ算法进行射线 路径和走时的计算对比,结果表明这两类算法数值精度相当,但辛几何算法的速度却快了３ 倍．本文还利用二阶Ｅｕｌｅｒ型辛差分格式对Ｍａｒｍｏｕｓｉ模型进行了射线追踪计算,结果显示所 得射线具有良好的光滑性和较强的阴影区穿透能力． 相似文献
16.
We present a new ray bending approach, referred to as the Eigenray method, for solving two‐point boundary‐value kinematic and dynamic ray tracing problems in 3D smooth heterogeneous general anisotropic elastic media. The proposed Eigenray method is aimed to provide reliable stationary ray path solutions and their dynamic characteristics, in cases where conventional initial‐value ray shooting methods, followed by numerical convergence techniques, become challenging. The kinematic ray bending solution corresponds to the vanishing first traveltime variation, leading to a stationary path between two fixed endpoints (Fermat's principle), and is governed by the nonlinear second‐order Euler–Lagrange equation. The solution is based on a finite‐element approach, applying the weak formulation that reduces the Euler–Lagrange second‐order ordinary differential equation to the first‐order weighted‐residual nonlinear algebraic equation set. For the kinematic finite‐element problem, the degrees of freedom are discretized nodal locations and directions along the ray trajectory, where the values between the nodes are accurately and naturally defined with the Hermite polynomial interpolation. The target function to be minimized includes two essential penalty (constraint) terms, related to the distribution of the nodes along the path and to the normalization of the ray direction. We distinguish between two target functions triggered by the two possible types of stationary rays: a minimum traveltime and a saddle‐point solution (due to caustics). The minimization process involves the computation of the global (all‐node) traveltime gradient vector and the traveltime Hessian matrix. The traveltime Hessian is used for the minimization process, analysing the type of the stationary ray, and for computing the geometric spreading of the entire resolved stationary ray path. The latter, however, is not a replacement for the dynamic ray tracing solution, since it does not deliver the geometric spreading for intermediate points along the ray, nor the analysis of caustics. Finally, we demonstrate the efficiency and accuracy of the proposed method along three canonical examples. 相似文献