纤维过滤介质非稳态过滤压力损失的研究.pdf

《纤维过滤介质非稳态过滤压力损失的研究.pdf》由会员分享，可在线阅读，更多相关《纤维过滤介质非稳态过滤压力损失的研究.pdf（5页珍藏版）》请在三一文库上搜索。

1、第3 0 卷第4 期纺织学报 V 0 1 3 0N o 4 2 0 0 9 年4 月J o u r n a lo fT e x t i l eR e s e a r c h A p r 2 0 0 9 文章编号：0 2 5 3 - 9 7 2 1 ( 2 0 0 9 ) 0 4 - 0 0 1 5 0 5 纤维过滤介质非稳态过滤压力损失的研究 付海明 ( 东华大学环境科学与工程学院，上海2 0 1 6 2 0 ) 摘要考虑粉尘颗粒沉积形成滤饼对过滤压力产生的影响，将滤饼层视为新增加的纤维捕集体，通过类比清洁 状态下过滤压力损失与含尘状态下纤维过滤压力损失的关系以及量纲分析的方法，推导非稳态过滤

2、状态下过滤压 力损失的数学模型，得出相关实验系数，且数学模型计算结果与实测结果基本吻合。通过测试滤料非稳态过滤的 压力损失，给出常用过滤介质非稳态过滤的压力损失增量计算的半经验公式，以便对过滤压力损失变化及运行时 间进行预测，及时清灰或更换过滤介质，延长滤袋使用寿命。 关键词纤维过滤介质；非稳态过滤；压力损失；数学模型 中圈分类号：X5 1 3文献标志码：A S t u d yo np r e s s u r ed r o po fn o n - s t a t i o n a r yf i l t r a t i o nw i t hf i b e rm e d i a F UH a i m

3、 i n g ( C o l l e g eo fE n v i r o n m e n t a lS c i e n c ea n dE n g i n e e r i n g ，D o n g h u aU n i v e r s i t y ，S h a n g h a i 2 0 1 6 2 0 。C h n a ) A b s t r a c tS i n c et h ed u s tc a k ef o r m e df r o mt h es e t t l i n go fd u s tp a r t i c l e se x e r ti m p a c to nt h e

4、f i l t r a t i o np r e s s u r e l o s s ，t h ed u s tc a k ei sr e g a r d e da san e wc o l l e c t o r T h ef i l t r a t i o np r e s s u r el o s su n d e rc l e a nc o n d i t i o ni s c o m p a r e dw i t ht h a tu n d e rd u s t yc o n d i t i o n ，a n dam a t hm o d e lo ff i l t r a t i

5、o np r e s s u r ed r o po fn o n s t a t i o n a r yi s d e v e l o p e d T h ee x p e r i m e n t a ld a t ar e v e a l t h a tt h et h e o r e t i c a lv a l u ec o m p u t e db yu s i n gt h em a t hm o d e li s b a s i c a l l yi na g r e e m e n t 诵t l lt h ee x p e r i m e n t a lv a l u e B

6、a e s e do nt h e t e s to fp r e s s u r el o s si nn o n - s t a t i o n a r y f i l t r a t i o n ，as e m i e m p i r i c a lf o r m u l af o rc o m p u t a t i o no ft h ep r e s s u r el o s si n c r e m e n ti nn o n s t a t i o n a r yt i t r a t i o n w h e nu s i n gc o m m o nh o m em a d

7、ef i l t r a t i o nm e d i as oa st oe s t i m a t et h ec h a n g eo ft h ef i l t r a t i o np r e s s u r el o s sa n d d u r a t i o no fo p e r a t i o n T i m e l yc l e a n i n go rr e p l a c i n gt h ef i l t r a t i o nm e d i ac a np r o l o n gt h es e r v i c el i f eo f t h e f i l t

8、r a t i o nb a g K e yw o r d s f i b e rf i l t e rm e d i a ；n o n s t a t i o n a r yf i l t r a t i o n ；p r e s s u r ed r o p ；m a t hm o d e l 研究非稳态过滤压力损失随时间变化的规律， 对降低过滤压力损失，优化过滤过程及减少能耗是 至关重要的n “ 。本文主要考虑粉尘颗粒沉积形成 的滤饼对过滤压力损失产生的影响，将滤饼层视为 新增加的纤维捕集体，通过清洁状态过滤压力损失 与含尘状态的类比分析，引入过滤压力损失系数及 纤维直径增量的物理概念，

