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

    重庆时时彩黑彩软件: 基于储层地质力学参数的微地震事件正演模拟方法.pdf

    关 键 词:
    基于 地质力学 参数 地震 事件 模拟 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201410131339.5

    申请日:

    2014.04.02

    公开号:

    CN103913772A

    公开日:

    2014.07.09

    当前法律状态:

    终止

    有效性:

    无权

    法律详情: 未缴年费专利权终止IPC(主分类):G01V 1/28申请日:20140402授权公告日:20160831终止日期:20170402|||著录事项变更IPC(主分类):G01V 1/28变更事项:发明人变更前:潘树林 尹成 吴波 高磊 朱仕军 丁峰 赵虎 廖细明变更后:潘树林 秦子雨 吴波 尹成 高磊 朱仕军 丁峰 赵虎 廖细明|||授权|||实质审查的生效IPC(主分类):G01V 1/28申请日:20140402|||公开
    IPC分类号: G01V1/28 主分类号: G01V1/28
    申请人: 西南石油大学
    发明人: 潘树林; 尹成; 吴波; 高磊; 朱仕军; 丁峰; 赵虎; 廖细明
    地址: 610500 四川省成都市新都区新都大道8号
    优先权:
    专利代理机构: 济南日新专利代理事务所 37224 代理人: 崔晓艳
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410131339.5

    授权公告号:

    ||||||||||||

    法律状态公告日:

    2018.04.17|||2017.03.15|||2016.08.31|||2014.08.06|||2014.07.09

    法律状态类型:

    专利权的终止|||著录事项变更|||授权|||实质审查的生效|||公开

    摘要

    本发明提供一种基于储层地质力学参数的微地震事件正演模拟方法,该基于储层地质力学参数的微地震事件正演模拟方法包括:为初始水力压裂的压力赋值,使用有限元素法进行模型模拟;根据时间变化进行每一特定时刻的岩石单元应力计算;根据应力计算结果,结合岩石破裂准则,判断当前岩石单元是否发生破裂;计算孔隙度参数;若岩石单元发生破裂,结合微震震源机制计算震源震级;以及完成所有时刻的模拟。该基于储层地质力学参数的微地震事件正演模拟方法计算结果能够真实的反映出地下储层的实际变化情况,可以用该方法进行地下裂缝及压裂过程中微震事件的预测。

    权利要求书

    权利要求书
    1.  基于储层地质力学参数的微地震事件正演模拟方法,其特征在于,该基于储层地质力学参数的微地震事件正演模拟方法包括:
    步骤1,为初始水力压裂的压力赋值,使用有限元素法进行模型模拟;
    步骤2,根据时间变化进行每一特定时刻的岩石单元应力计算;
    步骤3,根据应力计算结果,结合岩石破裂准则,判断当前岩石单元是否发生破裂;
    步骤4,计算孔隙度参数;
    步骤5,若岩石单元发生破裂,结合微震震源机制计算震源震级;以及
    步骤6,重复步骤2、步骤3、步骤4和步骤5,完成所有时刻的模拟。

    2.  根据权利要求1所述的基于储层地质力学参数的微地震事件正演模拟方法,其特征在于,该基于储层地质力学参数的微地震事件正演模拟方法还包括,在步骤1之前,建立地下储层的渗流-应力-损伤耦合模型,为模型赋储层地质力学参数。

    3.  根据权利要求1所述的基于储层地质力学参数的微地震事件正演模拟方法,其特征在于,在建立地下储层的渗流-应力-损伤耦合模型的步骤中,该渗流-应力-损伤耦合模型采用的储层参数控制方程为:
    σ·ij=2Gϵ·ij+(K-2G3)ϵ·kkδij+αp·δij+γ1T·δij---(1)]]>
    ζ·=αϵ·ii+βp·-γ2T·---(2)]]>式中σijσij和εij分别表示总应力和总应变,p和T分别表示孔隙压力和温度,α表示Biot系数,ζ表示孔隙流体含量,δij表示Kronecker符号,K表示体积模量,G表示剪切模量,γ1,γ2,β的值通过下面公式给出:
    β=α-φKs+φKf---(3)]]>
    γ1=Kαm    (4)
    γ2=ααm+(αf-αm)φ    (5)
    φ表示孔隙度,αm和αf分别表示岩石和流体的热膨胀系数,岩石的体积模量用Ks表示,流体的体积模量使用Kf表示;假设流体在岩石孔隙中流动遵循达西定律,热传导遵循傅里叶定律;则有下式:
    Jf=-ρfkη▿p---(6)]]>
    JT=-kT▿T---(7)]]>ρf表示流体密度,k表示渗透率,η表示粘滞系数,kT表示热传导系数,Jf和JT分别表示流体流量和热量;下式(8)(9)为空隙中流体流动和应力平衡方程;
    σij,j=0    (8)
    ∂ζ∂t=-1ρf▿Jf---(9)]]>
    将前面的控制方程(1)(2)带入平衡方程(8)(9),得到关于岩石形变和流体流动状态的场方程(10)(11);
    (K+G3)▿(▿.u)+G▿2u+m(α▿p+γ1▿T)=0---(10)]]>
    α(▿.u·)+βp·-kη▿2p-γ2T·=0---(11)]]>
    T·+v(▿T)-cT▿2T=0---(12)]]>
    u为位移量,当模型为二维时m=[1,1,0]T,模型为三维时m=[1,1,1,0,0,0]T,cT为热扩散系数,流体的速度和孔隙压力遵循达西定律v=-kη▿p.]]>

    4.  根据权利要求1所述的基于储层地质力学参数的微地震事件正演模拟方法,其特征在于,在步骤3,在模型中,储层岩石受应力作用后的过程可以描述为弹性形变和破裂两个阶段;受较小应力时,岩石发生弹性形变,这一阶段岩石内部不发生破裂;而当岩石受力增大达到岩石的破坏准则时,岩石开始破裂,随着破裂的不断发展,最终形成裂缝。

    5.  根据权利要求4所述的基于储层地质力学参数的微地震事件正演模拟方法,其特征在于,在步骤3,在弹性形变-破裂机制中,弹性模量的减小代表了岩石发生破裂的过程,
    E=(1-d)E0    (13)
    d表示岩石内部发生破裂的程度,可以称为岩石破裂系数,E0表示岩石的初始岩性模量,E为当前岩性模量;如果岩石受压力发生了破裂,岩石破裂系数可以通过岩石的剩余压力来表述:
    d=1-fcrE0ϵ‾---(14)]]>
    fcr为岩石的残余压力,为对应压力下发生的应变,则表示了初始压应力。

    6.  根据权利要求1所述的基于储层地质力学参数的微地震事件正演模拟方法,其特征在于,在步骤4,获得某一时刻各个地下单元的实际受力情况,并且对比实际受力大小和岩石单元保持不破裂所承受的最大应力,进行是否破裂的判断,如果未发生破裂,则岩石单元发生形变,孔隙度参数根据形变发生改变;如果岩石破裂,则孔隙度发生突变,转化为一个较大的数值,该数值为常数。

    7.  根据权利要求1所述的基于储层地质力学参数的微地震事件正演模拟方法,其特征在于,在步骤5,如果发生破裂则根据应力的变化情况模拟出微震事件的强度,微震事件震源位置为当前计算单元位置,震动时刻为当前时间。

    8.  根据权利要求1所述的基于储层地质力学参数的微地震事件正演模拟方法,其特征在于,该基于储层地质力学参数的微地震事件正演模拟方法还包括,在步骤5之后,将所有模拟结果输出,作为分析和解释储层压力的依据。

    说明书

    说明书基于储层地质力学参数的微地震事件正演模拟方法
    技术领域
    本发明涉及石油工程中水力压裂模拟及地震勘探方法中的微地震模拟,特别是涉及到一种基于储层地质力学参数的微地震事件正演模拟方法。
    背景技术
    水力压裂是提高致密油气、页岩气等非常规油气资源产量的重要手段。在压裂的过程中,监测压裂效果,对压裂方案提出完善,对于提高产量,降低成本有着重要的意义。微地震监测技术是目前被认可的用来实时监测压裂效果的唯一有效手段。在向低渗透率储层注水或注气的过程中,会引起流体压力前缘的移动和孔隙流体压力的变化,从而引发微地震事件。通过微地震监测,可以实时地调整、优化设计方案,有效地提高油气采收率,实现油田的整体开发的科学化管理。通过水力压裂下的微地震监测,可以对压裂范围、裂缝发育方向和大小进行有效的追踪,定位,客观的评价压裂工程的效果,尤其通过一系列裂缝属性分析,定位主裂缝的构造走向,以及次生裂缝的分布,近似计算出裂缝的长度、宽度、高度;通过微地震监测,可以对岩石内部流体前缘进行成像。油藏工程师通过对裂缝成像和驱动前缘波及状况的分析,对下一步的生产开发提供有效的指导。对于一次开发而言,准确的裂缝方位和分布范围将有助于从井中最有效地开发出移动油气。对储层的二次开发而言,有助于确定采油井位,提高采油率,节约成本,增加效益。
    目前对于微地震产生的机理了解不够,使得微地震资料在实际应用中效果有待进一步完善。常规的微地震正演技术是假设震源和传播介质已知,研究其走时以及振幅和频率等传播规律的方法,没有结合储层地质力学参数。为此我们发明了一种新的基于储层地质力学参数的微地震事件正演模拟方法,解决了以上技术问题。
    发明内容
    本发明的目的是提供一种可以进行地下裂缝及压裂过程中微震事件的预测的基于储层地质力学参数的微地震事件正演模拟方法。
    本发明的目的可通过如下技术措施来实现:基于储层地质力学参数的微地震事件正演模拟方法,该基于储层地质力学参数的微地震事件正演模拟方法包括:步骤1,为初始水力压裂的压力赋值,使用有限元素法进行模型模拟;步骤2,根据时间变化进行每一特定时刻的岩石单元应力计算;步骤3,根据应力计算结果,结合岩石破裂准则,判断当前岩石单元是否发生破裂;步骤4,计算孔隙度参数;步骤5,若岩石单元发生破裂,结合微震震源机制计算震源震级;以及步骤6,重复步骤2、步骤3、步骤4和步骤5,完成所有时刻的模拟。
    本发明的目的还可通过如下技术措施来实现:
    该基于储层地质力学参数的微地震事件正演模拟方法还包括,在步骤1之前,建立地下储层的渗流-应力-损伤耦合模型,为模型赋储层地质力学参数。
    在建立地下储层的渗流-应力-损伤耦合模型的步骤中,该渗流-应力-损伤耦合模型采用的储层参数控制方程为:
    σ·ij=2Gϵ·ij+(K-2G3)ϵ·kkδij+αp·δij+γ1T·δij---(1)]]>
    ζ·=αϵ·ii+βp·-γ2T·---(2)]]>
    式中σijσij开口εij分别表示总应力和总应变,p和T分别表示孔隙压力和温度,α表示Biot系数,ζ表示孔隙流体含量,δij表示Kronecker符号,K表示体积模量,G表示剪切模量,γ1,γ2,β的值通过下面公式给出:
    β=α-φKs+φKf---(3)]]>
    γ1=Kαm    (4)
    γ2=ααm+(αf-αm)φ    (5)
    φ表示孔隙度,αm和αf分别表示岩石和流体的热膨胀系数,岩石的体积 模量用Ks表示,流体的体积模量使用Kf表示;假设流体在岩石孔隙中流动遵循达西定律,热传导遵循傅里叶定律;则有下式:
    Jf=-ρfkη▿p---(6)]]>
    JT=-kT▿T---(7)]]>ρf表示流体密度,k表示渗透率,η表示粘滞系数,kT表示热传导系数,Jf和JT分别表示流体流量和热量;下式(8)(9)为空隙中流体流动和应力平衡方程;
    σij,j=0    (8)
    ∂ζ∂t=-1ρf▿Jf---(9)]]>
    将前面的控制方程(1)(2)带入平衡方程(8)(9),得到关于岩石形变和流体流动状态的场方程(10)(11);
    (K+G3)▿(▿.u)+G▿2u+m(α▿p+γ1▿T)=0---(10)]]>
    α(▿.u·)+βp·-kη▿2p-γ2T·=0---(11)]]>
    T·+v(▿T)-cT▿2T=0---(12)]]>
    u为位移量,当模型为二维时m=[1,1,0]T,模型为三维时m=[1,1,1,0,0,0]T,cT为热扩散系数,流体的速度和孔隙压力遵循达西定律
    在步骤3,在模型中,储层岩石受应力作用后的过程可以描述为弹性形变和破裂两个阶段;受较小应力时,岩石发生弹性形变,这一阶段岩石内部不发生破裂;而当岩石受力增大达到岩石的破坏准则时,岩石开始破裂,随着破裂的不断发展,最终形成裂缝。
    在步骤3,在弹性形变-破裂机制中,弹性模量的减小代表了岩石发生破 裂的过程,
    E=(1-d)E0    (13)
    d表示岩石内部发生破裂的程度,可以称为岩石破裂系数,E0表示岩石的初始岩性模量,E为当前岩性模量;如果岩石受压力发生了破裂,岩石破裂系数可以通过岩石的剩余压力来表述:
    d=1-fcrE0ϵ‾---(14)]]>
    fcr为岩石的残余压力,为对应压力下发生的应变,则表示了初始压应力。
    在步骤4,获得某一时刻各个地下单元的实际受力情况,并且对比实际受力大小和岩石单元保持不破裂所承受的最大应力,进行是否破裂的判断,如果未发生破裂,则岩石单元发生形变,孔隙度参数根据形变发生改变;如果岩石破裂,则孔隙度发生突变,转化为一个较大的数值,该数值为常数。
    在步骤5,如果发生破裂则根据应力的变化情况模拟出微震事件的强度,微震事件震源位置为当前计算单元位置,震动时刻为当前时间。
    该基于储层地质力学参数的微地震事件正演模拟方法还包括,在步骤5之后,将所有模拟结果输出,作为分析和解释储层压力的依据。
    本发明中的基于储层地质力学参数的微地震事件正演模拟方法,与常规的微地震模拟算法相比,考虑了压裂过程中储层地质力学参数的变化,因此计算结果能够真实的反映出地下储层的实际变化情况,可以用该方法进行地下裂缝及压裂过程中微震事件的预测。正演是模拟水力压裂实施过程中,储层受应力变化后的情况。使用该技术可以模拟出实际施工作业时,地下储层裂缝的产生、发育以及压裂过程中微震事件的震源位置、震源强度等信息。这些信息对于指导压裂作业有着重要的意义。
    附图说明
    图1为本发明的基于储层地质力学参数的微地震事件正演模拟方法的一具体实施例的流程图;
    图2为本发明的一具体实施例中模型弹性模量分布图;
    图3为本发明的一具体实施例中模型原始地应力分布图;
    图4为本发明的一具体实施例中正演模拟获得的模型微震事件示意图。
    具体实施方式
    为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合所附图式,作详细说明如下。
    如图1所示,图1为本发明的基于储层地质力学参数的微地震事件正演模拟方法的一具体实施例的流程图。
    在步骤101,建立地下储层的渗流-应力-损伤耦合模型,为模型赋储层地质力学参数。
    对水力压裂模拟必须结合流体压力,岩石受压力形变和岩石破裂过程的描述。地下裂缝的生成本质为岩石破裂,其产生和发展主要受岩石受到的张力影响。储层岩石的渗透率和应力是一个耦合的函数,当岩石破裂后,渗透率急剧增加;
    本发明模型采用的储层参数控制方程如下式:
    σ·ij=2Gϵ·ij+(K-2G3)ϵ·kkδij+αp·δij+γ1T·δij---(1)]]>
    ζ·=αϵ·ii+βp·-γ2T·---(2)]]>式中σijσij和εij分别表示总应力和总应变,p和T分别表示孔隙压力和温度,α表示Biot系数,ζ表示孔隙流体含量,δij表示Kronecker符号,K表示体积模量,G表示剪切模量,γ1,γ2,β的值通过下面公式给出:
    β=α-φKs+φKf---(3)]]>
    γ1=Kαm    (4)
    γ2=ααm+(αf-αm)φ    (5)
    φ表示孔隙度,αm和αf分别表示岩石和流体的热膨胀系数,岩石的体积 模量用Ks表示,流体的体积模量使用Kf表示。假设流体在岩石孔隙中流动遵循达西定律,热传导遵循傅里叶定律。则有下式:
    Jf=-ρfkη▿p---(6)]]>
    JT=-kT▿T---(7)]]>
    ρf表示流体密度,k表示渗透率,η表示粘滞系数,kT表示热传导系数,Jf和JT分别表示流体流量和热量。下式(8)(9)为空隙中流体流动和应力平衡方程。
    σij,j=0    (8)
    ∂ζ∂t=-1ρf▿Jf---(9)]]>
    将前面的控制方程(1)(2)带入平衡方程(8)(9),得到关于岩石形变和流体流动状态的场方程(10)(11)。
    (K+G3)▿(▿.u)+G▿2u+m(α▿p+γ1▿T)=0---(10)]]>
    α(▿.u·)+βp·-kη▿2p-γ2T·=0---(11)]]>
    T·+v(▿T)-cT▿2T=0---(12)]]>
    u为位移量,当模型为二维时m=[1,1,0]T,模型为三维时m=[1,1,1,0,0,0]T,cT为热扩散系数,流体的速度和孔隙压力遵循达西定律
    流程进入到步骤102。
    在步骤102,为初始水力压裂的压力赋值,对储层地质力学参数赋初值。使用方程(10)-(12)对个网格单元进行描述,使用有限元素法进行模型模拟。流程进入到步骤103。
    在步骤103,根据时间变化进行每一特定时刻的岩石单元应力计算。流 程进入到步骤104。
    在步骤104,根据应力计算结果,结合岩石破裂准则,判断当前岩石单元是否发生破裂。在模型中,储层岩石受应力作用后的过程可以描述为弹性形变和破裂两个阶段。受较小应力时,岩石发生弹性形变,这一阶段岩石内部不发生破裂。而当岩石受力增大达到岩石的破坏准则时,岩石开始破裂,随着破裂的不断发展,最终形成裂缝。在弹性形变-破裂机制中,弹性模量的减小代表了岩石发生破裂的过程。
    E=(1-d)E0    (13)
    d表示岩石内部发生破裂的程度,可以称为岩石破裂系数,E0表示岩石的初始岩性模量,E为当前岩性模量。如果岩石受压力发生了破裂,岩石破裂系数可以通过岩石的剩余压力来表述:
    d=1-fcrE0ϵ‾---(14)]]>
    fcr为岩石的残余压力,为对应压力下发生的应变,则表示了初始压应力。
    在进行储层地质力学参数描述时,流体的初始压力,及储层的初始地质力学参数如泊松比、弹性模量、初始地应力等都作为已知条件。需要求解的是流体压力改变后,储层内部应力及应变的变化情况。使用有限元素法求解方程(10)(11),并在求解的过程中加入摩尔-库伦岩石破裂判断准则,就可以获得在一定压裂条件下储层的响应。流程进入到步骤105。
    在步骤105,进行计算,获得某一时刻各个地下单元的实际受力情况,并且对比实际受力大小和岩石单元保持不破裂所承受的最大应力,进行是否破裂的判断,如果未发生破裂,则岩石单元发生形变,孔隙度参数根据形变发生改变。如果岩石破裂,则孔隙度发生突变,转化为一个较大的数值(该数值为常数)。流程进入到步骤106。
    在步骤106,若岩石单元发生破裂,结合震源机制来最终获得微地震事件信息。如果发生破裂则根据应力的变化情况模拟出微震事件的强度,微震事件震源位置为当前计算单元位置,震动时刻为当前时间。流程进入到步骤 107。
    在步骤107,重复步骤103到步骤106,完成所有时刻的模拟。流程进入到步骤108。
    在步骤108,将所有模拟结果输出,作为分析和解释储层压力的依据。流程结束。
    为了验证算法,建立了一个500米宽,500米长的一个二维模型,划分最小网格单元为边长5米的正方形。地质力学参数使用了某工区已知参数。图2为模型的弹性模量分布图,图3为地下应力的分布图。对该模型使用发明算法进行模拟?;竦昧送?所示的微震事件分布结果图。
    本发明提供一种油气田水力压裂过程中,基于储层地质力学参数的微地震事件模拟方法,主要应用于对水力压裂时产生微地震事件特征的分析和预测进一步压裂地下储层的变化及将要生成的微地震事件信息等。该方法主要流程是通过建立渗流-应力-损伤耦合的储层地质力学参数模型来描述储层的状态,通过对模型表述方程在边界条件下进行求解,从而实现对水力压裂过程地下储层变化的数值模拟,模拟将获得地下储层参数的变化,包括地质力学参数改变、地下裂缝分布,并且结合微地震震源机制原理,获得压裂过程中微地震事件产生的位置、时间和震级等信息,正演获得的微地震事件信息,对于处理和解释水力压裂过程中为了监控地下裂缝发育产生而采集的实际微地震资料具有重要意义。为了能够使用数学方程描述储层的受力及破裂情况,方法做以下假设:储层同时具有弹性和脆性,并且其状态使用弹性损伤力学进行描述;储层岩石受到的正应力和切应力遵循摩尔—库仑破裂准则;储层的渗透率是应力造成的弹性形变的函数,在岩石发生破裂或断层生成时,渗透率将急剧增加。在这些假设条件下,建立起符合研究要求的数学描述方程。在给定的边界条件下,就可以根据方程求解储层破裂情况,了解和掌握地下裂缝的分布。

    关于本文
    本文标题:基于储层地质力学参数的微地震事件正演模拟方法.pdf
    链接地址://www.4mum.com.cn/p-6130432.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
  • 天津11选5购买 山西11选5规则 河南十一选五选胆技巧 天津十一选五开将公告 山西11选5推荐导航 河南快赢481官网电话 大乐透后区生肖选号 股票配资名片 pc蛋蛋预测软件100%中 大乐透开奖结果 全球股票指数行情,东方财富网 QQ分分彩的漏洞 购买彩票合买 福建福彩七乐彩走势图 青海十一选五今天开奖 河北11选5技巧稳赚