CN103257114B - 基于光纤Bragg光栅传感网络的井下瓦斯检测方法 - Google Patents

基于光纤Bragg光栅传感网络的井下瓦斯检测方法 Download PDF

Info

Publication number
CN103257114B
CN103257114B CN201310167688.8A CN201310167688A CN103257114B CN 103257114 B CN103257114 B CN 103257114B CN 201310167688 A CN201310167688 A CN 201310167688A CN 103257114 B CN103257114 B CN 103257114B
Authority
CN
China
Prior art keywords
bragg grating
optical fiber
microprocessor
grating sensor
fiber bragg
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.)
Expired - Fee Related
Application number
CN201310167688.8A
Other languages
English (en)
Other versions
CN103257114A (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.)
Anhui University of Science and Technology
Original Assignee
Anhui 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 Anhui University of Science and Technology filed Critical Anhui University of Science and Technology
Priority to CN201310167688.8A priority Critical patent/CN103257114B/zh
Publication of CN103257114A publication Critical patent/CN103257114A/zh
Application granted granted Critical
Publication of CN103257114B publication Critical patent/CN103257114B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种基于光纤Bragg光栅传感网络的井下瓦斯检测方法,包括有光纤检测模块、数据采集中心、监控中心,光纤检测模块包括有光源、光纤Bragg光栅传感器,数据采集中心包括有依次连接的光电探测器、信号放大器、AD转换器,监控中心包括有依次连接的微处理器、计算机,光纤Bragg光栅传感器连接光电探测器,AD转换器连接微处理器;通过监测点的光信号被光纤Bragg光栅传感器的Bragg光栅滤波分成两束后被光电探测器分别接收,进行放大滤波后的数字信号量采用去极值数字滤波法进行数据处理转换后由微处理器分析处理并计算出瓦斯浓度。本发明很大程度上提升了系统的准确性和稳定性,确保监控系统的连续性、可靠性及有效性,最终实现“监控预防”的目的。

Description