9、推导得出非稳态纤维过 滤介质过滤压力损失计算数学模型，并将其理论计 算结果与实验研究进行对比分析，得出相关实验系 数，推导出非稳态过滤压力损失增量计算的半经验 公式，可进行非稳态过滤压力损失的预测，以便优化 过滤运行，节省能耗。 1 非稳态过滤压力损失模型 在清洁状态下，过滤压力损失可表达为m 8 1 ，工 一P 。A = z f J 。F o A d 并= 以F o A L ( 1 ) F 。= c 加d 哦丢p 2 ( 2 ) 收稿日期：2 0 0 8 0 4 2 1修回日期：2 0 0 8 一0 8 2 l 基金项目：上海市重点学科建设资助项目( B 6 0 4 ) 作者简介：付海明(

10、1 9 6 2 一) ，男，高级工程师，博士。主要研究领域包括暖通空调、空气品质及有害气体控制。E - m a i l ：f i l m d h u e d u c n o 万方数据 PDF Watermark Remover DEMO : Purchase from www.PDFWatermarkR to remove the watermark 1 6 纺织学报第3 0 卷 式中：l ，为单位过滤器中纤维的长度，I n ，Z ，= 4 a i r e d ；c 为清洁时过滤器压力损失系数；d ，为 清洁状态下纤维直径，m ；l D 。为空气密度，k g ，m 3 ；口为 过滤速度，m s

11、 ；P 。为稳态过滤压力损失，P a ；F 。为 稳态纤维捕集体过滤压力损失，P a 。 粉尘颗粒沉积在纤维表面，使纤维直径增加及 粉尘层内部结构孔隙率产生变化，随过滤的进一步 运行，粉尘沉积逐渐在纤维表面形成滤饼，沉积的粉 尘可视为新增加的捕集体，因而增加了对粉尘颗粒 的捕集效率，而新增的纤维捕集体又不断捕集颗粒， 形成新的纤维捕集体阳 ，如此循环，导致过滤器捕集 效率及过滤压力损失增加。 类比清洁状态下过滤压力损失与纤维压力损失 的关系，在非稳态过滤状态下，过滤压力损失与纤维 捕集体也应满足如下关系n 引： 广 一P 。A = f f IF 。A d x ( 3 ) J 0 F 。= c

12、D m d 缸I D 。寺F 2 ( 4 ) 式中：P 。为非稳态状态下过滤压力损失，P a ；F 。为 非稳态状态下纤维捕集体过滤压力损失，P a ；C D m 为 非清洁状态下压力损失系数；d h 为非清洁状态下纤 维捕集体直径，m 。 对比式( 2 ) 、( 4 ) ，非清洁状态下纤维捕集体压力 损失可表达为 ，。= ( C D m C ) ( d 缸，d ，) F o ( 5 ) 将式( 1 ) 、( 5 ) 代人式( 3 ) 可得 n ，工 P 。= 了 - 0l ( c D m ，c ) ( d 缸d f ) d 茗 ( 6 ) “J0 式( 6 ) 为非稳态过滤压力损失的机制模型

13、，确定 压力损失系数及纤维直径增量与其他因素的相关关 系，即确定非稳态过滤状态下的过滤压力损失。 根据大量实验数据回归分析可知：压力损失系 数与过滤厚度成非线形关系；纤维捕集体直径的增 量与粉尘沉积质量成非线形关系 I 卜挖 ，因此，假设 C D m C m = 口l 石+ 0 2 ( 7 ) d 缸d f = C 1 , 31 7 1 , ”2 + 口4 ( 8 ) 粉尘沉积量 m = T C o v t ( 9 ) 式中：叩为捕集效率；C 。为过滤粉尘浓度，k g m 3 ； t 为过滤时间，8 ；m 为单位面积纤维过滤介质沉积 量，k s m 2 ；口I 、口2 、口3 、口4 、r g

