CN1292263C - 一种用于地震勘探中射线追踪的方法 - Google Patents

Abstract

本发明属于石油地震勘探资料处理方法技术领域，尤其涉及一种最小旅行时地震射线追踪的方法。本发明的矩形网格三点Fermat射线追踪法是基于矩形网格三点扰动法的一种提高计算速度的方法。取矩形网格三个点，在Fermat最小旅行时原则下求取中间点的位置，而不象扰动法那样依次扰动。因此，计算速度比扰动法提高2倍多，同时不受扰动量大小选择的困扰。该方法继承了矩形网格三点扰动法的优点，对任意离散的速度场，总能找到最短时间路径，避免了射线盲区和追踪路径并非时间最短路径问题。在地震偏移、反演和成析成像方法应用中具有广泛的应用前景。

Description

一种用于地震勘探中射线追踪的方法

技术领域：

本发明属于石油地震勘探资料处理方法技术领域，尤其涉及一种最小旅行时地震射线追踪的方法。

背景技术：

地震勘探法，即按一定规律和方式排布的各地震波检测点接收由按一定规律或方式排布和移动的各放炮点在地表激发所产生的地震波的折射波和/或反射波，通过对检测到的由各地质构造层界面处对该地震波产生的折射波和/或反射波返回地表的时间和波形进行处理和分析，以便能准确了解和确定地质构造的形态和位置，是目前在石油和天然气等勘探中指导勘探钻井孔位部署的一种常用和有效的方法。由于地下各地质构造层在组成及密度、分布的均匀程度等诸方面的不同，放炮点激发产生的地震波在其中的传播速度相应也各不相同，各检测点检测到的由各地质构造层界面处产生而返回地表的相应折射波和/或反射波的时间早晚及波形等也相应有所差别。

地震波在地下传播的过程中，波穿过的路径和介质的速度决定了波的旅行时间。根据波动理论，波从激发点(震源点，如地震时的震中)出发，以介质速度在地下传播，经过一段时间，到达接收点(记录站点)，所走的路径总是旅行时间最短。这种最短路径(射线旅行时最短的路径)就是波的射线路径。如何从激发点到接收点确定波传播的路径和旅行时，就是射线追踪研究的主要内容。在地震偏移方法、层析成像和反演方法中，广泛采用射线追踪方法确定地震波旅行时，因此射线追踪方法是一个重要研究领域。现有技术中存在几种射线追踪的方法，如两点射线追踪法(Feregra，1980)，有限差分方法解程函方程法(Vidale，1988)，波前重建估算旅行时和振幅(Vinje等，1993)，波前追踪法(Bube等，2001)。

但，所有射线追踪技术主要面临的两个问题是：

1、射线在追踪过程中产生的盲区(即射线路径走不到的区域)和并非时间最短路径。射线的盲区主要发生在地层的低速度带，当射线遇到低速层时，绕过低速层，从高速度层穿过，因而产生了盲区。

2、有些射线追踪方法有可能找不到最短路径(旅行时间最小)，如四方网格打靶法。运算速度也是实际应用中必须考虑的问题。耗时过多的算法，尽管能克服射线的盲区问题，也无法在实际应用中得到广泛应用。

为了克服上述存在的问题，现有技术中提出了一种矩形网格三点扰动(中间点做水平位置的滑动)射线追踪法，该方法是针对速度模型离散的情况提出的，不论什么样的速度场分布或低速带的存在，该方法总能找到最短时间的路径，避免了射线盲区和追踪路径并非时间最短路径问题。其具体操作方法是：

1、对追踪区域作矩形网格化，每个网格上分布速度值。在震源点到接收点用直线联结，得到初始射线路径，初始路径与网格的交点分别为p₁、p₂、…p_N。计算路径p₁p₂…p_N总旅行时。

2、取初始路径与矩形网格三个交点p_i-1，p_i和p_i+1，固定两个端点p_i-1和p_i+1，对中间点p_i修正，给出一个扰动量(扰动的增量)对其扰动，计算相应的旅行时间。所有可能是最小旅行时的位置都进行这样的位置滑动，计算位置滑动后的旅行时，然后比较旅行时，确定最小旅行时，从而确定具有最小旅行时的位置q_i，用q_i代替p_i的位置。

3、取q_i，p_i+1和p_i+2三个点，固定两个端点q_i和p_i+2，重复2的扰动和旅行时计算。确定q_i+1代替p_i+1的位置。这样依次取三点修正直到修正到p_N为止，并计算新路径p₁q₁q₂…p_N的旅行时T。

4、把p₁q₁q₂…p_N的连线当作新的初始射线路径，重复2，3步。

5、反复重复1-4步，得到一系列弯曲射线路径和相应的旅行时序列T₀＞＝T₁＞＝…，如果某个旅行时T_M满足T_M-1-T_M＜＝ε，则追踪停止。

该技术特点是，能避免射线盲区和射线追踪路径并非时间最短路径问题，但此方法存在两点技术问题：

1、射线追踪速度慢：在具体操作中，在作扰动的过程中，必须对所有可能的点依次进行扰动，这样就增加了运算量，大大降低了射线追踪的速度。

2、扰动量的选择问题：在实际的操作和应用中，若扰动量选择过小，精度一般可以达到要求，但耗时明显增加。反之，若扰动量选取过大，则达不到最小旅行时的要求，因此需要对扰动量进行反复多次试验，以确定有效和合理的取值范围。以达到一方面即可以满足计算精度的要求，另一方面又不增加无谓的运算量。在实际的应用中，扰动量的合理选取十分困难。

发明内容：

现有技术中的矩形网格三点扰动射线追踪法是针对速度模型离散的情况提出的，不论什么样的速度场、低速带的影响等，总能找到最短时间路径，避免了射线盲区和追踪路径并非时间最短路径问题，是一个较为理想的射线追踪方法。但该方法存在两大难题，其一是速度太慢。在作扰动的过程中，必须对所有可能的点依次进行扰动。然后计算旅行时，进行最小旅行时的比较，确定最终的路径的位置。其二是扰动量大小的取舍问题。扰动量过大，速度能加快，但找不到最短时间路径；扰动量过小，计算精度达到后，最小旅行时不会有任何实际意义的改变，但速度明显减慢。只有一个合适的扰动量才能即满足计算精度，又不增加无谓的运算量。最令人困扰的还有这个扰动量随计算深度的不同，需要进行调整。深度增加时，扰动量需要适当地减小。在实际计算过程中，实现这种要求相当困难。

本发明是基于矩形网格三点扰动射线追踪法，利用Fermat最小旅行时原理，直接取得扰动中间点的位置，不需要对中间点进行依次扰动，从而提高了矩形网格三点扰动射线追踪法计算速度，同时由于在Fermat最小旅行时原则下计算最小旅行时的位置，因此避免了扰动量选取问题。

具体的方法如下：

一种用于地震勘探中射线追踪的方法，其特征在于：本方法将射线追踪场进行矩阵网格化后，利用矩形网格三点和FERMAT原理，逐一对射线路径点进行追踪，形成最短走时射线追踪路径；

其具体方法为：

A、对射线追踪区域进行网格化操作；

B、建立深度域的二维速度模型的离散速度场；

C、对离散速度场做平滑处理：每个深度界面的网格，沿水平轴方向取奇数个点，如5，7，9等进行数学平均或中值滤波的平滑处理；

D、设定初始射线路径和该初始路径与所述网格的交点p₁、p₂、…p_N；；

E、计算初始射线路径的旅行时间T₀；

F、逐点追踪射线路径：

(1)取点：取初始射线路径上三个顺序交点(p₁，p₂，p₃，)；

(2)计算：固定两端点p₁和p₃，对中间点p₂，用Fermat原理计算得到最小旅行时对应的中间点q₁，并用新中间点q₁代替原有的中间点p₂，使p₁q₁p₃为p₁到p₃的时间最短的弯曲射线路径；

(3)再取q₁，p₃，p₄三点，对其进行修正，重复(2)的操作，这样依次取三点修正直到修正到射线路径的终点p_N为止；

(4)经过步骤(1)-(3)后，形成一条新的射线路径p₁q₁q₂…p_N，

(5)计算旅行时间：根据p₁q₁q₂…p_N路径计算新路径的旅行时间T₁；

(6)阈值比较：将前后两次旅行时间(T₀和T₁)的差值和阈值进行比较；

(7)重复步骤：若阈值比较的结果不符合精度要求，则以p₁q₁q₂…p_N作为初始路径，重复上述步骤F中的(1)-(6)的追踪过程，得出新旅行时间T₂，并与T₁进行阈值比较，不满足精度要求则反复重复步骤F中的(1)-(6)的追踪过程，不断形成新的射线路径，逐点逼近最短走时的射线路径；

G、输出射线路径步骤：

(1)若步骤F中的(6)阈值比较的结果符合精度要求，即旅行时间的差值符合精度要求：|T_M-1-T_M|＜＝ε，追踪停止；

(2)输出：相应的射线路径即为最小旅行时间的地震射线追踪路径，输出此射线路径。

在具体的应用中，所述步骤A的网格化操作中横向网格间隔大于纵向网格间隔。

在具体的应用中，所述步骤A的网格化操作中网格间隔为横向网格间隔与纵向网格间隔的比例为1∶1。

在具体的操作中，所述步骤B中的建立深度域的二维速度模型的离散速度场为：从地震资料中，用速度分析方法得到纵波和对应深度的横波的叠加速度。用狄克斯公式将叠加速度换算为层速度，并将时间-速度域转换到深度-速度域，形成深度域的速度模型。在震源点-检波点之间对速度场进行内插和外插就生成深度域射线追踪需要的二维速度模型的离散速度场。

在具体的操作中，所述步骤C中的建立深度域的二维速度模型的离散速度场：操作过程中作平滑处理，取5-9道的平滑。

在具体的应用中，所述步骤F中的FERMAT原理，其核心公式：

$t\left(r\right)=\frac{\sqrt{{(x-r)}^2+h_1^2}}{v_1\left(r\right)}+\frac{\sqrt{r^2+h_2^2}}{v_2\left(r\right)}$

$\frac{\mathrm{dt}\left(r\right)}{\mathrm{dr}}=\frac{r-x}{v_1\sqrt{{(x-r)}^2+h_1^2}}+\frac{r}{v_2\sqrt{r^2+h_2^2}}=0$

其中x：网格交点P_i+2-P_i水平间距；

r：网格交点P_i+2-P_i+1水平间距；

h₁：网格交点P_i+1-P_i垂直间距；

h₂：网格交点P_i+2-P_i+1垂直间距；

v₁，v₂分别为P_i+1点上下网格速度；

t：从P_i～P_i+2的旅行时；

dt/dr：时间对r的求导。

在实际的应用中，所述步骤F中的逐点追踪递推公式：

x＝p₃-p₁、q₁＝p₁+(x-r)。

x＝p_i+2-q_i、q_i-1＝q_i+(x-r)，(i＝1，2，……，N-1)

其中r为Fermat原理求解的具有最小旅行时的位置。

式中p₁、p₂，…p_N系列为初始网格交点的水平坐标；q₁是对p₂点修改后的水平坐标；然后对p系列依次进行修改，并利用前面修改后的q值，形成新的q系列，也就得到了p₁至p_N新的射线路径。

根据上述的技术方案，本申请的矩形网格三点Fermat射线追踪法是基于矩形网格三点扰动法的一种提高计算速度的方法。取矩形网格三个点，在Fermat最小旅行时原则下求取中间点的位置，而不象扰动法那样依次扰动。因此，计算速度比扰动法提高2倍多，同时不受扰动过程中的扰动量大小选择的困扰。

该方法继承了矩形网格三点扰动法的优点，对任意离散的速度场，总能找到最短时间路径，避免了射线盲区和追踪路径并非时间最短路径问题。在地震偏移、反演和成析成像方法应用中具有广泛的应用前景。

附图说明：

图1为本发明从S点到R点，穿过界面的射线旅行示意图；

图2为本发明矩形网格三点Fermat射线追踪具体算法示意图；

图3为本发明矩形网格三点Fermat射线追踪方法的流程图；

图4为本发明中P波深度域的层速度模型；

图5为本发明50×50的网格三点Fermat射线追踪示意图；

图6为本发明矩形网格三点Fermat射线追踪法的PP波最小旅行时射线路径图；

图7为本发明和矩形网格三点扰动法射线追踪运行现场的操作界面示意图和耗时对比图；

图8为本发明和矩形网格三点扰动法射线追踪方法的最小旅行时路径叠合显示图。

具体的实施方式：