基于光纤Bragg光栅传感网络的井下瓦斯检测方法
技术领域
本发明涉及到井下瓦斯监控系统的传感技术领域,具体为一种基于光纤Bragg光栅传感网络的井下瓦斯检测方法。
背景技术
信息技术数据的采集主要依靠传感器和探测器进行,因此新型的传感技术在数据采集中尤其是在恶劣的井下环境中采集而来的数据在井下瓦斯监控系统中就显的尤为突出,作为煤矿安全监控系统中一项重要的指标,传感器网络自井下采集来的数据直接关系到监测系统的监控效率,所以传感器的稳定性、可靠性、以及传感器采集来的数据在瓦斯监控整个系统中就起到至关重要的作用。目前国内开发的的传感器虽然基本满足了井下瓦斯抽采系统的需求,但是面对恶劣的井下工作环境,现有的传感器无论是在敏感元件的质量上还是在传感器测量数据的准确性上都有待提高,长期工作稳定性差,调校频繁、复杂以及工作寿命短已经成为现有传感器无法回避的问题,所以本文提出一种光纤Bragg传感器代替传统的瓦斯传感器以应对现有系统对传感器使用寿命、稳定性、调节周期,精确度等提出的更高要求的现状。
为了保证煤矿安全监控系统的有效性和稳定性,研制开发一种基于光纤Bragg光栅传感网络的新型井下瓦斯监控系统很有必要,这样不仅可以保证井下数据的精确性,监测系统的稳定连续性,煤矿安全监控系统的正常运行,而且可以最大程度上避免因事故而带来的损失。
发明内容
本发明要解决的技术问题是针对现有矿井安全监控系统中传感系统的不足,提供一种采用差分吸收检测原理的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,以解决现有技术检测方法精度低,易受干扰,实时性差等的问题。
为了达到上述目的,本发明所采用的技术方案为:
基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:包括有瓦斯监测系统,所述的瓦斯监测系统包括有光纤检测模块、数据采集中心、监控中心,所述的光纤检测模块包括有光源、光纤Bragg光栅传感器,数据采集中心包括有依次连接的光电探测器、信号放大器、AD转换器,监控中心包括有依次连接的微处理器、计算机,所述的光纤Bragg光栅传感器的信号输出端连接光电探测器,AD转换器的信号输出端连接微处理器,微处理器还连接有数据存储器、声光报警器,计算机还连接有LCD显示模块;通过监测点的光信号被光纤Bragg光栅传感器的Bragg光栅滤波分成两束后被光电探测器分别接收,进行放大滤波后的数字信号量采用去极值数字滤波法进行数据处理转换后由微处理器分析处理并计算出瓦斯浓度,具体包括以下步骤:
(1)根据待测区域不同的地域情况,采用多个光纤Bragg光栅传感器分别检测不同待检测区域的待测参数,待测参数以光信号为载体返回控制室;
(2)所述监测点的光信号返回控制室被光纤Bragg光栅传感器的Bragg光栅滤波分成两束被光电探测器接收,电信号经过信号放大器处理后再经AD转换器转换为数字信号送入微处理器进行处理;
(3)在步骤(2)中所述微处理器内部操作系统中创建多个系统任务,用来保证系统各个模块功能的连续进行,所述的微处理器包含1个16通道差分输入多路选择器和1个10位SAR ADC并对接收来自A/D采样转换后的数字量采用去极值数字滤波法做了数据处理后,采用差分吸收法检测计算分析出结果,所述待检测区域的瓦斯浓度通过如下公式计算出来:
有比尔-朗伯定律:
I(λ)=K(λ)I0(λ)exp[-α(λ)LC+β(λ)+γ(λ)L+δ(λ)]
在此公式中,I(λ)为透射光的光强;I0(λ)为入射光的光强;α(λ)为待测气体的吸收系数;K(λ)为整个系统的光效率;β(λ)为瑞利散射系数;γ(λ)为米氏散射系数;C为待测气体浓度;L为待测气体的整个吸收光程;δ(λ)为气体密度造成波动的吸收系数;
波长相差较小分别为λ1、λ2且在吸收系数上有差别,对应待测气体的参考和工作波长,得到瓦斯气体浓度C:
C = 1 [ α ( λ 1 ) - α ( λ 2 ) L ln I ( λ 2 ) I ( λ 1 ) .
所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的光源采用InGaAsP/InP超辐射的发光二极管所提供的红外光波。
所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的光纤Bragg光栅传感器通过选定的用户分配各自唯一的地址码,同时进行扩频处理,数据接收时可以利用相关运算恢复初始指定信号并对井下多接触面多点进行分布式的定点测量。
所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的微处理器采用C8051F330。
本发明的有益效果在于:
本发明可以在井下复杂多变的环境下仍然可以在最大程度上保证监测系统的稳定性和监测数据的准确性,在面对各种干扰的情况下,基于光纤Bragg光栅传感网络瓦斯检测方法同时把编码技术和传感技术相结合并融合于一体,基本实现了在待测地域监测点的多点分布式测量,本系统同时还对环境中的例如温度、暗电流的造成的噪声等变化均有较强的抗干扰能力,很大程度上提升了系统的准确性和稳定性,确保监控系统的的连续性、可靠性及有效性,最终实现“监控预防”的目的。
附图说明
图1为本发明的系统图。
图2为本发明微处理器C8051F330建立的主程序流程图。
具体实施方式
如图1所示。基于光纤Bragg光栅传感网络的井下瓦斯检测方法,包括有瓦斯监测系统,瓦斯监测系统包括有光纤检测模块1、数据采集中心2、监控中心3,光纤检测模块1包括有光源4、光纤Bragg光栅传感器5,数据采集中心2包括有依次连接的光电探测器6、信号放大器7、AD转换器8,监控中心3包括有依次连接的微处理器9、计算机10,光纤Bragg光栅传感器5的信号输出端连接光电探测器6,AD转换器8的信号输出端连接微处理器9,微处理器9还连接有数据存储器11、声光报警器12,计算机10还连接有LCD显示模块13;通过监测点的光信号被光纤Bragg光栅传感器5的Bragg光栅滤波分成两束后被光电探测器分别接收,进行放大滤波后的数字信号量采用去极值数字滤波法进行数据处理转换后由微处理器分析处理并计算出瓦斯浓度,具体包括以下步骤:
(1)根据待测区域不同的地域情况,采用多个光纤Bragg光栅传感器5分别检测不同待检测区域的待测参数,待测参数以光信号为载体返回控制室;
(2)所述监测点的光信号返回控制室被光纤Bragg光栅传感器5的Bragg光栅滤波分成两束被光电探测器6接收,电信号经过信号放大器7处理后再经AD转换器8转换为数字信号送入微处理器9进行处理;
(3)在步骤(2)中所述微处理器9内部操作系统中创建多个系统任务,用来保证系统各个模块功能的连续进行,所述的微处理器9包含1个16通道差分输入多路选择器和1个10位SAR ADC并对接收来自A/D采样转换后的数字量采用去极值数字滤波法做了数据处理后,采用差分吸收法检测计算分析出结果,所述待检测区域的瓦斯浓度通过如下公式计算出来:
有比尔-朗伯定律:
I(λ)=K(λ)I0(λ)exp[-α(λ)LC+β(λ)+γ(λ)L+δ(λ)]
在此公式中,I(λ)为透射光的光强;I0(λ)为入射光的光强;α(λ)为待测气体的吸收系数;K(λ)为整个系统的光效率;β(λ)为瑞利散射系数;γ(λ)为米氏散射系数;C为待测气体浓度;L为待测气体的整个吸收光程;δ(λ)为气体密度造成波动的吸收系数;
波长相差较小分别为λ1、λ2且在吸收系数上有差别,对应待测气体的参考和工作波长,得到瓦斯气体浓度C:
C = 1 [ α ( λ 1 ) - α ( λ 2 ) L ln I ( λ 2 ) I ( λ 1 ) .
本发明采用美国Cygnal公司的C8051F330单片机作为系统的控制核心,一方面根据不同的待测区域不同的地域情况,采用多个光纤Bragg光栅传感器分别用来检测待检测区域的待测参数,来自InGaAsP/InP超辐射的发光二极管的红外光波经过待测地域后被接收转换为电信号,对电信号进行处理后给微处理器。另一方面微处理器先采用去极值数字滤波方法对采样转换来的数据进行处理,同时微处理器通过用户分配唯一的地址码和对的Bragg传感器进行扩频处理,以保证系统资源的合理利用以及增多传感器个数以加大区域监测范围,还可以确保当数据接收中断时可以利用相关运算恢复初始信号,监控主机采用以太网完成和数据库服务器数据之间的通讯;远程客户终端的WEB浏览服务由本地数据库的服务器提供,以便随时随地了解当前系统监控信息。
1. 差分吸收检测原理
当一束光通过某种介质时,就会有光的一部分会被这种介质吸收,少量的光波通过介质会被介质散射,其余的光波则沿着原来的方向继续的传播,所以有比尔-朗伯定律可知,介质的吸收长度、透射光的光强、光的吸收系数与物质浓度之间满足:
I=I0exp[-αLC+β+γL+δ]
式中,I,I0分别为透射和入射光强;L为待测气体的吸收光程、;C为待测气体的浓度;α为某单位波长下的浓度、单位长度介质的吸收系数;β为瑞利散射系数;δ为气体密度波动造成的吸收系数;γ为米氏散射系数。用光来检测气体时,应该还要包括比例系数K,为此公式可改写为
I=KI0exp[-αLC+β+γL+δ]
但是仅利用上式确定待测气体的浓度C时还是不足的,因为δ的反应的为平均变化数,随着时间的变化而不停的发生变化,所以它是随机量。所以用两个吸收系数有较大差距但是波长相近的两束波长为λ1、λ2的单色光在极短的时间通过将要测量的气体,有:
I(λ1)=K(λ1)I01)exp[-α(λ1)LC+β(λ1)+γ(λ1)L+δ(λ1)]
I(λ2)=K(λ2)I02)exp[-α(λ2)LC+β(λ2)+γ(λ2)L+δ(λ2)]
气体浓度关系表达式为:
C = 1 [ α ( λ 1 ) - α ( λ 2 ) L [ ln K ( λ 1 ) I ( λ 2 ) K ( λ 2 ) I ( λ 1 ) - ln I ( λ 1 ) I ( λ 2 ) - [ β ( λ 1 ) - β ( λ 2 ) ] L - [ γ ( λ 1 ) - γ ( λ 2 ) ] L - [ δ ( λ 1 ) - δ ( λ 2 ) ] ]
由于λ1、λ2数值相差较小,同时单色光也几乎是同一时间接近然后通过待测 气体的,则:
β(λ1)≈β(λ2),γ(λ1)≈γ(λ2),δ(λ1)≈δ(λ2),
所以 C = 1 [ α ( λ 1 ) - α ( λ 2 ) ] L ln K ( λ 1 ) I 0 ( λ 1 ) I ( λ 2 ) K ( λ 2 ) I 0 ( λ 2 ) I ( λ 1 )
对光学系统内部进行调节使的:
K(λ1)I01)=K(λ2)I02)
可得 C = 1 [ α ( λ 1 ) - α ( λ 2 ) ] L ln I ( λ 2 ) I ( λ 1 )
有光学上的相近条件可知:
[α(λ1)-α(λ2)]L≤1
在实际监控过程中,波长为λ1、λ2分别对应着不同吸收谱线,波长为λ2的光作为参考波长不能被待测气体吸收,检测过程中为空值,因此I(λ1)<I(λ2),I(λ2)/I(λ1)>1和I(λ1)/I(λ2)<1。对ln[I(λ2)/I(λ1)]进行泰勒展开:
ln I ( &lambda; 2 ) I ( &lambda; 1 ) = - ln [ + 1 ( I ( &lambda; 1 ) I ( &lambda; 2 ) - 1 ) ] &ap; I ( &lambda; 2 ) - ( &lambda; 1 ) I ( &lambda; 2 )
所以气体的浓度为
C = I ( &lambda; 2 ) - ( &lambda; 1 ) I ( &lambda; 2 ) 1 &alpha; ( &lambda; 1 ) - &alpha; ( &lambda; 2 )
在波长λ1、λ2、α1、α2已知情况下就可以通过I(λ2)-I(λ1)和I(λ2)关系式来测量待测气体浓度。
1)去极值的数字滤波方法
系统软件主要用来控制数据的分析与处理以及数据的保存,分析计算出的结果的显示以及需要保存的数据经由RS-232接口传递给上位机等功能更加强大的模块进行。启动脉冲来临前,首先进行初始化处理,对A/D、定时存储器等完成前期的初始化工作然后等待中断请求的到来,为了提高检测精度和消除随机噪声干扰,采用去极值的数字滤波方法对A/D采样转换来的数据进行处理。将事先标定好的线性系数存储在系统的ROM存储器中,当需要做数据处理时,根据当前分析所得数据获得相关的系数,求取当前瓦斯气体的浓度数据送显示器显示,当待测企图浓度超过设定的限定值时,处理器发出声光报警。软件部分采用模块化设计,程序的主流程图如图2所示
为了达到提高待测气体浓度数据精度的目的,对经过A/D采样转换后的数字量采用去极值数字滤波法进行数据处理。此方法首先需要建立一个具有一定容量的窗口平台,窗口内包含的数据量为n。将这n个数据放到一个数组当中送到微处理器进行循环分析处理得到n个结果,将n个结果进行比较获得该组检测数据当中待测值最小值与最大值,之后再计算出所有数据结果的累加和,在这个累加和的数值当中减去循环得到的最大值和最小值,这样的累加和的数值就组成变成了由原来的n个数据变成了n-2个数据,把得到n-2个数据进行累加再除以n-2得到的均值系统就当做一个监测数据。在一组数据的累加和中去除去最小值和最大值可以有效的消除一组监测数据中可能出现的变化较大的随机信号变量对最终测量结果的影响,从而达到了消除脉冲干扰,提高精确率的良好效果。根据误差理论的原则,相同的物理量在等同等精度测量初始条件下条件下,当测量次数越大时,其结果所得到的均值就越逼近真实值,当测量次数>10次之后,测量精度的提升将会放缓,这样就可以在确保检测精度和满足实时检测需要的这两个标准的要求下,可以采用选择连续采样12次的检测数据利用去极值数字滤波法处理后得到的数据做为当前瓦斯浓度检测的最终结果。

