水污染扩散

1 在线模拟示例

在线模拟示例

2 水污染扩散过程

地表水（surface water），是指陆地表面上动态水和静态水的总称，亦称“陆地水”，包括各种液态的和固态的水体，主要有河流、湖泊、沼泽、冰川、冰盖等。它是人类生活用水的重要来源之一，也是各国水资源的主要组成部分。

1.物理过程:物理过程主要是指污染物在水体中的混合稀释和自然沉淀过程。只改变进入水体污染物的物理性状、空间位置，而不改变其化学性质、不参与生物作用。
水体的混合稀释作用主要由下面三部分作用所致：紊动扩散、移流、离散。
（1）紊动扩散：由水流的紊动特性引起水中污染物自高浓度向低浓度区转移的扩散。
（2）移流：由于水流的推动使污染物的迁移随水流输移。
（3）离散：由于水流方向横断面上流速分布的不均匀（由河岸及河底阻力所致）而引起分散。
2.化学过程:污染物在水体中发生化学性质或形态、价态上的转化，使水质发生化学性质的变化。
主要包括酸碱中和、氧化-还原、分解-化合、吸附-解吸、胶溶-凝聚等过程
3.生物自净过程:是水体中的污染物经生物吸收、降解作用而发生消失或浓度降低的过程。
影响生物自净作用的关键是：溶解氧的含量，有机污染物的性质、浓度以及微生物的种类、数量等。

 3 一维模型

一维河流水质模型是目前应用最广泛的模型，它适用于污染物浓度仅在一个方向上有变化的场合，如宽度比较窄的河流。如果河段横截面、流速、流量、污染物的排入量和弥散系数不随时间变化，且污染物的衰减符合一级反应，河段不考虑源和汇。

4 二维模型

二维河流水质模型在如下条件下使用：
（1）平直、断面形状规则河段混合过程段；
（2）持久性污染物；    
（3）河流为恒定流动；
（4）连续稳定排放； 

5 部分代码

package com.planet.engine.water.param;

import lombok.Data;

@Data

public class BusD1Param {

private double ex = 0.045;

private double c0 = 10d;

private double k = 0.045d;

private double k1 = 0.194d;

private double bupper = 100d;

private double aupper = 300d;

}

@RestController

@RequestMapping("/api/water/1d")

@Api(value = "WaterD1SimuApi", description = "水污染扩散模拟-一维模型")