14、 l 、n 2 为与实验相关的系 数。由文献 1 3 可得出 7 = kd ? 2 口勺p 驴毋t ( 1 0 ) 式中：d ，为粉尘颗粒粒径，m ；口为纤维过滤介质填 充率；L 为纤维过滤介质厚度，m ；k 、茹，、茹：、舅，、扎、 菇，、髫。、茹，为与实验相关的系数。 将式( 9 ) 、( 1 0 ) 代入式( 7 ) 、( 8 ) 可得 ( c ，c ) ( d h d f ) = k I 诉晓矿f 7 菇+ k 2 诉戗J t 7 + J 3 髫+ k 4 ( 1 1 ) 式中：k 。、k ：、k ，、k 。、d 、：I 、卢、r 均为与实验相关的系 数。将式( 1 1 ) 代入式(

15、6 ) 可得 P m - 钮( C D I C 捌d 缸d f ) d 菇= P o 【( 篙“4 e oe “ 羔+k1 】 ( 1 2 ) + 厅t 。 4 】 ” 假设 ：黑岫 弱：篙“ K l = g oP o 诉 则 P 。= K ：P o + K oP o 诉a 俨广= 蜀P o + K la r o t 7 ( 1 3 ) 当过滤时间为0 时，过滤介质的压力损失应等 于过滤介质清洁时的压力损失，即式( 1 3 ) 应满足磁 为1 的条件，因此，可得 P 。= P o + K l 成r ( 1 4 ) 非稳态过滤压力损失增量的表达式为 A p 。= P 。一P o = K l 晓扩

16、r ( 1 5 ) 式中：A p 。为非稳态过滤压力损失增量，P a ；K 。为过 滤介质特性影响系数，可通过实验确定。它们与填 充率口、纤维直径D ，、过滤介质厚度、黏性系数 t 等因素相关。 2 过滤压力损失实验 参照工业用袋式除尘器搭建实验台，滤袋形状 为圆形布袋；过滤方式为下进风外滤袋式除尘器，系 统风量通过S A M C O i 高功能静音式变频器调节，风 速变化区域大，为2 5 0 。29 5 0m r a i n ，可实现无级变 速。实验台示意图见图1 。 万方数据 PDF Watermark Remover DEMO : Purchase from www.PDFWaterma

17、rkR to remove the watermark 第4 期 付海明：纤维过滤介质非稳态过滤压力损失的研究 1 7 注：1 一滤袋；2 一支撑架；3 一粉尘装置；4 一数字微压计；5 一整流栅；6 - - - S w e m a A i r3 0 0 风速仪；7 一变频器；8 一风机；9 一电源。 图1 实验台示意图 F i g 1 S k e t c hm a po fe x p e r i m e n t a ls e t 3 结果与讨论 含尘过滤状态下，对粉尘浓度分别为2 0 、4 0 、6 0 、 8 0g ，m 3 进行过滤测试，可得到不同纤维过滤介质在 含尘过滤状态下过滤速度对

18、过滤压力损失影响的变 化曲线。过滤介质为2 0 8 滤料含尘气流中各影响因 素对过滤压力损失影响的实验结果如图2 4 所示。 图3 中，实验粉尘浓度为2 0g ，m 3 ，过滤速度为 O 9 8m m i 1 。 罢 j K 辐 农 田 图2 过滤压力损失与过滤速度 F i g 2 C u r v eo ff i l t r a t i o np m * s s u r ed r o pw i t hv e l o c i t y 可以看出：随着过滤速度的增加，过滤压力损失 上升。随着粉尘浓度的增加，滤筒压力损失增大。 随着时间的增加，压力损失也随之增大。在过滤初 期，压力损失迅速上升，过滤一

19、段时间后，压力损失 曲线斜率逐渐发生变化。压力损失曲线在一段弯曲 后，又近似呈线性趋势增加，曲线斜率很小，即随着 时间增加，压力损失以极其缓慢的速度上升。 由非稳态过滤压力损失模型式( 1 5 ) 可知：影响 含尘过滤介质滤筒过滤压力损失的主要参数有入口 粉尘浓度、过滤速度和过滤时间。本文通过对过滤 罢 水 辑 R 田 入口粉尘浓度( g m 一3 ) 图3 过滤压力损失与粉尘浓度 F i g 3 C u r v eo ff i l t r a t i o np r e s s u r ed r o pw i t h d u s tc o n c e n t r a t i o n 苣 水 鞲

20、 穴 田 图4 过滤压力损失与过滤时间 F i g 4 C u r v eo ff i l t r a t i o np r e s s u r ed r o pw i t hf d t r a t i o nt i m e 介质2 0 8 、过滤介质7 2 9 及覆膜过滤介质7 2 9 实验数 据进行多元回归，找出其压力损失变化规律，如表1 所示。对表1 进行非样本实验数据检测，给出过滤 介质2 0 8 的回归公式计算值与实测值，如表2 所示。 通过检测得出：对于过滤介质2 0 8 计算值与实 测值偏差为0 5 2 一4 8 8 ；对于过滤介质7 2 9 计 万方数据 PDF Waterma

21、rk Remover DEMO : Purchase from www.PDFWatermarkR to remove the watermark 1 8 纺织学报第3 0 卷 算值与实测值偏差为9 1 8 一1 6 8 9 ；对于覆膜过 滤介质7 2 9 计算值与实测值偏差为6 8 5 1 4 6 4 。 质K I 为7 7 1 2 4 、口为0 3 2 6 、卢为0 6 8 5 、r 为0 3 2 3 ； 对于覆膜7 2 9 过滤介质K 。为11 6 0 9 6 、口为0 3 5 3 、 p 为0 9 0 3 、r 为0 3 5 3 。 表1 含尘过滤介质压力损失多元回归公式 呦1M 砌V

22、 啦她n 删细山甜p 嘲md 呷 4 结论 o ft l l t r a t l o nH l e d l ai nd u s tc o n d i t i o n 注：实验过程过滤风速范围为8 5 5 1 3 9 2m m i n 。 裹2 过滤介质2 0 8 回归公式计算值与实测值对比 T 曲2C a l c u l a t er e s u l to fn l t r a t t o l lm e d i a2 0 8w i t h e x p e r i m e n t a lr e s u l t 由于风机的压头不变，当过滤粉尘负荷不断增 加时，滤料粉尘层不断增厚，滤饼过滤压力损失不

23、断 增加，因而导致过滤风速下降。清洁状态下过滤风 速与含尘状态过滤风速不同，应采用测量纤维过滤 介质压力损失增量进行回归分析。如果采用含尘运 行过滤后过滤介质的总压力损失进行回归分析，各 实验系数是不同的，此时，过滤压力损失约与过滤速 度的二次方成正比。文献 1 4 理论推导得出的过滤 压力损失等于过滤速度的一次方与过滤速度二次方 之和的表达式存在不完善之处，其表达式中的实验 系数依然与过滤速度相关，因此，通过理论推导及实 验验证，非稳态过滤压力损失如果表达为过滤速度 的单因素函数，表达式应为 P 。= a f ) + 6 t ，“( 1 6 ) 式中口、b 、n 为与过滤速度不相关的实验系数