图1为本发明从S点到R点，穿过界面的射线旅行示意图。在矩形网格三点扰动追踪法的基础上，本发明提出矩形网格三点Fermat射线追踪法，同样取矩形网格和射线路径的三个交点x_i-1，x_i和x_i+1，对中间点x_i，不是给出一个扰动量对其扰动，而是根据Fermat最小旅行时原理，确定具有最小旅行时的位置(r)。这样一次就能够确定出最小旅行时的位置。不是最小旅行的中间点，不必依次扰动。以新的中间点位置替代x_i，顺序取下一个交点，反复上述过程，直到最后一个交点，得到新的射线路径，计算新路径的旅行时，以新的射线路径作为新的初始射线路径，重复上述过程，直到旅行时的改变量小于给定的精度，终止追踪计算。其计算速度上比扰动法快，而且继承了扰动法的优点，适应各种复杂速度场。

$t\left(r\right)=\frac{\sqrt{{(x-r)}^2+h_1^2}}{v_1\left(r\right)}+\frac{\sqrt{r^2+h_2^2}}{v_2\left(r\right)}$

$\frac{\mathrm{dt}\left(r\right)}{\mathrm{dr}}=0$

其中，x：网格交点P_i+2-P_i水平间距；r：网格交点P_i+2-P_i+1水平间距；

h₁：网格交点P_i+1-P_i垂直间距；h₂：网格交点P_i+2-P_i+1垂直间距；

v₁，v₂分别为P_i+1点上下网格速度；t：从P_i～P_i+2的旅行时；

dt/dr：时间对r的求导。

如图1所示，上式中x＝x_i+1-x_i-1是固定两端点S(x_i-1)和R(x_i+1)的水平距离；r＝x_i+1-x_i是中间点(x_i)到R(x_i+1)的水平距离；v1、v2分别是扰动点上下网格的速度；h1，h2是扰动点上下网格的高度。t(r)是从S经过中间点x_i，到达R的射线旅行时间，如果希望从S到R这段旅行时最小，要重新确定中间点x_i的位置，三点扰动法是用一个扰动量在界面上不断改变x_i位置，计算相应的旅行时t(r)，比较各x_i的旅行时，求出最小旅行时，从而确定x_i的位置作为最小旅行时的射线路径。发明的三点Fermat方法，不是对x_i进行依次扰动，而是对旅行时t(r)求关于r的导数，当导数为零时，求解出关于r的解，对应的x_i就是最小旅行时的射线路径。这就是Fermat原理。

图2为本发明矩形网格三点Fermat射线追踪具体算法示意图。

具体方法：图中，a为初始射线，b为最小时射线。

1、连接S，R的直线作为初始射线路径，求出此直线与网格线的交点，P₁，P₂，…，P_N，并计算旅行时T₀。

2、用P₁，P₂，P₃三点修正初始射线路径，其中P₁，P₃保持不动，用Fermat最小旅行时公式，直接计算得出Q₁点的位置，使P₁Q₁P₃为P₁到P₃的时间最短的弯曲射线路径。

3、取Q₁，P₃，P₄三点作修正，同第2步，这样依次取三点修正直到修正到P_N为止，并计算新路径P₁Q₁Q₂…P_N的旅行时T₁。

4、把P₁Q₁Q₂…P_N的连线当作新的初始射线路径，重复2，3步。

5、反复重复1-4步，得到一系列弯曲射线路径和相应的旅行时序列T₀＞＝T₁＞＝…，如果某个旅行时T_M满足T_M-1-T_M＜＝ε，则追踪停止。

图3为本发明矩形网格三点Fermat射线追踪方法的流程图，其具体步骤同上。

图4为本发明中P波深度域的层速度模型图。

取射线追踪的区域为深度0-5000米，最大震源点-检波点距离4000米。按照图3的计算流程，

1、用40×50米的网格对射线追踪的区域进行网格化；

2、从地震资料中，用速度分析方法可以得到纵波(P波)和横波(S波)的叠加速度。用狄克斯公式将叠加速度换算为层速度，并将时间-速度域转换到深度-速度域，形成深度域的速度模型。在震源点-检波点(接收点)之间对速度场进行内插和外插就生成深度域射线追踪需要的二维速度模型的离散速度场(见图4)。离散速度场的速度在纵、横向变化明显，速度层中夹有高速层，高速层下有低速度层，整个速度场比较复杂；

