• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
    • / 23
    • 下载费用:20 金币  

    重庆时时彩即开奖: 一种反褶积方法及装置.pdf

    摘要
    申请专利号:

    重庆时时彩单双窍门 www.4mum.com.cn CN201410132417.3

    申请日:

    2014.04.03

    公开号:

    CN103954992A

    公开日:

    2014.07.30

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01V 1/28申请日:20140403|||公开
    IPC分类号: G01V1/28 主分类号: G01V1/28
    申请人: 中国石油天然气股份有限公司
    发明人: 王万里; 魏新建; 何欣; 禄娟
    地址: 100007 北京市东城区东直门北大街9号
    优先权:
    专利代理机构: 北京三友知识产权代理有限公司 11127 代理人: 党晓林
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410132417.3

    授权公告号:

    ||||||

    法律状态公告日:

    2017.04.05|||2014.08.27|||2014.07.30

    法律状态类型:

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

    摘要

    本发明提供了一种反褶积方法及装置。该方法包括以下处理步骤:将采集的地震数据处理成叠后地震数据记录;构建包括反射系数序列在内的目标函数;在所述叠后地震数据记录中选取一道所需处理的地震资料,从所述构建的目标函数中计算出所述地震资料中反射系数序列的振幅和位置;读取下一道所需处理的地震资料,计算出所述下一道地震资料反射系数序列的振幅和位置;计算出所述叠后地震数据记录中所有地震道反射系数序列的振幅和位置,输出反射系数序列,完成所述叠后地震数据记录的反褶积处理。本发明提供的反褶积方法能进一步提高地震资料的分辨率。

    权利要求书

    权利要求书
    1.  一种反褶积方法,其特征在于,包括以下处理步骤:
    S1:将所需处理的地震数据处理成叠后地震数据记录;
    S2:构建包括反射系数序列在内的目标函数;
    S3:在所述叠后地震数据记录中选取一道所需处理的地震资料,从所述构建的目标函数中计算出所述地震资料的反射系数序列的振幅和位置;
    S4:读取下一道所需处理的地震资料,计算出所述下一道地震资料反射系数序列的振幅和位置;
    S5:计算出所述叠后地震数据记录中所有地震道反射系数序列的振幅和位置,输出反射系数序列,完成所述叠后地震数据记录的反褶积处理。

    2.  如权利要求1所述的一种反褶积方法,其特征在于,所述构建目标函数包括:
    S201:建立地震数据褶积模型;
    S202:将所述地震数据褶积模型转换成矩阵形式;
    S203:由所述矩阵形式构建目标函数;所述构建的目标函数为:
    Jα=||Aα-S||2
    上式中,Jα为关于反射系数序列α的范数,A为一个N×M的矩阵,其表达式为其中,ωt为地震子波,t=1,2,…N,N为所述叠后地震数据记录的时间窗长度,M为预先设置的所述时间窗长度N内反射系数序列的个数,τM第M个反射系数的时延,s为t时刻所述叠后地震数据记录的值。

    3.  如权利要求2所述的一种反褶积方法,其特征在于,所述构建的目标函数为:
    J=||Aα-S||2+μ||Cα-ε||2=Jα+μJc
    上式中,μ为阻抗权系数,C为左下三角全为1的矩阵,ε为波阻抗。

    4.  如权利要求3所述的一种反褶积方法,其特征在于,所述构建的目标函数为:
    J=Jα+μJc+βJα
    上式中,β为正则化的阻尼系数,其中,β=β0f0,β0为预设值,f0=max{(ATA+μCTCii)},i为矩阵(ATA+μCTC)的下标。

    5.  如权利要求1-4任意一项所述的一种反褶积方法,其特征在于,从所述构建的目标函数中计算出所述地震资料的反射系数序列的振幅包括步骤:
    S310:从所述叠后地震数据记录中估计出地震子波;
    S320:将所述估计出的地震子波代入所述目标函数中,计算出所述反射系数序列的振幅;
    其中,所述估计地震子波的方法包括步骤:
    S311:建立所述叠后地震数据记录的褶积模型;
    S312:将所述建立的褶积模型变换到复赛谱域;
    S313:设计低通滤波器,在复赛谱域将所述建立的褶积模型的地震子波和反射系数序列进行分离,得到复赛谱域的地震子波;
    S314:将所述复赛谱域的地震子波转换到时间域,得到时间域的地震子波。

    6.  如权利要求5所述的一种反褶积方法,其特征在于,所述地震子波的估计方法还包括步骤:
    315:对所述估计出的地震子波相位进行校正。

    7.  如权利要求5所述的一种反褶积方法,其特征在于,采用模拟退火算法计算反射系数序列中每一个反射系数的时间位置,其操作过程包括以下步骤:
    S301:在所述时间窗长度N内随机给定反射系数序列a(j)的初始时间位置,由所述初始时间位置计算所述目标函数Ja当前解,j=1,2…M;
    S302:在所述时间窗长度N内对所述反射系数序列a(j)进行扰动,将所述扰动后的反射系数序列a(j)代入所述目标函数Ja中,计算出所述目标函数Ja的新解;
    S303:将所述Jα的新解与当前解进行比较,根据Metropolis准则判断是否接受所述目标函数Jα的新解;
    如果判断结果为接受所述新解,则接受反射系数α(j)当前的时间位置;此时,所述新解作为目标函数Ja当前解;
    S304:重复S302、S303直到达到收敛条件后,停止迭代,得到所述反射系数α(j)的时间位置。

    8.  如权利要求7所述的一种反褶积方法,其特征在于,所述模拟退火算法还包括在所述S303之后:
    S3031:设置误差项σ,将所述新解与当前解的差值与所述误差项σ比较;如果所述新解与当前解的差值大于所述误差项σ,则不接受所述目标函数Jα的新解。

    9.  如权利要求7所述的一种反褶积方法,其特征在于,在达到收敛条件后,将所述模拟退火算法得出的反射系数序列α(j)的振幅值作为地震资料中反射系数序列的振幅。

    10.  一种反褶积装置,其特征在于,该装置包括:数据预处理???、目标函数建立???、目标函数计算???、资料读取???、反射系数序列输出???,其中:
    所述数据预处理???,用于将所需处理的地震数据处理成叠后地震数据记录;
    所述目标函数建立???,用于构建包括反射系数序列在内的目标函数;
    所述目标函数计算???,用于在所述叠后地震数据记录中选取一道所需处理的地震资料,从所述构建的目标函数中计算出所述地震资料的反射系数序列的振幅和位置;
    所述资料读取???,用于读取下一道所需处理的地震资料,计算出所述下一道地震资料反射系数序列的振幅和位置;
    所述反射系数序列输出???,用于计算出所述叠后地震数据记录中所有地震道反射系数序列的振幅和位置,输出反射系数序列,完成所述叠后地震数据记录的反褶积处理。

    11.  如权利要求10所述的一种反褶积装置,其特征在于,所述目标函数建立??楣菇ǖ哪勘旰?BR>Jα=||Aα-S||2
    上式中,Jα为关于反射系数序列α的范数形式的函数,A为一个N×M的矩阵,其表达式为其中,ωt为地震子波,t=1,2,…N,N为所述叠后地震数据记录的时间窗长度,M为预先设置的所述时间窗长度N内反射系数序列的个数,τM第M个反射系数的时延,S为t时刻所述叠后地震数据记录的值。

    12.  如权利要求11所述的一种反褶积装置,其特征在于,所述目标函数建立??楣菇ǖ哪勘旰?BR>J=||Aα-S||2+μ||Cα-ε||2=Jα+μJc
    上式中,μ为阻抗权系数,C为左下三角全为1的矩阵,ε为波阻抗。

    13.  如权利要求12所述的一种反褶积装置,其特征在于,所述目标函数建立??楣菇ǖ哪勘旰?BR>J=Jα+μJc+βJα
    上式中,β为正则化的阻尼系数,其中,β=β0f0,β0为预设值,f0=max{(ATA+μCTCii)},i为矩阵(ATA+μCTC)的下标。

    14.  如权利要求11所述的一种反褶积装置,其特征在于,所述目标函数计算??榘ㄕ穹扑隳??、位置计算???,其中:
    所述振幅计算???,用于从所述目标函数建立??橹屑扑愠鏊龅卣鹱柿现蟹瓷湎凳蛄械恼穹?;
    所述位置计算???,用于从所述目标函数建立??橹屑扑愠鏊龅卣鹱柿现蟹瓷湎凳蛄械氖奔湮恢?。

    15.  如权利要求14所述的一种反褶积装置,其特征在于,位置计算??榘ǔ跏寄??、扰动???、迭代???、收敛???,其中:
    所述初始???,用于在所述时间窗长度N内随机给定反射系数序列a(j)的初始时间位置,由所述初始时间位置计算所述目标函数Ja当前解,j=1,2…M;
    在所述扰动???,用于在所述时间窗长度N内对所述反射系数序列a(j)进行扰动,将所述扰动后的反射系数序列a(j)代入所述目标函数Ja中,计算出所述目标函数Ja的新解;
    所述迭代???,用于将所述Jα的新解与当前解进行比较,根据Metropolis准则判断是否接受所述目标函数Jα的新解;
    如果判断结果为接受所述新解,则接受反射系数α(j)当前的时间位置;此时,所述新解作为目标函数Ja当前解;
    所述收敛???,用于继续计算所述新解,并与所述当前解进行比较直到达到收敛条件后,停止迭代,得到所述反射系数α(j)的时间位置。

    关 键 词:
    一种 反褶积 方法 装置
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    0条评论

    还可以输入200字符

    暂无评论,赶快抢占沙发吧。

    关于本文
    本文标题:一种反褶积方法及装置.pdf
    链接地址://www.4mum.com.cn/p-6140788.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