24、，其值 可由实验确定。 非稳态过滤的过滤介质压力损失增量计算模 型为 A P 。= K lC o 矿f( 1 7 ) 式中：K 。为压力损失系数，随过滤介质材料不同而 变化，其数值可通过实验确定。它们与填充率口、纤 维直径D ，、过滤介质厚度、黏性系数口等因素相 关。对于国产2 0 8 过滤介质墨为7 7 8 2 1 、口为 0 3 3 1 、卢为0 7 7 9 、r 为0 3 4 8 ；对于国产7 2 9 过滤介 通过建立过滤介质非稳态过滤的机制模型，根 据对实验结果的分析，假设过滤压力损失系数与过 滤介质厚度呈非线形关系，在粉尘负荷状态下纤维 直径的增量与粉尘沉积量呈非线形关系，经过理论

25、推导及实测可得出以下结论。 1 ) 过滤压力损失与粉尘量及过滤时间呈非线 性变化关系，粉尘发尘量越大，过滤时间越长，其过 滤压力损失越大，随粉尘量及粉尘过滤时间增加，过 滤压力损失的增加逐渐平缓。过滤压力损失随过滤 平均速度的增加而增加，过滤压力损失随粉尘浓度 增加而增加。 2 ) 非稳态过滤。介质压力损失计算模型也可以 采用过滤速度单因素表达式P 。= 口口+ b y “，但此过 滤压力损失表达式不能预测过滤清灰的周期。 参考文献： Y O S H I Y U K IE n d o ，C H E ND a r e n ，D A V I D ，e ta 1 E f f e c t s o fp

