2016-2017年换热管破裂计算公式分析

2016-2017年换热管破裂计算公式分析 安全阀排放量计算公式分析 在安全阀的排放量计算方法里，在目前能见到的计算方法中，我认为有几个计算公式还有讨论的必要，因为就同一个工况而言，不同的公式相互矛盾，我先以换热管破裂的排放量为例，讨论如下: 一 换热管破裂的安全阀计算公式. 原来使用的公式1) SEI 420.5 ××?×ρ×W=246.310d(p) kg/h Vv 420.5××?ρ×W=16.810d(p/) kg/h LL 3,内径?高低压侧压差ρ，ρ气、液相重度 m p Mpa kg/m vL 日本三菱重工的公式2) 换热器低压侧的设计压力大于/的高压侧的设计压力时，由于整根1) 23 管子完全断裂极少发生，可不考虑换热器管子破裂，不需设安全阀; 换热器低压侧设计压力小于或等于/高压侧的操作压力时，整根管2) 23 子完全断裂虽很少发生，但要考虑换热器管子破裂，需设安全阀; 换热器高压侧的操作压力大于等于(),约(),，3) 7MPa.G1000psi.G或低压侧的介质是能闪蒸液体，或介质是含有蒸气、会汽化的液体;整根管子完全断裂的故障应考虑，而且与差压无关，需设安全阀。 换热管破裂引起的安全阀排放量，按下面式计算:3 0.5高压侧是气体的公式:××× W=2345a(P-P) (M/T) kg/h HVC 0.5高压侧是液体的公式:××()× W=9963a[P-PS] kg/h HVC 高压侧是闪蒸液体的公式: 0.5当(),×时:××P-1.1PP-P W=9963a[(P-P)] HDL1VCHVC 0.5其中)×(× P=[0.96-0.28P/P] P VCVCV 二根换热管的面积气体泄放温度W kg/h a cm2 T K 气体分子量高压侧的操作压力，低压侧的设计压力M P P MPaG HDL 高压侧的操作压力，破裂处的脉冲压力P P MPaA 1VC 高压侧介质在操作温度下的蒸汽压P MPaA V 1 高压侧介质的临界压力高压侧液体的比重P MPaA S C 在决定安全阀的尺寸时，要考虑过程热量输入的潜在能力，不能只考虑正常的热量输入。例如，炉子燃烧器的最大负荷往往可达炉子铭牌负荷的125 %，此时应按泄压工况下的最大蒸汽发生量减去正常冷凝量或蒸汽流出量选用安全阀; 无火灾危险处的非易燃液化气，由于阳光暴晒、或其它原因造成辐射、对流、传导等的热传递，使容器内贮存的物质被加热，迅速汽化超压。根据压力容器安全技术监察规程，这类容器的安全阀的安全排量不低于火灾危险处的非易燃液化气的30 %; 管子或容器在充满冷液体后被关闭，然后又受到伴热管、蛇管、外界 热流或明火加热; 一个换热器的冷侧被切断，而热侧仍有流体通过; 管子或容器在接近常温时充满液体后，然后又直接受到太阳辐射加 热; 对于用蒸汽加热的再沸器和类似的管式换热器，在决定调节阀故障时的换热工况时，假设管子是清洁无污垢的; 从“工艺管线设计安装手册”下册中查到的计算公式:3) 420.5对气体××××ρ G=10710d(P) kg/h v 420.5对液体×××?γ G=4.710d(p/) m3/h 高压侧压力(表)?高低压侧压差ρ气体密度P atm p atm v kg/m3 γ 液体比重 换热管内径 d m 对气体排放而言，要求高压侧的压力至少是低压侧的两倍，如果低于两倍要适当降低排放量。 兰化院编的书中的公式:4) 20.5对液体×××? W=5.6d(Gp) kg/h L 液体密度?高低压侧压差换热管的内径G kg/m3 p Mpa d mm L 公司的计算公式:5) Kvaerner 20.5 对气体××××ρ W=1.5670.7d(P ) kg/h VHH 此式的假设是，高压气体排除的速度是音速的倍。0.7 2 3高压侧的最大操作压力ρ高压流体的密度P bara kg/m HH 换热管的内径d mm 在以上公式中，它们的计算基准都是以两根换热管的面积来计算排放量，符合API的论述。第二组公式考虑了闪蒸的流体，其它公式都没有考虑。故先讨论其它四组公式，分为气体和液体的公式分别讨论。 二. 对气体的计算公式 有三组公式: 420.5 ××?×ρ×W=246.310d(p) kg/h Vv ,内径?高低压侧压差ρ，ρ气、液相重度 m p Mpa v kg/m3 L 20.5 ××ρ××W=1.5670.7d(P ) kg/h VHH 高压侧的最大操作压力ρ高压流体的密度P barA kg/m3 HH 换热管的内径d mm 420.5××××ρG=10710d(P) kg/h v 高压侧压力(表)ρ气体密度换热管内径P atm v kg/m3 d m 420.5第二式经过单位换算，可换算为:×××ρ与第×W=10910d(P )Vv三式基本一致～ 我也试验了前面介绍的第二组公式(日本三菱重工的公式)，结果发现它的计算结果比其它公式的结果都要大很多。 对第一式我分析它与第二和第三式的区别，是用压差代替高压侧的压力计算，而压差是一个变化的数，不是常数，所以此式的常数项，应该是个变数，不是常数。 三. 对液体的计算公式 也有下面三个公式: 420.5××?ρ×W=16.810d(p/) kg/h LL ,内径?高低压侧压差ρ液相重度 m p Mpa kg/m3 L 420.53×××?γG=4.710d(p/) m/h ?高低压侧压差γ 液体比重 换热管内径p atm d m 20.5××?×W=5.6d(Gp) kg/h LL 3 3液体密度?高低压侧压差换热管的内径GL kg/m p Mpa d mm 首先对公式进行单位分析，第一和第三式的单位也有问题，第一式的单位应 3该是×?的结果不会是，才对,第三式的形式就有问题，用m/hGpkg/hL 我觉得这个公式本身就有问题，故建议不采用第三式～ 在第一式和第二式之间的差别，仅仅是ρ与γ，和的差别，经MpaatmL 过换算可知它们之间的差别是，×不是。0.031416.80.0314=0.5284.7 用一个例子来核算，高低压差是或，换热管径是0.5 Mpa5 atm ，介质是水，ρ，γ。0.025mm=1000=1 L 用第一式计算结果是;用第二式计算结果是差别很大～ 2.348 65.68说明这两个公式不一致，它们的系数有问题～ 四. 研究API推荐的公式 经过长时间的调查，我们总算找到了上面列举的全部公式的来源，即找到了它们的“老祖”～那就是英国CRANE公司编写的“FLOW OF FLUIDS”一书。在这本书中列出了有关流体流动的全部公式，分为液体、可压缩气体来论述，和换热管破裂有关的公式，是通过孔板的计算公式，我把它们分别列出在后面: 对于液体的计算公式: 20.5×××?ρQ=236Cd(p/) gal/min 20.5 w=0.525×××?×ρCd(p)lb/sec 20.5 ×××?×ρW=1891Cd(p)lb/h 23?流过孔板的压差，ρ 液体比重 p lb/in lb/ft 换热管内径，对流体的流量系数d in C 对于可压缩气体的计算公式是: 20.5 w=0.525××××?CdY(p/V)lb/sec 20.5 ××××?W=1891CdY(p/V)lb/h 23?流过孔板的压差， 气体的比容 p lb/in V ft/ lb 换热管内径，对流体的流量系数d in C Y 可压缩流体经过孔板后的膨胀系数 对比这些公式和上面介绍的公式，我们可以发现最大的差别是，前面介绍过的全部公式，都是以高低压端的压差为基准的，而公式的来源是以经过孔板的压 4 力损失来计算的～改为以高低压端的压差或就以高压端的压力为基准计算，是方便多了，可是需要合理的推导才能实现。 根据美国《Hydrocarbon processing February 1992》载文“PRV sizing for exchanger tube rupture” 作者是 Wing.Y.Wong，他对换热管破裂工况的安全阀排放量计算有严格的介绍，把他的介绍整理如下。 对于换热管破裂工况的计算可按以下步骤进行: 1)先依据换热器管侧及壳侧的高低压的设计压力，判断需不需要计算 换热管破裂的工况，依据是两端的设计压力之比是否小于2/3; 2)用 流程 快递问题件怎么处理流程河南自建厂房流程下载关于规范招聘需求审批流程制作流程表下载邮件下载流程设计 模拟软件计算破裂后的混合相的相态，得出两相的流量和密 度; 3)依据API520用公式计算破裂后的临界压力， K2 K,1，，P,P，CFR1K，1 P------高压侧的操作压力 1 K------由模拟计算得到，此例K=1.1 P---临界压力 CFR 4)计算所需的理论排放量 如果通过破裂口流体没有相变，对可压缩流体而言，应用下式计算气 20.5 体流量:W=1891×Y×d×C×(dP/V)V 此时Y=1-0.317dP/P1 Y值的计算式是从图中简化计算出来的。 此时C=0.6 这个流量系数是从锐边孔板C值的图形中查出的。 2d=4×A/3.1416 V 0.5 整理后 W=1444.6A(1-0.317dP/P) (dP×LO)VV1V 通过破裂口流出的是液体时的计算公式: 20.5 W=1891 dC (dP×LO)LL 对于此例C=0.6 2d=4×A/3.1416 L 0.5 整理后得:W=1444.6A (dP/LO)LVL 5)计算两相流中的汽相比例 5 R=W/(W+W) VVL 6)计算所需的实际排放量 举例说明: 见附图 低压端的设计压力:150psig 高压端的操作压力:330psig=330+14.7=344.7psia 高压端的设计压力:330×110%=363psig 150/363=0.41,2/3所以应该计算换热管破裂的工况， 计算高压流体通过破裂的换热管时的临界压力: 1.1K221.1,1=344.7×()=201.54psia K,1，，P,P，CFR1K，11.1，1 计算高低压混合后的流体通过安全阀喷嘴时的临界压力: 1.1 21.1,1P=P×()=179.7×0.585=105.07psiaCFOREL1.1，1 P是安全阀的排放压力=150×110%+14.7=179.7psia REL 通过流程模拟计算，可以得到以下数据: 3 W=266lb/h LO=1.91lb/ftVCVC 3 W=734lb/h LO=29.88lb/ftLCLC R=W/(W+W)=266/(266+734)=26.6% VVL dP 0.5W=1444.6×A×(1-0.317)(dP/V) VVP1 dP=P-Pcf=344.7-201.54=143.16psi 1 dP0.50.5W=1444.6A(1-0.317)(143.16)(1.91) VVP1 =20742.8 A V 0.50.5W=1444.6A (dP/LO)=1444.6 A (143.16×1.91) LVLV =94481.9 AL 查换热器的管径是0.532in.ID 22 A=A+A=2×(3.1416 d/4)=0.4446inVL 可以解以上方程得出以下结果: 3W=5744lb/h Q=5744/1.91=3007ft/h VV 3W=15847lb/h Q=15847/29.88=530 ft/h LL 6 Y=1-0.317×dP/P1 PRV 火炬 低压侧 上水 回水 高压侧 低压侧设计压力:150 psig 高压侧操作压力:330 psig 流向 高压侧 低压侧 背压 P=201.54 psia P=105.07 psia CFRCFO (tube rupture) (PRV offive) P=344.7 psia P=179.7 psia P=19.7 psia OPRELBACK 换热器管壁 安全阀喷嘴 7 气体、液体、两相共同涉及的参数: K--绝热指数，模拟计算值 d--换热管内径，d=换热管外径-2×换热管壁厚，mm C,Y值的计算: 把换热器的破裂口看成什么类型的孔板，过去没有认真研究，就和自控专家确定按偏心锐孔板的公式计算，最近在蔡尔辅写的第423页上看到，有明确规定按同心锐孔板计算，所以要修改这些公式。 Y = 1 - (0.1926+0.574×B + 0.9675×B2 - 4.24×B3 + 3.62×B4)×(?P / P1)×(1/K) 改为: 234Y=1-(0.333+0.574B+0.9675B-4.24B+3.62B) ×(?P / P1)×(1/K) B--孔径比 B = 换热管内径/换热器壳体外径 ?P = (P高压侧设计压力 – P低压侧设计压力)/P×[2/(K+1)]k/(k-1) 1 P—高压侧设计压力，MPa.A 1 K--绝热指数，模拟计算值 如果: d <= 100mm C = 0.5875 + 0.3813×B^2.1 + 0.6898×B^8 - 0.1963×B^4 / (1 – B^4) - 0.3366×B^3 + (7.3 - 15.7×B + 170.8×B^2 - 399.7×B^3 + 332.2×B^4) / Re^0.75 如果: d > 100mm C = 0.5949 + 0.4078×B^2.1 + 0.0547×B^8 - 0.0955×B^4 / (1 – B^4) - 0.5608×B^3 + (-139.7 + 1328.8×B - 4228.2×B^2 + 5691.9×B^3 - 2710.4×B^4) / Re^0.75 同心锐孔板计算公式: 如果: d <= 58.62mm C=0.5959+0.0312×B^2.1-0.184×B^8+0.0029×B^2.5×(10^6/Re)^0.75-0.85598×B^3/ D+0.039×B^4/(1-B^4) 如果: d > 58.62mm C=0.5959+0.0312×B^2.1-0.184×B^8+0.0029×B^2.5×(10^6/Re)^0.75-0.85598×B^3/ D+2.286×B^4/((D×(1-B^4) B--孔径比，B = 换热管内径/换热器壳体外径 Re = d×ρ× u /μ d--内径, m ρ--换热管内介质密度, kg/m3 换热管内介质流速, m/s u-- μ--换热管内介质粘度, mPa.s 分析换热管破裂时的公式 1.对于气体 42CRANE公式: W=400×10×Y×C×d×?(?P×ρ) V 8 42原BDI的公式: W=246.3×10×d×?(?P×ρ) V 这个公式没有Y，C值，就是把问题简化了，把单管直径和换热器壳体直径的比视为趋向于零，则 Y=1-0.1926×?P/P×1/k 1 因为有2/3原则，P,2/3，?P,1/3P，若?P= P，则 低高高 Y=1-0.1926×1/1.3×1=0.852 0.750.75C=0.5875+7.3/Re，可视 7.3/Re趋向于零。 所以C=0.5875+0=0.5875，而 Y×C=0.852×0.5875=0.5 400×0.5=200，近似246，但问题是此式没有计算两根管的面积，也不符合API风格。 9 10