CN102749643B - 一种面波地震记录的频散响应获取方法及其装置 - Google Patents
一种面波地震记录的频散响应获取方法及其装置 Download PDFInfo
- Publication number
- CN102749643B CN102749643B CN201110101906.9A CN201110101906A CN102749643B CN 102749643 B CN102749643 B CN 102749643B CN 201110101906 A CN201110101906 A CN 201110101906A CN 102749643 B CN102749643 B CN 102749643B
- Authority
- CN
- China
- Prior art keywords
- partiald
- wave
- frequency
- module
- medium
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种波动方程正演的瑞利面波频散响应计算方法及其装置,通过使用高精度的高阶有限差分、合理设置自由地表边界条件、有效压制模型外边界的伪反射、使用Radon变换法来计算面波的相速度谱,即瑞利面波在频率域中的频散响应,本发明基于各向异性弹性波动方程通过时域有限差分正演模拟来计算瑞利面波频散响应特征,其地质模型可以是各向异性、非水平层状介质,自由地表可以是起伏自由地表,可以用来研究复杂介质的瑞利面波频散响应特征。
Description
技术领域
本发明涉及地球物理勘探中的地震波数值模拟领域,特别涉及一种面波地震记录的频散响应获取方法及其装置。
背景技术
弹性波理论证明,当介质为半无限弹性介质时,在自由地表和弹性介质分界面上将出现一种波,这种沿界面在弹性介质内部传播的不均匀波就是瑞利波。瑞利波由瑞利于1887年首先指出其存在而得名。瑞利波是纵波与横波在弹性介质分界面附近干涉的结果。在表层附近,质点的运动轨迹为椭圆。在均匀介质中,瑞利波的传播速度大约是横波传播速度的0.92倍。瑞利波在弹性介质中的传播深度大约为1个波长。瑞利面波的传播速度与频率、介质参数的分布有密切关系:在均匀半空间介质中,瑞利面波的传播速度与频率呈线性关系,无频散现象;当地下为多于两层的层状介质时,瑞利面波的传播速度与频率不再保持线性关系,呈现出显著的频散特征,而且呈现多组模式,即含有多条频散曲线。因此,频散响应能够一定程度上反映地下介质的分布特性,研究瑞利面波的特征最直观的方式就是研究瑞利面波的频散响应。已知地下介质参数的分布研究频散响应特征就是瑞利面波的正演问题,通过频散响应特征研究地下介质参数的分布就是瑞利面波的反演问题。在水平层状介质模型中,瑞利面波的频散响应特征表现为多条不同模式的频散曲线。利用瑞利面波的频散响应,反演地下或地球内部的结构,这在工程物探和天然地震领域得到了广泛的应用。
在探测深度介于天然地震和工程物探之间的石油地震勘探中,瑞利面波是一种很常见的干扰波,在地震剖面上常常表现为强能量、低频低速特性,在以反射波勘探为主的地震勘探中,一般在地震资料处理当作噪声干扰被压制或消除。在复杂的山地、戈壁、沙漠等地区,强能量的面波往往掩盖掉了近地表的大量反射波信息,这对地震资料的静校正和后续处理非常不利。实际上,石油地震勘探中的瑞利面波也和折射波、反射波一样包含有地下介质的信息,充分利用瑞利面波信息有可能为研究表层结构开辟一条新的途径。另外,基于模型的瑞利面波波动方程正演也可以用来消除或压制地震记录的面波干扰。因此,在石油地震勘探中,将瑞利面波当作一种有用信息,研究瑞利面波的特性,特别是研究瑞利面波的频散响应特征具有十分重要的意义。
在实现本发明的过程中,发明人发现现有技术至少存在以下缺点:
目前国内外对瑞利面波频散响应特征的研究都是基于水平层状介质中瑞利面波频散方程,该方程基于弹性波动方程推导,并假设地下介质由多层水平层状介质组成,自由地表为水平自由地表。该方法计算频散曲线快速,常被用来做面波正演和面波反演。但是,该方法计算的频散曲线只能反映相位变化,不能反映振幅强弱的变化;不能解决地下介质为各向异性、粘弹性、双相介质等更复杂介质;不能解决自由地表为起伏地表的情形;无法研究非水平层状介质问题。弹性波动方程理论上可以模拟任意复杂介质和复杂地形中的地震波传播,例如基于波动方程有限差分正演来研究纵波和横波等体波的传播规律,其相应的应用如地震叠前逆时偏移、全波形反演、基于模型正演去多次波等已经在石油地球物理勘探领域发挥了重要作用。由于瑞利面波相对于体波传播得更慢、其能量在深度方向上衰减得很快、其激发对自由地表要求严格等因素,基于波动方程使用有限差分正演来研究面波要比研究体波复杂得多。只有高精度模拟了面波传播,才可以利用面波记录来研究复杂介质中面波的频散响应特征。
目前对于面波的频散响应特征研究都是基于水平自由地表、水平层状各向同性介质,此时瑞利面波的频散方程具有解析解形式;而实际的地下介质绝大多数为起伏自由地表、非水平层状介质,频散方程很难求出解析解形式。为了研究复杂介质中瑞利面波的频散响应特征,有必要从最原始的弹性波动方程出发,来解决上面上述存在的问题。
发明内容
为了实现对复杂介质的瑞利面波频散响应特征的研究,解决地下介质为各向异性、粘弹性、双相介质等更复杂介质,解决自由地表为起伏地表的情形等非水平层状介质问题,本发明提供了一种波动方程正演来获取面波地震记录的频散响应的方法,通过使用高精度的高阶有限差分、合理设置自由地表边界条件、有效压制模型外边界的伪反射、使用相移法来计算面波的相速度谱,即瑞利面波在频率域中的频散响应。
本发明实施例提供了一种面波地震记录的频散响应获取方法,用于地球物理勘探中的地震波数值模拟领域,所述方法包括以下步骤:
步骤一,根据波动方程建立地下介质的正演模型,并将地下介质的介质参数转换为波动方程张量矩阵中的弹性系数,所述地下介质的介质参数包括垂向纵波速度、垂向横波速度、密度以及各向异性参数;
步骤二,对正演模型中的弹性波场分量和弹性系数进行二维网格离散化,所述弹性波场分量和弹性系数交错地位于网格节点上;根据波速分布的范围[Vmin,Vmax]和给定的峰值频率fpeak,设置正演模拟的时间迭代步长dt、离散网格大小和空间网格离散步长dx、dz;
步骤三,设定自由地表边界条件和不同介质间的边界条件,
所述介质密度ρ为零时,其对应的位移分量为零;
所述弹性系数C55所涉及的其中一个相邻介质的弹性系数为零时,其对应的剪应力为零;
步骤四,通过弹性波动方程进行时域有限差分正演模拟,在地面上接收包含面波的地震记录τzz(x,t;z=0),并将该地震记录变换到频率域
步骤五,对步骤四中得到的频率域进行瑞利面波速度谱计算,得到面波地震记录在频率域的速度谱s(1:nω,1:nυ),即瑞利面波在频率域中的频散响应。
进一步,所述步骤一中所建正演模型公式如下所示:
其中,介质参数包括:
Vp0——垂向纵波速度;
Vs0——垂向横波速度;
ρ——介质密度;
ε和δ——各向异性参数,
所述弹性系数同介质参数之间的关系为:
进一步,所述步骤三中的两种所述边界条件均采用介质平均法计算位于相邻两离散网格节点中间的弹性系数。
进一步,所述步骤四中对正演模型进行时间偏导数和空间位移偏导数的计算,所述时间偏导数采用二阶中心差分进行离散,所述空间位移偏导数采用十二阶中心差分进行离散;通过对正演模型进行计算,分别得出所述弹性波场分量的二阶时间差分精度和十二阶空间差分精度的数值计算迭代公式。
进一步,所述步骤四中,在正演模型的左边界区域、右边界区域和下边界区域设置完全匹配层,所述完全匹配层的厚度为8~12层,在对应的一阶空间偏导数处设置一个辅助变量,用于吸收来自外边界的反射。
进一步,所述步骤五中,根据正演模型介质中面波速度的分布范围确定频率域中需要分析的单频个数,并对每个频率通过相移法计算得出所对应的速度,通过对每个频率和速度进行循环,得到面波的地震记录在所述频率域的速度谱s(1:nω,1:nυ)。
优选地,所述步骤二中,一个波长内至少被离散成8个离散网格,一个地震子波周期至少被离散为10个时间采样点。
进一步,将频率域作振幅一致性处理,用于对低频能量和高频能量进行补偿,其实现公式如下:
另外本发明还提供了一种面波地震记录的频散响应获取装置,用于地球物理勘探中的地震波数值模拟领域,所述装置包括:
弹性系数正演模型模块,根据波动方程建立地下介质的正演模型,并将地下介质的介质参数转换为波动方程张量矩阵中的弹性系数,所述地下介质的介质参数包括垂向纵波速度、垂向横波速度、密度以及各向异性参数;
与所述弹性系数正演模型模块相连的所述正演模型网格离散模块,所述正演模型网格离散模块,用于对正演模型中的弹性波场分量和所述弹性系数正演模型模块转换得到的所述弹性系数进行二维网格离散化,并设置时间迭代步长和空间离散步长;
与所述正演模型网格离散模块相连的边界条件设置模块,所述边界条件设置模块,用于根据所述正演模型网格离散模块得到的所述二维网格设定自由地表边界条件和不同介质间的边界条件;
与所述边界条件设置模块相连的有限差分正演模拟模块,所述有限差分正演模拟模块,通过弹性波动方程进行时域有限差分正演模拟,在地面上接收包含面波的地震记录τzz(x,t;z=0),并进行地震记录的计算;
与所述有限差分正演模拟模块相连的面波记录变换域模块,所述面波记录变换域模块,用于将所述有限差分正演模拟模块计算的所述地震记录τzz(x,t;z=0)变换到频率域
及与所述面波记录变换域模块相连的面波速度谱计算模块,所述面波速度谱计算模块,通过相移法计算频率域面波速度谱,得到面波地震记录在频率域的速度谱s(1:nω,1:nυ),即瑞利面波在频率域中的频散响应;
所述边界条件设置模块包括:自由地表边界条件模块、不同介质间边界条件模块及正演模型边界匹配层模块,所述正演模型边界匹配模块用于在正演模型的左边界、右边界和下边界设置完全匹配层。
本发明实施例提供的技术方案的有益效果是:
常规方法基于瑞利面波频散方程来计算面波频散响应特征,其地质模型只能是各向同性、水平层状介质模型,自由地表只能是水平自由地表。本发明基于各向异性弹性波动方程通过时域有限差分正演模拟来计算瑞利面波频散响应特征,其地质模型可以是各向异性、非水平层状介质,自由地表可以是起伏自由地表,因此可以用来研究复杂介质的瑞利面波频散响应特征。针对水平层状介质模型的频散响应特征,如图6所示,常规方法与本发明的结果保持一致,说明了本发明的正确性、可信性。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1弹性波波动方程有限差分交错网格离散示意图;
图2本发明瑞利面波频散响应计算流程图;
图3本发明瑞利面波频散响应计算装置;
图4瑞利面波在0.5s处传播的波场快照;
图5瑞利面波在1.0s处传播的波场快照;
图6三层水平层状介质正演合成的面波记录;
图7水平层状介质模型的瑞利面波频散响应;
图8水平层状介质模型的频散曲线比较(虚线为本发明,实线为理论解);
图9各向异性介质模型;
图10各向异性介质模型的瑞利面波频散响应;
图11含起伏反射界面的层状介质模型;
图12含起伏反射界面的层状介质模型的面波频散响应;
图13含起伏自由地表的层状介质模型;
图14含起伏自由地表的层状介质模型的瑞利面波频散响应。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
图2为本发明计算方法的总体流程图,其具体实现原理如下:
在二维各向异性弹性介质中,一阶位移-应力弹性波动方程为:
式中:x为水平方向坐标变量,z为垂直方向坐标变量;ρ=ρ(x,z)为介质密度,vx=vx(x,z,t)为质点的水平位移,vz=vz(x,z,t)为垂直位移;τxx=τxx(x,z,t),τzz=τzz(x,z,t)为正应力,τxz=τxz(x,z,t)为剪应力,C11、C13、C33、C55为介质的弹性系数张量矩阵元素,各向异性介质的速度和各向异性的强度隐含在弹性参数中,因此需要转换。
根据Thomsen公式,密度、垂向纵波速度Vp0、垂向横波速度Vs0、纵波各向异性参数ε、各向异性参数δ与张量矩阵弹性系数之间的关系为:
接下来,采用交错网格离散方案,对弹性波场变量和弹性系数的空间位置在笛卡尔坐标系中进行网格离散化,并简记:
i□x+Δx,i+0.5□x+0.5Δx
k□z+Δz,k+0.5□z+0.5Δz
式中Δx、Δz分别为水平方向和深度方向的空间网格离散步长。由于弹性系数离散时位于整网格节点上,考虑到有可能涉及到错半个网格的网格节点上弹性系数的计算,在这里采用介质平均法来实现:
同时对弹性波场变量的迭代更新做以下约束:
如果ρ(i,k+0.5)=0,则vz(i,k+0.5)=0;
如果ρ(i+0.5,k)=0,则vx(i+0.5,k)=0;
如果C55(i,k)、C55(i+1,k)、C55(i,k+1)、C55(i+1,k+1)中只要其中任何一个为零,那么C55(i+0.5,k+0.5)强制为零,τxz(i+0.5,k+0.5)=0。
上述介质平均法和特殊规定就能够实现起伏自由地表边界条件和两种介质分界面处的边界条件。
在数值实现弹性波动方程时,涉及到时间偏导数和空间偏导数的计算。时间偏导数用二阶中心差分进行离散,
其中f代表弹性波场分量vx,vz,τxx,τzz,τxz。
空间偏导数用十二阶中心差分进行离散,
其中,h代表x或z,Lh[f]表示数值离散后的f对h的一阶导数。差分系数计算公式为
记fn(i,k)=f(nΔt,iΔx,kΔz)=f(t;x,z),以正演模型公式中第三个方程的数值离散为例,可表示为
同理可得正演模型公式中其它4个方程的二阶时间差分精度、十二阶空间差分精度的数值计算迭代公式。
在模型的左边界区域、右边界区域和下边界区域设置完全匹配层吸收来自外边界的反射,其规律是:为了吸收某一方向的伪反射,在对应的一阶空间偏导数处引入一个辅助变量。以正演模型公式中的第三个方程所在模型的左边界为例,在x轴设置完全匹配层的厚度为10层离散网格,并在关于x的空间偏导数项Lx[vx(i,k)]处引入一个辅助变量有:
(i=0,9,k=0,nz)
其中下一时刻的计算可以使用下面公式来进行迭代更新
(i=0,9,k=0,nz)
其中
在水平自由地表z=0处记录每个采样步长内的波场,最终我们获得了面波正演记录τzz(1:nt,1:nx)。
将该记录变换到频率域得到由于使用的子波为带限雷克子波,在面波记录的低频部分或高频部分能量比较弱,需要对低频能量或高频能量进行补偿,即振幅一致性处理,实现公式为:
确定正演模型介质中面波速度的分布范围上限υmin和下限υmax,然后等分成nυ份,则第iυ等份处的速度变量为
然后,对于每个频率,使用Radon变换法计算速度谱:
其中x=(ix-1)dx,(ix=0,1,…,nx-1)
对频率和速度进行循环,最终得到面波记录在频率域的速度谱s(1:nω,1:nυ)。在该速度谱中的振幅谱|s(1:nω,1:nυ)|中,观测能量团强弱分布,拾取能量团的峰值曲线,就可以得到面波的频散曲线
首先设计模型。该模型含有多种介质,每种介质包含垂向纵波速度Vp0、垂向横波速度Vs0和密度ρ、各向异性参数ε和δ等五个参数。为了正演计算方便,将这五个参数转换为波动方程张量矩阵中的弹性系数。根据速度分布的范围[Vmin,Vmax]和给定的峰值频率fpeak,来设置正演模拟的时间迭代步长dt和空间网格离散步长dx、dz。一个波长内至少被离散成8个网格,同时一个地震子波周期至少得离散为10个时间采样点。根据模型的大小,以及设定的时间采样步长dt和空间采样步长dx、dz,对模型进行网格离散化,离散后的网格大小为nx*nz。
然后基于弹性波动方程使用有限差分进行正演模拟,并记录面波。主要涉及四方面的内容:一是弹性波场变量和弹性参数在离散网格上的分布;二是自由地表边界条件的实现;三是有限差分的精度;四是吸收边界条件的实现。
如图1所示,弹性波场变量交错地分布于离散网格上。其中正应力分量τxx和τzz位于实线相交的网格节点上,剪应力分量τxz位于虚线相交的网格节点上,位移分量νx和νz位于实线和虚线相交的网格节点上。弹性系数ρ、C11、C13、C33和C55皆位于实线相交的网格节点上,即与正应力在相同的网格节点上。剪应力分量和位移分量所在位置如果涉及到弹性参数,就使用介质平均法来计算弹性参数。
采用介质平均法,使得剪应力波场分量位于不同介质的分界面处。同时规定凡是密度为零的地方,其对应的位移分量也为零;弹性系数C55涉及的相邻介质只要其中一个为零,其对应的剪应力也为零。起伏自由地表处的边界条件和不同介质分界处的边界条件得以自动实现。
为了保持计算效率和计算精度,时间偏导数的计算采用二阶中心差分近似,空间偏导数的计算采用十二阶中心差分近似。
为了压制有限边界模型在有限边界处引起的虚假发射,在模型的左边界、右边界和下边界设置完全匹配层来消除虚假反射。一般完全匹配层的厚度设置为10层。
通过弹性波动方程有限差分正演模拟,在地面上接收到了包含面波的地震记录τzz(x,t;z=0)。接下来,将该地震记录变换到频率域,对单个频率成分在水平方向上做归一化处理,即对振幅做归一化处理。
确定正演模型介质中面波速度的分布范围上限υmin和下限υmax,然后等分成nυ份;确定需要分析的单频个数。对于每个频率,使用相移法计算速度谱。对每个频率和速度进行循环,最终得到面波记录在频率域的速度谱s(1:nω,1:nυ)。
本发明还提供了一种面波地震记录的频散响应获取装置,所述装置包括:
弹性系数正演模型模块,根据波动方程建立地下介质的正演模型,并将地下介质的介质参数转换为波动方程张量矩阵中的弹性系数;
正演模型网格离散模块,用于对正演模型中的弹性波场分量和弹性系数进行二维网格离散化,并设置时间迭代步长和空间离散步长;
边界条件设置模块,用于设定自由地表边界条件和不同介质间的边界条件;
有限差分正演模拟模块,通过弹性波动方程进行时域有限差分正演模拟,在地面上接收包含面波的地震记录τzz(x,t;z=0),并进行地震记录的计算;
面波记录变换域模块,用于将地震记录τzz(x,t;z=0)变换到频率域
及面波速度谱计算模块,通过相移法计算频率域面波速度谱,得到面波地震记录在频率域的速度谱s(1:nω,1:nυ),即瑞利面波在频率域中的频散响应。
所述边界条件设置模块包括:自由地表边界条件模块、不同介质间边界条件模块及正演模型边界匹配层模块,所述正演模型边界匹配模块用于在正演模型的左边界、右边界和下边界设置完全匹配层。
另外如图3所示,本发明还提供了一种面波地震记录的频散响应获取装置,包括:
弹性系数正演模型模块,根据波动方程建立地下介质的正演模型,并将地下介质的介质参数转换为波动方程张量矩阵中的弹性系数;
正演模型网格离散模块,用于对正演模型中的弹性波场分量和弹性系数进行二维网格离散化,并设置时间迭代步长和空间离散步长;
边界条件设置模块,用于设定自由地表边界条件和不同介质间的边界条件;
弹性波时域有限差分正演模拟模块,通过弹性波动方程进行时域有限差分正演模拟,在地面上接收包含面波的地震记录τzz(x,t;z=0),并进行地震记录的计算;
面波记录变换域模块,用于将地震记录τzz(x,t;z=0)变换到频率域
及面波速度谱计算模块,通过相移法计算频率域面波速度谱,得到面波地震记录在频率域的速度谱s(1:nω,1:nυ),即瑞利面波在频率域中的频散响应。
其中边界条件设置模块包括:自由地表边界条件模块、不同介质间边界条件模块及正演模型边界匹配层模块,所述正演模型边界匹配模块用于在正演模型的左边界、右边界和下边界设置完全匹配层。
实施例1 水平层状介质水平自由地表的面波频散响应
设计一个含3层介质的水平层状介质模型,自由地表为水平自由地表。模型大小为为562米×100米,被离散为750*200的网格,其中网格水平间距为0.75米,网格垂直间距为0.5米;第一层垂向纵波速度为650m/s,横波速度为200m/s,密度为1.82×103Kg/m3,2个各向异性参数皆设置为0。在水平自由地表处垂直震源激发,水平自由地表接收,主频为30Hz。图4为瑞利面波在0.5s处传播的波场快照,图5为瑞利面波在1.0s处的波场快照。图6是正演模拟合成的面波记录。图7是该面波记录在频率域的频散响应。图8是常规方法和本发明计算的频散曲线比较,实现是常规方法,虚线是本发明。
实施例2 各向异性水平层状介质的面波频散响应
设计一个如图9所示的含3层介质的各向异性水平层状介质模型,自由地表为水平自由地表。与实施例1中弹性系数不同的是,各向异性参数ε为0.1、各向异性参数δ为0.5。网格大小等其它参数的设置同实施例1。图10为各向异性介质模型的瑞利面波频散响应。
实施例3 含起伏反射界面层状介质模型的面波频散响应
设计一个如图11所示的含起伏反射界面的层状介质模型,各向同性介质,自由地表为水平自由地表。共3层介质,弹性系数的具体设置同实施例1设置。图12是含起伏反射界面的层状介质模型的面波频散响应。
实施例4 含起伏自由地表层状介质模型的面波频散响应
设计一个如图13所示的含起伏自由地表的层状介质模型,各向同性介质,自由地表为起伏自由地表。共3层介质,弹性系数的具体设置同实施例1设置。图14是含起伏自由地表的层状介质模型的面波频散响应。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种面波地震记录的频散响应获取方法,用于地球物理勘探中的地震波数值模拟领域,其特征在于,所述方法包括以下步骤:
获取数据
步骤一,根据波动方程建立地下介质的正演模型,并将地下介质的介质参数转换为波动方程张量矩阵中的弹性系数,所述地下介质的介质参数包括垂向纵波速度、垂向横波速度、密度以及各向异性参数;
步骤二,对正演模型中的弹性波场分量和弹性系数进行二维网格离散化,所述弹性波场分量和弹性系数交错地位于网格节点上;根据波速分布的范围[Vmin,Vmax]和给定的峰值频率fpeak,设置正演模拟的时间迭代步长dt、离散网格大小和空间网格离散步长dx、dz;
步骤三,设定自由地表边界条件和不同介质间的边界条件,
所述介质密度ρ为零时,其对应的位移分量为零;
所述弹性系数C55所涉及的其中一个相邻介质的弹性系数为零时,其对应的剪应力为零;
步骤四,通过弹性波动方程进行时域有限差分正演模拟,在地面上接收包含面波的地震记录τzz(x,t;z=0),并将该地震记录变换到频率域
步骤五,对步骤四中得到的频率域进行瑞利面波速度谱计算,得到面波地震记录在频率域的速度谱s(1:nω,1:nυ),即瑞利面波在频率域中的频散响应。
2.根据权利要求1所述的面波地震记录的频散响应获取方法,其特征在于:
所述步骤一中所建正演模型公式如下所示:
其中,介质参数包括:
Vp0——垂向纵波速度;
Vs0——垂向横波速度;
ρ——介质密度;
ε和δ——各向异性参数,
所述弹性系数同介质参数之间的关系为:
3.根据权利要求1所述的面波地震记录的频散响应获取方法,其特征在于:
所述步骤三中的两种所述边界条件均采用介质平均法计算位于相邻两离散网格节点中间的弹性系数。
4.根据权利要求1所述的面波地震记录的频散响应获取方法,其特征在于:
所述步骤四中对正演模型进行时间偏导数和空间位移偏导数的计算,所述时间偏导数采用二阶中心差分进行离散,所述空间位移偏导数采用十二阶中心差分进行离散;通过对正演模型进行计算,分别得出所述弹性波场分量的二阶时间差分精度和十二阶空间差分精度的数值计算迭代公式。
5.根据权利要求4所述的面波地震记录的频散响应获取方法,其特征在于:
所述步骤四中,在正演模型的左边界区域、右边界区域和下边界区域设置完全匹配层,所述完全匹配层的厚度为8~12层,并在对应的一阶空间偏导数处设置一个辅助变量,用于吸收来自外边界的反射。
6.根据权利要求1所述的面波地震记录的频散响应获取方法,其特征在于:
所述步骤五中,根据正演模型介质中面波速度的分布范围确定频率域中需要分析的单频个数,并对每个频率通过Radon变换法计算得出所对应的速度,通过对每个频率和速度进行循环,得到面波的地震记录在所述频率域的速度谱s(1:nω,1:nυ)。
7.根据权利要求1-6所述的面波地震记录的频散响应获取方法,其特征在于:
所述步骤二中,一个波长内至少被离散成8个离散网格,一个地震子波周期至少被离散为10个时间采样点。
8.根据权利要求7所述的面波地震记录的频散响应获取方法,其特征在于:
将频率域作振幅一致性处理,用于对低频能量和高频能量进行补偿,其实现公式如下:
9.一种面波地震记录的频散响应获取装置,用于地球物理勘探中的地震波数值模拟领域,其特征在于,所述装置包括:
弹性系数正演模型模块,根据波动方程建立地下介质的正演模型,并将地下介质的介质参数转换为波动方程张量矩阵中的弹性系数,所述地下介质的介质参数包括垂向纵波速度、垂向横波速度、密度以及各向异性参数;
与所述弹性系数正演模型模块相连的所述正演模型网格离散模块,所述正演模型网格离散模块,用于对正演模型中的弹性波场分量和所述弹性系数正演模型模块转换得到的所述弹性系数进行二维网格离散化,并设置时间迭代步长和空间离散步长;
与所述正演模型网格离散模块相连的边界条件设置模块,所述边界条件设置模块,用于根据所述正演模型网格离散模块得到的所述二维网格设定自由地表边界条件和不同介质间的边界条件;
与所述边界条件设置模块相连的有限差分正演模拟模块,所述有限差分正演模拟模块,通过弹性波动方程进行时域有限差分正演模拟,在地面上接收包含面波的地震记录τzz(x,t;z=0),并进行地震记录的计算;
与所述有限差分正演模拟模块相连的面波记录变换域模块,所述面波记录变换域模块,用于将所述有限差分正演模拟模块计算的所述地震记录τzz(x,t;z=0)变换到频率域
及与所述面波记录变换域模块相连的面波速度谱计算模块,所述面波速度谱计算模块,通过相移法计算所述面波记录变换域模块变换得到的频率域面波速度谱,得到面波地震记录在频率域的速度谱s(1:nω,1:nυ),即瑞利面波在频率域中的频散响应;
所述边界条件设置模块包括:自由地表边界条件模块、不同介质间边界条件模块及正演模型边界匹配层模块;
所述正演模型边界匹配模块用于在正演模型的左边界、右边界和下边界设置完全匹配层。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110101906.9A CN102749643B (zh) | 2011-04-22 | 2011-04-22 | 一种面波地震记录的频散响应获取方法及其装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110101906.9A CN102749643B (zh) | 2011-04-22 | 2011-04-22 | 一种面波地震记录的频散响应获取方法及其装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102749643A CN102749643A (zh) | 2012-10-24 |
CN102749643B true CN102749643B (zh) | 2015-06-03 |
Family
ID=47029986
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110101906.9A Active CN102749643B (zh) | 2011-04-22 | 2011-04-22 | 一种面波地震记录的频散响应获取方法及其装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102749643B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104730579B (zh) * | 2013-12-18 | 2018-02-13 | 中国石油化工股份有限公司 | 一种基于表层横波速度反演的纵横波联合静校正方法 |
CN104360390B (zh) * | 2014-11-13 | 2017-02-01 | 甘肃铁道综合工程勘察院有限公司 | 基于cmpcc二维面波的吸收散射综合分析法 |
CN106054248B (zh) * | 2016-07-15 | 2017-07-18 | 河海大学 | 一种基于大面积致密储层地震岩石物理反演方法 |
CN106324668B (zh) * | 2016-08-16 | 2019-08-02 | 广东石油化工学院 | 一种基于双变地质建模技术的薄储层地震正演模拟方法 |
CN106772590B (zh) * | 2017-03-17 | 2018-10-12 | 中国地质科学院地球物理地球化学勘查研究所 | 一种剧烈起伏自由地表有限差分正演模拟系统及方法 |
CN112987090B (zh) * | 2019-12-02 | 2024-05-28 | 中国石油天然气集团有限公司 | 面波频散曲线拾取方法及装置 |
CN111143991B (zh) * | 2019-12-25 | 2023-06-30 | 国网辽宁省电力有限公司沈阳供电公司 | 一种介质包裹导线的横磁波传输模型及其构建方法 |
CN113589362B (zh) * | 2020-04-30 | 2024-03-19 | 中国石油化工股份有限公司 | 三维陆上耦合波正演模拟方法 |
CN112444875B (zh) * | 2020-10-28 | 2023-10-03 | 西安理工大学 | 一种获取场地卓越周期的精确解的方法 |
CN112698400B (zh) * | 2020-12-04 | 2023-06-23 | 中国科学院深圳先进技术研究院 | 反演方法、反演装置、计算机设备和计算机可读存储介质 |
CN112749442A (zh) * | 2020-12-25 | 2021-05-04 | 青岛黄海学院 | 一种汽车震源近地表面波模拟算法 |
CN113534255A (zh) * | 2021-07-07 | 2021-10-22 | 南方海洋科学与工程广东省实验室(湛江) | 一种任意形态不连续面自适应表达的方法 |
CN116882214B (zh) * | 2023-09-07 | 2023-12-26 | 东北石油大学三亚海洋油气研究院 | 基于dfl粘弹性方程的瑞雷波数值模拟方法及系统 |
CN116893447A (zh) * | 2023-09-11 | 2023-10-17 | 中国科学院武汉岩土力学研究所 | 一种基于瞬变时域地质雷达频散差分高效分析方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2018883C1 (ru) * | 1992-06-30 | 1994-08-30 | Борис Георгиевич Келехсаев | Устройство обработки сейсмических сигналов |
GB2383414A (en) * | 2001-12-22 | 2003-06-25 | Westerngeco Ltd | Processing seismic data |
GB2431868A (en) * | 2005-11-02 | 2007-05-09 | Wonderland Nursery Goods | Swivel limiting device and quick release coupling for stroller castor wheel |
JP2008138514A (ja) * | 2008-01-11 | 2008-06-19 | Arukoihara:Kk | 地盤調査方法および装置 |
CN101334483A (zh) * | 2008-06-13 | 2008-12-31 | 徐基祥 | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 |
WO2009081210A1 (en) * | 2007-12-20 | 2009-07-02 | Statoilhydro Asa | Method of and apparatus for exploring a region below a surface of the earth |
CN101907727A (zh) * | 2010-08-17 | 2010-12-08 | 中国科学院地质与地球物理研究所 | 一种面波多分量转换波静校正方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2365336A1 (fr) * | 2001-12-21 | 2003-06-21 | Universite De Sherbrooke | Methode et algorithme d'utilisation des ondes de surface |
-
2011
- 2011-04-22 CN CN201110101906.9A patent/CN102749643B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2018883C1 (ru) * | 1992-06-30 | 1994-08-30 | Борис Георгиевич Келехсаев | Устройство обработки сейсмических сигналов |
GB2383414A (en) * | 2001-12-22 | 2003-06-25 | Westerngeco Ltd | Processing seismic data |
GB2431868A (en) * | 2005-11-02 | 2007-05-09 | Wonderland Nursery Goods | Swivel limiting device and quick release coupling for stroller castor wheel |
WO2009081210A1 (en) * | 2007-12-20 | 2009-07-02 | Statoilhydro Asa | Method of and apparatus for exploring a region below a surface of the earth |
WO2009081150A1 (en) * | 2007-12-20 | 2009-07-02 | Statoilhydro Asa | Method of and apparatus for exploring a region below a surface of the earth |
JP2008138514A (ja) * | 2008-01-11 | 2008-06-19 | Arukoihara:Kk | 地盤調査方法および装置 |
CN101334483A (zh) * | 2008-06-13 | 2008-12-31 | 徐基祥 | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 |
CN101907727A (zh) * | 2010-08-17 | 2010-12-08 | 中国科学院地质与地球物理研究所 | 一种面波多分量转换波静校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102749643A (zh) | 2012-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102749643B (zh) | 一种面波地震记录的频散响应获取方法及其装置 | |
Robertsson et al. | Viscoelastic finite-difference modeling | |
CN103713315B (zh) | 一种地震各向异性参数全波形反演方法及装置 | |
Fornberg | The pseudospectral method: Comparisons with finite differences for the elastic wave equation | |
Wu et al. | Scattering characteristics of elastic waves by an elastic heterogeneity | |
CN104570082B (zh) | 一种基于格林函数表征的全波形反演梯度算子的提取方法 | |
Pan et al. | Love-wave waveform inversion in time domain for shallow shear-wave velocity | |
CN101334483B (zh) | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 | |
CN103257361B (zh) | 基于Zoeppritz方程近似式的油气预测方法及系统 | |
CN103389513B (zh) | 应用声波测井资料约束反演提高地震资料分辨率的方法 | |
CN102466816A (zh) | 一种叠前地震数据地层弹性常数参数反演的方法 | |
Yamanaka et al. | Effects of sedimentary layers on surface-wave propagation | |
CN103926619B (zh) | 一种三维vsp数据的逆时偏移方法 | |
CN106353797A (zh) | 一种高精度地震正演模拟方法 | |
CN108508482A (zh) | 一种地下裂缝地震散射响应特征模拟方法 | |
CN105277980A (zh) | 高精度空间和时间任意倍数可变网格有限差分正演方法 | |
CN103792573A (zh) | 一种基于频谱融合的地震波阻抗反演方法 | |
CN111007567A (zh) | 基于地震波形反演的砂泥岩薄互层预测方法及系统 | |
CN115600373A (zh) | 黏滞各向异性介质qP波模拟方法、系统、设备及应用 | |
Dhakal et al. | An evaluation of 3-D velocity models of the Kanto basin for long-period ground motion simulations | |
CN101937101B (zh) | 一种鉴定能否实施时移地震的方法 | |
Yu et al. | Comparison of different BEM+ Born series modeling schemes for wave propagation in complex geologic structures | |
Nakagawa et al. | Three-dimensional elastic wave scattering by a layer containing vertical periodic fractures | |
CN116050045A (zh) | 一种粘弹介质中地震波数据正演模拟方法以及设备 | |
CN104062680B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |