CN103644460A - 一种输气管线泄漏音波信号的滤波优选方法 - Google Patents
一种输气管线泄漏音波信号的滤波优选方法 Download PDFInfo
- Publication number
- CN103644460A CN103644460A CN201310468802.0A CN201310468802A CN103644460A CN 103644460 A CN103644460 A CN 103644460A CN 201310468802 A CN201310468802 A CN 201310468802A CN 103644460 A CN103644460 A CN 103644460A
- Authority
- CN
- China
- Prior art keywords
- sound wave
- wave signals
- filtering
- transmission line
- leaks
- 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.)
- Pending
Links
Images
Landscapes
- Examining Or Testing Airtightness (AREA)
Abstract
本发明公开了一种输气管线泄漏音波信号的滤波优选方法,该方法包括:泄漏音波信号下降沿长度为信号幅值突降至幅值恢复至零值的信号长度;拟合真实泄漏音波信号时域波形,得到拟合泄漏音波信号;添加噪声信号,获得伪泄漏音波信号;利用预选择的各种滤波方法对伪泄漏音波信号进行滤波处理,得到滤波后波形;计算并获得滤波后信号的评价指标,确定最优滤波方法,以真实泄漏音波信号验证最优滤波方法的有效性。本发明通过满足信噪比评价指标最大、均方根误差评价指标最小并且泄漏下降沿幅值损失评价指标最小,选择最优滤波方法,定量地评价了泄漏音波信号的滤波效果,增大了泄漏音波信号的信噪比,提高了音波法泄漏检测技术的适用性。
Description
技术领域
本发明属于信号处理与管道安全监控技术领域,尤其涉及一种输气管线泄漏音波信号的滤波优选方法。
背景技术
音波法泄漏检测技术使用的传感器为动态压力传感器,传感器以管内流体的压力脉动作为输出值,管道发生泄漏时,管内压力骤降,传感器输出波形形成一个泄漏下降沿,即泄漏发生时刻,信号压力幅值出现一个脉冲下降,随后又迅速恢复至零值附近波动,泄漏下降沿是进行泄漏判断的依据,但其持续时间短,波形变化剧烈,同时,泄漏音波信号在传播过程中存在衰减,加之背景噪声叠加干扰,使得泄漏音波信号的滤波过程变得非常困难。
管道管径、结构和运行环境的不同,导致不同输气管道的泄漏音波信号和背景噪声信号各不相同,一种滤波方法对不同工况下采集到的信号的滤波效果也不相同。
音波法泄漏检测技术的原理以及众多的干扰因素对信号滤波方法提出了较高的要求,单纯寻求一种固定的滤波方法,对不同运行工况的输气管线适用性有限,一种有效的解决方案是综合分析比较各种 滤波方法,某一实际工况,选出一种适应性最强、滤波效果最佳的方法。
根据调研结果,现阶段国内外涉及输气管道泄漏音波信号处理的专利主要有:
美国专利US6389881公开了一种利用音波技术进行管道泄漏检测的技术,该技术利用传感器采集管内动态压力,采用模式匹配滤波技术对信号进行滤波处理,排除噪声,降低干扰,提高了定位精度;
中国实用新型专利200820078616.0公开了一种音波泄漏检测装置,传感器采集到的音波信号经前置放大后进行滤波处理,滤波为带通滤波,通过分路数据处理器实现滤波过程;
中国专利200710177617.0公开了一种基于压力信号和音波信号的泄漏检测定位技术,利用神经网络实现音波信号的滤波、特征融合和决策融合处理,采用的滤波方法为中值滤波;
现有的专利较少涉及泄漏音波信号的滤波过程,只是简单提及几种滤波方法,对滤波方法的选择过程没有描述,具体表现为:
(1)对泄漏音波信号的滤波方法缺乏统一的认知,现有专利没有给出采用某一种滤波方法的具体依据;
(2)对泄漏音波信号的滤波只使用一种固定的滤波方法,不能保证滤波效果最优,不能确保滤波方法对不同运行工况下的输气管道均具有适用性。
发明内容
本发明的目的在于提供一种输气管线泄漏音波信号的滤波优选方法,旨在解决现阶段缺乏选择泄漏音波信号滤波方法的问题。
本发明是这样实现的,一种输气管线泄漏音波信号的滤波优选方法,该输气管线泄漏音波信号的滤波优选方法包括以下步骤:
步骤一,泄漏音波信号下降沿长度为信号幅值突降至幅值恢复至零值的信号长度;拟合真实泄漏音波信号时域波形,得到拟合泄漏音波信号;
步骤二,添加噪声信号,获得伪泄漏音波信号;
步骤三,利用预选择的各种滤波方法对伪泄漏音波信号进行滤波处理,得到滤波后波形;
步骤四,计算并获得滤波后信号的评价指标,确定最优滤波方法,以真实泄漏音波信号验证最优滤波方法的有效性。
进一步,在步骤一中,真实泄漏音波信号的拟合为多项式拟合,拟合公式:p(n)=a0+a1t+a2t2+…+antn,式中,p(n)为声压,单位kPa,t为泄漏时间,单位s,a0至an为多项式系数,n为阶数,取值为5~10。
进一步,在步骤二中,添加噪声信号的方法为:采集得到输气管线稳定运行时的音波信号,计算得到音波信号的均值和方差,在拟合泄漏音波信号基础上,添加与管线稳定运行音波信号具有相同均值和方差的高斯白噪声。
进一步,在步骤三中,预选择的滤波方法有:傅里叶变换滤波、小波滤波、中值滤波、维纳滤波、FIR滤波器滤波和IIR滤波器滤波。
进一步,预选择的滤波方法的特征参数为:
(1)傅里叶变换滤波通带频率0Hz~100Hz;
(2)小波滤波中小波基为sym8小波,分解层数为5层,阈值选取原则为rigrsure规则,阈值处理原则为二次多项式差值法;
(4)维纳滤波窗口长度为个采样点,其中fs为采样频率,单位Hz;
(5)FIR滤波器选用Kaiser窗,通带截止频率为30Hz,阻带起始频率为60Hz,通带波纹最大值为0.01,阻带波纹最大值为0.1;
(6)IIR滤波器选用ChebyshevⅡ型滤波器,通带截止频率为30Hz,通带内最大衰减为0.5dB,阻带截止频率为90Hz,阻带内最小衰减为80dB。
进一步,在步骤四中,综合分析比较各种方法的实际滤波效果,计算滤波后的评价指标,滤波信号评价指标为信噪比、均方根误差和信号泄漏下降沿幅值损失百分比。
进一步,信噪比计算公式为:均方根误差计算公式:信号泄漏下降沿幅值损失计算公式为:式中x(i)为拟合泄漏音波信号幅值,单位kPa, 为滤波后的伪泄漏音波信号幅值,单位kPa,minx(i)为拟合泄漏 音波信号下降沿幅值极小值,单位kPa,为滤波后的伪泄漏音波信号下降沿幅值极小值,单位kPa。
进一步,在步骤四中,最优滤波方法满足:信噪比评价指标最大、均方根误差评价指标最小并且泄漏下降沿幅值损失评价指标最小。
本发明提供的输气管线泄漏音波信号的滤波优选方法,通过满足信噪比评价指标最大、均方根误差评价指标最小并且泄漏下降沿幅值损失评价指标最小,选择最优滤波方法,填补了现阶段选择泄漏音波信号滤波方法的空白,定量地评价了泄漏音波信号的滤波效果,最大限度地增大了泄漏音波信号的信噪比,提高了音波法泄漏检测技术的适用性。本发明方法简单,操作方便,较好的解决了现阶段缺乏选择泄漏音波信号滤波方法的问题。
附图说明
图1是本发明实施例提供的输气管线泄漏音波信号的滤波优选方法的流程图;
图2是本发明实施例提供的泄漏音波信号滤波及滤波优选方法原理流程图;
图3是本发明实施例提供的拟合泄漏音波信号示意图;
图4是本发明实施例提供的伪泄漏音波信号示意图;
图5a是本发明实施例提供的经傅里叶滤波后的伪泄漏音波信号示意图;
图5b是本发明实施例提供的经小波滤波后的伪泄漏音波信号示 意图;
图5c是本发明实施例提供的经中值滤波后的伪泄漏音波信号示意图;
图5d是本发明实施例提供的经维纳滤波后的伪泄漏音波信号示意图;
图5e是本发明实施例提供的经FIR滤波器滤波后的伪泄漏音波信号示意图;
图5f是本发明实施例提供的经IIR滤波器滤波后的伪泄漏音波信号。示意图;
图6是本发明实施例提供的经小波滤波后的实际泄漏音波信号示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
下面结合附图及具体实施例对本发明的应用原理作进一步描述。
如图1所示,本发明实施例的输气管线泄漏音波信号的滤波优选方法包括以下步骤:
S101:在泄漏音波信号下降沿长度内,拟合真实泄漏音波信号时域波形,得到拟合泄漏音波信号;
S102:添加噪声信号,获得伪泄漏音波信号;
S103:利用预选择的各种滤波方法对伪泄漏音波信号进行滤波处理,得到滤波后波形;
S104:计算并获得滤波后信号的评价指标,确定最优滤波方法,以真实泄漏音波信号验证最优滤波方法的有效性;
如图2所示,本发明的具体实施流程为:
音波传感器采集输气管线稳定运行音波信号与泄漏音波信号,输气管线稳定运行时的音波信号视为背景噪声信号;
选择泄漏音波信号下降沿长度,即从信号幅值突降到幅值恢复至零值附近波动的长度作为窗口长度,采用多项式拟合真实泄漏音波信号时域波形,获得拟合泄漏音波信号,多项式拟合公式为:p(n)=a0+a1t+a2t2+…+antn,式中,p(n)为声压,单位kPa,t为泄漏时间,单位s,a0至an为多项式系数,n为阶数,取值为5~10,图3为根据真实信号拟合获得的拟合泄漏音波信号;
计算背景噪声信号均值和方差,在拟合得泄漏音波信号基础上,添加与管线背景噪声信号具有相同均值和方差的高斯白噪声,构造获得伪泄漏音波信号,图4为实际构造的伪泄漏音波信号;
预选择现阶段通用的傅里叶变换滤波、小波滤波、中值滤波、维纳滤波、FIR滤波器滤波和IIR滤波器滤波六种滤波方法分别对伪泄漏音波信号进行滤波,傅里叶变换滤波通带频率0~100Hz;小波滤波中小波基为sym8小波,分解层数为5层,阈值选取原则为rigrsure规则,阈值处理原则为二次多项式差值法;中值滤波窗口长度为个采样点,其中fs为采样频率,单位Hz;维纳滤波窗口长度为个采样点,其中fs为采样频率,单位Hz;FIR滤波器选用Kaiser窗,通带截止频率为30Hz,阻带起始频率为60Hz,通带波纹最大值为0.01,阻带波纹最大值为0.1;IIR滤波器选用ChebyshevⅡ型滤波器,通带截止频率为30Hz,通带内最大衰减为0.5dB,阻带截止频率为90Hz,阻带内最小衰减为80dB。图5a~5f为各种滤波方法的实际滤波波形图;
选择信噪比(SNR)、均方根误差(RMSE)和信号泄漏下降沿幅值损失百分比作为滤波效果的评价指标,信噪比计算公式为: 均方根误差计算公式: 信号泄漏下降沿幅值损失计算公式为:式中x(i)为拟合泄漏音波信号幅值,单位kPa,为滤波后的伪泄漏音波信号幅值,单位kPa,minx(i)为拟合泄漏音波信号下降沿幅值极小值,单位kPa,为滤波后的伪泄漏音波信号下降沿幅值极小值,单位kPa,表1为某一实际的泄漏音波信号,运用六种滤波方法进行滤波后的实际评价指标,
表1六种滤波方法滤波效果评价指标
对比表1中的各项数据,小波滤波的信噪比指标最大,均方根误差指标最小,下降沿幅值损失百分比指标最小,小波滤波具有最优滤波效果,选择小波滤波作为最优滤波方法;
利用小波滤波方法对实际泄漏音波信号进行滤波处理,验证滤波效果的实际有效性,图6为经小波滤波后的实际泄漏音波信号,滤波后,泄漏下降沿幅值与平稳噪声幅值比值增大,泄漏特征得到突出,泄漏下降沿噪声干扰减小,泄漏发生时刻清晰明了,小波滤波对实际泄漏音波信号具有良好的滤波效果,选择小波滤波作为此工况条件运行输气管道的最优滤波方法。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种输气管线泄漏音波信号的滤波优选方法,其特征在于,该输气管线泄漏音波信号的滤波优选方法包括以下步骤:
步骤一,泄漏音波信号下降沿长度为信号幅值突降至幅值恢复至零值的信号长度;拟合真实泄漏音波信号时域波形,得到拟合泄漏音波信号;
步骤二,添加噪声信号,获得伪泄漏音波信号;
步骤三,利用预选择的各种滤波方法对伪泄漏音波信号进行滤波处理,得到滤波后波形;
步骤四,计算并获得滤波后信号的评价指标,确定最优滤波方法,以真实泄漏音波信号验证最优滤波方法的有效性。
2.如权利要求1所述的输气管线泄漏音波信号的滤波优选方法,其特征在于,在步骤一中,真实泄漏音波信号的拟合为多项式拟合,拟合公式:p(n)=a0+a1t+a2t2+…+antn,式中,p(n)为声压,单位kPa,t为泄漏时间,单位s,a0至an为多项式系数,n为阶数,取值为5~10。
3.如权利要求1所述的输气管线泄漏音波信号的滤波优选方法,其特征在于,在步骤二中,添加噪声信号的方法为:采集得到输气管线稳定运行时的音波信号,计算得到音波信号的均值和方差,在拟合泄漏音波信号基础上,添加与管线稳定运行音波信号具有相同均值和方差的高斯白噪声。
4.如权利要求1所述的输气管线泄漏音波信号的滤波优选方法,其特征在于,在步骤三中,预选择的滤波方法有:傅里叶变换滤波、小波滤波、中值滤波、维纳滤波、FIR滤波器滤波和IIR滤波器滤波。
5.如权利要求4所述的输气管线泄漏音波信号的滤波优选方法,其特征在于,预选择的滤波方法的特征参数为:
(1)傅里叶变换滤波通带频率0Hz~100Hz;
(2)小波滤波中小波基为sym8小波,分解层数为5层,阈值选取原则为rigrsure规则,阈值处理原则为二次多项式差值法;
(5)FIR滤波器选用Kaiser窗,通带截止频率为30Hz,阻带起始频率为60Hz,通带波纹最大值为0.01,阻带波纹最大值为0.1;
(6)IIR滤波器选用ChebyshevⅡ型滤波器,通带截止频率为30Hz,通带内最大衰减为0.5dB,阻带截止频率为90Hz,阻带内最小衰减为80dB。
6.如权利要求1所述的输气管线泄漏音波信号的滤波优选方法,其特征在于,在步骤四中,综合分析比较各种方法的实际滤波效果,计算滤波后的评价指标,滤波信号评价指标为信噪比、均方根误差和信号泄漏下降沿幅值损失百分比。
8.如权利要求1所述的输气管线泄漏音波信号的滤波优选方法,其特征在于,在步骤四中,最优滤波方法满足:信噪比评价指标最大、均方根误差评价指标最小并且泄漏下降沿幅值损失评价指标最小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310468802.0A CN103644460A (zh) | 2013-10-09 | 2013-10-09 | 一种输气管线泄漏音波信号的滤波优选方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310468802.0A CN103644460A (zh) | 2013-10-09 | 2013-10-09 | 一种输气管线泄漏音波信号的滤波优选方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103644460A true CN103644460A (zh) | 2014-03-19 |
Family
ID=50249718
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310468802.0A Pending CN103644460A (zh) | 2013-10-09 | 2013-10-09 | 一种输气管线泄漏音波信号的滤波优选方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103644460A (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105534600A (zh) * | 2015-12-29 | 2016-05-04 | 高建萍 | 一种智能检验采血台 |
CN105597124A (zh) * | 2015-12-29 | 2016-05-25 | 杨旭平 | 一种防感染净化式呼吸智能护理装置 |
CN105662546A (zh) * | 2016-01-04 | 2016-06-15 | 中国人民解放军北京军区总医院 | 一种彩超骨科穿刺双引导控制*** |
CN106510991A (zh) * | 2016-12-21 | 2017-03-22 | 王秀峰 | 一种多功能轮椅 |
CN106678552A (zh) * | 2017-01-05 | 2017-05-17 | 北京埃德尔黛威新技术有限公司 | 一种新型渗漏预警方法 |
CN106764468A (zh) * | 2017-01-05 | 2017-05-31 | 北京埃德尔黛威新技术有限公司 | 一种渗漏预警***及自适应频谱消噪方法 |
CN106841215A (zh) * | 2017-02-23 | 2017-06-13 | 深圳市地铁集团有限公司 | 基于互联网土木工程结构表现缺陷评定双ccd检测*** |
CN107357788A (zh) * | 2017-07-13 | 2017-11-17 | 吉林农业科技学院 | 一种英语会议翻译*** |
CN108801363A (zh) * | 2018-07-31 | 2018-11-13 | 深圳大图科创技术开发有限公司 | 一种环境保护*** |
CN108872516A (zh) * | 2018-07-31 | 2018-11-23 | 梧州市兴能农业科技有限公司 | 一种水质监测预警*** |
CN110953485A (zh) * | 2019-12-06 | 2020-04-03 | 北京无线电计量测试研究所 | 一种燃气管线泄漏点定位方法和*** |
CN113217823A (zh) * | 2021-03-10 | 2021-08-06 | 北京航空工程技术研究中心 | 一种用于航煤管道泄漏监测和定位的***及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6389881B1 (en) * | 1999-05-27 | 2002-05-21 | Acoustic Systems, Inc. | Method and apparatus for pattern match filtering for real time acoustic pipeline leak detection and location |
CN101196872A (zh) * | 2007-11-19 | 2008-06-11 | 清华大学 | 基于压力和声波信息融合的泄漏检测定位方法 |
CN201188051Y (zh) * | 2008-01-18 | 2009-01-28 | 北京世纪华扬能源科技有限公司 | 音波泄漏检测定位装置 |
-
2013
- 2013-10-09 CN CN201310468802.0A patent/CN103644460A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6389881B1 (en) * | 1999-05-27 | 2002-05-21 | Acoustic Systems, Inc. | Method and apparatus for pattern match filtering for real time acoustic pipeline leak detection and location |
CN101196872A (zh) * | 2007-11-19 | 2008-06-11 | 清华大学 | 基于压力和声波信息融合的泄漏检测定位方法 |
CN201188051Y (zh) * | 2008-01-18 | 2009-01-28 | 北京世纪华扬能源科技有限公司 | 音波泄漏检测定位装置 |
Non-Patent Citations (1)
Title |
---|
方丽萍,钟莹莹,梁金禄,黄建德: ""输气管道泄漏音波信号的几种滤波算法"", 《钦州学院学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105534600A (zh) * | 2015-12-29 | 2016-05-04 | 高建萍 | 一种智能检验采血台 |
CN105597124A (zh) * | 2015-12-29 | 2016-05-25 | 杨旭平 | 一种防感染净化式呼吸智能护理装置 |
CN105662546A (zh) * | 2016-01-04 | 2016-06-15 | 中国人民解放军北京军区总医院 | 一种彩超骨科穿刺双引导控制*** |
CN106510991A (zh) * | 2016-12-21 | 2017-03-22 | 王秀峰 | 一种多功能轮椅 |
CN106678552A (zh) * | 2017-01-05 | 2017-05-17 | 北京埃德尔黛威新技术有限公司 | 一种新型渗漏预警方法 |
CN106764468A (zh) * | 2017-01-05 | 2017-05-31 | 北京埃德尔黛威新技术有限公司 | 一种渗漏预警***及自适应频谱消噪方法 |
CN106841215A (zh) * | 2017-02-23 | 2017-06-13 | 深圳市地铁集团有限公司 | 基于互联网土木工程结构表现缺陷评定双ccd检测*** |
CN107357788A (zh) * | 2017-07-13 | 2017-11-17 | 吉林农业科技学院 | 一种英语会议翻译*** |
CN108801363A (zh) * | 2018-07-31 | 2018-11-13 | 深圳大图科创技术开发有限公司 | 一种环境保护*** |
CN108872516A (zh) * | 2018-07-31 | 2018-11-23 | 梧州市兴能农业科技有限公司 | 一种水质监测预警*** |
CN110953485A (zh) * | 2019-12-06 | 2020-04-03 | 北京无线电计量测试研究所 | 一种燃气管线泄漏点定位方法和*** |
CN113217823A (zh) * | 2021-03-10 | 2021-08-06 | 北京航空工程技术研究中心 | 一种用于航煤管道泄漏监测和定位的***及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103644460A (zh) | 一种输气管线泄漏音波信号的滤波优选方法 | |
CN104595729B (zh) | 一种基于声波幅值的油气管道泄漏定位方法 | |
CN104595730B (zh) | 一种基于声波幅值衰减模型的油气管道泄漏定位方法 | |
CN104747912B (zh) | 流体输送管道泄漏声发射时频定位方法 | |
CN101592288B (zh) | 一种管道泄漏辨识方法 | |
Pal et al. | Detecting & locating leaks in water distribution polyethylene pipes | |
CN104535275B (zh) | 基于气泡声学的水下气体泄漏量的检测方法和检测装置 | |
CN101551064A (zh) | 供水管道泄漏检测定位的信号处理方法 | |
CN104819766A (zh) | 基于谐噪比的包络解调频带确定方法 | |
CN104089186B (zh) | 一种基于组合滤波和动态阈值的管道压力异常诊断方法 | |
CN110469782B (zh) | 一种基于自适应滤波器的管道泄漏定位装置 | |
CN110645483B (zh) | 一种基于频谱分析的城市埋地管道早期泄漏诊断方法 | |
CN108332063A (zh) | 一种基于互相关的管道泄漏定位方法 | |
CN103234121A (zh) | 基于音波信号的输气管道泄漏检测装置及检测方法 | |
CN104535281B (zh) | 基于压力相关法的机坪管网密闭性测试装置进行测试的方法 | |
CN116577037B (zh) | 一种基于非均匀频率语谱图的风管泄漏信号检测方法 | |
CN106289121A (zh) | 一种变径管等效管长的计算方法 | |
CN108252708B (zh) | 一种油井动液面识别方法 | |
CN105116445A (zh) | 一种水陆检波器地震数据合并处理的方法及装置 | |
CN107389787B (zh) | 一种基于频域的气侵检测信号处理方法 | |
CN107389786B (zh) | 一种基于时域的气侵检测信号处理方法 | |
CN104869519A (zh) | 一种测试麦克风本底噪声的方法和*** | |
CN108758358B (zh) | 一种基于声回波信号解调的管道泄漏监测方法及*** | |
CN106678553B (zh) | 一种泄漏动态压力波在管内气体中传播速度的计算方法 | |
Wu et al. | Novel signal denoising approach for acoustic leak detection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140319 |