CN109916348A - 一种基于立式储罐本体泄漏口处液体落地水平距离s计算模型 - Google Patents

一种基于立式储罐本体泄漏口处液体落地水平距离s计算模型 Download PDF

Info

Publication number
CN109916348A
CN109916348A CN201910216521.3A CN201910216521A CN109916348A CN 109916348 A CN109916348 A CN 109916348A CN 201910216521 A CN201910216521 A CN 201910216521A CN 109916348 A CN109916348 A CN 109916348A
Authority
CN
China
Prior art keywords
storage tank
liquid
tank
level
leakage
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.)
Granted
Application number
CN201910216521.3A
Other languages
English (en)
Other versions
CN109916348B (zh
Inventor
何娟霞
周琪勇
周冬梅
黄丽文
姚健庭
段青山
李昂昂
李春花
刘磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangxi University
Original Assignee
Guangxi University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangxi University filed Critical Guangxi University
Priority to CN201910216521.3A priority Critical patent/CN109916348B/zh
Publication of CN109916348A publication Critical patent/CN109916348A/zh
Application granted granted Critical
Publication of CN109916348B publication Critical patent/CN109916348B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Filling Or Discharging Of Gas Storage Vessels (AREA)

Abstract

本发明提供一种基于立式储罐本体泄漏口处液体落地水平距离S计算模型，涉及立式液体储罐本体泄漏安全技术领域，包括建立立式液体储罐本体泄漏模型，确定储罐本体任意一泄漏口处，其泄漏口距离水平地面的高度及该泄漏口流出液体落至地面的水平泄漏距离之间的关系；通过已知参数储罐泄漏口距水平地面高度h1，快速计算出任意泄漏口处泄漏液体落地水平距离S，进而推算泄漏口处液体落地最远水平距离Smax，可确定易燃易爆、有毒泄漏液体波及范围，为储罐防火堤设置、易燃易爆有毒危险化学品的探测器安装提供科学依据，同时也可为事故预防提供有效的理论计算参考数据，为事故救援提供准确安全防护距离，科学制定事故应急救援方案、开展有效救援，减少事故损失。

Description

一种基于立式储罐本体泄漏口处液体落地水平距离S计算 模型
技术领域
本发明涉及立式液体储罐泄漏安全技术领域，具体是一种基于立式储罐本体泄漏口处液体落地水平距离S计算模型。
背景技术
企业中储存易燃易爆有毒性质的液体储罐一旦发生泄漏，将会造成严重的后果，易燃易爆有毒液体泄漏后形成的扩散气体，在空间借助风力作用扩散，与空气混合稀释，一旦达到爆炸极限，形成爆炸性混合物，遇明火、高热就极有可能引起燃烧爆炸；且有毒液体泄漏扩散后形成的毒性气体，人体吸入过量后，会影响人的正常生理功能，甚至会造成工作场所人员中毒死亡。因此，具有易燃易爆有毒性质的液体储罐区域都应被重点监控。
目前，迫切需要寻找出一种可预先快速计算出液体储罐泄漏情况的方法，能够为事故预防提供有效的理论计算参考数据，同时可以为储罐防火堤设置、易燃易爆有毒危险化学品的探测器安装等防护问题上提供有力的科学依据，用于各种储存液体储罐泄漏安全问题分析，在工业安全分析方面有广阔的应用前景。
发明内容
基于上述存在问题，本发明的目的在于，提供一种基于立式储罐本体泄漏口处液体落地水平距离S计算模型，本发明通过已知参数储罐泄漏口距水平地面高度h1，快速计算出任意泄漏口处泄漏出液体落至地面的水平距离S，进而推算泄漏口处液体落地最远水平距离Smax，可用于各种储存液体储罐泄漏安全问题分析，确定易燃易爆、有毒泄漏液体波及范围，为储罐防火堤设置、易燃易爆有毒危险化学品的探测器安装提供科学依据，同时也可为事故预防提供有效的理论计算参考数据，为事故救援提供准确安全防护距离，科学制定事故应急救援方案、开展有效救援，减少救援事故损失，在工业安全分析方面有广阔的应用前景。
本发明提供的方案为：
一种基于立式储罐本体泄漏口处液体落地水平距离S计算模型，包括建立立式液体储罐泄漏模型，确定储罐任意一泄漏口处，泄漏口距离水平地面的高度及该泄漏口处流出液体落地水平距离之间关系；通过已知参数储罐泄漏口距水平地面高度h1，快速计算出任意泄漏口处泄漏液体落地水平距离S；其计算关系公式为：
其中，公式中h1为储罐泄漏口距水平地面高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；C0为泄漏系数。
进一步地，由储罐泄漏口距水平地面高度h1计算泄漏口处泄漏液体落地水平距离S模型建立，具体包括如下步骤：
步骤一：建立储罐泄漏口处泄漏出流液体速度v模型
液体泄漏时，储罐内液面下降的静压能和势能，转化为泄漏口液体流出的动能，根据机械能守恒定律，可得到所述储罐内液面的下降速度为：
由质量守恒定律可知储罐内部的泄漏液面下降重量与储罐泄漏口处泄漏出的液体重量相同，即：A1v1＝Av (2)，
结合前式(1)、(2)可得储罐泄漏口处泄漏出流液体速度为：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；h为储罐原始液体高度，m；h1为储罐泄漏口距水平地面的高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s；
步骤二：建立储罐内原始液位h下降到泄漏口位置需要的泄漏时间t关系模型
当液体从储罐泄漏口泄漏后落至地面，此时泄漏口处的液体受重力和储罐内部原始液位高度的压力作用，泄漏液体做类平抛运动；(若不考虑其他影响因素，如风速、风向、恶劣天气状况、孔口拔塞等影响)，根据平抛运动模型，则液体由泄漏口落到地面的时间为：
泄漏口处液体落地水平泄漏距离为：
由公式(4)、(5)可得储罐内原始液位下降到泄漏口位置需要的泄漏时间为：
其中，h1为储罐泄漏口距水平地面的高度，m；v为储罐泄漏口处泄漏出流液体速度，m/s；C0为泄漏系数；A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；h为储罐原始液体高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s；
步骤三：建立储罐内水平液面下降高度Δh与泄漏口处液体落地水平距离S模型
储罐发生泄漏后罐内水平液面下降高度为：
将公式(7)带入到公式(6)，即得到储罐内水平液面下降高度Δh与所述泄漏口处液体落地水平泄漏距离S的关系：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；h1为储罐泄漏口距水平地面的高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s；
根据本发明推导出来的计算模型可计算不同泄漏场景下储罐内液面下降高度Δh与所述泄漏口处液体落地水平距离S的关系；且所述水平距离S与储罐内水平液面下降高度Δh的关系曲线如图1所示；
步骤四：建立储罐泄漏口距水平地面的高度h1与泄漏口处液体落地水平距离S模型
根据求根公式，继续推导公式(8)，且所述泄漏口处液体落地水平距离S的物理意义知其数值为正值，即得到所述泄漏口处液体落地水平距离S与储罐泄漏口距水平地面高度h1的关系为：
其中，h1为储罐泄漏口距水平地面的高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；C0为泄漏系数。
进一步地，所述公式(1)储罐内液面的下降速度v1的模型建立过程包括如下步骤：
步骤一：构建泄漏口上方液体高度hL与泄漏后储罐内液面下降流速v1的模型
根据现有标准《化工企业定量风险评价导则》(AQ/T 3046—2013)，当泄漏液体经储罐的泄漏孔流出时，其瞬时质量流量为：
其瞬时体积流量为：
由不可压缩液体连续性方程，可知：Qv＝Av(12)，
由公式(10)、(11)、(12)可得：
结合质量守恒定律及上式(13)，得到储罐内液面下降速度v1为：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；P为储罐内液体压力，Pa；P0为环境压力，Pa；hL为泄漏口上方液体高度，m；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；
步骤二：构建储罐内液面下降速度变化率a1模型
由于泄漏口上方储罐内液体高度为：hL＝h-h1-Δh (15)，
储罐内泄漏口上方压力为：P'＝P-P0 (16)，
根据(14)、(15)、(16)可得到：
进而推导出
引入a1为储罐内液面下降速度变化率，m/s2；对公式(18)进行微分后，令Δh对v1求导，得到：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；P为储罐内液体压力，Pa；P0为环境压力，Pa；hL为泄漏口上方液体高度，m；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；h1为储罐泄漏口距水平地面的高度，m；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；
步骤三：构建储罐内液面的下降速度v1的模型
对连续泄漏一段时间内储罐内液面下降速度变化率a1进行不定积分，可得储罐内液面下降速度v1：
其中，利用边界条件求积分常数C1，当t＝0时，v1为最大值，且Δh＝0，故有：
将公式(21)带入(20)中，即得到储罐内液面的下降速度：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；P’为储罐内泄漏口上方压力，Pa；h为储罐原始液体高度，m；h1为储罐泄漏口距水平地面高度，m；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；C1为积分常数；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s。
进一步地，所述泄漏系数C0是：液体从泄漏口流出液体的雷诺数和泄漏口直径的复杂函数得到的无量纲常数；按照如下情况取对应值：
⑴当雷诺数＞100时，泄漏口为圆形或多边形，则C0取0.65；泄漏口为三角形，则C0取0.60；泄漏口为长方形，则C0取0.55；
⑵当雷诺数≤100时，泄漏口为圆形或多边形，则C0取0.50；泄漏口为三角形，则C0取0.45；泄漏口为长方形，则C0取0.40；
进一步地，当储罐本体发生泄漏时，即刚开始发生泄漏时S为最大值；当储罐本体泄漏即将结束时，即储罐内水平液面下降到泄漏口位置，此时S∞0m，直至罐内液面高度低于泄漏口所在的高度位置，此时S＝0m，泄漏停止。
本发明具有以下有益效果：
1.本发明能够根据储罐内不同水平液面下降高度Δh得到其对应的泄漏口处泄漏液体落地水平距离S。能够根据此模型绘制出泄漏液体落地水平距离S随储罐内水平液面下降高度Δh关系变化曲线。
2.本发明能够根据储罐不同泄漏口距水平地面高度h1得到其所对应的泄漏口处泄漏液体落地水平距离S，同时可绘制出泄漏液体落地水平距离S与泄漏口距水平地面高度h1关系变化曲线。
3.将本发明运用到工业生产中，可用于各种储存液体储罐泄漏安全问题分析，可为事故预防提供有效的理论计算参考数据，为降低储罐区域的风险提供理论依据。
4.本发明方法适用于储存任意液体的立式储罐，尤其是化工企业涉及到生产、储存危化品的液体立式储罐，通过该方法，当储罐本体泄漏时，可确定液面下方任意泄漏口处泄漏液体落地水平距离S，包括最远水平距离Smax，该距离有助于确定具有易燃易爆、有毒泄漏液体波及范围，为救援防护安全距离的确定提供准确数据，从而科学制定事故应急救援方案、开展有效救援，减少救援事故损失。
附图说明
附图1为本发明的所述泄漏液体落地水平距离S与储罐内水平液面下降高度Δh的关系曲线。
附图2为本发明的所述泄漏液体落地水平距离S与储罐泄漏口距水平地面高度h1的关系曲线。
具体实施方式
下面通过实施例对本发明作进一步说明，但不作为对本发明的限制。
以中石化油库的汽油储罐为例，计算该汽油储罐泄漏口泄漏液体落地水平距离S：
(一)首先，收集所述汽油储罐相关参数来源于中石化某油库储罐信息，所述汽油储罐参数如下：储罐体积V为10000m3、内径D为30m、储罐高度H为19.341m、储罐内原有液体高度按h＝12m(储罐充装系数按0.84计算)；其具体相关参数见表1，
表1汽油储罐相关参数表
储罐类型 储存物料 储罐体积/m<sup>3</sup> 储罐充装系数 储罐内径/m 储罐高度/m 储罐内原始液体高度/m
常压立式储罐 汽油 10000 0.84 30 19.341 12
依据《化工企业定量风险评价导则》(AQ/T3046-2013)8.1.1中规定了泄漏场景，泄漏口径大小可分为完全破裂以及孔泄漏两大类，有代表性的泄漏场景及孔径范围见表2，
表2泄漏场景及泄漏口径取值
依据《事故调查与分析技术》(第二版)6.2.1.2中泄漏量的计算内容：若常压立式储罐发生泄漏的裂口是规则的，则裂口形状有圆形、多边形、三角形、长方形。液体泄漏系数C0取值见表3。
表3液体泄漏系数C0
(二)建立汽油储罐的泄漏模型
步骤一：构建泄漏口上方液体高度hL与泄漏后储罐内液面下降流速v1的模型
当汽油储罐内汽油发生泄漏时，储罐内液面下降的静压能和势能，转化为泄漏口液体流出的动能，符合机械能守恒定律。当泄漏液体经储罐的泄漏孔流出时，其瞬时质量流量为：
其瞬时体积流量为：
由不可压缩液体连续性方程，可知：Qv＝Av (12)，
由公式(10)、(11)、(12)可得：
结合质量守恒定律及上式(13)，得到储罐内液面下降速度v1为：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；P为储罐内液体压力，Pa；P0为环境压力，Pa；hL为泄漏口上方液体高度，m；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；
步骤二：构建储罐内液面下降速度变化率a1模型
由于泄漏口上方储罐内液体高度为：hL＝h-h1-Δh (15)，
储罐内泄漏口上方压力为：P'＝P-P0 (16)，
根据(14)、(15)、(16)可得到：
进而推导出
引入a1为储罐内液面下降速度变化率，m/s2；对公式(18)进行微分后，令Δh对v1求导，得到：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；P为储罐内液体压力，Pa；P0为环境压力，Pa；hL为泄漏口上方液体高度，m；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；h1为储罐泄漏口距水平地面的高度，m；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；
步骤三：构建储罐内液面的下降速度v1的模型
对连续泄漏一段时间内储罐内液面下降速度变化率a1进行不定积分，可得储罐内液面下降速度v1：
其中，利用边界条件求积分常数C1，当t＝0时，v1为最大值，且Δh＝0，故有：
将公式(21)带入(20)中，即得到储罐内液面的下降速度：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；P’储罐内泄漏口上方压力，Pa；h为储罐原始液体高度，m；h1为储罐泄漏口距水平地面高度，m；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；C1为积分常数；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s
步骤四：建立储罐泄漏口处泄漏出流液体速度v模型
当汽油储罐内汽油发生泄漏时，由质量守恒定律可知储罐内部的泄漏液面下降重量与储罐泄漏口处泄漏出的液体重量相同，即：A1v1＝Av (2)，
结合前式(1)、(2)可得储罐泄漏口处泄漏出流液体速度为：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；h为储罐原始液体高度，m；h1为储罐泄漏口距水平地面的高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s；
步骤五：建立储罐内原始液位h下降到泄漏口位置需要的泄漏时间t关系模型
当液体从储罐泄漏口泄漏后落至地面，若不考虑其他影响因素(风速、风向、恶劣天气状况、孔口拔塞等影响)，此时泄漏口处的液体受重力和储罐内部原始液位高度的压力作用，泄漏液体做类平抛运动；根据平抛运动模型，则液体由泄漏口落到地面的时间为：
泄漏口处的液体落至地面的水平泄漏距离为：
由公式(4)、(5)可得储罐内原始液位下降到泄漏口位置需要的泄漏时间为：
其中，h1为储罐泄漏口距水平地面的高度，m；v为储罐泄漏口处泄漏出流液体速度，m/s；C0为泄漏系数；A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；h为储罐原始液体高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s；
步骤六：建立储罐内水平液面下降高度Δh与泄漏口处液体落地水平距离S模型
储罐发生泄漏后罐内水平液面下降高度为：
将公式(7)带入到公式(6)，即得到储罐内水平液面下降高度Δh与所述泄漏口处液体落地水平距离S的关系：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；h1为储罐泄漏口距水平地面高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s；
根据已知条件储罐内原始液面高度h、泄漏口距水平地面高度h1、泄漏口形状及孔径，由本发明计算公式(8)计算储罐内不同的液面下降高度Δh所对应的泄漏口处液体落地水平距离S，实施例1至11为不同情况下的计算数值，其各实施例的计算结果见表4，其关系曲线如图1所示。
表4实施例计算结果
通过表4计算结果知，当储罐内液面高度h、距水平地面高度h1、泄漏口形状及孔径确定后，泄漏口处液体落地水平距离S随着储罐内液面下降高度Δh的增大而呈现减小的趋势，见图1。可知，其他同等条件下，泄漏刚刚开始时，即储罐内液面下降高度Δh为零时，泄漏口处液体落地水平距离S值最大，为5.3842m。
步骤七：建立储罐泄漏口距水平地面的高度h1与泄漏口处液体落地水平距离S模型
根据求根公式，继续推导公式(8)，且所述泄漏口处液体落地水平距离S的物理意义知其数值为正值，即得到所述泄漏口处液体落地水平距离S与储罐泄漏口距水平地面高度h1的关系为：
其中，h1为储罐泄漏口距水平地面的高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；C0为泄漏系数。
通过模拟储罐不同泄漏口距水平地面高度h1，计算出任意泄漏口处泄漏液体落地水平距离S；根据已知条件储罐内原始液面高度h，泄漏刚刚开始时，即储罐内液面下降高度Δh为零时，实施例12至33为在不同泄漏口距水平地面高度h1情况下，泄漏口处液体落地水平距离S的计算数值，其各实施例的计算结果见表5，其关系曲线如图2所示。
表5实施例计算结果
通过表5计算结果知，当储罐内液面高度h、泄漏口形状及孔径确定，泄漏口距水平地面高度h1不定时，泄漏口处液体落地水平最远距离S见图2。由图2及表5中可知，其他同等情况下，当泄漏口距水平地面高度h1为9.0m时，泄漏口处液体落地水平距离为最大值Smax为6.1740m。因此，该值Smax即为泄漏液体的最大波及范围，在设置储罐防火堤或为事故救援提供准确的安全防护距离，有利于开展有效救援，减少救援事故损失，在工业安全分析方面有广阔的应用前景。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围内。

Claims (4)

1.一种基于立式储罐本体泄漏口处液体落地水平距离S计算模型，其特征在于，包括建立立式液体储罐本体泄漏模型，确定储罐本体任意一泄漏口处，泄漏口距离水平地面的高度及该泄漏口流出液体的落地水平距离之间的关系；通过已知参数储罐泄漏口距水平地面的高度h1，计算出任意泄漏口处泄漏液体落地水平距离S；其计算关系公式为：
其中，公式中h1为储罐泄漏口距水平地面的高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；C0为泄漏系数。
2.根据权利要求1所述基于立式储罐本体泄漏口处液体落地水平距离S计算模型，其特征在于，由储罐泄漏口距水平地面的高度h1计算泄漏液体落地水平距离S模型建立，具体包括如下步骤：
步骤一：建立储罐泄漏口处泄漏出流液体速度v模型
液体泄漏时，储罐内液面下降的静压能和势能，转化为泄漏口液体流出的动能，根据机械能守恒定律，可得到所述储罐内液面的下降速度为：
由质量守恒定律可知储罐内部的泄漏液面下降重量与储罐泄漏口处泄漏出的液体重量相同，即：A1v1＝Av(2)，
结合前式(1)、(2)可得储罐泄漏口处泄漏出流液体速度为：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；h为储罐原始液体高度，m；h1为储罐泄漏口距水平地面的高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s；
步骤二：建立储罐内原始液位h下降到泄漏口位置需要的泄漏时间t关系模型
当液体从储罐泄漏口泄漏后落至地面，此时泄漏口处的液体受重力和储罐内部原始液位高度的压力作用，泄漏液体做类平抛运动；根据平抛运动模型，则液体由泄漏口落到地面的时间为：
泄漏口处的液体落至地面的水平距离为：
由公式(4)、(5)可得储罐内原始液位下降到泄漏口位置需要的泄漏时间为：
其中，h1为储罐泄漏口距水平地面的高度，m；v为储罐泄漏口处泄漏出流液体速度，m/s；C0为泄漏系数；A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；h为储罐原始液体高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s；
步骤三：建立储罐内水平液面下降高度Δh与泄漏口处液体落地水平距离S模型
储罐发生泄漏后罐内水平液面下降高度为：
将公式(7)带入到公式(6)，即得到储罐内水平液面下降高度Δh与所述泄漏口处液体落地水平距离S的关系：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；h1为储罐泄漏口距水平地面的高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s；
步骤四：建立储罐泄漏口距水平地面高度h1与泄漏口处液体落地水平距离S模型
根据求根公式，继续推导公式(8)，即得到所述泄漏口处液体落地水平距离S与储罐泄漏口距水平地面高度h1的关系：
其中，h1为储罐泄漏口距水平地面的高度，m；P’为储罐内泄漏口上方压力，Pa；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；C0为泄漏系数。
3.根据权利要求2所述基于立式储罐本体泄漏口处液体落地水平距离S计算模型，其特征在于，所述公式(1)储罐内液面的下降速度v1的模型建立过程包括如下步骤：
步骤一：构建泄漏口上方液体高度hL与泄漏后储罐内液面下降流速v1的模型
当泄漏液体经储罐泄漏孔流出时，其瞬时质量流量为：
其瞬时体积流量为：
由不可压缩液体连续性方程，可知：Qv＝Av(12)，
由公式(10)、(11)、(12)可得：
结合质量守恒定律及上式(13)，得到储罐内液面下降速度v1为：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；P为储罐内液体压力，Pa；P0为环境压力，Pa；hL为泄漏口上方液体高度，m；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；
步骤二：构建储罐内液面下降速度变化率a1模型
由于泄漏口上方储罐内液体高度为：hL＝h-h1-Δh(15)，
储罐内泄漏口上方压力为：P'＝P-P0(16)，
根据(14)、(15)、(16)可得到：
进而推导出
引入a1为储罐内液面下降速度变化率，m/s2；对公式(18)进行微分后，令Δh对v1求导，得到：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；P为储罐内液体压力，Pa；P0为环境压力，Pa；hL为泄漏口上方液体高度，m；h为储罐原始液体高度，m；Δh为储罐内液面下降高度，m；h1为储罐泄漏口距水平地面高度，m；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；
步骤三：构建储罐内液面下降速度v1的模型
对连续泄漏一段时间内储罐内液面下降速度变化率a1进行不定积分，可得储罐内液面下降速度v1：
其中，利用边界条件求积分常数C1，当t＝0时，v1为最大值，且Δh＝0，故有：
将公式(21)带入(20)中，即得到储罐内液面下降速度：
其中，A为储罐泄漏口面积，m2；A1为储罐水平横截面积，m2；P’为储罐内泄漏口上方压力，Pa；h为储罐原始液体高度，m；h1为储罐泄漏口距水平地面的高度，m；ρ为储罐内液体密度，kg/m3；g为重力加速度，m/s2；C0为泄漏系数；C1为积分常数；t为储罐原始液位下降到泄漏口位置需要的泄漏时间，s。
4.根据权利要求1所述基于立式储罐本体泄漏口处液体落地水平距离S计算模型，其特征在于，当储罐本体发生泄漏时，即刚开始发生泄漏时S为最大值；当储罐本体泄漏即将结束时，即储罐内水平液面下降到泄漏口位置，此时S∞0m，直至罐内液面高度低于泄漏口所在的高度位置，此时S＝0m，泄漏停止。
CN201910216521.3A 2019-03-21 2019-03-21 一种基于立式储罐本体泄漏口处液体落地水平距离s计算模型 Active CN109916348B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910216521.3A CN109916348B (zh) 2019-03-21 2019-03-21 一种基于立式储罐本体泄漏口处液体落地水平距离s计算模型

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910216521.3A CN109916348B (zh) 2019-03-21 2019-03-21 一种基于立式储罐本体泄漏口处液体落地水平距离s计算模型

Publications (2)

Publication Number Publication Date
CN109916348A true CN109916348A (zh) 2019-06-21
CN109916348B CN109916348B (zh) 2020-08-21

Family

ID=66966056

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910216521.3A Active CN109916348B (zh) 2019-03-21 2019-03-21 一种基于立式储罐本体泄漏口处液体落地水平距离s计算模型

Country Status (1)

Country Link
CN (1) CN109916348B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114722466A (zh) * 2022-03-31 2022-07-08 兵器工业卫生研究所 一种化学品围堰设计方法
CN115342994A (zh) * 2022-05-19 2022-11-15 北京大学 一种地下储氢气泄漏速度测量方法及测量装置
CN115342989A (zh) * 2022-07-21 2022-11-15 北京大学 一种氢气泄漏量测量方法及测量装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0473499A1 (fr) * 1990-08-30 1992-03-04 Sollac Tube protecteur de jet d'une installation de coulée continue et procédé en comportant application
JPH09318422A (ja) * 1996-05-31 1997-12-12 Kansei Corp 磁気式液面検出装置
JP2009002897A (ja) * 2007-06-25 2009-01-08 Showa Kiki Kogyo Co Ltd 液体貯蔵タンクの漏洩検知システムと漏洩検知方法
CN104297006A (zh) * 2014-09-25 2015-01-21 中国石油化工股份有限公司 Lng泄漏水幕抑制实验装置
CN207329366U (zh) * 2017-07-22 2018-05-08 江苏法利沃环保科技有限公司 带泄漏检测仪检测口的双层油罐
CN109033553A (zh) * 2018-07-04 2018-12-18 广西大学 一种基于常压立式储罐本体连续实时泄漏量的计算模型

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0473499A1 (fr) * 1990-08-30 1992-03-04 Sollac Tube protecteur de jet d'une installation de coulée continue et procédé en comportant application
JPH09318422A (ja) * 1996-05-31 1997-12-12 Kansei Corp 磁気式液面検出装置
JP2009002897A (ja) * 2007-06-25 2009-01-08 Showa Kiki Kogyo Co Ltd 液体貯蔵タンクの漏洩検知システムと漏洩検知方法
CN104297006A (zh) * 2014-09-25 2015-01-21 中国石油化工股份有限公司 Lng泄漏水幕抑制实验装置
CN207329366U (zh) * 2017-07-22 2018-05-08 江苏法利沃环保科技有限公司 带泄漏检测仪检测口的双层油罐
CN109033553A (zh) * 2018-07-04 2018-12-18 广西大学 一种基于常压立式储罐本体连续实时泄漏量的计算模型

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
滕天奎: "关于"液体从器壁流出离壁距离"的讨论", 《大学物理》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114722466A (zh) * 2022-03-31 2022-07-08 兵器工业卫生研究所 一种化学品围堰设计方法
CN115342994A (zh) * 2022-05-19 2022-11-15 北京大学 一种地下储氢气泄漏速度测量方法及测量装置
CN115342994B (zh) * 2022-05-19 2024-05-28 北京大学 一种地下储氢气泄漏速度测量方法及测量装置
CN115342989A (zh) * 2022-07-21 2022-11-15 北京大学 一种氢气泄漏量测量方法及测量装置
CN115342989B (zh) * 2022-07-21 2024-05-28 北京大学 一种氢气泄漏量测量方法及测量装置

Also Published As

Publication number Publication date
CN109916348B (zh) 2020-08-21

Similar Documents

Publication Publication Date Title
CN109916348A (zh) 一种基于立式储罐本体泄漏口处液体落地水平距离s计算模型
CN108038275A (zh) 一种天然气管道泄漏声场的数值仿真及特性分析方法
CN103605822B (zh) 石油化工码头液态化学品管道泄漏危险性动态分级方法
CN105805559A (zh) 液相管道泄漏不规则孔等效圆孔实验装置及实验方法
CN103527934B (zh) 一种管道泄漏实验装置及实验方法
CN109033553B (zh) 一种基于常压立式储罐本体连续实时泄漏量的计算模型
CN108171373A (zh) 一种化工园区毒气泄露逃生路径规划方法
CN108108525A (zh) 基于地理信息系统的气体泄漏事故模拟推演方法及装置
He et al. Simulation and application of a detecting rapid response model for the leakage of flammable liquid storage tank
CN112989717B (zh) 一种加氢站火灾事故热辐射评估方法及系统
CN102096768B (zh) 危险化学品泄漏污染评估方法
CN102880734A (zh) 基于cfd的飞机尾喷流大气扩散建模方法
CN109884250A (zh) 一种长贮库房内偏二甲肼扩散浓度分布的快速估算方法
CN108345759A (zh) 液化烃储存装置泄漏扩散区域浓度分布的测算方法
CN105550505A (zh) 一种确定气体扩散范围的方法及系统
Havens Review of dense gas dispersion field experiments
CN109297636A (zh) 储罐泄漏上风侧探测快速报警响应判断优化安装计算模型
CN111368445B (zh) 一种基于天然气热辐射的安全距离确定方法及系统
CN105868426B (zh) 一种基于烟团模型的重气持续泄露模拟方法
Das et al. Improving flammable mass estimation for vapor cloud explosion modeling in an offshore QRA
CN216247116U (zh) 一种集输立管两相流型在线分析实验系统
CN111611752B (zh) 一种非液体泄漏影响范围实时预测方法
CN110442091B (zh) 一种化工厂电子化应急预案及事故处置辅助方法
Li et al. Pressure distributions on prism-shaped buildings in experimentally simulated downburst
CN104573253A (zh) 压缩天然气加气站灾难后果预测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant