CN105988135A - 基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法 - Google Patents
基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法 Download PDFInfo
- Publication number
- CN105988135A CN105988135A CN201510073617.0A CN201510073617A CN105988135A CN 105988135 A CN105988135 A CN 105988135A CN 201510073617 A CN201510073617 A CN 201510073617A CN 105988135 A CN105988135 A CN 105988135A
- Authority
- CN
- China
- Prior art keywords
- ray
- dimensional
- gaussian beam
- zone
- projection
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法,该基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法包括:步骤1,建立起伏地表的三维地质模型;步骤2,进行运动学射线追踪;步骤3,进行动力学射线追踪;步骤4,在三维介质中用投影菲涅尔带椭圆来约束高斯束在射线终点处能量的分布;以及步骤5,采用高斯波包法合成地震记录。本方法中的基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法在中深部地层具有更好的保幅性,对起伏地表和复杂模型具有更好的适应性。
Description
技术领域
本发明涉及地震资料处理,特别是涉及到一种基于投影菲涅尔带的三维起伏地表高斯束演模拟方法。
背景技术
近年来,国内外许多学者对起伏地表下的正演模拟和高斯束正演模拟进行了很多研究。首先,在不同介质中应用方面:Cerveny(1982)在引入高斯束理论后,将其应用于单层二维非均匀介质、侧向变速的层状介质、三维弹性非均匀介质以及三维弹性侧向变速层状介质的正演模拟中,Muler(1988)将其应用于二维非均匀吸收介质的正演模拟中,徐盛岩等(1988)将高斯束理论用于二维粘弹性介质地震波场的模拟。其次,在高斯束参数选择方面:Cerveny(1982)给出了最初的动力学追踪初始参数;Muler(1984)给出了二维介质中几种初始参数的选择,并讨论了其对高斯束性质的影响;George(1987)对束参数的选择给出了一些改进,确保了合成记录的稳定性。Cruz(2012)首先提出在二维介质中利用投影菲涅尔带半径来约束高斯束的有效半宽度,并给出了相应的初始参数,但没有给出投影菲涅尔带半径的计算方法。近年来,基于起伏地表的三维高斯束正演方法还处于起步阶段。为此我们发明了一种新的基于投影菲涅尔带的三维起伏地表高斯束演模拟方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种针对起伏地表条件下的三维高斯束正演模拟的基于投影菲涅尔带的三维起伏地表高斯束演模拟方法。
本发明的目的可通过如下技术措施来实现:基于投影菲涅尔带的三维起伏地表高斯束演模拟方法,该基于投影菲涅尔带的三维起伏地表高斯束演模拟方法包括:步骤1,建立起伏地表的三维地质模型;步骤2,进行 运动学射线追踪;步骤3,进行动力学射线追踪;步骤4,在三维介质中用投影菲涅尔带椭圆来约束高斯束在射线终点处能量的分布;以及步骤5,采用高斯波包法合成地震记录。
本发明的目的还可通过如下技术措施来实现:
在步骤2中,运动学射线追踪在三维坐标系下满足如下方程:
其中:τ为旅行时,v为中心射线的P波速度,vi(i=x,y,z)为中心射线处速度函数的偏导数,α和β分别为其坐标系下的倾角与方位角;在层内采用四阶龙格-库塔法进行求解各时间步的坐标,在界面上利用二分法求出交点的近似坐标,并利用矢量Snell定律计算反射或透射方向。
在步骤3中,将动力学追踪方程:
其中:
化为线性动力学追踪方程:
dXi/ds=HXi (i=1,2) (3)
其中: Xi为X的列向量, 0为2×2零矩阵,I为2×2单位矩阵;(s,m,n)为三维中心射线坐标,V为速度二阶偏导数矩 阵,Q,P为动力射线追踪方程的参数;
在给定四阶单位矩阵的初始条件下,得到X的通解:
X=Π(s)·C (4)
其中:Π(s)为4×4传播矩阵,其列向量为方程(3)四组线性独立的解,C为一个2×4复值初始参数矩阵;
把Π(s)和C分块如下:
其中:Πij和Ci(i,j=1,2)为2×2子矩阵;
这样,根据 得到高斯束的动力学参数矩阵:
由于射线经过界面时传播矩阵Π(s)发生突变,需要在界面处重新计算Π(s);
在射线与界面交点O点处边界条件为:
E11(O)=-sinisv-2(O)[(1+cos2is)v,z1-εcosissinisv,z3]
E12(O)=E21(O)=-sinisv-2(O)v,z2 (9)
其中:S0为炮点,O为射线与界面交点,为下一个界面交点,GT为G的转置,为方向指数;is为入射角,iR为反射角或透射角;D为界面的曲率矩阵;z1、z2和z3为局部笛卡尔坐标系中的三个分量;κ为射线中心坐标系中en和局部笛卡尔坐标系中的夹角;当发生透射时取等式(9)上面的符号,发生反射时取等式(9)下面的符号。
征在于,在步骤4中,在三维介质中用投影菲涅尔带椭圆来约束高斯束在射线终点处能量的分布;
当中心射线终点处高斯束的有效半宽度椭圆与射线的投影菲涅尔带椭圆一致时,可得:
其中:eig为特征值,Gr为中心射线终点,Hp为投影菲涅尔带矩阵,其值可以通过经典射线传播矩阵Π(s)、面面传播矩阵Τ(s)以及两者转换关系求出;
根据传播矩阵Π(s)的辛属性及数学推导可得到:
其中:λ1、λ2(λ1≥λ2)为Π11的特征值,η1、η2(η1≥η2)为Π12的特征值,ζ1、ζ2(ζ1≥ξ2)为投影菲涅尔带矩阵Hp的特征值。
在步骤5中,记以初始角(αi,βj)出发在检波点R处t时刻高斯波包为g(R,t,αi,βj),其近似表达式为:
其中:fm、γ、ν为高斯子波参数,R'为R在中心射线上的投影, qT=(q1,q2),q1,q2为射线在R点的坐标,A为振幅,在层状介质中 ;ρ为密度,v为速度,Ri为经过第i个界面反射系数,N表示射线经过界面的个数,“~”表示生成射线一侧的量值;Φ=ω/2π·|det(Q(R'))|1/2·{-det[M(R')-Re(M(R'))]},为能量叠加的权系数;
那么,检波点R处离散的能量叠加表达式为:
其中:I,J分别表示射线在倾角和方位角上离散的个数,Δα,Δβ为射线间隔。
在步骤1中,结合地震和测井资料及其解释结果包括层位、断层等, 建立模型框架,之后加上各层的速度参数,建立起伏地表的三维地质模型。
本发明中的基于投影菲涅尔带的三维起伏地表高斯束演模拟方法,引入了三维投影菲涅尔带的思想,发展了基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法,在中深部地层具有更好的保幅性,对起伏地表和复杂模型具有更好的适应性。相对于常规方法,在投影法菲涅尔带约束下高斯束的传播更稳定,能量聚集性更好;本方法使高斯束能量分布在投影菲涅尔带内,赋予了能量分布范围明确的物理意义,使得基于渐进射线理论的高斯束更符合波动理论;与常规初始参数正演结果相比,本方法在中深部地层具有更好的保幅性,对起伏地表和复杂模型具有更好的适应性。
附图说明
图1为本发明的基于投影菲涅尔带的三维起伏地表高斯束演模拟方法的一具体实施例的流程图;
图2为本发明的一具体实施例中用于测试的简单三维水平地表层状模型;
图3为本发明的一具体实施例中利用两种初始参数对三维水平地表层状模型进行三维高斯束正演模拟得到的结果的示意图;
图4为本发明的一具体实施例中用于测试的经典起伏地表台地模型;
图5为本发明的一具体实施例中起伏地表台地模型盲区示意图;
图6为本发明的一具体实施例中利用两种初始参数对起伏地表台地模型进行三维高斯束正演模拟得到的结果的示意图;
图7为本发明的一具体实施例中用于测试的起伏地表复杂模型;
图8为本发明的一具体实施例中利用两种初始参数对起伏地表复杂模型进行三维高斯束正演模拟得到的结果的示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合所附图式,作详细说明如下:
如图1所示,图1为本发明的基于投影菲涅尔带的三维起伏地表高斯 束演模拟方法的流程图。
在步骤101,建立起伏地表的三维地质模型。结合地震和测井资料及其解释结果包括层位、断层等,建立模型框架,之后加上各层的速度参数,建立起伏地表的三维地质模型。流程进入到步骤102。
在步骤102,进行运动学射线追踪。
运动学射线追踪在三维坐标系下满足如下方程:
其中:τ为旅行时,v为中心射线的P波速度,vi(i=x,y,z)为中心射线处速度函数的偏导数,α和β分别为其坐标系下的倾角与方位角;在层内采用四阶龙格-库塔法进行求解各时间步的坐标,在界面上利用二分法求出交点的近似坐标,并利用矢量Snell定律计算反射或透射方向。流程进入到步骤103。
在步骤103,进行动力学射线追踪。
将动力学追踪方程:
其中:
化为线性动力学追踪方程:
dXi/ds=HXi (i=1,2) (3)
其中: Xi为X的列向量, 0为2×2零矩阵,I为2×2单位矩阵;(s,m,n)为三维中心射线坐标,V为速度二阶偏导数矩阵,Q,P为动力射线追踪方程的参数;
在给定四阶单位矩阵的初始条件下,得到X的通解:
X=Π(s)·C (4)
其中:Π(s)为4×4传播矩阵,其列向量为方程(3)四组线性独立的解,C为一个2×4复值初始参数矩阵;
把Π(s)和C分块如下:
其中:Πij和Ci(i,j=1,2)为2×2子矩阵;
这样,根据 得到高斯束的动力学参数矩阵:
由于射线经过界面时传播矩阵Π(s)发生突变,需要在界面处重新计算Π(s);
在射线与界面交点O点处边界条件为:
E11(O)=-sinisv-2(O)[(1+cos2is)v,z1-εcosissinisv,z3]
E12(O)=E21(O)=-sinisv-2(O)v,z2 (9)
其中:S0为炮点,O为射线与界面交点,为下一个界面交点,GT为G的转置,为方向指数;is为入射角,iR为反射角或透射角;D为界面的曲率矩阵;z1、z2和z3为局部笛卡尔坐标系中的三个分量;κ为射线中心坐标系中en和局部笛卡尔坐标系中的夹角;
当发生透射时取等式(9)上面的符号,发生反射时取等式(9)下面的符号。流程进入到步骤104。
在步骤104,在三维介质中用投影菲涅尔带椭圆来约束高斯束在射线终点处能量的分布。
当中心射线终点处高斯束的有效半宽度椭圆与射线的投影菲涅尔带椭圆一致时,可得:
其中:eig为特征值,G为中心射线终点,Hp为投影菲涅尔带矩阵,其值可以通过经典射线传播矩阵Π(s)、面面传播矩阵Τ(s)以及两者转换关系求出。
根据传播矩阵Π(s)的辛属性及数学推导可得到:
其中:λ1、λ2(λ1≥λ2)为Π11的特征值,η1、η2(η1≥η2)为Π12的特征值,ζ1、ζ2(ζ1≥ξ2)为投影菲涅尔带矩阵Hp的特征值。流程进入到步骤105。
在步骤105,高斯波包法合成地震记录。
记以初始角(αi,βj)出发在检波点R处t时刻高斯波包为g(R,t,αi,βj),其近似表达式为:
其中:fm、γ、ν为高斯子波参数,R'为R在中心射线上的投影, qT=(q1,q2),q1,q2为射线在R点的坐标,A为振幅,在层状介质中 ;ρ为密度,v为速度,Ri为经过第i个界面反射系数,N表示射线经过界面的个数,“~”表示生成射线一侧的量值;Φ=ω/2π·|det(Q(R'))|1/2·{-det[M(R')-Re(M(R'))]},为能量叠加的权系数;
那么,检波点R处离散的能量叠加表达式为:
其中:I,J分别表示射线在倾角和方位角上离散的个数,Δα,Δβ为射线间隔。
在应用本发明的一具体实施例中,包括:
1)选取如图2所示的简单三维水平地表层状模型,模型尺寸为4000m×4000m×4000m,各层速度如图2所示,Δα=2°,Δβ=10°;
2)观测***为:三炮四线制(如图2,实线为炮线,虚线为检波点线),道间距25m,线间距500m,炮间距50m,炮线距500m,采样间隔4ms,记录时间4s,全排列接收。经试算得到的射线路径如图1中红色线条所示;
3)利用两种初始参数进行三维高斯束正演模拟得到的单炮记录如图3所示;
在应用本发明的另一具体实施例中,包括:
1)选取如图4所示的经典起伏地表台地模型,模型尺寸为4000m×4000m×2500m,上层速度为2000m/s,下层速度为3000m/s,Δα=2°,Δβ=2°;
2)观测***与图2所示模型相同,试算得到的射线路径如图4所示,得到射线追踪的盲区如图5;
3)利用两种初始参数进行三维高斯束正演模拟得到的单炮记录如图6所示;
在应用本发明的一具体实施例中,包括:
1)选取如图6所示的经典起伏地表台地模型,模型尺寸为4000m×4000m×4500m,接收时间5s,速度分布如图7;
2)观测***与图2所示模型相同,试算得到的射线路径如图7所示;
3)利用两种初始参数进行三维高斯束正演模拟得到的单炮记录如图8所示;
对比图3中(a)、(b)两者炮记录,(a)为常规初始参数模拟的炮记录,(b)为投影菲涅尔带约束下模拟的炮记录,可知,(1)两者同相轴的形态和位置是一致的,说明两者的运动学特征一致;(2)在纵向上,a和b同相轴主要能量分布范围(图3中矩形框)从浅层到深层都逐渐增大,但b比a增加的更慢些,这是因为投影菲涅尔带椭圆和常规有效半宽度椭圆都随射线路径增大而增大,但前者增大的速度没有后者速度快;(3)在横向上, 同一深度b中同相轴主要能量分布范围比a小,并且深度越大,两者差异越大;经过此简单模型的试算和分析,验证了基于投影菲涅尔带的高斯束正演方法的正确性,尤其是本方法还具有一定的保幅性。
图6(a)为常规初始参数模拟的炮记录,(b)为投影菲涅尔带约束下模拟的炮记录,由图5和图6可知,(1)对比两种初始参数合成的地震记录,可知本方法对起伏地表角点模型是适应的;(2)在普通射线方法盲区内(如图5中的空白区域),两种初始参数试算的炮记录都有一定的绕射能量(图6中箭头),但b中绕射能量范围(约25道)比a中(约40道)的小,并且没有a中的能量干涉现象(图6a矩形框),究其原因可能是在焦点处投影菲涅尔带范围小于常规初始参数计算的有效范围,这样使得绕射能量分布在合理的范围内。经过起伏地表台地模型的计算结果证实了本方法的对起伏地表和焦点绕射一定的适应性和保幅性。
图8(a)为常规初始参数模拟的炮记录,(b)为投影菲涅尔带约束下模拟的炮记录,对比图8中(a)、(b)两者炮记录可知,本方法对复杂构造也具有较好的适应性,对绕射能量具有很好的控制作用,因而模拟结果较常规三维高斯束正演模拟方法具有更好的保幅性。
Claims (6)
1.基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法,其特征在于,该基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法包括:
步骤1,建立起伏地表的三维地质模型;
步骤2,进行运动学射线追踪;
步骤3,进行动力学射线追踪;
步骤4,在三维介质中用投影菲涅尔带椭圆来约束高斯束在射线终点处能量的分布;以及
步骤5,采用高斯波包法合成地震记录。
2.根据权利要求1所述的基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法,其特征在于,在步骤2中,运动学射线追踪在三维坐标系下满足如下方程:
其中:τ为旅行时,v为中心射线的P波速度,vi(i=x,y,z)为中心射线处速度函数的偏导数,α和β分别为其坐标系下的倾角与方位角;在层内采用四阶龙格-库塔法进行求解各时间步的坐标,在界面上利用二分法求出交点的近似坐标,并利用矢量Snell定律计算反射或透射方向。
3.根据权利要求2所述的基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法,其特征在于,在步骤3中,将动力学追踪方程:
其中:M=PQ-1
化为线性动力学追踪方程:
dXi/ds=HXi(i=1,2) (3)
其中:Xi为X的列向量,0为2×2零矩阵,I为2×2单位矩阵;(s,m,n)为三维中心射线坐标,V为速度二阶偏导数矩阵,Q,P为动力射线追踪方程的参数;
在给定四阶单位矩阵的初始条件下,得到X的通解:
X=Π(s)·C (4)
其中:Π(s)为4×4传播矩阵,其列向量为方程(3)四组线性独立的解,C为一个2×4复值初始参数矩阵;
把Π(s)和C分块如下:
其中:Πij和Ci(i,j=1,2)为2×2子矩阵;
这样,根据得到高斯束的动力学参数矩阵:
由于射线经过界面时传播矩阵Π(s)发生突变,需要在界面处重新计算Π(s);
在射线与界面交点O点处边界条件为:
E11(O)=-sinisv-2(O)[(1+cos2is)v,z1-εcosissinisv,z3]
E12(O)=E21(O)=-sinisv-2(O)v,z2
其中:S0为炮点,O为射线与界面交点,为下一个界面交点,GT为G的转置,为方向指数;is为入射角,iR为反射角或透射角;D为界面的曲率矩阵;z1、z2和z3为局部笛卡尔坐标系中的三个分量;κ为射线中心坐标系中en和局部笛卡尔坐标系中的夹角;当发生透射时取等式(9)上面的符号,发生反射时取等式(9)下面的符号。
4.根据权利要求3所述的基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法,其特征在于,在步骤4中,在三维介质中用投影菲涅尔带椭圆来约束高斯束在射线终点处能量的分布;
当中心射线终点处高斯束的有效半宽度椭圆与射线的投影菲涅尔带椭圆一致时,可得:
其中:eig为特征值,Gr为中心射线终点,Hp为投影菲涅尔带矩阵,其值可以通过经典射线传播矩阵Π(s)、面面传播矩阵Τ(s)以及两者转换关系求出;
根据传播矩阵Π(s)的辛属性及数学推导可得到:
其中:λ1、λ2(λ1≥λ2)为Π11的特征值,η1、η2(η1≥η2)为Π12的特征值,ζ1、ζ2(ζ1≥ξ2)为投影菲涅尔带矩阵Hp的特征值。
5.根据权利要求4所述的基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法,其特征在于,在步骤5中,记以初始角(αi,βj)出发在检波点R处t时刻高斯波包为g(R,t,αi,βj),其近似表达式为:
g(R,t,αi,βj)=2πfm|ΦA|exp{-[2πfm(t-θ)γ]2+(2πfmG/γ)2-2πfmG}
(12)
×cos(2πf*(t-θ)+ν-arg(ΦA)+π/2)
其中:fm、γ、ν为高斯子波参数,R'为R在中心射线上的投影, qT=(q1,q2),q1,q2为射线在R点的坐标,A为振幅,在层状介质中
ρ为密度,v为速度,Ri为经过第i个界面反射系数,N表示射线经过界面的个数,“~”表示生成射线一侧的量值;Φ=ω/2π·|det(Q(R'))|1/2·{-det[M(R')-Re(M(R'))]},为能量叠加的权系数;
那么,检波点R处离散的能量叠加表达式为:
其中:I,J分别表示射线在倾角和方位角上离散的个数,Δα,Δβ为射线间隔。
6.根据权利要求1所述的基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法,其特征在于,在步骤1中,结合地震和测井资料及其解释结果包括层位、断层等,建立模型框架,之后加上各层的速度参数,建立起伏地表的三维地质模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510073617.0A CN105988135B (zh) | 2015-02-11 | 2015-02-11 | 基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510073617.0A CN105988135B (zh) | 2015-02-11 | 2015-02-11 | 基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105988135A true CN105988135A (zh) | 2016-10-05 |
CN105988135B CN105988135B (zh) | 2018-09-14 |
Family
ID=57041162
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510073617.0A Active CN105988135B (zh) | 2015-02-11 | 2015-02-11 | 基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105988135B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107918144A (zh) * | 2016-10-09 | 2018-04-17 | 中国石油化工股份有限公司 | 各向异性介质初至波射线追踪方法及*** |
CN109143333A (zh) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | 基于三角剖分模型的正演方法及计算机可读存储介质 |
CN109655882A (zh) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | 基于高斯束波场模拟的地震波正演方法和装置 |
CN113108897A (zh) * | 2021-04-23 | 2021-07-13 | 自然资源部第三海洋研究所 | 一种基于非均匀风关声源的海洋环境噪声场预报方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2262723C1 (ru) * | 2004-07-08 | 2005-10-20 | Завалишин Борис Родионович | Способ определения эффективных скоростей распространения сейсмических волн |
US20120010820A1 (en) * | 2010-07-08 | 2012-01-12 | Winbow Graham A | Fresnel Zone Fat Ray Tomography |
WO2014195435A1 (en) * | 2013-06-07 | 2014-12-11 | Cgg Services Sa | Regularization of multi-component seismic data |
CN104237940A (zh) * | 2014-09-29 | 2014-12-24 | 中国石油天然气股份有限公司 | 一种基于动力学特征的绕射波成像方法及装置 |
-
2015
- 2015-02-11 CN CN201510073617.0A patent/CN105988135B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2262723C1 (ru) * | 2004-07-08 | 2005-10-20 | Завалишин Борис Родионович | Способ определения эффективных скоростей распространения сейсмических волн |
US20120010820A1 (en) * | 2010-07-08 | 2012-01-12 | Winbow Graham A | Fresnel Zone Fat Ray Tomography |
WO2014195435A1 (en) * | 2013-06-07 | 2014-12-11 | Cgg Services Sa | Regularization of multi-component seismic data |
CN104237940A (zh) * | 2014-09-29 | 2014-12-24 | 中国石油天然气股份有限公司 | 一种基于动力学特征的绕射波成像方法及装置 |
Non-Patent Citations (2)
Title |
---|
JIANPING HUANG 等: "3D Gaussian Beam Forward Modelling Based on Projected Fresnel Zone for Irregular Topography", 《JOURNAL OF EARTH SCIENCE RESEARCH》 * |
刘学才 等: "三维高斯射线束法合成三分量VSP记录", 《石油地球物理勘探》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107918144A (zh) * | 2016-10-09 | 2018-04-17 | 中国石油化工股份有限公司 | 各向异性介质初至波射线追踪方法及*** |
CN107918144B (zh) * | 2016-10-09 | 2019-10-11 | 中国石油化工股份有限公司 | 各向异性介质初至波射线追踪方法及*** |
CN109143333A (zh) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | 基于三角剖分模型的正演方法及计算机可读存储介质 |
CN109655882A (zh) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | 基于高斯束波场模拟的地震波正演方法和装置 |
CN113108897A (zh) * | 2021-04-23 | 2021-07-13 | 自然资源部第三海洋研究所 | 一种基于非均匀风关声源的海洋环境噪声场预报方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105988135B (zh) | 2018-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102759746B (zh) | 一种变偏移距垂直地震剖面数据反演各向异性参数方法 | |
CN103605151B (zh) | 基于相位测量的分布式群波浅层微震定位方法 | |
CN105093274B (zh) | 一种水力压裂裂缝震源机制的反演方法及*** | |
CN105988135A (zh) | 基于投影菲涅尔带的三维起伏地表高斯束正演模拟方法 | |
CN105093292A (zh) | 一种地震成像的数据处理方法和装置 | |
CN105242318B (zh) | 一种确定砂体连通关系的方法及装置 | |
CN106353793A (zh) | 一种基于走时增量双线性插值射线追踪的井间地震层析反演方法 | |
CN105445789A (zh) | 基于多次反射折射波约束的三维菲涅尔体旅行时层析成像方法 | |
CN104237937B (zh) | 叠前地震反演方法及其*** | |
CN107843922A (zh) | 一种基于地震初至波和反射波走时联合的层析成像方法 | |
CN106646645A (zh) | 一种新的重力正演加速方法 | |
CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN105549077B (zh) | 基于多级多尺度网格相似性系数计算的微震震源定位方法 | |
CN108414983B (zh) | 一种基于逆时射线追踪方法的微地震定位技术 | |
CN101923175B (zh) | 波动方程偏移直接产生角道集方法 | |
CN109375252A (zh) | 考虑不同发震构造最大可信地震的地震动参数评价方法 | |
CN105093301B (zh) | 共成像点反射角角道集的生成方法及装置 | |
CN109375253A (zh) | 基于全部发震构造最大可信地震的地震动参数评价方法 | |
CN107290722A (zh) | 微震源的定位方法和装置 | |
CN105425299B (zh) | 确定地层裂缝分布的方法和装置 | |
CN104360396B (zh) | 一种海上井间tti介质三种初至波走时层析成像方法 | |
CN106483559A (zh) | 一种地下速度模型的构建方法 | |
CN106199704A (zh) | 一种三维三分量海底电缆地震资料速度建模方法 | |
CN102636809A (zh) | 一种传播角度域共成像点道集的生成方法 | |
CN105093279A (zh) | 针对山前带的三维地震初至波菲涅尔体层析反演方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |