计算机测量与控制．２０２２．３０（１） 犆狅犿狆
狌狋犲狉犕犲犪狊狌狉犲犿犲狀狋牔犆狅狀狋狉狅犾
·７８ ·
收稿日期：２０２１１１２３； 修回日期：
２０２１１２１６。基金项目：山西省自然科学基金项目（２０１９０３Ｄ１２１１３３；２０１９０１Ｄ２１１２５１
）。作者简介：白文斌（１９９７）
，男，山西晋中人，硕士，主要从事信息探测与处理方向的研究。通讯作者：李 凯（１９７９）
，男，山西太原人，硕士生导师，副教授，主要从事信号与信息处理、信息探测、电磁检测、传感器方向的研究。 引用格式：白文斌，李 凯．三线圈电感式油液金属磨粒检测系统［Ｊ］．计算机测量与控制，２０ ２２，３０（１）：７８８５．
文章编号：１６７１４５９８（２０２２）０１００７８０８ ＤＯＩ：１０．１６５２６／ｊ
．ｃｎｋｉ．１１－４７６２／ｔｐ．２０２２．０１．０１３ 中图分类号：ＴＰ３文献标识码：Ａ
三线圈电感式油液金属磨粒检测系统
白文斌，李 凯
（中北大学信息探测与处理山西省重点实验室，太原 ０３００５１
）摘要：机械传动金属磨粒可以反馈机械设备故障的特性，针对机械传动系统故障在线检测需求，提出了一种三线圈电感式金属磨粒检测系统；通过建立传感器数学模型，对影响传感器金属探测灵敏度的参数进行仿真，结合传感器实际使用情况求出最优解；通过设计励磁信号源模块、调幅模块、信号调理模块等电路并结合上位机对油液中的金属磨粒进行检测；经孔径为８ｍｍ的流道的实验测试，实现了５００μｍ铁磨粒及１０００μｍ铜磨粒的检测精度；该系统为发动机中金属磨粒的检测提供了技术支持，对故障预防与诊断具有重要意义。
关键词：油液；电感式；磨粒传感器；电磁检测；解调
犜犺狉犲犲犆狅犻犾犐狀犱狌犮狋犻狏犲犗犻犾犕犲狋犪犾犃犫狉犪狊犻狏犲犘犪狉狋犻犮犾犲犇犲狋犲犮狋犻狅狀犛狔
狊狋犲犿ＢＡＩＷｅｎｂｉｎ，ＬＩＫａｉ
（ＳｈａｎｘｉＰｒｏｖｉｎｃｉａｌＫｅｙＬａｂｏｒａｔｏｒｙｏｆＩｎｆｏｒｍａｔｉｏｎＤｅｔｅｃｔｉｏｎａｎｄＰｒｏｃｅｓｓｉｎｇ
，ＮｏｒｔｈＵｎｉｖｅｒｓｉｔｙｏｆＣｈｉｎａ，Ｔａｉｙ
ｕａｎ ０３００５１，Ｃｈｉｎａ）犃犫狊狋狉犪犮狋：Ｔｈｅｍｅｔａｌｗｅａｒｄｅｂｒｉｓｉｎｍｅｃｈａｎｉｃａｌｔｒａｎｓｍｉｓｓｉｏｎｃａｎｆｅｅｄｂａｃｋｔｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｍｅｃｈａｎｉｃａｌｅ ｑｕｉｐｍｅｎｔｆａｉｌｕｒｅ，ａｉｍｉｎｇａｔｔｈｅｒｅｑｕｉｒｅｍｅｎｔｏｆｏｎｌｉｎｅｆａｕｌｔｄｅｔｅｃｔｉｏｎｏｆｍｅｃｈａｎｉｃａｌｔｒａｎｓｍｉｓｓｉｏｎｓｙｓｔｅｍ，ａｔｈｒｅｅ－ｃｏｉｌｉｎｄｕｃｔｉｖｅｍｅｔａｌａｂｒａｓｉｖｅｐａｒｔｉｃｌｅｄｅｔｅｃｔｉｏｎｓｙｓｔｅｍｉｓｐｒｏｐｏｓｅｄ．Ｂｙｅｓｔａｂｌｉｓｈｉｎｇｔｈｅｍａｔｈｅｍａｔｉｃａｌｍｏｄｅｌｏｆｔｈｅｓｅｎｓｏｒ，ｔｈｅｐａｒａｍｅｔｅｒｓａｆｆｅｃｔｉｎｇｔｈｅｍｅｔａｌｄｅｔｅｃｔｉｏｎｓｅｎｓｉｔｉｖｉｔｙｏｆｔｈｅｓｅｎｓｏｒａｒｅｓｉｍｕｌａｔｅｄ，ａｎｄｔｈｅｏｐ
ｔｉｍａｌｓｏｌｕｔｉｏｎｉｓｏｂｔａｉｎｅｄｃｏｍｂｉｎｅｄｗｉｔｈｔｈｅａｃｔｕａｌｕｓｅｏｆｔｈｅｓｅｎｓｏｒ．Ｔｈｒｏｕｇｈｔｈｅｄｅｓｉｇｎｏｆｅｘｃｉｔａｔｉｏｎｓｉｇｎａｌｓｏｕｒｃｅｍｏｄｕｌｅ，ａｍ ｐｌｉｔｕｄｅｍｏｄｕｌａｔｉｏｎｍｏｄｕｌｅ，ｓｉｇｎａｌｃｏｎｄｉｔｉｏｎｉｎｇｍｏｄｕｌｅａｎｄｏｔｈｅｒｃｉｒｃｕｉｔｓ，ｃｏｍｂｉｎｅｄｗｉｔｈｔｈｅｈｏｓｔｃｏｍｐｕｔｅｒ，ｔｈｅｍｅｔａｌａｂｒａｓｉｖｅｐａｒｔｉｃｌｅｓｉｎｔｈｅｏｉｌａｒｅｄｅｔｅｃｔｅｄ．Ｔｈｒｏｕｇｈｔｈｅｅｘｐｅｒｉｍｅｎｔａｌｔｅｓｔｏｆｔｈｅｒｕｎｎｅｒｗｉｔｈａｎａｐｅｒｔｕｒｅｏｆ８ｍｍ，ｔｈｉｓｓｙｓｔｅｍｃａｎａｃｈｉｅｖｅｔｈｅｄｅｔｅｃｔｉｏｎａｃｃｕｒａｃｙｏｆ５００μｍｉｒｏｎａｂｒａｓｉｖｅｇｒａｉｎｓａｎｄ１０００μｍｃｏｐｐｅｒａｂｒａｓｉｖｅｇｒａｉｎｓ．Ｔｈｅｓｙｓｔｅｍｐｒｏｖｉｄｅｓｔｅｃｈｎｉｃａｌｓｕｐｐｏｒｔｆｏｒｔｈｅｄｅｔｅｃｔｉｏｎｏｆｍｅｔａｌａｂｒａｓｉｖｅｐａｒｔｉｃｌｅｓｉｎｅｎｇｉｎｅ，ａｎｄｉｓｏｆｇｒｅａｔｓｉｇｎｉｆｉｃａｎｃｅｆｏｒｆａｕｌｔｐｒｅｖｅｎｔｉｏｎａｎｄｄｉａｇ
ｎｏｓｉｓ．犓犲狔
狑狅狉犱狊：ｏｉｌ；ｉｎｄｕｃｔｉｖｅ；ａｂｒａｓｉｖｅｓｅｎｓｏｒ；ｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｄｅｔｅｃｔｉｏｎ；ｄｅｍｏｄｕｌａｔｉｏｎ０ 引言
现代工业、民用、军事为了准确控制机械设备的运行，液压系统和润滑系统广泛应用于机械自动化等领域，对机械设备在运转过程的在线故障诊断需求逐
渐提高［１］
。机械设备在运行过程中由于相互摩擦、碰撞、污染物侵入等原因产生磨粒［
２］
，其在线、离线检测是判断机械设备磨损的重要手段［３］
，也是预测设备
使用寿命、故障检测的重要依据，成为了一个重要的
机械设备故障监测研究方向［４］。
投稿网址：ｗｗｗ．ｊｓｊｃｌｙ
雨水弃流井
红外线加热灯ｋｚ．ｃｏｍ
第１期
白文斌，等：
三线圈电感式油液金属磨粒检测系统
·７９ ·
油液在线检测传感器可分为光学［５］、超声［６
］、电容［７］
和电感传感器。基于金属磨粒通过油路时产生的
磁场扰动具有较好的金属类别表征特性，电感检测法在油液中金属磨粒检测领域受到高度重视。国外有关金属磨粒电感探测方法研究较早，主要包括加拿大ＧａｓＴＯＰＳ公司研发的ＭｅｔａｌＳＣＡＮ传感器与英国
Ｋｉｔｔｉｗａｋｅ公司研制的ＦＧ型电感传感器［８］
。
针对微流道油液传感器，大连海事大学做了很多研究。为提高微流道油液传感器检测精度，文献［９
］通过在传感器内部填充超顺磁性材料提升磁场强度，对于颗粒检测信噪比提升了２０％左右，与未填充磁性纳米材料的传感器相比，具有更高的检测信噪比以及更低的检测下限。为提高微流道油液传感器吞吐量，文献［１０］基于相分复用原理搭建的双通道检测系统，在提高检测流量的同时，达到了检测７４～８８μ
ｍ铁磨粒的效果。基于对非铁磁性颗粒检测更为敏感的电阻检测法，文献［１１］提出了一种双螺线管线圈结构的电感电阻微流控传感器，有效弥补电感传感器对非铁磁性颗粒检测能力低的缺陷，完成了５００μ
ｍ流道中２０μｍ铁磨粒和６０μｍ铜磨粒的检测。基于此检测方法，文献［１２］提出一种平面线圈与铁氧体磁芯组合的新型结构，增大吞吐量的同时提高了信噪比。为实现微流道油液传感器多污染物的检测，基于电感与电容检测法，张洪朋团队共提出了４类多参数测量传感器，可检测润滑油中存在金属磨粒、气泡与水滴。其一，文献［１３］将两个单层平面线圈正对排布在微通道两侧，通过调节激励源实现电感电容模式的切换。在３００μｍ流道中，电容模式可对１８０μｍ水滴和２４０μｍ气泡进行区分检测，电感模式可对８０μ
ｍ铁颗粒和１５０μｍ铜颗粒进行区分检测。其二，为提升检测分辨率，文献［１４］采用流道垂直嵌入两个贴合的平面线圈内孔的传感器，在３００μｍ流道中成功检测到４０μｍ铁颗粒、１１０μｍ
、１００μｍ水滴及１８０μｍ气泡。为提升高吞吐量流道的检测精度，文献［１５］将油液通道改为环状结构增大流量同时在油液通道内的毛细玻璃管加入硅钢片，在缩短检测时间的同时可完成４０μｍ铁磁性颗粒、１３０μｍ铜颗粒的检测，切换为电容检测模式时，可完成２００μｍ水滴和２７０μｍ气泡的检测。在此基础上，文献［１６］将平面线圈优化为螺线管线圈结构，保持高通量的同时提升了检测精度，可检测到３０μｍ铁颗粒、１１０μｍ
铜颗粒、１００μｍ水滴和１８０μｍ气泡。其三，文献
［１７
］提出一种包含电容式传感器和电感式传感器的集成式磨粒检测装置，该装置能够检测并区分８０μｍ气泡，３０μｍ铁颗粒和４５μｍ铜颗粒。目前关于微流道金属磨粒检测技术的研究较多且可达到的检测精度较高，但是微流道金属磨粒检测精度是在牺牲传感器流量的情况下实现的，为了保证检测精度，常将微流道传感器的孔径设置过小而受到限制，无法满足大流量机械设备油路中的金属磨粒检测，而且由于孔径过小，使得油液在流道口处受到较大的阻力，管路承受的压力增大，进而产生杂质堆积，造成流道堵塞，同时，由于传感器通量较小，通常只能安装在润滑油旁路，导致传感器仅对部分润滑油进行取样检测，检测结果随机性较大。
目前电感式磨粒传感器的主要技术瓶颈是传感器的检测精度受限于流道孔径。高通量传感器可完成多
种机械设备润滑油路的高吞吐量金属磨粒检测，但存
在检测精度低等问题［
１８］
口吃矫正器。现阶段关于电磁感应传感器技术主要包括结构设计技术、电路补偿技术。结构
设计主要包括传感器最优结构参数设计［
１９］
、一激多感及多激一感等新型传感器探测结构设计、传感器内部磁场增强结构设计等，电路补偿技术主要包括谐振电路、相分复用电路、激励线圈不平衡补偿电
路［
２０
］等。针对目前高通量传感器存在检测精度低的问题，本文设计了基于ＦＰＧＡ的电感式金属磨粒检测系统，通过对双激励传感器结构建立数学模型，对影响传感器灵敏度参数进行仿真，结合实际情况maop
并求出最优解。经实验验证：本系统可在油液流量为５～１０Ｌ／ｍｉｎ的工况下实现５００μｍ铁磨粒及１０００μｍ铜磨粒的检测精度。
１ 电感式金属磨粒检测原理
电感式磨粒检测传感器结构如图１所示，线圈骨架采用不影响传感器磁场分布的酚醛塑酯材料设计而成。激励线圈位于骨架两侧且对称放置，线圈绕向既可同向也可反向，骨架中部为感应线圈。两激励线圈同向绕制时，由相位相差１８０°的信号源驱动，而激励线圈反向绕制时，仅由单个信号源驱动即可。当信号源驱动激励线圈后，传感器内部产生激励磁场，两参数相同的激励线圈绕制时采用反向绕制，在同一信
投稿网址：ｗｗｗ．ｊｓｊｃｌｙ
ｋｚ．ｃｏｍ
计算机测量与控制 第３０
卷
·８０ ·号源的驱动下，位于传感器骨架中点位置的磁场强度接近于零。当有金属磨粒通过传感器时，基于磨粒自身的磁化或涡流效应，扰动中心零磁场状态，感应线圈内部产生感应电动势。通过分析传感器输出电压信号的变化趋势及幅值，即可判断金属磨粒的数量、类别、大小等信息
。
图１ 电感式磨粒检测传感器结构
金属磨粒通过传感器时感应电动势的数学模型见式（１
）。犈＝－３．７０８π狉４μ０犖（狀－犿）犐犿［
狉２＋（狀－犿）２］５２μ１狉３
１狏２狉２＋１．２３６μ１狉２１（１） 令：
犓＝－３．７０８π狉４
μ０犖（狀－犿）犐犿［狉２＋（
狀－犿）２］５
２
（２
隔离dcdc电源） 式（
１）可简化为：犈ｍａｘ＝犓μ１狉３
１狏２狉２＋１．２３６μ
１狉２
１（３）式中，狉１为球形磨粒半径，μ１为磨粒相对磁导率，狏为油液流速，狋为磨粒经过第一个线圈时间，狀为激励线圈到感应线圈中心距离，狉为螺线管内径，μ０为真空磁导率，犿为初级线圈长度，犖为激励线圈匝数，犐为电流大小，犓值中包含影响传感器灵敏度的传感器结构参数。
２ 检测系统框架设计
金属磨粒检测系统整体设计方案如图２所示，主要由磨粒传感器模块（
激励线圈Ａ、激励线圈Ｂ、感应线圈）、信号源模块、调幅模块、电源模块、信号调理模块、通信模块构成。磨粒传感器模块用于获取油液中的磨粒信息，电源模块为其他模块供电，信号源模块输出的正弦交流电压通过调幅模块实现对磨粒传感器的激励，激励线圈产生极性相反的动态交变磁场保证感应线圈处于零磁场状态，当磨粒通过传感器时，信号调理模块对感应线圈输出的电压信号进行滤波、放大、解调等信号调理，通信模块将调理后的电压信号传输至计算机进行数字信号处理。同时在传感
器模块与整体电路模块加入磁场防干扰结构减弱外部磁场干扰，增强系统稳定性
。
图２ 检测系统框架
２ １ 励磁信号源
传感器输出的应信号特性受信号源影响，励磁信号源中的高次谐波会在感应线圈上产生同频率的感应信号，干扰信号解调，为此设计的励磁信号源应具有频率、相位、幅值稳定可调、波形失真小等特点。
本方案采用直接数字合成器技术（ＤＤＳ）
［２１］
，ＤＤＳ以数控振荡器的方式，产生频率、相位可调的信号源，主要由基准时钟模块、相位累加器、相位调制器、正弦ＲＯＭ查表、ＤＡＣ、低通滤波电路构成。
输出信号源的频率及相位分别由频率控制字和相位控制字计算控制。ＤＤＳ系统的核心是相位累加器，时钟脉冲每触发一次，累加器将频率控制数据与相位寄存器输出的累加相位数据相加后的结果输入至相位寄存器的数据输入端。相位寄存器将累加后的数据重新输入至累加器，以便加法器在下一时钟的作用下继续与频率控制数据相加。当相位累加器累加满时，产生溢出，完成一个频率周期的动作，相位累加器的溢出频率就是ＤＤＳ输出的信号频率。相位寄存器的输出与相位控制字相加，结果作为正弦查表的地址，根据正弦ＲＯＭ表中的地址，将相应的正弦矢量值发送到８位ＤＡＣ的输入端，通过ＤＡＣ可将ＦＰＧＡ输出的８位数字信号的波形幅值转换成所要求的合成离散的模拟信号。
除此之外，为满足感应信号解调条件，还需输出一路励磁信号作为解调模块的参考频率信号。基于ＤＤＳ模块，可获得频率、相位、幅值稳定可调的正弦稳流激励信号，根据感应信号数学模型、励磁线圈匝数条件并结合实际传感器输出信号特点设置励磁信
投稿网址：ｗｗｗ．ｊｓｊｃｌｙ
ｋｚ．ｃｏｍ
霍尔式角度传感器第１期
白文斌，等： 三线圈电感式油液金属磨粒检测系统
·８１ ·
号源激励频率。
２ ２ 低通滤波电路
ＤＡＣ输出端为离散的模拟信号，需要将离散的模拟信号变为平滑的模拟信号，同时输出信号中含有高次谐波成分，为避免高次谐波影响后续信号解调，设计无源ＬＣ低通滤波电路，用于衰减和滤除不需要的取样分量，以便输出频谱纯净的正弦波信号。低通滤波器的截止频率为４０ＭＨｚ，频率响应曲线如图３所示
。
图３ 频率响应曲线图
２ ３ 调幅电路
ＤＡＣ芯片输出为差分输出，为此基于高性能１４５ＭＨｚ带宽的运放ＡＤ８０６５设计调幅电路将差分输出变为单端输出，但经由ＡＤ８０６５芯片设计的单端输出电路后输出的电压幅值无法满足传感器所需的激励幅值的要求，由于绕制工艺的限制，相同参数的传感器所需最佳激励幅值也会有差异，基于检测系统普适性要求，同时考虑到传感器输出的感应电动势大小受励磁信号源输出幅度的影
响，需根据感应信号数学模型、励磁线圈匝数等条件并结合实际传感器输出信号特点设置励磁信号幅值大小，因此设计调幅电路如图４所示，调幅电路最终输出范围为－５～５Ｖ。
２ ４ 解调电路
当信号源激励线圈，金属磨粒通过传感器时，感
应线圈输出的信号为振幅调制信号［２２］
，需使用解调
方法进行频谱搬移，将磨粒信号从载波中解调出来。解调方法包括包络检波、相干解调等，传统的二极管包络检波后的信号会出现较大失真，因此采用相干解调方法
［２３］
。
乘法器电路如图５所示。相干解调需满足调制信号与激励信号同频同相，在进行相干解调时，
调制信
图４ 调幅电路
号与ＤＤＳ输出的一路励磁信号作为解调模块的参考频率信号相乘后，经过低通滤波器取出低频分量，此刻输出信号为单个磨粒通过传感器时的感应信号
。
图５ 乘法器电路
３ 仿真
３ １ 模型参数的确定
由感应电动势的数学模型可知，传感器输出的感应电动势即传感器的探测灵敏度主要受传感器结构参数与磨粒参数的影响，传感器结构参数主要包括传感器孔径、激励线圈长度、激励线圈匝数、激励线圈到感应线圈中心距离等，磨粒参数主要包括磨粒半径、磨粒相对磁导率、磨粒流速等。
为了使磨粒通过传感器时输出的感应电动势能够有效表达金属磨粒的特性，即传感器的金属磨粒探测灵敏度达到最优，使用ＭＡＴＬＡＢ对影响传感器灵敏度因素的变量进行仿真，求得影响传感器灵敏度因素的最优解，传感器参数表如表１所示。
表１ 传感器参数
参数数值参数数值狉１
／ｍ２５０×１０－６犖
／匝数３００μ
０／（Ｈ／ｍ）４π×１０－７
狉／ｍ０．００５μ
１／（Ｈ／ｍ）６００μ０
狀／ｍ０．０３５狏／（ｍ·ｓ－１
）５犿／ｍ０．０２犳
／ｋＨｚ１００
犐／Ａ
１
投稿网址：ｗｗｗ．ｊｓｊｃｌｙ
ｋｚ．ｃｏｍ
计算机测量与控制 第３０
卷
·８２ ·３ ２ 磨粒半径与感应电动势关系
当磨粒材料相同粒径不同时，磨粒扰动中心零磁场状态产生的电感变化量不同，不同的电感变化量影响传感器输出的感应电动势，传感器可通过感应电动势的大小区分材料相同但粒径不同的磨粒，为明确磨粒粒径大小与传感器输出感应电动势的变化关系，根据感应电动势的数学模型，对磨粒半径与感应电动势关系进行仿真，仿真采用的磨粒材料为铁磁性材料，磁导率为μ１，磨粒半径的变化范围是０～５００μｍ，磨粒流速为５ｍ／ｓ
。磨粒半径与感应电动势的关系仿真结果如图６所示，从仿真结果可以看出，当磨粒半径为零时，传感器输出的感应电动势为零，仿真处于理想情况，在忽略其他因素干扰的条件下，当磨粒粒径为零（即没有磨粒通过传感器）时，感应电动势输出即为零，将磨粒粒径为零代入感应电动势的数学模型，感应电动势输出结果也为零，在实际传感器运行过程中，当没有磨粒通过传感时，传感器的输出信号应包含环境噪声及电路噪声，但此类噪声对磨粒通过传感器时输出感应电动势干扰较小，仿真结果同样与实际情况相符。
当粒径为微米级别的磨粒通过传感器时，传感器输出的感应电动势量级为纳伏级别，由于微米级别的磨粒通过传感器时，感应线圈上引起的电感变化范围小，导致输出的感应电动势也十分微弱，因此在实际应用中，需要对传感器输出的感应电动势进行放大、滤波等信号调理。从图中可以看出，随着磨粒半径逐渐增大，感应电动势也逐渐增大，二者之间近似呈三次方的变化关系，与数学模型中磨粒半径与感应电动势的变化关系保持一致。不同粒径的磨粒产生的感应电动势相差较大，仿真证明，传感器可以对微米级别粒径的磨粒进行检测，同时根据感应电动势的幅值大小对磨粒粒径进行区分。
３ ３ 磨粒磁导率与输出电压关系
当磨粒粒径相同材料不同时，磨粒扰动中心零磁场状态产生的电感变化量不同，进而影响传感器输出的感应电动势，不同的感应电动势的幅值可代表粒径相同但磁导率不同的磨粒，根据感应电动势的数学模型，对磨粒磁导率与感应电动势关系进行仿真，仿真采用的磨粒粒径为５００μｍ，磨粒流速为５ｍ／ｓ，磨粒磁导率的变化范围是０～１０００Ｈ／ｍ
。
图６ 磨粒半径与感应电动势关系图
磨粒磁导率与感应电动势关系仿真结果如图７所示，从仿真图可以看出，随着磨粒磁导率的增加，感
应电动势逐渐增大，同一尺寸不同磁导率的金属磨粒输出的感应电动势不同，当磨粒通过传感器时，传感器输出的感应电动势量级为微伏级，当磨粒磁导率从０～３００Ｈ／ｍ变化时，不同磁导率的磨粒产生的感应电动势相差较大，当磁导率超过９００Ｈ／ｍ时，感应电动势变化趋于平缓
。
图７ 磨粒磁导率与输出电压关系图３ ４ 传感器参数与传感器灵敏度关系
传感器检测灵敏度受机械设备实际运行工况的影响，在机械设备的实际运转中，其运动部件在轴承、气缸以及导板等设备内部不同位置产生相对运动，发生磨损，由上述磨损原因产生的金属磨粒具有较大区别，具体表现为金属磨粒的成分与尺寸，即磨粒磁导率与磨粒粒径不同。当油路孔径一定时，磨粒流速是由润滑油液流量决定的，当润滑油裹挟金属磨粒通过
投稿网址：ｗｗｗ．ｊｓｊｃｌｙ
ｋｚ．ｃｏｍ
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议