CN109444955A - 三维地震射线追踪的双线性走时扰动插值方法 - Google Patents

三维地震射线追踪的双线性走时扰动插值方法 Download PDF

Info

Publication number
CN109444955A
CN109444955A CN201910020530.5A CN201910020530A CN109444955A CN 109444955 A CN109444955 A CN 109444955A CN 201910020530 A CN201910020530 A CN 201910020530A CN 109444955 A CN109444955 A CN 109444955A
Authority
CN
China
Prior art keywords
point
travel time
ray
wave front
wave
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910020530.5A
Other languages
English (en)
Other versions
CN109444955B (zh
Inventor
张建中
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ocean University of China
Original Assignee
Ocean University of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Ocean University of China filed Critical Ocean University of China
Priority to CN201910020530.5A priority Critical patent/CN109444955B/zh
Publication of CN109444955A publication Critical patent/CN109444955A/zh
Application granted granted Critical
Publication of CN109444955B publication Critical patent/CN109444955B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/67—Wave propagation modeling
    • G01V2210/671—Raytracing

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Image Generation (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种三维地震射线追踪的双线性走时扰动插值方法，属于地震勘探技术领域。本发明包括如下步骤：建立三维速度模型，设置炮点与检波点位置；从震源点开始扩展波前；计算与震源相邻的网格节点上的波传播时间，组成波前窄带；在波前窄带中按照规则找出次级震源点；计算与次级震源点相邻的网格节点处的波传播时间；不断向外扩展波前，直至计算出模型全部网格节点的波前走时；从检波点处出发，反向计算射线与离散单元边界的交点；将该交点作为新的检波点，重复上一步骤直至追踪至炮点所在单元；顺次连接炮点、射线与各单元边界的交点和检波点，就得到连接炮点与检波点的直达波射线路径；实现了三维复杂介质中地震走时和射线追踪的高精度计算。

Description

三维地震射线追踪的双线性走时扰动插值方法
技术领域
本发明涉及一种三维地震射线追踪的双线性走时扰动插值方法，属于地震勘探技术领域。
背景技术
地震射线追踪是研究地震波的一种快速有效技术，被广泛应用在地震层析成像、偏移成像、地震定位以及地震数据采集设计等领域。线性走时插值法包括二维走时插值法和三维走时插值法，是地震走时计算和和射线追踪的主要方法之一。该方法假设走时沿离散单元边界线性变化，单元边界上任意点的走时可通过相邻离散网格节点上的走时线性插值表示。
实际上，单元边界上的走时呈非线性变化，因而当离散单元较大时，线性走时插值法的线性走时假设会导致较大的走时和射线路径计算误差。为解决这一问题，提出了计算三维模型地震走时的基于规则六面体单元的线性走时扰动插值方法(Zhang et al.，2015)。该方法将单元边界的实际走时分解为等效匀速介质中的参考走时和走时扰动，假设走时扰动呈线性变化，参考走时沿单元边界非线性变化，保留了实际走时沿单元边界非线性变化的特征，具有较高的走时计算精度。但该方法未涉及射线路径的追踪，并且规则六面体单元不利于模拟地形和地层界面的起伏以及各地层速度的变化。
由于在二维不规则四边形单元中，单元边界为直线，单元边界上任意一点的走时扰动可通过相邻两个网格节点处走时扰动的线性插值求得；类似地，在三维规则六面体单元中，单元边界为与某一坐标轴垂直的水平平面或垂直平面。但是，在不规则六面体单元中，单元边界上的四个节点可能不在同一平面上，即单元边界面可能为曲面，建立在单元边界为平面基础上的走时扰动插值方程不适用于这种情况。因此，需要推导适用于三维不规则单元的走时扰动插值方法。
发明内容
针对现有技术存在的上述缺陷，本发明提出了一种三维地震射线追踪的双线性走时扰动插值方法，用于计算复杂三维模型的地震波前走时并追踪射线路径。
本发明所述的三维地震射线追踪的双线性走时扰动插值方法，包括如下步骤：
步骤S1：建立三维速度模型：根据各个地层分界面和速度变化确定各个地层中离散网格的大小，将三维模型离散为一系列不规则六面体单元，设置炮点与检波点位置；
步骤S2：从震源点开始，采用群波前扩展算法扩展波前：将震源点处波前走时值记为0，并且记K＝2，表示此处已获得波前走时；其余网格节点的波前走时值均记为无穷大，并令K＝0，表示此处网格节点未计算波前走时；
步骤S3：计算与炮点相邻的各个网格节点上的波传播时间：上述节点组成波前窄带L，令K＝1，表示该网格节点已计算出波传播时间，但还未确定波前时间；
步骤S4：在波前窄带L中按照一定规则找出次级震源点，各次级震源间的波传播能量不会相互干扰，将满足条件的网格节点从波前窄带L中移至次级震源集合B中，并令K＝2；
步骤S5：使用新提出的双线性走时扰动插值方法更新与次级震源集合B中的网格节点相邻的节点处的波传播时间，相邻节点需满足K≠2；若K＝0，则将该节点移至波前窄带L中并令K＝1；
步骤S6：若波前窄带集合L非空，则跳转至步骤S4；若集合L为空，则波前扩展结束；
步骤S7：在计算出模型中全部网格节点的波前走时后，根据费马原理和互换原理，从接收点出发，应用新提出的双线性走时扰动插值方法反向计算射线与检波点所在的离散单元各个边界的交点，并选取具有最小走时的交点作为射线与该离散单元边界的交点；
步骤S8：将该交点作为新的接收点：重复步骤S7计算射线与下一个离散单元边界的交点，若该离散单元为炮点所在单元则跳至步骤S9；
步骤S9：顺次连接炮点、射线与各单元边界的交点和接收点，便得到具有最小走时的直达波射线路径。
优选地，所述步骤S1中，模型离散过程具体按照以下步骤进行：
步骤S11：在x方向和y方向分别用一系列等间隔的垂向平面剖分模型；
步骤S12：由浅至深分别将每个地层划分成多个薄层：在同一地层中，薄层数目相等，并且在同一x、y坐标处，各薄层在z方向上的厚度相等；不同地层的薄层数目可以不同，这样就把三维速度模型离散成一系列不规则六面体单元。
优选地，所述步骤S4中，次级震源点集合B的选择需遵循如下规则：
B＝{(i,j,k)∈L,ti,j,k≤tL,min+ΔtL} (1)
其中，i，j，k分别表示离散网格节点在x，y，z三个方向的编号；L表示当前波前窄带；ti，j，k表示节点(i，j，k)处的波传播时间；tL，min表示波前窄带L中的最小波传播时间；ΔtL为波前窄带中波传播能量不会相互影响的两个相邻网格节点的波传播时间之差。
优选地，所述步骤S5中，新提出的双线性走时扰动插值方法计算步骤如下：
步骤S51：计算由炮点出发到达待计算点的射线与离散单元边界的交点：假设射线由炮点S(xS，yS，zS)出发，穿过离散单元左边界到达待计算点F(xF，yF，zF)；离散单元左边界顶点坐标为Ei(i＝1，2，3，4)，且顶点处波前走时已知为ti(i＝1，2，3，4)；射线与离散单元左边界的交点E(xE，yE，zE)可表示为：
其中，
ax＝x2-x1，bx＝x3-x1，cx＝x1-x2-x3+x4，dx＝x1 (3)
ay＝y2-y1，by＝y3-y1，cy＝y1-y2-y3+y4，dy＝y1 (4)
az＝z2-z1，bz＝z3-z1，cz＝z1-z2-z3+z4，dz＝z1 (5)
公式(2)中，m，n为点Ei(i＝1，2，3，4)坐标对E点坐标的贡献权重，可由下式表示：
式中，
at＝t2-t1-s1(h2-h1)，bt＝t3-t1-s1(h3-h1)，ct＝t1-t2-t3+t4-s1(h1-h2-h3+h4) (13)
式中，hi(i＝1，2，3，4)为点S至点Ei(i＝1，2，3，4)的直线距离；s1为炮点至单元左边界的参考慢度，有s1＝(t1/h1+t1/h1+t1/h1+t1/h1)/4；
步骤S52：利用与E点相邻的网格节点处的走时扰动，通过插值计算E点处的走时扰动，计算公式如下：
ΔtE＝atm+btn+ctmn+dt (14)
式中，m，n可通过公式(6)-(7)计算得到，at，bt，ct可由公式(13)计算得到，dt＝t1-s1h1；
步骤S53：计算走时待计算点F处的走时：假设离散单元的平均慢度为s2，射线路径在离散单元中为直线，则点F处的走时可表示为：
本发明的有益效果是：本发明所述的三维地震射线追踪的双线性走时扰动插值方法，克服了基于三维规则单元的走时扰动插值方法难以适应复杂的地层起伏和速度变化的问题，实现了三维复杂介质中地震走时和射线追踪的高精度计算；本发明易于实现，对复杂模型具有很强的适应性，且计算精度和计算效率高。。
附图说明
图1为本发明的流程框图；
图2为本发明的模型离散及波前扩展示意图；
图3为本发明的离散单元走时计算示意图；
图4(a)为新提出的双线性走时扰动插值方法计算的波前走时相对误差分布图(y＝1000m切片)；
图4(b)为三维线性走时插值法计算的波前走时相对误差分布图(y＝1000m切片)；
图5(a)为新提出的双线性走时扰动插值方法计算的射线路径与理论路径对比图；
图5(b)为三维线性走时插值法计算的射线路径与理论路径对比图；
图6为检波点接收到的射线走时的相对误差图。
具体实施方式
为了使本发明目的、技术方案更加清楚明白，下面结合实施例，对本发明作进一步详细说明。
实施例1：
如图1至图6所示，本发明所述的三维地震射线追踪的双线性走时扰动插值方法，包括如下步骤：
步骤S1：建立三维速度模型，根据各个地层分界面和速度变化确定各个地层中离散网格的大小，将三维模型离散为一系列不规则六面体单元，设置炮点与检波点位置；
步骤S2：如图2所示，从震源点开始，采用群波前扩展算法(GMM)扩展波前。将震源点处波前走时值记为0，并且记K＝2，表示该处已获得波前走时。其余网格节点的波前走时值均记为无穷大，并令K＝0，表示该处网格节点未计算波前走时；
步骤S3：计算与炮点相邻的各个网格节点上的波传播时间，这些节点组成波前窄带L，令K＝1，表示该网格节点已计算出波传播时间，但还未确定波前时间；
步骤S4：在波前窄带L中按照一定规则找出次级震源点，各次级震源间的波传播能量不会相互干扰，将满足条件的网格节点从波前窄带L中移至次级震源集合B中，并令K＝2；
步骤S5：使用新提出的双线性走时扰动插值方法更新与次级震源集合B中的网格节点相邻的节点处的波传播时间，相邻节点需满足K≠2。若K＝0，则将该节点移至波前窄带L中并令K＝1；
步骤S6：若波前窄带集合L非空，则跳转至步骤S4；若集合L为空，则波前扩展结束；
步骤S7：在计算出模型中全部网格节点的波前走时后，根据费马原理和互换原理，从接收点出发，应用新提出的双线性走时扰动插值方法反向计算射线与检波点所在的离散单元各个边界的交点，并选取具有最小走时的交点作为射线与该离散单元边界的交点；
步骤S8：将该交点作为新的接收点。重复步骤S7计算射线与下一个离散单元边界的交点，若该离散单元为炮点所在单元则跳至步骤S9；
步骤S9：顺次连接炮点、射线与各单元边界的交点和接收点，便得到具有最小走时的直达波射线路径。
进一步的，所述步骤S1中，如图2所示的模型离散过程具体按照以下步骤进行：
步骤S11：在x方向和y方向分别用一系列等间隔的垂向平面剖分模型；
步骤S12：由浅至深分别将每个地层划分成多个薄层。在同一地层中，薄层数目相等，并且在同一x、y坐标处，各薄层在z方向上的厚度相等。不同地层的薄层数目可以不同。这样,就把三维速度模型离散成一系列不规则六面体单元。
进一步的，所述步骤S4中，次级震源点集合B的选择需遵循如下规则：
B＝{(i,j,k)∈L,ti,j,k≤tL,min+ΔtL} (1)
其中，i,j,k分别表示离散网格节点在x,y,z三个方向的编号；L表示当前波前窄带；ti,j,k表示节点(i,j,k)处的波传播时间；tL，min表示波前窄带L中的最小波传播时间；ΔtL为波前窄带中波传播能量不会相互影响的两个相邻网格节点的波传播时间之差。
进一步的，所述步骤S5中，新提出的双线性走时扰动插值方法计算步骤如下：
步骤S51：计算由炮点出发到达待计算点的射线与离散单元边界的交点。为方便理解，假设射线由炮点S(xS,yS,zS)出发，穿过如图3所示的离散单元左边界到达待计算点F(xF,yF,zF)。离散单元左边界顶点坐标为Ei(i＝1,2,3,4)，且顶点处波前走时已知为ti(i＝1,2,3,4)。射线与离散单元左边界的交点E(xE,yE,zE)可表示为：
其中，
ax＝x2-x1，bx＝x3-x1，cx＝x1-x2-x3+x4，dx＝x1 (3)
ay＝y2-y1，by＝y3-y1，cy＝y1-y2-y3+y4，dy＝y1 (4)
az＝z2-z1，bz＝z3-z1，cz＝z1-z2-z3+z4，dz＝z1 (5)
公式(2)中，m,n为点Ei(i＝1,2,3,4)坐标对E点坐标的贡献权重，可由下式表示：
式中，
at＝t2-t1-s1(h2-h1)，bt＝t3-t1-s1(h3-h1)，ct＝t1-t2-t3+t4-s1(h1-h2-h3+h4) (13)
式中，hi(i＝1,2,3,4)为点S至点Ei(i＝1,2,3,4)的直线距离；s1为炮点至单元左边界的参考慢度，有s1＝(t1/h1+t1/h1+t1/h1+t1/h1)/4。
步骤S52：利用与E点相邻的网格节点处的走时扰动，通过插值计算E点处的走时扰动，计算公式如下：
ΔtE＝atm+btn+ctmn+dt (14)
式中，m,n可通过公式(4)计算得到，at,bt,ct可由公式(13)计算得到，dt＝t1-s1h1。
步骤S53：计算走时待计算点F处的走时。假设离散单元的平均慢度为s2,射线路径在离散单元中为直线，则点F处的走时可表示为
实施例2：
为了进一步说明本方法的实现思路及实现过程并证明方法的准确性，使用三维模型进行测试，并和三维线性走时插值方法的计算结果进行比较。因匀速模型的理论射线为直线，因此选用匀速模型，便于说明射线追踪和走时计算误差。
S1：建立速度为2000m/s的匀速模型，模型尺寸为3500×3500×3500m3，离散网格大小为50×50×50m3，炮点位于(1000m,1000m,3200m)，检波点共9个，以400m间距均匀分布在检波点1(3000m,200m,0m)与检波点9(3000m,3400m,0m)连接的直线上；
S2：应用群波前扩展算法从炮点开始扩展波前，利用新提出的双线性走时扰动插值方法计算模型中各个节点处的波前走时，y＝1000m切片处的波前走时相对误差分布如图4(a)所示,最大相对误差为1.0×10-11。图4(b)为三维线性走时插值方法的计算误差，最大相对误差达1.0×10-1。可以看出使用新提出的双线性走时扰动插值方法计算走时的相对误差远小于线性走时插值方法的走时相对误差。
S3：从接收点出发使用新提出的双线性走时扰动插值方法反向逐步计算射线与相应单元边界的交点，直至追踪至炮点为止。顺次连接炮点、射线与各单元边界的交点、接收点，所得到的直达波射线路径与理论射线路径的比较如图5(a)所示。图5(b)为使用三维线性走时插值方法的射线路径追踪与理论射线路径的比较，图6所示是各个检波点接收到的射线走时与理论走时的相对误差。可以看出，使用新提出的双线性走时扰动插值方法计算的射线路径和射线走时其精度均远高于三维线性走时插值方法。
本发明可广泛运用于地震勘探场合。
以上所述仅为本发明的较佳实施例而己，并不以本发明为限制，凡在本发明的精神和原则之内所作的均等修改、等同替换和改进等，均应包含在本发明的专利涵盖范围内。

Claims (4)

1.一种三维地震射线追踪的双线性走时扰动插值方法，其特征在于，包括如下步骤：
步骤S1：建立三维速度模型：根据各个地层分界面和速度变化确定各个地层中离散网格的大小，将三维模型离散为一系列不规则六面体单元，设置炮点与检波点位置；
步骤S2：从震源点开始，采用群波前扩展算法扩展波前：将震源点处波前走时值记为0，并且记K＝2，表示此处已获得波前走时；其余网格节点的波前走时值均记为无穷大，并令K＝0，表示此处网格节点未计算波前走时；
步骤S3：计算与炮点相邻的各个网格节点上的波传播时间：上述节点组成波前窄带L，令K＝1，表示该网格节点已计算出波传播时间，但还未确定波前时间；
步骤S4：在波前窄带L中按照一定规则找出次级震源点，各次级震源间的波传播能量不会相互干扰，将满足条件的网格节点从波前窄带L中移至次级震源集合B中，并令K＝2；
步骤S5：使用新提出的双线性走时扰动插值方法更新与次级震源集合B中的网格节点相邻的节点处的波传播时间，相邻节点需满足K≠2；若K＝0，则将该节点移至波前窄带L中并令K＝1；
步骤S6：若波前窄带集合L非空，则跳转至步骤S4；若集合L为空，则波前扩展结束；
步骤S7：在计算出模型中全部网格节点的波前走时后，根据费马原理和互换原理，从接收点出发，应用新提出的双线性走时扰动插值方法反向计算射线与检波点所在的离散单元各个边界的交点，并选取具有最小走时的交点作为射线与该离散单元边界的交点；
步骤S8：将该交点作为新的接收点：重复步骤S7计算射线与下一个离散单元边界的交点，若该离散单元为炮点所在单元则跳至步骤S9；
步骤S9：顺次连接炮点、射线与各单元边界的交点和接收点，便得到具有最小走时的直达波射线路径。
2.根据权利要求1所述的三维地震射线追踪的双线性走时扰动插值方法，其特征在于，
所述步骤S1中，模型离散过程具体按照以下步骤进行：
步骤S11：在x方向和y方向分别用一系列等间隔的垂向平面剖分模型；
步骤S12：由浅至深分别将每个地层划分成多个薄层：在同一地层中，薄层数目相等，并且在同一x、y坐标处，各薄层在z方向上的厚度相等；不同地层的薄层数目可以不同，这样就把三维速度模型离散成一系列不规则六面体单元。
3.根据权利要求1所述的三维地震射线追踪的双线性走时扰动插值方法，其特征在于，所述步骤S4中，次级震源点集合B的选择需遵循如下规则：
B＝{(i,j,k)∈L,ti,j,k≤tL,min+ΔtL} (1)
其中，i，j，k分别表示离散网格节点在x，y，z三个方向的编号；L表示当前波前窄带；ti，j，k表示节点(i，j，k)处的波传播时间；tL，min表示波前窄带L中的最小波传播时间；ΔtL为波前窄带中波传播能量不会相互影响的两个相邻网格节点的波传播时间之差。
4.根据权利要求3所述的三维地震射线追踪的双线性走时扰动插值方法，其特征在于，所述步骤S5中，新提出的双线性走时扰动插值方法计算步骤如下：
步骤S51：计算由炮点出发到达待计算点的射线与离散单元边界的交点：假设射线由炮点S(xS，yS，zS)出发，穿过离散单元左边界到达待计算点F(xF，yF，zF)；离散单元左边界顶点坐标为Ei(i＝1，2，3，4)，且顶点处波前走时已知为ti(i＝1，2，3，4)；射线与离散单元左边界的交点E(xE，yE，zE)可表示为：
其中，
ax＝x2-x1，bx＝x3-x1，cx＝x1-x2-x3+x4，dx＝x1 (3)
ay＝y2-y1，by＝y3-y1，cy＝y1-y2-y3+y4，dy＝y1 (4)
az＝z2-z1，bz＝z3-z1，cz＝z1-z2-z3+z4，dz＝z1 (5)
公式(2)中，m，n为点Ei(i＝1，2，3，4)坐标对E点坐标的贡献权重，可由下式表示：
式中，
at＝t2-t1-s1(h2-h1)，bt＝t3-t1-s1(h3-h1)，ct＝t1-t2-t3+t4-s1(h1-h2-h3+h4)(13)
式中，hi(i＝1，2，3，4)为点S至点Ei(i＝1，2，3，4)的直线距离；s1为炮点至单元左边界的参考慢度，有s1＝(t1/h1+t1/h1+t1/h1+t1/h1)/4；
步骤S52：利用与E点相邻的网格节点处的走时扰动，通过插值计算E点处的走时扰动，计算公式如下：
ΔtE＝atm+btn+ctmn+dt (14)
式中，m，n可通过公式(6)-(7)计算得到，at，bt，ct可由公式(13)计算得到，dt＝t1-s1h1；
步骤S53：计算走时待计算点F处的走时：假设离散单元的平均慢度为s2，射线路径在离散单元中为直线，则点F处的走时可表示为：
CN201910020530.5A 2019-01-09 2019-01-09 三维地震射线追踪的双线性走时扰动插值方法 Expired - Fee Related CN109444955B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910020530.5A CN109444955B (zh) 2019-01-09 2019-01-09 三维地震射线追踪的双线性走时扰动插值方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910020530.5A CN109444955B (zh) 2019-01-09 2019-01-09 三维地震射线追踪的双线性走时扰动插值方法

Publications (2)

Publication Number Publication Date
CN109444955A true CN109444955A (zh) 2019-03-08
CN109444955B CN109444955B (zh) 2020-05-19

Family

ID=65544058

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910020530.5A Expired - Fee Related CN109444955B (zh) 2019-01-09 2019-01-09 三维地震射线追踪的双线性走时扰动插值方法

Country Status (1)

Country Link
CN (1) CN109444955B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110660135A (zh) * 2019-09-20 2020-01-07 西南石油大学 一种利用三角网格化射线实现波前构建的方法
CN110972259A (zh) * 2019-03-29 2020-04-07 国家无线电监测中心检测中心 台站定位装置
CN112596103A (zh) * 2020-11-24 2021-04-02 中国地质科学院地球物理地球化学勘查研究所 射线追踪方法、装置和电子设备
CN114879249A (zh) * 2022-04-13 2022-08-09 中国海洋大学 基于四面体单元走时扰动插值的地震波前走时计算方法
CN114879251A (zh) * 2022-06-23 2022-08-09 中国地震局地球物理研究所 一种基于射线追踪技术的地震定位图形方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0541451A1 (fr) * 1991-11-06 1993-05-12 TOTAL Société anonyme dite : Procédé de détermination de la trajectoire d'un rai à travers un milieu géologique tridimensionnel
US5570321A (en) * 1994-03-03 1996-10-29 Atlantic Richfield Company Seismic velocity model optimization method using simulated annearling to determine prestack travel-times
US5724310A (en) * 1996-10-15 1998-03-03 Western Atlas International, Inc. Traveltime generation in the presence of velocity discontinuities
CN105093279A (zh) * 2014-05-16 2015-11-25 中国石油化工股份有限公司 针对山前带的三维地震初至波菲涅尔体层析反演方法
CN107505654A (zh) * 2017-06-23 2017-12-22 中国海洋大学 基于地震记录积分的全波形反演方法
CN107505651A (zh) * 2017-06-26 2017-12-22 中国海洋大学 地震初至波和反射波联合斜率层析成像方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0541451A1 (fr) * 1991-11-06 1993-05-12 TOTAL Société anonyme dite : Procédé de détermination de la trajectoire d'un rai à travers un milieu géologique tridimensionnel
US5570321A (en) * 1994-03-03 1996-10-29 Atlantic Richfield Company Seismic velocity model optimization method using simulated annearling to determine prestack travel-times
US5724310A (en) * 1996-10-15 1998-03-03 Western Atlas International, Inc. Traveltime generation in the presence of velocity discontinuities
CN105093279A (zh) * 2014-05-16 2015-11-25 中国石油化工股份有限公司 针对山前带的三维地震初至波菲涅尔体层析反演方法
CN107505654A (zh) * 2017-06-23 2017-12-22 中国海洋大学 基于地震记录积分的全波形反演方法
CN107505651A (zh) * 2017-06-26 2017-12-22 中国海洋大学 地震初至波和反射波联合斜率层析成像方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JIANZHONG ZHANG ET AL.: "Fast and accurate 3-D ray tracing using bilinear traveltime interpolation and the wave front group marching", 《GEOPHYSICAL JOURNAL INTERNATIONAL》 *
李同宇 等: "地震射线追踪的线性走时扰动插值算法", 《石油地球物理勘探》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110972259A (zh) * 2019-03-29 2020-04-07 国家无线电监测中心检测中心 台站定位装置
CN110660135A (zh) * 2019-09-20 2020-01-07 西南石油大学 一种利用三角网格化射线实现波前构建的方法
CN112596103A (zh) * 2020-11-24 2021-04-02 中国地质科学院地球物理地球化学勘查研究所 射线追踪方法、装置和电子设备
CN114879249A (zh) * 2022-04-13 2022-08-09 中国海洋大学 基于四面体单元走时扰动插值的地震波前走时计算方法
CN114879251A (zh) * 2022-06-23 2022-08-09 中国地震局地球物理研究所 一种基于射线追踪技术的地震定位图形方法及系统
CN114879251B (zh) * 2022-06-23 2023-02-21 中国地震局地球物理研究所 一种基于射线追踪技术的地震定位图形方法及系统

Also Published As

Publication number Publication date
CN109444955B (zh) 2020-05-19

Similar Documents

Publication Publication Date Title
CN109444955B (zh) 三维地震射线追踪的双线性走时扰动插值方法
CN105277978B (zh) 一种确定近地表速度模型的方法及装置
CN107843922A (zh) 一种基于地震初至波和反射波走时联合的层析成像方法
CN108984818A (zh) 固定翼时间域航空电磁数据拟三维空间约束整体反演方法
WO2019242045A9 (zh) 一种虚拟震源二维波前构建地震波走时计算方法
CN107817516B (zh) 基于初至波信息的近地表建模方法及系统
CN105549077B (zh) 基于多级多尺度网格相似性系数计算的微震震源定位方法
CN105445789A (zh) 基于多次反射折射波约束的三维菲涅尔体旅行时层析成像方法
CN108303736B (zh) 各向异性ti介质最短路径射线追踪正演方法
CN109459787B (zh) 基于地震槽波全波形反演的煤矿井下构造成像方法及系统
CN104360396B (zh) 一种海上井间tti介质三种初至波走时层析成像方法
Zhang et al. A method of shortest path raytracing with dynamic networks
CN116774292B (zh) 一种地震波走时确定方法、系统、电子设备及存储介质
CN111158059A (zh) 基于三次b样条函数的重力反演方法
CN111948708A (zh) 一种浸入边界起伏地表地震波场正演模拟方法
CN105425286A (zh) 地震走时获取方法及基于其的井间地震走时层析成像方法
CN105093279A (zh) 针对山前带的三维地震初至波菲涅尔体层析反演方法
CN113960532B (zh) 一种基于假想源的二次定位计算的微地震定位方法
CN103698810A (zh) 混合网最小走时射线追踪层析成像方法
CN109085642B (zh) 一种各向异性介质微地震事件定位方法
CN111665556A (zh) 地层声波传播速度模型构建方法
CN113221395A (zh) 基于分层介质的地震走时参数化网格模型构建方法及应用
Manchuk et al. Implementation aspects of sequential Gaussian simulation on irregular points
Yu et al. A minimum traveltime ray tracing global algorithm on a triangular net for propagating plane waves
CN109212594B (zh) 一种各向异性介质纵横波联合定位方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200519

Termination date: 20220109