共查询到19条相似文献，搜索用时 125 毫秒
1.
《物探化探计算技术》2021,(2)
保幅高斯束逆时偏移是一种兼具高斯束偏移灵活、高效和逆时偏移高精度特性的深度域成像算法,不仅能够消除偏移过程中由于入射角变化和几何扩散引起的振幅误差,而且具有面向目标成像的能力。其以Kirchhoff积分为基础,采用高斯束的加权积分来构建格林函数,进而实现正反向波场的延拓,最后通过相应的成像条件进行成像。这里从基于相速度的各向异性射线追踪算法出发,采用互相关成像条件,将保幅高斯束逆时偏移推广到VTI介质。随后通过模型试算对方法的正确性和有效性进行了验证。 相似文献
2.
栗学磊 《吉林大学学报(地球科学版)》2010,(Z1)
由于具有计算速度快的优点,射线追踪方法在地震勘探中一直起着非常重要的作用。V.Cerveny在上世纪80年代初提出了高斯波束方法,该方法同时考虑了波的动力学和运动学特征,并且无需两点试射追踪,因而运算速度快且精度高;另外高斯波束对异常区波场有较好的效果,比如焦散点、阴影区、临界和超临界反射。研究的主要内容是2.5维简单结构模型的高斯波束法。利用2.5维模型进行数值计算相当于在二维空间内计算三维问题,可以节约很大一部分计算。 相似文献
3.
基于波前构建的射线追踪:一种Java实现 总被引：3，自引：2，他引：1
针对波前构建法射线追踪在实现过程中存在的困难和问题,提出了一种基于Java语言的实现方案。与常规的Fortran实现方案相比,Java实现避免了由射线插值和网格变换所带来的额外数据重排及数据点反复标记问题。因此,在计算过程中可以随时对数据点进行插入操作而不用考虑再定位问题。鉴于在关于波前构建法的文献中很少涉及算法和计算流程的设计问题,首先对波前构建法的基本方程、数值计算方案、计算流程及具体算法进行了简要的回顾和讨论,然后讨论了Java类和Java函数设计。我们的经验和数值计算结果表明,用Java语言实现的波前构建过程能够降低编程工作的难度、提高射线追踪的计算效率并且为今后的跨平台开发提供良好的基础。 相似文献
4.
基于二维三次卷积插值算法的波前构建射线追踪 总被引：1，自引：0，他引：1
针对波前构建法在实现波前射线路径追踪过程中所遇到的速度以及速度导数的插值问题,提出一种更精确的二维三次卷积插值算法,将该算法应用到波前构建法中,得到比较精确的射线路径。实例证明,所采用的插值算法在应用到波前构建的过程中,在求解非网格节点的速度及导数时与双三次多项式插值算法相比能够更加快速准确地获得射线路径。 相似文献
5.
《物探化探计算技术》2017,(3)
常规的逆时偏移算法都是基于双程波动方程的数值解法,实际上利用Kirchhoff积分进行波场延拓,也可以实现逆时偏移。以Kirchhoff积分为基础,采用高斯束叠加积分表示的格林函数来表示正反向延拓波场,将正反向延拓波场进行保幅的互相关得到地下成像。保幅高斯束逆时偏移方法结合了射线类偏移方法的高计算效率和逆时波动方程偏移的高精度,能很好地处理焦散点、陡倾角等问题,并且具有面向目标成像的能力,并消除了偏移过程中的几何扩散效应和入射角变化引起的振幅误差。 相似文献
6.
相对于常规射线层析速度建模,基于波动理论的层析速度建模考虑了波的带限特性,反演分辨率更高。波动理论层析的核心在于波路径(有限频核函数)的计算。文中详细介绍了一种基于高斯束算子计算有限频核函数的方法,分析了初始束宽度和高斯束出射角度间隔对计算精度的影响;并针对高斯束近源处误差较大的缺陷,提出了改进的束参数以提高近源精度;详细分析了初始束宽度和角度间隔对改进高斯束方法的影响及改进高斯束的聚焦特性;数值算例验证了该方法在缓变介质中计算有限频核函数的可行性,计算效率较高且可处理回折波的核函数。 相似文献
7.
青藏高原东北缘是青藏高原隆升的前缘地带，其构造变形特征一直是研究的热点。前人在该区域进行过面波成像研究，但不同结果存在较明显的差异，可能与不同成像方法的分辨率有关，此次研究分别使用基于射线的方法和基于程函方程的方法进行了瑞利波相速度成像。笔者对比了20 s、30 s和60 s的成像结果，两种方法的结果表现出相似的速度特征。总体的速度特征与研究区域的主要构造单元分布相吻合，青藏高原东北缘表现出低速异常，鄂尔多斯块体表现出高速异常。在银川河套地堑，基于程函方程的成像在20 s和30 s得到更明显的低速异常，低速异常体的分布与地堑的轮廓吻合更好，说明基于程函方程的方法对数据利用更充分。 相似文献
8.
9.
10.
基于弹性波动方程,推导了二维非均匀介质条件下纵横波解耦且无转换波的一阶速度一应力弹性波方程,使得弹性波在延拓过程中实现矢量偏移保证偏移质量,且不产生其他类型的转换干扰波。从激发时间成像条件以及反射系数理论出发,并与纵横波解耦且无转换波的一阶速度一应力弹性波方程结合,给出了适用于P＋s波震源模拟数据的成像条件以及偏移方法。文中给出了两个典型模型的理论记录的偏移结果,试验表明,笔者提出的方法能够对P＋s波震源的产生的波场准确归位,提高资料处理精度。 相似文献
11.
二维介质速度成像的新方法是基于地震波的旅行时间，由射线路径沿线的速度分布特征所决定的，将二维介质中的速度分布用两个分别以二维空间坐标为自变量的勒让德多项式的线性组合来描述。根据走时积分方程，即可建立地震波走时与其射线路径沿线的速度分布之间的对应关系，进而重建出地震波射线覆盖区域的速度结构。该成像技术是一个迭代过程，且适用于折射波及反射波。 相似文献
12.
井间地震层析成像分辨率研究 总被引：12，自引：1，他引：11
通过大量数值模型试验,并根据射线密度、射线正交性和点扩散函数,对井间地震层析成像的分辨率问题进行了研究。结果表明,井间地震层析成像的分辨率在空间上是可能性变的,垂直方向的分辨率好于水平方向,而且高速异常体和低速异常体的图像形态和分辨率是非一致的。高速异常体波形CT的图像水平方向最小分辨率距离为1/2波长,垂直分辨率最小分辨率距离为3/10波长;而低速异常体的波形CT图像的水平方向的最小分辨率约为2/5波长,垂向最小分辨率约1/5个波长;走时CT的图像水平方向最小分辨率为3个像元,垂直方向的最小分辨率距离为2个像元;当井深与井间距之比为2时,图像分辨率最高。 相似文献
13.
Analytical Study of Underground Explosion-Induced Ground Motion 总被引：1，自引：0，他引：1
A theoretical approach is presented to study the ground motion induced by an underground tunnel explosion. The ground motion is caused by two coupled stress waves, i.e., the reflected body wave and the secondary surface wave or Raleigh wave. Based on the principle of conservation of momentum at the wavefronts, the reflected body waves along the ground surface are derived. The interaction of the body wavefront and the ground surface induces the secondary surface wave which transfers outwards on the ground. The particle velocity and particle acceleration on the ground surface are subsequently derived. The analytical results are compared with results from numerical simulations and empirical formulae with different material damping ratios. The effects of the loading density and the material damping on the ground motion are investigated. Finally, the limitations of the proposed theoretical approach for ground motion prediction are discussed. 相似文献
14.
地震波走时计算的逆风差分算法 总被引：5，自引：0，他引：5
快速精确地旅行时计算在地震资料的叠前偏移与层析反演中起重要作用。利用逆足差分格式求解程函方程，首先在波阵面上寻找全局极小点，然后计算该点周围各点的地震波初至走时，同时，采用局部算法以确保计算的稳定性。该算法对于任意复杂地质模型都具有较高的精度。 相似文献
15.
Anomalous crustal and upper mantle structure of northern Juan de Fuca plate is revealed from wide-angle seismic and gravity modelling. A 2-D velocity model is produced for refraction line II of the 1980 Vancouver Island Seismic Project (VISP80). The refraction data were recorded on three ocean bottom seismometers (OBSs) deployed at the ends and middle of a 110 km line oriented parallel to the North American continental margin. The velocity model is constructed via ray tracing and conforms to first-arrival amplitude observations and travel time picks of direct, converted and reflected phases. Between sub-sediment depths of 3 to 11 km, depths normally associated with the lower crust and upper oceanic mantle, the final model shows that compressional-wave velocities decrease significantly from southeast to northwest along the profile. At sub-sediment depths of 11 km at the northwestern end of the profile, P-wave velocities are as low as 7.2 km/s. A complementary 2-D gravity model using the geometry of the velocity model and velocity–density relationships characteristic of oceanic crust is produced. The high densities required to match the gravity field indicate the presence of peridotites containing 25–30% serpentine by volume, rather than excess gabbroic crust, within the deep low velocity zone. Anomalous travel time delays and unusual reflection characteristics observed from proximal seismic refraction and reflection experiments suggest a broader zone of partially serpentinized peridotites coincident with the trace of a pseudofault. We propose that partial serpentinization of the upper mantle is a consequence of slow spreading at the tip of a propagating rift. 相似文献
16.
声波速度-频移衰减联合层析成像检测技术 总被引：1，自引：1，他引：0
混凝土构件中存在的诸如松散、空洞、缝隙夹层等缺陷,不仅造成声波速度的变化,而且对声波具有较强的吸收衰减作用,利用声波速度CT和衰减CT联合检测技术是检测混凝土构件工程质量的有效方法。介绍了声波速度CT检测原理和SIRT算法、接收波信号频移与介质衰减吸收系数的关系,导出了适用于CT成像反演算法的介质吸收因数沿射线路径的积分与频移参数的关系表达式。对存在松散和空洞缺陷的两种混凝土构件模型的声波检测结果,分别进行了速度CT和频移衰减CT处理,处理结果表明,声波速度和频移衰减联合层析成像检测技术能够获得更加准确可靠的检测效果。 相似文献
17.
A Ram 《Journal of Earth System Science》1978,87(11):263-270
A set of computer programs was written which enables one to determine automatically, the apparent velocity and azimuth of any portion of the seismic wave-train recorded at various seismic arrays. All that is required to use these programs is that the raw data be available on digital tape and that the user have access to any modern computer. The programs are written in Fortran IV and make use of the adaptive processing method whereby the arrival times of the wavefront at each sensor are accurately determined by cross-correlating the observed wavelet of interest with the corresponding wave on the beam trace. The new arrival times are then used to create a new and improved beam and the whole operation is repeated in an iterative manner until convergence takes place. 相似文献
18.
19.
CHANG Xu LIU YikeInstitute of Geology Geophysics Chinese Academy of Sciences All Datun Ro Chaoyang District Beijing E-mail: changxu@mail.c-geos.ac.cn 《《地质学报》英文版》2000,74(1):102-112
In this paper, 3-D velocity images of the crust and upper mantle beneath the Nanbei tectonic zone of China are constructed using P-wave travel time residuals of earthquakes, with the data supplied by China's seismic networks and the International Seismic Centre.During the model parameterization in the tomographic inversion the interpolation function of grid node velocities is used as the space function of velocity, and the velocity function is allowed to be discontinuous. The difficulties caused by large memory demand and high computing cost in solving the system of equations are avoided by utilizing the least squares QR decomposition algorithm. Thus, the stability of the algorithm is ensured. Though resolution images are not derived in the inversion process based on the resolution theory of Backus and Gilbert, a covariance resolution method is put forward by the authors. The resolution analyses have proved the reliability of the image results. 3-D ray tracing is conducted to obtain the ray paths in a s 相似文献