CN113686251A - 一种综采工作面设备上窜下滑偏移测量方法及系统 - Google Patents

一种综采工作面设备上窜下滑偏移测量方法及系统 Download PDF

Info

Publication number
CN113686251A
CN113686251A CN202110955008.3A CN202110955008A CN113686251A CN 113686251 A CN113686251 A CN 113686251A CN 202110955008 A CN202110955008 A CN 202110955008A CN 113686251 A CN113686251 A CN 113686251A
Authority
CN
China
Prior art keywords
point cloud
cloud data
distance
offset
reference plane
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110955008.3A
Other languages
English (en)
Other versions
CN113686251B (zh
Inventor
臧万顺
冯龙
张强
苏金鹏
田莹
沈伟挺
刘文卓
赵进
姜玉燕
王博申
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN202110955008.3A priority Critical patent/CN113686251B/zh
Publication of CN113686251A publication Critical patent/CN113686251A/zh
Priority to PCT/CN2021/135365 priority patent/WO2023019808A1/zh
Application granted granted Critical
Publication of CN113686251B publication Critical patent/CN113686251B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/02Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
    • G01B11/03Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness by measuring coordinates of points
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/02Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
    • G01B11/026Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness by measuring distance between sensor and object

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本申请实施例公开了一种综采工作面设备上窜下滑偏移测量方法,方法包括:接收第一激光雷达采集的第一巷道壁对应的第一点云数据;以及,接收第二激光雷达采集的第二巷道壁对应的第二点云数据;分别对第一点云数据与第二点云数据进行滤波处理以及曲面拟合处理,得到第一基准平面与第二基准平面;确定第一激光雷达与第一基准平面之间的第一距离,以及第二激光雷达与第二基准平面之间的第二距离;基于第一距离与预存的第一标准距离,确定第一偏移量,以及基于第二距离与预存的第二标准距离,确定第二偏移量;根据第一偏移量以及第二偏移量,确定综采设备发生的偏移类型。用以解决现有的综采设备偏移量测量方法的测量精度低的技术问题。

Description

一种综采工作面设备上窜下滑偏移测量方法及系统
技术领域
本申请涉及偏移测量领域,尤其涉及一种综采工作面设备上窜下滑偏移测量方法及系统。
背景技术
煤矿综采工作面的煤层都是天然形成的,煤层厚度分布不均,且具有一定的坡度。根据坡度的不同,可将煤层分为近水平层、缓倾斜层、倾斜层、急倾斜层。在不同工作面的煤炭采集过程中,刮板输送机等综采设备在重力、推力等因素的作用下,会在采煤工作面倾斜的方向上,产生“向上”或“向下”的分力,造成刮板运输机机头处搭接错位,成为引发事故发生的导火索。危险程度与上窜下滑的偏移量成正比,且在刮板输送机连续推进时,由于液压支架与刮板输送机连接销耳间隙窜动量的累积,以及传感器误差,导致刮板运输机上窜下滑偏移加剧。如果不能及时加以检查以及调整,将导致刮板输送机的端头驱动部损坏、液压支架挤架或倒架、推杆损坏、电缆损坏、生产中断等问题继而导致停产。由于在综采设备推进过程中,发生上窜或者下滑偏移是随时的,偏移量是实时发生变化的,因此,需要实现对于设备上窜下滑偏移的实时测试,以便进行实时调整。
目前,对于综采设备偏移量测量方法主要通过激光测距、超声波测距传感器等设备直接测量综采设备与巷道的实时距离,一方面,由于这些测距方法多数是将采集的信息传输到计算机进行处理,测量效率低。另一方面,综采工作面在有设备工作时会出现粉尘,一般的光原理测距仪器会把粉尘扫描进去,有一定的概率将扫描到的粉尘当作巷道的侧壁,导致测量精度大大降低,且煤矿巷道的侧壁是凹凸不平的,如果只测量巷道侧壁上一个点的距离并不能代表综采设备与巷道的真实距离,也会导致偏移量的测量出现偏差。
发明内容
本申请实施例提供了一种综采工作面设备上窜下滑偏移测量方法及系统,用于解决如下技术问题:现有的综采设备偏移量测量方法的测量精度低。
本申请实施例采用下述技术方案:
本申请实施例提供了一种综采工作面设备上窜下滑偏移测量方法,方法包括:接收第一激光雷达采集的第一巷道壁对应的第一点云数据;以及,接收第二激光雷达采集的第二巷道壁对应的第二点云数据;分别对所述第一点云数据与所述第二点云数据进行滤波处理以及曲面拟合处理,以得到第一基准平面与第二基准平面;确定所述第一激光雷达与所述第一基准平面之间的第一距离,以及所述第二激光雷达与所述第二基准平面之间的第二距离;基于所述第一距离与预存的第一标准距离,确定第一偏移量,以及基于所述第二距离与预存的第二标准距离,确定第二偏移量;根据所述第一偏移量以及所述第二偏移量,确定所述综采设备发生的偏移类型。
在一种可行的实施方式中,分别对所述第一点云数据与所述第二点云数据进行滤波处理以及曲面拟合处理,以得到第一基准平面与第二基准平面,具体包括:对所述第一点云数据进行高斯滤波,去除所述第一点云数据中的孤立点;以及,对所述第二点云数据进行高斯滤波,去除所述第二点云数据中的孤立点;通过预设算法,对进行高斯滤波后的所述第一点云数据进行曲面拟合,得到所述第一基准平面;以及,对进行高斯滤波后的所述第二点云数据进行曲面拟合,得到所述第二基准平面。
在一种可行的实施方式中,对所述第一点云数据进行高斯滤波,去除所述第一点云数据中的孤立点;以及,对所述第二点云数据进行高斯滤波,去除所述第二点云数据中的孤立点,具体包括:根据
Figure BDA0003219899730000031
得到二维高斯滤波权函数h(x,y);其中,
Figure BDA0003219899730000032
x代表偏离二维高斯滤波权函数中心的点云数据的x坐标,y代表偏离二维高斯滤波权函数中心的点云数据的y坐标,λxc以及λyc代表低通高斯滤波器的截止波长;根据w(x,y)=∫∫z(x-ξ,y-η)h(ξ,η)dξdη,得到二维高斯滤波函数w(x,y);其中,ξ、η为卷积积分所需的微分变量,z(x-ξ,y-η)为原始点云数据;将所述二维高斯滤波函数w(x,y)离散化,得到二维离散高斯滤波过程公式:
Figure BDA0003219899730000033
其中,g、k为计算高斯评定基准面w所需的离散计算系数,g的范围为g1~g2,k的范围为k1~k2;其中,i=g1,...,Lx-g2,j=k1,...,Ly-k2,Lx和Ly为采样数据点;Δx、Δy为采样间隔;通过所述二维离散高斯滤波过程公式w(xi,yi),对所述第一点云数据以及所述第二点云数据进行二维高斯滤波。
在一种可行的实施方式中,通过预设算法,对进行高斯滤波后的所述第一点云数据进行曲面拟合,得到所述第一基准平面;以及,对进行高斯滤波后的所述第二点云数据进行曲面拟合,得到所述第二基准平面,具体包括:对于滤波后的点云数据点列Qj,k(j=1,2,...,n;k=1,2,...e),分别按照下标j和k的顺序,构造单增的参数序列
Figure BDA0003219899730000034
以及{pk};根据所述参数序列
Figure BDA0003219899730000035
以及
Figure BDA0003219899730000036
构造B样条基函数
Figure BDA0003219899730000037
根据对n行点云数据的一元函数拟合公式
Figure BDA0003219899730000038
得到每一行点云数据的r个空间点列
Figure BDA0003219899730000039
根据对r列空间点列
Figure BDA00032198997300000310
的一元函数拟合的公式
Figure BDA00032198997300000311
得到中间参数lzk,再对所述中间参数lzk进行求和计算,得到控制顶点ljk;根据
Figure BDA00032198997300000312
分别对所述第一点云数据以及所述第二点云数据进行曲面拟合,得到所述第一基准平面以及所述第二基准平面。
在一种可行的实施方式中,确定所述第一激光雷达与所述第一基准平面之间的第一距离,以及所述第二激光雷达与所述第二基准平面之间的第二距离,具体包括:以所述第一激光雷达为原点,建立第一三维坐标系;基于所述第一三维坐标系的原点坐标,以及所述第一基准平面上的所有点云数据在所述第一三维坐标系中的三维坐标,确定所述第一激光雷达与所述第一基准平面上的所有点云数据的距离,并将得到的所述距离中的最小值确定为所述第一激光雷达与所述第一基准平面之间的第一距离;以所述第二激光雷达为原点,建立第二三维坐标系;基于所述第二三维坐标系的原点坐标,以及所述第二基准平面上的所有点云数据在所述第二三维坐标系中的三维坐标,确定所述第二激光雷达与所述第二基准平面上的所有点云数据的距离,并将得到的所述距离中的最小值确定为所述第二激光雷达与所述第二基准平面之间的第二距离。
在一种可行的实施方式中,基于所述第一距离与预存的第一标准距离,确定第一偏移量,以及基于所述第二距离与预存的第二标准距离,确定第二偏移量,具体包括:在所述综采设备运行之前,确定所述第一激光雷达与所述第一基准平面的第一标准距离,并确定所述第二激光雷达与所述第二基准平面的第二标准距离;将所述第一标准距离以及所述第二标准距离保存在内存中;在所述综采设备运行过程中,将所述第一距离与所述第一标准距离进行对比并求第一差值,将所述第一差值的绝对值确定为所述第一偏移量;将所述第二距离与所述第二标准距离进行对比并求第二差值,将所述第二差值的绝对值确定为所述第二偏移量。
在一种可行的实施方式中,基于所述第一偏移量以及所述第二偏移量,确定所述综采设备发生的偏移类型,具体包括:在所述第一距离大于第一标准距离且所述第二距离小于第二标准距离的情况下,确定所述综采设备的偏移类型为下滑偏移;在所述第一距离小于第一标准距离且所述第二距离大于第二标准距离的情况下,确定所述综采设备的偏移类型为上窜偏移;在所述第一距离等于第一标准距离且所述第二距离等于第二标准距离的情况下,确定所述综采设备没有发生偏移。
在一种可行的实施方式中,在根据所述第一偏移量以及所述第二偏移量,确定所述综采设备发生的偏移类型之后,所述方法还包括:在所述综采设备发生偏移的情况下,通过无线通讯模块将所述第一偏移量、所述第二偏移量以及所述偏移类型发送到所述综采设备的控制中心;在所述综采设备的偏移类型为下滑偏移的情况下,通过所述控制中心控制所述综采设备向所述第一巷道壁移动所述第一偏移量;在所述综采设备的偏移类型为上窜偏移的情况下,通过所述控制中心控制所述综采设备向所述第二巷道壁移动所述第二偏移量;将所述第一偏移量、所述第二偏移量、所述偏移类型以及所述控制中心对所述综采设备调整的结果发送到显示设备,并在所述显示设备中进行显示。
本申请实施例还提供了一种综采工作面设备上窜下滑偏移测量系统,系统包括:激光雷达,包括第一激光雷达以及第二激光雷达,用于实时采集所述综采设备第一巷道壁的第一点云数据以及所述综采设备第二巷道壁的第二点云数据;处理器,用于分别对所述第一点云数据与所述第二点云数据进行滤波处理以及曲面拟合处理,得到第一基准平面与第二基准平面;并确定所述第一激光雷达与所述第一基准平面的第一距离,以及所述第二激光雷达与所述第二基准平面的第二距离;所述处理器还用于基于所述第一距离与预存的第一标准距离,确定第一偏移量;基于所述第二距离与预存的第二标准距离,确定第二偏移量;以及基于所述第一偏移量以及所述第二偏移量,确定所述综采设备发生的偏移类型;综采设备控制中心,用于根据所述第一偏移量、所述第二偏移量以及所述偏移类型对所述综采设备作出调整。
在一种可行的实施方式中,所述第一激光雷达安装于所述综采设备的第一防护板上,所述第二激光雷达安装于所述综采设备的第二防护板上;所述第一防护板与所述第二防护板分别安装于所述综采设备靠近巷道壁的两端,用于保护所述综采设备。
本申请实施例采用的上述至少一个技术方案能够达到以下有益效果:
1.通过激光雷达采集综采设备两侧巷道壁的三维点云数据,并传输至ARM处理器进行处理,得到综采设备距离两侧巷道壁的距离,激光雷达与ARM处理器组合的形式对于测量综采设备偏移量的精度更高,且无需将数据传输至计算机处理,数据处理速度更快,测量效率更高,能够适应综采工作面的工作智能化要求。
2.ARM处理器对激光雷达采集的点云数据进行高斯滤波,将综采工作面产生的粉尘过滤掉,减少粉尘造成的误差。然后再通过曲面拟合,将凹凸不平的煤壁拟合为一张光滑的基准平面,然后通过三维坐标计算激光雷达与基准平面的距离,此方法会避免因煤壁凹凸不平而造成的影响真实距离的测量问题,可以得到较为真实的平面。排除了粉尘和煤壁不平等干扰因素,激光雷达与煤壁距离的测量准确度也更高。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本申请实施例提供的一种综采工作面设备上窜下滑偏移测量方法流程图;
图2为本申请实施例提供的一种综采工作面设备结构示意图;
1、巷道壁;2、第一激光雷达;3第一防护板、;4第一液压支架、;5、第二液压支架;6、第三液压支架;7、第二防护板;8、第二激光雷达;9、刮板运输机;10、底调液压缸;
图3为本申请实施例提供的一种综采工作面设备上窜下滑偏移测量系统结构示意图;
图4为本申请实施例提供的一种综采工作面设备上窜下滑偏移测量设备结构示意图;
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本说明书实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1为本申请实施例提供的一种综采工作面设备上窜下滑偏移测量方法流程图,如图1所示,方法具体包括S101-S106:
S101、具体地,处理器采用ARM处理器。ARM处理器通过安装在综采设备第一防护板上的第一激光雷达,实时采集综采设备第一巷道壁对应的第一点云数据,通过安装在综采设备第二防护板上的第二激光雷达,实时采集综采设备第二巷道壁对应的第二点云数据。
具体地,综采设备用于在采煤工作面进行煤炭采集。在综采设备的左右两侧分别安装一个激光雷达,用于测量综采设备两侧与左右巷道壁的距离。ARM处理器实时接收两个激光雷达传回的点云数据。
作为一种可行的实施方式,如图2所示,综采设备主要包括刮板运输机9、底调液压缸10、第一液压支架4、第二液压支架5以及第三液压支架6。在综采设备的两侧分别安装第一防护板3以及第二防护板7,并在第一防护板3上安装第一激光雷达2,在第二防护板7上安装第二激光雷达8。第一激光雷达2实时采集第一巷道壁的三维点云数据,第二激光雷达实时采集第二巷道壁的三维点云数据。需要说明的是,与第一激光雷达距离更近的一侧巷道壁即为第一巷道壁,与第二激光雷达距离更近的一侧巷道壁即为第二巷道壁。图2中的左右关系并不代表现实中的左右关系,仅用于指示各部件的位置关系。
在一个实施例中,若2为第一激光雷达,则1为第一巷道壁,若2为第二激光雷达,则1为第二巷道壁。在综采设备的两端安装防护板,将激光雷达安装于防护板上,在对综采设备提供保护的同时,方便确定综采设备与巷道壁的距离。
S102、ARM处理器分别对第一点云数据与第二点云数据进行滤波处理以及曲面拟合处理,以得到第一基准平面与第二基准平面。
具体地,ARM处理器在接收到两个激光雷达采集的两侧巷道壁的点云数据后,分别对两侧巷道壁的点云数据进行高斯滤波,去除点云数据中的孤立点。通过预设的曲面拟合算法,对进行高斯滤波后的点云数据进行曲面拟合,得到对应的基准平面。
作为一种可行的实施方式,可以根据
Figure BDA0003219899730000081
得到二维高斯滤波权函数h(x,y);其中,
Figure BDA0003219899730000082
x代表偏离二维高斯滤波权函数中心的点云数据的x坐标,y代表偏离二维高斯滤波权函数中心的点云数据的y坐标,λxc以及λyc代表低通高斯滤波器的截止波长;
进一步地,根据w(x,y)=∫∫z(x-ξ,y-η)h(ξ,η)dξdη,得到二维高斯滤波函数w(x,y);其中,ξ、η为卷积积分所需的微分变量,z(x-ξ,y-η)为原始点云数据;
进一步地,由于采集到的点云数据为离散点,因此,将二维高斯滤波函数w(x,y)离散化,得到二维离散高斯滤波过程公式:
Figure BDA0003219899730000083
其中,g、k为计算高斯评定基准面w所需的离散计算系数,g的范围为g1~g2,k的范围为k1~k2;i=g1,...,Lx-g2,j=k1,...,Ly-k2,Lx和Ly为采样数据点;Δx、Δy为采样间隔;通过公式w(xi,yi),对两侧巷道壁的点云数据进行二维高斯滤波,保留有效三维点云数据,将采集到的粉尘的点云数据过滤掉。
作为一种可行的实施方式,对于高斯滤波后的点云数据点列Qj,k(j=1,2,...,n;k=1,2,...e),分别按照下标j和k的顺序,构造单增的参数序列
Figure BDA0003219899730000091
以及{pk};然后根据参数序列
Figure BDA0003219899730000092
以及{pk},构造B样条基函数
Figure BDA0003219899730000093
{Ak,x(p)};
进一步地,先根据
Figure BDA0003219899730000094
对n行点云数据分别作一元函数拟合,在上述公式中代入{Az,x(pk)}以及Qj,k,得到每一行点云数据的r个空间点列
Figure BDA0003219899730000095
进一步地,根据
Figure BDA0003219899730000096
对r列空间点列
Figure BDA0003219899730000097
进行一元函数拟合,在上述公式中代入
Figure BDA0003219899730000098
以及{Az,c(mj)},得到中间参数lzk,再对中间参数lzk进行求和计算,得到点云数据的控制顶点ljk
最后,根据
Figure BDA0003219899730000099
分别对两侧巷道壁的点云数据进行曲面拟合,得到第一基准平面以及第二基准平面。
S103、ARM处理器确定第一激光雷达2与第一基准平面之间的第一距离,以及第二激光雷达8与第二基准平面之间的第二距离。
具体地,分别以第一激光雷达2为原点,建立第一三维坐标系,以及以第二激光雷达8为原点,建立第二三维坐标系。在第一三维坐标系中,确定第一基准平面上的所有拟合后的点云数据的三维坐标,在第二三维坐标系中,确定第二基准平面上的所有拟合后的点云数据的三维坐标。
在一个实施例中,第一三维坐标系与第二三维坐标系的坐标轴可以与大地坐标系的坐标轴相同,也可以指向任意三个互相垂直的方向,可按照实际需求进行设置,本申请实施例对此不作限定。
进一步地,基于第一三维坐标系的原点坐标(0,0,0),以及第一基准平面上的所有点云数据在第一三维坐标系中的三维坐标,计算第一激光雷达2与第一基准平面上的所有点云数据的坐标距离,得到的所有坐标距离中的最小值即为第一激光雷达2与第一基准平面之间的第一距离。基于第二三维坐标系的原点坐标(0,0,0),以及第二基准平面上的所有点云数据在第二三维坐标系中的三维坐标,计算第二激光雷达8与第二基准平面上的所有点云数据的坐标距离,得到的所有坐标距离中的最小值即为第二激光雷达8与第二基准平面之间的第二距离。若需要得到激光雷达与基准平面的实际距离,则用计算出的坐标距离乘以三维坐标系的单位长度即可。
S104、ARM处理器基于第一距离与预存的第一标准距离,确定第一偏移量,以及基于第二距离与预存的第二标准距离,确定第二偏移量。
具体地,在综采设备开始运行之前,ARM处理器先根据S101-S103所示的过程确定第一激光雷达2与第一基准平面的第一标准距离,以及确定第二激光雷达8与第二基准平面的第二标准距离,并将第一标准距离和第二标准距离保存在ARM处理器的内存中。
进一步地,在综采设备运行过程中,ARM处理器实时将计算出的第一距离与内存中的第一标准距离进行对比并计算两者的第一差值,第一差值的绝对值即为第一偏移量。同时,ARM处理器也实时将计算出的第二距离与第二标准距离进行对比并求两者的第二差值,第二差值的绝对值即为第二偏移量。
在一个实施例中,若计算出的第一距离为12cm,第一标准距离为10cm,则第一偏移量为2cm。
S105、ARM处理器根据第一偏移量以及第二偏移量,确定综采设备发生的偏移类型。
具体地,在第一距离大于第一标准距离且第二距离小于第二标准距离的情况下,综采设备的偏移类型为下滑偏移;在第一距离小于第一标准距离且第二距离大于第二标准距离的情况下,确定综采设备的偏移类型为上窜偏移;在第一距离等于第一标准距离且第二距离等于第二标准距离的情况下,确定综采设备没有发生偏移。
在一个实施例中,如图2所示,若计算出的第一距离为a,第二距离为b,预存的第一标准距离为A,第二标准距离为B。那么若a>A且b<B,则判断综采设备发生了下滑偏移。若a<A且b>B,判断设备发生上窜偏移。同样,第一距离也可以为b,第二距离也可以为a,因此,上窜偏移和下滑偏移具体是向哪个方向偏移,取决于工作人员将哪个方向的激光雷达设为第一激光雷达。
S106、ARM处理器在综采设备发生偏移的情况下,将第一偏移量、第二偏移量以及偏移类型发送到综采设备的控制中心,以使控制中心实时调整综采设备。
具体地,在综采设备发生偏移的情况下,通过无线通讯模块将第一偏移量、第二偏移量以及偏移类型发送到综采设备的控制中心。在综采设备的偏移类型为下滑偏移的情况下,控制中心控制综采设备向第一巷道壁移动第一偏移量。在综采设备的偏移类型为上窜偏移的情况下,控制中心控制综采设备向第二巷道壁移动第二偏移量。
在一个实施例中,若控制中心接收到的信息为综采设备产生下滑偏移,且第一偏移量为3cm,则控制综采设备向第一巷道壁移动3cm。若控制中心接收到的信息为综采设备产生上窜偏移,且第二偏移量为3cm,则控制综采设备向第二巷道壁移动3cm。
需要说明的是,一般情况下,综采设备若发生偏移,其第一偏移量与第二偏移量是相同的,例如综采设备的左侧向左偏移了3cm,右侧也会向左偏移3cm。本申请为避免特殊情况发生,在调整综采设备时,调整下滑偏移以第一偏移量为准,调整上窜偏移以第二偏移量为准。
进一步地,ARM处理器将实时计算出的第一偏移量、第二偏移量、偏移类型以及控制中心对综采设备调整的结果发送到显示设备,并在显示设备中进行显示,以便工作人员查看综采设备状态。
图3为本申请实施例提供的一种综采工作面设备上窜下滑偏移测量系统结构示意图,如图3所示,系统包括:
第一激光雷达310以及第二激光雷达320,用于实时采集综采设备第一巷道壁的第一点云数据以及综采设备第二巷道壁的第二点云数据。
作为一种可行的实施方式,第一激光雷达310安装于综采设备的第一防护板上,第二激光雷达320安装于综采设备的第二防护板上。第一防护板与第二防护板分别安装于综采设备靠近巷道壁的两端,用于保护综采设备不被巷道壁磨损。
ARM处理器330,用于分别对第一点云数据与第二点云数据进行滤波处理以及曲面拟合处理,得到第一基准平面与第二基准平面;并确定第一激光雷达与第一基准平面的第一距离,以及第二激光雷达与第二基准平面的第二距离;ARM处理器还用于基于第一距离与预存的第一标准距离,确定第一偏移量;基于第二距离与预存的第二标准距离,确定第二偏移量;以及基于第一偏移量以及第二偏移量,确定综采设备发生的偏移类型。
其中,分别对第一点云数据与第二点云数据进行滤波处理以及曲面拟合处理,得到第一基准平面与第二基准平面,具体包括:对第一点云数据进行高斯滤波,去除第一点云数据中的孤立点;以及,对第二点云数据进行高斯滤波,去除第二点云数据中的孤立点;通过预设算法,对进行高斯滤波后的第一点云数据进行曲面拟合,得到第一基准平面;以及,对进行高斯滤波后的第二点云数据进行曲面拟合,得到第二基准平面。
其中,对第一点云数据进行高斯滤波,去除第一点云数据中的孤立点;以及,对第二点云数据进行高斯滤波,去除第二点云数据中的孤立点,具体包括:根据
Figure BDA0003219899730000121
得到二维高斯滤波权函数h(x,y);其中,
Figure BDA0003219899730000122
x代表偏离二维高斯滤波权函数中心的点云数据的x坐标,y代表偏离二维高斯滤波权函数中心的点云数据的y坐标,λxc以及λyc代表低通高斯滤波器的截止波长;根据w(x,y)=∫∫z(x-ξ,y-η)h(ξ,η)dξdη,得到二维高斯滤波函数w(x,y);其中,ξ、η为卷积积分所需的微分变量,z(x-ξ,y-η)为原始点云数据;将二维高斯滤波函数w(x,y)离散化,得到二维离散高斯滤波过程公式:
Figure BDA0003219899730000131
其中,g、k为计算高斯评定基准面w所需的离散计算系数,g的范围为g1~g2,k的范围为k1~k2;其中,i=g1,...,Lx-g2,j=k1,...,Ly-k2,Lx和Ly为采样数据点;Δx、Δy为采样间隔;通过二维离散高斯滤波过程公式w(xi,yi),对第一点云数据以及第二点云数据进行二维高斯滤波。
通过预设算法,对进行高斯滤波后的第一点云数据进行曲面拟合,得到第一基准平面;以及,对进行高斯滤波后的第二点云数据进行曲面拟合,得到第二基准平面,具体包括:对于滤波后的点云数据点列Qj,k(j=1,2,...,n;k=1,2,...e),分别按照下标j和k的顺序,构造单增的参数序列
Figure BDA0003219899730000132
以及{pk};根据参数序列
Figure BDA0003219899730000133
以及{pk},构造B样条基函数
Figure BDA0003219899730000134
根据对n行点云数据的一元函数拟合公式
Figure BDA0003219899730000135
得到每一行点云数据的r个空间点列
Figure BDA0003219899730000136
根据对r列空间点列
Figure BDA0003219899730000137
的一元函数拟合的公式
Figure BDA0003219899730000138
得到中间参数lzk,再对中间参数lzk进行求和计算,得到控制顶点ljk;根据
Figure BDA0003219899730000139
分别对第一点云数据以及第二点云数据进行曲面拟合,得到第一基准平面以及第二基准平面。
其中,确定第一激光雷达与第一基准平面之间的第一距离,以及第二激光雷达与第二基准平面之间的第二距离,具体包括:以第一激光雷达为原点,建立第一三维坐标系;基于第一三维坐标系的原点坐标,以及第一基准平面上的所有点云数据在第一三维坐标系中的三维坐标,确定第一激光雷达与第一基准平面上的所有点云数据的距离,并将得到的距离中的最小值确定为第一激光雷达与第一基准平面之间的第一距离;以第二激光雷达为原点,建立第二三维坐标系;基于第二三维坐标系的原点坐标,以及第二基准平面上的所有点云数据在第二三维坐标系中的三维坐标,确定第二激光雷达与第二基准平面上的所有点云数据的距离,并将得到的距离中的最小值确定为第二激光雷达与第二基准平面之间的第二距离。
其中,基于第一距离与预存的第一标准距离,确定第一偏移量,以及基于第二距离与预存的第二标准距离,确定第二偏移量,具体包括:在综采设备运行之前,确定第一激光雷达与第一基准平面的第一标准距离,并确定第二激光雷达与第二基准平面的第二标准距离;将第一标准距离以及第二标准距离保存在内存中;在综采设备运行过程中,将第一距离与第一标准距离进行对比并求第一差值,将第一差值的绝对值确定为第一偏移量;将第二距离与第二标准距离进行对比并求第二差值,将第二差值的绝对值确定为第二偏移量。
其中,基于第一偏移量以及第二偏移量,确定综采设备发生的偏移类型,具体包括:在第一距离大于第一标准距离且第二距离小于第二标准距离的情况下,确定综采设备的偏移类型为下滑偏移;在第一距离小于第一标准距离且第二距离大于第二标准距离的情况下,确定综采设备的偏移类型为上窜偏移;在第一距离等于第一标准距离且第二距离等于第二标准距离的情况下,确定综采设备没有发生偏移。
作为一种可行的实施方式,ARM处理器330包括数据处理模块以及通讯模块,激光雷达采用对应通讯模块与ARM处理器330进行数据通信,增强数据传输中的抗干扰能力。
综采设备控制中心340,用于根据第一偏移量、第二偏移量以及偏移类型对综采设备作出调整。
图4为本申请实施例提供的一种综采工作面设备上窜下滑偏移测量设备结构示意图,如图4所示,设备包括:
至少一个处理器;以及,与至少一个处理器通信连接的存储器;其中,存储器存储有可被至少一个处理器执行的指令,指令被至少一个处理器执行,以使至少一个处理器能够执行S101-S106中的任一步骤。
本申请中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置、设备、非易失性计算机存储介质实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
上述对本申请特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请的实施例可以有各种更改和变化。凡在本申请实施例的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。

Claims (10)

1.一种综采工作面设备上窜下滑偏移测量方法,其特征在于,所述方法包括:
接收第一激光雷达采集的第一巷道壁对应的第一点云数据;以及,接收第二激光雷达采集的第二巷道壁对应的第二点云数据;
分别对所述第一点云数据与所述第二点云数据进行滤波处理以及曲面拟合处理,以得到第一基准平面与第二基准平面;
确定所述第一激光雷达与所述第一基准平面之间的第一距离,以及所述第二激光雷达与所述第二基准平面之间的第二距离;
基于所述第一距离与预存的第一标准距离,确定第一偏移量,以及基于所述第二距离与预存的第二标准距离,确定第二偏移量;
根据所述第一偏移量以及所述第二偏移量,确定所述综采设备发生的偏移类型。
2.根据权利要求1所述的一种综采工作面设备上窜下滑偏移测量方法,其特征在于,分别对所述第一点云数据与所述第二点云数据进行滤波处理以及曲面拟合处理,以得到第一基准平面与第二基准平面,具体包括:
对所述第一点云数据进行高斯滤波,去除所述第一点云数据中的孤立点;以及,对所述第二点云数据进行高斯滤波,去除所述第二点云数据中的孤立点;
通过预设算法,对进行高斯滤波后的所述第一点云数据进行曲面拟合,得到所述第一基准平面;以及,对进行高斯滤波后的所述第二点云数据进行曲面拟合,得到所述第二基准平面。
3.根据权利要求2所述的一种综采工作面设备上窜下滑偏移测量方法,其特征在于,对所述第一点云数据进行高斯滤波,去除所述第一点云数据中的孤立点;以及,对所述第二点云数据进行高斯滤波,去除所述第二点云数据中的孤立点,具体包括:
根据
Figure FDA0003219899720000021
得到二维高斯滤波权函数h(x,y);其中,
Figure FDA0003219899720000022
x代表偏离二维高斯滤波权函数中心的点云数据的x坐标,y代表偏离二维高斯滤波权函数中心的点云数据的y坐标,λxc以及λyc代表低通高斯滤波器的截止波长;
根据w(x,y)=∫∫z(x-ξ,y-η)h(ξ,η)dξdη,得到二维高斯滤波函数w(x,y);其中,ξ、η为卷积积分所需的微分变量,z(x-ξ,y-η)为原始点云数据;
将所述二维高斯滤波函数w(x,y)离散化,得到二维离散高斯滤波过程公式:
Figure FDA0003219899720000023
其中,g、k为计算高斯评定基准面w所需的离散计算系数,g的范围为g1~g2,k的范围为k1~k2;其中,i=g1,...,Lx-g2,j=k1,...,Ly-k2,Lx和Ly为采样数据点;Δx、Δy为采样间隔;
通过所述二维离散高斯滤波过程公式w(xi,yi),对所述第一点云数据以及所述第二点云数据进行二维高斯滤波。
4.根据权利要求2所述的一种综采工作面设备上窜下滑偏移测量方法,其特征在于,通过预设算法,对进行高斯滤波后的所述第一点云数据进行曲面拟合,得到所述第一基准平面;以及,对进行高斯滤波后的所述第二点云数据进行曲面拟合,得到所述第二基准平面,具体包括:
对于滤波后的点云数据点列Qj,k(j=1,2,...,n;k=1,2,...e),分别按照下标j和k的顺序,构造单增的参数序列{mj}以及{pk};
根据所述参数序列{mj}以及{pk},构造B样条基函数{Aj,c(m)}、{Ak,x(p)};
根据对n行点云数据的一元函数拟合公式
Figure FDA0003219899720000024
得到每一行点云数据的r个空间点列
Figure FDA0003219899720000025
根据对r列空间点列
Figure FDA0003219899720000031
的一元函数拟合的公式
Figure FDA0003219899720000032
得到中间参数lzk,再对所述中间参数lzk进行求和计算,得到控制顶点ljk
根据
Figure FDA0003219899720000033
分别对所述第一点云数据以及所述第二点云数据进行曲面拟合,得到所述第一基准平面以及所述第二基准平面。
5.根据权利要求1所述的一种综采工作面设备上窜下滑偏移测量方法,其特征在于,确定所述第一激光雷达与所述第一基准平面之间的第一距离,以及所述第二激光雷达与所述第二基准平面之间的第二距离,具体包括:
以所述第一激光雷达为原点,建立第一三维坐标系;
基于所述第一三维坐标系的原点坐标,以及所述第一基准平面上的所有点云数据在所述第一三维坐标系中的三维坐标,确定所述第一激光雷达与所述第一基准平面上的所有点云数据的距离,并将得到的所述距离中的最小值确定为所述第一激光雷达与所述第一基准平面之间的第一距离;
以所述第二激光雷达为原点,建立第二三维坐标系;
基于所述第二三维坐标系的原点坐标,以及所述第二基准平面上的所有点云数据在所述第二三维坐标系中的三维坐标,确定所述第二激光雷达与所述第二基准平面上的所有点云数据的距离,并将得到的所述距离中的最小值确定为所述第二激光雷达与所述第二基准平面之间的第二距离。
6.根据权利要求1所述的一种综采工作面设备上窜下滑偏移测量方法,其特征在于,基于所述第一距离与预存的第一标准距离,确定第一偏移量,以及基于所述第二距离与预存的第二标准距离,确定第二偏移量,具体包括:
在所述综采设备运行之前,确定所述第一激光雷达与所述第一基准平面的第一标准距离,并确定所述第二激光雷达与所述第二基准平面的第二标准距离;
将所述第一标准距离以及所述第二标准距离保存在内存中;
在所述综采设备运行过程中,将所述第一距离与所述第一标准距离进行对比并求第一差值,将所述第一差值的绝对值确定为所述第一偏移量;
将所述第二距离与所述第二标准距离进行对比并求第二差值,将所述第二差值的绝对值确定为所述第二偏移量。
7.根据权利要求6所述的一种综采工作面设备上窜下滑偏移测量方法,其特征在于,基于所述第一偏移量以及所述第二偏移量,确定所述综采设备发生的偏移类型,具体包括:
在所述第一距离大于第一标准距离且所述第二距离小于第二标准距离的情况下,确定所述综采设备的偏移类型为下滑偏移;
在所述第一距离小于第一标准距离且所述第二距离大于第二标准距离的情况下,确定所述综采设备的偏移类型为上窜偏移;
在所述第一距离等于第一标准距离且所述第二距离等于第二标准距离的情况下,确定所述综采设备没有发生偏移。
8.根据权利要求1所述的一种综采工作面设备上窜下滑偏移测量方法,其特征在于,在根据所述第一偏移量以及所述第二偏移量,确定所述综采设备发生的偏移类型之后,所述方法还包括:
在所述综采设备发生偏移的情况下,通过无线通讯模块将所述第一偏移量、所述第二偏移量以及所述偏移类型发送到所述综采设备的控制中心;
在所述综采设备的偏移类型为下滑偏移的情况下,通过所述控制中心控制所述综采设备向所述第一巷道壁移动所述第一偏移量;
在所述综采设备的偏移类型为上窜偏移的情况下,通过所述控制中心控制所述综采设备向所述第二巷道壁移动所述第二偏移量;
将所述第一偏移量、所述第二偏移量、所述偏移类型以及所述控制中心对所述综采设备调整的结果发送到显示设备,并在所述显示设备中进行显示。
9.一种综采工作面设备上窜下滑偏移测量系统,其特征在于,所述系统包括:
激光雷达,包括第一激光雷达以及第二激光雷达,分别用于实时采集所述第一巷道壁的第一点云数据,以及所述第二巷道壁的第二点云数据;
处理器,用于分别对所述第一点云数据与所述第二点云数据进行滤波处理以及曲面拟合处理,得到第一基准平面与第二基准平面;并确定所述第一激光雷达与所述第一基准平面的第一距离,以及所述第二激光雷达与所述第二基准平面的第二距离;所述处理器还用于基于所述第一距离与预存的第一标准距离,确定第一偏移量;基于所述第二距离与预存的第二标准距离,确定第二偏移量;以及基于所述第一偏移量以及所述第二偏移量,确定所述综采设备发生的偏移类型;
综采设备控制中心,用于根据所述第一偏移量、所述第二偏移量以及所述偏移类型对所述综采设备作出调整。
10.根据权利要求9所述的一种综采工作面设备上窜下滑偏移测量系统,其特征在于,
所述第一激光雷达安装于所述综采设备的第一防护板上,所述第二激光雷达安装于所述综采设备的第二防护板上;
所述第一防护板与所述第二防护板分别安装于所述综采设备靠近巷道壁的两端,用于保护所述综采设备。
CN202110955008.3A 2021-08-19 2021-08-19 一种综采工作面设备上窜下滑偏移测量方法及系统 Active CN113686251B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202110955008.3A CN113686251B (zh) 2021-08-19 2021-08-19 一种综采工作面设备上窜下滑偏移测量方法及系统
PCT/CN2021/135365 WO2023019808A1 (zh) 2021-08-19 2021-12-03 一种综采工作面设备上窜下滑偏移测量方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110955008.3A CN113686251B (zh) 2021-08-19 2021-08-19 一种综采工作面设备上窜下滑偏移测量方法及系统

Publications (2)

Publication Number Publication Date
CN113686251A true CN113686251A (zh) 2021-11-23
CN113686251B CN113686251B (zh) 2022-12-13

Family

ID=78580708

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110955008.3A Active CN113686251B (zh) 2021-08-19 2021-08-19 一种综采工作面设备上窜下滑偏移测量方法及系统

Country Status (2)

Country Link
CN (1) CN113686251B (zh)
WO (1) WO2023019808A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023019808A1 (zh) * 2021-08-19 2023-02-23 山东科技大学 一种综采工作面设备上窜下滑偏移测量方法及系统

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0479273A2 (en) * 1990-10-03 1992-04-08 Aisin Seiki Kabushiki Kaisha Laser radar for a vehicle lateral guidance system
CN101961698A (zh) * 2010-08-04 2011-02-02 中国科学院自动化研究所 一种嵌入式喷枪位姿实时测量装置及方法
JP2011191238A (ja) * 2010-03-16 2011-09-29 Daihatsu Motor Co Ltd 運転支援装置
WO2017032121A1 (zh) * 2015-08-25 2017-03-02 太原理工大学 一种煤矿综采工作面大型装备集中控制平台的实现方法
CN109255837A (zh) * 2018-08-06 2019-01-22 上海大学 一种用于激光雷达点云数据处理的高效b样条曲面的构造方法
CN109373921A (zh) * 2018-10-26 2019-02-22 南京航空航天大学 一种隧道监测方法及装置
CN110658528A (zh) * 2019-10-12 2020-01-07 山东科技大学 基于激光雷达的综采工作面成套设备偏移监测方法
CN110703266A (zh) * 2019-10-12 2020-01-17 山东科技大学 一种掘进机精准定位及导航系统
CN110986788A (zh) * 2019-11-15 2020-04-10 华南农业大学 一种基于三维点云牲畜表型体尺数据的自动测量方法
CN111612902A (zh) * 2020-04-20 2020-09-01 杭州鼎控自动化技术有限公司 一种基于雷达点云数据的煤矿巷道三维模型构建方法
WO2021082971A1 (zh) * 2019-10-28 2021-05-06 北京天地玛珂电液控制系统有限公司 一种综采工作面数字化开采模型构建系统及方法
CN113128248A (zh) * 2019-12-26 2021-07-16 深圳一清创新科技有限公司 障碍物检测方法、装置、计算机设备和存储介质

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001271583A (ja) * 2000-03-24 2001-10-05 Toda Constr Co Ltd トンネル内壁測定システムおよびトンネル内壁測定方法
CN103899358B (zh) * 2014-03-19 2016-03-02 北京天地玛珂电液控制系统有限公司 一种运输机上窜下滑自动控制系统
CN113686251B (zh) * 2021-08-19 2022-12-13 山东科技大学 一种综采工作面设备上窜下滑偏移测量方法及系统

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0479273A2 (en) * 1990-10-03 1992-04-08 Aisin Seiki Kabushiki Kaisha Laser radar for a vehicle lateral guidance system
JP2011191238A (ja) * 2010-03-16 2011-09-29 Daihatsu Motor Co Ltd 運転支援装置
CN101961698A (zh) * 2010-08-04 2011-02-02 中国科学院自动化研究所 一种嵌入式喷枪位姿实时测量装置及方法
WO2017032121A1 (zh) * 2015-08-25 2017-03-02 太原理工大学 一种煤矿综采工作面大型装备集中控制平台的实现方法
CN109255837A (zh) * 2018-08-06 2019-01-22 上海大学 一种用于激光雷达点云数据处理的高效b样条曲面的构造方法
CN109373921A (zh) * 2018-10-26 2019-02-22 南京航空航天大学 一种隧道监测方法及装置
CN110658528A (zh) * 2019-10-12 2020-01-07 山东科技大学 基于激光雷达的综采工作面成套设备偏移监测方法
CN110703266A (zh) * 2019-10-12 2020-01-17 山东科技大学 一种掘进机精准定位及导航系统
WO2021082971A1 (zh) * 2019-10-28 2021-05-06 北京天地玛珂电液控制系统有限公司 一种综采工作面数字化开采模型构建系统及方法
CN110986788A (zh) * 2019-11-15 2020-04-10 华南农业大学 一种基于三维点云牲畜表型体尺数据的自动测量方法
CN113128248A (zh) * 2019-12-26 2021-07-16 深圳一清创新科技有限公司 障碍物检测方法、装置、计算机设备和存储介质
CN111612902A (zh) * 2020-04-20 2020-09-01 杭州鼎控自动化技术有限公司 一种基于雷达点云数据的煤矿巷道三维模型构建方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
孔祥斌: "综采工作面刮板输送机上窜下滑的问题研究", 《中国设备工程》 *
张晓等: "综采工作面设备上窜下滑原因分析及预防措施", 《科技创新与应用》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023019808A1 (zh) * 2021-08-19 2023-02-23 山东科技大学 一种综采工作面设备上窜下滑偏移测量方法及系统

Also Published As

Publication number Publication date
WO2023019808A1 (zh) 2023-02-23
CN113686251B (zh) 2022-12-13

Similar Documents

Publication Publication Date Title
US11773544B2 (en) Determining milled volume or milled area of a milled surface
CN111442759B (zh) 一种综采工作面设备位姿统一监测系统
Simmons et al. Obstacle avoidance and safeguarding for a lunar rover
US8462208B2 (en) System and method for tridimensional cartography of a structural surface
AU2016202377B2 (en) Improved mining machine and method
CN110703266A (zh) 一种掘进机精准定位及导航系统
CN113686251B (zh) 一种综采工作面设备上窜下滑偏移测量方法及系统
CN110658528A (zh) 基于激光雷达的综采工作面成套设备偏移监测方法
CN110067566A (zh) 盾构纠偏力矩的预测方法及系统
JP7221549B2 (ja) 情報処理装置および移動ロボット
KR102095799B1 (ko) 전자기파 분석을 통한 지하구조물과 공동형상 맵핑장치 및 그 맵핑방법
CN108827220B (zh) 一种基于捷联惯导的煤矿综采工作面直线度检测方法
CN108267116A (zh) 一种基于路面三维数据的路面错台检测方法及装置
KR102619960B1 (ko) 땅 밀림 모니터링과 암석 균열 측정 및 분석 시스템
CN113944494B (zh) 一种基于超声无线测距的液压支架自动调直方法及系统
CN212363185U (zh) 一种综采工作面设备位姿监测系统
CN113772363A (zh) 一种刮板输送机位姿模型建立方法和系统
CN114593730A (zh) 一种定位导航系统
CN114662259B (zh) 一种基于高精度惯导的采煤机三维定位方法及系统
Dong et al. A Weld Line Detection Method Based on 3D Point Cloud for Automatic NDT
CN115235412A (zh) 一种三维形变监测装置及监测方法
CN115291232A (zh) 一种用于井下水仓自动清理机器人的位姿识别方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant