天然气脱蜡旋风分离器分离效率的模拟.pdf

化 工 进 展 2 0 1 4年第 3 3卷第 1 期 C HE MI C AL I ND US T R Y A ND E NGINE E R I NG P R O GR E S S 4 3 天然气脱蜡旋风分离器分离效率的模拟 袁惠新，方毅，付双成 常州大学机械工程学院，江苏 常州 2 1 3 0 1 6 摘要旋风分离器分离效率高，不易堵塞，用于天然气脱蜡效果显著。通过 C F D软件 F l u e n t 模拟 C Y G ． S型天 然气脱蜡旋风分离器的两相流场，得到了旋风分离器内的压力、切向速度、轴向速度分布。对比了不同入 口速 度下的模拟与理论计算的分割粒径 X 发现具有很好的吻合度，两相模拟有一定的可靠性。结果表明在旋风 分离器锥段底部靠近壁面处的石蜡液滴质量浓度较 高；随着进 口流量的增加，旋风分 离器分离效率提 高，当进 口流量为 1 0 0 0 m ／ h时， X 5 o 可以达到 5 ． 3 m；大粒径液滴的分离效果明显，但在所研究的进 口流量范围内，进口 流量的变化不能明显地影响粒径小于5 g m液滴的分离效率；柱段和锥段长度的增加使得旋风分 离器的整体长度 增加，延长了液滴在旋风分离器内的停留时间，提 高了旋风分离器的分离效率。 关键词旋风分离器；数值模拟；分离效率；进 口流量；柱段和锥段长度 中图分类号T Q 0 5 1 ． 8 文献标志码A 文章编号1 0 0 06 6 1 3 2 0 1 40 1 0 0 4 30 7 DoI 1 0 ． 3 9 6 9 ． i s s n ． 1 0 0 0 6 6 1 3 ．2 0 1 4 ． 0 1 ． 0 0 8 Nume r i c a l s i m ul a t i o n o f t he s e pa r a t i o n e ffi c i e n c y o f na t ur a l g a s de wa x i ng c y c l o ne s e pa r a t o r Y U A N Hu i x i n，F ANG Y i ，FUS h u a n g c h e n g S c h o o l o f Me c h a n i c a l E n g i n e e r i n g，Ch a n g z h o u Un i v e r s i t y ，Ch a n g z h o u 2 1 3 01 6，J i a n g s u，C h i n a Ab s t r a c t T h e c y c l o n e s e p a r a t o r h a s h i g h s e p a r a t i o n e ffic i e n c y a n d i s n o t e a s y t o b e b l o c k e d ， i t h a s g o o d e f f e c t o n t h e n a t u r a l g a s d e wa x i n g ． Nu me r i c a l s i mu l a t i o n o f t h e t wo p h a s e fl o w o f t h e CYG S n a tur a l g a s d e wa x i n g c y c l o n e s e p a r a t o r wa s s tud i e d b y u s i n g CF D s o f t wa r e F l u e n t ，a n d g e t t h e i n t e r n a l p r e s s u r e ，t a n g e n t a n d a x i a l v e l o c i t y d i s t r i b u t i o n s ． T h e s i mu l a t i o n a n d c a l c u l a t e d c u t s i z e X 5 0 i n d i f f e r e n t i n l e t v e l o c i ty a r e c o n t r a s t e d，t h e r e s u l t s h a s a g o o d d e g r e e o f a g r e e me n t ，S O t h e t wo - p h a s e s i mu l a t i o n a r e r e l i a b l e ． T h e r e s u l t s s h o w t h a t t h e ma s s c o n c e n t r a t i o n o f p a r a ffi n d r o p l e t i S h i g h e r a t t h e b o t t o m o f t h e c o n e s e c t i o n ． W i t h t h e i n c r e a s i n g o f t h e i n l e t fl o w，t h e s e p a r a t i o n e ffic i e n c y o f c y c l o n e s e p a r a t o r wa s i mp r o v e d ， u n d e r t h e 1 0 0 0 m ／ h i n l e t flo w ． t h e X 5 0 c a n r e a c h 5 ． 3 g m． Ho we v e r ， i n t h e s t u d i e d r a n g e o f i n l e t flo w，t h e i n c r e a s e r a t e o f t h e d r o p l e s s t h a n 5 g m i s n o t o b v i o u s ． T h e l e n g t h e n e d o f c o l u r n fi a n d t a p e r s e c t i o n i n c r e a s e s t h e l e n g t h o f c y c l o n e s e p a r a t o r ，e x t e n t s t h e r e s i d e n c e t i me o f t h e d r o p l e t s i n t h e c y c l o n e s e p a r a t o r ， a n d i mp r o v e s t h e s e p a r a t i o n e ffi c i e n c y o f t h e c y c l o n e s e p a r a t o r Ke y wor ds c y c l o n e s e pa r a t o r ；n ume r i c a l s i mu l a t i o n； s e pa r a t i o n e ffi c i e n c y； i n l e t flo w ； c o l umn a n d t a p e r s e c t i o n l e n g t h 石蜡是由正构烷烃 、带支链 或异构 烷烃和 环烷烃的组成 的复杂混合物 ，有很多品种 的熔点范 围在常温范围 2 3 6 ℃[ ，以固态和半固态的 的形式存在于天然气中。由于石蜡具有黏结和流动 性差的特性， 极易在管道内壁沉积， 甚至堵塞设备， 所 以对于天然气中石蜡的分离，传统的过滤分离会 非常容易堵塞过滤器 ，严重影响生产效率，存在安 收稿日期2 0 1 3 0 7 ． 1 6 ；修改稿 日期2 0 1 3 0 8 - 2 0 。 第一作者及联系人袁惠新 1 9 5 7 一 ，男，教授，博士生导师，从事 分离技术与设备研究。E ma i l y u a n h u i x i n 2 0 0 0 1 6 3 ． c o rn。 学兔兔 w w w .x u e t u t u .c o m 学兔兔 w w w .x u e t u t u .c o m 第 1 期 袁惠新等天然气脱蜡旋风分离器分离效率的模拟 表 1 旋风分离器的主要结构参数 结构参数 尺寸／ m m 旋流器直径 D 柱段长度 三 l 锥段长度 2 矩形进料口尺寸a x b 溢流管插入深度 L n 溢流口直径 D 。 底流口直径 D 分了 3 2 6 9 8个网格 。 1 ． 2 计算方法 在流场 的数值模拟 中对 于各个通用方程 的离 散， 目前 比较常用的差分格式主要有一阶迎风格式、 二阶迎风格式、Q UI C K格式和中心差分格式等 。当 流动具有强旋转特性时， QUI C K差分格式明显具有 较高的精度，能最准确地预报旋流器内的流场【 9 J 。 采用 F l u e n t P r e s s u r e B a s e d隐式求解器； 计算模型选 用 R S M 模型，控制方程中的压力一 速度耦合项选用 协调一致的 S I MP L E C算法。 1 ． 3 边界条件及部分物性参数 1 ． 3 ． 1 流体相 的边界条件 1 进 口边界旋风分离器入 口为流场 的速 度入 口，速度方 向垂直于入 口界面指向入 口内侧 ， 压力为该处静压。 2 出 口边界旋风分离器的溢流 口，为 自由 出口。底流 口设为壁面 ，为液滴捕集面之一。 3 壁面在旋风分离器的近壁区， 湍流输运 性质会发生急剧的变化 ，应用壁面函数法进行近壁 处理 ，解决计算精确度的问题。 4 基础物性旋风分离器 内的流体为天然 气 ，分散相为石蜡液滴，其组成及相关的物性参数 如表 2所示。 1 - 3 ． 2 分散相边界条件 液 滴粒 径 分 布通 过 对 液 滴尺 寸 进 行 R o s i n ． R a mml e r 分布拟合[ ] ，所有的颗粒尺寸范围被划分 为离散的尺寸组，每个尺寸组由该组射流源中的单 个颗粒所代表，在计算运动轨迹的时候所采用的就 表 2 物料及相关物性参数 嵝 g _ K ．_皿 图2 R o s i n ． R a mml e r 分布 是这个平均直径 。图 2 为模拟采用的液滴的粒径 分布。 溢流 口设置为逃逸 出口，当液滴运动到溢流 口 的时候会逃逸，此时停止对颗粒的追踪，终止对该 颗粒轨道的计算 。 壁面和底流 口均设置为捕集面 ，液滴接触壁面 或到达底流 口就认为液滴得到分离，终止对该颗粒 轨道 的计算 。 2 结果与分析 2 ． 1 压 力分布 图 3 a 为旋风分离器在进 口流量 Q i 为 5 0 0 m3 ／ h 的纵截面压力分布云图；图 3 b 为距离底流 口 9 0 0 mm 处横截面上的压力分布。压力呈对称分布 ，由 外至 内逐渐减小，壁面处压力最高，轴心处压力最 低，有负压产生。压力分布比较稳定，可以看出流 场比较稳定，有利于两相分离。 2 ． 2 速度分布 2 ． 2 ． 1 切 向速度 切 向速度反映了流体在旋风分离器内旋转的快 慢 ， 决定 了离心力 的大小，是使液滴沉降的驱动力。 图 4 a 为旋风分离器在进 口流量 Q i 为 5 0 0 m3 ／ h 时的纵截面切向速度分布云图；图 4 b 为距离底 流口9 0 0 1T U T I 处横截面上的切向速度分布。 模拟发现切 向速度满足一般旋风分离器 的准 自 由涡规律从边壁开始随着半径的减小，切 向速度 缓慢增大，在某一径 向位置达到最大值，然后急剧 降低，在整个直径范围内，呈现“ M” 形驼峰分布。 2 ． 2 ． 2 轴向速度 轴 向速度反映了流体 向两个轴 向出口的流动情 况。图5 a 为旋风分离器在进口流量 Q i 为 5 0 0 m V h 时的纵截面上的轴向速度分布云图和零轴速包络 面；图 5 b 为距离底流 口9 0 0 n l l T l 处横截面上的轴 O 0 2 嗡 0 0 ， 枷 m 主蓦 学兔兔 w w w .x u e t u t u .c o m 学兔兔 w w w .x u e t u t u .c o m 学兔兔 w w w .x u e t u t u .c o m 4 8 化 工 进 展 2 0 1 4年第 3 3卷 随着进 口流量 的增加，级效率 曲线 向左偏 移，分割粒径 奶0 减小，分级效率提高。这是因为 随着进口流量 的增加， 旋风分离器内流体的切向 速度增大，产生更大的离心分离作用，使得液滴更 快靠近边壁而被分离出。 在所研究的进 口流量范围内，对于小于 5 g m 的液滴，流量 的增加没有显著增加其分离效率。这 是由于粒度小的液滴沉 降速度小， 来不及完成沉降， 被直接夹带出溢流口，未实现分离。 当进 口流量 大于 8 0 0 m3 ／ h后， 级效率 曲线变 化不是很显著。这是 由于在此时旋流器内的湍流强 度增大，不利于液滴的分离，同时，湍流加速了液 滴的破碎，使液滴粒径变小，增加了分离难度。 2 ． 4 ． 2 段长度 1 和锥段长度 2 对分离效率 的影响 柱段长度和锥段长度增加，增加 了旋风分离器 的长度，延长了液滴在旋风分离器内的停留时间。 流体在旋风分离器内的停留时间决定于流体在旋风 分离器内运动的速度和经过的路径长短，或流体的 轴向速度 和路径的轴 向距离即旋风分离器的长 度 。旋风分离器内的物料的平均停留时间可以表 达为式 3 。 ∞ f 生 垒 ] 3 ∞ 一 W j nD 由式 3 可知，停 留时间与进 口流量 、旋 风分离器直径 D和柱段长径比 L ／ D 以及锥段长径 比 ￡ 2 等有关。加大 D、三 1 ／ D或 2 ，或者减小 都能延长停留时间，使更小的液滴完成沉降，从 而提高分离效率。但是，加大 D 或减小 Q i 都会显 著地降低离心加速度， 使得液滴沉降的驱动力不够。 因此，只有增加旋风分离器的长度才是增加停留时 间的有效手段 。 通过对不 同柱段和锥段长度 的旋风分离器 的 D P M 两相模拟， 增加旋风分离器的长度可有效提高 旋风分离器的分离效率。 保持锥段长度 2 不变，图 1 0为不同柱段长径 比 1 的旋风分离器的级效率曲线。 旋风分离器的 级效率 曲线随着柱段长径 比 ／ D的增大向左偏移， 5 o减小，分离效率提高。 保持柱段长度 1 不变，随着锥段长度 2 的增 加， 如图 l 1 ，级效率曲线向左偏移， X 5 0 减小， 分离 效率增大。 在实际工程 中，更长的旋风分离器的设计和制 造成本更高，需要综合考虑合适的长径 比。 图 1 0 不同柱段长度 L 1 的级效率曲线 Q i 5 o o m3 h 图 1 1 不同锥段长度L 2 的级效率曲线 Q i 5 0 0 m h 3 结 论 I 使用 R S M 湍流模型模拟天然气脱蜡旋风 分离器，所得的流场分布规律符合一般旋风分离器 的流场分布规律 旋风分离器内部压力呈对称分布， 由壁面向中心逐渐减小，在 中心轴处压力最低并有 负压产生；切向速度符合准 自由涡分布规律，呈现 “ M” 形驼峰分布；轴向速度分布反映了向下的外旋 流和 向上的内旋流的流动规律。 2 使用 DP M 模型模拟旋风分离器 内的两相 流动 ，随着进 口流量的增加，旋风分离器 的级效率 曲线 向左移动，奶。 减小，分离效率提高，与理论计 算相符。当入 口流量为 1 0 0 0 m3 ／ } l 时，X 5 0 可以达到 5 ． 3 g m。大粒径液滴的分离效果明显，但在所研 究 的处理量范围内，处理量 的变化不能明显地影响小 于 5 la m液滴的分离效率。 3 柱段和锥段长度 的增加增长 了旋风分离 器的的总长度 ，延长了液滴在旋风分离器内的停留 时间， 旋风分离器的级效率曲线向左偏移， X 。 减小 ， 分离效率增大 。 学兔兔 w w w .x u e t u t u .c o m 第 1 期 袁惠新等天然气脱蜡旋风分离器分离效率的模拟 4 9 参考文献 冯异勇．石蜡沉积及清除数值模拟研究[ D】 ．四川 西南石油学院， 2 0 0 2 ． 袁惠新．天然气净化装置中国，1 0 2 7 8 6 9 9 5 A 【 P 1 ． 2 0 1 2 ． 0 8 ． 3 1 ． H s i e h K T，R a j a m R K _ Ma th e ma fi a c mo d e l o f t h e h y d r o c y c l o n e b a s e d o n p h y s ic s o f fl u i d fl o w n． A I C h E J o u r n a l ， 1 9 9 1 ， 3 7 5 7 3 5 7 4 6 ． 刘子红，肖波， 杨家宽．旋风分离器两相流研究综述『 J 】 ． 过滤与分 离，2 0 0 3 ，1 3 1 4 2 ． 4 5 ． 王海刚，刘石． 不同湍流模型在旋风分离器三维数值模拟中的应 用和 比较f J ] ． 热能动力工程，2 0 0 3 ，1 8 4 3 3 7 3 4 2 ． 葛坡 ，袁惠新，付双成． 对称多入口型旋风分离器的数值模拟[ J ] ． 化工进展，2 0 1 2 ，3 1 2 2 9 6 ． 2 9 9 ． 牟春宇． 基于 C F D 的旋风除尘器分离效率的数值模型研究【 D ] ． 无锡江南大学，2 0 0 8 ． 许妍霞． 水力旋流分离过程数值模拟与分析【 D ] ．上海华东理工 大学，2 0 1 2 ． 刘淑艳，张雅 ，王保国． 用 R S M 模拟旋风分离器内的三维湍流流 场[ J ] ．北京理工大学学报，2 0 0 5 ，2 5 5 3 7 7 ． 3 7 9 ． 张宪瑞 ，闰鹏 ，毕荣山． 气液喷射器的不同液滴粒径分布特性研 究[ J 1 ． 科学技术与工程，2 0 1 2 ，1 2 6 3 9 3 3 ． 3 9 3 6 ． 魏耀东，宋健斐，陈建义，等．旋风分离器内气相流场的相似模 化分析 I 流动参数『 J 】 ． 化工学报，2 0 1 0 ，6 1 9 2 2 6 5 ． 2 2 7 3 ． 宋健斐，杨光福，陈建义，等．旋风分离器内气相流场的相似模 化分析 I I 尺寸参数f J ] ． 化工学报，2 0 1 0 ，6 1 9 2 2 7 4 ． 2 2 7 9 ． 袁惠新，冯晶．分离工程【 M】 ．北京中国石化出版社，2 0 0 1 ． 痧 痧 痧 、 痧 ， 、 上接第 1 5页 [ 3 2 ] 戴海风． C O 2 水合物浆流动传热特性的研究 D] ．上海 上海理工大 R e n e w a b l e a n d S u s t a i n a b l e E n e r g yR e v i ews ，2 0 1 3 ，2 2 1 0 8 1 2 0 ． 学，2 0 l 3 ． [ 4 l 】 马志伟，张鹏，王如竹．水合物浆体在直管中的流动换热特性研 【 3 3 】 S a r i O，H u J ，E i c h e r S ，e t a 1 ． I n t e r n a t i o n a l R e f r i g e r a t i o n a n d A i r 究[ J ] ．工程热物理学报，2 0 1 0 ，3 1 8 1 3 9 8 ． 1 4 0 2 ． C o n d i t i o n i n g C o n f e r e n c e [ C ] ． P u r d u e P u r d u e Un i v e r s i t y ， J u l y 1 4 - 1 7 ， 【 4 2 1 S t a ma t i o u E，Me e w i s s e J W ，K a wa j i M ．I c e s l u r r y g e n e r a t i o n 2 0 0 8 i n v o 1 v i n g mo v i n gp a r t s [ J ] ．／ n t． R e f r i g ． ，2 0 0 5 ，2 8 1 6 0 ． 7 2 ． [ 3 4 】 K i t a n o v s k i A，P o r e d o s A C o n c e n t r a t i o n d i s t r i b u t i o n a n d v i s c o s i ty o f 【 4 3 ] Cl a r k e M A ，Bi s h n o i P R． De t e r mi n a t i 0 n o f t h e i n t r i n s i c r a t e c o n s t a n t i c e 。 s l u r r y i n h e t e r o g e n e o u s fl o w [ J ] ． ／ n t ． R ef r ig ． ，2 0 0 2 ，2 5 6 and a c t i v a t i o n e n e r g y o f C O 2 g a s h y d r a t e d e c o m p o s i t i 0 n u s i n g i n - s i t u 8 2 7 - 8 3 5 l p a r t i c l e s i z e ana l y s i s [ J ] ． ； h e m． En g． S c i ． ， 2 0 0 4，5 9 1 4 2 98 3 2 9 9 3 ． 【 3 5 ] O y a ma H，E b i n u ma T，S h i m a d a W，e t a 1 ． 4 t h I n t ． C o n e o n G a s [ 4 4 1 C l a r k e M A，B i s h n 0 i P R ．De t e r mi n a t i o n o f t h e i n t r i n s i c k i n e t i c s o f H y d r a t e s [ C ] ． Y o k o h a ma ，J a p a n ，2 0 0 2 ． C O2 g a s h y dra t e f o r ma t io n u s i n g i n s i tu p a r t i c l e s i z e ana l y s is f J ] ． [ 3 6 ] 肖睿，宋文吉，黄冲，等． T B AB包络化合物浆的管内流动再层流 C h e m． E n g ． S c i ． ，2 0 0 5 ，6 0 3 6 9 5 7 0 9 ． 化现象[ J ] ．工程热物理学报 ，2 0 0 9 ，3 0 6 9 7 1 - 9 7 3 ， 『 4 5 1 My r i a m D，Mi c h e l C，J e an． Mi c h e 1 H． R h e o l o g i c a l s t u d y o f T B AB [ 3 7 】 Ma z w，Z h a n g P ，Wang R Z ，e t a 1 ． F o r c e d fl o w a n d c o n v e c t i v e h y d r a t e s l u r r i e s a s s e c o n d a r y t w o - p h a s e r e f r i g e r ant s [ J ] ． ／ n t ． 辔． ， me l t i n g h e a t tr a n s f e r o f c l a t h r a t e h y d r a t e s l u r r y i n t u b e s [ J ] ． ／ n t ． o r ． 2 0 0 5 ，2 8 5 6 6 3 6 7 1 ． He a t Ma s s T r a n sf． ，2 0 1 0 ，5 3 1 9 2 0 3 7 4 5 - 3 7 5 7 ． 『 4 6 ] De l a h a y e A，F o u ma i s o n L，Gu i l p a r t J ． C h ar a c t e r i z a t io n o f i c e a n d 【 3 8 ] 宋文吉，肖睿 ，董凯军，等． T B A B包络化合物浆水平管内的传热 1 1 { F h y dra t e s l u r r y c r y s t a l s i z e d i s t r i b u t i o n b y mi c r o s c o p i c 特性研 究『 J ] ． 工程热物理学报，2 0 0 9 ，3 0 9 1 5 0 2 1 5 0 4 ． o b s e r v a t i o n me t h o d [ J ] ． ／ n t ． R ef r ig ． ，2 0 1 0 ，3 3 8 1 6 3 9 1 6 4 7 ． 【 3 9 】 Ma Z W，Z h ang P P r e s s u r e d r 0 p and h e a t tr ans f e r c h a r a c t e r i s t i c s o f 【 4 7 ] I n a d a T，Mo d a k P R ． Gr o w t h c o n tr o l o f i c e c r y s t a l s b y p o l y v i n y l c l a t h r a t e h y dra t e s l u r r y i n a p l a t e h e a t e x c h a n g e r [ J ] ． ／ n t ． R ef r i g ． ， a l c o h o 1 and ant i fre e z e p r o t e i n i n i c e s l u r r i e s [ J ] ． C h e m ． E n g ． S c i ． ， 2 0 1 1 ，3 4 3 7 9 6 8 0 6 ． 2 0 0 6 ， 6 1 1 O 3 1 4 9 31 5 8 ． [ 4 0 】 Z h a i XQ，Wa n gxL ，WangT ，e t a 1 ． A r e v i e wo np h ase c h ang e c o l d [ 4 8 】 H u o Z ，F r e e r E，L a mar M，e t a 1 ． Hy dra t e p l u g p r e v e n t i o n b y s t o r a g e i n a i r - c o n d i t i o n i n g s y s t e mMa t e r i a l s a n d a p p l i c a t i o n s [ J ] ． ant i - a g g l o me r a t i o n [ J ] ． C h e m． Eng ． S c L ，2 0 0 1 ，5 6 1 7 4 9 7 9 - 4 9 9 1 ． 。 ． 、 移、 、 ， 、 、 移 、 I 上接第 2 0页 [ 1 3 】 [ 1 4 】 [ 1 5 】 [ 1 6 ] [ 1 7 】 【 l 8 】 [ 1 9 1 [ 2 O 】 刘新 民，王兵，荆海洋，等． 用 R1 4 2 b裂解制备偏氟乙烯的方法 和设各中国，1 0 1 0 0 3 4 6 0 『 P 】 ． 2 0 0 7 ． 0 7 2 5 ． F r a n c i s H ， Wa l k e r At t i l a E ． De h y d r o h a lo g e n a t i o n o f l ， l ， 1 一 tr i h a l o e t h ane s [ J ] ．J o u r n a l o f O r g a n i c C h e m i s t r y， 1 9 6 5，3 0 3 2 8 4 3 2 8 5 ． M u l l e r H， Erni g G ， Ho f ma n n H． S t u d i e s o n h e t e r o g en e o u s ly c a t a l y z e d d e h y d r o c h l o r i n a fi o n o f 1 ． 1 - d i f l u o rn- l - c h l o r o e th ane [ J ] ．C h e m I n g ． T e c h ． ， l 9 8 4 ， 5 6 8 6 2 6 ． 6 2 8 ． E l s h e l k h ． Ca t a l y t i c p r o c e s s f o r t h e d e h y d r o h a lo g e n a t io n o f 1 - c h l o r o - 1 ， 1 ． d i fl u o r o e t h ane E P ，0 4 0 7 7 1 1 f P ] ． 1 9 9 1 0 1 - 1 6 ． M a h e r M i c h a e 1 ．P r o d u c t i o n o f v iny l i d e n e fl u o r i d e b y r e a c t i o n o f 1 - c h l o r o 1 ． 1 - d i fl u o r o e t h ane a n d h y d r o g e n c h l o r i d e i n g a s e o u s p h a s e a n di nth ep r e s e n c e o f a c a ta l y s t F R ，2 6 9 0 6 8 7 [ P ] ． 1 9 9 3 1 l O 5 ． W a c k e r C G V e r f a h r e n z u r d e h y dr o h a l o g e n i e r u n g v o n a s y mme t r i s c h s u b s t i t u i e r t e n p o l y h a l o g e n a e t h ane n DE ，1 9 2 5 5 6 8 [ P 】 ． 1 9 7 0 - 1 1 - 2 6 ． R u s s e l l M，Wi l l j ．姗 S ． Se l e c t i v e d e h y dro h a l o g e n a t i o n o f fl u o roh a l o a l k ane s u s in g a c o p p e r c a t a l y s t US， 2 7 7 4 7 9 9 [ P ] ． 1 9 5 6 ． 1 】 一 l 8 ． 王建 明．催化裂解生产低碳烯烃技术和工业应用进展[ J 1 ．化工进 展 ，2 O1 l ，3 0 5 9 1 1 - 91 7 ． [ 2 l 】 郭天祥．催化裂解多产低碳烯烃催化剂的研究[ D ] ． 北京北京化 工大 学 ，2 0 1 1 ． [ 2 2 1 S c h e r e r O，S t e i n me t z A， Ku h n H，e t a 1 ． P r o c e s s f o r t h e m a n u f a c t u r e o f flu o r e t h y l e n e s DE，1 0 6 8 6 9 5 『 P 1 ． 1 9 5 8 1 1 1 2 ． [ 2 3 1 C l a u s L ，G u e n t h e r E ． V e r f a h r e n z u r h e r s t e l l u n g y o n fl u o r a l k e n e m DE， 1 9 0 6 2 0 6 『 P 1 ． 1 9 7 0 1 0 0 8 ． [ 2 4 】 B a u e r W V，S z e M C ．S t e am p y r o l y s i s o f h y d r o c a r b o n s U S ， 4 3 4 2 6 4 2 f P 1 ． 1 9 8 2 0 8 0 3 ． [ 2 5 ] 蔡宗铉，李源镐，丁相文，等．裂解管和使用裂解管的裂解方法 中国，1 4 5 7 3 5 6 『 P 1 ． 2 0 0 3 - 1 1 - 1 9 ． [ 2 6 ] 王军，王汉利，郝宏伟，等． 一种高纯偏氟乙烯单体的生产方法 中国 ，1 O l 7 0 4 7 0 9 [ P ] ． 2 0 1 0 0 5 ． 1 2 ． [ 2 7 】 周兴贵，粱聪强，吴君毅．一种偏氟乙烯单体的制备方法中国， 1 0 2 8 3 8 4 4 7 [ P ] ． 2 0 1 2 ． 1 2 2 6 ． [ 2 8 1 郑根土． 偏氟乙烯生产过程仿真研究和优化[ J ] ．化工生产与技术， 2 0 1 l ， l 8 5 l 】 一 l 4． [ 2 9 】 唐杰．偏氟乙烯分离工艺的过程模拟分析[ 】 】 ． 上海 华东理工大 学 ，2 0 1 2 ． Ⅲ Ⅲ 嘲 嘲 学兔兔 w w w .x u e t u t u .c o m