public class WaterD1SimuApi extends MyBaseApi {

@Autowired

private WaterDiffService waterDiffService;

@Autowired

private SimuHistService simuHistService;

@LimitRequest

@ResponseBody

@RequestMapping(value = "/mixing-length", method = RequestMethod.POST)

@ApiOperation(value = "一维-混合过程段长度")

public Result d1SimuMixLength(@RequestBody WaterD1MixParam params) {

if (params.getProcParams().getGridLen() != 0

& params.getProcParams().getSimuLength() / params.getProcParams().getGridLen() > 1000) {

throw new DataOptionException("云主机资源有限，请缩短模拟时间和步长距离!");

}

D1MixWrapData result = this.waterDiffService.mixLength(params);

SecurityUser securityUser = this.getLoginUser();

String paramsStr = this.getReqParams(params);

this.simuHistService.save(securityUser, SimuType.WATER_1D_MIXING_LENGTH, paramsStr, StatusType.SUCCESS.code(), result);

return this.getResult(ResultCode.SUCCESS, result);

}

@LimitRequest

@ResponseBody

@RequestMapping(value = "/", method = RequestMethod.POST)

@ApiOperation(value = "一维模式-通用水质模拟模型")

public Result d1Simu(@RequestBody WaterD1Param params) {

if (params.getProcParams().getGridLen() != 0

& params.getProcParams().getSimuLength() / params.getProcParams().getGridLen() > 1000) {

throw new DataOptionException("云主机资源有限，请缩短模拟时间和格网距离!");

}

String result = this.waterDiffService.d1Diff(params);

SecurityUser securityUser = this.getLoginUser();

String paramsStr = this.getReqParams(params);

this.simuHistService.save(securityUser, SimuType.WATER_1D, paramsStr, StatusType.SUCCESS.code(), result);

return this.getResult(ResultCode.SUCCESS, result);

}

@LimitRequest

@ResponseBody

@RequestMapping(value = "/std/continue", method = RequestMethod.POST)

@ApiOperation(value = "一维模式-水质模拟模型(持续排放)")

public Result d1SimuStd(@RequestBody WaterD1Param params) {

if (params.getProcParams().getGridLen() != 0

& params.getProcParams().getSimuLength() / params.getProcParams().getGridLen() > 1000) {

throw new DataOptionException("云主机资源有限，请缩短模拟时间和格网距离!");

}

String result = this.waterDiffService.d1DiffStd(params);

SecurityUser securityUser = this.getLoginUser();

String paramsStr = this.getReqParams(params);

this.simuHistService.save(securityUser, SimuType.WATER_1D_CONTINUE, paramsStr, StatusType.SUCCESS.code(), result);

return this.getResult(ResultCode.SUCCESS, result);

}

@LimitRequest

@ResponseBody

@RequestMapping(value = "/std/instant", method = RequestMethod.POST)

@ApiOperation(value = "一维模式-水质模拟模型(瞬时排放)")

public Result d1SimuInstantStd(@RequestBody WaterD1InstParam params) {

if (params.getProcParams().getGridLen() != 0

& params.getProcParams().getSimuLength() / params.getProcParams().getGridLen() > 1000) {

throw new DataOptionException("云主机资源有限，请缩短模拟时间和格网距离!");

}

D1InstWrapData result = this.waterDiffService.d1SimuInstantStd(params);

SecurityUser securityUser = this.getLoginUser();

String paramsStr = this.getReqParams(params);

this.simuHistService.save(securityUser, SimuType.WATER_1D_INSTANT,

paramsStr, StatusType.SUCCESS.code(), result);

return this.getResult(ResultCode.SUCCESS, result);

}

@LimitRequest

@ResponseBody

@RequestMapping(value = "/s-p/bod", method = RequestMethod.POST)

@ApiOperation(value = "一维模式-S-P水质模拟模型(BOD浓度)")

public Result d1SpBodSimu(@RequestBody WaterD1SpParam params) {

if (params.getProcParams().getGridLen() != 0

& params.getProcParams().getSimuLength() / params.getProcParams().getGridLen() > 1000) {

throw new DataOptionException("云主机资源有限，请缩短模拟时间和步长距离!");

}

String result = this.waterDiffService.spBodDiff(params);

SecurityUser securityUser = this.getLoginUser();

String paramsStr = this.getReqParams(params);

this.simuHistService.save(securityUser, SimuType.WATER_1D_SP, paramsStr,

StatusType.SUCCESS.code(), result);

return this.getResult(ResultCode.SUCCESS, result);

}

@LimitRequest

@ResponseBody

@RequestMapping(value = "/s-p/odis", method = RequestMethod.POST)

@ApiOperation(value = "一维模式-S-P水质模拟模型(溶氧量)")

public Result d1SpOdDisSimu(@RequestBody WaterD1SpOdParam params) {

if (params.getProcParams().getGridLen() != 0

& params.getProcParams().getSimuLength() / params.getProcParams().getGridLen() > 1000) {

throw new DataOptionException("云主机资源有限，请缩短模拟时间和格网距离!");

}

D1SpWrapData result = this.waterDiffService.spOdDisDiff(params);

SecurityUser securityUser = this.getLoginUser();

String paramsStr = this.getReqParams(params);

this.simuHistService.save(securityUser, SimuType.WATER_1D_SP_ODIS, paramsStr,

StatusType.SUCCESS.code(), result);

return this.getResult(ResultCode.SUCCESS, result);

}

@LimitRequest

@ResponseBody

@RequestMapping(value = "/s-p/odif", method = RequestMethod.POST)

@ApiOperation(value = "一维模式-S-P水质模拟模型(亏氧量)")

public Result d1SpOdDifSimu(@RequestBody WaterD1SpOdParam params) {

if (params.getProcParams().getGridLen() != 0

& params.getProcParams().getSimuLength() / params.getProcParams().getGridLen() > 1000) {

throw new DataOptionException("云主机资源有限，请缩短模拟时间和格网距离!");

}

D1SpWrapData result = this.waterDiffService.spOdDifDiff(params);

SecurityUser securityUser = this.getLoginUser();

String paramsStr = this.getReqParams(params);

this.simuHistService.save(securityUser, SimuType.WATER_1D_SP_ODIF, paramsStr,

StatusType.SUCCESS.code(), result);

return this.getResult(ResultCode.SUCCESS, result);

}

}

public static void main0(String[] args) {

File file = new File("D:\fenghe-desolveriver-line.shp");

SimpleFeatureCollection featureCollection = null;

try {

ShapefileDataStore dataStore = new ShapefileDataStore(file.toURI().toURL());

dataStore.setCharset(StandardCharsets.UTF_8);

System.out.println(dataStore.getCharset() + "获取编码11111111111111111");

SimpleFeatureSource featureSource = dataStore.getFeatureSource();

featureCollection = featureSource.getFeatures();

SimpleFeatureIterator iterator1 = (SimpleFeatureIterator) featureCollection.features();

SimpleFeature simpleFeature = iterator1.next();

MultiLineString geom = (MultiLineString) simpleFeature.getDefaultGeometry();

LineString river = (LineString) geom.getGeometryN(0);

Coordinate endPt = river.getCoordinates()[river.getCoordinates().length - 1];

iterator1.close();

ShapefileDataStore ds = new ShapefileDataStore(s1.toURI().toURL());

ds.setCharset(Charset.forName("GBK"));

System.out.println(ds.getCharset() + "获取编码22222222");

SimpleFeatureType sft = fs.getSchema();

ds.createSchema(sft);

FeatureWriter<SimpleFeatureType, SimpleFeature> writer = ds.getFeatureWriter(ds.getTypeNames()[0], Transaction.AUTO_COMMIT);

while (iterator.hasNext()) {

SimpleFeature feature = iterator.next();

SimpleFeature featureBuf = writer.next();

featureBuf.setAttributes(feature.getAttributes());

Geometry geo = (Geometry) feature.getAttribute("the_geom");

featureBuf.setAttribute("the_geom", geo);

}

writer.write();

writer.close();

iterator.close();

} catch (Exception e) {

e.printStackTrace();

}

}

相关知识

水污染防治工作简报
中华人民共和国水污染防治法
水污染治理标语
..::刘洪岩：《水污染防治法》新解
城市湿地水生植物治理水污染应用分析
室内环境警示：警惕家庭装修造成饮用水污染
水污染对土壤有哪些影响
大型水生植物在水污染治理中的应用
湘教版地理选修6水污染及其防治2(编辑修改稿)
涟源市人大调研六亩塘镇水污染防治情况

网址: 水污染扩散 https://www.huajiangbk.com/newsview1352899.html