发明名称 一种加氢反应流出物系统失效控制参数的建模方法
摘要 本发明公开了一种加氢反应流出物系统失效控制参数的建模方法。失效控制参数建模包括Kp值、NH4HS浓度、理论注水量、NH4HS沉积温度、NH4Cl沉积温度、空冷器平均流速等共六个部分。通过DCS控制系统数据库读取加氢反应流出物系统运行参数,结合原料化验分析数据、加氢反应原则流程图,针对具体加氢反应流出物系统变工况运行中的失效控制参数进行建模,确定Kp值、NH4HS浓度、理论注水量、NH4HS沉积温度、NH4Cl沉积温度、空冷器平均流速,可供现场操作人员对加氢反应流出物系统失效控制参数进行有效地监控,避免加氢反应流出物系统失效引发的非计划停工事故,确保加氢反应流出物系统的长周期、安全、稳定运行。
申请公布号 CN101225326B 申请公布日期 2011.03.23
申请号 CN200710160373.5 申请日期 2007.12.18
申请人 浙江理工大学 发明人 偶国富;郑智剑;王艳萍;曹海彬;金浩哲
分类号 C10G45/72(2006.01)I;C10G47/36(2006.01)I;C10G49/26(2006.01)I;G05B19/418(2006.01)I 主分类号 C10G45/72(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 林怀禹
主权项 1.一种加氢反应流出物系统失效控制参数的建模方法,其特征在于:失效控制参数包括六个部分:Kp值、NH<sub>4</sub>HS浓度、理论注水量、NH<sub>4</sub>HS沉积温度、NH<sub>4</sub>Cl沉积温度、空冷器平均流速;其中:1)Kp值:首先,通过DCS控制系统数据库读取加氢反应流出物系统原料油进料量实际值W<sub>f</sub>、循环氢流量mX<sub>1</sub>和新氢流量mX<sub>2</sub>,结合原料油进料量设计值<img file="FSB00000342690900011.GIF" wi="88" he="56" />空冷器入口混合物流量设计值mH<sub>1</sub>,确定空冷器入口混合物流量实际值mH:<maths num="0001"><![CDATA[<math><mrow><mi>mH</mi><mo>=</mo><msub><mi>mH</mi><mn>1</mn></msub><mo>&times;</mo><mfrac><msub><mi>W</mi><mi>f</mi></msub><msub><mi>W</mi><msub><mi>f</mi><mn>1</mn></msub></msub></mfrac></mrow></math>]]></maths>其次,输入加氢反应流出物系统原料化验分析数据,结合原料油进料量实际值W<sub>f</sub>和循环氢流量mX<sub>1</sub>分别确定H<sub>2</sub>S和NH<sub>3</sub>的摩尔流量:<maths num="0002"><![CDATA[<math><mrow><msub><mi>mNH</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>W</mi><mi>f</mi></msub><mo>&times;</mo><msub><mi>F</mi><mi>n</mi></msub></mrow><mrow><mo>(</mo><mi>mwN</mi><mo>)</mo></mrow></mfrac></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>mH</mi><mn>2</mn></msub><mi>S</mi><mo>=</mo><mfrac><mrow><msub><mi>W</mi><mi>f</mi></msub><mo>&times;</mo><msub><mi>F</mi><mi>s</mi></msub></mrow><mrow><mo>(</mo><mi>mwS</mi><mo>)</mo></mrow></mfrac><mo>+</mo><msub><mi>mX</mi><mn>1</mn></msub><mo>&times;</mo><msubsup><mi>F</mi><mi>s</mi><mo>&prime;</mo></msubsup></mrow></math>]]></maths>式中:F<sub>n</sub>为原料中N的含量;F<sub>s</sub>为原料中S的含量;F′<sub>s</sub>为脱硫前循环氢中的H<sub>2</sub>S含量;最后,分析空冷器入口混合物流量实际值mH、H<sub>2</sub>S和NH<sub>3</sub>的摩尔流量确定加氢反应流出物系统的Kp值:<maths num="0004"><![CDATA[<math><mrow><mi>Kp</mi><mo>=</mo><mrow><mo>(</mo><msub><mi>H</mi><mn>2</mn></msub><mi>S</mi><mo>)</mo></mrow><mi>mol</mi><mo>%</mo><mo>&times;</mo><mrow><mo>(</mo><msub><mi>NH</mi><mn>3</mn></msub><mo>)</mo></mrow><mi>mol</mi><mo>%</mo><mo>=</mo><mfrac><mrow><msub><mi>mH</mi><mn>2</mn></msub><mi>S</mi></mrow><mi>mH</mi></mfrac><mo>&times;</mo><mn>100</mn><mo>&times;</mo><mfrac><msub><mi>mNH</mi><mn>3</mn></msub><mi>mH</mi></mfrac><mo>&times;</mo><mn>100</mn></mrow></math>]]></maths>2)NH<sub>4</sub>HS浓度:通过DCS控制系统数据库读取加氢反应流出物系统的实际注水量WW′<sub>r</sub>,结合加氢反应流出物系统原料化验分析数据,确定NH<sub>4</sub>HS浓度:<maths num="0005"><![CDATA[<math><mrow><msub><mi>Wt</mi><mrow><msub><mi>NH</mi><mn>4</mn></msub><mi>HS</mi></mrow></msub><mo>%</mo><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><msub><mi>mwNH</mi><mn>4</mn></msub><mi>HS</mi><mo>)</mo></mrow><mo>&times;</mo><msub><mi>W</mi><mi>f</mi></msub><mo>&times;</mo><msub><mi>F</mi><mi>n</mi></msub></mrow><mrow><mrow><mo>(</mo><msub><mi>mwNH</mi><mn>4</mn></msub><mi>HS</mi><mo>)</mo></mrow><mo>&times;</mo><msub><mi>W</mi><mi>f</mi></msub><mo>&times;</mo><msub><mi>F</mi><mi>n</mi></msub><mo>+</mo><msub><msup><mi>WW</mi><mo>&prime;</mo></msup><mi>r</mi></msub><mo>&times;</mo><mrow><mo>(</mo><mi>mwN</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>3)理论注水量:通过DCS控制系统数据库读取加氢反应流出物系统注水点处的平均温度T,查饱和水蒸汽的绝对压力表通过插值可得温度T时饱和水蒸汽的绝对压力P<sub>satstm</sub>;其次,根据注水点处工艺介质的相分率Z<sub>c</sub>和Kp值建模过程中空冷器入口混合物流量实际值mH,确定注水点处气相摩尔流量H<sub>c</sub>:H<sub>c</sub>=z<sub>c</sub>×mH最后,结合温度T时饱和水蒸汽的绝对压力P<sub>satstm</sub>、注水点处工艺压力P<sub>system</sub>及工艺压力系数F<sub>c</sub>,确定理论注水量WW<sub>r</sub>:<maths num="0006"><![CDATA[<math><mrow><msub><mi>WW</mi><mi>r</mi></msub><mo>=</mo><msub><mi>F</mi><mi>c</mi></msub><mo>&times;</mo><msub><mi>H</mi><mi>c</mi></msub><mo>&times;</mo><mfrac><mrow><msub><mi>P</mi><mi>satstm</mi></msub><mo>/</mo><msub><mi>P</mi><mi>system</mi></msub></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mi>satstm</mi></msub><mo>/</mo><msub><mi>P</mi><mi>system</mi></msub><mo>)</mo></mrow></mfrac><mo>&times;</mo><mn>1.25</mn><mo>&times;</mo><mrow><mo>(</mo><msub><mi>mwH</mi><mn>2</mn></msub><mi>O</mi><mo>)</mo></mrow></mrow></math>]]></maths>4)NH<sub>4</sub>HS沉积温度:根据理论注水量建模过程中计算得出的注水点处气相摩尔流量H<sub>c</sub>、注水点处工艺压力P<sub>system</sub>、Kp值建模过程中计算得出的H<sub>2</sub>S和NH<sub>3</sub>的摩尔流量,确定NH<sub>4</sub>HS沉积温度系数K<sub>1</sub>:<maths num="0007"><![CDATA[<math><mrow><msub><mi>K</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>mNH</mi><mn>3</mn></msub><mo>&times;</mo><msub><mi>mH</mi><mn>2</mn></msub><mi>S</mi></mrow><msup><mrow><mo>(</mo><msub><mi>H</mi><mi>c</mi></msub><mo>&times;</mo><mn>6.89</mn><mo>)</mo></mrow><mn>2</mn></msup></mfrac><mo>&times;</mo><msubsup><mi>P</mi><mi>system</mi><mn>2</mn></msubsup></mrow></math>]]></maths>通过NH<sub>4</sub>HS沉积温度曲线确定NH<sub>4</sub>HS沉积温度;5)NH<sub>4</sub>Cl沉积温度:根据理论注水量建模过程中计算得出的注水点处气相摩尔流量H<sub>c</sub>,结合原料化验分析数据、注水点处工艺压力P<sub>system</sub>、原料油进料量实际值W<sub>f</sub>确定NH<sub>4</sub>Cl沉积温度系数K<sub>2</sub>:<maths num="0008"><![CDATA[<math><mrow><mi>mHCl</mi><mo>=</mo><mfrac><mrow><msub><mi>W</mi><mi>f</mi></msub><mo>&times;</mo><msub><mi>F</mi><mi>e</mi></msub></mrow><mrow><mo>(</mo><mi>mwCl</mi><mo>)</mo></mrow></mfrac></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><msub><mi>K</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>mNH</mi><mn>3</mn></msub><mo>&times;</mo><mi>mHCl</mi></mrow><msup><mrow><mo>(</mo><msub><mi>H</mi><mi>c</mi></msub><mo>&times;</mo><mn>6.89</mn><mo>)</mo></mrow><mn>2</mn></msup></mfrac><mo>&times;</mo><msubsup><mi>P</mi><mi>system</mi><mn>2</mn></msubsup></mrow></math>]]></maths>式中:F<sub>e</sub>为原料中Cl的含量;通过NH<sub>4</sub>Cl沉积温度曲线确定NH<sub>4</sub>Cl沉积温度;6)空冷器平均流速:根据空冷器入口混合物流量实际值mH、原料油进料量实际值W<sub>c</sub>、实际注水量WW′<sub>r</sub>、新氢流量mX<sub>2</sub>、循环氢流量mX<sub>1</sub>及原料化验分析数据,分析确定空冷器平均流速v:<maths num="0010"><![CDATA[<math><mrow><mi>v</mi><mo>=</mo><mfrac><mrow><msub><mi>W</mi><mi>f</mi></msub><mo>+</mo><msup><msub><mi>WW</mi><mi>r</mi></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>mX</mi><mn>1</mn></msub><mo>&times;</mo><mrow><mo>(</mo><msub><mi>mwX</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>mX</mi><mn>2</mn></msub><mo>&times;</mo><mrow><mo>(</mo><msub><mi>mwX</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow><mrow><mi>&rho;</mi><mo>&times;</mo><mi>m</mi><mo>&times;</mo><mi>n</mi><mo>&times;</mo><mi>&pi;</mi><mo>&times;</mo><msup><mrow><mo>(</mo><mi>d</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>&times;</mo><mn>3600</mn></mrow></mfrac></mrow></math>]]></maths>式中:ρ为空冷器入口混合物的密度;m为空冷器的个数;n为每台空冷器管束数量;d为管束内径。
地址 310018 浙江省杭州市江干区经济技术开发区白杨街道2号大街5号