• 浪潮孙丕恕从信息时代到智能时代 人工智能价值将爆发式释放 2019-12-21
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
    • / 13
    • 下载费用:30 金币  

    重庆时时彩直播视频软件: 基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法.pdf

    关 键 词:
    基于 风力发电机 转子 电量 信息 碳刷 滑环 火花 监测 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201611143936.5

    申请日:

    2016.12.13

    公开号:

    CN106772036A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G01R 31/34申请日:20161213|||公开
    IPC分类号: G01R31/34 主分类号: G01R31/34
    申请人: 浙江运达风电股份有限公司
    发明人: 许国东; 邱纪星; 陶钧炳; 韩小良
    地址: 311106 浙江省杭州市余杭经济技术开发区(钱江经济开发区)顺风路558号
    优先权:
    专利代理机构: 杭州斯可睿专利事务所有限公司 33241 代理人: 王利强
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201611143936.5

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

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

    摘要

    一种基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法,包括如下步骤:1)机组并网运行时,实时采集发电机转子三相电流iA、iB、iC,利用Morlet小波,分别对其进行变换;2)提取小波脊线;3)根据小波脊线得出瞬时频率信息,检测瞬时频率中是否存在特征频率,即50Hz;4)对转子三相电流信号进行滤波、加窗、插值之后,进行FFT变换,分别比较当前在转子基频和50Hz处频谱幅值与正常工况下在转子基频和50Hz处频谱幅值的差异,如果两处幅值都有减小的趋势,判定为发生打火。本发明直接采集转子三相电流,数据获取方便,对双馈风电机组碳刷和滑环间打火的早期诊断具有一定的工程应用价值。

    权利要求书

    1.一种基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法,其特征在
    于:所述监测方法包括如下步骤:
    1)机组并网运行时,实时采集发电机转子三相电流iA、iB、iC。利用Morlet小波,分别对
    其进行变换;
    2)提取小波脊线;
    3)根据小波脊线得出瞬时频率信息,检测瞬时频率中是否存在特征频率,即50Hz;
    4)对转子三相电流信号进行滤波、加窗、插值之后,进行FFT变换,分别比较当前在转子
    基频和50Hz处频谱幅值与正常工况下在转子基频和50Hz处频谱幅值的差异,如果两处幅值
    都有减小的趋势,判定为发生打火。
    2.如权利要求1所述的基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方
    法,其特征在于:在步骤1)中,对转子三相电流进行处理的过程如下:
    1.1)设ψ(t)为平方可积函数,即ψ(t)∈L2(R),如果Fourier变换满足容许性条件:
    <mrow> <msub> <mi>C</mi> <mi>&psi;</mi> </msub> <mo>=</mo> <munder> <mo>&Integral;</mo> <mi>R</mi> </munder> <mfrac> <mover> <mrow> <msup> <mi>&psi;</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> </mrow> <mo>~</mo> </mover> <mrow> <mo>|</mo> <mi>&omega;</mi> <mo>|</mo> </mrow> </mfrac> <mi>d</mi> <mi>&omega;</mi> <mo>&lt;</mo> <mi>&infin;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
    则ψ(t)为母小波,将其按照尺度参数a和平移参数b进行伸缩和平移,生成的函数族
    {ψa,b(t)}为分析小波。
    <mrow> <msub> <mi>&psi;</mi> <mrow> <mi>a</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>|</mo> <mi>a</mi> <msup> <mo>|</mo> <mrow> <mo>-</mo> <mn>0.5</mn> </mrow> </msup> <mi>&psi;</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>t</mi> <mo>-</mo> <mi>b</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mo>,</mo> <mi>a</mi> <mo>,</mo> <mi>b</mi> <mo>&Element;</mo> <mi>R</mi> <mo>,</mo> <mi>a</mi> <mo>&NotEqual;</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
    1.2)变换中所用的Morlet小波时域形式如下。
    <mrow> <mi>&psi;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>&pi;</mi> <mrow> <mo>-</mo> <mn>1</mn> <mo>/</mo> <mn>4</mn> </mrow> </msup> <mrow> <mo>(</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>i&omega;</mi> <mn>0</mn> </msub> <mi>t</mi> </mrow> </msup> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msubsup> <mi>&omega;</mi> <mi>0</mi> <mi>2</mi> </msubsup> <mo>/</mo> <mn>2</mn> </mrow> </msup> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> <mo>/</mo> <mn>2</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
    1.3)根据步骤1.1)、1.2)中,可以得到信号各个频率的相位和幅值信息,对任意的(a,
    b),相位可表示为:

    对于式(4),任意的(a,b),其在t上只有一个驻点,满足:

    3.如权利要求2所述的基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方
    法,其特征在于:在步骤3)中,由驻点的性质知φ′a,b(ts)=0,而小波脊线是由相平面上满
    足ts(a,b)=b的点,两者结合即可得到脊线。


    因此知道小波脊线ar(b),即可求出瞬时频率其中小波脊线的提取可借助小波曲
    线。
    4.如权利要求1~3之一所述的基于双馈风力发电机转子侧电量信息的碳刷和滑环火
    花监测方法,其特征在于:在步骤3)中,小波脊线与瞬时频率的关系如下:

    式(2)中,ar(b)为小波脊线,为瞬时频率。
    5.如权利要求1~3之一所述的基于双馈风力发电机转子侧电量信息的碳刷和滑环火
    花监测方法,其特征在于:在步骤4)中,对转子三相电流进行如下处理:
    4.1)对转子相电流进行低通滤波、加窗函数,Hanning窗的时域形式如下:
    <mrow> <mi>w</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mfrac> <mrow> <mi>&pi;</mi> <mi>t</mi> </mrow> <mi>T</mi> </mfrac> <mo>)</mo> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mi>t</mi> <mo>|</mo> <mo>&le;</mo> <mi>T</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mi>t</mi> <mo>|</mo> <mo>&gt;</mo> <mi>T</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
    4.2)进行分段线性插值,各插值点用折线连接,再将插值之后的信号做FFT变换;
    4.3)将打火和正常运行时的处理结果作比较,观测两者在转子基频和50Hz处频谱的差
    异。

    说明书

    基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法

    技术领域

    本发明属于风电机组的在线故障诊断技术领域,涉及双馈风力发电机碳刷和滑环
    火花监测方法。

    背景技术

    对于MW级的双馈风电机组,励磁电流较大,通常同一相分布着多个碳刷,也使相互
    之间的均流效果难以得到保证。实际运行中,碳刷和滑环间电阻的变化,极易导致电刷电流
    的二次分配,使部分碳刷和滑环高温过热,再加上周围空气湿度大,两者之间发生打火的几
    率会显著增加,轻则损坏碳刷,重则造成发电机滑环及转轴的损坏,影响发电机的正常运
    行。

    实际中导致碳刷与滑环间出现电火花的原因很多,比如:整机的振动过大,引起电
    刷振动过大;滑环表面有灼伤或表面粗糙使两者接触不良;碳刷运行过热和周围环境湿度
    大;恒压弹簧压力不够导致碳刷在刷握内摇摆等。若能够对火花进行诊断,及时更换碳刷,
    将有非常重要的过程意义。

    目前双馈风力发电机组碳刷和滑环火花的监测可通过光电管监测法、电刷电流、
    接触电压的峰值等手段,但是这些方法大都需要增加设备,不仅成本大,而且可靠性不高、
    适应性较差。近来各研究机构及风机制造商也对两者间的火花监测进行了相关研究,以及
    公布了一些文章和专利,但大都从电刷电流的角度出发。如“双馈异步发电机电刷滑环电阻
    变化理论与仿真研究”,该方法通过对电刷与滑环的等效电阻进行分析,测量每个碳刷的电
    流并进行傅里叶变换,提取特征分量。该方法需要对电刷电流进行处理,在现有风机硬件基
    础上难以实现。

    发明内容

    为了克服现有的双馈风力发电机碳刷和滑环火花监测方式的成本大、可靠性较
    差、适应性较差的不足,本发明提供一种简单易行、准确性良好、适用性较好的基于双馈风
    力发电机转子侧电量信息的碳刷和滑环火花监测方法。

    本发明解决其技术问题所采用的技术方案是:

    一种基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法,所述监测
    方法包括如下步骤:

    1)机组并网运行时,实时采集发电机转子三相电流iA、iB、iC。利用Morlet小波,分
    别对其进行变换;

    2)提取小波脊线;

    3)根据小波脊线得出瞬时频率信息,检测瞬时频率中是否存在特征频率,即50Hz;

    4)对转子三相电流信号进行滤波、加窗、插值之后,进行FFT变换,分别比较当前在
    转子基频和50Hz处频谱幅值与正常工况下在转子基频和50Hz处频谱幅值的差异,如果两处
    幅值都有减小的趋势,判定为发生打火。

    进一步,所述步骤1)中,需实时采集发电机转子侧的电流,实际中可采用变流器采
    集的三相转子电流。

    再进一步,在步骤1)中,对转子三相电流进行处理的过程如下:

    1.1)设ψ(t)为平方可积函数,即ψ(t)∈L2(R),如果Fourier变换满足容许性
    条件:


    则ψ(t)为母小波,将其按照尺度参数a和平移参数b进行伸缩和平移,生成的函数
    族{ψa,b(t)}为分析小波。


    1.2)变换中所用的Morlet小波时域形式如下。


    1.3)根据步骤1.1)、1.2)中,可以得到信号各个频率的相位和幅值信息,对任意的
    (a,b),相位可表示为:


    对于式(4),任意的(a,b),其在t上只有一个驻点,满足:

    φ'a,b(ts)=0&φ”a,b(ts)≠0(5)。

    在步骤3)中,由驻点的性质知φ'a,b(ts)=0,而小波脊线是由相平面上满足ts(a,
    b)=b的点,两者结合即可得到脊线。



    因此知道小波脊线ar(b),即可求出瞬时频率其中小波脊线的提取可借助小
    波曲线。

    在步骤3)中,小波脊线与瞬时频率的关系如下:


    式(2)中,ar(b)为小波脊线,为瞬时频率。

    在步骤4)中,对转子三相电流进行如下处理:

    4.1)对转子相电流进行低通滤波、加窗函数,Hanning窗的时域形式如下:


    4.2)进行分段线性插值,各插值点用折线连接,再将插值之后的信号做FFT变换;

    4.3)将打火和正常运行时的处理结果作比较,观测两者在转子基频和50Hz处频谱
    的差异。

    本发明的有益效果主要表现在:1)本发明涉及的信息的提取无需增加新的设备,
    硬件成本较低;

    2)通过发明所述的方法可以有效的提取特征量,且与正常运行时的数据对比,可
    以明显观察到打火时的区别。

    附图说明

    图1为本发明流程图。

    图2为滑环和碳刷电阻理论示意图,其中,(a)是等效电路,(b)是碳刷与滑环滑动
    接触示意图。

    图3为实验机组滑环电阻周期性变化图。

    图4为打火时转子A相和B相电流。

    图5打火相瞬时频率信息。

    图6为打火时转子A相电流FFT分析。

    图7为正常运行时转子A相电流FFT分析。

    具体实施方式

    下面结合附图对本发明作进一步描述。

    参照图1~图7,一种基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测
    方法,包括如下步骤:

    1)机组并网运行时,实时采集发电机转子三相电流iA、iB、iC。碳刷和滑环之间发生
    打火时,会在转子电流中有所体现,现从两者间的电阻理论出发具体介绍,主要考虑接触电
    阻和滑环的周期性电阻。

    碳刷和滑环的接触面存在附件电阻,表达式如下:


    式中,ρl、ρr为碳刷、滑环的电阻率;ε为材料的变形系数;Hl、Hr为碳刷、滑环接触材
    料的布氏硬度;n、F分别是接触点数和接触压力;ρa膜层的电阻率。

    滑环的引线触头与碳刷接触部分时刻变化,故滑环电阻也做周期性变化,分析原
    理如图2所示。从图中可以看到电刷电流经2路流向滑环,一路是L1,另一路是L2,两段电阻
    可表示如下:


    式中,ρ为电阻率,A为滑环径向截面积。


    式中ω可由频率f来表示,即ω(t)=f×t×360。再有R1和R2的并联关系可得总的
    滑动电阻,其周期性变化如图3所示。

    Rtot(t)=Rring×f×(t-f×t2)Rring=2×π×r×ρ/A(13)

    进一步设转子A相电压u(t)=cos(ωt),对(4)式按泰勒变换转化成正余弦形式,
    可求出转子A相电流与碳刷和滑环电阻的关系,两者电阻的变化会引起转子电流在ω0、ω0+
    ωm频率处的变化。


    利用Morlet小波对其变换;

    2)提取小波脊线;

    3)根据小波脊线得出瞬时频率的变化,检测瞬时频率图中是否存在特征频率,即
    50Hz。

    4)为进一步区分,打火与正常运行时,两者在特征频率及基波频率处的区别,对转
    子三相电流信号进行滤波、加窗、插值之后,进行FFT变换,分别比较当前在转子基频和50Hz
    处频谱幅值与正常工况下在转子基频和50Hz处频谱幅值的差异,如果两处幅值都有减小的
    趋势,判定为发生打火。

    即当前转子基频和50Hz处频谱处的幅值,与正常工况时转子基频和50Hz处频谱处
    的幅值相比较,如果同时满足当前转子基频的幅值都比正常工况时基频的幅值小,当前转
    子50Hz的幅值都比正常工况时50Hz的幅值小,判定为发生打火。

    进一步,在步骤1)中,对转子三相电流进行处理的过程如下:

    1.1)设ψ(t)为平方可积函数,即ψ(t)∈L2(R),如果Fourier变换满足容许性
    条件:


    则ψ(t)为母小波,将其按照尺度参数a和平移参数b进行伸缩和平移,生成的函数
    族{ψa,b(t)}为分析小波。


    1.2)变换中所用的Morlet小波时域形式如下。


    1.3)根据步骤1.1)、1.2)中,可以得到信号各个频率的相位和幅值信息,对任意的
    (a,b),相位可表示为:


    对于式(9),任意的(a,b),其在t上只有一个驻点,满足:

    φ'a,b(ts)=0&φ”a,b(ts)≠0(5)

    1.4)由驻点的性质知φ'a,b(ts)=0,而小波脊线是由相平面上满足ts(a,b)=b的
    点。两者结合即可得到脊线。



    因此知道小波脊线ar(b),即可求出瞬时频率其中小波脊线的提取可借助小
    波曲线。

    在步骤3)中,小波脊线与瞬时频率的关系如下:


    式(2)中,ar(b)为小波脊线,为瞬时频率。

    再进一步,在步骤4)中,对转子三相电流进行如下处理:

    4.1)对转子相电流进行低通滤波、加窗函数,Hanning窗的时域形式如下:


    4.2)进行分段线性插值,各插值点用折线连接,再将插值之后的信号做FFT变换;

    4.3)将打火和正常运行时的处理结果作比较,观测两者在转子基频和50Hz处频谱
    的差异。

    本发明涉及的信息的提取无需增加新的设备,硬件成本较低;

    因故障特征量本身很小,提取困难。文中采用瞬时频率的方法进行提取。并且通过
    分析及实验结果可以明显看到打火时,转子电流基波幅值及特征频率处的幅值较正常运行
    时偏小。

    以某2.0MW双馈风力发电机组为研究对象,转子每相6个碳刷,为使打火易于发生
    使转子K相2个碳刷运行,实验时转速为1750rpm附近,输出有功1.5MW,转子相电流有效值
    410A。图2和图3分别是滑环和碳刷间电阻理论的等效示意图和实验机组的滑环周期性电
    阻,实验时在滑环室观察到的火花现象以及实验结束后,转子K相碳刷表面留下的电蚀现
    象。

    需要说明的是,实际实验时打火一般难以实现,故将转子K相2个碳刷运行,以增加
    电刷电流密度,增大打火的几率。

    记录打火前后转子电流数据,数据包时间长度为10s,图4记录的是该时间段内转
    子A和B相的电流,其中A相对应K相,B相对应L相,期间A相发生打火。采用Morlet小波分析,
    其瞬时频率如图5所示,图中平均频率对应转子基频8.3Hz附近,在6.3s左右出现50Hz分量。

    为进一步分析,A相打火与正常运行时的区别,分别记录两种工况的数据进行分
    析。图6和图7分别对应打火时和正常运行时A相电流的FFT分析结果,所用方法如上述所述。
    打火期间,对应转子基频8.3Hz频谱幅值523.3,特征频率50Hz处幅值8.014;而正常工况时,
    转子基频处幅值532.6,50Hz处幅值8.394。从结果看出不管是转子基频还是50Hz处,在打火
    发生时两者幅值明显减小,与根据滑环和碳刷电阻理论分析的结果相一致。

    根据实验结果,可以根据50Hz处幅值和转子基频幅值的变化判断碳刷和滑环轻微
    打火的发生,并且本发明的判定条件是基于大量实验数据的,并不局限于该实例。

    本发明通过双馈风电机组转子电流来判断碳刷和滑环之间的打火,简单有效,成
    本相对较低,是一种有效可靠的判断方法。

    最后需要说明的是,以上实施例仅用以说明本发明的技术方案而非限制,本领域
    技术人员应当理解,可以在形式上和细节上对本发明做出改变,而不偏离发明权利要求书
    所限定的范围。

    还需要说明的是,本文中所公开的实例描述的各参数,仅是为了更好的描述本发
    明,专业人员可以意识到,通过修改本发明的参数值可达到同样的诊断效果,但是这种实现
    不应认为超出本发明的范围。

    对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。
    本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实
    现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理
    和新颖特点相一致的最宽的范围。

       内容来自专利网重庆时时彩单双窍门 www.4mum.com.cn转载请标明出处

    关于本文
    本文标题:基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法.pdf
    链接地址://www.4mum.com.cn/p-6001410.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 浪潮孙丕恕从信息时代到智能时代 人工智能价值将爆发式释放 2019-12-21
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
  • 湖北11选5走势图表 青海快3开奖时间表 幸运农场高手指导 网球比分怎么算 手机捕鱼游戏编程原理 pc蛋蛋走势图 福彩开奖272 雷速体育直播篮球 南粤风彩36选7 河南十一选五中奖奖金 北京时时彩开奖号码 麻将技巧视频 6肖复式5肖多少组 吉林11选5彩票走势图 亚上彩群 辽宁快乐12组选走势图