3、给出初始射线路径，计算初始射线路径与网格交点，并计算初始射线路径的旅行时；

4、循环交点，取射线上三个顺序的交点，固定两端点，对中间点用Fermat原理计算最小旅行时对应的中间点，用新的中间点代替老的中间点；

5、判断交点是否到最后一个点？否，用新的点和顺序向下取射线交点，组成新的三个点，重复步骤步3。是，计算新路径的旅行时；

6、判断旅行时是否满足计算精度？否，以新的路径当作新的初始射线路径，重复步骤3、4和5。是，终止射线追踪。

图5为本发明50×50的网格三点Fermat射线追踪示意图；

用50×50的网格，网格间距(50米)，速度场在网格作横向和纵向线性变化。固定两点(震源点-检波点)，进行初始射线追踪和最小旅行时追踪。图5是震源点-检波点的两点射线追踪的射线图。直线是初始射线，曲线是最小旅行时射线。弯曲射线的旅行时最小，而直射线并非是最小旅行时的路径。由于速度场线性变化，所以弯曲射线十分圆滑。其中，a为初始射线，b为最小时射线。

图6为本发明矩形网格三点Fermat射线追踪法的PP波最小旅行时射线路径图。从深度点2000米至4800米，固定震源点-反射点-检波点，矩形三点Fermat射线追踪PP波射线路径图。从图6可见，射线除了明显弯曲外，没有出现盲区，所有射线都能够通过复杂速度模型，以最小旅行时到达地面接收点。

图7为本发明和矩形网格三点扰动法射线追踪运行现场的操作界面。矩形网格三点扰动法射线追踪运行现场(上)，输入扰动量为0.05，运行开始时间为8点58分41秒，深度点为2500米，PP波(PP波为纵波)的最小旅行时为1814毫秒，PS波(PS波为转换波)的最小旅行时为2522毫秒，终止时间为8点59分24秒，计算耗时43秒；本发明矩形网格三点Fermat射线追踪法运行现场(下)，运行开始时间为9点4分9秒，深度点为2500米，PP波的最小旅行时为1812毫秒，PS波的最小旅行时为2520毫秒，终止时间为9点4分29秒，计算耗时20秒。

在SGI工作站上，用矩形网格(网格间隔40×50)三点扰动法，扰动量选0.05(扰动量小于0.05时，最小旅行时的改变量不到4毫秒，但耗时明显增加。扰动量过大，如取0.1，则达不到最小旅行时的要求)，因此，扰动量0.05是经过试验后确定的量，可以即满足计算精度，又不增加无谓的运算量。在震源点，P波从三个角度入射到地下界面，在深度点2500米的界面反射P波和转换S波到地面检波点，PP波最小旅行时为1814毫秒，PS波最小旅行时为2522毫秒，计算总耗时为43秒，运算现场如图7上。同样的计算精度(甚至比扰动法精确2毫秒，图7下)、同样的深度点和震源点-检波点，矩形网格三点Fermat射线追踪法的计算总耗时为20秒。提高速度大约2.15倍。两种方法的PP波和PS波最小旅行时的射线路径见图8，实线是矩形网格三点扰动法的纵波射线追踪的结果，空心圆为矩形网格三点Fermat法纵波射线追踪法的结果；小方框是矩形网格三点扰动法的转换波射线追踪的结果，十子为矩形网格三点Fermat法转换波射线追踪法的结果。PP和PS波的下行入射P波路径相同，在界面上遵守Snall定律，反射P波和转换S波上行至检波点位置。两种方法的射线路径基本重叠。从最小旅行时的结果也可以看出，2毫秒的差别在射线路径上基本反映不出多少差别。但两种方法的运算速度差别很大，在深度2500米处，两者的计算速度差2.15倍，在浅层，如1000米处，计算结果表明，两者的计算速度差2.67倍。即在相同的扰动量条件下，界面愈浅，两种方法的计算速度差别愈大。

对于矩形网格三点扰动算法，随着深度的增加，固定一个扰动量还不能满足最小旅行时的要求。需要根据深度的增加，适当减小扰动量，这样使这种方法的耗时进一步增加。因此矩形网格三点Fermat射线追踪法比三点扰动法运算速度提高2倍多，在实际的应用中得到广泛的应用。

图8为两种射线追踪方法的最小旅行时路径叠合显示图。