Claims (3)

1.基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:包括有瓦斯监测系统,所述的瓦斯监测系统包括有光纤检测模块、数据采集中心、监控中心,所述的光纤检测模块包括有光源、光纤Bragg光栅传感器,数据采集中心包括有依次连接的光电探测器、信号放大器、AD转换器,监控中心包括有依次连接的微处理器、计算机,所述的光纤Bragg光栅传感器的信号输出端连接光电探测器,AD转换器的信号输出端连接微处理器,微处理器还连接有数据存储器、声光报警器,计算机还连接有LCD显示模块;通过监测点的光信号被光纤Bragg光栅传感器的Bragg光栅滤波分成两束后被光电探测器分别接收,进行放大滤波后的数字信号量采用去极值数字滤波法进行数据处理转换后由微处理器分析处理并计算出瓦斯浓度,所述的光纤Bragg光栅传感器通过选定的用户分配各自唯一的地址码,同时进行扩频处理,数据接收时利用相关运算恢复初始指定信号并对井下多接触面多点进行分布式的定点测量;具体包括以下步骤:
(1)根据待测区域不同的地域情况,采用多个光纤Bragg光栅传感器分别检测不同待检测区域的待测参数,待测参数以光信号为载体返回控制室;
(2)所述监测点的光信号返回控制室被光纤Bragg光栅传感器的Bragg光栅滤波分成两束被光电探测器接收,电信号经过信号放大器处理后再经AD转换器转换为数字信号送入微处理器进行处理;
(3)在步骤(2)中所述微处理器内部操作系统中创建多个系统任务,用来保证系统各个模块功能的连续进行,所述的微处理器包含1个16通道差分输入多路选择器和1个10位SAR ADC并对接收来自A/D采样转换后的数字量采用去极值数字滤波法做了数据处理后,采用差分吸收法检测计算分析出结果,所述待检测区域的瓦斯浓度通过如下公式计算出来:
比尔-朗伯定律:
在此公式中,为透射光的光强;为入射光的光强;为待测气体的吸收系数;为整个系统的光效率;为瑞利散射系数;为米氏散射系数;C为待测气体浓度;L为待测气体的整个吸收光程;为气体密度造成波动的吸收系数;
波长相差较小分别为且在吸收系数上有差别,对应待测气体的参考和工作波长,得到瓦斯气体浓度C:
2.根据权利要求1所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的光源采用InGaAsP/InP超辐射的发光二极管所提供的红外光波。
3.根据权利要求1所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的微处理器采用C8051F330。
CN201310167688.8A 2013-05-08 2013-05-08 基于光纤Bragg光栅传感网络的井下瓦斯检测方法 Expired - Fee Related CN103257114B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310167688.8A CN103257114B (zh) 2013-05-08 2013-05-08 基于光纤Bragg光栅传感网络的井下瓦斯检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310167688.8A CN103257114B (zh) 2013-05-08 2013-05-08 基于光纤Bragg光栅传感网络的井下瓦斯检测方法

