环境影响评价程序1
——BOD-DO偶合矩阵模型一、程序参数说明：
1、1、待用户输入参数
N ——断面数；
Q0 ——上游（起始）河段的流量，m3/s；
L0 ——上游（起始）河段的BOD浓度，mg/L；
O0 ——上游（起始）河段的DO浓度，mg/ L；
Q ——各断面进入河流的污水流量，m3/s；
L ——各断面进入河流的BOD浓度，mg/L；
O ——各断面进入河流的DO浓度，mg/L；
Q3 ——各断面处的取水量，m3/s；
k1 ——各断面下游河段的BOD衰减速率常数（耗氧系数），d-1；
k2 ——各断面下游河段的大气复氧速率常数（复氧系数），d-1；
t ——各断面下游河段内的流行时间，d；
Os ——饱和溶解氧值，mg/L；
2、2、其他参数
Q1(i) ——上游进入断面i的流量，m3/s；
L1(i) ——上游进入断面i的BOD浓度，mg/L；
O1(i) ——上游进入断面i的DO浓度，mg/L；
Q2(i) ——由断面i输出到下游的流量，m3/s；
L2(i) ——由断面i输出到下游的BOD浓度，mg/L；
O2(i) ——由断面i输出到下游的DO浓度，mg/L；
aa(i) ——α(i)；
bb(i) ——β(i)；
r(i) ——γ(i)；
s(i) ——δ(i)；
二、程序
% 多河段BOD-DO耦合矩阵模型
N=input('N='); %输入断面数
Q=input('Q='); %输入从i至i+5断面支流流入的污水量
Q3=input('Q3='); %输入i至i+5断面支流流出的污水量
L=input('L=');
L=L';
O=input('O=');
O=O';
kd=input('kd='); %输入i至i+5相邻断面断面之间的有机物降解速
率常数
ka=input('ka='); %输入i至i+5相邻断面之间的大气复氧数率常数t=input('t='); %输入i至i+5相邻断面之间的有机物降解时间Q2(1)=input('Q20='); %输入从上游进入i断面的水流量
L20=input('L20='); %输入从i断面输入水的初始COD浓度
O20=input('O20='); %输入从主干道流入起始断面的水中的溶氧量的值
Os=input('Os='); %输入各断面的饱和溶氧量
for i=1:N
Q2(i+1)=Q2(i)+Q(i+1)-Q3(i+1);
Q1(i+1)=Q2(i);
end
for i=1:N;
aa(i)=exp(-kd(i)*t(i));
a(i)=aa(i)*(Q1(i+1)-Q3(i+1))/Q2(i+1);
b(i+1)=Q(i+1)/Q2(i+1);
r(i)=exp(-ka(i)*t(i));
c(i)=(Q1(i+1)-Q3(i+1))*r(i)/Q2(i+1);
bb(i)=kd(i)*(aa(i)-r(i))/(ka(i)-kd(i));
d(i)=bb(i)*(Q1(i+1)-Q3(i+1))/Q2(i+1);
s(i)=Os*(1-r(i));
f(i)=s(i)*(Q1(i+1)-Q3(i+1))/Q2(i+1);
end
aa
r
c
f=f';
g(1)=a(1)*L20;
g(2:N)=0;
g=g';
h(1)=c(1)*O20-d(1)*L20;
h(2:N)=0;
h=h';
A=eye(N);
for i=1:N-1
A(i+1,i)=-a(i+1);
end
for i=1:N
B(i,i)=b(i+1);
-0.0027 -0.0035 -0.0032 0 0
-0.0029 -0.0041 -0.0049 -0.0016 0
m = [1.8736 , 1.6580 , 1.4541 , 1.2962 , 1.1450]T
n = [ 6.1348, 6.4043 , 6.5532 , 6.6722 , 6.7594]T
L2 = [ 4.1464 , 7.5552 , 14.7068 , 18.0360 , 23.1092]T
O2 =[ 6.1348 , 6.2173 , 5.9537, 5.1114 , 4.4118]T
环境影响评价程序2
——多箱模型
一、程序参数说明：
1、待用户输入参数
L ——研究区域的长度，m
b ——研究区域的宽度，m
H ——混合层高度，m
m ——将高度方向离散的子高度个数
n ——将长度方向离散的子长度个数
u ——各高度层上的平均流速，m/s
E ——高度方向上相邻两层的湍流扩散系数，m2/s
C0——各高度层的本底浓度mg/m3
Q ——底层子箱的强源, mg/(m2·s)
2、其他参数
h ——子高度△h
l ——子长度△L
二、程序：
%大气质量模拟多箱模型
m=input('m='); %高度方向的子箱数
b=input('b='); %宽度，m
L=input('L='); %总长度，m
n=input('n='); %长度方向上的子箱个数
l=L/n; %子长度，m
H=input('H='); %子高度,m
h=H/n; %子长度,m
Q=input('Q='); %各子箱的强源,kg/s
Q=Q./(b*l); %各子箱源强2，kg/(m2.s)
E=input('E='); %高度方向上相邻两层的湍流扩散系数,m2/s u=input('u='); %高度方向上第i层的平均风速,m/s
C0=input('C0='); %高度方向第i层的污染物本底浓度 mg/m3 for i=1:m
a(i)=u(i)*h;
end
for i=1:m-1
e(i)=E(i)*l/h;
end
A(1,1)=a(1)+e(1);
A(m,m)=a(m)+e(m-1);
A(1,2)=-e(1);
A(m,m-1)=-e(m-1);
for i=2:m-1
A(i,i)=a(i)+e(i-1)+e(i);
A(i,i-1)=-e(i-1);
A(i,i+1)=-e(i);
end
a=a';
C0=C0';
D1=a.*C0;
D1(1)=D1(1)+Q(1)*l;
C1=inv(A)*D1;
D2=a.*D1;
D2(1)=D2(1)+Q(2)*l;
C2=inv(A)*D2;
D3=a.*D2;
D3(1)=D3(1)+Q(3)*l;
C3=inv(A)*D3;
D4=a.*D3;
D4(1)=D4(1)+Q(4)*l;
C4=inv(A)*D4;
CC=[C1,C2,C3,C4];
A %输出系数矩阵A、D
D
C %输出子箱浓度矩阵
三、程序运行
程序运行如下：
请输入将高度方向离散的子高度个数m=4
请输入研究区域的宽度，m b=1000
请输入研究区域的长度，m L=4000
请输入将长度方向离散的子长度个数n=4
请输入混合层高度, m H=120
输入底层子箱的强源, mg/(m2·s) Q=[0.1, 0.3, 0.4, 0.2]
输入10米处的风速,（m/s）u10=4
输入风速垂直分布指数p=0.07
输入高度方向上相邻两层的湍流扩散系数，m2/s E=[0.25,0.25,0.25,0.25]
输入各高度层的本底浓度mg/m3 C0=[0.1 0.1 0.2 0.3]
系数矩阵A=
129.4478 -8.3333 0 0
-8.3333 149.8239 -8.3333 0
0 -8.3333 154.7829 -8.3333
0 0 -8.3333 149.7719
系数矩阵D=
1.0e+003 *
0.1121 0.8484 1.1311 0.5641
0.0133 0.0198 0.0528 0.1139
0.0276 0.0279 0.0300 0.0353
0.0424 0.0417 0.0411 0.0408
子箱浓度矩阵C=
0.2946 0.2906 0.2884 0.2889
0.2023 0.2237 0.2557 0.2988
0.1488 0.5111 0.8555 1.0230 0.8757 6.5868 8.7933 4.4239。