图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软件可进行区域性个人风险的计算。在这里不做详细介绍。 危险品园区的区域风险分析(5):http://www.751com.cn/guanli/lunwen_9676.html