• 浪潮孙丕恕从信息时代到智能时代 人工智能价值将爆发式释放 2019-12-21
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
    • / 11
    • 下载费用:30 金币  

    重庆时时彩500如何刷水: 多波联合振幅随入射角变化反演方法.pdf

    关 键 词:
    联合 振幅 入射角 变化 反演 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201410250358.X

    申请日:

    2014.06.06

    公开号:

    CN103984018A

    公开日:

    2014.08.13

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||著录事项变更IPC(主分类):G01V 1/28变更事项:发明人变更前:石瑛 芦俊变更后:芦俊 石瑛|||实质审查的生效IPC(主分类):G01V 1/28申请日:20140606|||公开
    IPC分类号: G01V1/28 主分类号: G01V1/28
    申请人: 中国地质大学(北京)
    发明人: 石瑛; 芦俊
    地址: 100083 北京市海淀区学院路29号
    优先权:
    专利代理机构: 北京国昊天诚知识产权代理有限公司 11315 代理人: 许志勇;王中
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410250358.X

    授权公告号:

    |||||||||

    法律状态公告日:

    2017.01.11|||2016.06.22|||2014.09.10|||2014.08.13

    法律状态类型:

    授权|||著录事项变更|||实质审查的生效|||公开

    摘要

    本发明适用于地震波技术领域,提供了一种多波联合振幅随入射角变化反演方法,适于地震波的多波叠前联合反演。所述方法包括:接收PP波与PS波的转换波道集,并提取PP波及压缩后的PS波的不同入射角的子波;根据不同入射角的子波,建立第一与第二矩阵;根据薄互层初始模型,建立第一与第二雅各比矩阵及第一与第二反射系数矩阵;根据所述PP波与PS波的转换波道集,建立第三与第四矩阵,并计算第一与第二目标函数;根据第一、第二、第三与第四矩阵、第一与第二雅各比矩阵、第一与第二反射系数矩阵,得到修正值,并修正薄互层初始模型,且调整所述目标函数。通过应用该方法,解决了现有技术存在的多波叠前反演方法会存在较大的误差的问题。

    权利要求书

    权利要求书
    1.  一种多波联合振幅随入射角变化反演方法,适于地震波的多波叠前联合反演,其特征在于,所述目标函数的计算方法包括:
    1)接收反射纵波PP波的转换波道集与反射转换横波PS波的转换波道集,并提取所述PP波的不同入射角的子波以及压缩后的所述PS波的不同入射角的子波,其中压缩后的所述PS波为压缩所述PS波至PP波的预设时;
    2)根据所述不同入射角的子波,建立第一矩阵与第二矩阵;
    根据薄互层初始模型,建立第一雅各比矩阵与第二雅各比矩阵及对应所述薄互层初始模型的第一反射系数矩阵与第二反射系数矩阵;
    3)根据所述PP波的转换波道集与所述PS波的转换波道集,建立第三矩阵与第四矩阵,并计算第一目标函数与第二目标函数;
    4)根据所述第一矩阵、所述第二矩阵、所述第三矩阵、所述第四矩阵、所述第一雅各比矩阵、所述第二雅各比矩阵、所述第一反射系数矩阵与所述第二反射系数矩阵,得到修正值,并修正所述薄互层初始模型,且调整所述第一目标函数与所述第二目标函数。

    2.  根据权利要求1所述的目标函数的计算方法,其特征在于,所述第一矩阵、所述第二矩阵、所述第三矩阵、所述第四矩阵、所述第一雅各比矩阵、所述第二雅各比矩阵、所述第一反射系数矩阵、所述第二反射系数矩阵与所述修正值组合成如下公式:
    [(Wpp*Gpp)T(Wpp*Gpp)+(Wps*Gps)T(Wps*Gps)]△V,=(Wpp*Gpp)T(Dpp-Wpp*Rpp0)+(Wps*Gps)T(Dps-Wps*Rps0)
    其中,Wpp为所述第一矩阵,Wps为所述第二矩阵,Gpp为所述第一雅各比矩阵,Gps为所述第二雅各比矩阵,Rpp0为所述第一反射系数矩阵,Rps0为所述第二反射系数矩阵,(Dpp-Wpp*Rpp0)为所述第三矩阵,(Dps-Wps*Rps0)为所述第四矩阵,△V为所述修正值。

    3.  根据权利要求2所述的目标函数的计算方法,其特征在于,所述公 式更加入阻尼矩阵;
    其中所述阻尼矩阵为λI,其中I为单位向量矩阵,λ取极小值。

    4.  根据权利要求2所述的目标函数的计算方法,其特征在于,所述第一目标函数与第二目标函数分别为如下公式:
    Q(V)=‖Spp-Dpp‖2+‖Sps-Dps‖2,
    T=λ1·(|Spp-Dpp|+|Sps-Dps|)+λ2·(|vpi-vpw|+|vsi-vsw|+|ρi-ρw|),]]>
    其中,Dpp、Dps为实际地震记录,Spp=Wpp*Rpp和Sps=Wps*Rps为模型响应,i为反演得到的参数,w为测井模型外推得到的参数,λ1、λ2分别对应反演中地震道与测井模型约束作用的权重。

    5.  根据权利要求1所述的目标函数的计算方法,其特征在于,所述调整所述第一目标函数与所述第二目标函数包括迭代所述步骤1)、2)、3)、4)至所述第一目标函数与所述第二目标函数的值最小为止。

    6.  根据权利要求1所述的目标函数的计算方法,其特征在于,所述薄互层初始模型、所述PP波的转换波道集与PS波的转换波道集为标准地震数据记录格式。

    7.  根据权利要求1所述的目标函数的计算方法,其特征在于,所述不同入射角的子波为文本格式。

    说明书

    说明书多波联合振幅随入射角变化反演方法
    技术领域
    本发明涉及地震波技术领域,尤其涉及一种多波联合振幅随入射角变化反演方法。
    背景技术
    目前从许多的文献来看,多波振幅随入射角变化(Amplitude Various Angle,AVA)叠前联合反演方法一般都是针对油气田的单界面反演,采用了阿基-理查德(Aki-Richards)近似公式以及与之相关的改造公式来模拟单界面的反射系数。这类反演方法的特点是假设界面上下的介质的波阻抗是弱反差的。
    但对于薄互层来说,薄互层的波阻抗一般与围岩相差太大,并且薄互层是薄层,对于地震波频带来说,薄互层顶、底板的反射是无法区分的,这也使得上述多波叠前AVA联合反演方法用于薄互层会存在较大的误差。
    发明内容
    本发明的主要目的在于提供一种多波联合振幅随入射角变化反演方法,以解决现有技术存在的多波叠前AVA联合反演方法用于薄互层会存在较大的误差的问题。
    为解决上述问题,本发明实施例提供一种目标函数的计算方法,适于地震波的多波叠前联合反演。所述方法包括:1)接收PP波的转换波道集与PS波的转换波道集,并提取所述PP波的不同入射角的子波以及压缩后的所述PS波的不同入射角的子波,其中压缩后的所述PS波为压缩所述PS波至PP波的预设时;2)根据所述不同入射角的子波,建立第一矩阵与第二矩阵;根据一薄互层初始模型,建立第一雅各比矩阵与第二雅各比矩阵及对应所述薄互层初始模型的第一反射系数矩阵与第二反射系数矩阵;3)根据所述PP波 的转换波道集与PS波的转换波道集,建立第三矩阵与第四矩阵,并计算第一目标函数与第二目标函数;4)根据所述第一矩阵、所述第二矩阵、所述第三矩阵、所述第四矩阵、所述第一雅各比矩阵、所述第二雅各比矩阵、所述第一反射系数矩阵与所述第二反射系数矩阵,得到修正值,并修正所述薄互层初始模型,且调整所述第一目标函数与所述第二目标函数。
    其中,所述第一矩阵、所述第二矩阵、所述第三矩阵、所述第四矩阵、所述第一雅各比矩阵、所述第二雅各比矩阵、所述第一反射系数矩阵、所述第二反射系数矩阵与所述修正值组合成如下公式:[(Wpp*Gpp)T(Wpp*Gpp)+(Wps*Gps)T(Wps*Gps)]ΔV=(Wpp*Gpp)T(Dpp-Wpp*Rpp0)+(Wps*Gps)T(Dps-Wps*Rps0)]]>其中,Wpp为所述第一矩阵,Wps为所述第二矩阵,Gpp为所述第三矩阵,Gps为所述第四矩阵,Rpp0为所述第一雅各比矩阵,Rps0为所述第二雅各比矩阵,(Dpp-Wpp*Rpp0)为所述第一反射系数矩阵,Rps0为所述第二反射系数矩阵,△V为所述修正值。
    其中,所述公式更加入阻尼矩阵;其中所述阻尼矩阵为λI,其中I为单位向量矩阵,λ取极小值。
    其中,所述第一目标函数与第二目标函数分别为如下公式:
    Q(V)=‖Spp-Dpp‖2+‖Sps-Dps‖2,
    T=λ1·(|Spp-Dpp|+|Sps-Dps|)+λ2·(|vpi-vpw|+|vsi-vsw|+|ρi-ρw|),]]>
    其中,Dpp、Dps为实际地震记录,Spp=Wpp*Rpp和Sps=Wps*Rps为模型响应,i为反演得到的参数,w为测井模型外推得到的参数,λ1、λ2分别对应反演中地震道与测井模型约束作用的权重。
    其中,所述调整所述第一目标函数与所述第二目标函数包括迭代所述步骤1)、2)、3)、4)至所述第一目标函数与所述第二目标函数的值最小为止。
    其中,所述薄互层初始模型、所述PP波的转换波道集与PS波的转换波 道集为标准地震数据记录格式。
    其中,所述不同入射角的子波为文本格式。
    根据本发明的技术方案,采用多波地震数据以及测井共同约束反演的结果,基于广义反演理论同时反演薄互层的纵波速度、横波速度以及密度。其中,PP波与PS波的反射系数采用策普里兹方程组精确求解。由于反射系数的求解过程无任何近似,所以本发明的目标函数的计算方法即使在无井约束的情况下也能对薄互层的储层参数精确反演。
    附图说明
    此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
    图1是根据本发明实施例的模拟得到的PP波的AVA道集的示意图;
    图2是根据本发明实施例的模拟得到的PS波的AVA道集的示意图;
    图3是根据本发明实施例的多波联合AVA反演方法的流程图。
    具体实施方式
    本发明的主要思想在于,采用多波地震数据以及测井共同约束反演的结果,基于广义反演理论同时反演薄互层的纵波速度、横波速度以及密度。其中,PP波与PS波的反射系数采用策普里兹(Zoeppritz)方程组精确求解。由于反射系数的求解过程无任何近似,所以本发明的目标函数的计算方法即使在无井约束的情况下也能对薄互层的储层参数精确反演。
    为使本发明的目的、技术方案和优点更加清楚,以下结合附图及具体实施例,对本发明作进一步地详细说明。
    首先,在各向同性介质中,当一个平面纵波倾斜入射到两种介质分界面时,会产生四种波,即反射纵波、反射横波、透射纵波和透射横波,且满足斯涅尔定律(Snell's Law),如公式(1.1)所示:
    Vp1sinα=Vs1sinβ=Vp2sinα=Vs2sinβ.---(1.1)]]>
    根据边界条件解波动方程并引入纵波反射系数RPP、转换横波反射系数RPS、纵波透射系数TPP和转换横波透射TPS,可以得出四个波的位移振幅应当满足的策普里兹(Zoeppritz)方程组,如公式(1.2)所示:
    sinαRPP+cosβRPS-sinαTPP+cosβTPS=-sinα,cosαRPP+sinβRPS-cosαTPP+sinβTPS=cosα,cos2βRPP-VS1VP1sin2βRPS-ρ2VP2ρ1VS2cos2βTPPρ2VP2ρ1VS2sin2-βTPS=-cos2βVS22VP1sin2αRPP+VS1cos2βRPS+ρ2VS22ρ1VP2sin2αTPP-ρ2VS2ρ1cosβ2TPS=VS12VP1sin2α---(1.2)]]>
    在公式(1.2)中,VP1和VS1分别为上层介质纵、横波速度;VP2和VS2为下层介质的纵、横波速度;ρ1和ρ2为上、下层介质的密度;α和α′为纵波入射角和透射角;β和β′为横波反射角和透射角。对上述公式(1.2)求解则可得到P-P波、P-SV波的反射系数RPP、RPS。
    根据公式(1.2)可以分析各向同性介质条件下薄互层的多波AVA地震响应特征。如表1所示为薄互层模型参数表,从浅至深的地层分别为:第四系地层、砂岩(Sand)、薄互层(Coal)、砂岩、砂岩,其中薄互层厚度(Thickness)为3米。
    Vp(m/s)Vs(m/s)ρ(g/cc)层(m)第四系地层28009332.28100砂岩320016002.35100薄互层200013332.03砂岩330016502.3580砂岩340017002.5270
    表1薄互层模型参数表
    图1是根据本发明实施例的模拟得到的PP波的AVA道集的示意图。图2是根据本发明实施例的模拟得到的PS波的AVA道集的示意图。如图1及2所示,其中PS波道集已经被压缩至PP波的T0时间,子波为30Hz的雷克子波。
    从图1及2可以看出:由于薄互层为薄层,薄互层的PP波与PS波地震反射无法区分顶底板的位置,顶、底板反射形成了复合波形态;PP波与PS波的AVA特征有较大区别,对于第四系地层与砂岩界面来说,PP波随着入射角增大,反射能量减弱,但PS波反之;但对于薄互层反射来说,PP波与PS波的反射能量均随着入射角增大而增大;PP波与PS波的极性大多数相反。
    另外,反演采用双目标函数分别为公式(1.3)及(1.4),如下所示:
    Q(V)=‖Spp-Dpp‖2+‖Sps-Dps‖2→minimum,                  (1.3)
    T=λ1·(|Spp-Dpp|+|Sps-Dps|)+λ2·(|vpi-vpw|+|vsi-vsw|+|ρi-ρw|)→minimum,---(1.4)]]>
    在公式(2.3)中,Dpp、Dps为实际地震记录,Spp=Wpp*Rpp和Sps=Wps*Rps为模型响应;在公式(2.4)中,下标“i”表示反演得到的参数,下标“w”表示测井模型外推得到的参数;λ1、λ2分别对应反演中地震道与测井模型约束作用的权重,λ1+λ2=1,当各取50%时,测井数据与地震数据在反演过程中的约束权重各占一半;当λ1=1时表明地震数据完全约束反演结果,测井数据不起约束作用;当λ2=1时,表明测井数据完全约束反演结果,地震数据不起约束作用,相当于将测井模型进行外推。若某一数据的可靠性较低,可以相应降低其约束的权重。
    在公式(1.3)中,对纵波和转换波进行泰勒展开成公式(1.5),如下所示:
    R(V0+△V0)=R(V0)+G(V0)△V0,                   (1.5)
    其中,V0为纵波速度、横波速度和密度的初始猜测值,可表示为向量,如公式(1.6)所示:
    (vp vs ρ),                                 (1.6)
    △V0为修正量,如公式(1.7)所示:
    (△vp △vs △ρ),                          (1.7)
    G(V0)称为f(V0)在V0处的梯度,即雅各比矩阵。
    反射系数线性近似公式代入反演目标(公式(1.3)),可得公式(1.8),如下所示:
    Q(V)=‖Wpp*(Rpp0+Gpp△V)-Dpp‖2+‖Wps*(Rps0+Gps△V)-Dps‖2=[Wpp*(Rpp0+Gpp△V)-Dpp]T[Wpp*(Rpp0+Gpp△V)-Dpp],    (1.8)+[Wps*(Rps0+Gps△V)-Dps]T[Wps*(Rps0+Gps△V)-Dps]
    令可得公式(1.9),如下所示:
    ∂Q(V)/∂ΔV=(Wpp*Rpp0-Dpp)T(Wpp*Gpp)+(Wpp*Gpp)T(Wpp*Rpp0-Dpp)+2(Wpp*Gpp)T(Wpp*Gpp)ΔV+(Wps*Rps0-Dps)T(Wps*Gps)+(Wps*Gps)T(Wps*Rps0-Dps)+2(Wps*Gps)T(Wps*Gps)ΔV=0,---(1.9)]]>
    所以根据公式(1.9),可以得出公式(1.10),如下所示:
    [(Wpp*Gpp)T(Wpp*Gpp)+(Wps*Gps)T(Wps*Gps)]△V,   (1.10)=(Wpp*Gpp)T(Dpp-Wpp*Rpp0)+(Wps*Gps)T(Dps-Wps*Rps0)
    其中:

    以上,大略说明了本发明之实施例所需要运用到的相关公式,以下将提 供对应的实施例来进行说明。根据本发明的实施例,提供了一种目标函数的计算方法,其适于地震波的多波叠前联合反演。
    图3是根据本发明实施例的多波联合AVA反演方法的流程图,如图3所示,该方法包括:
    步骤S302,接收PP波的转换波道集与PS波的转换波道集,并提取所述PP波的不同入射角的子波以及压缩后的所述PS波的不同入射角的子波,其中压缩后的所述PS波为压缩所述PS波至PP波的预设时。其中,所述PP波的预设时例如为PP波的T0时。
    步骤S304,根据所述不同入射角的子波,建立第一矩阵与第二矩阵。进一步来说,例如建立公式(1.10)中的Wpp与Wps,其中Wpp为所述第一矩阵,Wps为所述第二矩阵。并且,所述不同入射角的子波例如为文本格式(TXT)。
    步骤S306,根据薄互层初始模型,建立第一雅各比矩阵与第二雅各比矩阵及对应所述薄互层初始模型的第一反射系数矩阵与第二反射系数矩阵。进一步来说,例如建立公式(1.10)中的Gpp、Gps、Rpp0与Rps0,其中Gpp为所述第一雅各比矩阵,Gps为所述第二雅各比矩阵,Rpp0为所述第一反射系数矩阵,Rps0为所述第二反射系数矩阵。并且,所述薄互层初始模型例如为标准地震数据记录格式(SEGY)。
    步骤S308,根据所述PP波的转换波道集与所述PS波的转换波道集,建立第三矩阵与第四矩阵,并计算第一目标函数与第二目标函数。进一步来说,例如建立公式(1.10)中的(Dpp-Wpp*Rpp0)与(Dps-Wps*Rps0),其中(Dpp-Wpp*Rpp0)为所述第三矩阵,(Dps-Wps*Rps0)为所述第四矩阵,而第一目标函数例如为公式(1.3),第二目标函数例如为公式(1.4)。并且,所述PP波的转换波道集与所述PS波的转换波道集例如也为标准地震数据记录格式(SEGY)。
    步骤S310,根据所述第一矩阵、所述第二矩阵、所述第三矩阵、所述第四矩阵、所述第一雅各比矩阵、所述第二雅各比矩阵、所述第一反射系数矩阵与所述第二反射系数矩阵,得到修正值,并修正所述薄互层初始模型,且调整所述第一目标函数与所述第二目标函数。其中,例如根据公式(1.10)求解,可得到△V,其中△V例如为修正值。
    另外,[(Wpp*Gpp)T(Wpp*Gpp)+(Wps*Gps)T(Wps*Gps)]在解矩阵方程需要加入阻尼矩阵,以克服奇异,其中阻尼矩阵为λI,I为单位向量矩阵,λ取极小值,例如0.000001等。由于奇异出现的情况不可预料,有时λ取值很难决断。一般来说,若λ取值误差较大,一般不会影响反演结果,但是个别情况会导致整个反演结果发生整体漂移。这个漂移量可以用井数据来校正。
    综上所述,根据本发明的技术方案,采用多波地震数据以及测井共同约束反演的结果,基于广义反演理论同时反演薄互层的纵波速度、横波速度以及密度。其中,PP波与PS波的反射系数采用策普里兹(Zoeppritz)方程组精确求解。由于反射系数的求解过程无任何近似,所以本发明的目标函数的计算方法即使在无井约束的情况下也能对薄互层的储层参数精确反演。
    以上所述仅为本发明的实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内?!  ∧谌堇醋宰ɡ鴚ww.www.4mum.com.cn转载请标明出处

    关于本文
    本文标题:多波联合振幅随入射角变化反演方法.pdf
    链接地址://www.4mum.com.cn/p-6140458.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 浪潮孙丕恕从信息时代到智能时代 人工智能价值将爆发式释放 2019-12-21
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 股票推荐每日一股9月6日 国彩合买靠谱吗 1000炮打鱼机李逵劈鱼 快乐十分彩票投注技巧 海南七星彩规律图长条2079 2012年奥运会足球直播 广东好彩1开奖时间 江西多乐彩一定牛 竞彩胜分差怎么算的 神手麻将下载 上海快三实时开奖视频 体球比分直播网 大乐购群 辽宁35选7开奖结果2元 黑龙江十一选五任三遗漏 黑龙江时时彩