Publications (2)

Publication Number Publication Date
CN103257114A CN103257114A (zh) 2013-08-21
CN103257114B true CN103257114B (zh) 2015-04-22

Family

ID=48961161

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310167688.8A Expired - Fee Related CN103257114B (zh) 2013-05-08 2013-05-08 基于光纤Bragg光栅传感网络的井下瓦斯检测方法

Country Status (1)

Country Link
CN (1) CN103257114B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103926209B (zh) * 2013-12-26 2017-01-18 安徽理工大学 一种基于doas双采样技术的煤层瓦斯联测系统
CN111380670A (zh) * 2020-02-28 2020-07-07 浙江工业大学 一种体状及微棱镜镜片的棱镜度智能检测系统
CN111397852A (zh) * 2020-02-28 2020-07-10 浙江工业大学 一种体状及微棱镜镜片的棱镜度智能检测方法
CN114200875A (zh) * 2020-07-02 2022-03-18 广东电网有限责任公司江门供电局 一种超辐射双光源外破隐患巡查系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101030822B (zh) * 2006-02-27 2012-02-01 华为技术有限公司 光网络单元和光线路终端及无源光纤接入网及其传输方法
CN102207458A (zh) * 2011-04-08 2011-10-05 东北大学 基于光纤光栅的近红外吸收浓度传感器信号解调方法

