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

    重庆时时彩杀号: 一种基于叠前地震数据的裂缝储层快速预测方法.pdf

    关 键 词:
    一种 基于 地震 数据 裂缝 快速 预测 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201611193625.X

    申请日:

    2016.12.21

    公开号:

    CN106772601A

    公开日:

    2017.05.31

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01V 1/30申请日:20161221|||公开
    IPC分类号: G01V1/30 主分类号: G01V1/30
    申请人: 成都理工大学
    发明人: 熊晓军; 简世凯; 李翔; 刘阳; 龚思宇; 崔泽飞
    地址: 610059 四川省成都市二仙桥东三路1号
    优先权:
    专利代理机构: 北京轻创知识产权代理有限公司 11212 代理人: 谈杰
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201611193625.X

    授权公告号:

    ||||||

    法律状态公告日:

    2018.10.23|||2017.06.23|||2017.05.31

    法律状态类型:

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

    摘要

    本发明公开了一种基于叠前地震数据的裂缝储层快速预测方法,根据叠前地震资料的品质,采用去噪??拉平??切除的组合优化对地震资料进行处理,在优化处理后的道集上,拟合某一时间点各道的振幅值,将该时间点各道的拟合值减去看作均匀介质的各道近道集振幅值的平均值,得到各道的差值;求其标准差作为异常值,可用该异常值有效地表征裂缝储层发育强度;以采样率为步长,移动时间点,得到整个地震数据的异常值数据体。本发明快捷有效;不依赖于各向异性椭圆拟合的裂缝储层预测分析技术,解决宽窄方位角地震数据裂缝储层预测的问题。

    权利要求书

    1.一种基于叠前地震数据的裂缝储层快速预测方法,其特征在于,所述基于叠前地震
    数据的裂缝储层快速预测方法根据叠前地震资料的品质,采用去噪-拉平-切除的组合优化
    对地震资料进行处理,在优化处理后的道集上,拟合某一时间点各道的振幅值,将该时间点
    各道的拟合值减去看作均匀介质的各道近道集振幅值的平均值,得到各道的差值;求其标
    准差作为异常值,可用该异常值有效地表征裂缝储层发育强度;异常值看作是裂缝和流体
    的综合影响的结果;以采样率为步长,移动时间点,得到整个地震数据的异常值数据体。
    2.如权利要求1所述的基于叠前地震数据的裂缝储层快速预测方法,其特征在于,所述
    采用去噪-拉平-切除的组合优化对地震资料进行处理具体包括:
    步骤一,采用四阶多项式拟合进行去噪处理;
    步骤二,采用非地表一致性剩余静校正实现道集拉平处理;
    步骤三,叠前远道集在经过去噪、拉平处理后,若仍然存在杂乱,道集数据缺损、噪音干
    扰严重现象,需要对数据采取切除30°以外的道集。
    3.如权利要求2所述的基于叠前地震数据的裂缝储层快速预测方法,其特征在于,所述
    采用四阶多项式拟合进行去噪处理包括:
    (1)基于AVO曲线拟合,采用四阶多项式拟合方法拟合叠前道集某一时刻的数据点振幅
    值,得到拟合曲线;
    (2)选择初始参数,将拟合曲线分别向上向下移动,对比同一偏移距的数据,去除高于
    向上移动值和低于向下移动值的采样点值,用相同偏移距的多项式拟合值代替;
    (3)检验拟合的AVO曲线,是否与正演道集AVO曲线趋势一致,一致则输出去噪道集,否
    则调整参数值大小,重复(1)和(2),输出去噪后道集。
    4.如权利要求2所述的基于叠前地震数据的裂缝储层快速预测方法,其特征在于,所述
    采用非地表一致性剩余静校正实现道集拉平处理包括:
    (1)将一定偏移距内的道集进行叠加作为初始道S(t);在一定时窗t内,分别计算各道
    与初始道的相关系数R(t,I),采用如下公式求取,其中t1、t2为时窗的起止时间,X(t-τ,I)为
    输入道,I为道号,τ=0,1,2……n,n为互相关函数的长度;
    <mrow> <mi>R</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>,</mo> <mi>I</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> </mrow> <msub> <mi>t</mi> <mn>2</mn> </msub> </munderover> <mi>S</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>X</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&tau;</mi> <mo>,</mo> <mi>I</mi> <mo>)</mo> </mrow> </mrow>
    (2)选择相关系数最大的一道作为模型道,从某道开始依次在指定可移动时窗T范围
    内,以某个时窗t滑动求取与相邻道的相关系数,所用公式与(1)中所述相同;
    (3)将相关系数最大时所滑动的时窗作为当前道的剩余时差,在指定时间段T内,各道
    减去看作剩余时差的时间之后,输出拉平道集。
    5.如权利要求2所述的基于叠前地震数据的裂缝储层快速预测方法,其特征在于,所述
    统计计算AVO异常值的方法包括:
    第一步,对优化处理之后的叠前道集一条测线同一时间点上的振幅值进行多项式拟
    合,采用最小二乘法原理来实现;
    第二步,根据均匀介质与各向异性介质在入射角0°-6°时其值基本相同,首先根据第一
    步得到某一时刻各道的拟合值,同时根据深度和偏移距的关系计算近道集角度,并求取各
    道入射角0°-6°的近道集振幅值的平均值,将道集各道的拟合值减去看作均匀介质的各道近
    道集振幅值的平均值,得到各道的差值;
    第三步,利用如下公式:
    <mrow> <mi>S</mi> <mo>=</mo> <msqrt> <mrow> <mfrac> <mn>1</mn> <mi>M</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>b</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>,</mo> <mi>M</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>3</mn> <mo>,</mo> <mn>4.......</mn> <mo>;</mo> </mrow>
    其中bi为原始值,bii为拟合值,M为叠前道集某一道的总道数,i为道号,S为某道的标准
    差;统计各道的差值,计算标准差S作为各道此时间点的异常值,异常值看作是裂缝和流体
    综合影响的结果,用该异常值有效地表征裂缝储层发育强度;
    第四步,以采样率为步长,移动时间点,循环第一步-第三步,得到整条线的地震数据的
    异常值;
    第五步,改变测线,循环第一步-第四步,得到整个地震数据的异常数据体,通过剖面和
    平面成图,得到裂缝储层预测的图形效果,实现裂缝储层预测。

    说明书

    一种基于叠前地震数据的裂缝储层快速预测方法

    技术领域

    本发明属于油气地球物理技术领域,尤其涉及一种基于叠前地震数据的裂缝储层
    快速预测方法。

    背景技术

    根据采集的地震数据,现有的基于叠前地震数据裂缝预测技术,可分为针对宽方
    位角叠前地震数据的裂缝预测方法和针对窄方位角叠前地震数据的裂缝预测方法两类:
    (1)针对宽方位角叠前地震裂缝预测,国内外学者大多采用的是基于P波各向异性地震属性
    椭圆拟合技术,如Grechka和Tsvankin等(1999年)提出对水平层状介质中关于横向速度变
    化的动校正时差椭圆拟合,来预测随深度变化的裂缝倾向;曲寿利等(2001年)提出了根据
    阻抗随方位角变化来检测裂缝的方法;Gray(2004年)提出AVO梯度随全方位角变化的方法
    来预测裂缝(其过程与拟合椭圆相近);孙炜等(2014年)提出了一种改进的方位各向异性裂
    缝预测方法。以上方法实现过程均较复杂,且在处理窄方位角地震数据的裂缝预测时,拟合
    椭圆受到限制;(2)针对窄方位角地震资料的裂缝预测,国内学者也提出了一些技术方法,
    如田立新等(2010年)通过远近偏移距属性差异地震正演和方位地震属性差异的地震正演,
    模拟得到针对窄方位地震采集区块的叠前裂缝储层预测方法;苏世龙等(2015年)通过采取
    限炮检距接收及一些特殊的数据规则化技术改善窄方位资料的性质从而进行叠前裂缝预
    测。这些针对窄方位角地震数据裂缝预测的技术方法,采取分方位角处理,过程复杂,理论
    不够完善。此外,多波技术也用于裂缝预测中,多波资料相对于裂缝的敏感程度更强,但是
    采集多波分量数据的施工成本太高,资料处理难度大,在实际生产中多波采集较少,因此多
    波技术受到很大的限制,不利于普遍适用。

    综上所述,现有的基于叠前地震数据裂缝预测方法存在实现过程复杂,拟合椭圆
    受到限制,施工成本太高,资料处理难度大。

    发明内容

    本发明的目的在于提供一种基于叠前地震数据的裂缝储层快速预测方法,旨在解
    决现有的基于叠前地震数据裂缝预测方法存在实现过程复杂,拟合椭圆受到限制,施工成
    本太高,资料处理难度大的问题。

    本发明是这样实现的,一种基于叠前地震数据的裂缝储层快速预测方法,所述基
    于叠前地震数据的裂缝储层快速预测方法根据叠前地震资料的品质,采用去噪-拉平-切除
    的组合优化对地震资料进行处理,在优化处理后的道集上,拟合某一时间点各道的振幅值,
    将该时间点各道的拟合值减去看作均匀介质的各道近道集振幅值的平均值,得到各道的差
    值;求其标准差作为异常值,可用该异常值有效地表征裂缝储层发育强度;异常值看作是裂
    缝和流体的综合影响的结果;以采样率为步长,移动时间点,得到整个地震数据的异常值数
    据体。

    进一步,所述采用去噪-拉平-切除的组合优化对地震资料进行处理具体包括:

    步骤一,采用四阶多项式拟合进行去噪处理;

    步骤二,采用非地表一致性剩余静校正实现道集拉平处理;

    步骤三,叠前远道集在经过去噪、拉平处理后,若仍然存在杂乱,道集数据缺损、噪
    音干扰严重等现象,需要对数据采取切除30°以外的道集。

    进一步,所述采用四阶多项式拟合进行去噪处理包括:

    (1)基于AVO曲线拟合,采用四阶多项式拟合方法拟合叠前道集某一时刻的数据点
    振幅值,得到拟合曲线;

    (2)选择初始参数,将拟合曲线分别向上向下移动,对比同一偏移距的数据,去除
    高于向上移动值和低于向下移动值的采样点值,用相同偏移距的多项式拟合值代替;

    (3)检验拟合的AVO曲线,是否与正演道集AVO曲线趋势一致,一致则输出去噪道
    集,否则调整参数值大小,重复(1)和(2),输出去噪后道集。

    进一步,所述采用非地表一致性剩余静校正实现道集拉平处理包括:

    (1)将一定偏移距内的道集进行叠加作为初始道S(t);在一定时窗t内,分别计算
    各道与初始道的相关系数R(t,I),采用如下公式求取,其中t1、t2为时窗的起止时间,X(t-τ,
    I)为输入道,I为道号,τ=0,1,2……n,n为互相关函数的长度;


    (2)选择相关系数最大的一道作为模型道,从某道开始依次在指定可移动时窗T范
    围内,以某个时窗t滑动求取与相邻道的相关系数,所用公式与(1)中所述相同;

    (3)将相关系数最大时所滑动的时窗作为当前道的剩余时差,在指定时间段T内,
    各道减去看作剩余时差的时间之后,输出拉平道集。

    进一步,所述统计计算AVO异常值的方法包括:

    第一步,对优化处理之后的叠前道集一条测线同一时间点上的振幅值进行多项式
    拟合,采用最小二乘法原理来实现;

    第二步,根据均匀介质与各向异性介质在入射角0°-6°时其值基本相同,首先根据
    第一步得到某一时刻各道(这里各道指的是叠前道集包括CRP、CDP、CMP)的拟合值,同时根
    据深度和偏移距的关系计算近道集角度,并求取各道近道集(入射角0°-6°)振幅值的平均
    值,将道集各道的拟合值减去看作均匀介质的各道近道集振幅值的平均值,得到各道的差
    值;

    第三步,利用如下公式:


    其中bi为原始值,bii为拟合值,M为叠前道集某一道的总道数,i为道号,S为某道的
    标准差。统计各道的差值,计算标准差S作为各道此时间点的异常值,异常值可看作是裂缝
    和流体综合影响的结果,可用该异常值有效地表征裂缝储层发育强度;

    第四步,以采样率为步长,移动时间点,循环第一步-第三步,得到整条线的地震数
    据的异常值。

    第五步,改变测线,循环第一步-第四步,得到整个地震数据的异常数据体,通过剖
    面和平面成图,得到裂缝储层预测的图形效果,实现裂缝储层预测。

    进一步,多项式拟合,最小二乘法原理包括:

    1)在某时刻,若观测到的数据为一个入射角xi对应一个振幅值yi如下表所示:

    xi
    x1x2……xn
    yi
    y1y2…….yn

    2)假设拟合函数不是线性函数,用如下的公式来进行拟合:

    y=a0+a1x+a2x2+a3x3+.......+anxn,n=0,1,2,3........;

    3)设拟合值和真实值之间的残差平方和为δi:


    4)当残差平方最小时,拟合曲线与真实值之间最为接近。所以,其关键问题就是求
    解多项式的系数。转化问题,关于多项式系数的函数,为得到最小值则求取其导数,令
    经过一系列的推导可以得到如下公式:


    求出ai,然后就可以进行拟合曲线。

    本发明提供的基于叠前地震数据的裂缝储层快速预测方法,适用范围较广、不受
    宽窄方位角地震资料的限制、实现过程快捷;基于各向异性理论,考虑在均匀介质(不含裂
    缝和流体的介质)中,根据Zoeppritz方程及其简化式可知随入射角的变化,振幅大小基本
    上不变;而在各向异性介质(含裂缝和流体)中,振幅随入射角的变化而变化;并验证分析两
    者的AVO曲线随入射角变化的特征,得知裂缝方向无论如何展布,AVO都会表现一定的差异,
    常规的裂缝预测方法基本忽视了这一点;本发明将这种差异定性的作为裂缝储层预测的基
    础,综合考虑了裂缝及流体的综合影响;快速有效地解决宽窄方位角地震数据裂缝储层预
    测的问题;是在优化处理后的道集上,拟合某一时间点各道的振幅值,将该时间点各道的拟
    合值减去看作均匀介质的各道近道集(约0°-6°之间)振幅值的平均值,得到各道的差值;求
    其标准差作为异常值,可用该异常值有效地表征裂缝储层发育强度。如先计算1000ms处的
    异常值之后,若采样率为2ms,则以采样率为步长,移动时间点到1002ms,计算得到异常值,
    依次循环到整个记录时长,得到整个时长的异常值。

    本发明对优化处理之后的叠前道集同一时间点上的振幅值进行多项式拟合;将道
    集各道的拟合值减去看作均匀介质的各道近道集(约0°-6°之间)振幅值的平均值,得到各
    道的差值;统计各道的差值计算标准差作为各道此时间点的异常值,此方法相对于常规的
    基于P波各向异性地震属性椭圆拟合技术的计算量少,计算效率高;且不依赖于各向异性椭
    圆拟合的裂缝储层预测分析技术,解决宽窄方位角地震数据裂缝储层预测的问题。假设优
    化处理步骤一样,基于椭圆拟合的方法,假设也用异常值来拟合椭圆,需要全方位角度的计
    算,还要进行拟合椭圆,来预测裂缝,而本发明的方法,求得异常值,即可作为裂缝预测的依
    据。

    本发明的方法对某工区进行实际资料(CRP道集)分析,在井的周围都有裂缝的发
    育,与实际钻井吻合(应用本发明裂缝储层预测的符合率达到83%);有效地验证了本发明
    的实用性。

    附图说明

    图1是本发明实施例提供的基于叠前地震数据的裂缝储层快速预测方法流程图。

    具体实施方式

    为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明
    进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于
    限定本发明。

    下面结合附图对本发明的应用原理作详细的描述。

    如图1所示,本发明实施例提供的基于叠前地震数据的裂缝储层快速预测方法包
    括以下步骤:

    S101:根据叠前地震资料的品质,首先采用“去噪-拉平-切除”的组合优化技术对
    地震资料进行处理;

    S102:对优化处理之后的叠前道集同一时间点上的振幅值进行二次以上的多项式
    拟合;

    S103:根据均匀介质与各向异性介质在入射角非常小(约0°-6°之间)时其值基本
    相同,将道集各道的拟合值减去看作均匀介质的各道近道集(入射角约0°-6°之间)振幅值
    的平均值,得到各道的差值;

    S104:统计各道拟合值与原始值的标准差作为各道此时间点的异常值,可用该异
    常值有效地表征裂缝储层发育强度;

    S105:以采样率为步长,移动时间点,循环步骤得到整个地震数据的异常值数据
    体。

    下面结合具体实施例对本发明的应用原理作进一步的描述。

    本发明实施例提供的基于叠前地震数据的裂缝储层快速预测方法,具体实施步骤
    如下:

    (1)叠前道集优化处理:

    根据叠前地震资料的品质,首先采用“去噪-拉平-切除”的一个组合优化技术对地
    震资料进行处理,其具体步骤如下:

    步骤1:采用四阶多项式拟合进行去噪处理:

    ①基于AVO曲线拟合,通过四阶多项式拟合叠前道集某一时刻的数据点,得到拟合
    曲线;

    ②选择初始参数,将拟合曲线上下移动,对比同一时刻的数据采样点值,将上下移
    动值之外的值去除,用拟合值代替;

    ③检验拟合的AVO曲线,是否与正演道集AVO曲线趋势一致,一致则输出去噪道集,
    否则调整参数值大小,重复①和②,输出去噪后道集。

    步骤2:采用非地表一致性剩余静校正实现道集拉平处理:

    ①将一定偏移距内的道集进行叠加作为初始道;在一定时窗内,分别计算各道与
    初始道的相关系数,选择相关系数最大的一道作为模型道;

    ②从某道开始依次在指定可移动时窗范围内们,以时窗t滑动求取与相邻道的相
    关系数;

    ③将相关系数最大的滑动数作为当前的剩余时差,在指定时间段T内,各道减去剩
    余时差之后,输出拉平道集。

    步骤3:叠前远道集在经过去噪、拉平处理后,若仍然存在杂乱,道集数据缺损、噪
    音干扰严重等现象,需要对数据采取切除30°以外的道集,这样能保证叠前道集品质及有效
    信息,且减少了后续分析的时间。

    若不存在上述问题,则省略此步骤。

    (2)统计计算AVO异常值

    步骤1:对优化处理之后的叠前道集同一时间点上的振幅值进行多项式拟合,采用
    最小二乘法原理来实现;

    步骤2:根据均匀介质与各向异性介质在入射角非常小(约0°-6°之间)时其值基本
    相同,将道集各道的拟合值减去看作均匀介质的各道近道集(入射角约0°-6°之间)振幅值
    的平均值,得到各道的差值;

    步骤3:利用如下公式(1):

    其中M=1,2,3,4.......(1)

    (bi为原始值,bii为拟合值)统计各道的差值,计算标准差作为各道此时间点的异
    常值,异常值可看作是裂缝和流体综合影响的结果,可用该异常值有效地表征裂缝储层发
    育强度;

    步骤4:以采样率为步长,移动时间点,循环步骤1、步骤2、步骤3,得到整个地震数
    据的异常值数据体,通过剖面和平面成图,得到裂缝储层预测的图形效果,从而实现裂缝储
    层预测。

    进一步,多项式拟合,最小二乘法原理包括:

    1)在某时刻,若观测到的数据为一个入射角xi对应一个振幅值yi如下表所示:

    xi
    x1x2……xn
    yi
    y1y2…….yn

    2)假设拟合函数不是线性函数,用如下的公式来进行拟合:

    y=a0+a1x+a2x2+a3x3+.......+anxn,n=0,1,2,3........;

    3)设拟合值和真实值之间的残差平方和为δi:


    4)当残差平方最小时,拟合曲线与真实值之间最为接近。所以,其关键问题就是求
    解多项式的系数。转化问题,关于多项式系数的函数,为得到最小值则求取其导数,令
    经过一系列的推导可以得到如下公式:


    求出ai,然后就可以进行拟合曲线。

    本发明对优化处理之后的叠前道集同一时间点上的振幅值进行多项式拟合;将道
    集各道的拟合值减去看作均匀介质的各道近道集(约0°-6°之间)振幅值的平均值,得到各
    道的差值;统计各道的差值计算标准差作为各道此时间点的异常值,异常值可看作是裂缝
    和流体的综合影响的结果,此方法相对于常规的基于P波各向异性地震属性椭圆拟合技术
    的计算量少,计算效率高;且不依赖于各向异性椭圆拟合的裂缝储层预测分析技术,解决宽
    窄方位角地震数据裂缝储层预测的问题。假设优化处理步骤一样,基于椭圆拟合的方法,假
    设也用异常值来拟合椭圆,需要全方位角度的计算,还要进行拟合椭圆,来预测裂缝,而本
    发明的方法,求得异常值,即可作为裂缝预测的依据。

    本发明的方法对某工区进行实际资料(CRP道集)分析,在井的周围都有裂缝的发
    育,与实际钻井吻合(应用本发明裂缝储层预测的符合率达到83%);有效地验证了本发明
    的实用性。

    以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精
    神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的?;し段е?。

    关于本文
    本文标题:一种基于叠前地震数据的裂缝储层快速预测方法.pdf
    链接地址://www.4mum.com.cn/p-6001522.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
  • 新时时跟号玩法 pk10人工精准计划软件手机 赛车pk拾免费软件 大乐透常用复式投注及中奖金额表 冠亚和值计划软件下载 北京pk10 稳赚 20选5稳赚不赔技巧 北赛车pk10官网 全网36码 专业版 北单上下单双玩法 后三组六9码稳赚方案 幸运pk10快艇在线直播 北京pk10赛车开结果 幸运飞艇计划软件 app 快3稳赚