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

    重庆时时彩我怎么老输: 基于同位多相协同克里金的地震弹性参数相控建模方法.pdf

    关 键 词:
    基于 同位 多相 协同 克里金 地震 弹性 参数 建模 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201710098888.0

    申请日:

    2017.02.23

    公开号:

    CN106772587A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G01V 1/28申请日:20170223|||公开
    IPC分类号: G01V1/28 主分类号: G01V1/28
    申请人: 河海大学
    发明人: 张宏兵; 韩飞龙; 尚作萍; 郭强; 曹呈浩
    地址: 211100 江苏省南京市江宁开发区佛城西路8号
    优先权:
    专利代理机构: 南京经纬专利商标代理有限公司 32200 代理人: 徐莹
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201710098888.0

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

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

    摘要

    本发明公开了基于同位多相协同克里金的地震弹性参数相控建模方法,包括步骤:在目标层内采用等比例网格剖分方法对建模网格进行剖分,形成若干个平面,作为弹性参数建模的目标位置;分别提取测井弹性参数、地震属性参数和沉积相信息作为建模数据,选择主变量,选择第一及第二协同变量;将所得主变量、第一协同变量及第二协同变量作为计算参数,利用同位多相协同克里金方法进行插值计算,获得平面上所有待估计点的参数值作为多源数据融合建模结果;使用各向异性扩散法进行概化处理,得到概化处理后的地震弹性参数建模结果。本发明能够消除计算异常点和边界噪点,提高了多信息参数综合建模的精度并融入了更多真实的地质信息,具有较好适用性。

    权利要求书

    1.基于同位多相协同克里金的地震弹性参数相控建模方法,其特征在于,包括以下步
    骤:
    步骤S1、在目标层内采用等比例网格剖分方法对建模网格进行剖分,在相同比例厚度
    的位置上形成若干个平面,作为弹性参数建模的目标位置;
    步骤S2、分别提取平面上已知点测井弹性参数以及分布于整个区域的地震属性参数和
    沉积相信息作为建模数据,即选择测井参数的低频参数背景值作为作为主变量,选择地震
    属性参数作为第一协同变量,及利用各向异性扩散法对沉积相信息进行概化处理得到沉积
    相数字模型作为第二协同变量;
    步骤S3、将步骤S2所得主变量、第一协同变量及第二协同变量作为计算参数,利用同位
    多相协同克里金方法进行插值计算,获得平面上所有待估计点的参数值作为弹性参数的多
    源数据融合建模结果;
    步骤S4、使用各向异性扩散法对步骤S3所得多源数据融合建模结果进行概化处理,得
    到概化处理后的地震弹性参数建模结果。
    2.根据权利要求1所述的基于同位多相协同克里金的地震弹性参数相控建模方法,其
    特征在于,所述步骤S2中选择测井参数的低频参数背景值作为主变量,具体包括:
    使用多项式拟合方法对提取的测井弹性参数进行平滑处理,根据不同的多项式拟合次
    数得到若干组平滑后的低频参数背景值;
    选择大于预设低频参数背景值变化范围的一组低频参数背景值作为主变量。
    3.根据权利要求1所述的基于同位多相协同克里金的地震弹性参数相控建模方法,其
    特征在于,所述步骤S2中,根据主变量与地震属性参数间的相关系数选择地震属性参数作
    为第一协同变量。
    4.根据权利要求1所述的基于同位多相协同克里金的地震弹性参数相控建模方法,其
    特征在于,所述步骤S2中利用各向异性扩散法对沉积相信息进行概化处理,采用公式:
    <mrow> <msubsup> <mi>Z</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> <mrow> <mo>(</mo> <mi>m</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>AZ</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> </msubsup> <mo>+</mo> <mfrac> <mi>B</mi> <mn>8</mn> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>,</mo> <mi>l</mi> <mo>=</mo> <mo>-</mo> <mn>1</mn> </mrow> <mn>1</mn> </munderover> <msub> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mi>l</mi> </mrow> </msub> <mo>&dtri;</mo> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mi>l</mi> </mrow> </msub> </mrow>
    其中,Z为待处理点参数值,▽Z为各方向临近点梯度值,A、B、C为各方向的权重系数,i、
    j、k为所处理点的相对位置,m为点的迭代次数。
    5.根据权利要求1所述的基于同位多相协同克里金的地震弹性参数相控建模方法,其
    特征在于,所述步骤S3中利用同位多相协同克里金方法插值进行计算,具体包括:
    步骤S31、利用同位协同多相克里金方法,建立计算平面上待估计点的参数值模型;
    步骤S32、根据无偏最优条件得到克里金方程组,及将克里金方程组转换为矩阵;
    步骤S33、计算不同位置待估计点的各参数间的协方差和互协方差,以得到步骤S32所
    得矩阵中所有变量值;
    步骤S34、通过求解步骤S32所得矩阵得到各变量的权重系数,将其代入步骤S31所建立
    待估计点的参数值模型计算得到待估计点的参数值。
    6.根据权利要求1所述的基于同位多相协同克里金的地震弹性参数相控建模方法,其
    特征在于,所述步骤S3中获得弹性参数的多源数据融合建模结果为:
    <mrow> <msup> <mi>Z</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> <mi>Z</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <msub> <mi>Y</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&beta;</mi> <mn>2</mn> </msub> <msub> <mi>Y</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow>
    其中,待估计点位置的参数估计值为Z*(u0),主变量的测井弹性参数为Z(ui),u为参数
    点位置,作为第一协同变量的地震属性参数为Y1,作为第二协同变量的沉积相数字模型为
    Y2,权重系数分别为αi、β1和β2。
    7.根据权利要求1所述的基于同位多相协同克里金的地震弹性参数相控建模方法,其
    特征在于,所述步骤S4中,使用各向异性扩散法概化处理具体包括:
    步骤S41、设定待估计点的参数变化阈值,将大于或小于设定参数变化阈值的待估计点
    参数值视为异常点去除;
    步骤S42、利用各向异性扩散法概化处理剩余待估计点的参数值,通过反复迭代得到概
    化处理后的地震弹性参数建模结果。

    说明书

    基于同位多相协同克里金的地震弹性参数相控建模方法

    技术领域

    本发明涉及一种基于同位多相协同克里金的地震弹性参数相控建模方法,属于地
    球物理中参数建模的技术领域。

    背景技术

    地震弹性参数建模是一个多信息综合分析的过程,在复杂的地质环境下,保持波
    速建模的稳定性和高精度非常困难。由于整个地质系统的变化受到自然环境、层位埋深及
    土壤和岩层厚度的影响,同时地质构造的复杂性会导致速度模型呈现极大的非均质性和不
    确定性。因此需要一种能够综合多种地质情况的建模方法,来表现这种非均质性并降低不
    确定性。要得到精细的地震反演或深度偏移结果,同样依赖于地震波速建模构建的初始模
    型。尽管线性反演方法可以利用井位处的测井信息构建合适的边界条件,使其与地震信息
    相结合,但初始模型的准确性会在很大程度上影响反演的精度。地震波走时速度分析方法
    通过对初始模型的不断修改使其收敛,但对强横向变化的适应性较差,并且在计算过程中
    容易产生较大的累积误差。

    海上碳酸盐岩区域地质条件复杂,非均质性强,作为主变量的测井数据又相对稀
    少,为了提高弹性参数建模精度和降低不确定性,需要在有限的观测数据中加入有效的约
    束条件。研究区域所能得到的约束条件往往是地震属性参数和沉积相信息,而这些信息都
    能在一定程度上反应主变量的变化特征,故需要选取合适的属性参数以约束建模过程?;?br />于协同克里金的弹性参数建模方法除主变量外,同时使用了一个次级变量作为约束条件,
    有效表达了地质条件的先验信息,降低了建模的不确定性。由于使用协同克里金法所得到
    的估计值综合表达了周围样点的参数值,样点权重系数不仅与主变量空间协方差函数有
    关,还受到一种协同变量的约束控制,因此在一定程度上提高了建模的稳定性。显然,当继
    续增加约束时,带有多个次级变量的协同克里金法能够进一步提高建模的准确性,但约束
    条件的提取与使用仍有很多困难。

    使用常规地震属性参数对不同主变量进行约束,基于协同克里金方法所得到的建
    模结果往往不具有清晰的边界且存在大量噪点,不能真实反应地下介质弹性参数的分布情
    况,特别是在介质属性横向连续性较差的区域,难以分辨真实的局部扰动和由于计算所产
    生的微小误差,对后续的反演及相关工作造成困难。

    发明内容

    本发明所要解决的技术问题在于克服现有技术的不足,提供一种基于同位多相协
    同克里金的地震弹性参数相控建模方法,解决基于协同克里金方法所得到的建模结果往往
    不具有清晰的边界且存在大量噪点,不能真实反应地下介质弹性参数的分布情况的问题,
    可以消除计算异常点和边界噪点,提高了多信息参数综合建模的精度并融入了更多真实的
    地质信息,具有较好适用性。

    本发明具体采用以下技术方案解决上述技术问题:

    基于同位多相协同克里金的地震弹性参数相控建模方法,包括以下步骤:

    步骤S1、在目标层内采用等比例网格剖分方法对建模网格进行剖分,在相同比例
    厚度的位置上形成若干个平面,作为弹性参数建模的目标位置;

    步骤S2、分别提取平面上已知点的测井弹性参数以及分布于整个区域的地震属性
    参数和沉积相信息作为建模数据,即选择测井参数的低频参数背景值作为作为主变量,选
    择地震属性参数作为第一协同变量,及利用各向异性扩散法对沉积相信息进行概化处理得
    到沉积相数字模型作为第二协同变量;

    步骤S3、将步骤S2所得主变量、第一协同变量及第二协同变量作为计算参数,利用
    同位多相协同克里金方法进行插值计算,获得平面上所有待估计点的参数值作为弹性参数
    的多源数据融合建模结果;

    步骤S4、使用各向异性扩散法对步骤S3所得多源数据融合建模结果进行概化处
    理,得到概化处理后的地震弹性参数建模结果。

    进一步地,作为本发明的一种优选技术方案:所述步骤S2中选择测井参数的低频
    参数背景值作为作为主变量,具体包括:

    使用多项式拟合方法对透气的测井弹性参数进行平滑处理,根据不同的多项式拟
    合次数得到若干组平滑后的低频参数背景值;

    选择大于预设低频参数背景值变化范围的一组低频参数背景值作为主变量。

    进一步地,作为本发明的一种优选技术方案:所述步骤S2中,根据主变量与地震属
    性参数间的相关系数选择地震属性参数作为第一协同变量。

    进一步地,作为本发明的一种优选技术方案:所述步骤S2中利用各向异性扩散法
    对沉积相信息进行概化处理,采用公式:


    其中,Z为待处理点参数值,为各方向临近点梯度值,A、B、C为各方向的权重系
    数。

    进一步地,作为本发明的一种优选技术方案:所述步骤S3中利用同位多相协同克
    里金方法插值进行计算,具体包括:

    步骤S31、利用同位协同多相克里金方法,建立计算平面上待估计点的参数值模
    型;

    步骤S32、根据无偏最优条件得到克里金方程组,及将克里金方程组转换为矩阵;

    步骤S33、计算不同位置待估计点参数间的协方差和互协方差,以得到步骤S32所
    得矩阵中所有变量值;

    步骤S34、通过求解步骤S32所得矩阵得到各变量的权重系数,将其代入步骤S31所
    建立待估计点的参数值模型计算得到待估计点的参数值。

    进一步地,作为本发明的一种优选技术方案:所述步骤S3中获得弹性参数的多源
    数据融合建模结果为:


    其中,待估计点位置的参数估计值为Z*(u0),主变量的测井弹性参数为Z(ui),u为
    参数点位置,作为第一协同变量的地震属性参数为Y1,作为第二协同变量的沉积相数字模
    型为Y2,权重系数分别为αi、β1和β2。

    进一步地,作为本发明的一种优选技术方案:所述步骤S4中,使用各向异性扩散法
    概化处理具体包括:

    步骤S41、设定待估计点的参数变化阈值,将大于或小于设定的参数变化阈值的待
    估计点参数值视为异常点去除;

    步骤S42、利用各向异性扩散法概化处理剩余待估计点参数值,通过反复迭代得到
    概化处理后的地震弹性参数建模结果。

    本发明采用上述技术方案,能产生如下技术效果:

    本发明通过对多个变量参数值的提取和基于同位协同克里金方法对多参数进行
    的弹性参数融合建模,获得了建模过程中合理优化参数的方法,实现了利用沉积相数字模
    型作为协同变量的相控建模,提高了地球物理领域弹性参数的建模精度,使用各向异性扩
    散法去除计算噪点的同时保留了地质体所包含的局部信息。最终实现了多参数融合的相控
    建模,提高了基于常规资料的弹性参数建模的准确性和完整性??梢杂行扑阋斐5?br />和边界噪点,提高了多信息参数综合建模的精度并融入了更多真实的地质信息,是具有较
    好适用性的基于同位多相协同克里金的弹性参数相控建模方法。

    附图说明

    图1为本发明的基于同位多相协同克里金的弹性参数相控建模方法流程图。

    图2为本发明的等比例网格剖分示意图。

    图3为本发明示例所使用的测井原始参数及平滑后的低频参数背景值示意图。

    图4为本发明示例所使用的沉积相数字模型示意图。

    图5为本发明示例所使用的基于同位多相协同克里金的弹性参数相控建模方法对
    测井横波波速参数进行建模所得到的结果示意图。

    具体实施方式

    下面结合说明书附图对本发明的实施方式进行描述。

    如图1所示,本发明提出的基于同位多相协同克里金的弹性参数相控建模方法,主
    要包含建模网格剖分、多源信息提取、基于同位多相克里金的插值计算和基于各向异性扩
    散法的概化处理四大部分,具体如下:

    步骤S1、在目标层内采用等比例网格剖分方法对建模网格进行剖分,如图2所示,
    在相同比例厚度的位置上形成多个平面,作为弹性参数建模的目标位置;若目标位置没有
    数据点则使用竖直方向的最邻近点作为目标位置的值。其中,平面上存在若干个已知点和
    待估计点,所述已知点可以直接获得测井弹性参数。

    步骤S2、多源信息提取,分别提取平面上所有已知点的测井弹性参数以及分布于
    整个区域的地震属性参数和沉积相信息作为建模数据,即选择测井参数的低频参数背景值
    作为作为主变量,选择地震属性参数作为第一协同变量,及利用各向异性扩散法对沉积相
    信息进行概化处理得到沉积相数字模型作为第二协同变量;分析数据体相关性并对所提取
    数据进行初步处理,实现多源数据融合建模流程,提高地球物理弹性参数建模精度。该过程
    包括以下步骤:

    步骤S21、提取已知点的所有测井弹性参数低频背景值,尝试不同的拟合次数并结
    合建模要求与测井原始数据进行比对,获得能够满足建模精度要求的数据频率,如图3所
    示。

    其具体为:

    步骤S211、使用多项式拟合方法对已知点的所有测井弹性参数数据进行平滑处
    理,根据不同的多项式拟合次数得到多组平滑后的低频参数背景值。

    步骤S212、绘制低频参数背景值曲线图,并将预设低频参数背景值范围设置为40m
    范围变化,由于建模结果所需频率实际低于测井采集的原始数据,因此根据弹性参数变化
    特点选择能够反映大于40m范围变化的一组低频背景值作为主变量,在深度域数据体中使
    低频背景值能够反映出大段岩层并且在薄层区域反映出一定的变化趋势。

    步骤S22、提取地震属性参数,通过统计方法获得能够在侧面反映主变量变化的参
    数类型,计算整个区域不同地震属性数据体与测井低频背景值的相关系数,根据计算结果
    选择普遍适用的地震属性作为第一协同变量,实现建模参数的优化。其包括以下步骤:

    步骤S221、绘制各地震属性不同井位处的曲线图,统计不同井位处地震属性的变
    化趋势,宏观上了解在不同沉积相带中地震属性参数的变化特征,将井位处地震属性特征
    与沉积相信息对比分析。

    步骤S222、统计全部已知井位处由步骤S21所得到的主变量与各地震属性间的相
    关系数,结合步骤S221地震属性与沉积相信息的吻合程度,选择相关系数较大的地震属性
    参数作为第一协同变量。

    步骤S23、将沉积相信息通过分类标号的方式转换为沉积相数字模型,使沉积相作
    为控制条件加入到建模计算过程中,实现真正的相控建模。具体过程如下:

    步骤S231、根据地质调查结果中的沉积相描述,使用数字表示沉积相变化,数字大
    小与远离陆地距离有关,具体的:1、2表示台内礁滩相;3表示台地潮坪相;4、5表示台缘斜坡
    相,利用定量的数字表示代替了常规的文字描述,如此得到的数字化沉积相信息就可以直
    接约束建模过程;使得平面上每一个点的沉积相数字表示可以将沉积相信息加入到建模的
    计算过程中。

    步骤S232、使用各向异性扩散法对所得到的数字化沉积相信息进行概化处理,得
    到沉积相数字模型作为第二协同变量;

    所述各向异性扩散法的计算公式如下:


    其中,Z为待处理点参数值,为各方向临近点梯度值,A、B、C为各方向的权重系
    数,i、j、k为所处理点的相对位置,m为点的迭代次数。通过求解公式(1)可以得到概化处理
    后的沉积相模型。

    其中,通过不断修改待处理点位置,对整个数字化沉积相信息数据体进行概化处
    理,最终得到所处理平面上的沉积相数字模型,如图4所示。

    步骤S3、将步骤S2所得主变量、第一协同变量及第二协同变量作为计算参数,利用
    同位多相协同克里金方法进行插值计算,获得平面上所有待估计点的参数值作为弹性参数
    的多源数据融合建模结果。

    其中,平面上所有待估计点的参数值的获取过程如下:

    步骤S31、以步骤S2中所述的主变量、第一协同变量和第二协同变量为计算参数,
    根据同位协同多相克里金方法进行插值计算,建立计算平面上待估计点的参数值模型:


    其中,待估计点位置的参数估计值为Z*(u0),测井主变量参数为Z(ui),u为参数点
    位置,提取的作为第一协同变量的地震属性参数为Y1,作为第二协同变量的沉积相数字模
    型为Y2,权重系数分别为αi、β1和β2。

    步骤S32、根据无偏最优条件,通过构建拉格朗日方程使估计克里金方差最小,可
    以得到克里金方程组:


    其中C表示参数间的协方差或互协方差,为了更清楚的表示参数间的关系,将方程
    组转换为矩阵形式:


    其中M表示协方差矩阵,Λ表示主变量权重系数向量,ε表示误差参数,M的表达式
    为:


    步骤S33、计算不同位置的待估计点测井参数、地震属性和沉积相之间的协方差和
    互协方差,可以得到矩阵M中所有变量值;

    步骤S34、通过求解公式(4)可以得到各变量的权重系数,将其代入公式(2)可以计
    算得到待估计点参数值。

    上述步骤S31至S34可计算出平面上的一个待估计点参数值,为了获取全部平面上
    所有待估计点的参数值,则修改待估计点位置,将步骤S34中计算得到的参数值视为已知
    点,重复步骤S3,计算下一个待估计点参数值,当平面上所有点全部计算完成后结束。

    步骤S4、使用各向异性扩散法对步骤3所得多源数据融合建模结果进行概化处理,
    得到概化处理后的地震弹性参数建模结果。具体包括:

    步骤S41、设定待估计点的参数变化阈值,将大于或小于设定参数变化阈值的待估
    计点的参数值视为异常点去除,即统计参数变化范围将负值及超出最大值数量级范围内的
    参数值视为异常点,随后将步骤S3多源数据融合建模结果中的异常点去除。

    步骤S42、利用各向异性扩散法概化处理对去除异常点的多元数据融合建模结果
    处理,即对剩余的待估计点参数值概化处理,通过公式(1)反复迭代可以得到概化处理后的
    地震弹性参数建模结果,如图5所示,所得结果保留了完整的参数变化边界和局部变化,并
    且在很大程度上消除了干扰建模清晰度的噪点,从图中可以看出采用本发明的方法可以获
    得更合理的弹性参数建模结果。

    因此,本发明的方法实现了利用沉积相数字模型作为协同变量的相控建模,提高
    了地球物理领域弹性参数的建模精度,使用各向异性扩散法去除计算噪点的同时保留了地
    质体所包含的局部信息,能够消除计算异常点和边界噪点,提高了多信息参数综合建模的
    精度并融入了更多真实的地质信息,具有较好适用性。

    上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施
    方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下
    做出各种变化。

    关于本文
    本文标题:基于同位多相协同克里金的地震弹性参数相控建模方法.pdf
    链接地址://www.4mum.com.cn/p-6001419.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精准计划 欢乐生肖规则 山东时时走势图 最新双色球规律与技巧 安徽快三计划软件手机版 极速赛车6码必中规律 北京pk赛车计划软件最准 水果机 什么是黑彩 德国pk10赛车计划 电子投注是否就是买了 新疆时时开奖结果表 四川时时下载 北京pk10定位胆计划 大乐透近100走势图