26、 a r t i c l e sp o l y d i s p e r s i t ya n ds h a p e sf a c t o r sd u r i n gd u s t c a k el o a d i n gO i la i rf i l t e r s J P o w e r sT e c h n o l o g y ，1 9 9 8 ， 9 8 ：2 4 l 一2 4 9 1 H O M A SD ，C O N r A LP ，R E N A U D I NV ，e ta 1 M o d d l i J l g p r e s s u r ed r o pi nH E P A

27、f i l t e r sd u r i n gc l y n a I I l i cf i l t r a t i o n J A e r o s o lS e i ，1 9 9 9 ，3 0 ( 2 ) ：2 3 5 2 4 6 S T O C K L M A Y E RC h H O F U N G E RW S i m u l a t i o no ft h e f i l t r a t i o nb e h a v i o u ro fd u s tf j d t e r s J S i m u l a t i o nP a r e i t i e e a n dT h e o r

28、 y ，1 9 9 8 ，6 2 8 1 2 9 6 T H O M A SD ，P E N I C O TP ，C O N T A LP ，e ta 1 C l o g g i n go f f i b r o u sf i l t e r sb ys o l i da e r o s o lp a r t i c l e se x p e r i m e n t a la n d m o d e l i n gs t u d y J C h e m i c a lE n g i n e e r i n gS c i e n c e ，2 0 0 1 ， 5 6 ：3 5 4 9 3 5 6

29、 1 D E C L 怠NC W A L S H R e c e n ta d v a n c ei nt h eu n d e r s t a n d i n g o ff i b r o u sf i l t e rb e h e a v e ru n d e rs o l i dp a r t i c l el o a d JJ F i l t r a t i o na n dS e p a r a t i o n ，1 9 9 6 ( 7 ) ：5 0 1 5 0 6 S I 0 C K L M A Y E RC h H O F U N G E RW S i m 蛳O no f t h

30、 e r e g e n e r a t i o no fd u s tr t l t e r s J 1 l a t h e m a t i e sa n dC o m p u t e r s i nS i m u l a t i o n ，1 9 9 8 。4 6 ：6 0 l 一6 0 9 C H I K A 0g a m a o k a P e r f o r m a n c eo f 蚰a i rf i l t e ra td u s t - l o a d e dc o n d i t i o n I V l A d v a n c e si nA e r o s o lF i l

31、 t r a t i o n N e w Y o r k ：L e n i sP u b l i s h e s ，3 2 3 3 3 5 ? C H I K A 0K a n a o k a S O T O J Ir n r a g i P l s S i l l ed r o po fa i r ( 下转第2 3 页) 1 J 1 J 1 J 1 J 1 J 1 J 1 J 1 J l 2 3 4 5 6 7 8 r L r L r L r L r L r L r L r L 万方数据 PDF Watermark Remover DEMO : Purchase from www.PDFWa

32、termarkR to remove the watermark 第4 期哈丽丹买买提等：纤维素氨基甲酸酯溶液黏度与其结构及组成的关系 。2 3 7 8 r e g e n e r a t e dc e l l u l o s ef i b r e s J J o u r n a lo fC h i n aT e x t i l e U n i v e r s i t y ，1 9 9 8 ，2 4 ( 4 ) ：1 1 1 一1 1 4 哈丽丹买买提，努尔买买提，吾满江艾力植物纤维 素聚合度的测定【J 合成纤维工业，2 0 0 6 ，2 9 ( 1 ) ： 4 0 4 2 M A M A T

33、H a l i d a n 。N U R M M A M A T ，A I L IW u m a n j i a n g D e t e r m i n a t i o no fp o l y m e r i z a t i o no fc e l l u l o s eb yv i s c o s i t yt t j C h i n aS y n t h e t i cF i b e rI n d u s t r y ，2 0 0 6 ，2 9 ( 1 ) ：4 0 4 2 哈丽丹买买提，夏木西卡玛尔，吾满江艾力植物纤 维素的限制性酶解特性 J 化学与生物工程，2 0 0 6 。 2 3

34、( 6 ) ：2 5 2 6 M A M A TH a l i d a n ，X A M U X I K A M A R 。A I UW u m a n j i a n g R e s t r i c t i v ee n z y m o o l y s i sp r o p e r t i e so fp l a n tc e n u l o w i t h 9 1 0 c e l h l a s e J C h e I I I i s t r y B i o e n g i n e e r i n g ，2 0 0 6 ，2 3 ( 6 ) ： 2 5 2 6 哈丽丹买买提，努尔买买提，吾