图中pp-r：三点扰动法的纵波射线；pp-f：三点Fermat法的纵波射线；ps-r：三点扰动法的转换波射线；ps-f：三点Fermat法的转换波射线。

Claims (1)

1、一种用于地震勘探中射线追踪的方法，其特征在于：

所述方法包括以下步骤：

A、对射线追踪区域进行网格化操作：所述网格化操作中横向网格间隔与纵向网格间隔的比例为1∶1；

B、建立深度域的二维速度模型的离散速度场：从地震资料中，用速度分析方法得到纵波和对应深度的横波的叠加速度；再用狄克斯公式将叠加速度换算为层速度，并将时间-速度域转换到深度-速度域，形成深度域的速度模型；然后在震源点-检波点之间对速度场进行内插和外插就生成深度域射线追踪需要的二维速度模型的离散速度场；

C、对离散速度场做平滑处理：对于每个深度界面的网格，沿水平轴方向取奇数个点进行数学平均或中值滤波的平滑处理：对离散速度场作平滑处理中取5-9道的平滑；

D、设定初始射线路径和确定该初始路径与所述网格的交点p₁、p₂、…p_N；

E、计算初始射线路径的旅行时间T₀；

F、逐点追踪射线路径：

(1)取点：取初始射线路径上三个顺序交点(p₁，p₂，p₃，)；

(2)计算：固定两端点p₁和p₃，对中间点p₂，用Fermat原理计算得到最小旅行时对应的中间点q₁，并用新中间点q₁代替原有的中间点p₂，使p₁q₁p₃成为p₁到p₃时间最短的弯曲射线路径：

所述步骤F中的Fermat原理，其应用公式为：

$t\left(r\right)=\frac{\sqrt{{(x-r)}^2+h_1^2}}{v_1\left(r\right)}+\frac{\sqrt{r^2+h_2^2}}{v_2\left(r\right)}$

$\frac{\mathrm{dt}\left(r\right)}{\mathrm{dr}}=\frac{r-x}{v_1\sqrt{{(x-r)}^2+h_1^2}}+\frac{r}{v_2\sqrt{r^2+h_2^2}}=0$

其中，x：网格交点P_i+2-P_i水平间距；

      r：网格交点P_i+2-P_i+1水平间距；

      h₁：网格交点P_i+1-P_i垂直间距；

      h₂：网格交点P_i+2-P_i+1垂直间距；

    v₁，v₂分别为P_i+1点上下网格速度；

    t：从P_i～P_i+2的旅行时；

    dt/dr：时间对r的求导；

(3)再取q₁，p₃，p₄三点，对p₃进行修正，重复(2)的操作，这样依次取三点修正直到修正到射线路径的终点p_N为止；

(4)经过步骤(1)-(3)后，形成一条新的射线路径p₁q₁q₂...p_N；

(5)计算旅行时间：根据p₁q₁q₂...p_N路径计算新路径的旅行时间T₁；

(6)阈值比较：将前后两次旅行时间(T₀和T₁)的差值和阈值进行比较；

(7)重复步骤：若阈值比较的结果不符合精度要求，则以p₁q₁q₂...p_N作为初始路径，重复上述步骤F中的(1)-(6)的追踪过程，得出新旅行时间T₂，并与T₁进行阈值比较，不满足精度要求则反复重复步骤F中的(1)-(6)的追踪过程，不断形成新的射线路径，逐点逼近最短走时的射线路径；

所述步骤F中的逐点追踪递推公式：

x＝p₃-p₁、q₁＝p₁+(x-r)；

x＝p_i+2-q_i、q_i+1＝q_i+(x-r)，(i＝1，2，……，N-1)；

其中r为Fermat原理求解的具有最小旅行时的位置；

p₁、p₂，…p_N系列为初始网格交点的水平坐标；

q₁是对p₂点修改后的水平坐标；然后对p系列依次进行修改，

并利用前面修改后的q值，形成新的q系列，也就得到了从p₁至p_N新的射线路径；

G、输出射线路径步骤：

(1)若步骤F中的(6)阈值比较的结果符合精度要求，即旅行时间的差值符合精度要求：|T_M-1-T_M|＜＝ε，追踪停止；

(2)输出：相应的射线路径即为最小旅行时间的地震射线追踪路径，输出此射线路径。

Images