共查询到19条相似文献，搜索用时 768 毫秒
1.
作为一种模拟地震波走时的传统途径，射线追踪不仅对地震波理论研究具有重要意义。而且也可直接应用于地震波反演及偏移成像等过程。在一般条件下，射线追踪只能用数值方法来实现，但当介质上具有垂向不均匀性时，射线路径及其走时可以由解析式给出。根据适当的解析公式，可建立反射点和射线路径以及走时之间的关系，并用VC^ 6．0编制了基于图形界面的计算程序。这些工作对于逐步开发地震波的偏移和反偏移程序具有基础性的意义。 相似文献
2.
3.
二维介质速度成像的新方法是基于地震波的旅行时间，由射线路径沿线的速度分布特征所决定的，将二维介质中的速度分布用两个分别以二维空间坐标为自变量的勒让德多项式的线性组合来描述。根据走时积分方程，即可建立地震波走时与其射线路径沿线的速度分布之间的对应关系，进而重建出地震波射线覆盖区域的速度结构。该成像技术是一个迭代过程，且适用于折射波及反射波。 相似文献
4.
地震波走时是克希霍夫叠前深度偏移的重要参数。三维克希霍夫偏移成像由于需频繁读取大数据量的走时表文件,所以计算效率不高。这里提出一种通过压缩三维射线追踪走时表,来提高克希霍夫偏移计算效率的方法：首先规划一个具有规则网格控制点并覆盖所有走时表点集的最小长方体区域,然后以三维三次B样条函数为插值基函数进行最小二乘法曲面体拟合,求出规则网格控制点的数值并以数组形式存储入内存,采用稀疏化存储进一步节省了内存空间。在偏移成像时,再由这些规则网格控制点的数值,使用线性插值公式解编出走时表。实际资料算例验证了该走时表压缩方法不仅近似精度高,计算稳定度高,计算效率高,而且由于省去了频繁进行大数据量走时表文件的读写操作,所以克希霍夫偏移的计算效率提高了二倍以上。 相似文献
5.
这里将单程波真振幅方程与分步傅里叶算子（SSF）相结合,同时还结合了保幅算法和分步傅里叶算法的优点,因此该方法具有计算量小,占内存少,能处理横向变化的速度等优点。并且克服了傅里叶有限差分方法偏移后的振幅都有很大的偏差的不足。与目前广泛应用的常规的分步有限差分叠前深度偏移相比,具有成像精度高,保持地震波动力学特征等优点。在Marmousi模型上成功地进行了真振幅分步傅里叠前深度偏移处理,取得了理想的成像效果。 相似文献
6.
在隧道或巷道工程地震超前探测中较常用偏移成像技术,计算地震射线走时是该技术的核心部分。由于在近似地下全空间区域内成像,隧道或巷道、空洞、采空区等非穿透体对地震波走时计算有较大影响。为此,文中发展了基于FMM(fast marching method)的含非穿透体的走时算法。该算法采用非穿透体区域标记法,当FMM窄带区在计算到非穿透体时会自动避开或绕过,使得波前推进更加符合实际传播情况。这种算法也适用于起伏地表模型,只需要将起伏地表以上区域也作为非穿透体来对待。因而,新算法可以同时处理含有起伏地表的模型。改进的算法与常规算法相比只是增加了标记点,保持了FMM的计算精度和效率。理论模型试验表明,改进的算法能够较准确地计算走时,对复杂异常体的适应性较强,而且有很好的稳定性。 相似文献
7.
8.
先后推导了水平层与倾斜层2-D横向各向同性介质地震波反射波走时曲线方程,并就似SH波进行了模型计算,发现2-D情况下,似SH波的反射波走时曲线基本上符合双曲线形式。 相似文献
9.
为了实现适应崎岖海底、大陡坡、海底火山等复杂海底地质条件且灵活、稳定、精度高的射线追踪方法,并基于该方法详细分析复杂海底条件下各种地震波型的运动学特征,综合多种算法的优势,实现了一种快速推进迎风插值射线追踪方法。首先,采用混合网格法剖分复杂海底地质模型;其次,通过融入迎风差分思想的线性插值策略来构建精度高、无条件稳定且灵活的局部走时和射线路径计算公式;然后,综合应用这些公式和多级次快速推进法,灵活计算整个模型各种类型地震波的走时,并基于逆向追踪方法计算射线路径;最后,对该方法进行了精度分析,发现其能够获得相对高的走时和射线路径计算精度,且反射波的计算精度远高于入射波。此外,计算实例分析显示,初至波中富含折射波和陡倾构造的反射波在很大偏移距还能被接收,崎岖海底各种波型传播路径复杂;基于此提出加大采集排列长度和采用直达波走时可提高复杂构造成像质量等复杂海底地震数据采集与处理方面的思考与建议。 相似文献
10.
在利用有限差分等基于网格的数值分析方法解地震波走时所满足的程函方程时,由于速度模型的网格化离散等原因,会使走时在各网格节点之间不具有计算射线路径所要求的光滑性,即走时在邻近网格节点之间不具有连续的一阶导数。因此,直接利用网格节点走时计算射线路径会使最终的射线路径不光滑。为解决这个问题,已有研究者提出了基于B样条插值的逆向梯度方案（法）。然而,在速度发生突变时,B样条逆向梯度法所计算出的射线路径会具有较大的误差。针对这个问题,首先采用适合于解最小零偏差逼近及最佳平方逼近问题的Chebyshev多项式取代B样条对来自于分区多级计算方案的网格节点走时进行最佳逼近,得到在最小平方意义下的最优走时公式;然后采用与B样条逆向梯度法类似的计算过程得到光滑的射线路径。数值实验表明,利用Chebyshev多项式逼近走时可以得到具有很高精度的多次反射射线路径,在多次波偏移成像研究中具有潜在的价值。 相似文献
11.
旅行时线性插值（LTI）射线追踪算法是基于线性假设的,在向前处理过程中仅用按行或按列扫描的方法来计算节点旅行时没有考虑逆向传播射线,导致其计算精度与网格剖分大小有关,在处理复杂介质时会使得追踪出来的射线路径不一定满足最短旅行时。因此,笔者提出了两点改进措施：在向前处理时需采用全方位循环的方法来计算节点最小旅行时;在网格边界加入次生节点。模型试算结果表明：采用全方位循环的LTI法考虑了来自各个方向的射线,可提高其对复杂模型的适应能力;在节点间距相同的情况下,网格边界插入次生节点的LTI法较传统的LTI法计算精度至少可以提高一个数量级,同时,计算速度也更快;随着节点间距剖分的越精细,计算耗时下降也越明显,计算速度较传统的方法可提高n~10n倍。 相似文献
12.
基于连续地球模型的ＢＧ理论，系统地研究了地球物理反问题，描述了地球物理反问题高度的非唯一性，给出了构造和评价解估计的一系列准则和方法。根据ＢＧ方法的一般原理，针对地震勘探等一类具体问题，利用地震走时，实现了地震波速度重建。 相似文献
13.
We present a probabilistic analysis of seismic travel-time equations using the Bayesian Method. The assessment of models and data is crucial in 3D seismic travel-time tomography, and a method quantitatively assess the quality of both the data and the model is necessary in order to attain the most realistic results. The Bayesian method that we propose here is more effective than the frequentist approach, both in analysis time and uncertainty minimization, when processing large sets of tomographic data. 相似文献
14.
代数重建技术ART(Algebra Reconstruction Technology)是一种层析处理数学方法,具有计算速度快、图形简洁直观等优点,被应用于电磁波CT方法的数据处理中。但该处理方法易发散,结果难以控制,因此它的敛散性成为提高其应用的主要问题。这里通过在ART算法中加入阻尼系数的做法保证了迭代的收敛,并对阻尼系数的取值范围进行了探讨,成功地用于电磁波数据处理中。 相似文献
15.
16.
Yin Changchun 《吉林大学学报(地球科学版)》1994,(4)
研究发现，两种目前比较常用的可控源音频磁大地电流法（ＣＳＡＭＴ）的正演计算方法，即简化积分法和二次插值法之间计算结果差异较大。为了正确地评价这两种方法的计算精度，本文提出了利用直接积分进行，ＣＳＡＭＴ一维正演的算法，即通过将任意角度频率测深￣［１］正演结果沿发射源方向直接进行数值积分，得到有限长度发射源的ＣＳＡＭＴ正演结果。大量计算和精度评价发现，在这些算法中，简化积分法和本文提出的直接积分法计算结果吻合很好，是精度较高的算法，而二次插值法与前两种算法结果差别较大，是精度较低的算法。 相似文献
17.
TEM中心回线法计算考虑关断时间的全区视电阻率 总被引：2，自引：0，他引：2
详细推导了TEM法中心回线装置发射电流为斜阶跃波形条件下全区视电阻率迭代反演计算方法。对正演过程中对关断时间的积分采用了解析表示。对反演过程中求导运算,通过变量变换,将对电阻率的求导变换为对时间的求导,从而也给出简单的解析关系。避开了数值积分运算和差分替代导数的近似,提高了计算精度和速度。为克服迭代反演过程中要多次计算概率积分的耗时问题,采用标准数据库的方法一劳永逸地加以解决,也大大提高了计算精度和速度。对实测资料进行了计算,并通过算例对一些问题进行了讨论。 相似文献
18.
考虑冻土的非线性断裂力学特征,基于现有J积分测试方法,提出了一种新的冻土非线性断裂韧度的测试方法--修正因子法。在传统方法中J积分可以写成如下表达：J =Je +Jp,其中Je为J积分的线弹性分量,Jp为J积分的塑性分量。Je的计算相对简单;但计算Jp相对困难,原因是其中的参数Up难于确定。为了使计算简单,可将塑性分量改写为Jp =qJe,J积分则改写为J=(1+q)Je,q称为塑性修正因子。该方法通过循环加、卸载试验得到非线性载荷-位移曲线（P-Δ曲线）,并由载荷与位移的增量确定修正因子,在获得线弹性分量Je基础上,获得冻土非线性断裂韧度。采用该方法进行了冻土非线性断裂韧度测试,并将测试结果与相关文献结果进行比较,证明了该方法的合理性和有效性。 相似文献
19.
体积分方程法模拟复电阻率三维体电磁响应 总被引：1，自引：0，他引：1
利用体积分方程法计算了均匀半空间中复电阻率(激电)效应和电磁效应同时存在时的三维体响应。在计算中,对于需要计算三重积分的张量格林函数电荷项一次部分,应用一种差分近似的方法求解,这种方法在保证计算精度的同时更加便于计算机实现;采用二次剖分的算法解决了计算张量格林函数时的奇异值问题;计算含有贝塞尔函数的积分项时,利用一种结合连分式展开的高斯求积代替常规的快速汉克尔变换方法。验证了计算结果并分析了三维体复电阻率模型(Cole-Cole)参数对正演结果的影响,为三维体复电阻率及其参数反演提供了正演依据。 相似文献