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

    重庆时时彩后一毒胆: 一种考虑尾水管水压脉动情况的水电孤岛频率抑制方法.pdf

    关 键 词:
    一种 考虑 水管 水压 脉动 情况 水电 孤岛 频率 抑制 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201410019214.3

    申请日:

    2014.01.16

    公开号:

    CN103779870A

    公开日:

    2014.05.07

    当前法律状态:

    终止

    有效性:

    无权

    法律详情: 未缴年费专利权终止IPC(主分类):H02J 3/36申请日:20140116授权公告日:20160203终止日期:20170116|||授权|||实质审查的生效IPC(主分类):H02J 3/36申请日:20140116|||公开
    IPC分类号: H02J3/36; H02J1/00; G06F17/50 主分类号: H02J3/36
    申请人: 华北电力大学; 中国电力科学研究院; 国网上海市电力公司
    发明人: 刘崇茹; 张梦琳; 张健; 韩玉蓉; 王建明; 魏佛送; 王歆; 李庚银
    地址: 102206 北京市昌平区回龙观朱辛庄2号
    优先权:
    专利代理机构: 北京麟保德和知识产权代理事务所(普通合伙) 11428 代理人: 周恺丰
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410019214.3

    授权公告号:

    |||||||||

    法律状态公告日:

    2018.03.06|||2016.02.03|||2014.06.11|||2014.05.07

    法律状态类型:

    专利权的终止|||授权|||实质审查的生效|||公开

    摘要

    本发明公开了电力系统运行和控制技术领域中的一种考虑尾水管水压脉动情况的水电孤岛频率抑制方法。包括根据水力系统和电力系统的耦合数学模型,计算水电孤岛直流送端系统的固有振荡频率;对水轮机模型进行改进,得到考虑尾水管水压脉动的详细水轮机模型;令详细水轮机模型谐振函数中的特征频率和尾水管压力脉动频率都等于系统的固有振荡频率,采集由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据;对由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据进行抑制,得到抑制后的水电孤岛直流送端系统振荡频率数据。本发明设计的详细水轮机模型和频率限制器能抑制考虑实际尾水管水压脉动引起的孤岛送端频率振荡,适用范围更广。

    权利要求书

    权利要求书
    1.  一种考虑尾水管水压脉动情况的水电孤岛频率抑制方法,其特征是所述方法包括:
    步骤1:根据水力系统和电力系统的耦合数学模型,计算水电孤岛直流送端系统的固有振荡频率fF;
    所述水力系统和电力系统的耦合数学模型包括水力系统模型和电力系统模型;
    所述水力系统模型包括水轮机模型和调速器模型;
    所述电力系统模型包括发电机模型和励磁系统模型;
    步骤2:对水轮机模型进行改进,得到考虑尾水管水压脉动的详细水轮机模型;
    所述考虑尾水管水压脉动的详细水轮机模型为:
    ΔPm(s)=KmH1-Zptanh(Teps)1+0.5tanh(Teps)ΔG(s)+1.5Hw1+0.5Tws;]]>
    其中,KmH为容量转换系数;
    ZP为引水管的规格化水力阻抗且ZP=Tw/Tep;
    Tw=Lv0/gH0;
    L为引水管管道长度;
    v0为初始水流速度;
    H0为初始水头;
    g为重力加速度常数;
    Tep为引水管的弹性时间;
    tanh(·)为双曲正切函数;
    ΔG(s)为导叶开度偏差;
    G(s)为导叶开度的函数;
    s为拉普拉斯算子;

    Y(t)为谐振函数且Y(t)=L-1[G×(sωc)1+2ζ(sωc)+(sωc)2];]]>
    ωc=2πfc且fc为特征频率;
    G为增益系数;
    ζ为阻尼系数;
    A为尾水管压力脉动的幅值;
    f为尾水管压力脉动频率;
    为正弦波的初始相位;
    *为卷积运算;
    步骤3:令谐振函数中的特征频率fc和尾水管压力脉动频率f都等于系统的固有振荡频率fF,采集由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据f0;
    步骤4:对由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据f0 进行抑制,得到抑制后的水电孤岛直流送端系统振荡频率数据f2和水电孤岛直流送端系统的实际传输功率。

    2.  根据权利要求1所述的方法,其特征是所述步骤4具体包括:
    子步骤A1:根据公式消除由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据中的直流分量;
    其中,f0为由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据;
    f1为消除直流分量后的水电孤岛直流送端系统振荡频率数据;
    G1为抑制环节增益系数;
    T1为微分环节时间常数;
    s为拉普拉斯算子;
    子步骤A2:根据公式消除f1中的所有高次谐波;
    其中,f2为消除所有高次谐波和直流分量后的水电孤岛直流送端系统振荡频率数据;
    z为阻尼系数;
    ω0为截止频率;
    子步骤A3:通过PI控制器将频率调整量转换为功率调整量;
    其中,频率调整量的计算公式为
    f2为消除所有高次谐波和直流分量后的水电孤岛直流送端系统振荡频率数据;
    Kp为比例系数;
    T2为积分时间常数;
    s为拉普拉斯算子;
    子步骤A4:将功率调整量与水电孤岛直流送端系统的额定直流功率相加,得到水电孤岛直流送端系统的实际传输功率。

    说明书

    说明书一种考虑尾水管水压脉动情况的水电孤岛频率抑制方法
    技术领域
    本发明属于电力系统运行和控制技术领域,尤其涉及一种考虑尾水管水压脉动情况的水电孤岛频率抑制方法。
    背景技术
    伴随着西电东送战略的推进,一大批水电基地正在或即将在我国西南地区兴建,该类新建水电站装机容量大,受电端一般位于东部沿海城市,水电站提供的电能传输距离较远。而与传统的交流输电相比,高压直流输电具有输送灵活、损耗小、能够实现快速控制等优点,所以很多水电站都采用高压直流输电方式提高外送电力。例如,向家坝-上海、云南-广东以及正在建设中的溪洛渡-浙西±800Kv直流工程都是水电机组采用高压直流输电的典型代表。但值得注意的一点是,该类水电站大都地处偏僻山区,常通过超高压线路与送端主网单点联系,在一定条件下可能形成孤岛,有些甚至直接采用孤岛运行方式,不与送端主网相联,这就形成了水电直流送端孤岛系统。
    水电直流送端孤岛系统自身网架结构单薄,维持频率稳定性能力弱,再加上水轮机尾水管压力脉动与直流系统之间的耦合作用,如果控制措施不当,极有可能引发水轮机组参与下的低频振荡现象。因此,研究高压直流输电送端水电机组孤岛运行方式下的频率抑制方法显得尤为重要。
    目前,在研究水力系统和电力系统耦合作用的方面,由于缺乏实际的尾水管压力脉动数据,建模与仿真中均用正弦或余弦激励信号来模拟尾水管水压脉动,这大大简化了实际的脉动情况,不能真实地展现脉动过程。在频率抑制方 法上,也是以采用直流频率限制器为主,但这种直流频率限制器仅能抑制用正弦或余弦激励信号模拟的尾水管水压脉动引发的孤岛频率振荡,不具有普遍适用性。
    故本发明设计了一种谐振函数来模拟尾水管的实际脉动情况,并建立了既考虑尾水管水压脉动又考虑弹性水柱的水轮机模型,由此得到由尾水管水压脉动引起的HVDC(高压直流)送端水电孤岛频率振荡波形,通过新设计的直流频率限制器来抑制由尾水管水压脉动引起的孤岛送端频率振荡。
    发明内容
    本发明的目的在于,提供一种考虑尾水管水压脉动情况的水电孤岛频率抑制方法,用于解决传统尾水管水压脉动建模过于简单和孤岛频率抑制方法过于单一的问题。
    为了实现上述目的,本发明提出的技术方案是,一种考虑尾水管水压脉动情况的水电孤岛频率抑制方法,其特征是所述方法包括:
    步骤1:根据水力系统和电力系统的耦合数学模型,计算水电孤岛直流送端系统的固有振荡频率fF;
    所述水力系统和电力系统的耦合数学模型包括水力系统模型和电力系统模型;
    所述水力系统模型包括水轮机模型和调速器模型;
    所述电力系统模型包括发电机模型和励磁系统模型;
    步骤2:对水轮机模型进行改进,得到考虑尾水管水压脉动的详细水轮机模型;
    所述考虑尾水管水压脉动的详细水轮机模型为:
    ΔPm(s)=KmH1-Zptanh(Teps)1+0.5tanh(Teps)ΔG(s)+1.5Hw1+0.5Tws;]]>
    其中,KmH为容量转换系数;
    ZP为引水管的规格化水力阻抗且ZP=Tw/Tep;
    Tw=Lv0/gH0;
    L为引水管管道长度;
    v0为初始水流速度;
    H0为初始水头;
    g为重力加速度常数;
    Tep为引水管的弹性时间;
    tanh(·)为双曲正切函数;
    ΔG(s)为导叶开度偏差;
    G(s)为导叶开度的函数;
    s为拉普拉斯算子;

    Y(t)为谐振函数且Y(t)=L-1[G×(sωc)1+2ζ(sωc)+(sωc)2];]]>
    ωc=2πfc且fc为特征频率;
    G为增益系数;
    ζ为阻尼系数;
    A为尾水管压力脉动的幅值;
    f为尾水管压力脉动频率;
    为正弦波的初始相位;
    *为卷积运算;
    步骤3:令谐振函数中的特征频率fc和尾水管压力脉动频率f都等于系统的固有振荡频率fF,采集由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据f0;
    步骤4:对由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据f0进行抑制,得到抑制后的水电孤岛直流送端系统振荡频率数据f2和水电孤岛直流送端系统的实际传输功率。
    所述步骤4具体包括:
    子步骤A1:根据公式消除由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据中的直流分量;
    其中,f0为由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据;
    f1为消除直流分量后的水电孤岛直流送端系统振荡频率数据;
    G1为抑制环节增益系数;
    T1为微分环节时间常数;
    s为拉普拉斯算子;
    子步骤A2:根据公式消除f1中的所有高次谐波;
    其中,f2为消除所有高次谐波和直流分量后的水电孤岛直流送端系统振荡 频率数据;
    z为阻尼系数;
    ω0为截止频率;
    子步骤A3:通过PI控制器将频率调整量转换为功率调整量;
    其中,频率调整量的计算公式为
    f2为消除所有高次谐波和直流分量后的水电孤岛直流送端系统振荡频率数据;
    Kp为比例系数;
    T2为积分时间常数;
    s为拉普拉斯算子;
    子步骤A4:将功率调整量与水电孤岛直流送端系统的额定直流功率相加,得到水电孤岛直流送端系统的实际传输功率。
    本发明设计了一个全新的数学模型来模拟尾水管的脉动情况,并且建立了既考虑尾水管水压脉动又考虑弹性水柱的详细水轮机模型,这对于通过时域仿真抑制共振型低频振荡提供了更真实的数学模型,为找到更好的抑制尾水管水压脉动引起的孤岛送端频率波动的方法也提供了理论基??;同时在设计直流频率限制器上也做出了改进,由于模型的改变,新设计的直流频率限制器能抑制考虑实际尾水管水压脉动引起的孤岛送端频率振荡,适用范围更广。
    附图说明
    图1是水电孤岛高压直流送端系统模型示意图;
    图2是水力系统和电力系统的耦合数学模型示意图;
    图3是包括阀控和伺服系统的调速器PID数学模型示意图;
    图4是利用matlab中的prony工具包拟合的频率曲线图;
    图5是加入谐振函数后模拟尾水管水压脉动情况引起的送端频率波动曲线图;
    图6是利用正弦波模拟尾水管水压脉动情况引起的送端频率波动曲线图;
    图7是根据本发明提供的抑制方法设计频率限制器模型示意图;
    图8是投入频率限制器后的送端频率曲线。
    具体实施方式
    下面结合附图,对优选实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
    图1是水电孤岛高压直流送端系统模型示意图。本实施例以图1所示在电力系统模型在仿真软件pscad上的仿真系统为例,说明本发明在HVDC送端水电机组孤岛运行方式下的对频率的抑制作用。
    仿真系统中直流主电路模型为±500kV,单极12脉动,直流传输功率1000MV,控制系统基本模型采用CIGRE直流输电标准测试系统的控制器;由于交流受端系统强度远远大于孤岛送端系统,因而受端交流系统??槟诓课刃У奈尴薮蟮缭茨P?。
    仿真模型的整流侧是包括4台450MW的水轮机、调速器、发电机和励磁系统的水力系统和电力系统的耦合模型,如图2所示。其中,G′为导叶开度,ω为角速度,ωref是基准角速度,Pm为水轮机输出机械功率,Vf为励磁输出电压,Vt为机端电压。
    基于该仿真系统,本发明考虑尾水管水压脉动情况的水电孤岛频率抑制方法的实施过程是:
    步骤1:根据水力系统和电力系统的耦合数学模型,计算水电孤岛直流送 端系统的固有振荡频率fF。
    发电机选用三阶实用模型,假定:(1)忽略定子绕组电磁暂态过程(2)阻尼绕组的作用计入阻尼系数D,数学模型为:
    Δδ=ω0ΔωTjΔω=Δpm-Δpe-DΔωTd0ΔEq=ΔEfd-ΔEq---(1)]]>
    其中Δδ和Δω分别为转子角和角速度增量,ΔE'q、ΔEq和ΔEfd分别为q轴暂态电势、空载电势和励磁电势增量,Δpm和Δpe分别为机械功率和电磁功率增量,Tj和T'd0分别为机组惯性时间常数和d轴开路暂态时间常数,ω′0为发电机的额定角速度。
    励磁系统模型的传递函数为:
    ΔEfd=Ge(s)(ΔV′f-ΔVt)              (2)
    Ge(s)=KA(1+TFs)(1+TAs)(1+Tes)(1+TFs)+KAKFs---(3)]]>
    其中:Vt为机端电压,ΔVt为机端控制电压扰动值,V′f为励磁附加控制信号,ΔV′f为励磁附加控制信号扰动值。KA和TA分别为电压调节器的放大倍数和时间常数。Te为励磁机时间常数,KF和TF分别为转子电压软反馈的放大倍数和时间常数,s为拉普拉斯算子。
    调速器模型选用包括阀控和伺服系统的PID数学模型,传递函数如图3所示。其中,K′p为比例系数,KI为积分系数,K0为微分系数,TB是伺服系统时间常数,Rp是永态转差率,TC是导叶伺服增益,TD是导叶伺服电机时间常数。μ为导叶开度,ω为角速度,ωref是基准角速度。
    水轮机模型先选用计及弹性水柱的非线性水轮机模型,在仿真平台上建立水力系统和电力系统的耦合模型,仿真得到孤岛送端频率的波形,通过matlab中的prony工具包采样波形中的点来计算系统的固有振荡频率fF。计及弹性水柱的非线性水轮机模型是
    ΔPm(s)=KmH1-ZPtanh(Teps)1+0.5ZPtanh(Teps)ΔG(s)---(4)]]>
    其中,KmH为容量转换系数,ZP为引水管的规格化水力阻抗且ZP=Tw/Tep,Tw=Lv0/gH0,L为引水管管道长度,v0为初始水流速度,H0为初始水头,g为重力加速度常数,Tep为引水管的弹性时间,tanh(·)为双曲正切函数,ΔG(s)为导叶开度偏差,G(s)为导叶开度的函数,s为拉普拉斯算子。取不同的展开级数n时,可得到不同的近似传递函数。对于本发明,n=1即满足研究要求。
    将仿真时长设为70s,将水轮机参数Tw=1秒、Te=1秒。仿真得到水电孤岛直流送端系统的固有振荡频率fF的波形,用matlab中的prony工具包对3个周期选取400点进行拟合,拟合图形如图4所示。经过prony计算可得到水电孤岛直流送端系统的固有振荡频率fF=0.11Hz。
    步骤2:对水轮机模型进行改进,得到考虑尾水管水压脉动的详细水轮机模型。
    将尾水管压力脉动的模型加入到计及弹性水柱的水轮机模型中形成考虑尾水管水压脉动并计及弹性水柱的非线性水轮机模型,其公式如下:
    ΔPm(s)=KmH1-Zptanh(Teps)1+0.5tanh(Teps)ΔG(s)+1.5Hw1+0.5Tws---(5)]]>
    其中,KmH为容量转换系数,ZP为引水管的规格化水力阻抗且ZP=Tw/Tep,Tw=Lv0/gH0,L为引水管管道长度,v0为初始水流速度,H0为初始水头,g为重 力加速度常数,Tep为引水管的弹性时间,tanh(·)为双曲正切函数,ΔG(s)为导叶开度偏差,G(s)为导叶开度的函数,s为拉普拉斯算子。
    Hw为增加了谐振函数的尾水管水压脉动的扰动量,其公式为:

    Y(t)为谐振函数,且
    Y(t)=L-1[G×(sωc)1+2ζ(sωc)+(sωc)2]---(7)]]>
    在时域中解得
    Y(t)=Be-ζωct---(8)]]>
    B=c[cosh(ζ2-1·ωc·t)-ζ·sinh(ζ2-1·ωc·t)ζ2-1]---(9)]]>
    其中,ωc=2πfc且fc为特征频率,G为增益系数,ζ为阻尼系数,A为尾水管压力脉动的幅值,f为尾水管压力脉动频率,为正弦波的初始相位,*为卷积运算。
    步骤3:令谐振函数中的特征频率fc和尾水管压力脉动频率f都等于系统的固有振荡频率fF,通过仿真得到由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据f0。
    其中,仿真中令A=0.02p.u,fc=fF=0.11Hz(fF为搭建的仿真系统的固有振荡频率),Tw=1秒、Te=1秒,ωc=2πfc,G=1,ζ=-0.001。仿真时长100s,扰动持续50s后消失。因为根据电力系统强迫功率振荡理论,持续的周期性小扰 动会引起电力系统强迫振荡,扰动频率与系统固有振荡频率相同时会引起系统谐振,导致功率大幅振荡,对于送端孤岛高压直流输电系统来说,会导致频率大幅振荡。所以fc、f和fF三者频率都设为0.11Hz。仿真时长设为100s,仿真所得到的模拟尾水管水压脉动引起的送端频率振荡如图5所示。将送端孤岛的振荡频率与仅仅使用的正弦波扰动得到的波形(如图6所示)进行对比,对比图5和图6可知,加入谐振函数更准确的模拟了由尾水管水压脉动所引起的送端孤岛频率振荡的情况,对抑制措施的提出具有更好的指导作用。而只采用正弦波模拟尾水管水压脉动没有考虑到尾水管水压脉动的自激过程,与实际情况差别过大。
    步骤4:对由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据f0进行抑制,得到抑制后的水电孤岛直流送端系统振荡频率数据f2和水电孤岛直流送端系统的实际传输功率。
    根据图5的送端频率出现的低频振荡现象,加入抑制措施,如图7所示。先根据公式消除由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据中的直流分量。f0为由尾水管水压脉动引起的水电孤岛直流送端系统振荡频率数据,f1为消除直流分量后的水电孤岛直流送端系统振荡频率数据,G1为抑制环节增益系数,T1为微分环节时间常数,s为拉普拉斯算子。
    再根据公式消除f1中的所有高次谐波。f2为消除所有高次谐波和直流分量后的水电孤岛直流送端系统振荡频率数据,也即抑制后的水电孤岛直流送端系统振荡频率数据,z为阻尼系数,ω0为截止频率。
    再通过PI控制器将频率调整量转换为功率调整量。如图7所示,将抑制后 的频率通过PI环节得到频率调整量。其中,频率调整量的计算公式为f2为抑制后的水电孤岛直流送端系统振荡频率数据,Kp为比例系数,T2为积分时间常数,s为拉普拉斯算子。
    在本实施例中,取G1=100,T1=0.01,z=5,ω0=0.15,Kp=30,T2=15。加上频率限制器后,PI环中的限幅取±500MV,将送端频率有效的抑制在49.8Hz-50.2Hz之间,如图8所示,抑制效果很好。由于通过PI环节后,频率调整量被转换为功率调整量Δp,因此再将功率调整量Δp与水电孤岛直流送端系统的额定直流功率相加,得到水电孤岛直流送端系统的实际传输功率。
    从上述实施例可以看出,改进后的频率限制器可以有效的抑制考虑实际尾水管水压脉动引起的HVDC送端孤岛频率振荡的问题,而且控制器操作简单,适用范围更广。
    以上所述,仅为本发明较佳的具体实施方式,但本发明的?;し段Р⒉痪窒抻诖?,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的?;し段е?。因此,本发明的?;し段вΩ靡匀ɡ蟮谋;し段?。

    关于本文
    本文标题:一种考虑尾水管水压脉动情况的水电孤岛频率抑制方法.pdf
    链接地址://www.4mum.com.cn/p-6185469.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 中国福彩时时彩计划 双色球有没有历史记录 天津快乐十分走势图 易发棋牌正宗官方下载 杰克棋牌完整版 3d今日预测号码预测分析 湖北11选5胆拖计算器 腾讯棋牌欢乐麻将 后三组三 重庆百变王牌中奖技巧 2006排列三走势图 北京pk10前二8码万能码 腾游棋牌游戏 2011年上证指数走势图 大智慧手机炒股