• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 16
    • 下载费用:30 金币  

    重庆时时彩庄闲怎么看: 一种宽频地震多域联合AVO反演方法.pdf

    关 键 词:
    一种 宽频 地震 联合 AVO 反演 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201610986167.9

    申请日:

    2016.11.09

    公开号:

    CN106772615A

    公开日:

    2017.05.31

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01V 1/48申请日:20161109|||公开
    IPC分类号: G01V1/48; G01V1/50 主分类号: G01V1/48
    申请人: 中国石油大学(华东)
    发明人: 宗兆云; 印兴耀; 吴国忱; 李坤
    地址: 266580 山东省青岛市黄岛经济开发区长江西路66号
    优先权:
    专利代理机构: 北京汇泽知识产权代理有限公司 11228 代理人: 胡小龙
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610986167.9

    授权公告号:

    ||||||

    法律状态公告日:

    2018.11.13|||2017.06.23|||2017.05.31

    法律状态类型:

    授权|||实质审查的生效|||公开

    摘要

    本发明涉及地震勘探领域,具体公开了一种宽频地震多域联合AVO反演方法。该方法包括两个步骤:第一步,拉普拉斯??傅里叶域AVO反演进行弹性参数低频模型预测;第二步,时间域AVO反演进行弹性参数全频带信息预测:利用第一步得到的弹性参数低频模型作为时间域AVO反演的初始模型,实现时间域弹性参数全频带反演预测。本发明利用宽频地震资料,通过在拉普拉斯??傅里叶域反演实现弹性参数低频模型预测,可以克服常规地震AVO反演借助测井信息建模在复杂地质条件下难以得到合理结果的不足,实现模型参数低频模型合理稳定估计,并利用该低频模型作为反演初始模型,在时间域实现弹性参数全频带反演预测。

    权利要求书

    1.一种宽频地震多域联合AVO反演方法,其特征在于,其包括两个步骤:
    第一步,拉普拉斯-傅里叶域AVO反演进行弹性参数低频模型预测:
    地震信号s(t)的拉普拉斯-傅里叶域变换S(ξ)为:
    <mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>&xi;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>S</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>-</mo> <mi>i</mi> <mi>&sigma;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <mrow> <mo>+</mo> <mi>&infin;</mi> </mrow> </msubsup> <mi>s</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mo>&lsqb;</mo> <mo>-</mo> <mi>i</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>-</mo> <mi>i</mi> <mi>&sigma;</mi> <mo>)</mo> </mrow> <mi>t</mi> <mo>&rsqb;</mo> <mi>d</mi> <mi>t</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
    方程(1)中,ξ=ω-iσ为复频率,ω为频率,σ为拉普拉斯域阻尼因子,i为虚部,t为时
    间,
    对于叠前地震资料,s(t)表示为:
    <mrow> <mi>s</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <mrow> <mo>+</mo> <mi>&infin;</mi> </mrow> </msubsup> <mi>w</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&tau;</mi> <mo>)</mo> </mrow> <mi>r</mi> <mrow> <mo>(</mo> <mi>&tau;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&tau;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
    方程(2)中,w(t)为地震子波,r(t)为地震反射系数表示为:
    <mrow> <mi>r</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msup> <mi>sec</mi> <mn>2</mn> </msup> <mi>&theta;</mi> </mrow> <mn>2</mn> </mfrac> <mfrac> <mrow> <msub> <mi>&Delta;V</mi> <mi>p</mi> </msub> </mrow> <msub> <mi>V</mi> <mi>p</mi> </msub> </mfrac> <mo>-</mo> <mn>4</mn> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>V</mi> <mi>s</mi> </msub> <msub> <mi>V</mi> <mi>p</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <mi>sin</mi> <mn>2</mn> </msup> <mi>&theta;</mi> <mfrac> <mrow> <msub> <mi>&Delta;V</mi> <mi>s</mi> </msub> </mrow> <msub> <mi>V</mi> <mi>s</mi> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>4</mn> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>V</mi> <mi>s</mi> </msub> <msub> <mi>V</mi> <mi>p</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <mi>sin</mi> <mn>2</mn> </msup> <mi>&theta;</mi> <mo>)</mo> </mrow> <mfrac> <mrow> <mi>&Delta;</mi> <mi>&rho;</mi> </mrow> <mi>&rho;</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
    方程(3)中,θ为入射角,和分别为纵波速度、横波速度和密度反射系数,
    将方程(2)和方程(3)带入方程(1)并进行拉普拉斯-傅里叶域变换得:
    S'=F×C'σ×R (4)
    方程(4)中,
    <mrow> <msup> <mi>S</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>S</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>S</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>S</mi> <mi>k</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mrow> <mi>m</mi> <mi>k</mi> <mo>*</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> </mrow>

    <mrow> <msubsup> <mi>C</mi> <mi>&sigma;</mi> <mo>&prime;</mo> </msubsup> <mo>=</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>C</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>C</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>C</mi> <mi>k</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mrow> <mi>n</mi> <mi>k</mi> <mo>*</mo> <mi>n</mi> </mrow> </msub> <mo>,</mo> </mrow>
    <mrow> <mi>R</mi> <mo>=</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>r</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>r</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>r</mi> <mi>n</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mrow> <mi>n</mi> <mo>*</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> </mrow>
    其中,



    且m,n和k分别是子波矩阵、反射系数矩阵和地震数据的行数,
    利用方程(4)作为反演正演算子,建立拉普拉斯-傅里叶域AVO反演目标函数K为:
    <mrow> <mi>K</mi> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>S</mi> <mo>&prime;</mo> </msup> <mo>-</mo> <msubsup> <mi>FC</mi> <mi>&sigma;</mi> <mo>&prime;</mo> </msubsup> <mi>R</mi> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msup> <mi>S</mi> <mo>&prime;</mo> </msup> <mo>-</mo> <msubsup> <mi>FC</mi> <mi>&sigma;</mi> <mo>&prime;</mo> </msubsup> <mi>R</mi> <mo>)</mo> </mrow> <mo>+</mo> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mi>l</mi> <mi>n</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msubsup> <mi>r</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>/</mo> <msubsup> <mi>&sigma;</mi> <mi>m</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msub> <mi>&lambda;</mi> <mi>l</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>-</mo> <msub> <mi>&beta;R</mi> <mrow> <mi>r</mi> <mi>l</mi> </mrow> </msub> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>-</mo> <msub> <mi>&beta;R</mi> <mrow> <mi>r</mi> <mi>l</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
    其中,σm和σn分别为拉普拉斯-傅里叶域地震噪音和弹性参数协方差,
    Rl=[Vp,Vs,ρ],λl为不同弹性参数权重因子,β为α

    采用迭代加权最小二乘算法实现对方程(8)优化求解,进而得到弹性参数低频模型;
    第二步,时间域AVO反演进行弹性参数全频带信息预测:
    利用第一步得到的弹性参数低频模型作为时间域AVO反演的初始模型,实现时间域弹
    性参数全频带反演预测,具体包括如下步骤:
    时间域AVO反演所需的正演算子为:
    d=Gm (9)
    方程(9)中,d,G和m分别为时间域叠前地震数据,子波矩阵和弹性参数矩阵,
    根据方程(9),建立时间域AVO反演目标函数为:
    <mrow> <mtable> <mtr> <mtd> <mrow> <mi>F</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <mi>G</mi> <mi>m</mi> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <mi>G</mi> <mi>m</mi> <mo>)</mo> </mrow> <mo>+</mo> <mn>2</mn> <msup> <msub> <mi>&sigma;</mi> <mi>n</mi> </msub> <mn>2</mn> </msup> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <mi>ln</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msubsup> <mi>m</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>/</mo> <msubsup> <mi>&sigma;</mi> <msub> <mi>m</mi> <mi>i</mi> </msub> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msub> <mi>&lambda;</mi> <mi>M</mi> </msub> <msup> <mrow> <mo>(</mo> <msub> <mi>&eta;</mi> <mi>M</mi> </msub> <mo>-</mo> <msub> <mi>PM</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>&eta;</mi> <mi>M</mi> </msub> <mo>-</mo> <msub> <mi>PM</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&lambda;</mi> <mi>&mu;</mi> </msub> <msup> <mrow> <mo>(</mo> <msub> <mi>&eta;</mi> <mi>&mu;</mi> </msub> <mo>-</mo> <msub> <mi>P&mu;</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>&eta;</mi> <mi>&mu;</mi> </msub> <mo>-</mo> <msub> <mi>P&mu;</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&lambda;</mi> <mi>D</mi> </msub> <msup> <mrow> <mo>(</mo> <msub> <mi>&eta;</mi> <mi>D</mi> </msub> <mo>-</mo> <msub> <mi>PD</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>&eta;</mi> <mi>D</mi> </msub> <mo>-</mo> <msub> <mi>PD</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
    方程(10)中,Mr,μr和Dr分别为时间域纵波速度反射系数矩阵、横波速度反射系数矩阵、
    密度反射系数矩阵,
    λM,λμ,λD分别为时间域纵波速度、横波速度和密度权重系数,
    P为
    <mrow> <msub> <mi>&eta;</mi> <mi>M</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>l</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>/</mo> <msub> <mi>M</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
    ηM=1/2*ln(μ/μ0),
    ηM=1/2*ln(D/D0),
    其中,M,μ和D分别为时间域纵波速度、横波速度和密度矩阵,M0,μ0和D0分别为时间域纵
    横波速度和密度初始值,
    采用迭代加权最小二乘算法实现对方程(10)优化求解,得到全频带弹性参数。

    说明书

    一种宽频地震多域联合AVO反演方法

    技术领域

    本发明涉及地震勘探领域,特别涉及一种宽频地震多域联合AVO反演方法。

    背景技术

    AVO反演是获取地下介质弹性参数的重要手段。常规地震AVO反演主要利用地震资
    料随偏移距的变化信息,建立地震反射系数与地下介质弹性参数及偏移距或入射角之间的
    定量关系,通过选择合适的反演优化算法,实现弹性介质参数反演预测。常规AVO反演主要
    包括子波提取、低频模型构建、反演目标函弹性参数建立、优化反演几个步骤。其中,低频模
    型构建是否合理是评价反演合理性的重要因素。

    常规地震资料频带一般在8-100Hz范围内,欠缺有效的低频信息,因此需要借助其
    他地球物理物理资料(如测井资料等)的低频信息来建立低频模型约束AVO反演,以便获取
    全频带参数信息。常规地震勘探低频模型主要利用目标区已知测井资料,在地质或层位约
    束下通过插值或外推建立,这种方法在地质条件较为简单,研究目标区横向物性变化较小
    的情况下可以获取较好的建模结果,但是在地质条件复杂,横向变化较大的目标区很难建
    立合理的初始模型。

    发明内容

    本发明要解决的技术问题是:利用宽频地震资料,通过在拉普拉斯-傅里叶域反演
    实现弹性参数低频模型预测,并利用该低频模型作为反演初始模型,在时间域实现弹性参
    数全频带反演预测。

    为了解决上述技术问题,本发明所采用的技术方案是:一种宽频地震多域联合AVO
    反演方法,其包括两个步骤:

    第一步,拉普拉斯-傅里叶域AVO反演进行弹性参数低频模型预测:

    地震信号s(t)的拉普拉斯-傅里叶域变换S(ξ)为:


    方程(1)中,ξ=ω-iσ为复频率,ω为频率,σ为拉普拉斯域阻尼因子,i为虚部,t为
    时间,

    对于叠前地震资料,s(t)表示为:


    方程(2)中,w(t)为地震子波,r(t)为地震反射系数表示为:


    方程(3)中,θ为入射角,和分别为纵波速度、横波速度和密度反射
    系数,

    将方程(2)和方程(3)带入方程(1)并进行拉普拉斯-傅里叶域变换得:

    S'=F×C'σ×R (4)

    方程(4)中,





    其中,




    且m,n和k分别是子波矩阵、反射系数矩阵和地震数据的行数,

    利用方程(4)作为反演正演算子,建立拉普拉斯-傅里叶域AVO反演目标函数K为:


    其中,σm和σn分别为拉普拉斯-傅里叶域地震噪音和弹性参数协方差,

    Rl=[Vp,Vs,ρ],λl为不同弹性参数权重因子,β为
    α为

    采用迭代加权最小二乘算法实现对方程(8)优化求解,进而得到弹性参数低频模
    型;

    第二步,时间域AVO反演进行弹性参数全频带信息预测:

    利用第一步得到的弹性参数低频模型作为时间域AVO反演的初始模型,实现时间
    域弹性参数全频带反演预测,具体包括如下步骤:

    时间域AVO反演所需的正演算子为:

    d=Gm (9)

    方程(9)中,d,G和m分别为时间域叠前地震数据,子波矩阵和弹性参数矩阵,

    根据方程(9),建立时间域AVO反演目标函数为:


    方程(10)中,Mr,μr和Dr分别为时间域纵波速度反射系数矩阵、横波速度反射系数
    矩阵、密度反射系数矩阵,

    λM,λμ,λD分别为时间域纵波速度、横波速度和密度权重系数,

    P为


    ηM=1/2*ln(μ/μ0),

    ηM=1/2*ln(D/D0),

    其中,M,μ和D分别为时间域纵波速度、横波速度和密度矩阵,M0,μ0和D0分别为时间
    域纵横波速度和密度初始值,

    采用迭代加权最小二乘算法实现对方程(10)优化求解,得到全频带弹性参数。

    本发明技术方案带来的有益效果是:克服常规地震AVO反演借助测井信息建模在
    复杂地质条件下难以得到合理结果的不足,实现模型参数低频模型合理稳定估计,并利用
    该低频模型作为反演初始模型,在时间域实现弹性参数全频带反演预测。

    附图说明

    图1本发明一种宽频地震多域联合AVO反演方法主要步骤流程图。

    图2A为小角度地震剖面示意图。

    图2B为中角度地震剖面示意图。

    图2C为大角度地震剖面示意图。

    图3A为拉普拉斯-傅里叶域AVO纵波速度反演结果示意图。

    图3B为拉普拉斯-傅里叶域AVO横波速度反演结果示意图。

    图3C为拉普拉斯-傅里叶域AVO密度反演结果示意图。

    图4A为时间域AVO纵波速度反演结果示意图。

    图4B为时间域AVO横波速度反演结果示意图。

    图4C为时间域AVO密度反演结果示意图。

    图5A为井位置处纵波速度对比结果示意图。

    图5B为井位置处横波速度对比结果示意图。

    图5C为井位置处密度对比结果示意图。

    具体实施方式

    下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以
    更好的理解本发明并能予以实施,但所举实施例不作为对本发明的限定。

    宽频地震资料有效频带低频端可低至2Hz,给AVO反演低频模型建立带了新的机
    遇,即可以充分利用宽频地震有效的低频信息进行低频模型预测。本发明充分利用拉普拉
    斯-傅里叶域反演在挖掘低频信息方面的优势,提出一种宽频地震多域联合AVO反演方法。

    如图1所示,本发明提出的一种宽频地震多域联合AVO反演方法主要包括两个步
    骤:

    第一步,拉普拉斯-傅里叶域AVO反演进行弹性参数低频模型预测:

    地震信号s(t)的拉普拉斯-傅里叶域变换S(ξ)为:


    其中,ξ=ω-iσ为复频率,ω为频率,σ为拉普拉斯域阻尼因子,i为虚部,t为时间;

    对于叠前地震资料,s(t)可以表示为:


    其中,w(t)为地震子波,r(t)为地震反射系数,可以表示为:


    其中,θ为入射角,和分别为纵波速度、横波速度和密度反射系数;

    将方程(2)和(3)带入方程(1)并进行拉普拉斯-傅里叶域变换得:

    S'=F×C'σ×R (4)

    其中,





    其中,




    且m,n和k分别是子波矩阵、反射系数矩阵和地震数据的行数;

    利用方程(4)作为反演正演算子,建立拉普拉斯-傅里叶域AVO反演目标函数K为:


    其中,σm和σn分别为拉普拉斯-傅里叶域地震噪音和弹性参数协方差;

    Rl=[Vp,Vs,ρ],λl为不同弹性参数权重因子,β为
    α为

    采用迭代加权最小二乘算法即可实现对方程(8)优化求解,进而得到弹性参数低
    频模型。

    第二步,时间域AVO反演进行弹性参数全频带信息预测:

    利用第一步得到的弹性参数低频模型作为时间域AVO反演的初始模型,可以实现
    时间域弹性参数全频带反演预测,具体包括如下步骤:

    时间域AVO反演所需的正演算子为:

    d=Gm (9)

    其中,d,G和m分别为时间域叠前地震数据,子波矩阵和弹性参数矩阵。

    根据方程(9),考虑弹性参数先验信息及模型约束可建立时间域AVO反演目标函数
    为:


    其中,Mr,μr和Dr分别为时间域纵波速度反射系数矩阵、横波速度反射系数矩阵、密
    度反射系数矩阵。

    λM,λμ,λD分别为时间域纵波速度、横波速度和密度权重系数,

    P为


    ηM=1/2*ln(μ/μ0),

    ηM=1/2*ln(D/D0),

    其中,M,μ和D分别为时间域纵波速度、横波速度和密度矩阵,M0,μ0和D0分别为时间
    域纵横波速度和密度初始值;

    同样采用迭代加权最小二乘算法即可实现对方程(10)优化求解,即可得到全频带
    弹性参数。

    为验证本发明方法的有效性,采用分角度部分叠加地震资料(如图2A、图2B、图2C
    所示)进行多域联合AVO反演方法,其中,图2A为小角度地震剖面示意图,图2B为中角度地震
    剖面示意图,图2C为大角度地震剖面示意图。

    利用拉普拉斯-傅里叶域AVO反演得到的弹性参数低频模型如图3A、图3B、图3C所
    示,图3A为拉普拉斯-傅里叶域AVO纵波速度反演结果,图3B为拉普拉斯-傅里叶域AVO横波
    速度反演结果,图3C为拉普拉斯-傅里叶域AVO密度反演结果。

    利用该低频模型作为时间域反演初始模型,可以得到时间域反演结果如图4A、图
    4B、图4C所示。图4A为时间域AVO纵波速度反演结果,图4B为时间域AVO横波速度反演结果,
    图4C为时间域AVO密度反演结果。

    抽取反演结果井旁道对比结果如图5A、图5B、图5C所示,图5A为井位置处纵波速度
    对比结果,图5B为井位置处横波速度对比结果,图5C为井位置处密度对比结果,其中黑色细
    线为真实时间域测井曲线,黑色粗线为拉普拉斯-傅里叶域AVO反演结果,虚线为时间域反
    演结果。

    从图5A、图5B、图5C中可知,拉普拉斯-傅里叶域AVO反演得到的弹性参数能够较好
    的反演参数的低频背景,时间域反演结果和实际时间域测井曲线有较高的吻合度,验证了
    本发明方法的有效性和合理性。

    本发明利用宽频地震资料,通过在拉普拉斯-傅里叶域反演实现弹性参数低频模
    型预测,可以克服常规地震AVO反演借助测井信息建模在复杂地质条件下难以得到合理结
    果的不足,实现模型参数低频模型合理稳定估计,并利用该低频模型作为反演初始模型,在
    时间域实现弹性参数全频带反演预测。

    以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的?;し?br />围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明
    的?;し段е?。本发明的?;し段б匀ɡ笫槲?。

    关于本文
    本文标题:一种宽频地震多域联合AVO反演方法.pdf
    链接地址://www.4mum.com.cn/p-6001533.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

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


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 类似彩名堂的计划软件 上海时时开奖信息查询 香港36码特围 刘伯温的六肖期期准 快速时时走势图 七乐彩开奖直播电视台 时时彩 稳赚验证 曾道长六肖期期准14期开奖 快3稳赚投注技巧 分分彩9码平刷稳赚方案 重庆时时前三星走势图 手机重庆时彩计划软件下载 安徽时时计划软件手机版下载手机版下载手机版下载 幸运飞艇稳赚技巧 全天pk10计划两期稳定版 3d福彩五码复式计划大师