CN107966729A - 一种三维tti介质射线追踪方法及系统 - Google Patents
一种三维tti介质射线追踪方法及系统 Download PDFInfo
- Publication number
- CN107966729A CN107966729A CN201610912944.5A CN201610912944A CN107966729A CN 107966729 A CN107966729 A CN 107966729A CN 201610912944 A CN201610912944 A CN 201610912944A CN 107966729 A CN107966729 A CN 107966729A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- mfrac
- mtr
- mtd
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012512 characterization method Methods 0.000 claims abstract description 8
- 230000000644 propagated effect Effects 0.000 claims description 15
- 238000006073 displacement reaction Methods 0.000 claims description 11
- 230000004048 modification Effects 0.000 claims description 6
- 238000012986 modification Methods 0.000 claims description 6
- 238000010304 firing Methods 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims 1
- 238000003384 imaging method Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 8
- 230000005012 migration Effects 0.000 description 8
- 238000013508 migration Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000000694 effects Effects 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 125000006850 spacer group Chemical group 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
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
-
- 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
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
-
- 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
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)
- Measurement Of Radiation (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种三维TTI介质射线追踪方法及系统，该方法包括基于声波近似方程推导出三维TTI介质下的程函方程；将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式；基于解析后的程函方程建立三维TTI介质下的射线追踪方程；基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的射线旅行时。该方法能够获得精确的旅行时，利用射线追踪的结果对地震资料进行处理，能够有效解决目前地震资料处理方法精度低的问题。
Description
技术领域
本发明涉及地震资料处理领域，尤其涉及一种三维TTI介质射线追踪方法及系统。
背景技术
在地震资料处理领域，对于横向各向同性(Transverse Isotropic，TI)介质，理想情况下，若介质对称轴是垂直的，则该横向各向同性介质为VTI(Transversely IsotropyWith A Vertical Axis of Symmetry)介质，若介质对称轴是水平的，则该横向各向同性介质为HTI(Transversely Isotropy With A Horizontal Axis of Symmetry)介质。但是在实际情况中，TI介质的对称轴往往是倾斜的，此时称该横向各向同性介质称为TTI(Transversely Isotropy With A Tilt Axis of Symmetry)介质。
目前的地震资料处理中，采用基于各向同性介质和VTI介质的理论假设前提条件对地下介质中地震波的传播进行简化表达，从而导致地震偏移成像数据处理中目的层的构造成像位置错误，甚至造成无论在深度上还是纵向上皆与测井数据不吻合的现象，给后期的圈闭规模计算、井位预测埋下隐患。
因此，亟需一种算法稳健、精度较高的地震资料处理方法。
发明内容
本发明所要解决的技术问题是提供一种算法稳健、精度较高的地震资料处理方法。
为了解决上述技术问题，本发明提供了一种三维TTI介质射线追踪方法，该方法包括了以下步骤：步骤一、基于声波近似方程推导出三维TTI介质下的程函方程；步骤二、将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式；步骤三、基于解析后的程函方程建立三维TTI介质下的射线追踪方程；步骤四、基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的射线旅行时。
优选地，所述三维TTI介质下的程函方程为：
其中，τ表示射线旅行时，表示表示射线在z轴上的位移量，表示z轴方向上的射线参数的分量，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，a0,a1,a2,a3,a4均为方程的系数。
优选地，将所述三维TTI介质下的程函方程中的系数a0,a1,a2,a3,a4分别解析为以下参数方程：
其中，v表示动校正速度，η为非椭圆率，ε、δ表示Thomsen各向异性参数，θ表示TTI介质对称轴倾角，φ表示TTI介质对称轴方位角。
优选地，根据如下表达式建立所述三维TTI介质下的射线追踪方程：
其中，dx表示单位时间步长射线追踪的路径长度在x轴方向的投影，ds表示单位时间步长射线所追踪的路径长度，分别表示射线追踪过程中x,y,z轴方向射线的角度信息，表示程函方程，px,py,pz分别表示射线参数p在x,y,z轴上的分量，分别表示射线追踪单位时间步长对应的射线参数在x,y,z三个方向的改变量，分别表示射线在x,y,z轴上的位移量。
优选地，在步骤四之前还包括对所述射线追踪方程进行初始化：设置初始射线追踪的相角以及所述沿TTI介质对称轴方向传播的纵波地震波速度；根据所述初始射线追踪的相角以及所述沿TTI介质对称轴方向传播的纵波地震波速度获取相速度。
优选地，根据下式获取所述相速度V(θ′,θ,φ,φ′)：
其中，θ为TTI介质对称轴倾角，θ′为TTI介质初始出射角，φ为TTI介质对称轴方位角，φ′为射线方位角，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，相角γ表示为
系数E＝sinθ′sinθcos(φ′-φ)+cosθ′cosθ，
F＝(sinθ′cosθcos(φ′-φ)-cosθ′sinθ)2+sin2θ′sin2(φ′-φ)，
D＝[(1+2ε)F-E2]2+4(1+2δ)E2F；
ε表示Thomsen各向异性参数。
优选地，所述记录射线追踪到目标位置的射线旅行时，包括：设置追踪步长；对于地下任一目标位置，记录射线追踪的起始位置到所述目标位置射线追踪使用的总的步长数目作为射线旅行时。
本发明还提供了一种三维TTI介质射线追踪系统，所述系统包括：程函方程建模模块，其基于声波近似方程推导出三维TTI介质下的程函方程；程函方程解析模块，其将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式；射线追踪方程建模模块，其基于所述程函方程建立三维TTI介质下的射线追踪方程；射线追踪模块，其基于上述解析后的三维TTI介质下的程函方程以及射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的旅行时。
优选地，所述程函方程建模模块基于所述程函方程建立三维TTI介质下的射线追踪方程为：
其中，τ表示射线旅行时，表示表示射线在z轴上的位移量，表示z轴方向上的射线参数的分量，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，a0,a1,a2,a3,a4均为方程的系数,如以下参数方程所示：
其中，v表示动校正速度，η为非椭圆率，ε、δ表示Thomsen各向异性参数，θ表示TTI介质对称轴的倾角，φ表示TTI介质对称轴方位角。
优选地，所述射线追踪方程建模模块建立三维TTI介质下的射线追踪方程为：
其中，dx表示单位时间步长射线追踪的路径长度在x轴方向的投影，ds表示单位时间步长射线所追踪的路径长度，分别表示射线追踪过程中x,y,z轴方向射线的角度信息，表示程函方程，px,py,pz分别表示射线参数p在x,y,z轴上的分量，分别表示射线追踪单位时间步长对应的射线参数在x,y,z三个方向的改变量，分别表示射线在x,y,z轴上的位移量。
与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
本发明实施例从声波近似方程出发，推导出三维TTI介质下的程函方程，并将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式，进而基于解析后的程函方程建立三维TTI介质下的射线追踪方程并根据射线追踪方程进行射线追踪，该方法能够获得精确的旅行时，利用射线追踪的结果对地震资料进行处理，能够有效解决目前地震资料处理方法精度低的问题。
本发明的其他优点、目标，和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书，权利要求书，以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本申请的技术方案或现有技术的进一步理解，并且构成说明书的一部分。其中，表达本申请实施例的附图与本申请的实施例一起用于解释本申请的技术方案，但并不构成对本申请技术方案的限制。
图1为根据本发明一实施例的三维TTI介质射线追踪方法的流程示意图；
图2a-图2h为根据本发明一实施例的SEG管道模型示意图；
图3a-图3f为根据本发明一实施例的盐丘模型示意图；
图4为根据本发明另一实施例的三维TTI介质射线追踪系统的结构示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。需要说明的是，只要不构成冲突，本发明中的各个实施例以及各实施例中的各个特征可以相互结合，所形成的技术方案均在本发明的保护范围之内。
图1为根据本发明一实施例的三维TTI介质射线追踪方法的流程示意图，如图所示，包括以下步骤：
步骤S110、基于声波近似方程推导出三维TTI介质下的程函方程。
步骤S120、将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式。
步骤S130、基于解析后的程函方程建立三维TTI介质下的射线追踪方程。
步骤S140、基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的旅行时。
具体的，在步骤S110中，基于声波近似方程推导出的TTI介质下的程函方程，如表达式(1)所示：
其中，表示程函方程，τ表示射线旅行时，表示表示射线在z轴上的位移量，表示z轴方向上的射线参数的分量，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，a0,a1,a2,a3,a4均为方程的系数。
接下来，在步骤S120中，将上述程函方程中的a0,a1,a2,a3,a4五个参数分别解析为下列参数方程：
其中，v表示动校正速度，η为非椭圆率，ε、δ表示Thomsen各向异性参数，θ表示TTI介质对称轴倾角，φ表示TTI介质对称轴方位角。
进一步的，在上述程函方程和参数方程中，令其中pz，py，px分别表示射线参数p在x,y,z轴方向的分量，分别表示射线在x,y,z轴上的位移量。
对表达式(1)所示的程函方程进行化简，过程如下：
在中，将相应的a0，a1，a2，a3，a4五个参数方程中的用进行替换。
接下来，将a0，a1，a2，a3，a4五个参数的表达式的各项用b系列参数进行化简，其具体方法如下：
对参数a0表达式中各项用b系列参数表示为：
b02＝v2(1+2η)cos2θ(cosφpx-sinφpy)2
b04＝v2(1+2η)(sinφpx+cosφpy)2
则将a0化简为：a0＝b01-b02+b03-b04-b05+1，
同理，对参数a1表达式中各项用b系列参数表示为：
b12＝sin2θcosφpx(2vtsinφpy-1)
b13＝vt(px)2(cos2θcos2φ+sin2φ)
b14＝py(vtpy(cos2θsin2φ+cos2φ)+sin2θsinφ)
b15＝vt(2v2η(b12+b13+b14)-vt)+v2(1+2η)
则将a1化简为：a1＝b11*b15。
同理，对参数a2表达式中各项用b系列参数表示为：
b21＝-4sin2θ(3cos(2θ)+2)cos(2φ)((px)2-(py)2)
b22＝8sin2θ(3cos(2θ)+2)sin(2φ)pxpy
b23＝(2cos(2θ)+3cos(4θ)+3)((px)2+(py)2)
则将a2化简为：
a3，a4暂可不用化简，即：
在步骤S130中，基于解析后的程函方程建立三维TTI介质下的射线追踪方程的，如表达式(2)所示：
其中，dx表示单位时间步长射线追踪的路径长度在x轴方向的投影，ds表示单位时间步长射线所追踪的路径长度，分别表示射线追踪过程中x,y,z轴方向射线的角度信息，表示程函方程，px,py,pz分别表示射线参数p在x,y,z轴上的分量，分别表示射线追踪单位时间步长对应的射线参数在x,y,z轴三个方向的改变量，分别表示射线在x,y,z轴上的位移量。
在求解上述方程中的时，的求解步骤如下：
因为在常规射线追踪系统中
而依次求解方程中即可得到值，同理可求得
的具体计算过程如下：
进一步的，依次求解方程中即：
同理，具体计算过程如下：
具体计算过程如下：
在求解上述方程中的时，求解步骤如下：
因为在常规射线追踪系统
而
依次求解方程中即可得到值，同理可求得
的具体计算过程如下：
的具体计算过程如下：
的具体计算过程如下：
接下来，在得到三维TTI介质下的射线追踪方程之后，在步骤S140中，基于射线追踪方程进行射线追踪，具体包括基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的射线旅行时。
在步骤S140之前，还包括对射线追踪方程进行初始化。
具体的，设置初始射线追踪的相角以及TTI介质对称轴方向传播的纵波地震波速度vt，根据初始射线追踪的相角以及TTI介质对称轴方向传播的纵波地震波速度获取相速度，相速度如表达式(3)所示：
其中，θ为TTI介质对称轴倾角，θ′为TTI介质初始出射角，φ为TTI介质对称轴方位角，φ′为射线方位角，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，相角γ表示为：
系数E＝sinθ′sinθcos(φ′-φ)+cosθ′cosθ，
F＝(sinθ′cosθcos(φ′-φ)-cosθ′sinθ)2+sin2θ′sin2(φ′-φ)，
D＝[(1+2ε)F-E2]2+4(1+2δ)E2F；
ε表示Thomsen各向异性参数。
进一步地，射线旅行时的记录方法如下：设置固定追踪步长，对于地下任一目标位置，记录射线追踪起始位置到目标位置射线追踪总共使用的步长数目来记录旅行时。
本发明实施例从声波近似方程出发，推导出三维TTI介质下的程函方程，并将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式，进而基于解析后的程函方程建立三维TTI介质下的射线追踪方程并根据射线追踪方程进行射线追踪，该方法能够获得精确的旅行时，利用射线追踪的结果对地震资料进行处理，能够有效解决目前地震资料处理方法精度低的问题。
在一个实际的示例中，将本发明实施例提供的三维TTI介质射线追踪方法应用于TTI高斯束偏移成像，具体为通过加窗局部倾斜叠加并计算由震源和接收点出射的各向异性高斯束。基于三维TTI介质射线追踪方法获取射线参数、射线路径和射线旅行时，并根据射线旅行时获得高斯束的走时。根据射线参数、射线路径和高斯束旅行时计算动力学参量Q(s)、P(s)。根据动力学参量Q(s)、P(s)以及高斯束函数计算各条射线周围设定范围内不同频率的波场值。累加上述不同频率的波场值，并根据累加的结果进行成像。
图2a-图2h为根据本发明一实施例的SEG管道模型示意图，通过SEG管道模型可以看到利用本发明实施例的射线追踪方法进行高斯束偏移成像的效果。如图所示，一方面各向异性地层下面的地层分界面被拉平，与理论参数模型所呈现图像一致，而各向同性介质射线追踪应用于高斯束偏移成像，各向异性地层下面的地层分界面未拉平。另一方面，通过成像道集对比，各向同性成像道集道集不平，而TTI成像道集道集拉平。高斯束偏移成像与理论参数模型所呈现图像一致，成像精度高。
在另一个实际的示例中，利用本发明实施例提供的三维TTI介质射线追踪方法进行TTI参数建模，并将该TTI参数模型用于偏移成像。
具体为，利用初始参数模型进行TTI偏移，提供成像道集与偏移剖面，如果成像道集同相轴拉平，则输出初始模型参数及成像道集和偏移剖面数据，如果成像道集同相轴未被拉平，则在偏移剖面上进行属性参数拾取，在成像道集中进行剩余深度拾取，利用拾取得到的属性参数在初始参数场上进行TTI射线追踪，提供射线旅行时及路径信息，与剩余深度一起作为矩阵建立模块的输入数据构建矩阵，最终通过求解矩阵得到更新后的参数场数据。
进一步的，在进行深度拾取时，对于不同的反演参数需要的道集偏移距和角度范围是有区别的。具体的，如果单独反演vt，仅仅需要角度小于30°的道集或者偏移距深度比值小于1范围内的数据。如果单独反演ε，则需要角度小于45°的道集或者偏移距深度比值大1范围内的数据。如果同时反演三个各向异性参数vt,ε,δ，则需要对所有数据进行剩余深度拾取。对于不同的反演参数，矩阵求解过程中需要求解的对象不同。
将上述成像道集参数和偏移剖面参数用于各向异性深度偏移成像。
图3a-图3f为根据本发明一实施例的盐丘模型示意图，通过盐丘模型可以看到利用本发明实施例的射线追踪方法进行TTI参数建模偏移成像的效果。图3a-图3f依次为初始偏移剖面、参数建模之后的偏移剖面、理论值对应的偏移剖面、初始成像道集、参数建模之后的成像道集、理论值对应的成像道集，如图所示，得到的图像与理论模型精度相差不大，成像精度高。
图4为根据本发明另一实施例的三维TTI介质射线追踪系统的结构示意图，如图所示，该系统包括，
程函方程建模模块41，其执行步骤S110中的操作，程函方程解析模块42，其执行步骤S120中的操作，射线追踪方程建模模块43，其执行步骤S130中的操作，射线追踪模块44，其执行步骤S140中的操作，此处不再赘述。
虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。
Claims (10)
1.一种三维TTI介质射线追踪方法，其特征在于，所述方法包括以下步骤：
步骤一、基于声波近似方程推导出三维TTI介质下的程函方程；
步骤二、将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式；
步骤三、基于解析后的程函方程建立三维TTI介质下的射线追踪方程；
步骤四、基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的射线旅行时。
2.如权利要求1所述的方法，其特征在于，所述三维TTI介质下的程函方程为：
其中，τ表示射线旅行时，表示表示射线在z轴上的位移量，表示z轴方向上的射线参数的分量，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，a0,a1,a2,a3,a4均为方程的系数。
3.如权利要求2所述的方法，其特征在于，在步骤二中，将所述三维TTI介质下的程函方程中的系数a0,a1,a2,a3,a4分别解析为以下参数方程：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mn>2</mn>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>v</mi>
<mi>t</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>&eta;sin</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&theta;cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>2</mn>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>v</mi>
<mi>t</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>&eta;sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msubsup>
<mi>v</mi>
<mi>t</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
</mfrac>
<mi>sin</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mi>&eta;</mi>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&theta;cos</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mo>+</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>(</mo>
<mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&theta;sin</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mo>+</mo>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>&eta;</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mi>&eta;</mi>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>4</mn>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>8</mn>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mn>3</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mrow>
<mn>4</mn>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mn>3</mn>
<mo>)</mo>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mn>4</mn>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<msubsup>
<mi>v</mi>
<mi>t</mi>
<mn>2</mn>
</msubsup>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mi>&eta;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&eta;sin</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&theta;cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<msubsup>
<mi>v</mi>
<mi>t</mi>
<mn>2</mn>
</msubsup>
</mfrac>
</mrow>
其中，v表示动校正速度，η为非椭圆率，ε、δ表示Thomsen各向异性参数，θ表示TTI介质对称轴倾角，φ表示TTI介质对称轴方位角。
4.如权利要求3所述的方法，其特征在于，根据如下表达式建立所述三维TTI介质下的射线追踪方程：
其中，dx表示单位时间步长射线追踪的路径长度在x轴方向的投影，ds表示单位时间步长射线所追踪的路径长度，分别表示射线追踪过程中x,y,z轴方向射线的角度信息，表示程函方程，px,py,pz分别表示射线参数p在x,y,z轴上的分量，分别表示射线追踪单位时间步长对应的射线参数在x,y,z三个方向的改变量，分别表示射线在x,y,z轴上的位移量。
5.如权利要求1至4中任一项所述的方法，其特征在于，在步骤四之前还包括对所述射线追踪方程进行初始化：
设置初始射线追踪的相角以及所述沿TTI介质对称轴方向传播的纵波地震波速度；
根据所述初始射线追踪的相角以及所述沿TTI介质对称轴方向传播的纵波地震波速度获取相速度。
6.如权利要求5所述的方法，其特征在于，根据下式获取所述相速度V(θ′,θ,φ,φ′)：
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>&theta;</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mi>&theta;</mi>
<mo>,</mo>
<mi>&phi;</mi>
<mo>,</mo>
<msup>
<mi>&phi;</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<msqrt>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>&epsiv;</mi>
<mo>)</mo>
<mi>F</mi>
<mo>+</mo>
<msup>
<mi>E</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msqrt>
<mi>D</mi>
</msqrt>
</mrow>
</msqrt>
</mrow>
其中，θ为TTI介质对称轴倾角，θ′为TTI介质初始出射角，φ为TTI介质对称轴方位角，φ′为射线方位角，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，相角γ表示为
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&gamma;</mi>
<mo>=</mo>
<msup>
<mi>sin&theta;</mi>
<mo>&prime;</mo>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>&phi;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>cos&theta;</mi>
<mo>&prime;</mo>
</msup>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&gamma;</mi>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>sin&theta;</mi>
<mo>&prime;</mo>
</msup>
<mi>cos</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mo>-</mo>
<msup>
<mi>cos&theta;</mi>
<mo>&prime;</mo>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&theta;</mi>
<mo>&prime;</mo>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>&phi;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
系数E＝sinθ′sinθcos(φ′-φ)+cosθ′cosθ，
F＝(sinθ′cosθcos(φ′-φ)-cosθ′sinθ)2+sin2θ′sin2(φ′-φ)，
D＝[(1+2ε)F-E2]2+4(1+2δ)E2F；
ε表示Thomsen各向异性参数。
7.如权利要求1至6中任意一项所述的方法，其特征在于，所述记录射线追踪到目标位置的射线旅行时，包括：
设置追踪步长；
对于地下任一目标位置，记录射线追踪的起始位置到所述目标位置射线追踪使用的总的步长数目作为射线旅行时。
8.一种三维TTI介质射线追踪系统，其特征在于，所述系统包括：
程函方程建模模块，其基于声波近似方程推导出三维TTI介质下的程函方程；
程函方程解析模块，其将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式；
射线追踪方程建模模块，其基于解析后的程函方程建立三维TTI介质下的射线追踪方程；
射线追踪模块，基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的射线旅行时。
9.如权利要求8所述的系统，其特征在于，所述程函方程建模模块基于所述程函方程建立三维TTI介质下的射线追踪方程为：
其中，τ表示射线旅行时，表示表示射线在z轴上的位移量，表示z轴方向上的射线参数的分量，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，a0,a1,a2,a3,a4均为方程的系数,如以下参数方程所示：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mn>2</mn>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>v</mi>
<mi>t</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>&eta;sin</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&theta;cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>2</mn>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>v</mi>
<mi>t</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>&eta;sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msubsup>
<mi>v</mi>
<mi>t</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
</mfrac>
<mi>sin</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mi>&eta;</mi>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&theta;cos</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mo>+</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>(</mo>
<mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&theta;sin</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mo>+</mo>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>&eta;</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mi>&eta;</mi>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>4</mn>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>8</mn>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mn>3</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mrow>
<mn>4</mn>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mn>3</mn>
<mo>)</mo>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mn>4</mn>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<msubsup>
<mi>v</mi>
<mi>t</mi>
<mn>2</mn>
</msubsup>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mi>&eta;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&eta;sin</mi>
<mn>2</mn>
</msup>
<msup>
<mi>&theta;cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<msubsup>
<mi>v</mi>
<mi>t</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>,</mo>
</mrow>
其中，v表示动校正速度，η为非椭圆率，ε、δ表示Thomsen各向异性参数，θ表示TTI介质对称轴的倾角，φ表示TTI介质对称轴方位角。
10.如权利要求9所述的系统，其特征在于，所述射线追踪方程建模模块建立三维TTI介质下的射线追踪方程为：
其中，dx表示单位时间步长射线追踪的路径长度在x轴方向的投影，ds表示单位时间步长射线所追踪的路径长度，分别表示射线追踪过程中x,y,z轴方向射线的角度信息，表示程函方程，px,py,pz分别表示射线参数p在x,y,z轴上的分量，分别表示射线追踪单位时间步长对应的射线参数在x,y,z三个方向的改变量，分别表示射线在x,y,z轴上的位移量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610912944.5A CN107966729B (zh) | 2016-10-19 | 2016-10-19 | 一种三维tti介质射线追踪方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610912944.5A CN107966729B (zh) | 2016-10-19 | 2016-10-19 | 一种三维tti介质射线追踪方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107966729A true CN107966729A (zh) | 2018-04-27 |
| CN107966729B CN107966729B (zh) | 2019-07-19 |
Family
ID=61996348
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610912944.5A Active CN107966729B (zh) | 2016-10-19 | 2016-10-19 | 一种三维tti介质射线追踪方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107966729B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108983290A (zh) * | 2018-09-12 | 2018-12-11 | 中国石油大学(华东) | 一种三维横向各向同性介质中旅行时确定方法及系统 |
| CN110873891A (zh) * | 2018-08-30 | 2020-03-10 | 中国石油化工股份有限公司 | 基于相速度的正交各向异性介质射线追踪方法及其系统 |
| CN112505765A (zh) * | 2020-11-18 | 2021-03-16 | 东华理工大学 | Lax Friedrichs扫描地震波旅行时间的方法 |
| CN115793058A (zh) * | 2023-02-08 | 2023-03-14 | 成都理工大学 | 局部路径频变复旅行时的计算方法、装置、设备及介质 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070064530A1 (en) * | 2005-08-26 | 2007-03-22 | Ian Moore | Method for processing a record of seismic traces |
| US20070168167A1 (en) * | 2006-01-06 | 2007-07-19 | Baker Hughes Incorporated | Traveltime calculation in three dimensional transversely isotropic (3d tti) media by the fast marching method |
| US20090257308A1 (en) * | 2008-04-11 | 2009-10-15 | Dimitri Bevc | Migration velocity analysis methods |
| CN102239429A (zh) * | 2008-12-03 | 2011-11-09 | 雪佛龙美国公司 | Tti地球模型的多各向异性参数反演 |
| CN104360396A (zh) * | 2014-12-04 | 2015-02-18 | 中国海洋石油总公司 | 一种海上井间tti介质三种初至波走时层析成像方法 |
| CN104422955A (zh) * | 2013-08-22 | 2015-03-18 | 中国石油化工股份有限公司 | 一种利用旅行时变化量进行各向异性参数提取的方法 |
| CN105549081A (zh) * | 2016-01-29 | 2016-05-04 | 中国石油大学(华东) | 各向异性介质共炮域高斯束偏移成像方法 |
-
2016
- 2016-10-19 CN CN201610912944.5A patent/CN107966729B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070064530A1 (en) * | 2005-08-26 | 2007-03-22 | Ian Moore | Method for processing a record of seismic traces |
| US20070168167A1 (en) * | 2006-01-06 | 2007-07-19 | Baker Hughes Incorporated | Traveltime calculation in three dimensional transversely isotropic (3d tti) media by the fast marching method |
| US20090257308A1 (en) * | 2008-04-11 | 2009-10-15 | Dimitri Bevc | Migration velocity analysis methods |
| CN102239429A (zh) * | 2008-12-03 | 2011-11-09 | 雪佛龙美国公司 | Tti地球模型的多各向异性参数反演 |
| CN104422955A (zh) * | 2013-08-22 | 2015-03-18 | 中国石油化工股份有限公司 | 一种利用旅行时变化量进行各向异性参数提取的方法 |
| CN104360396A (zh) * | 2014-12-04 | 2015-02-18 | 中国海洋石油总公司 | 一种海上井间tti介质三种初至波走时层析成像方法 |
| CN105549081A (zh) * | 2016-01-29 | 2016-05-04 | 中国石油大学(华东) | 各向异性介质共炮域高斯束偏移成像方法 |
Non-Patent Citations (2)
| Title |
|---|
| 李晓玲 等: ""起伏层状TI介质中多次波射线追踪"", 《石油地球物理勘探》 * |
| 白海军 等: ""基于波前构建法的TTI介质射线追踪"", 《石油地球物理勘探》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110873891A (zh) * | 2018-08-30 | 2020-03-10 | 中国石油化工股份有限公司 | 基于相速度的正交各向异性介质射线追踪方法及其系统 |
| CN110873891B (zh) * | 2018-08-30 | 2021-12-24 | 中国石油化工股份有限公司 | 基于相速度的正交各向异性介质射线追踪方法及其系统 |
| CN108983290A (zh) * | 2018-09-12 | 2018-12-11 | 中国石油大学(华东) | 一种三维横向各向同性介质中旅行时确定方法及系统 |
| CN112505765A (zh) * | 2020-11-18 | 2021-03-16 | 东华理工大学 | Lax Friedrichs扫描地震波旅行时间的方法 |
| CN112505765B (zh) * | 2020-11-18 | 2023-05-09 | 东华理工大学 | Lax Friedrichs扫描地震波旅行时间的方法 |
| CN115793058A (zh) * | 2023-02-08 | 2023-03-14 | 成都理工大学 | 局部路径频变复旅行时的计算方法、装置、设备及介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107966729B (zh) | 2019-07-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107966729A (zh) | 一种三维tti介质射线追踪方法及系统 | |
| CN104133245B (zh) | 一种地震资料的静校正方法及系统 | |
| CN104267429B (zh) | 确定地层压力的方法及装置 | |
| CN105510880A (zh) | 一种基于双差法的微地震震源定位方法 | |
| CN104570069B (zh) | 一种从起伏地表直接成像的地震成像方法和装置 | |
| CN103984022B (zh) | 大斜度井纵、横波测井速度各向异性同时校正方法 | |
| CN104730579B (zh) | 一种基于表层横波速度反演的纵横波联合静校正方法 | |
| CN104422963B (zh) | 一种变偏移距vsp资料时差校正方法 | |
| CN107490808B (zh) | 一种高可靠性地震勘探观测系统的建立方法 | |
| CN105277978A (zh) | 一种确定近地表速度模型的方法及装置 | |
| CN106526677A (zh) | 一种海上自适应压制鬼波的宽频逆时偏移成像方法 | |
| CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
| CN106291687A (zh) | 各向异性多波高斯束叠前深度偏移成像方法 | |
| CN106556861A (zh) | 一种基于全方位地震资料的方位avo反演方法 | |
| CN106896403A (zh) | 弹性高斯束偏移成像方法和系统 | |
| CN107817516A (zh) | 基于初至波信息的近地表建模方法及系统 | |
| CN107703545A (zh) | 一种三分量地震槽波波场分离方法及系统 | |
| CN106970417A (zh) | 椭圆展开转换波速度分析方法与系统 | |
| CN106249297A (zh) | 基于信号估计的水力压裂微地震震源定位方法及系统 | |
| CN104482942B (zh) | 一种基于惯性系的最优两位置对准方法 | |
| CN105607119B (zh) | 近地表模型构建方法与静校正量求取方法 | |
| CN103149588A (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
| CN109917461A (zh) | 一种提高侏罗系地层及基岩成像质量的方法及装置 | |
| CN104155690B (zh) | 基于椭球展开的三维地震数据叠加速度求取方法 | |
| CN107918152B (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 |