35、满江艾力纤维素氨 基甲酸脂的合成 J 高分子材料科学与工程，2 0 0 5 ， 2 1 ( 4 ) ：1 0 2 1 0 4 M A M A TH a l i d a n ，N U R M M A M A T 舡U W u m a n j i a n g S y n t h e s i so fc e U u l o 靶e a r b a m a t e J P o l y m e rM a t e r i l s S c i e n c ea n dE n g i n e e r i n g ，2 0 0 5 ，2 1 ( 4 ) ：1 0 2 1 0 4 哈丽丹买买提，吾满江艾力纤维素氨基

36、甲酸酯溶 液制备工艺研究 J 纺织学报，2 0 0 8 ，2 9 ( 1 ) ：5 8 M A M A TH a h d a n 。A I uW u m a n j i a n g S t u d yo np r e p a r a t i o n o fe e l h l o s ec a r b a m as o l u t i o n J J o u r n a lo fT e x t i l e R e s e a r c h 。2 0 0 8 ，2 9 ( 1 ) ：5 8 删州州o 删州_ ooqh”o_“洲o“o州删1州o帝 ( 上接第1 8 页) f i l t e rw i

37、t hd u s tl o a d J JA e r o s o lS c i ，1 9 9 0 ，2 1 ( 1 ) ： S e p a r a t i o n ，1 9 8 6 ( 1 ) ：1 1 3 1 1 8 1 2 7 1 3 7 1 3 付海明，张鹏峰，亢燕铭空气过滤捕集效率影响因素 9 G A B E LIT a r d o s S e p a r a t i o no fa i r b e m ed u s ti nd e e p - b e d分析及多元关联式的确定 J 洁净与空调技术，2 0 0 6 f i l t r a t i o n M A d v a n c e

38、si nA e r o s o lF i l t r a t i o n N e w , Y o r k ：( 2 ) ：2 5 2 8 L e w i spublishes，2000：241257FUH a i m i n g ，Z H A N GP e n g f e n g ，K A N GY a n m i n g A n a l y s e s 1 0 V I N C E N TJN o v i c k ，J A N EFK l t l a e n P r e d i c t i n gp r e s s u r ea n dd e f i n i t i o n so fc a l

39、 c u l a t i n g c o l l e c t i n ge f f i c i e n c y o f r e s p o n s ec h a r a c t e r i s t i c sa c l r 0 6 sp a r t i c l e l o a d e df i l t e r s M f l i e r s “ J C o n t a m i n a t i o na n dA i r - c o n d i t i o n l n gT e c h n o l o g e ， A d v a n c ei nA e r o s o lF i l t r a

40、t i o n N e wY o r k ：L e w i sp u b l i s h e s ， 2 0 0 6 ( 2 ) ：2 5 2 8 2 0 0 0 ：3 3 7 3 4 7 1 4 向晓东，朱东升，彭伟功非稳态过滤除尘压力损失的 1 1 C H E N GY u h s i a n g ，T S MC h u e n j i n n F a c t o r si n f l u e n c i n g研究 J 工业安全与环保，2 0 0 2 ( 9 ) ：5 7 p r e s s u r ed r o pt h r e u g h tad u s tc a k ed u r

41、i n gf i l t r a t i o n J X I A N GX i a o d o n g ，Z H UD o n g s h e n g ，P E N GW e i g o n g S t u d y A e r o s o ls c i e n c ea n dt e c h n o l o g y ，1 9 9 8 ，2 9 ( 4 ) ：3 1 5 3 2 8 o nt h ep r e s s u r ed r o po fd y n a m i cf i l t r a t i o no ff i l t e r s J 1 2 B O T H W E L LE A n

42、a n a l y s i so ff a b r i cd u s tf i l t r a t i o nI ：m o d e lI n d u s t r i a lS a f e t ya n dE n v i r o n m e n t a lP r o t e c t i o n ，2 0 0 2 ( 9 ) ： o b s e r v a t i o n so fd u s tc a k ef o r m a t i o n J F i l t r a t i o na n d5 7 万方数据 PDF Watermark Remover DEMO : Purchase from www.PDFWatermarkR to remove the watermark