图3.2.-3 区域风险评价网格示意图
3.2.3 个人风险评价模型
首先计算区域内每一危险源对每一网格中心产生的个人风险,然后对每一危险源产生的风险进行叠加,得到每一网格中心总的个人风险值。任意点处(x,y)通过将网格中心个人风险的离散结果进行内插得到。最后将个人风险值相等的点连接起来,便得到园区内不同水平个人风险的等值线[20]。
园区内任意网格中心的个人风险评价模型可表述为式(3):
(3)
式中:
IR(x,y),为危险源在空间地理坐标(x,y)处产生的个人风险;
f,为园区危险源发生事故的概率;
vs(x,y),为第s个危险源发生事故在位置(x,y)处引起的个体死亡概率;
M,为园区内危险源个数。
3.2.4 结合matlab计算个人风险
多数评价项目缺乏 AutoCAD 地图,在缺少地区坐标信息的情况使用 GIS 系统可以很好的解决该问题。通过建立重点位置坐标图并结合 matlab 在数值和图形处理的强大功能,能快速得出区域各点的个人风险值并绘制出个人风险等值线。
Matlab 个人风险程序编制思路如下:
网格化区域地图,每个网格中心点风险视为矩阵的一个元素。每个危险源导致的个人风险对应于一个矩阵,使用矩阵相加得到多个危险源的叠加风险。利用 surf 函数绘制出等值线图,或在需要时绘制三文立体图形。
算法程序如下:
m=input('输入行数:')
n=input('输入列数:')
c=input('输入网格边长,整数 m:')
risk1=input('输入个人风险值1×10-4:')
risk2=input('输入个人风险值1×10-5:')
risk3=input('输入个人风险值1×10-6:')
r1=input('输入个人风险值1×10-4半径,整数 m:')
r2=input('输入个人风险值1×10-5半径,整数 m:')
r3=input('输入个人风险值1×10-6半径,整数 m:')
x1=input('输入危险源 x 轴坐标:')
y1=input('输入危险源 y 轴坐标:')
A=zeros(m,n) %先将矩阵所有元素赋值 0
for f=1:n
x=(2*f-1)*c %取网格中心点横坐标
for g=1:m
y=(2*g-1)*c %取网格中心点纵坐标
r=sqrt((x-x1)^2+(y-y1)^2); %网格中心位置到危险源的距离
if r<=r1 %位于1×10-4区域内
A(g,f)=risk1
elseif r1<r<=r2 %位于1×10-5区域内
A(g,f)=risk2
elseif r2<r<=r3 %位于1×10-6区域内
A(g,f)=risk3
else %其它区域
A(g,f)=0
end
end
end
重复计算各个危险源风险矩阵,之后使用矩阵相加得总风险矩阵 M。利用 surf 函数绘制风险等值线:
f=c/2:c:(2n-1)*c/2;
g=c/2:c: (2m-1)*c/2;
surf(f,g,M)
3.2.5 利用CASST-QRA软件计算个人风险
中国安全生产科学研究院科研人员在国家“八五”、 “九五”和“十五”科技攻关成果的基础上研制开发的CASST-QRA软件可进行区域性个人风险的计算。在这里不做详细介绍。
- 上一篇:泸州老窖白酒行业网络营销模式分析
- 下一篇:《和讯网》营销策略的研究
-
-
-
-
-
-
-
g-C3N4光催化剂的制备和光催化性能研究
现代简约美式风格在室内家装中的运用
巴金《激流三部曲》高觉新的悲剧命运
浅析中国古代宗法制度
NFC协议物理层的软件实现+文献综述
上市公司股权结构对经营绩效的影响研究
江苏省某高中学生体质现状的调查研究
高警觉工作人群的元情绪...
中国传统元素在游戏角色...
C++最短路径算法研究和程序设计