CN108957111A - 一种可滤除交流采样中残余直流分量的馈线终端计量方法 - Google Patents

一种可滤除交流采样中残余直流分量的馈线终端计量方法 Download PDF

Info

Publication number
CN108957111A
CN108957111A CN201810597330.1A CN201810597330A CN108957111A CN 108957111 A CN108957111 A CN 108957111A CN 201810597330 A CN201810597330 A CN 201810597330A CN 108957111 A CN108957111 A CN 108957111A
Authority
CN
China
Prior art keywords
component
sampling
value
filter
residual
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.)
Withdrawn
Application number
CN201810597330.1A
Other languages
English (en)
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.)
Qingdao Topscomm Communication Co Ltd
Original Assignee
Qingdao Topscomm Communication Co Ltd
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 Qingdao Topscomm Communication Co Ltd filed Critical Qingdao Topscomm Communication Co Ltd
Priority to CN201810597330.1A priority Critical patent/CN108957111A/zh
Publication of CN108957111A publication Critical patent/CN108957111A/zh
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/25Arrangements for measuring currents or voltages or for indicating presence or sign thereof using digital measurement techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/25Arrangements for measuring currents or voltages or for indicating presence or sign thereof using digital measurement techniques
    • G01R19/2506Arrangements for conditioning or analysing measured signals, e.g. for indicating peak values ; Details concerning sampling, digitizing or waveform capturing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Current Or Voltage (AREA)

Abstract

本发明公开了一种可滤除交流采样中残余直流分量的馈线终端计量方法,包括以下步骤,步骤一,建立直流滤波器差分方程;步骤二,根据滤波器原理及差分方程推导出滤波算法;步骤三,采用单片机进行数字化处理,根据采样定理对所需要采集的交流模拟量进行等间隔采样,从而获得一个离散信号序列;步骤四,在单片机的AD中断中部署交流采样及信号处理、滤波算法,当采集的数据量累积到一定数量时,程序前端会根据累积的缓存数据进行周期任务或秒任务计算,计算电压和电流的有效值、电网频率、功率参数及相关电气参数;在测量前完成对交采得到的每相的电压、电流进行残余直流滤波,以消除直流分量对测量的影响,实现AD中断计量。

Description

一种可滤除交流采样中残余直流分量的馈线终端计量方法
技术领域
本发明涉及配电网自动化系统领域,涉及一种可滤除交流采样中残余直流分量的馈线终端计量方法。
背景技术
随着科学技术的发展,伴随着配电网自动化技术的不断革新,不论是从硬件设计还是从算法创新,多种多样的配电终端设备都处在不断地更新换代中。近几年,随着国家电网智能化水平的不断提升,人们对配电设备的要求越来越高,除了硬件上的不断完善,配电设备在数据采集及控制算法上的创新与优化也越来越受到人们的重视。在智能电表及配电终端领域,馈线终端(Feeder Terminal Unit,FTU)作为一种集采集、监测和控制于一体的配电自动化设备,被广泛地应用于配电网自动化系统的升级和改造中来。数据采集作为FTU的一项基本功能,其采集的精度将直接影响到配电网参数的测量,甚至影响FTU对配电线路上相关设备的控制。因此,提高数据采集的精度显得尤为重要。
目前,针对馈电线路上三相交流电的数据采集,一般的做法是利用交流采样技术将电网二次侧的电压和电流经高精度的PT或CT变换成小电流,再通过FTU中的采样电路进行AD采样并数字化,然后通过MCU进行数字化处理。这种交流采样方法是按照一定的采样频率对施加于电网二次侧的电压和电流的瞬时值进行采样,具有实时性好、方便处理的优点。但是这种数据采集方式要求采样电路有很高的稳定性,并且要求互感器有较高的精度,在不同的采集环境下或在有外界电磁干扰的情况下,都有可能会对采集的数据造成影响,使采样数据中引入误差,从而影响对配电网参数的测量精度或影响FTU对馈电线路上相关设备的控制判断等。
残余直流分量作为一种在交采前端由传感器或板级电路引入的采样误差,其对数据采集精度的影响不容忽视,需要采取一定的手段将其滤除。而且在实际应用中,残余直流不是一成不变的,而是随着数据的采集而不断变化的。在进行数字化信号处理时,信号中的直流分量通常需要去除。在配电网自动化应用中,FTU在进行数据采集时,第一步是采集经过PT或CT转换而来的小电流,经采样电路进入单片机。但是,在交采的前端会因为传感器、采样电路或外界电磁干扰等因素的存在,引入残余直流分量,使交采数据存在误差。由于采用模拟电路去除直流分量则比较复杂,目前,通常的做法是在交采后端也即单片机中通过算法进行数字化滤波,去除直流分量。常规去直流的方法是对单片机采集到的数字量进行累加,然后将累加值进行移位操作即可获得直流分量,如累加8192个瞬时值数据点,则直流分量可由累加总和右移13位得到。但是根据该算法原理便可看出该方法所得到的直流分量为常数,而且需要有足够的采样点进行计算得到。并且在实际配电网应用中,每个瞬时采样点对应的直流分量并不一致,如果按照一个常数去滤除会引入新的误差。现有的方法无法解决上述问题,因此,需要发展一种自适应的滤波算法,能够根据当前采样值进行计算,得到当前值对应的直流分量并滤除,提高采样精度,同时,发展和实现配电终端的软件计量方法也急需解决。
发明内容
本发明针对上述问题,克服现有技术的不足,提出一种可滤除交流采样中残余直流分量的馈线终端计量方法,通过FTU内置单片机的滤波算法根据采集到的电压、电流瞬时值实现自适应的滤除残余直流分量,使交采得到的电压、电流值更加准确,确保配电网参数的测量有更高的精度,同时减少硬件设备的投资、省略实际的滤波电路,能够实现馈线终端的软件计量,充分发挥计算机算法具有的计算高精度的优点。
本发明方法通过FTU内置单片机的滤波算法可根据单片机采集到的已量化的瞬时值实现自适应的数字化滤波,消除在电网信号采集中可能出现的直流分量干扰,为测量电网参数之前提供一个准确的电压、电流瞬时值,从而提升测量精度。该方法以直流滤波器为基础,采用差分方程计算每个瞬时采样点所对应的直流分量,然后进行数字化滤波,通过选择不同的滤波器系数实现在不同情况下对不同数据的残余直流滤波,实现馈线终端的软件计量。
为了实现上述发明目的,本发明采取如下技术方案:
一种可滤除交流采样中残余直流分量的馈线终端计量方法,包括以下步骤,
步骤一,建立直流滤波器差分方程;
步骤二,根据滤波器原理及差分方程推导出滤波算法;
步骤三,采用单片机进行数字化处理,根据采样定理对所需要采集的交流模拟量进行等间隔采样,从而获得一个离散信号序列;
步骤四,在单片机的AD中断中部署交流采样及信号处理、滤波算法,当采集的数据量累积到一定数量时,程序前端会根据累积的缓存数据进行周期任务或秒任务计算,计算电压和电流的有效值、电网频率、功率参数及相关电气参数;在测量前完成对交采得到的每相的电压、电流进行残余直流滤波,以消除直流分量对测量的影响,实现AD中断计量。
步骤一中建立直流滤波器差分方程为,
y[n]=(x[n]-y[n-1])×K+y[n-1],K=1/2m
其中,y[n]为当前采样点对应的直流分量输出,y[n-1]为上一采样点对应的直流分量输出,x[n]为当前瞬时采样值输入,K为滤波器系数,m值则根据实际情况进行选择。m值越大,K值越小,直流分量收敛越慢,但是波动较小,稳定性较高;m值越小,K值越大,直流分量收敛的越快,但是波动比较大,稳定性较差。
步骤二中根据滤波器原理及差分方程推导出滤波算法为,
yn=((xn-yn-1)>>m+yn-1
xn=xn-yn
其中,m为滤波器系数中的指数值,根据实际应用选择合适的m值。yn和yn-1分别代表当前和上一次计算得到的直流分量,xn为本次采集瞬时值。
步骤四中AD中断计量处理过程为,读取三相电压值和三相电流值,对每一相进行残余直流滤波,电压值、电流值数据量累积,周期任务或秒任务计算。采集的数据量累积到一定数量,该数量一般为1个工频周期或1秒对应的数据量。
本发明的有益效果是:本发明一种可滤除交流采样中残余直流分量的馈线终端计量方法,能够根据瞬时采样值自适应地进行直流滤波,消除残余直流对交采精度的影响,使交采得到的电压、电流值更加准确,提高FTU对配电网参数的测量精度,确保FTU对配电网线路上相关设备控制的准确度,减少硬件设备的投资、省略实际的滤波电路,能够实现馈线终端的软件计量,充分发挥计算机算法具有的计算高精度的优点。
附图说明
图1为本发明AD中断计量流程图。
图2为含有直流分量的滤波之前波形。
图3为本发明K=1/214时直流滤波直流分量收敛波形。
图4为本发明K=1/214时直流滤波后波形。
图5为本发明K=1/212时直流滤波直流分量收敛波形。
图6为本发明K=1/212时直流滤波后波形。
图7为本发明K=1/28时直流滤波直流分量收敛波形。
图8为本发明K=1/28时直流滤波后波形。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明,以具体阐述本发明的技术方案。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
一种可滤除交流采样中残余直流分量的馈线终端计量方法,包括以下步骤,
步骤一,建立直流滤波器差分方程,建立直流滤波器差分方程为,
y[n]=(x[n]-y[n-1])×K+y[n-1],K=1/2m
其中,y[n]为当前采样点对应的直流分量输出,y[n-1]为上一采样点对应的直流分量输出,x[n]为当前瞬时采样值输入,K为滤波器系数,m值则根据实际情况进行选择。m值越大,K值越小,直流分量收敛越慢,但是波动较小,稳定性较高;m值越小,K值越大,直流分量收敛的越快,但是波动比较大,稳定性较差。
步骤二,根据滤波器原理及差分方程推导出滤波算法为,
yn=((xn-yn-1)>>m+yn-1
xn=xn-yn
其中,m为滤波器系数中的指数值,根据实际应用选择合适的m值。yn和yn-1分别代表当前和上一次计算得到的直流分量,xn为本次采集瞬时值。
步骤三,采用单片机进行数字化处理,根据采样定理对所需要采集的交流模拟量进行等间隔采样,从而获得一个离散信号序列;
步骤四,在单片机的AD中断中部署交流采样及信号处理、滤波算法,当采集的数据量累积到一定数量时,程序前端会根据累积的缓存数据进行周期任务或秒任务计算,计算电压和电流的有效值、电网频率、功率参数及相关电气参数;在测量前完成对交采得到的每相的电压、电流进行残余直流滤波,以消除直流分量对测量的影响,实现AD中断计量;AD中断计量处理过程为,读取三相电压值和三相电流值,对每一相进行残余直流滤波,电压值、电流值数据量累积,周期任务或秒任务计算。采集的数据量累积到一定数量,该数量一般为1个工频周期或1秒对应的数据量。
本实施例中:利用MATLAB仿真软件对直流滤波算法进行仿真验证。构造正弦波形,设置直流分量为1。如图2为含有直流分量的滤波之前波形,其中横坐标为时间,纵坐标为幅值。本实施例分别针对K=1/214、K=1/212以及K=1/28三种情况,仿真利用该滤波算法进行直流滤波时直流分量收敛的波形以及滤波后波形,结果分别如图3为本发明K=1/214时直流滤波直流分量收敛波形;图4为本发明K=1/214时直流滤波后波形;图5为本发明K=1/212时直流滤波直流分量收敛波形;图6为本发明K=1/212时直流滤波后波形;图7为本发明K=1/28时直流滤波直流分量收敛波形;图8为本发明K=1/28时直流滤波后波形。
比较图3、图5、图7可以看出,随着K值的变化,滤波效果有明显的不同。结合K值越小,直流分量收敛越慢,但收敛后波形波动小,稳定性较好。因此,K=1/214时直流分量收敛最慢,在8s的时候还没有完全收敛到1,而当K=1/212时,直流分量在3s-4s之间就已经收敛到1,随着K的进一步增大,当K=1/28时,直流分量收敛的更快,但相对来说其波动性也更大,收敛稳定性不好。
比较图4、图6、图8滤波后波形可以看出,利用该算法进行直流滤波均能达到滤除直流分量的目的,区别在于随着K值的不同,滤波后的波形稳定在零值附近的时间不同。K值越大,稳定时间越快,但稳定后的误差也越大,如图8所示,稳定后的波形幅值上下限并没有达到1和-1。
分析图3至图8可以看到,当K=1/212时,直流分量具有较快的收敛速度,同时滤波后的波形稳定后幅值没有衰减,具有较好的滤波效果。因此,在使用该算法进行直流滤波时,根据采集数据的频率以及采样参数选择合适的K值,能够保证滤波后的波形在具有较好稳定性的同时,同时具有较快的收敛速度。
综上所述,本发明通过FTU内置单片机的滤波算法根据采集到的电压、电流瞬时值实现自适应的滤除残余直流,消除残余直流对交采精度的影响,使交采得到的电压、电流值更加准确,确保配电网参数的测量有更高的精度,同时减少硬件设备的投资、省略实际的滤波电路,能够实现馈线终端的软件计量,充分发挥计算机算法具有的计算高精度的优点。
以上实施例是对本发明的具体实施方式的说明,而非对本发明的限制,有关技术领域的技术人员在不脱离本发明的精神和范围的情况下,还可以做出各种变换和变化而得到相对应的等同的技术方案,因此所有等同的技术方案均应该归入本发明的专利保护范围。

Claims (4)

1.一种可滤除交流采样中残余直流分量的馈线终端计量方法,其特征在于:包括以下步骤,
步骤一,建立直流滤波器差分方程;
步骤二,根据滤波器原理及差分方程推导出滤波算法;
步骤三,采用单片机进行数字化处理,根据采样定理对所需要采集的交流模拟量进行等间隔采样,从而获得一个离散信号序列;
步骤四,在单片机的AD中断中部署交流采样及信号处理、滤波算法,当采集的数据量累积到一定数量时,程序前端会根据累积的缓存数据进行周期任务或秒任务计算,计算电压和电流的有效值、电网频率、功率参数及相关电气参数;在测量前完成对交采得到的每相的电压、电流进行残余直流滤波,以消除直流分量对测量的影响,实现AD中断计量。
2.根据权利要求1所述的一种可滤除交流采样中残余直流分量的馈线终端计量方法,其特征在于:步骤一中建立直流滤波器差分方程为,
y[n]=(x[n]-y[n-1])×K+y[n-1],K=1/2m
其中,y[n]为当前采样点对应的直流分量输出,y[n-1]为上一采样点对应的直流分量输出,x[n]为当前瞬时采样值输入,K为滤波器系数,m值则根据实际情况进行选择;m值越大,K值越小,直流分量收敛越慢,但是波动较小,稳定性较高;m值越小,K值越大,直流分量收敛的越快,但是波动比较大,稳定性较差。
3.根据权利要求1所述的一种可滤除交流采样中残余直流分量的馈线终端计量方法,其特征在于:步骤二中根据滤波器原理及差分方程推导出滤波算法为,
yn=((xn-yn-1)>>m+yn-1
xn=xn-yn
其中,m为滤波器系数中的指数值,根据实际应用选择合适的m值;yn和yn-1分别代表当前和上一次计算得到的直流分量,xn为本次采集瞬时值。
4.根据权利要求1所述的一种可滤除交流采样中残余直流分量的馈线终端计量方法,其特征在于:步骤四中AD中断计量处理过程为,读取三相电压值和三相电流值,对每一相进行残余直流滤波,电压值、电流值数据量累积,周期任务或秒任务计算;采集的数据量累积到一定数量,该数量一般为1个工频周期或1秒对应的数据量。
CN201810597330.1A 2018-06-06 2018-06-06 一种可滤除交流采样中残余直流分量的馈线终端计量方法 Withdrawn CN108957111A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810597330.1A CN108957111A (zh) 2018-06-06 2018-06-06 一种可滤除交流采样中残余直流分量的馈线终端计量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810597330.1A CN108957111A (zh) 2018-06-06 2018-06-06 一种可滤除交流采样中残余直流分量的馈线终端计量方法

Publications (1)

Publication Number Publication Date
CN108957111A true CN108957111A (zh) 2018-12-07

Family

ID=64488156

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810597330.1A Withdrawn CN108957111A (zh) 2018-06-06 2018-06-06 一种可滤除交流采样中残余直流分量的馈线终端计量方法

Country Status (1)

Country Link
CN (1) CN108957111A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111061422A (zh) * 2019-12-13 2020-04-24 明见(厦门)技术有限公司 一种逼近式直流分量提取方法、终端设备及存储介质
WO2020125494A1 (zh) * 2018-12-17 2020-06-25 中国科学院深圳先进技术研究院 一种生物信号的分析处理方法、装置及设备
CN112505612A (zh) * 2020-10-26 2021-03-16 青岛鼎信通讯股份有限公司 一种应用于电能表检定装置的标准表电能脉冲输出方法
CN112817045A (zh) * 2021-01-28 2021-05-18 中国科学院地质与地球物理研究所 一种时域脱敏的直流偏置实时去除装置及方法
CN114076678A (zh) * 2020-08-19 2022-02-22 计算系统有限公司 从振动波形中移除渐近衰减的直流偏置影响

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101505060A (zh) * 2009-03-18 2009-08-12 天津大学 电网中利用窄带滤波实现多采样率信号的频率转换方法
CN102005744A (zh) * 2010-12-15 2011-04-06 国网电力科学研究院 基于差流过零识别ct电流缓慢衰减形成差流的方法
CN103926449A (zh) * 2014-04-24 2014-07-16 江苏省电力公司常州供电公司 电力电缆接地电流的自适应监测方法
CN105738686A (zh) * 2016-04-12 2016-07-06 辽宁大学 一种电流检测方法
CN106372369A (zh) * 2016-10-08 2017-02-01 国家电网公司 用于逆变器输出直流分量、谐波分量检测的瞬时计算方法
CN107102193A (zh) * 2016-02-23 2017-08-29 上海贝岭股份有限公司 电参量计量过程中的信号噪声处理装置
CN107656134A (zh) * 2017-10-13 2018-02-02 国网安徽省电力公司经济技术研究院 一种自适应滤除直流分量的故障电流检测方法及装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101505060A (zh) * 2009-03-18 2009-08-12 天津大学 电网中利用窄带滤波实现多采样率信号的频率转换方法
CN102005744A (zh) * 2010-12-15 2011-04-06 国网电力科学研究院 基于差流过零识别ct电流缓慢衰减形成差流的方法
CN103926449A (zh) * 2014-04-24 2014-07-16 江苏省电力公司常州供电公司 电力电缆接地电流的自适应监测方法
CN107102193A (zh) * 2016-02-23 2017-08-29 上海贝岭股份有限公司 电参量计量过程中的信号噪声处理装置
CN105738686A (zh) * 2016-04-12 2016-07-06 辽宁大学 一种电流检测方法
CN106372369A (zh) * 2016-10-08 2017-02-01 国家电网公司 用于逆变器输出直流分量、谐波分量检测的瞬时计算方法
CN107656134A (zh) * 2017-10-13 2018-02-02 国网安徽省电力公司经济技术研究院 一种自适应滤除直流分量的故障电流检测方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
寇英刚 等: "馈线终端交流采样中的残余直流滤波算法研究", 《电测与仪表》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020125494A1 (zh) * 2018-12-17 2020-06-25 中国科学院深圳先进技术研究院 一种生物信号的分析处理方法、装置及设备
CN111061422A (zh) * 2019-12-13 2020-04-24 明见(厦门)技术有限公司 一种逼近式直流分量提取方法、终端设备及存储介质
CN111061422B (zh) * 2019-12-13 2022-09-02 明见(厦门)技术有限公司 一种逼近式直流分量提取方法、终端设备及存储介质
CN114076678A (zh) * 2020-08-19 2022-02-22 计算系统有限公司 从振动波形中移除渐近衰减的直流偏置影响
CN114076678B (zh) * 2020-08-19 2024-04-26 计算系统有限公司 从振动波形中移除渐近衰减的直流偏置影响
CN112505612A (zh) * 2020-10-26 2021-03-16 青岛鼎信通讯股份有限公司 一种应用于电能表检定装置的标准表电能脉冲输出方法
CN112817045A (zh) * 2021-01-28 2021-05-18 中国科学院地质与地球物理研究所 一种时域脱敏的直流偏置实时去除装置及方法
CN112817045B (zh) * 2021-01-28 2021-08-27 中国科学院地质与地球物理研究所 一种时域脱敏的直流偏置实时去除装置及方法

Similar Documents

Publication Publication Date Title
CN108957111A (zh) 一种可滤除交流采样中残余直流分量的馈线终端计量方法
CN105044456B (zh) 一种基于正交子带的电网瞬时频率测量与跟踪方法
CN101441231B (zh) 一种谐波电能计量误差分析装置
Bayindir et al. Development of a real time energy monitoring platform
CN104375007B (zh) 一种基于相位识别的土壤电阻率测量方法及装置
CN105606900B (zh) 一种基于方波信号的单相谐波阻抗测量方法
CN110470926A (zh) 一种交流电网背景谐波阻抗的实用估算方法
CN105911495A (zh) 一种带有iec61850数字接口的模拟标准源及其检测方法
CN102253281A (zh) 一种基于软件抗谐波干扰并自动换相测量电网频率的方法
CN108333426A (zh) 基于傅氏算法的电力系统频率测量方法
CN102981083A (zh) 自校准电能质量监测装置
CN109541450A (zh) 一种基于matlab的直流开关振荡回路参数测量系统及计算方法
CN105021869A (zh) 一种低压配电短路电流峰值预测方法
CN110231581A (zh) 蓄电池内阻测量电路及测量方法
CN106324353A (zh) 一种感应滤波换流变压器谐波阻抗测量方法及装置
CN102522742B (zh) 基于单点量测信息的外网戴维南等值参数的估计方法
CN107238750B (zh) 一种基于运放电路的电网波形侦测方法
CN109581045A (zh) 一种满足iec标准框架的间谐波功率计量方法
CN110873822A (zh) 一种馈线终端交流采样中频率计算的方法
CN105866576A (zh) 智能变电站二次侧电能计量误差影响的模拟检测系统及其检测分析方法
CN105375470A (zh) 一种利用磁场数据逆推架空输电线路三相电流的方法
CN112731177A (zh) 干电池电量检测方法及装置
CN106526365A (zh) 单火线用电数据测量方法及单火线智能断路器
CN202305668U (zh) 基于can总线的新能源发电电能质量监测装置
CN100508388C (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
DD01 Delivery of document by public notice

Addressee: Sun Siqi

Document name: Notice of first office action

DD01 Delivery of document by public notice
WW01 Invention patent application withdrawn after publication

Application publication date: 20181207

WW01 Invention patent application withdrawn after publication