公开/公告号CN109061732A
专利类型发明专利
公开/公告日2018-12-21
原文格式PDF
申请/专利权人 中国石油大学(华东);
申请/专利号CN201811173341.3
申请日2018-10-09
分类号
代理机构北京汇泽知识产权代理有限公司;
代理人程殿军
地址 266580 山东省青岛市经济技术开发区长江西路66号
入库时间 2023-06-19 07:52:37
法律状态公告日
法律状态信息
法律状态
2020-04-17
授权
授权
2019-01-15
实质审查的生效 IPC(主分类):G01V1/30 申请日:20181009
实质审查的生效
2018-12-21
公开
公开
技术领域
本发明涉及地震波散射技术领域,尤其涉及基于Rayleigh散射理论的近场散射解析方法。
背景技术
地震波散射目前主要是基于Rayleigh远场散射理论,其主要适用于深层介质的情况。现在对于浅层介质的研究也越来越多,因此浅层介质散射情况也变得尤为重要。
常规的地震波散射理论通常利用的是Rayleigh远场散射来进行研究的,针对地震波散射问题,不少地球物理工作者尝试对远场散射理论进行进一步的深入研究来改善和丰富散射理论,然而该理论忽略了近场散射的情况,只考虑了平面波入射而没有考虑球面波入射。
在浅层介质情况下,地震波是以球面波形式入射,该种情况应该采用近场散射来求解该问题,而常规的远场散射理论不再适用。
总之,现有的Rayleigh散射是远场散射的形式,并且是平面简谐波入射,该种形式只适用于深层介质的情况,在使用中具有一定的约束,而对于浅层散射的情况不能适用。在浅层介质的情况下,地震波是以球面波的形式入射的,远场散射的形式不能够应用到浅层介质的情况。
发明内容
本发明要解决的技术问题是:提出了地震波近场散射的解析方法,解决了常规Rayleigh散射难以适用于浅层介质散射情况的问题。
为了解决上述技术问题,本发明所采用的技术方案是:一种Rayleigh散射的近场散射解析方法,其包括如下步骤:
步骤1,通过对位移场方程(1)提取近场位移得到一个方程(2),
方程(1)、(2)中,α0,β0分别是纵波速度和横波速度,γi,γj,γk分别是方向余弦,其中i、j、k分别表示的是x、y、z轴方向,ρ0是介质的密度,δ是Dirac函数,r是散射点与接收点的距离,G表示格林函数,F,M分别是等价的单力和等价的力矩张量,Fi,Mjk两者的下标都表示受力的方向,Ui表示位移场,Ui下标i表示位移方向,x表示的是当前位置,x0表示的是散射体的中心位置;
步骤2,设当球面波到达散射点的时候是沿着x1方向入射的,将x、y、z三个方向分别记为x1,x2,x3,则波函数为:
方程(3)中,ω表示波的频率,t表示传播时间,u表示位移,下标表示位移方向;
步骤3,已知等效体力为:
通过积分得到
计算得到等价单力和等价力矩张量为:
方程(6)、(7)中,
步骤4,将方程(6)和(7)代入到方程(2)得到近场散射方程(8):
方程(8)中,r表示的是散射点与接收点之间的距离,α0为介质的纵波速度,β0为介质的横波速度,V表示介质的体积,
步骤5,设采用球坐标系对方程(8)进行简化,假设x1方向为极轴方向,得到:
γ1=cosθ
γ2=sinθcosφ
γ3=sinθsinφ>
方程(9)中,θ为散射射线与x1方向的夹角,φ为散射射线在yoz平面内的投影与x2方向的夹角;
步骤6,将方程(9)代入到方程(8)中得到了P波入射情况下的Rayleigh近场散射的散射场表达式为:
在方程(10)中,
本发明技术方案带来的有益效果是:
1.本发明,考虑了近场散射和球面波入射的情况,适用于浅层介质的散射问题的研究,得到了Rayleigh散射的近场散射表达式;
2.与常规的Rayleigh散射方程相比,本发明提出了地震波近场散射的解析方法,解决了常规Rayleigh散射难以适用于浅层介质散射情况的问题,丰富了地震波散射理论。
附图说明
图1等效力
图2等效力
图3等效力
图4等效力
图5等效力
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好的理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
本发明一种Rayleigh散射的近场散射解析方法包括如下步骤:
步骤1,通过对定量地震学(Aki)中的位移场方程(1)提取近场位移得到一个新的方程(2),
方程(1)、(2)中,α0,β0分别是纵波速度和横波速度,γi,γj,γk分别是方向余弦,其中i、j、k分别表示的是x、y、z轴方向,ρ0是介质的密度,δ是Dirac函数,r是散射点与接收点的距离,G表示格林函数,F,M分别是等价的单力和等价的力矩张量,Fi,Mjk两者的下标都表示受力的方向,Ui表示位移场,Ui下标i表示位移方向,并且所有参数的下标都表示受力的方向。x表示的是当前位置,x0表示的是散射体的中心位置。
步骤2,设当球面波到达散射点的时候是沿着x1方向(为了便于研究,将x、y、z三个方向分别记为x1,x2,x3)入射的,则波函数为:
其中ω表示波的频率,t表示传播时间,u表示位移,下标表示位移方向。
步骤3,计算等价的单力和等价力矩张量,
已知等效体力为:
通过积分我们可以得到
因此通过计算得到等价单力和等价力矩张量为:
其中
步骤4,将方程(6)和(7)代入到方程(2)可以得到一个新的近场散射方程:
其中,r表示的是散射点与接收点之间的距离,α0为介质的纵波速度,β0为介质的横波速度,V表示介质的体积,
步骤5,为了方便研究,我们采用球坐标系对上述公式进行简化,如果我们假设x1方向为极轴方向,因此可以得到
γ1=cosθ
γ2=sinθcosφ
γ3=sinθsinφ>
其中θ为散射射线与x1方向的夹角,φ为散射射线在yoz平面内的投影与x2方向的夹角。
步骤6,将方程(9)代入到方程(8)中得到了P波入射情况下的Rayleigh近场散射的散射场表达式为,
在方程(10)中,
如图1-5所示,对本发明方法的得到的方程(10)进行简单的分析和讨论。
图1表示的是等效体力
图2表示的是等效体力
图3表示的是等效体力
图4表示的是等效体力
图5表示的是等效体力
总结说来:图1到图5分别表示的是不同等效体力形成的空间散射图样,其不同的散射图样是方程(10)中
应用本发明方法,通过不同等效体力的空间散射图样的研究,对于近场散射有了一个更加清晰的认识和了解,对于近场散射的研究和利用具有重要作用。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。
机译: 散射型近场光学探针的制造方法,散射型近场光学探针和散射型近场光学探针的制造方法
机译: 散射式近场显微镜和散射式近场光谱系统
机译: 散射式近场显微镜和散射式近场光谱系统