Also Published As

Publication number Publication date
CN103257114A (zh) 2013-08-21

Similar Documents

Publication Publication Date Title
CN103196465B (zh) 一种相敏光时域反射仪传感信号噪声分离及信号提取方法
CN104864979B (zh) 一种分布式拉曼光纤测温系统测量误差的修正方法
CN102879081B (zh) 一种分布式光纤振动系统中的数据处理方法
CN103439232A (zh) 一种减光式烟尘颗粒浓度测量方法及其装置
CN104964699B (zh) 基于φ‑OTDR光纤分布式扰动传感器的扰动判断方法和装置
CN103257114B (zh) 基于光纤Bragg光栅传感网络的井下瓦斯检测方法
CN102292621A (zh) 分布式光纤感测中的改进
CN102636459A (zh) 一种前向散射与透射结合的能见度测量仪及其测量方法
CN204556499U (zh) 调谐二极管吸收光谱的多通道高速数据采集和处理系统
CN103727968A (zh) 一种分布式光纤传感装置及同时测量温度、应变、振动的方法
CN103439233A (zh) 烟道粉尘浓度检测系统
CN105953941A (zh) 一种基于拉曼散射的分布式光纤测温方法及装置
CN105157874A (zh) 一种分布式边界安防监测系统及方法
CN107545684B (zh) 一种分布式光纤振动传感系统中的多点扰动事件报警方法
CN107727122B (zh) 双端探测的联合拉曼和布里渊散射的分布式光纤传感装置
CN108344432B (zh) 获取布里渊光纤分布式传感系统的传感信息的方法
CN109444895B (zh) 分布式振动传感器消除干涉衰落的振动信息定位方法
CN202631153U (zh) 带有自动补偿功能的单端口分布式光纤温度传感器
CN115615332A (zh) 一种光纤复合架空地线的状态在线监测系统
CN115752696A (zh) 基于前向传输的分布式光纤振动传感系统及信号处理方法
CN205068674U (zh) 一种智能型分布式光纤感温油田油井火灾预警系统
CN211696444U (zh) 一种可实现长距离输电线路融冰监测的光纤监测系统
US20180340831A1 (en) Wavelength determination using an optical filter having complementary transmission and reflection coefficients
CN204879459U (zh) 一种光纤法集油管线泄漏检测装置
CN207515900U (zh) 基于脉冲编码和边沿滤波法的光纤光栅振动测量系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150422

Termination date: 20160508