地下水质量评价的计算程序

c zhiliang.for，根据水样分析报告数据计算“综（复）合评价分值F（i）”，

c 遵照中华人民共和国国家标准《地下水质量标准》（GBT14848-93）

c 程序中有m个水样，每个水样最多可含n＝39项［综（复）合评价分值所要求的］检测组分

c 当某水样实测项目个数nt小于39项时，其所有单项评价分值的算术平均值FI用分母nt计算

c 某水样实测项目nt通过计数器把非-1项累加求得，其间剔除了-1项—即未测项

c 某水样实测项目所有单项评价分值的最大值fmax通过循环比较获得

c x（m），y（m）-水样点坐标，id（m）-水样点统一编号，c（n，m）-m个水样点n项检测组分的测试值mg/L，缺测项赋-1占位，c0（n，5）-n项检测组分（GBT14848-93）的一至五级水标准值mg/L，F（m）-m个水样每个水样的综（复）合评价分值悄搏，Fi（n，m）-m个水样每个水样所有检测组分的单项评价分值，Fs（m）-m个水样每个水样所有检测组分单项评价分值的算术平均值，fmax（m）-m个水样每个水样所有单项评价分值中挑出的最大值，nt（m）m个水样每个水样所有实测项目的个数，

c jmax（m）-每个水样挑出的单项评价分值中最大值是哪一项 parameter（m＝9，n＝39）

dimension x（m），y（m），id（m），c（n，m），c0（n，5），f（m），fi（n，m），fs（m）

dimension fmax（m），nt（m），jmax（m）

c先读进m个水样的编号、坐标、“色（度）、嗅和味、浑浊旦此度（度）、肉眼可见物、pH值”5个单项评分值Fi（n，m），缺测项-1占位open（1，file＝’fi1-5.dat’）

read（1，*）（id（j），j＝1，m）

read（1，*）（x（j），j＝1，m）

read（1，*）（y（j），j＝1，m）

do 10 i＝1，5

10 read（1，*）（fi（i，j），j＝1，m）

close（1）

c 再读进m个水样第6～39项组分的测试值mg/L，缺测项-1占位

open（1，file＝’c6-39.dat’）

do 20 i＝6，n

20 read（1，*）（c（i，j），j＝1，m）

close（1）

c 读进第6～39项检测组分（GBT14848-93）的一至五级水标准值mg/L—c0（n，5）

open（1，file＝’c0.dat’）

do 30 i＝6，n

30 read（1，*）ia，（c0（i，j），j＝1，5）

close（1）

c 第100句循环体对第j个水样点第i（i＝1—5）项组分的单项评分值Fi（1-5，m）统计有值项个数nt（j）

c 找最大项fmax（j）及其项目序号jmax（j）、累加有值项分值之和fs（j）

do 100 j＝1，m

nt（j）＝0

fmax（j）＝0

f（j）＝0

fs（j）＝0

do 90 i＝1，5

c 如果第i项组分是缺测项，则跳出i循环计算下一项组分

if（fi（i，j）.eq.-1）goto 90

nt（j）＝nt（j）+1

fs（j）＝fs（j）+fi（i，j）

if（fi（i，j）.le.fmax（j））goto 90

fmax（j）＝fi（i，j）

jmax（j）＝i

90 continue

100 continue

c 第170句循环体对第j个水样点第i（i＝6—39）项组分的单项评分值Fi（1-5，m）统计有值项个数nt（j）

c 找最大项fmax（j）及其项目序号jmax（j）、累加有值项分值之和fs（j）

c 第200句循环体对第j个水样求nt（j）项组分分值和fs（j）的均值fs（j）、综（复）合评价分值f（j）

do 200 j＝1，m

do 170 i＝6，模运迅n

fi（i，j）＝-1

c 如果第i项组分是缺测项则分值保留-1，且跳出i循环计算下一项组分：

if（c（i，j）.eq.-1）goto 170

nt（j）＝nt（j）+1

c 第j个水样点第i项组分为1级，分值为零：

if（c（i，j）.gt.c0（i，1））goto 110

fi（i，j）＝0

goto 150

c 第j个水样点第i项组分为1级，分值为1

110 if（c（i，j）.gt.c0（i，2））goto 120

fi（i，j）＝1

goto 150

120 if（c（i，j）.gt.c0（i，3））goto 130

fi（i，j）＝3

goto 150

130 if（c（i，j）.gt.c0（i，4））goto 140

fi（i，j）＝6

goto 150

140 fi（i，j）＝10

150 fs（j）＝fs（j）+fi（i，j）

if（fi（i，j）.le.fmax（j））goto 170

fmax（j）＝fi（i，j）

jmax（j）＝i

170 continue

fs（j）＝fs（j）/nt（j）

f（j）＝sqrt（fs（j）**2/2+fmax（j）**2/2）

200 continue

open（1，file＝’fenzhi.dat’）

c 输出：水样点坐标x（i），y（i），综（复）合评价分值f（i），水样编号id（i）

do 210 j＝1，m

210 write（1，220）x（j），y（j），f（j），id（j）

close（1）

220 format（2f8.2，1x，f8.2，1x，i4）

c 输出：水样点单项分值fi（j），水样编号id（j），每个水样实测项目个数nt（j）

c 每个水样单项分值中最大值fmax（j），每个水样的单项分值最大值是哪一项jmax（j）

open（1，file＝’fi.dat’）

write（1，240）（id（j），j＝1，m）

do 230 i＝1，n

230 write（1，250）i，（fi（i，j），j＝1，m）

write（1，260）（nt（j），j＝1，m）

write（1，270）（fmax（j），j＝1，m）

write（1，280）（jmax（j），j＝1，m）

close（1）

240 format（’水样编号’，9i5）

250 format（i8，9f5.0）

260 format（’测试项数’，9i5）

270 format（’最大分值’，9f5.0）

280 format（’最大项号’，9i5）

stop

end

一御档物、地下水环境质量标准是什么

1、地下水环境质量标准 ：

（1）Ⅰ类主要反映地下水化学组分的天然低背景含量。适用于各种用途；

（2）Ⅱ类主要反映地下水化学组分的天然背景含量。适用于各种用途；

（3）Ⅲ类以人体健康基准值为依据。主要适用于集中式生活饮用水水源及工、农业用水；

（4）Ⅳ类以农业和工业用水要求为依据。除适用于农业和部分工业用水外，适当处理后可作生活饮用水；

（5）Ⅴ类不宜饮用，其他用水可根据使用目的选用。

2、法律依据：《中华人民共和国水污染防治法》第十三条

国务院环境保护主管部门会同国务院水行政主管部门和有关省、自治区、直辖市人民政府，可以根据国家确定的重要江河、湖泊流域水体的使用功能以及有关地镇液区的经济、技术条件，确定该重要江河、湖泊流域的省界水体适用的水环境质量标准，报国务院批准后施行。

第十五条

国务院环境保护主管部门和省、自治区、直辖市人民政府，应当根据水污染防治的要求和国家或者地方的经济、蠢桐技术条件，适时修订水环境质量标准和水污染物排放标准。

二、地下水环境质量有哪些调查程序

地下水环境质量调查程序：

1、收集已有的定位水质监测资料，确定重点调查地区，制定调查计划；

2、进行现场查勘，了解污染源的分布情况，估算废污水排放量和有机农药使用量，对污染严重的河段和水井进行取样分析；

3、将水质调查资料与定位监测资料相结合，对水体水质概况进行评价，提出控制污染的建议。

---