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

    重庆时时彩开奖结果视频: 一种基于任意广角波动方程的成像方法及装置.pdf

    关 键 词:
    一种 基于 任意 广角 波动 方程 成像 方法 装置
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201310226448.0

    申请日:

    2013.06.07

    公开号:

    CN103308941A

    公开日:

    2013.09.18

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01V 1/28申请日:20130607|||公开
    IPC分类号: G01V1/28 主分类号: G01V1/28
    申请人: 中国石油天然气集团公司; 中国石油大学(北京)
    发明人: 周辉; 林鹤; 陈汉明; 王颖; 李彦奇
    地址: 100007 北京市东城区东直门北大街9号
    优先权:
    专利代理机构: 北京三友知识产权代理有限公司 11127 代理人: 郭智
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201310226448.0

    授权公告号:

    ||||||

    法律状态公告日:

    2015.10.28|||2013.10.23|||2013.09.18

    法律状态类型:

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

    摘要

    本发明实施例提供一种基于任意广角波动方程的成像方法及装置,所述方法包括:获取建立的速度模型和给定的震源波??;利用正演模拟得到速度模型的合成地震记录;对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理;从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度;采用有限差分方法求解最优化参考速度的任意广角波动方程;采用互相关成像条件进行偏移成像,以提高任意广角波动方程的偏移成像精度??梢杂行У匮怪迫我夤憬遣ǘ匠痰吧疃绕频馁渴挪ǔ上裨胍?,提高了成像质量;所提出的参考速度优化方法,提高了任意广角波动方程叠前深度偏移成像精度,因此可以对更陡倾角的地层成像。

    权利要求书

    权利要求书
    1.   一种基于任意广角波动方程的成像方法,其特征在于,所述基于任意广角波动方程的成像的方法,包括:
    获取建立的速度模型和给定的震源波??;
    利用正演模拟得到速度模型的合成地震记录;
    对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理;
    从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度;
    采用有限差分方法求解最优化参考速度的任意广角波动方程;
    采用互相关成像条件进行偏移成像。

    2.   如权利要求1所述基于任意广角波动方程的成像方法,其特征在于,所述利用正演模拟得到速度模型的合成地震记录,包括:
    采用交错网格有限差分方法求解时空域各向同性均匀介质的一阶速度—应力声波波动方程,从而得到所述速度模型的合成地震记录。

    3.   如权利要求1所述基于任意广角波动方程的成像方法,其特征在于,所述对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理,包括:
    利用频率波数域的扇形滤波器对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理。

    4.   如权利要求1所述基于任意广角波动方程的成像方法,其特征在于,所述从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度,包括:
    从频散关系出发,分析任意广角波动方程的成像精度,利用非线性牛顿法得到任意广角波动方程的最优化参考速度。

    5.   如权利要求1所述基于任意广角波动方程的成像方法,其特征在于,所述采用有限差分方法求解最优化参考速度的任意广角波动方程,包括:
    采用有限差分方法求解最优化参考速度的任意广角波动方程,并在边界处采用完美匹配层PML吸收边界条件。

    6.   一种基于任意广角波动方程的成像装置,其特征在于,所述基于任意广角波动方程的成像的装置,包括:
    获取单元,用于获取建立的速度模型和给定的震源波??;
    正演模拟单元,用于利用正演模拟得到速度模型的合成地震记录;
    倾角滤波单元,用于对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理;
    分析单元,用于从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度;
    方程求解单元,用于采用有限差分方法求解最优化参考速度的任意广角波动方程;
    成像单元,用于采用互相关成像条件进行偏移成像。

    7.   如权利要求6所述基于任意广角波动方程的成像装置,其特征在于,
    所述正演模拟单元,进一步用于采用交错网格有限差分方法求解时空域各向同性均匀介质的一阶速度—应力声波波动方程,从而得到所述速度模型的合成地震记录。

    8.   如权利要求6所述基于任意广角波动方程的成像装置,其特征在于,
    所述倾角滤波单元,进一步用于利用频率波数域的扇形滤波器对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理。

    9.   如权利要求6所述基于任意广角波动方程的成像装置,其特征在于,
    所述分析单元,进一步用于从频散关系出发,分析任意广角波动方程的成像精度,利用非线性牛顿法得到任意广角波动方程的最优化参考速度。

    10.   如权利要求6所述基于任意广角波动方程的成像装置,其特征在于,
    所述方程求解单元,进一步用于采用有限差分方法求解最优化参考速度的任意广角波动方程,并在边界处采用完美匹配层PML吸收边界条件。

    说明书

    说明书一种基于任意广角波动方程的成像方法及装置
    技术领域
    本发明涉及勘探地球物理领域,具体地说涉及一种基于任意广角波动方程的成像方法及装置。
    背景技术
    勘探目标的复杂化和勘探精度的提高,要求地震偏移成像方法可以对高陡倾角构造和横向速度变化较剧烈的地质体准确成像,这在很大程度上促进了叠前深度偏移方法的发展。
    目前,叠前深度偏移方法主要有基于射线理论的Kirchhoff积分类偏移成像方法和基于波场延拓的波动方程偏移成像方法?;掷喾椒ǖ募负我庖灞冉现惫?,具有灵活的目标处理能力和较高的计算效率,目前已成为工业界普遍采用的叠前深度偏移成像方法。但是高频近似假设使其在地下介质比较复杂的情况下得不到满意的成像结果?;诓ǔ⊙油氐牟ǘ匠唐品椒梢院芎玫乇A舻卣鸩ǖ脑硕Ш投ρ卣?,在对复杂构造区域成像时,具有一定优势。
    基于波场延拓的波动方程偏移方法主要包含单程波方程和双程波方程两类:基于双程波方程的逆时偏移方法(RTM)完全遵循全波方程,不存在偏移倾角限制,可以处理速度场的剧烈变化,在保幅处理方面也存在一定优势。但是其在计算量和数据读写方面要求很高,因此至今仍未得到工业界大规模广泛应用。与逆时偏移相比,单程波方法的计算效率较高,但是因为存在对平方根算子的近似,所以对成像倾角和速度的横向变化剧烈程度有一定限制。
    针对传统偏移方法中存在的一些问题,Guddati(2005)提出基于任意广角波动方程(AWWE)的偏移方法。与逆时偏移方法相比,该偏移方法可以避免低频成像噪音和大存储量问题;与传统的单程波偏移方法相比,该偏移方法更适用于速度横向变化较大的复杂构造区以及陡倾角界面的成像。但是因为AWWE存在对平方根算子的近似,所以它的成像精度有待提高。
    发明内容
    本发明实施例提供一种基于任意广角波动方程的成像方法及装置,所提出的倾角滤波方法和参考速度优化方法可以提高任意广角波动方程的偏移成像精度。
    一方面,本发明实施例提供了一种基于任意广角波动方程的成像方法,所述基于任意广角波动方程的成像的方法,包括:
    获取建立的速度模型和给定的震源波??;
    利用正演模拟得到速度模型的合成地震记录;
    对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理;
    从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度;
    采用有限差分方法求解最优化参考速度的任意广角波动方程;
    采用互相关成像条件进行偏移成像。
    可选的,在本发明一实施例中,所述利用正演模拟得到速度模型的合成地震记录,包括:采用交错网格有限差分方法求解时空域各向同性均匀介质的一阶速度—应力声波波动方程,从而得到所述速度模型的合成地震记录。
    可选的,在本发明一实施例中,所述对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理,包括:利用频率波数域的扇形滤波器对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理。
    可选的,在本发明一实施例中,所述从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度,包括:从频散关系出发,分析任意广角波动方程的成像精度,利用非线性牛顿法得到任意广角波动方程的最优化参考速度。
    可选的,在本发明一实施例中,所述采用有限差分方法求解最优化参考速度的任意广角波动方程,包括:采用有限差分方法求解最优化参考速度的任意广角波动方程,并在边界处采用完美匹配层PML吸收边界条件。
    另一方面,本发明实施例提供了一种基于任意广角波动方程的成像装置,所述基于任意广角波动方程的成像的装置,包括:
    获取单元,用于获取建立的速度模型和给定的震源波??;
    正演模拟单元,用于利用正演模拟得到速度模型的合成地震记录;
    倾角滤波单元,用于对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理;
    分析单元,用于从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度;
    方程求解单元,用于采用有限差分方法求解最优化参考速度的任意广角波动方程;
    成像单元,用于采用互相关成像条件进行偏移成像。
    可选的,在本发明一实施例中,所述正演模拟单元,进一步用于采用交错网格有限差分方法求解时空域各向同性均匀介质的一阶速度—应力声波波动方程,从而得到所述速度模型的合成地震记录。
    可选的,在本发明一实施例中,所述倾角滤波单元,进一步用于利用频率波数域扇形滤波器对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理。
    可选的,在本发明一实施例中,所述分析单元,进一步用于从频散关系出发,分析任意广角波动方程的成像精度,利用非线性牛顿法得到任意广角波动方程的最优化参考速度。
    可选的,在本发明一实施例中,所述方程求解单元,进一步用于采用有限差分方法求解最优化参考速度的任意广角波动方程,并在边界处采用完美匹配层PML吸收边界条件。
    上述技术方案具有如下有益效果:因为采用所述基于任意广角波动方程的成像的方法,包括:获取建立的速度模型和给定的震源波??;利用正演模拟得到速度模型的合成地震记录;对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理;从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度;采用有限差分方法求解最优化参考速度的任意广角波动方程;采用互相关成像条件进行偏移成像的技术手段,所以达到了如下的技术效果:(1)提出的倾角滤波方法可以有效地压制任意广角波动方程叠前深度偏移的倏逝波成像噪音,提高成像质量;(2)优化了任意广角波动方程的参考速度参数,提高了任意广角波动方程叠前深度偏移成像精度,因此使基于任意广角波动方程的叠前深度偏移方法可以对更陡倾角的地层成像。
    附图说明
    为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
    图1为本发明实施例提供了一种基于任意广角波动方程的成像方法流程图;
    图2为本发明实施例一种基于任意广角波动方程的成像装置结构示意图;
    图3(a)为本发明应用实例AWWE的频散关系曲线实部之间的关系示意图;
    图3(b)为本发明应用实例AWWE的频散关系曲线虚部之间的关系示意图;
    图4为本发明应用实例频率波数域扇形滤波器的频率波数响应示意图;
    图5(a)为本发明应用实例优化二参数的任意广角波动方程的频散关系曲线示意图;
    图5(b)为本发明应用实例优化三参数的任意广角波动方程的频散关系曲线示意图;
    图6为本发明应用实例复杂介质速度模型示意图;
    图7为本发明应用实例复杂介质速度模型的多炮偏移成像结果示意图。
    具体实施方式
    下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例?;诒痉⒚髦械氖凳├?,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明?;さ姆段?。
    如图1所示,为本发明实施例提供了一种基于任意广角波动方程的成像方法流程图,所述基于任意广角波动方程的成像的方法,包括:
    101、获取建立的速度模型和给定的震源波??;
    102、利用正演模拟得到速度模型的合成地震记录;
    103、对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理;
    104、从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度;
    105、采用有限差分方法求解最优化参考速度的任意广角波动方程;
    106、采用互相关成像条件进行偏移成像。
    可选的,所述利用正演模拟得到速度模型的合成地震记录,包括:采用交错网格有限差分方法求解时空域各向同性均匀介质的一阶速度—应力声波波动方程,从而得到所述速度模型的合成地震记录。
    可选的,所述对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理,包括:利用频率波数域的扇形滤波器对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理。
    可选的,所述从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度,包括:从频散关系出发,分析任意广角波动方程的成像精度,利用非线性牛顿法得到任意广角波动方程的最优化参考速度。
    可选的,所述采用有限差分方法求解最优化参考速度的任意广角波动方程,包括:采用有限差分方法求解最优化参考速度的任意广角波动方程,并在边界处采用完美匹配层PML(Perfectly Matched Layer)吸收边界条件。
    如图2所示,为本发明实施例一种基于任意广角波动方程的成像装置结构示意图,所述基于任意广角波动方程的成像的装置,包括:
    获取单元21,用于获取建立的速度模型和给定的震源波??;
    正演模拟单元22,用于利用正演模拟得到速度模型的合成地震记录;
    倾角滤波单元23,用于对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理;
    分析单元24,用于从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度;
    方程求解单元25,用于采用有限差分方法求解最优化参考速度的任意广角波动方程;
    成像单元26,用于采用互相关成像条件进行偏移成像。
    可选的,所述正演模拟单元22,进一步用于采用交错网格有限差分方法求解时空域各向同性均匀介质的一阶速度—应力声波波动方程,从而得到所述速度模型的合成地震记录。
    可选的,所述倾角滤波单元23,进一步用于利用频率波数域的扇形滤波器对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理。
    可选的,所述分析单元24,进一步用于从频散关系出发,分析任意广角波动方程的成像精度,利用非线性牛顿法得到任意广角波动方程的最优化参考速度。
    可选的,所述方程求解单元25,进一步用于采用有限差分方法求解最优化参考速度的任意广角波动方程,并在边界处采用完美匹配层PML吸收边界条件。
    本发明实施例上述技术方案具有如下有益效果:因为采用所述基于任意广角波动方程的成像的方法,包括:获取建立的速度模型和给定的震源波??;利用正演模拟得到速度模型的合成地震记录;对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理;从频散关系出发,分析任意广角波动方程的成像精度,得到任意广角波动方程的最优化参考速度;采用有限差分方法求解最优化参考速度的任意广角波动方程;采用互相关成像条件进行偏移成像的技术手段,所以达到了如下的技术效果:(1)提出的倾角滤波方法可以有效地压制任意广角波动方程叠前深度偏移的倏逝波成像噪音,提高成像质量;(2)优化了任意广角波动方程的参考速度参数,提高了任意广角波动方程叠前深度偏移成像精度,因此使基于任意广角波动方程的叠前深度偏移方法可以对更陡倾角的地层成像。
    以下举应用实例进行详细说明:
    本发明应用实例的目的在于提供一种提高任意广角波动方程成像精度的方法。提出了一种基于倾角滤波的成像噪音压制方法,提高了偏移成像质量。从频散关系出发,用非线性牛顿法求解单程波算子的最优化问题,使任意广角波动方程的频散关系与准确的平方根算子达到最佳匹配,提高了方程的成像倾角,使之适合更陡倾角界面的成像。
    第一步,建立速度模型,给定震源波场。
    根据设计要求建立一个理论速度模型,设定好地下介质速度。根据震源子波给定震源波场。
    第二步,正演模拟得到速度模型的合成地震记录。
    采用交错网格有限差分方法求解时空域各向同性均匀介质的一阶速度—应力声波波动方程,从而得到所设计速度模型的合成地震记录。均匀各项同性介质一阶速度—应力声波方程的形式为
    <mrow><MFRAC><MROW><MO>&amp;PartialD;</MO> <MI>p</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MI>t</MI> </MROW></MFRAC><MO>=</MO> <MO>-</MO> <MI>K</MI> <MROW><MO>(</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>v</MI> <MI>x</MI> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MI>x</MI> </MROW></MFRAC><MO>+</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>v</MI> <MI>z</MI> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MI>z</MI> </MROW></MFRAC><MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR><MATHS num="0002"><MATH><![CDATA[ <mrow><MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>v</MI> <MI>x</MI> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MI>t</MI> </MROW></MFRAC><MO>=</MO> <MO>-</MO> <MFRAC><MN>1</MN> <MI>ρ</MI> </MFRAC><MFRAC><MROW><MO>&amp;PartialD;</MO> <MI>p</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MI>x</MI> </MROW></MFRAC><MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>1</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR><MATHS num="0003"><MATH><![CDATA[ <mrow><MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>v</MI> <MI>z</MI> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MI>t</MI> </MROW></MFRAC><MO>=</MO> <MO>-</MO> <MFRAC><MN>1</MN> <MI>ρ</MI> </MFRAC><MFRAC><MROW><MO>&amp;PartialD;</MO> <MI>p</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MI>z</MI> </MROW></MFRAC></MROW>]]&gt;</MATH></MATHS> <BR>其中K,ρ为介质的体积模量和密度,它们与声波传播速度V之间的关系为K=ρV2,P,vx,vz为介质中压力和质点振动速度在两个坐标方向的分量,x,z分别为水平和垂直坐标,t为时间变量。 <BR>第三步,对正演模拟的合成地震记录和给定的震源波场做倾角滤波处理。 <BR>如图3(a)所示,为本发明应用实例AWWE的频散关系曲线实部之间的关系示意图,如图3(b)所示,为本发明应用实例AWWE的频散关系曲线虚部之间的关系示意图,其中exact表示准确的频散关系,AWWE表示任意广角波动方程的频散关系。通过分析图3(a)、图3(b)的频散关系曲线和一系列的推导论证,得出偏移成像噪音是由于不准确的倏逝波近似产生的。我们根据倏逝波的慢度特征,在频率—波数域设计滤波器,将频率波数域的扇形滤波器作用于接收到的记录波场和给定的震源波场,这样就可以滤除不准确倏逝波。频率波数域的扇形滤波器的表达式为 <BR><MATHS num="0004"><MATH><![CDATA[ <mrow><MI>H</MI> <MROW><MO>(</MO> <MI>ω</MI> <MO>,</MO> <MI>k</MI> <MO>)</MO> </MROW><MO>=</MO> <MFENCED close="" open="{"><MTABLE><MTR><MTD><MN>1</MN> <MO>,</MO> </MTD><MTD><MO>|</MO> <MI>k</MI> <MO>|</MO> <MO>≤</MO> <MSUB><MI>k</MI> <MN>1</MN> </MSUB></MTD></MTR><MTR><MTD><MI>g</MI> <MROW><MO>(</MO> <MI>k</MI> <MO>)</MO> </MROW><MO>,</MO> </MTD><MTD><MSUB><MI>k</MI> <MN>1</MN> </MSUB><MO>&lt;</MO> <MO>|</MO> <MI>k</MI> <MO>|</MO> <MO>&lt;</MO> <MSUB><MI>k</MI> <MN>2</MN> </MSUB></MTD></MTR><MTR><MTD><MN>0</MN> <MO>,</MO> </MTD><MTD><MO>|</MO> <MI>k</MI> <MO>|</MO> <MO>&amp;GreaterEqual;</MO> <MSUB><MI>k</MI> <MN>2</MN> </MSUB></MTD></MTR></MTABLE></MFENCED><MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>2</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>其中,k1=ω/c?10dk,k2=ω/c+10dk,H(ω,k)为频率波数域的扇形滤波器的频率—波数响应函数,k为圆波数,c为检波器所在位置处的地层速度,g(k)为余弦镶边函数,k1,k2为镶边函数的左右边界,镶边函数的长度取为二十倍采样波数间隔dk。如图4所示,为本发明应用实例频率波数域的扇形滤波器的频率波数响应示意图。为了避免吉普斯效应,我们对频率波数域的扇形滤波器进行镶边,采用的余弦镶边函数为 <BR><MATHS num="0005"><MATH><![CDATA[ <mrow><MI>g</MI> <MROW><MO>(</MO> <MI>ω</MI> <MO>,</MO> <MI>k</MI> <MO>)</MO> </MROW><MO>=</MO> <MFRAC><MN>1</MN> <MN>2</MN> </MFRAC><MROW><MO>(</MO> <MN>1</MN> <MO>+</MO> <MI>cos</MI> <MFRAC><MROW><MO>|</MO> <MI>k</MI> <MO>|</MO> <MO>-</MO> <MSUB><MI>k</MI> <MN>1</MN> </MSUB></MROW><MROW><MSUB><MI>k</MI> <MN>2</MN> </MSUB><MO>-</MO> <MSUB><MI>k</MI> <MN>1</MN> </MSUB></MROW></MFRAC><MI>π</MI> <MO>)</MO> </MROW><MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>3</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>第四步,从频散关系出发,分析任意广角波动方程的成像精度,利用非线性牛顿法求解任意广角波动方程参考速度的最优化问题,得到任意广角波动方程的最优化参考速度。 <BR>时空域上行任意广角波动方程为 <BR><MATHS num="0006"><MATH><![CDATA[ <mrow><MFRAC><MROW><MSUP><MO>&amp;PartialD;</MO> <MN>2</MN> </MSUP><MI>u</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MI>z</MI> <MO>&amp;PartialD;</MO> <MI>t</MI> </MROW></MFRAC><MO>-</MO> <MFRAC><MN>1</MN> <MI>c</MI> </MFRAC><MROW><MO>(</MO> <MSUB><MI>Λ</MI> <MN>1</MN> </MSUB><MO>+</MO> <MSUB><MI>Λ</MI> <MN>2</MN> </MSUB><MO>)</MO> </MROW><MFRAC><MROW><MSUP><MO>&amp;PartialD;</MO> <MN>2</MN> </MSUP><MI>u</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MSUP><MI>t</MI> <MN>2</MN> </MSUP></MROW></MFRAC><MO>+</MO> <MI>c</MI> <MSUB><MI>Λ</MI> <MN>2</MN> </MSUB><MFRAC><MROW><MSUP><MO>&amp;PartialD;</MO> <MN>2</MN> </MSUP><MI>u</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MSUP><MI>x</MI> <MN>2</MN> </MSUP></MROW></MFRAC><MO>=</MO> <MN>0</MN> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>4</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>其中,d=(1&nbsp;0&nbsp;0&nbsp;…&nbsp;0)T,u=(u&nbsp;u1&nbsp;u2&nbsp;…&nbsp;un?1)T,T表示转置,u表示波场值,u1,u2,…,un?1是推导(4)式时引入的辅助变量;n?1表示辅助变量的个数,x,z分别表示空间水平和深度坐标,t表示时间变量,c是地震波在地下介质中的传播速度, <BR> <BR> <BR>c1,c2,…,cn?1,cn是推导(4)式引入的参考速度。 <BR>波动方程的理论频散关系为 <BR><MATHS num="0007"><MATH><![CDATA[ <mrow><MSUBSUP><MI>k</MI> <MI>x</MI> <MN>2</MN> </MSUBSUP><MO>+</MO> <MSUBSUP><MI>k</MI> <MI>z</MI> <MN>2</MN> </MSUBSUP><MO>=</MO> <MFRAC><MSUP><MI>ω</MI> <MN>2</MN> </MSUP><MSUP><MI>c</MI> <MN>2</MN> </MSUP></MFRAC><MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>7</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>其中kx是水平波数,kz是垂直波数,ω是角频率。令X=ckx/ω=sinθ,R=ckz/ω,则方程(7)可以简化为 <BR><MATHS num="0008"><MATH><![CDATA[ <mrow><MI>R</MI> <MO>=</MO> <MSQRT><MN>1</MN> <MO>-</MO> <MSUP><MI>X</MI> <MN>2</MN> </MSUP></MSQRT><MO>=</MO> <MI>cos</MI> <MI>θ</MI> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>8</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>其中θ表示波的传播方向与界面法线之间的夹角。对方程(4)中的x、z、t做傅里叶变换,得到频率—波数域中的任意广角波动方程 <BR><MATHS num="0009"><MATH><![CDATA[ <mrow><MSUB><MI>k</MI> <MI>z</MI> </MSUB><MI>du</MI> <MO>-</MO> <MFRAC><MI>ω</MI> <MI>c</MI> </MFRAC><MO>{</MO> <MSUB><MI>Λ</MI> <MN>1</MN> </MSUB><MO>+</MO> <MO>[</MO> <MN>1</MN> <MO>-</MO> <MSUP><MROW><MO>(</MO> <MFRAC><MROW><MI>c</MI> <MSUB><MI>k</MI> <MI>x</MI> </MSUB></MROW><MI>ω</MI> </MFRAC><MO>)</MO> </MROW><MN>2</MN> </MSUP><MO>]</MO> <MSUB><MI>Λ</MI> <MN>2</MN> </MSUB><MO>}</MO> <MI>u</MI> <MO>=</MO> <MN>0</MN> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>9</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>整理方程(9)可以得到任意广角波动方程对平方根算子的近似函数Fn(θ,c1,…,cn)。当选取两个参考速度,并归一化速度后得到 <BR><MATHS num="0010"><MATH><![CDATA[ <mrow><MSUB><MI>F</MI> <MN>2</MN> </MSUB><MROW><MO>(</MO> <MI>θ</MI> <MO>,</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>,</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB><MO>)</MO> </MROW><MO>=</MO> <MROW><MO>(</MO> <MFRAC><MN>1</MN> <MROW><MN>2</MN> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MROW></MFRAC><MO>+</MO> <MFRAC><MSUB><MI>c</MI> <MN>1</MN> </MSUB><MN>2</MN> </MFRAC><MO>)</MO> </MROW><MO>-</MO> <MFRAC><MSUB><MI>c</MI> <MN>1</MN> </MSUB><MN>2</MN> </MFRAC><MSUP><MI>sin</MI> <MN>2</MN> </MSUP><MI>θ</MI> <MO>-</MO> <MFRAC><MSUP><MROW><MO>[</MO> <MROW><MO>(</MO> <MFRAC><MN>1</MN> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MFRAC><MO>-</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>)</MO> </MROW><MO>+</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MSUP><MI>sin</MI> <MN>2</MN> </MSUP><MI>θ</MI> <MO>]</MO> </MROW><MN>2</MN> </MSUP><MROW><MN>2</MN> <MO>[</MO> <MROW><MO>(</MO> <MFRAC><MN>1</MN> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MFRAC><MO>+</MO> <MFRAC><MN>1</MN> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MFRAC><MO>)</MO> </MROW><MO>+</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>+</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB><MO>]</MO> <MO>-</MO> <MN>2</MN> <MROW><MO>(</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>+</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB><MO>)</MO> </MROW><MSUP><MI>sin</MI> <MN>2</MN> </MSUP><MI>θ</MI> </MROW></MFRAC><MO>.</MO> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>10</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>从平方根算子表达的频散关系出发,对参考速度进行优化。优化的准则是在最小平方意义下,使任意广角波动方程对平方根算子的近似与准确的平方根算子达到最佳的匹配。下面在两个参考速度的情况下,详细地说明参数的优化过程。(8)式与(10)式之间的相对近似误差为 <BR><MATHS num="0011"><MATH><![CDATA[ <mrow><MSUB><MI>E</MI> <MI>r</MI> </MSUB><MO>=</MO> <MSUB><MI>F</MI> <MN>2</MN> </MSUB><MO>-</MO> <MI>R</MI> <MO>=</MO> <MROW><MO>(</MO> <MFRAC><MN>1</MN> <MROW><MN>2</MN> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MROW></MFRAC><MO>+</MO> <MFRAC><MSUB><MI>c</MI> <MN>1</MN> </MSUB><MN>2</MN> </MFRAC><MO>)</MO> </MROW><MO>-</MO> <MFRAC><MSUB><MI>c</MI> <MN>1</MN> </MSUB><MN>2</MN> </MFRAC><MSUP><MI>sin</MI> <MN>2</MN> </MSUP><MI>θ</MI> <MO>-</MO> <MFRAC><MSUP><MROW><MO>[</MO> <MROW><MO>(</MO> <MFRAC><MN>1</MN> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MFRAC><MO>-</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>)</MO> </MROW><MO>+</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MSUP><MI>sin</MI> <MN>2</MN> </MSUP><MI>θ</MI> <MO>]</MO> </MROW><MN>2</MN> </MSUP><MROW><MN>2</MN> <MO>[</MO> <MROW><MO>(</MO> <MFRAC><MN>1</MN> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MFRAC><MO>+</MO> <MFRAC><MN>1</MN> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MFRAC><MO>)</MO> </MROW><MO>+</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>+</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB><MO>]</MO> <MO>-</MO> <MN>2</MN> <MROW><MO>(</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>+</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB><MO>)</MO> </MROW><MSUP><MI>sin</MI> <MN>2</MN> </MSUP><MI>θ</MI> </MROW></MFRAC><MO>-</MO> <MI>cos</MI> <MI>θ</MI> <MO>.</MO> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>11</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>定义目标函数 <BR><MATHS num="0012"><MATH><![CDATA[ <mrow><MI>Ir</MI> <MROW><MO>(</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>,</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB><MO>)</MO> </MROW><MO>=</MO> <MSUBSUP><MO>&amp;Integral;</MO> <MSUB><MI>θ</MI> <MN>1</MN> </MSUB><MSUB><MI>θ</MI> <MN>2</MN> </MSUB></MSUBSUP><MSUBSUP><MI>E</MI> <MI>r</MI> <MN>2</MN> </MSUBSUP><MROW><MO>(</MO> <MI>θ</MI> <MO>)</MO> </MROW><MI>dθ</MI> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>12</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>其中θ1和θ2是所选择的角度范围,通常选取θ1=0°,θ2=90°。通过求取目标函数Ir的最小值,得到参考速度c1和c2的值。在最小平方意义约束下,得到 <BR><MATHS num="0013"><MATH><![CDATA[ <mrow><MFENCED close="" open="{"><MTABLE><MTR><MTD><MFRAC><MROW><MO>&amp;PartialD;</MO> <MI>Ir</MI> <MROW><MO>(</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>,</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB><MO>)</MO> </MROW></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MROW></MFRAC><MO>=</MO> <MN>0</MN> </MTD></MTR><MTR><MTD><MFRAC><MROW><MO>&amp;PartialD;</MO> <MI>Ir</MI> <MROW><MO>(</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>,</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB><MO>)</MO> </MROW></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MROW></MFRAC><MO>=</MO> <MN>0</MN> </MTD></MTR></MTABLE></MFENCED><MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>13</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>采用牛顿法求解非线性方程组(13)。首先给定c1和c2的初值和然后采用(14)式更新参考速度。 <BR><MATHS num="0014"><MATH><![CDATA[ <mrow><MFENCED close="" open="{"><MTABLE><MTR><MTD><MSUBSUP><MI>c</MI> <MN>1</MN> <MROW><MO>(</MO> <MI>k</MI> <MO>+</MO> <MN>1</MN> <MO>)</MO> </MROW></MSUBSUP><MO>=</MO> <MSUBSUP><MI>c</MI> <MN>1</MN> <MROW><MO>(</MO> <MI>k</MI> <MO>)</MO> </MROW></MSUBSUP><MO>+</MO> <MI>Δ</MI> <MSUP><MSUB><MI>c</MI> <MN>1</MN> </MSUB><MROW><MO>(</MO> <MI>k</MI> <MO>)</MO> </MROW></MSUP></MTD></MTR><MTR><MTD><MSUBSUP><MI>c</MI> <MN>2</MN> <MROW><MO>(</MO> <MI>k</MI> <MO>+</MO> <MN>1</MN> <MO>)</MO> </MROW></MSUBSUP><MO>=</MO> <MSUBSUP><MI>c</MI> <MN>2</MN> <MROW><MO>(</MO> <MI>k</MI> <MO>)</MO> </MROW></MSUBSUP><MO>+</MO> <MI>Δ</MI> <MSUP><MSUB><MI>c</MI> <MN>2</MN> </MSUB><MROW><MO>(</MO> <MI>k</MI> <MO>)</MO> </MROW></MSUP></MTD></MTR></MTABLE></MFENCED><MO>,</MO> <MI>k</MI> <MO>=</MO> <MN>0,1,2</MN> <MO>.</MO> <MO>.</MO> <MO>.</MO> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>14</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>其中表示第k次的参考速度迭代值,Δc1(k),Δc2(k)表示第k次的参考速度更正值,是第k+1次的参考速度迭代值,参考速度的更正值Δc1(k),Δc2(k)由(15)式求出。 <BR> <BR>其中Δc(k)=[Δc1(k),Δc2(k)]是参考速度更正向量,g(c(k))=[g(c1),g(c2)](k)是目标函数的梯度向量,<MATHS num="0015"><MATH><![CDATA[ <mrow> <MSUP><MI>H</MI> <MROW><MO>(</MO> <MI>k</MI> <MO>)</MO> </MROW></MSUP><MO>=</MO> <MSUP><MFENCED close="]" open="["><MTABLE><MTR><MTD><MSUB><MI>h</MI> <MN>11</MN> </MSUB></MTD><MTD><MSUB><MI>h</MI> <MN>12</MN> </MSUB></MTD></MTR><MTR><MTD><MSUB><MI>h</MI> <MN>21</MN> </MSUB></MTD><MTD><MSUB><MI>h</MI> <MN>22</MN> </MSUB></MTD></MTR></MTABLE></MFENCED><MROW><MO>(</MO> <MI>k</MI> <MO>)</MO> </MROW></MSUP></MROW>]]&gt;</MATH></MATHS>是目标函数的Hessian矩阵,梯度向量和Hessian矩阵中各元素由(16)—(20)式定义,即 <BR><MATHS num="0016"><MATH><![CDATA[ <mrow><MI>g</MI> <MROW><MO>(</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>)</MO> </MROW><MO>=</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MI>Ir</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MROW></MFRAC><MO>=</MO> <MSUBSUP><MO>&amp;Integral;</MO> <MSUB><MI>θ</MI> <MN>1</MN> </MSUB><MSUB><MI>θ</MI> <MN>2</MN> </MSUB></MSUBSUP><MN>2</MN> <MSUB><MI>E</MI> <MI>r</MI> </MSUB><MROW><MO>(</MO> <MI>θ</MI> <MO>)</MO> </MROW><MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>E</MI> <MI>r</MI> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MROW></MFRAC><MI>dθ</MI> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>16</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR><MATHS num="0017"><MATH><![CDATA[ <mrow><MI>g</MI> <MROW><MO>(</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB><MO>)</MO> </MROW><MO>=</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MI>Ir</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MROW></MFRAC><MO>=</MO> <MSUBSUP><MO>&amp;Integral;</MO> <MSUB><MI>θ</MI> <MN>1</MN> </MSUB><MSUB><MI>θ</MI> <MN>2</MN> </MSUB></MSUBSUP><MN>2</MN> <MSUB><MI>E</MI> <MI>r</MI> </MSUB><MROW><MO>(</MO> <MI>θ</MI> <MO>)</MO> </MROW><MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>E</MI> <MI>r</MI> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MROW></MFRAC><MI>dθ</MI> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>17</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR><MATHS num="0018"><MATH><![CDATA[ <mrow><MSUB><MI>h</MI> <MN>11</MN> </MSUB><MO>=</MO> <MFRAC><MROW><MSUP><MO>&amp;PartialD;</MO> <MN>2</MN> </MSUP><MI>Ir</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MSUP><MSUB><MI>c</MI> <MN>1</MN> </MSUB><MN>2</MN> </MSUP></MROW></MFRAC><MO>=</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>g</MI> <MN>1</MN> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MROW></MFRAC><MO>=</MO> <MN>2</MN> <MSUBSUP><MO>&amp;Integral;</MO> <MSUB><MI>θ</MI> <MN>1</MN> </MSUB><MSUB><MI>θ</MI> <MN>2</MN> </MSUB></MSUBSUP><MO>[</MO> <MSUP><MROW><MO>(</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>E</MI> <MI>r</MI> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MROW></MFRAC><MO>)</MO> </MROW><MN>2</MN> </MSUP><MO>+</MO> <MSUB><MI>E</MI> <MI>r</MI> </MSUB><MROW><MO>(</MO> <MI>θ</MI> <MO>)</MO> </MROW><MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUBSUP><MI>E</MI> <MI>r</MI> <MN>2</MN> </MSUBSUP></MROW><MROW><MO>&amp;PartialD;</MO> <MSUP><MSUB><MI>c</MI> <MN>1</MN> </MSUB><MN>2</MN> </MSUP></MROW></MFRAC><MO>]</MO> <MI>dθ</MI> <MO>,</MO> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>18</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR><MATHS num="0019"><MATH><![CDATA[ <mrow><MSUB><MI>h</MI> <MN>12</MN> </MSUB><MO>=</MO> <MSUB><MI>h</MI> <MN>21</MN> </MSUB><MO>=</MO> <MFRAC><MROW><MSUP><MO>&amp;PartialD;</MO> <MN>2</MN> </MSUP><MI>Ir</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MROW></MFRAC><MO>=</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>g</MI> <MN>1</MN> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MROW></MFRAC><MO>=</MO> <MN>2</MN> <MSUBSUP><MO>&amp;Integral;</MO> <MSUB><MI>θ</MI> <MN>1</MN> </MSUB><MSUB><MI>θ</MI> <MN>2</MN> </MSUB></MSUBSUP><MO>[</MO> <MROW><MO>(</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>E</MI> <MI>r</MI> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB></MROW></MFRAC><MO>)</MO> </MROW><MROW><MO>(</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>E</MI> <MI>r</MI> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MROW></MFRAC><MO>)</MO> </MROW><MO>+</MO> <MSUB><MI>E</MI> <MI>r</MI> </MSUB><MROW><MO>(</MO> <MI>θ</MI> <MO>)</MO> </MROW><MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUBSUP><MI>E</MI> <MI>r</MI> <MN>2</MN> </MSUBSUP></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>1</MN> </MSUB><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MROW></MFRAC><MO>]</MO> <MI>dθ</MI> <MO>,</MO> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>19</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR><MATHS num="0020"><MATH><![CDATA[ <mrow><MSUB><MI>h</MI> <MN>22</MN> </MSUB><MO>=</MO> <MFRAC><MROW><MSUP><MO>&amp;PartialD;</MO> <MN>2</MN> </MSUP><MI>Ir</MI> </MROW><MROW><MO>&amp;PartialD;</MO> <MSUP><MSUB><MI>c</MI> <MN>2</MN> </MSUB><MN>2</MN> </MSUP></MROW></MFRAC><MO>=</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>g</MI> <MN>2</MN> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MROW></MFRAC><MO>=</MO> <MN>2</MN> <MSUBSUP><MO>&amp;Integral;</MO> <MSUB><MI>θ</MI> <MN>1</MN> </MSUB><MSUB><MI>θ</MI> <MN>2</MN> </MSUB></MSUBSUP><MO>[</MO> <MSUP><MROW><MO>(</MO> <MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>E</MI> <MI>r</MI> </MSUB></MROW><MROW><MO>&amp;PartialD;</MO> <MSUB><MI>c</MI> <MN>2</MN> </MSUB></MROW></MFRAC><MO>)</MO> </MROW><MN>2</MN> </MSUP><MO>+</MO> <MSUB><MI>E</MI> <MI>r</MI> </MSUB><MROW><MO>(</MO> <MI>θ</MI> <MO>)</MO> </MROW><MFRAC><MROW><MO>&amp;PartialD;</MO> <MSUBSUP><MI>E</MI> <MI>r</MI> <MN>2</MN> </MSUBSUP></MROW><MROW><MO>&amp;PartialD;</MO> <MSUP><MSUB><MI>c</MI> <MN>2</MN> </MSUB><MN>2</MN> </MSUP></MROW></MFRAC><MO>]</MO> <MI>dθ</MI> <MO>.</MO> <MO>-</MO> <MO>-</MO> <MO>-</MO> <MROW><MO>(</MO> <MN>20</MN> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>把(11)式代入到(16)—(20)式中,再根据(15)式就可以得到每一次的参考速度更正值。 <BR>当参考速度的更正值小于给定的值(本发明选取10?6)时,终止迭代,并输出此时的参考速度。输出的参考速度可以使任意广角波动方程的精度在所选择的角度范围内达到全局最优。如图5(a)所示,为本发明应用实例优化二参数的任意广角波动方程的频散关系曲线示意图,如图5(b)所示,为本发明应用实例优化三参数的任意广角波动方程的频散关系曲线示意图,可以看出方程的成像倾角得到了大大地提高。 <BR>第五步,采用有限差分方法求解最优化参考速度的任意广角波动方程。 <BR>将第四步得到的参考速度带入到方程(4)中,然后采用有限差分方法求解方程(4),边界处采用PML吸收边界条件,同时将震源波场和记录波场沿着深度方向做波场延拓。 <BR>第六步,采用互相关成像条件进行偏移成像。如图6所示,为本发明应用实例复杂介质速度模型示意图,如图7所示,为本发明应用实例复杂介质速度模型的多炮偏移成像结果示意图,可以看出主要的速度界面都得到了较好的成像。 <BR>本发明应用实例与现有技术相比,具有如下的有益效果:(1)提出的倾角滤波方法可以有效地压制任意广角波动方程叠前深度偏移的倏逝波成像噪音,提高成像质量;(2)优化了任意广角波动方程的参考速度参数,提高了任意广角波动方程叠前深度偏移成像精度,因此使基于任意广角波动方程的叠前深度偏移方法可以对更陡倾角的地层成像。 <BR>本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrative&nbsp;logical&nbsp;block),单元,和步骤可以通过电子硬件、电脑软件,或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability),上述的各种说明性部件(illustrative&nbsp;components),单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例?;さ姆段?。 <BR>本发明实施例中所描述的各种说明性的逻辑块,或单元都可以通过通用处理器,数字信号处理器,专用集成电路(ASIC),现场可编程门阵列或其它可编程逻辑装置,离散门或晶体管逻辑,离散硬件部件,或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器,可选地,该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现,例如数字信号处理器和微处理器,多个微处理器,一个或多个微处理器联合一个数字信号处理器核,或任何其它类似的配置来实现。 <BR>本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件???、或者这两者的结合。软件??榭梢源娲⒂赗AM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD?ROM或本领域中其它任意形式的存储媒介中。示例性地,存储媒介可以与处理器连接,以使得处理器可以从存储媒介中读取信息,并可以向存储媒介存写信息??裳〉?,存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中,ASIC可以设置于用户终端中??裳〉?,处理器和存储媒介也可以设置于用户终端中的不同的部件中。 <BR>在一个或多个示例性的设计中,本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD?ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外,任何连接都可以被适当地定义为电脑可读媒介,例如,如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘,磁盘通常以磁性复制数据,而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。 <BR>以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的?;し段?,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的?;し段е?。</p> </div> </div> </div> </div> <div class="mt10 works-comment"> <div class="works-comment-hd"> <span class="font-tahoma">关于本文</div> <div style="line-height: 25px; padding: 10px 20px;"> 本文标题:一种基于任意广角波动方程的成像方法及装置.pdf<br /> 链接地址:<a href="//www.4mum.com.cn/p-5775527.html">//www.4mum.com.cn/p-5775527.html</a><br /> </div> </div> </div> <div class="boxright" id="boxright" > <div class="fr detail-aside" id="Div11" style="width:270px;"> <div class="box hot-keywords mt10" style="overflow: hidden;width: 268px; border:solid 1px #dedede;" id="relatebox0"> <div class="boxHd" > <div class="fl nt-ico mr5 ml13 ico" style="margin-top: 3px;"> </div> <h2 class="fl font-normal font16 font-yahei" style="font-size: 16px; font-weight: 100; margin-left: 0px; margin-top: 3px; font-family: 微软雅黑"> 当前资源信息</h2> </div> <div id="Div2" class="author-works-list bgF" style="overflow: hidden; padding:10px 10px; "> <table><tr><td> <dt class="author-avatar-box fl"><a class="author-avatar" title="t****" href="u-242.html"> <img src="//www.4mum.com.cn/FlexPaper/ZoomImage/UploadPhoto/tx_20181218222850.png" style="border-radius:50px 50px; " onerror="this.src='images/noavatar_small.gif'" alt="t****"></a> </dt></td><td> <div class="author-name fl w100 ellipsis"> <a href="u-242.html" target="_blank"> t****</a><img style="height:15px; width:20px; overflow:hidden; margin-right:10px;background:url(images/bg_index_ie6_781d95ab.png) no-repeat -185px 4px;display:none" src="images/s.gif" alt="企业认证" title="企业认证"/></div> <div class="author-level-bar"> <span class="mr5 author-grade author-grade5" title="会员等级"></span> </div></td></tr></table> <div class="fl" style="width: 240px;"> <p class="kh_cpZl"> 编号: cj20190318102822324534</p> <p class="kh_cpZl"> 类型: 共享资源</p> <p class="kh_cpZl"> 格式: PDF</p> <p class="kh_cpZl"> 大?。?5.10MB</p> <p class="kh_cpZl"> 上传时间: 2019-03-18</p> </div> </div> </div> <div class="box hot-keywords mt10" style="overflow: hidden;width: 268px; border:solid 1px #dedede;" id="relatebox"> <div class="boxHd" style="padding-bottom: 0px;"> <div class="fl keywords-ico mr5 ml13 ico"> </div> <h2 class="fl font-normal font16 font-yahei" style="font-size: 16px; font-weight: 100; margin-left: 0px; margin-top: 3px; font-family: 微软雅黑"> 相关资源</h2> </div> <div id="author-works-list" class="author-works-list bgF"> <li class="vprdn"> <img alt="枪支在位检测装置.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326895.html" title="枪支在位检测装置.pdf">枪支在位检测装置.pdf</a> </li> <li class="vprdn"> <img alt="便捷装拆式竖向送枪支架.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326824.html" title="便捷装拆式竖向送枪支架.pdf">便捷装拆式竖向送枪支架.pdf</a> </li> <li class="vprdn"> <img alt="一种制导弹药电性能测试仿真试验平台.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326533.html" title="一种制导弹药电性能测试仿真试验平台.pdf">一种制导弹药电性能测试仿真试验平台.pdf</a> </li> <li class="vprdn"> <img alt="一种制导弹药舵机放大器仿真测试装置.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326506.html" title="一种制导弹药舵机放大器仿真测试装置.pdf">一种制导弹药舵机放大器仿真测试装置.pdf</a> </li> <li class="vprdn"> <img alt="北斗卫星枪支动态定位系统.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326502.html" title="北斗卫星枪支动态定位系统.pdf">北斗卫星枪支动态定位系统.pdf</a> </li> <li class="vprdn"> <img alt="制导弹药多路时序状态信号测量电路及其测量方法.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326480.html" title="制导弹药多路时序状态信号测量电路及其测量方法.pdf">制导弹药多路时序状态信号测量电路及其测量方法.pdf</a> </li> <li class="vprdn"> <img alt="RF集成电路的测试方法和系统.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-6421345.html" title="RF集成电路的测试方法和系统.pdf">RF集成电路的测试方法和系统.pdf</a> </li> <li class="vprdn"> <img alt="电力变压器绕组故障模拟试验方法.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-6421258.html" title="电力变压器绕组故障模拟试验方法.pdf">电力变压器绕组故障模拟试验方法.pdf</a> </li> <li class="vprdn"> <img alt="磁传感器.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-6421149.html" title="磁传感器.pdf">磁传感器.pdf</a> </li> <li class="vprdn"> <img alt="一种利用二次反射波测量超声波渡越时间的电路.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-6421130.html" title="一种利用二次反射波测量超声波渡越时间的电路.pdf">一种利用二次反射波测量超声波渡越时间的电路.pdf</a> </li> <li class="vprdn"> <img alt="接收器与通信信号的接收方法.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-6420702.html" title="接收器与通信信号的接收方法.pdf">接收器与通信信号的接收方法.pdf</a> </li></div> </div> <div class="box hot-keywords mt10" id="box3" style="overflow: hidden;width: 268px; border:solid 1px #dedede;"> <div class="boxHd" style="border: none;padding-bottom: 0px;"> <div class="fl keywords-ico mr5 ml13 ico"> </div> <h2 class="fl font-normal font16 font-yahei" style="font-size: 16px; font-weight: 100; margin-left: 0px; margin-top: 3px; font-family: 微软雅黑"> 相关搜索</h2> </div> <div class="hot-keywords-list"> <a target="_blank" href="search.html?p=0&q=%e4%b8%80%e7%a7%8d" class="tag-item ico" title="一种" hidefocus="true"><span class="ico"><em> 一种</em></span></a> <a target="_blank" href="search.html?p=0&q=%e5%9f%ba%e4%ba%8e" class="tag-item ico" title="基于" hidefocus="true"><span class="ico"><em> 基于</em></span></a> <a target="_blank" href="search.html?p=0&q=%e4%bb%bb%e6%84%8f" class="tag-item ico" title="任意" hidefocus="true"><span class="ico"><em> 任意</em></span></a> <a target="_blank" href="search.html?p=0&q=%e5%b9%bf%e8%a7%92" class="tag-item ico" title="广角" hidefocus="true"><span class="ico"><em> 广角</em></span></a> <a target="_blank" href="search.html?p=0&q=%e6%b3%a2%e5%8a%a8" class="tag-item ico" title="波动" hidefocus="true"><span class="ico"><em> 波动</em></span></a> <a target="_blank" href="search.html?p=0&q=%e6%96%b9%e7%a8%8b" class="tag-item ico" title="方程" hidefocus="true"><span class="ico"><em> 方程</em></span></a> <a target="_blank" href="search.html?p=0&q=%e6%88%90%e5%83%8f" class="tag-item ico" title="成像" hidefocus="true"><span class="ico"><em> 成像</em></span></a> <a target="_blank" href="search.html?p=0&q=%e6%96%b9%e6%b3%95" class="tag-item ico" title="方法" hidefocus="true"><span class="ico"><em> 方法</em></span></a> <a target="_blank" href="search.html?p=0&q=%e8%a3%85%e7%bd%ae" class="tag-item ico" title="装置" hidefocus="true"><span class="ico"><em> 装置</em></span></a> <input name="ctl00$Content$hiddenCategoryID" type="hidden" id="Content_hiddenCategoryID" value="117" /> </div> </div> </div> </div> <!--end boxright--> <!--end layout01--> </div> <!--end mainpart--> </div> <div class="aspNetHidden"><input type="hidden" name="__VIEWSTATEGENERATOR" id="__VIEWSTATEGENERATOR" value="541A68AF" /></div></form> <div style="display: none" id="Div3"> </div> <div id="jubao" class="jubao" style="display: none;"> <div class="jubaoTitle"> 举报类型</div> <ul class="vprdn"> <li class="vprdn"> <input type="radio" id="jubaoSel12" name="jubaoSel12" value="12" />&nbsp;&nbsp; <label for="jubaoSel12">广告或垃圾信息</label> </li> <li class="vprdn"> <input type="radio" id="jubaoSel13" name="jubaoSel13" value="13" />&nbsp;&nbsp; <label for="jubaoSel13">色情、淫秽、低俗信息</label> </li> <li class="vprdn"> <input type="radio" id="jubaoSel14" name="jubaoSel14" value="14" />&nbsp;&nbsp; <label for="jubaoSel14">反政府、反人类、反社会等反动信息</label> </li> <li class="vprdn"> <input type="radio" id="jubaoSel15" name="jubaoSel15" value="15" />&nbsp;&nbsp; <label for="jubaoSel15">散布赌博、暴力、凶杀、恐怖或者教唆犯罪等信息</label> </li> <li class="vprdn"> <input type="radio" id="jubaoSel16" name="jubaoSel16" value="16" />&nbsp;&nbsp; <label for="jubaoSel16">侮辱、诽谤等人身攻击信息</label> </li> <li class="vprdn"> <input type="radio" id="jubaoSel17" name="jubaoSel17" value="17" />&nbsp;&nbsp; <label for="jubaoSel17">散布谣言、扰乱社会秩序,破坏社会稳定等信息</label> </li> <li style="text-align: center; height: 50px; margin-top: 10px;"> <input type="hidden" id="hidTypeID" /> <input type="button" onclick="saveJuBao()" value="提 交" />&nbsp;&nbsp;&nbsp;&nbsp;<input type="button" onclick="$('#jubao').hide();" value="取 消" /></li> </ul> </div> <div id="jubaoResponse" class="jubao" style="display: none;"> <div class="jubaoTitle"> 提示</div> <ul class="vprdn"> <li class="vprdn"> <h4 style="text-align: center; color: rgb(83, 136, 35);"> </h4> </li> <li style="text-align: center; height: 50px; margin-top: 10px;"> <input type="button" onclick="$('#jubaoResponse').hide();" value="关 闭" /> </ul> </div></div> <div class="tempdiv" style="display:none;line-height:0px;height:0px; overflow:hidden;"> </div> <!-- JiaThis Button END --> <span id="LabelScript"></span> </div> </div> <!--foot--><div class="bg_100 foot_nav_bg" style=" min-width:1200px;"> <div class="foot_nav"> <a href="//www.4mum.com.cn/h-33.html" target="_blank" rel="nofollow">关于我们</a> - <a href="//www.4mum.com.cn/h-34.html" target="_blank" rel="nofollow">网站声明</a> - <a href="//www.4mum.com.cn/h-35.html" target="_blank" rel="nofollow"> 网站地图</a> - <a href="//www.4mum.com.cn/sitemap.html" target="_blank" rel="nofollow"> 资源地图</a> - <a href="//www.4mum.com.cn/friend.aspx" target="_blank" rel="nofollow">友情链接</a> - <a rel="nofollow">网站客服</a> - <a href="//www.4mum.com.cn/h-93.html" target="_blank" rel="nofollow">联系我们</a> </div></div><div class="bg_100 siteInner_bg" style=" min-width:1200px;"> <div class="siteInner"> <p style="text-align: center; line-height: 1.75em;"><span style="font-size: 14px;"><span style="font-size: 14px; color: rgb(102, 102, 102); font-family: 微软雅黑, Arial, &#39;Times New Roman&#39;; line-height: 20px; text-align: center;">copyright@ 2017-2018 www.4mum.com.cn网站版权所有</span><br style="color: rgb(102, 102, 102); font-family: 微软雅黑, Arial, &#39;Times New Roman&#39;; font-size: 12px; line-height: 20px; text-align: center; white-space: normal; "/><span style="font-size: 14px; color: rgb(102, 102, 102); font-family: 微软雅黑, Arial, &#39;Times New Roman&#39;; line-height: 20px; text-align: center;">经营许可证编号:粤ICP备17046363号-1&nbsp;</span></span> &nbsp;</p><p><br/></p> </div></div><!--foot end--> <!-- 代码部分begin --><div class="QQ_S" style="height: 172px;position:fixed;right: 0px;bottom: 20px; top:auto;"> <div class="Q_top" onclick="HideFoot()"> <span class="signi"></span>收起</div> <div class="Q_botm"> <div class="Q_pic"> <div class="Q_pic_hide"> <a target="_blank" title="在线客服" rel="nofollow"><span class="hide_pic"></span>在线客服</a> </div> </div> <div class="Q_anser"> <div class="Q_anser_hide"><a target="_blank" title="意见反馈"> <span class="hide_pic1"></span>意见反馈 </a> </div> </div> <div class="Q_sign"> <div class="Q_sign_hide" onclick="backToTop();"><a href="javascript:void(0)" title="返回顶部"><span class="hide_pic2"></span>返回顶部 </a></div> </div> </div></div><div class="QQ_S1"> <div class="Q_top1" onclick="ShowFoot()"> <span class="signj"></span>展开</div> <div class="Q_botm1"> <div class="Q_pic1"> <div class="Q_pic1_hide"> <a href="//www.4mum.com.cn/" title="重庆时时彩单双窍门">重庆时时彩单双窍门</a></div> </div> <div class="Q_sign1"> <div class="Q_sign1_hide" onclick="backToTop();"><a href="javascript:void(0)">返回顶部</a></div> </div> </div></div> <!-- 代码部分end --> <li><a href="//www.4mum.com.cn/77d944/69034815.html ">四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 </a> 2019-05-13</li> <li><a href="//www.4mum.com.cn/euy29/3433514.html ">银保监会新规剑指大企业多头融资和过度融资 </a> 2019-05-12</li> <li><a href="//www.4mum.com.cn/a62999/3431813.html ">韩国再提4国联合申办世界杯 中国网友无视:我们自己来 </a> 2019-05-11</li> <li><a href="//www.4mum.com.cn/4kf294/98434412.html ">中国人为什么一定要买房? </a> 2019-05-11</li> <li><a href="//www.4mum.com.cn/33r392/76278011.html ">十九大精神进校园:风正扬帆当有为&nbsp;勇做时代弄潮儿 </a> 2019-05-10</li> <li><a href="//www.4mum.com.cn/1jd106/47417610.html ">粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 </a> 2019-05-09</li> <li><a href="//www.4mum.com.cn/997397/2262689.html ">太原设禁鸣路段 设备在测试中 </a> 2019-05-09</li> <li><a href="//www.4mum.com.cn/br3559/3646538.html ">拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 </a> 2019-05-08</li> <li><a href="//www.4mum.com.cn/9h7915/8581547.html ">“港独”没出路!“梁天琦们”该醒醒了 </a> 2019-05-07</li> <li><a href="//www.4mum.com.cn/3n3511/2345366.html ">陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 </a> 2019-05-06</li> <li><a href="//www.4mum.com.cn/b0059/6332815.html ">人民日报客户端辟谣:“合成军装照”产品请放心使用 </a> 2019-05-05</li> <li><a href="//www.4mum.com.cn/v5j888/7405464.html ">【十九大·理论新视野】为什么要“建设现代化经济体系”?&nbsp;&nbsp; </a> 2019-05-04</li> <li><a href="//www.4mum.com.cn/pl1926/4722623.html ">聚焦2017年乌鲁木齐市老城区改造提升工程 </a> 2019-05-04</li> <li><a href="//www.4mum.com.cn/14g861/9039202.html ">【专家谈】上合组织——构建区域命运共同体的有力实践者 </a> 2019-05-03</li> <li><a href="//www.4mum.com.cn/7jn906/5662951.html ">【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! </a> 2019-05-03</li> <script>(function(){ var src = (document.location.protocol == "") ? "//js.passport.qihucdn.com/11.0.1.js?9ed1f3a8f9c3ff069b7b95c01474c743":"https://jspassport.ssl.qhimg.com/11.0.1.js?9ed1f3a8f9c3ff069b7b95c01474c743"; document.write('<script src="' + src + '" id="sozz"><\/script>'); })(); </script> <script> (function(){ var bp = document.createElement('script'); var curProtocol = window.location.protocol.split(':')[0]; if (curProtocol === 'https') { bp.src = 'https://zz.bdstatic.com/linksubmit/push.js'; } else { bp.src = '//push.zhanzhang.baidu.com/push.js'; } var s = document.getElementsByTagName("script")[0]; s.parentNode.insertBefore(bp, s); })(); </script> <a href="//www.467045.live">红警ol</a> <a href="//www.138381.live">欢乐捕鱼人弹头怎么买</a> <a href="//www.8118743.live">广西快乐十分软件</a> <a href="//www.472012.live">丰禾棋牌官网</a> <a href="//www.968441.live">湖南幸运赛车爱彩人</a> <a href="//www.206500.live">2018年一码三中三</a> <a href="//www.5897214.live">彩票大赢家游戏</a> <a href="//www.967856.live">陕西快乐10分选号规律</a> <a href="//649519.live">决策天机炒股软件</a> <a href="//www.200490.live">福彩喜乐彩</a> <a href="//www.752633.live">四川金7乐二维码</a> <a href="//www.1240419.live">双色球免费预测</a> <a href="//www.5175569.live">黑龙江22选5走势图2008年</a> <a href="//www.1987656.live">甘肃十一选五开奖结</a> <a href="//www.188197.live">怎么买新时时彩</a> <a href="//www.259520.live">青海快3第26期开奖结果查询</a> </body></html>