技术领域
[0001] 本
发明涉及
地震勘探领域,尤其涉及一种基于弹性波最小二乘逆时偏移成像的方法及系统。
背景技术
[0002] 地震偏移的目的是得到能够用于后续地质解释的地下成像,且随着勘探程度的不断深入,
地震勘探对地下成像的
精度要求越来越高。
[0003] 叠前逆时偏移成像精度较高,但仅针对单一纵波的成像也可能形成地下介质成像盲区。鉴于叠前逆时偏移成像的
缺陷,弹性波逆时偏移成像被提出,且由于弹性波逆时偏移成像可形成多波模式的成像数据,使得其能够提供更为丰富的地下构造信息,即,能够提供鲁棒性更好精度更高的地下成像。但是,现有的弹性波逆时偏移成像的缺陷是:由于地球介质是弹性的,且地震资料为矢量多分量形式,而现有的弹性波逆时偏移成像中所采用的大多为标量
波动方程,难以对矢量的多分量地震数据进行偏移成像。
发明内容
[0004] 本发明要解决的技术问题是提供一种基于弹性波最小二乘逆时偏移成像的方法及系统,基于矢量弹性波方程,直接对多分量地震资料成像,使纵波和横波的反射同相轴准确归位,以提高地下成像精度。
[0005] 为解决上述技术问题,本发明提供的技术方案是:
[0006] 一方面,本发明提供一种基于弹性波最小二乘逆时偏移成像的方法,包括,[0007] 数据接收步骤,接收多分量地震数据;
[0008] 梯度构建步骤,根据预先构建的弹性波方程构建弹性波最小二乘逆时偏移的梯度;
[0009] 数据模拟步骤,根据梯度,采用逆时反偏移方式获取弹性波波场的模拟数据;
[0010] 残差计算步骤,根据多分量地震数据和模拟数据,以及根据预先构建的目标函数,获取数据残差;
[0011] 成像获取步骤,根据数据残差,采用预先构建的共轭梯度,以
迭代计算方式获取成像数据;
[0012] 成像输出步骤,输出最小二乘逆时偏移成像剖面。
[0013] 进一步地,梯度构建步骤,具体为,
[0014] 根据预先构建一阶速度-应
力方程组获取正传
震源波场;
[0015] 采用预设的存储规则,保存正传震源波场在有效边界内的速度分量波场;
[0016] 根据速度分量波场,以及根据预先构建的伴随反传方程组获取反传震源波场;
[0017] 根据伴随反传方程组,获取反传检波点波场;
[0018] 根据反传震源波场与反传检波点波场构建弹性波最小二乘逆时偏移的梯度。
[0020]
[0021] 伴随反传方程组为,
[0022]
[0023] 其中,vx为质点在
水平方向的振动速度,vz为质点在垂直方向的振动速度,τxx为水平方向的正应力,τzz为垂直方向的正应力,τxz为切应力,ρ为介质的
密度,λ和μ为介质的
拉梅常数,t为时间变量, 为伴随水平速度, 为伴随垂直速度, 为水平方向的伴随正应力, 为垂直方向的伴随正应力, 为伴随切应力,Δρ为密度的变化量,Δλ和Δμ为拉梅系数的变化量,Δvx和Δvz是模拟速度数据,vxobs和vzobs是野外观测的速度数据。
[0024] 进一步地,数据模拟步骤,具体为,根据梯度,通过求解预先构建的逆时反偏移方程组获取弹性波波场的模拟数据,且逆时反偏移方程组为,
[0025]
[0026] 其中,逆时反偏移方程组反映的是模拟数据与扰动模型参数间的关系。
[0027] 进一步地,采用时间二阶,空间十二阶的交错网格有限差分法对一阶速度-应力方程组、伴随反传方程组和逆时反偏移方程组进行求解,Δτxx为水平方向的模拟正应力,Δτzz为垂直方向的模拟正应力,Δτxz为模拟切应力。
[0028] 进一步地,目标函数为,E=||Δvx-vxobs||2+||Δvz-vzobs||2,Δvx和Δvz是模拟速度数据,vxobs和vzobs是野外观测的速度数据。
[0029] 进一步地,弹性波最小二乘逆时偏移的梯度的表达式为,
[0030]
[0031] 其中, 为密度的相对变化量, 为纵
波速度的相对变化量,为横波速度的相对变化量。
[0032] 进一步地,共轭梯度βk为 其中,k是迭代次数,gk是梯度方向,且 dk是共轭梯度方向,且
E是目标函数, 和Rρ分别是纵波速度、横波速度
和密度的相对变化量。
[0033] 另一方面,本发明还提供了一种基于弹性波最小二乘逆时偏移成像的系统,包括,[0034] 数据接收单元,用于接收多分量地震数据;
[0035] 梯度构建单元,用于根据预先构建的弹性波方程构建弹性波最小二乘逆时偏移的梯度;
[0036] 数据模拟单元,用于根据梯度,采用逆时反偏移方式获取弹性波波场的模拟数据;
[0037] 残差计算单元,用于根据多分量地震数据和模拟数据,以及根据预先构建的目标函数,获取数据残差;
[0038] 成像获取单元,用于根据数据残差,采用预先构建的共轭梯度,以迭代计算方式获取成像数据;
[0039] 成像输出单元,用于输出最小二乘逆时偏移成像剖面。
[0040] 进一步地,梯度构建单元,具体用于,
[0041] 根据预先构建一阶速度-应力方程组获取正传震源波场;
[0042] 采用预设的存储规则,保存正传震源波场在有效边界内的速度分量波场;
[0043] 根据速度分量波场,以及根据预先构建的伴随反传方程组获取反传震源波场;
[0044] 根据伴随反传方程组,获取反传检波点波场;
[0045] 根据反传震源波场与反传检波点波场构建弹性波最小二乘逆时偏移的梯度。
[0046] 本发明提供的基于弹性波最小二乘逆时偏移成像的方法及系统,基于矢量弹性波方程,直接对多分量地震资料成像,使纵波和横波的反射同相轴准确归位,以提高地下成像精度,鲁棒性好。
附图说明
[0047] 图1是本发明
实施例提供的基于弹性波最小二乘逆时偏移成像的方法的
流程图;
[0048] 图2是本发明实施例提供的基于弹性波最小二乘逆时偏移成像的方法的又一流程图;
[0049] 图3是本发明实施例提供的正传震源波场的边界存储示意图;
[0050] 图4是本发明实施例提供的正传的震源波场的垂直速度分量波场;
[0051] 图5是本发明实施例提供的重建的震源波场的垂直速度分量波场;
[0052] 图6是本发明实施例提供的震源波场的垂直速度分量波场的重建误差分析图;
[0053] 图7是本发明实施例提供的成像结果分析图;
[0054] 图8是本发明实施例提供的又一成像结果分析图;
[0055] 图9是本发明实施例提供的基于弹性波最小二乘逆时偏移成像的系统的
框图。
具体实施方式
[0056] 下面通过具体的实施例进一步说明本发明,但是,应当理解为,这些实施例仅仅是用于更详细具体地说明之用,而不应理解为用于以任何形式限制本发明。
[0057] 实施例一
[0058] 结合图1,本实施例提供的基于弹性波最小二乘逆时偏移成像的方法,包括,[0059] 数据接收步骤S1,接收多分量地震数据;
[0060] 梯度构建步骤S2,根据预先构建的弹性波方程构建弹性波最小二乘逆时偏移的梯度;
[0061] 数据模拟步骤S3,根据梯度,采用逆时反偏移方式获取弹性波波场的模拟数据;
[0062] 残差计算步骤S4,根据多分量地震数据和模拟数据,以及根据预先构建的目标函数,获取数据残差;
[0063] 成像获取步骤S5,根据数据残差,采用预先构建的共轭梯度,以迭代计算方式获取成像数据;
[0064] 成像输出步骤S6,输出最小二乘逆时偏移成像剖面。
[0065] 本发明实施例提供的基于弹性波最小二乘逆时偏移成像的方法,基于矢量弹性波方程,直接对多分量地震资料成像,使纵波和横波的反射同相轴准确归位,以提高地下成像精度,鲁棒性好。
[0066] 优选地,如图2所示地,梯度构建步骤S2,具体为,
[0067] 根据预先构建阶速度-应力方程组获取正传震源波场;
[0068] 采用预设的存储规则,保存正传震源波场在有效边界内的速度分量波场;
[0069] 根据速度分量波场,以及根据预先构建的伴随反传方程组获取反传震源波场;
[0070] 根据伴随反传方程组,获取反传检波点波场;
[0071] 根据反传震源波场与反传检波点波场构建弹性波最小二乘逆时偏移的梯度。
[0072] 本实施例中,采用如图3所示的存储规则对正传震源波场进行存储。此外,图4是本实施例提供的正传的震源波场的垂直速度分量波场,以及图5是本实施例提供的重建的震源波场的垂直速度分量波场,根据图4和图5,可以获知,本实施例提供的方法能够很好地对震源波场进行重建与恢复,此外,更加具体地,针对图5中存在的重建误差,进行了更加具体的分析,且图6是图5中对应的x=0.8km处的重建误差分析图,且图6中的第一列图为正传的垂直速度分量波场的单道,第二列图为重建的垂直速度分量波场的单道,第三列图为二者的绝对误差,从图6中可见,误差极小,也就是说,本实施例提供的方法能够很好地对震源波场进行重建与恢复。
[0073] 此外,需要说明的是,图7是层状介质图,图7中的第一列的上下两图分别为真实的纵波速度相对变化量与真实的横波速度相对变化量,第二列的上下两图分别为采用本实施例的方法成像获取的纵波速度的相对变化量(即, )和横波速度的相对变化量(即, )。图8是地下形态图,图8中的第一列的上下两图分别为真实的纵波速度相对变化量与真实的横波速度相对变化量,第二列的上下两图分别为采用本实施例的方法成像获取的纵波速度的相对变化量(即, )和横波速度的相对变化量(即, )。图7和图8具体的阐述了本实施例提供的方法的有效性与实用性。
[0074] 具体地,一阶速度-应力方程组为,
[0075]
[0076] 伴随反传方程组为,
[0077]
[0078] 其中,vx为质点在水平方向的振动速度,vz为质点在垂直方向的振动速度,τxx为水平方向的正应力,τzz为垂直方向的正应力,τxz为切应力,ρ为介质的密度,λ和μ为介质的拉梅常数,t为时间变量, 为伴随水平速度, 为伴随垂直速度, 为水平方向的伴随正应力, 为垂直方向的伴随正应力, 为伴随切应力,Δρ为密度的变化量,Δλ和Δμ为拉梅系数的变化量,Δvx和Δvz是模拟速度数据,vxobs和vzobs是野外观测的速度数据。
[0079] 需要说明的是,本实施例中,一阶速度-应力方程组是描述弹性波场在地下的传播方程,且一阶速度-应力方程组必须满足伴随反传方程组。此外,本实施例中,采用时间二阶,空间十二阶的交错网格有限差分法对一阶速度-应力方程组及其伴随方程组进行求解。此外,需要说明的是,通过求解以上二组方程组能够得到正传的震源波场和伴随反传的检波点波场,且在获得正传的震源波场和伴随反传的检波点波场后,可应用梯度的表达式获取目标函数随成像结果的梯度。
[0080] 更加具体地,弹性波最小二乘逆时偏移的梯度的表达式为,
[0081]
[0082] 其中, 为密度的相对变化量, 为纵波速度的相对变化量,为横波速度的相对变化量。
[0083] 此外,需要说明的是,本实施例中的弹性波最小二乘逆时偏移的梯度是由中间梯度推导获得,且中间梯度的表达式为:
[0084]
[0085] 进一步优选地,数据模拟步骤,具体为,根据梯度,通过求解预先构建的逆时反偏移方程组获取弹性波波场的模拟数据,且逆时反偏移方程组为,
[0086]
[0087] 其中,逆时反偏移方程组反映的是模拟数据与扰动模型参数间的关系,Δτxx为水平方向的模拟正应力,Δτzz为垂直方向的模拟正应力,Δτxz为模拟切应力。
[0088] 需要说明的是,本实施例中,根据扰动理论,从弹性波的一阶速度-应力方程组能够推导出逆时反偏移方程组,它描述了模拟数据和扰动模型参数之间的关系。具体地,本实施例中,采用时间二阶,空间十二阶的交错网格有限差分法对逆时反偏移方程组进行求解。此外,逆时反偏移方程组的意义在于:该方程组建立了模拟数据和扰动模型参数(它的作用等同于反射率成像结果)之间的关系,基于该方程组,能够使用模拟数据去拟合观测数据,且当模拟数据和观测数据之间的二范数残差值最小时,此时,扰动模型参数(它的作用等同于反射率成像结果)即为最优的成像结果。
[0089] 进一步优选地,目标函数为,E=||Δvx-vxobs||2+||Δvz-vzobs||2,Δvx和Δvz是模拟速度数据,vxobs和vzobs是野外观测的速度数据。
[0090] 进一步优选地,共轭梯度βk为其中,k是迭代次数,gk是梯度方向,且 dk是共轭梯度方向,且
E是目标函数, 和Rρ分别是纵波速度、横波速度和
密度的相对变化量。
[0091] 需要说明的是,本实施例中,
现有技术中常用的直接采用梯度下降方法,收敛速度很慢,使用本实施例的共轭梯度进行计算,则可实现在仅增加一点点计算量的情况下获取明显的速度提升,实用性高。
[0092] 实施例二
[0093] 结合图9,本,实施例提供的基于弹性波最小二乘逆时偏移成像的系统,包括,[0094] 数据接收单元1,用于接收多分量地震数据;
[0095] 梯度构建单元2,用于根据预先构建的弹性波方程构建弹性波最小二乘逆时偏移的梯度;
[0096] 数据模拟单元3,用于根据梯度,采用逆时反偏移方式获取弹性波波场的模拟数据;
[0097] 残差计算单元4,用于根据多分量地震数据和模拟数据,以及根据预先构建的目标函数,获取数据残差;
[0098] 成像获取单元5,用于根据数据残差,采用预先构建的共轭梯度,以迭代计算方式获取成像数据;
[0099] 成像输出单元6,用于输出最小二乘逆时偏移成像剖面。
[0100] 本发明实施例提供的基于弹性波最小二乘逆时偏移成像的方法,基于矢量弹性波方程,直接对多分量地震资料成像,使纵波和横波的反射同相轴准确归位,以提高地下成像精度,鲁棒性好。
[0101] 优选地,梯度构建单元2,具体用于,
[0102] 根据预先构建阶速度-应力方程组获取正传震源波场;
[0103] 采用预设的存储规则,保存正传震源波场在有效边界内的速度分量波场;
[0104] 根据速度分量波场,以及根据预先构建的伴随反传方程组获取反传震源波场;
[0105] 根据伴随反传方程组,获取反传检波点波场;
[0106] 根据反传震源波场与反传检波点波场构建弹性波最小二乘逆时偏移的梯度。
[0107] 本实施例中,采用如图3所示的存储规则对正传震源波场进行存储。此外,图4是本实施例提供的正传的震源波场的垂直速度分量波场,以及图5是本实施例提供的重建的震源波场的垂直速度分量波场,根据图4和图5,可以获知,本实施例提供的系统能够很好地对震源波场进行重建与恢复,此外,更加具体地,针对图5中存在的重建误差,进行了更加具体的分析,且图6是图5中对应的x=0.8km处的重建误差分析图,且图6中的第一列图为正传的垂直速度分量波场的单道,第二列图为重建的垂直速度分量波场的单道,第三列图为二者的绝对误差,从图6中可见,误差极小,也就是说,本实施例提供的系统能够很好地对震源波场进行重建与恢复。
[0108] 此外,需要说明的是,图7是层状介质图,图7中的第一列的上下两图分别为真实的纵波速度相对变化量与真实的横波速度相对变化量,第二列的上下两图分别为采用本实施例的系统成像获取的纵波速度的相对变化量(即, )和横波速度的相对变化量(即, )。图8是地下形态图,图8中的第一列的上下两图分别为真实的纵波速度相对变化量与真实的横波速度相对变化量,第二列的上下两图分别为采用本实施例的系统成像获取的纵波速度的相对变化量(即, )和横波速度的相对变化量(即, )。图7和图8具体的阐述了本实施例提供的系统的有效性与实用性。
[0109] 尽管本发明已进行了一定程度的描述,明显地,在不脱离本发明的精神和范围的条件下,可进行各个条件的适当变化。可以理解,本发明不限于所述实施方案,而归于
权利要求的范围,其包括所述每个因素的等同替换。