CN109754182A - A kind of calculation method and system of contaminated site soil remediation amount - Google Patents
A kind of calculation method and system of contaminated site soil remediation amount Download PDFInfo
- Publication number
- CN109754182A CN109754182A CN201811648129.8A CN201811648129A CN109754182A CN 109754182 A CN109754182 A CN 109754182A CN 201811648129 A CN201811648129 A CN 201811648129A CN 109754182 A CN109754182 A CN 109754182A
- Authority
- CN
- China
- Prior art keywords
- soil
- contaminated site
- area
- repaired
- sample data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Processing Of Solid Wastes (AREA)
Abstract
The present invention provides the calculation methods and system of a kind of contaminated site soil remediation amount, comprising: according to the soil sample data of contaminated site, obtains the pollution risk distribution map of different soil;The region for selecting pollution risk grade to reach predetermined level simplifies area's boundary profile to be repaired as area to be repaired, and according to preset simplified tolerance;According to soil sample data, piecemeal zoning is carried out using contaminant transportation depth of the Thiessen polygon to different soil, obtains contaminant transportation depth block diagram;The pollution risk distribution map of identical soil layer and contaminant transportation depth block diagram are overlapped, contaminated site soil remediation amount is calculated according to the area to be repaired on stacking chart, and obtain the inflection point coordinate in area to be repaired, constructed for soil remediation.Operation of the present invention is simple, is convenient for practical application, can scientifically and rationally calculate contaminated site soil remediation amount, reduce place soil repairing and treating difficulty, improve construction efficiency.
Description
Technical field
The present invention relates to space enrironment technique for investigation field, the calculation method of espespecially a kind of contaminated site soil remediation amount and
System.
Background technique
China's industrial and mineral area environment quality allows of no optimist.It was shown according to 2014 " national Soil Pollution Investigation bulletin ",
Polluted enterprises land used, Industrial Wasteland, industrial park, solid waste centralized processing place and minery soil point exceeding standard rate
Respectively 36.3%, 34.9%, 29.4%, 21.3% and 33.4%.And environmental pollution improvement's total investment in 2016 accounts for the country
Total output value specific gravity is relatively low, is 1.24% (data are derived from " China Statistical Yearbook in 2017 "), this is required using limited
Investment carries out pollution control.Environment remediation cost is reduced, it is extremely urgent to improve control quality.
It is a kind of to reduce later period environment remediation treatment cost by the space enrironment technique for investigation of exploiting economy adaptive type
Feasible method.Identification to contaminated site soil remediation side amount is the key that space enrironment technique for investigation.Currently, soil remediation
The identification just measured is portrayed based on the three-dimensional simulation of pollutant transition process in the soil, and common soil pollutant three-dimensional is moved
Condition of shifting one's love can by EVS (EarthVolumetricStudio, visualize geology modeling software) environment visualization system and
The three-dimensional irregular figure that Voxler software is realized, however depicted is not particularly suited for actual soil remediation construction work
Journey carries out soil according to the three-dimensional irregular figure and excavates work, then needs constantly to adjust and excavate opening for shape and each layer
Depth is dug, and to each layer of cutting depth without related science foundation, in addition repairs the more curves in boundary, also can artificially cause to repair
Engineering construction difficulty increases, project amount is bigger than normal.
Summary of the invention
The object of the present invention is to provide the calculation method and system of a kind of contaminated site soil remediation amount, by will be complicated
Pollutant distribution threedimensional model is simplified, and is allowed to be suitble to practice of construction demand, has stronger operability, reduce place
Soil remediation difficulty of governance, improves construction efficiency.
Technical solution provided by the invention is as follows:
A kind of calculation method of contaminated site soil remediation amount, comprising: according to the soil sample data of contaminated site, obtain
The pollution risk distribution map of the different soil of the contaminated site;According to the pollution risk distribution map, pollution risk etc. is selected
Grade reaches the region of predetermined level as area to be repaired, and carries out letter to area's boundary profile to be repaired according to preset simplified tolerance
Change;According to the soil sample data, piecemeal zoning is carried out using contaminant transportation depth of the Thiessen polygon to different soil,
Obtain contaminant transportation depth block diagram;By the pollution risk distribution map of identical soil layer and the contaminant transportation depth point
Block figure is overlapped, and calculates contaminated site soil remediation amount according to the area to be repaired on stacking chart, and obtain the area to be repaired
Inflection point coordinate, for subsequent soil remediation construct.
It is further preferred that the soil sample data according to contaminated site, obtain the difference soil of the contaminated site
The pollution risk distribution map of layer specifically includes: according to place geotechnical stratified situation, carrying out to the soil sample data of contaminated site
Processing, obtains each soil layer sample data;According to each soil layer sample data, other are gone out using inverse distance weighting interpolation and is not surveyed
The data for measuring position, to obtain the pollutant concentration of each position of each soil layer;According to the pollution of each position of each soil layer
Object concentration classifies to pollution risk referring to soil background and soil pollutant concentration limit, obtains corresponding soil layer
Pollution risk distribution map.
It is further preferred that the soil sample data to contaminated site are handled, each soil layer sample number is obtained
According to including: from the soil sample data of contaminated site, by all soil sample data for belonging to the same sampled point of same soil layer
Averaging is summed up, gained mean value corresponds to the sample data of sampled point as corresponding soil layer;The all of same soil layer will be belonged to
The sample data of sampled point constitutes corresponding soil layer sample data.
It is further preferred that described simplify area's boundary profile to be repaired according to preset simplified tolerance, it is specific to wrap
It includes: based on Douglas-Pu Ke algorithm, extra folding being removed to area's boundary profile to be repaired according to maximum allowable offset is preset
Point obtains first and simplifies boundary profile;Based on Zhou-Jones algorithm, simplify boundary to described first according to default minimum area
Profile removes extra break.
It is further preferred that it is described preset maximum allowable offset be 2-4 meter, the default minimum area for 20-40 put down
Square rice.
It is further preferred that obtaining the difference of the contaminated site in the soil sample data according to contaminated site
It include: to carry out soil sample in multiple sampled points of contaminated site respectively before the pollution risk distribution map of soil layer, by obtaining
The soil sample obtained carries out pollutant concentration test, obtains soil sample data.
The present invention also provides a kind of computing systems of contaminated site soil remediation amount, comprising: pollution risk evaluation module is used
In the soil sample data according to contaminated site, the pollution risk distribution map of the different soil of the contaminated site is obtained;It is to be repaired
Multiple area simplifies module, for according to the pollution risk distribution map, the region for selecting pollution risk grade to reach predetermined level to be made
For area to be repaired, and area's boundary profile to be repaired is simplified according to preset simplified tolerance;Transfer Depth evaluation module is used
According to the soil sample data, piecemeal zoning is carried out using contaminant transportation depth of the Thiessen polygon to different soil,
Obtain contaminant transportation depth block diagram;Reparation amount computing module, by the pollution risk distribution map of identical soil layer and described
Contaminant transportation depth block diagram is overlapped, and calculates contaminated site soil remediation amount according to the area to be repaired on stacking chart, and
The inflection point coordinate in the area to be repaired is obtained, is constructed for subsequent soil remediation.
It is further preferred that the pollution risk evaluation module, is further used for according to place geotechnical stratified situation, to dirt
The soil sample data in dye place are handled, and each soil layer sample data is obtained;According to each soil layer sample data, using anti-
Distance weighting method interpolation goes out the data of other unmeasured positions, to obtain the pollutant concentration of each position of each soil layer;Root
According to the pollutant concentration of each position of each soil layer, referring to soil background and soil pollutant concentration limit to pollution risk into
Row classification obtains the pollution risk distribution map of corresponding soil layer.
It is further preferred that the area to be repaired simplifies module, it is further used for pressing based on Douglas-Pu Ke algorithm
Extra break is removed to area's boundary profile to be repaired according to maximum allowable offset is preset, first is obtained and simplifies boundary profile;Base
In Zhou-Jones algorithm, simplifies boundary profile to described first according to default minimum area and remove extra break.
It is further preferred that further include: sample data obtains module, for clicking through respectively in multiple samplings of contaminated site
Row soil sample obtains soil sample data by carrying out pollutant concentration test to the soil sample of acquisition.
The calculation method and system of a kind of contaminated site soil remediation amount provided through the invention, can bring with down toward
It is few a kind of the utility model has the advantages that
1, for the present invention by simplifying complicated pollutant distribution threedimensional model, division obtains simple contaminated soil
Block is allowed to be suitble to practice of construction demand, easy to operate, is convenient for practical application, reduces place soil repairing and treating difficulty, mention
High construction efficiency.
2, the present invention fully considers practice of construction demand, obtains to contaminated soil block inflection point coordinate, convenient for practical
The design and progress of soil remediation construction work.
3, repair coverage boundary design is closing straightway by the present invention, and solving same court restoring area boundary is curve
The problem of, the implementation of recovery project is facilitated, more there is operability.
Detailed description of the invention
Below by clearly understandable mode, preferred embodiment is described with reference to the drawings, a kind of contaminated site soil is repaired
The calculation method of complexor and above-mentioned characteristic, technical characteristic, advantage and its implementation of system are further described.
Fig. 1 is a kind of flow chart of one embodiment of the calculation method of contaminated site soil remediation amount of the invention;
Fig. 2 is a kind of flow chart of another embodiment of the calculation method of contaminated site soil remediation amount of the invention;
Fig. 3 is a kind of structural representation of one embodiment of the computing system of contaminated site soil remediation amount of the invention
Figure;
Fig. 4 is a kind of structural representation of another embodiment of the computing system of contaminated site soil remediation amount of the invention
Figure;
Fig. 5 is the first soil in a kind of another embodiment of the calculation method of contaminated site soil remediation amount of the invention
Layer pollution risk distribution map (left side) and the second soil layer pollution risk distribution map (right side);
Fig. 6 is the first soil in a kind of another embodiment of the calculation method of contaminated site soil remediation amount of the invention
Layer simplification figure and the second soil layer simplification figure.
Drawing reference numeral explanation:
100. sample data obtain module, 110. pollution risk evaluation modules, 120. Transfer Depth evaluation modules, 130. to
It repairs area and simplifies module, 140. reparation amount computing modules.
Specific embodiment
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, Detailed description of the invention will be compareed below
A specific embodiment of the invention.It should be evident that drawings in the following description are only some embodiments of the invention, for
For those of ordinary skill in the art, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing, and obtain other embodiments.
To make simplified form, part related to the present invention is only schematically shown in each figure, they are not represented
Its practical structures as product.In addition, there is identical structure or function in some figures so that simplified form is easy to understand
Component only symbolically depicts one of those, or has only marked one of those.Herein, "one" is not only indicated
" only this ", can also indicate the situation of " more than one ".
In one embodiment of the invention, as shown in Figure 1, a kind of calculation method of contaminated site soil remediation amount, packet
It includes:
Step S100 obtains the pollution wind of the different soil of the contaminated site according to the soil sample data of contaminated site
Dangerous distribution map.
Specifically, soil sample data are to carry out experimental analysis by soil sampling on the spot to the soil of sampling and obtain.Soil
Earth sample data reflects pollutant concentration contained by soil.Consider human and material resources cost, the sampling point of soil sampling is limited on the spot
, according to soil sample data, the pollutant concentration of unmeasured point is obtained by space interpolation, to obtain the dirt of different soil
Dye distribution;Definition according to the pollutant concentration of each position of each soil layer and country or professional standard to pollution risk grade,
Obtain the pollution risk distribution map of each soil layer.
Step S200 according to the pollution risk distribution map, select pollution risk grade reach the region of predetermined level as
Area to be repaired, and area's boundary profile to be repaired is simplified according to preset simplified tolerance.
Specifically, pollution risk grade is such as divided into three-level, it is respectively high, medium, low;If predetermined level is
It is medium, then select pollution risk grade to reach medium or more region, i.e., high, medium region is as area to be repaired;Such as
Fruit predetermined level be it is high, then selecting pollution risk grade is high region as area to be repaired.To soil remediation amount shadow
In the case where ringing less, area's boundary profile to be repaired is simplified, in favor of the construction of subsequent soil remediation.Road can be used
Lattice Lars-Pu Ke algorithm, and/or Wang-M ü ller algorithm, and/or Zhou-Jones algorithm to area's boundary profile to be repaired into
Row simplifies.Douglas-Pu Ke algorithm retains the key point for constituting the basic configuration of boundary profile, and removes every other point;
Each break will be deleted at a distance from Trendline less than the break of tolerance to the vertical range of Trendline on the algorithm measurement boundary
It removes.Wang-M ü ller algorithm finds out borderline bending by shape recognition techniques and analyzes its feature, then eliminates nothing and fastens
The bending wanted.Zhou-Jones algorithm identifies effective triangle area of each break on boundary, then uses a series of measurements
These triangles are weighted to compare the flatness of each area, the degree of bias and convexity, instruct to move by the area after weighting
Except corresponding break to simplify facial contour, while retention characteristic as much as possible.
Step S300 is according to the soil sample data, using Thiessen polygon to the contaminant transportation depth of different soil
Piecemeal zoning is carried out, contaminant transportation depth block diagram is obtained.
Specifically, each sample has corresponding Transfer Depth on practical contaminated site, and sample is opposite
It is discrete, and the Transfer Depth of each sample is different, so using Thiessen polygon to the pollutant of different soil in place
Transfer Depth carries out piecemeal zoning, and rough operation principle and process is as follows:
Thiessen polygon is a kind of subdivision method to investigation range plane, its main feature is that any position in polygon from
The distance of sampled point is nearest in the polygon, remote with a distance from sampled point in adjacent polygons, and each polygon includes and only
Include a sampled point.Thiessen polygon is constructed as follows:
1. according to sequential scan input point (i.e. sampled point) from left to right, from top to bottom.
Meet delaunay criterion 2. marking off in all the points (not comprising other points in i.e. each triangle circumscribed circle)
Irregular triangle network.
3. doing perpendicular bisector to each side of all triangles, these perpendicular bisectors form the side of Thiessen polygon, each to hang down
The intersection point of straight bisector determines the position of Thiessen polygon break.
Step S400 carries out the pollution risk distribution map of identical soil layer and the contaminant transportation depth block diagram
Superposition calculates contaminated site soil remediation amount according to the area to be repaired on stacking chart, and the inflection point for obtaining the area to be repaired is sat
Mark is constructed for subsequent soil remediation.
Specifically, the pollution risk distribution map of identical soil layer and contaminant transportation depth block diagram are overlapped, obtain
The stacking chart of each soil layer.Area to be repaired on the stacking chart of each soil layer includes several independent regions, and each region is wrapped again
Containing several piecemeals, the piecemeal in each region is obtained according to Thiessen polygon again, and it is deep that each piecemeal corresponds to a kind of contaminant transportation
Degree, obtains the volume of piecemeal according to the area of each piecemeal and contaminant transportation depth, by all areas to be repaired of different soil
The volume of interior all piecemeals sums up, and obtains contaminated site soil remediation amount.The inflection point coordinate in area to be repaired is after simplifying
Area's boundary profile to be repaired on break.
The present embodiment simplifies area's boundary profile to be repaired in conjunction with practice of construction demand, such as the side in area to be repaired
Boundary line is changed to broken line by curve, in favor of practice of construction;By pollution risk distribution map and contaminant transportation depth block diagram into
Row superposition calculates contaminated site soil remediation amount according to the area to be repaired on stacking chart, can more reasonably estimate contaminated site
Soil remediation amount calculates the reparation amount more acurrate.
In another embodiment of the present invention, as shown in Fig. 2, a kind of calculation method of contaminated site soil remediation amount,
Include:
Step S000 respectively contaminated site multiple sampled points carry out soil sample, by the soil sample to acquisition into
The test of row pollutant concentration, obtains soil sample data.
Step S110 is handled the soil sample data of contaminated site according to place geotechnical stratified situation, is obtained every
One soil layer sample data;
Wherein, from the soil sample data of contaminated site, by all soil-likes for belonging to the same sampled point of same soil layer
Notebook data sums up averaging, and gained mean value corresponds to the sample data of sampled point as corresponding soil layer;Same soil layer will be belonged to
The sample datas of all sampled points constitute corresponding soil layer sample data.
Step S120 goes out the number of other unmeasured positions using inverse distance weighting interpolation according to each soil layer sample data
According to obtain the pollutant concentration of each position of each soil layer;
For step S130 according to the pollutant concentration of each position of each soil layer, reference soil background and soil pollutant are dense
Degree limit value classifies to pollution risk, obtains the pollution risk distribution map of corresponding soil layer.
Step S210 according to the pollution risk distribution map, select pollution risk grade reach the region of predetermined level as
Area to be repaired;
Step S220 is based on Douglas-Pu Ke algorithm, takes turns according to maximum allowable offset is preset to area boundary to be repaired
Exterior feature removes extra break, obtains first and simplifies boundary profile;
Step S230 is based on Zhou-Jones algorithm, simplifies boundary profile removal to described first according to default minimum area
Extra break;
It is described preset maximum allowable offset be 2-4 meter, the default minimum area be 20-40 square metres;
Step S300 is according to the soil sample data, using Thiessen polygon to the contaminant transportation depth of different soil
Piecemeal zoning is carried out, contaminant transportation depth block diagram is obtained;
Step S400 carries out the pollution risk distribution map of identical soil layer and the contaminant transportation depth block diagram
Superposition calculates contaminated site soil remediation amount according to the area to be repaired on stacking chart, and the inflection point for obtaining the area to be repaired is sat
Mark is constructed for subsequent soil remediation.
Specifically, Qinghai chemical plant mainly produces chlorate by taking the region of Qinghai as an example, stopped production in 2010.It should
There is residential block except east in region, east is farmland, and there are the lands used such as residential block, school in west side.Plant area's former address area is about
19000m2.Since factory's production technology characteristic and supervision gaps and omissions, chromium residue become the main source of plant area's soil pollution.
It is in place geotechnical engineering investigation report statistics indicate that, from top to bottom, soil layer is divided into three layers: 1. laminin is banketed (thickness 0.7-
0.9 meter) 2. 3. layer cobble (5.2-5.6 meters of thickness), stratum are in samely layer Loessial silty clay (6.8-7.2 meters of thickness)
On looks unit, engineering characteristic is without significant changes.
Collect soil sample in the venue by engineering drilling, amount to 23 sampled points and each sampled point 13 meters with
Interior vertical depth carries out multidraw, for example, the sampling of different depth is carried out with 0.5 meter of spacing, until heavy metal in soil
Concentration limits in safety limit, such as lower than Portable heavy metal analyzer (XRF) detection.The pollutant concentration that will test out is exceeded
Sample send to lab analysis, obtain the pollutant concentration of each sample, these data constitute soil sample data.Soil-like
Notebook data shows contaminant transportation to the second soil layer.
According to place geotechnical stratified situation, soil sample data are handled, the same sampled point of same soil layer will be belonged to
Soil sample data sum up averaging, gained mean value corresponds to the sample data of sampled point as corresponding soil layer;It will belong to
The sample data of all sampled points of same soil layer constitutes corresponding soil layer sample data.
The sample data of each each sampled point of soil layer can be calculated by following formula:
WhereinFor the xth soil layer soil sample pollutant measured data mean value (x=1,2,3......, n),
CiFor same i-th of the sample data of soil layer (i=1,2,3......, n) of same soil sampling point.
Each soil layer sample data is obtained into the data of unmeasured position through inverse distance weighting interpolation.According to sample point
The data of data and unmeasured position obtain the pollutant concentration of each position of same soil layer, are portrayed according to the pollutant concentration
Pollutant layer distributed situation, and referring to soil background and " soil environment quality-construction land soil pollution risk management and control mark
It is quasi- " the second class of pollutant land used control value classifies (as shown in Figure 5) to pollution risk, obtain the pollution of each soil layer
Risk distribution figure.
According to pollution risk distribution map, selecting pollution risk grade is high region as area to be repaired.In conjunction with reality
Construction demand simplifies area's boundary profile to be repaired.First using Douglas-Pu Ke algorithm to area's boundary profile to be repaired
A removal is carried out, the algorithm principle is as follows:
It is connected to the first and last point of each curve and obtains a Trendline, asks all the points on curve vertical with the Trendline
Distance, and maximum range value dmax is found out, with dmax and tolerance D (D is the maximum allowable offset in preset simplified tolerance)
Compare: if dmax < D, the intermediate point on this curve is all cast out;If dmax >=D, the corresponding coordinate points of reservation dmax, and with
The point is boundary, and curve is divided into two parts, to constitute two new trend lines, then measures remaining break again to this two lines
Vertical range, until whole process will continue to that all breaks whole at a distance from Trendline less than tolerance is deleted, thus
It obtains first and simplifies boundary profile.
Simplify boundary profile to first, using Zhou-Jones algorithm according to the default minimal face in preset simplified tolerance
Product removes extra break, area's boundary profile to be repaired after being simplified.It is obtained by Zhou-Jones algorithm every on boundary
Effective triangle area of a break, if effective triangle area of break is less than the default minimum in preset simplified tolerance
Area then removes the break.
By Douglas-Pu Ke algorithm and Zhou-Jones algorithm, area's boundary profile to be repaired is simplified, such as
Boundary profile will be repaired, broken line is changed to by curve, more meets practice of construction demand.
The simplification in area to be repaired is handled by the cartographic generaliztion tool box of ArcGIS10.2, and sets simplified tolerance most
Big allowable offset is 2 meters, minimum area is 20 square metres.Using the element tool box of ArcGIS10.2, will be rolled in simplification figure
Point is extracted as point shape file;Using the regular toolbox of ArcGIS10.2, fold point shape file is handled, and deletes point
Concentrate equitant point element;(x, y) is carried out to treated point shape file using the element tool box of ArcGIS10.2 to sit
Assignment is marked, each inflection point coordinate of simplification figure is obtained.
Piecemeal zoning is carried out to different soil contaminant transportation depth using Thiessen polygon, obtains contaminant transportation depth
Distribution map, wherein the corresponding Thiessen polygon of each delimitation block in the first soil layer, the second soil layer, a kind of corresponding pollutant
Transfer Depth.
The pollution risk distribution map of above-mentioned identical soil layer and contaminant transportation depth block diagram are overlapped, by different soil
The volume of all piecemeals in all areas to be repaired of floor sums up, and obtains contaminated site soil remediation amount.As shown in fig. 6,
Finally obtaining the total side's amount of first layer soil block 10,11,13,15,19,26,31,37,38,40,42 to be repaired is
3631.77m3;The total side's amount of second layer soil block 26,28,32,39,44,53,55,58 to be repaired is 10253.92m3。
The amount of contaminated site soil remediation side described in the present embodiment calculates optimization method, and complicated pollutant distribution is three-dimensional
Model is simplified, and division obtains simple contaminated soil piecemeal, to calculate corresponding soils reparation side amount;Fully consider knot
Practice of construction demand is closed, contaminated soil block inflection point coordinate is obtained, convenient for the design of practical soil remediation construction work
With progress.
In another embodiment of the present invention, as shown in figure 3, a kind of computing system of contaminated site soil remediation amount,
Include:
Pollution risk evaluation module 110 obtains the contaminated site for the soil sample data according to contaminated site
The pollution risk distribution map of different soil;
Area to be repaired simplifies module 120, for selecting pollution risk grade to reach pre- according to the pollution risk distribution map
If the region of grade simplifies area's boundary profile to be repaired as area to be repaired, and according to preset simplified tolerance;
Transfer Depth evaluation module 130 is used for according to the soil sample data, using Thiessen polygon to different soil
Contaminant transportation depth carry out piecemeal zoning, obtain contaminant transportation depth block diagram;
Reparation amount computing module 140, for by the pollution risk distribution map of identical soil layer and the contaminant transportation
Depth block diagram is overlapped, and calculates contaminated site soil remediation amount according to the area to be repaired on stacking chart, and obtain it is described to
The inflection point coordinate for repairing area is constructed for subsequent soil remediation.
Specifically, soil sample data are to carry out experimental analysis by soil sampling on the spot to the soil of sampling and obtain.Soil
Earth sample data reflects pollutant concentration contained by soil.Consider human and material resources cost, the sampling point of soil sampling is limited on the spot
, according to soil sample data, the pollutant concentration of unmeasured point is obtained by space interpolation, to obtain the dirt of different soil
Dye distribution;Definition according to the pollutant concentration of each position of each soil layer and country or professional standard to pollution risk grade,
Obtain the pollution risk distribution map of each soil layer.
For example pollution risk grade is divided into three-level, it is respectively high, medium, low;If predetermined level be it is medium,
Selection pollution risk grade reaches medium or more region, i.e., high, medium region is as area to be repaired;If default etc.
Grade be it is high, then selecting pollution risk grade is high region as area to be repaired.It is influenced less on soil remediation amount
In the case of, area's boundary profile to be repaired is simplified, in favor of the construction of subsequent soil remediation.Douglas-can be used
General gram of algorithm, and/or Wang-M ü ller algorithm, and/or Zhou-Jones algorithm simplify area's boundary profile to be repaired.
Douglas-Pu Ke algorithm retains the key point for constituting the basic configuration of boundary profile, and removes every other point;The algorithm is surveyed
Measure boundary on each break arrive Trendline vertical range, at a distance from Trendline be less than tolerance break will be deleted.Wang-
M ü ller algorithm finds out borderline bending by shape recognition techniques and analyzes its feature, then eliminates inessential folding
It is curved.Zhou-Jones algorithm identifies effective triangle area of each break on boundary, then using a series of measurements to these
Triangle is weighted to compare the flatness of each area, the degree of bias and convexity, is removed by the area guidance after weighting corresponding
Break is to simplify facial contour, while retention characteristic as much as possible.
Each sample has a corresponding Transfer Depth on practical contaminated site, and sample relative discrete, and
The Transfer Depth of each sample is different, so using Thiessen polygon to the contaminant transportation depth of different soil in place
Piecemeal zoning is carried out, rough operation principle and process is as follows:
Thiessen polygon is a kind of subdivision method to investigation range plane, its main feature is that any position in polygon from
The distance of sampled point is nearest in the polygon, remote with a distance from sampled point in adjacent polygons, and each polygon includes and only
Include a sampled point.Thiessen polygon is constructed as follows:
1. according to sequential scan input point (i.e. sampled point) from left to right, from top to bottom.
Meet delaunay criterion 2. marking off in all the points (not comprising other points in i.e. each triangle circumscribed circle)
Irregular triangle network.
3. doing perpendicular bisector to each side of all triangles, these perpendicular bisectors form the side of Thiessen polygon, each to hang down
The intersection point of straight bisector determines the position of Thiessen polygon break.
The pollution risk distribution map of identical soil layer and contaminant transportation depth block diagram are overlapped, each soil layer is obtained
Stacking chart.Area to be repaired on the stacking chart of each soil layer includes several independent regions, and each region includes several points again
Block, the piecemeal in each region are obtained according to Thiessen polygon again, and each piecemeal corresponds to a kind of contaminant transportation depth, according to every
The area and contaminant transportation depth of a piecemeal obtain the volume of piecemeal, by all points in all areas to be repaired of different soil
The volume of block sums up, and obtains contaminated site soil remediation amount.The inflection point coordinate in area to be repaired is simplified area to be repaired
Break on boundary profile.
The present embodiment simplifies area's boundary profile to be repaired in conjunction with practice of construction demand, such as the side in area to be repaired
Boundary line is changed to broken line by curve, in favor of practice of construction;By pollution risk distribution map and contaminant transportation depth block diagram into
Row superposition calculates contaminated site soil remediation amount according to the area to be repaired on stacking chart, can more reasonably estimate contaminated site
Soil remediation amount calculates the reparation amount more acurrate.
In another embodiment of the present invention, as shown in figure 4, a kind of computing system of contaminated site soil remediation amount,
Include:
Sample data obtains module 100, for carrying out soil sample in multiple sampled points of contaminated site respectively, by right
The soil sample of acquisition carries out pollutant concentration test, obtains soil sample data.
Pollution risk evaluation module 110 is used for according to place geotechnical stratified situation, to the soil sample data of contaminated site
It is handled, obtains each soil layer sample data;Wherein, from the soil sample data of contaminated site, by it is all belong to it is same
The soil sample data of the same sampled point of soil layer sum up averaging, and gained mean value corresponds to the sample of sampled point as corresponding soil layer
Notebook data;The sample data that all sampled points of same soil layer will be belonged to constitutes corresponding soil layer sample data;And according to each
Soil layer sample data, the data of other unmeasured positions are gone out using inverse distance weighting interpolation, so that it is each to obtain each soil layer
The pollutant concentration of position;And the pollutant concentration according to each position of each soil layer, it is dirty referring to soil background and soil
Dye object concentration limit classifies to pollution risk, obtains the pollution risk distribution map of corresponding soil layer.
Area to be repaired simplifies module 120, for selecting pollution risk grade to reach pre- according to the pollution risk distribution map
If the region of grade is as area to be repaired;And it is based on Douglas-Pu Ke algorithm, according to presetting maximum allowable offset pair
Area's boundary profile to be repaired removes extra break, obtains first and simplifies boundary profile;Based on Zhou-Jones algorithm, according to pre-
If minimum area, which simplifies boundary profile to described first, removes extra break;Wherein, presetting maximum allowable offset is 2-4
Rice, presetting minimum area is 20-40 square metres.
Transfer Depth evaluation module 130 is used for according to the soil sample data, using Thiessen polygon to different soil
Contaminant transportation depth carry out piecemeal zoning, obtain contaminant transportation depth block diagram;
Reparation amount computing module 140, by the pollution risk distribution map of identical soil layer and the contaminant transportation depth
Block diagram is overlapped, and calculates contaminated site soil remediation amount according to the area to be repaired on stacking chart, and obtain described to be repaired
The inflection point coordinate in area is constructed for subsequent soil remediation.
Specifically, Qinghai chemical plant mainly produces chlorate by taking the region of Qinghai as an example, stopped production in 2010.It should
There is residential block except east in region, east is farmland, and there are the lands used such as residential block, school in west side.Plant area's former address area is about
19000m2.Since factory's production technology characteristic and supervision gaps and omissions, chromium residue become the main source of plant area's soil pollution.
It is in place geotechnical engineering investigation report statistics indicate that, from top to bottom, soil layer is divided into three layers: 1. laminin is banketed (thickness 0.7-
0.9 meter) 2. 3. layer cobble (5.2-5.6 meters of thickness), stratum are in samely layer Loessial silty clay (6.8-7.2 meters of thickness)
On looks unit, engineering characteristic is without significant changes.
Collect soil sample in the venue by engineering drilling, amount to 23 sampled points and each sampled point 13 meters with
Interior vertical depth carries out multidraw, for example, the sampling of different depth is carried out with 0.5 meter of spacing, until heavy metal in soil
Concentration limits in safety limit, such as lower than Portable heavy metal analyzer (XRF) detection.The pollutant concentration that will test out is exceeded
Sample send to lab analysis, obtain the pollutant concentration of each sample, these data constitute soil sample data.Soil-like
Notebook data shows contaminant transportation to the second soil layer.
According to place geotechnical stratified situation, soil sample data are handled, the same sampled point of same soil layer will be belonged to
Soil sample data sum up averaging, gained mean value corresponds to the sample data of sampled point as corresponding soil layer;It will belong to
The sample data of all sampled points of same soil layer constitutes corresponding soil layer sample data.
The sample data of each each sampled point of soil layer can be calculated by following formula:
WhereinFor the xth soil layer soil sample pollutant measured data mean value (x=1,2,3......, n),
CiFor same i-th of the sample data of soil layer (i=1,2,3......, n) of same soil sampling point.
Each soil layer sample data is obtained into the data of unmeasured position through inverse distance weighting interpolation.According to sample point
The data of data and unmeasured position obtain the pollutant concentration of each position of same soil layer, are portrayed according to the pollutant concentration
Pollutant layer distributed situation, and referring to soil background and " soil environment quality-construction land soil pollution risk management and control mark
It is quasi- " the second class of pollutant land used control value classifies (as shown in Figure 5) to pollution risk, obtain the pollution of each soil layer
Risk distribution figure.
According to pollution risk distribution map, selecting pollution risk grade is high region as area to be repaired.In conjunction with reality
Construction demand simplifies area's boundary profile to be repaired.First using Douglas-Pu Ke algorithm to area's boundary profile to be repaired
A removal is carried out, the algorithm principle is as follows:
It is connected to the first and last point of each curve and obtains a Trendline, asks all the points on curve vertical with the Trendline
Distance, and maximum range value dmax is found out, with dmax and tolerance D (D is the maximum allowable offset in preset simplified tolerance)
Compare: if dmax < D, the intermediate point on this curve is all cast out;If dmax >=D, the corresponding coordinate points of reservation dmax, and with
The point is boundary, and curve is divided into two parts, to constitute two new trend lines, then measures remaining break again to this two lines
Vertical range, until whole process will continue to that all breaks whole at a distance from Trendline less than tolerance is deleted, thus
It obtains first and simplifies boundary profile.
Simplify boundary profile to first, using Zhou-Jones algorithm according to the default minimal face in preset simplified tolerance
Product removes extra break, area's boundary profile to be repaired after being simplified.It is obtained by Zhou-Jones algorithm every on boundary
Effective triangle area of a break, if effective triangle area of break is less than the default minimum in preset simplified tolerance
Area then removes the break.
By Douglas-Pu Ke algorithm and Zhou-Jones algorithm, area's boundary profile to be repaired is simplified, such as
Boundary profile will be repaired, broken line is changed to by curve, more meets practice of construction demand.
The simplification in area to be repaired is handled by the cartographic generaliztion tool box of ArcGIS10.2, and sets simplified tolerance most
Big allowable offset is 2 meters, minimum area is 20 square metres.Using the element tool box of ArcGIS10.2, will be rolled in simplification figure
Point is extracted as point shape file;Using the regular toolbox of ArcGIS10.2, fold point shape file is handled, and deletes point
Concentrate equitant point element;(x, y) is carried out to treated point shape file using the element tool box of ArcGIS10.2 to sit
Assignment is marked, each inflection point coordinate of simplification figure is obtained.
Piecemeal zoning is carried out to different soil contaminant transportation depth using Thiessen polygon, obtains contaminant transportation depth
Distribution map, wherein the corresponding Thiessen polygon of each delimitation block in the first soil layer, the second soil layer, a kind of corresponding pollutant
Transfer Depth.
The pollution risk distribution map of above-mentioned identical soil layer and contaminant transportation depth block diagram are overlapped, by different soil
The volume of all piecemeals in all areas to be repaired of floor sums up, and obtains contaminated site soil remediation amount.As shown in fig. 6,
Finally obtaining the total side's amount of first layer soil block 10,11,13,15,19,26,31,37,38,40,42 to be repaired is
3631.77m3;The total side's amount of second layer soil block 26,28,32,39,44,53,55,58 to be repaired is 10253.92m3。
The amount of contaminated site soil remediation side described in the present embodiment calculates optimization method, and complicated pollutant distribution is three-dimensional
Model is simplified, and division obtains simple contaminated soil piecemeal, to calculate corresponding soils reparation side amount;Fully consider knot
Practice of construction demand is closed, contaminated soil block inflection point coordinate is obtained, convenient for the design of practical soil remediation construction work
With progress.
It should be noted that above-described embodiment can be freely combined as needed.The above is only of the invention preferred
Embodiment, it is noted that for those skilled in the art, in the premise for not departing from the principle of the invention
Under, several improvements and modifications can also be made, these modifications and embellishments should also be considered as the scope of protection of the present invention.
Claims (10)
1. a kind of calculation method of contaminated site soil remediation amount characterized by comprising
According to the soil sample data of contaminated site, the pollution risk distribution map of the different soil of the contaminated site is obtained;
According to the pollution risk distribution map, pollution risk grade is selected to reach the region of predetermined level as area to be repaired, and
Area's boundary profile to be repaired is simplified according to preset simplified tolerance;
According to the soil sample data, piecemeal area is carried out using contaminant transportation depth of the Thiessen polygon to different soil
It draws, obtains contaminant transportation depth block diagram;
The pollution risk distribution map of identical soil layer and the contaminant transportation depth block diagram are overlapped, according to superposition
Area to be repaired on figure calculates contaminated site soil remediation amount, and obtains the inflection point coordinate in the area to be repaired, for subsequent
Soil remediation construction.
2. the calculation method of contaminated site soil remediation amount according to claim 1, which is characterized in that described according to pollution
The soil sample data in place, the pollution risk distribution map for obtaining the different soil of the contaminated site specifically include:
According to place geotechnical stratified situation, the soil sample data of contaminated site are handled, obtain each soil layer sample number
According to;
According to each soil layer sample data, the data of other unmeasured positions are gone out using inverse distance weighting interpolation, to obtain
The pollutant concentration of each each position of soil layer;
According to the pollutant concentration of each position of each soil layer, referring to soil background and soil pollutant concentration limit to pollution
Risk is classified, and the pollution risk distribution map of corresponding soil layer is obtained.
3. the calculation method of contaminated site soil remediation amount according to claim 2, which is characterized in that described to Pollution Field
The soil sample data on ground are handled, and are obtained each soil layer sample data and are included:
From the soil sample data of contaminated site, all soil sample data for belonging to the same sampled point of same soil layer are carried out
Adduction is averaging, and gained mean value corresponds to the sample data of sampled point as corresponding soil layer;All samplings of same soil layer will be belonged to
The sample data of point constitutes corresponding soil layer sample data.
4. the calculation method of contaminated site soil remediation amount according to claim 1, which is characterized in that described according to default
Simplification tolerance area's boundary profile to be repaired is simplified, specifically include:
It is extra according to presetting maximum allowable offset and being removed to area's boundary profile to be repaired based on Douglas-Pu Ke algorithm
Break obtains first and simplifies boundary profile;
Based on Zhou-Jones algorithm, simplifies boundary profile to described first according to default minimum area and remove extra break.
5. the calculation method of contaminated site soil remediation amount according to claim 4 characterized by comprising
It is described preset maximum allowable offset be 2-4 meter, the default minimum area be 20-40 square metres.
6. the calculation method of contaminated site soil remediation amount according to claim 1, which is characterized in that described according to dirt
The soil sample data for contaminating place, the pollution risk distribution map for obtaining the different soil of the contaminated site include: before
Soil sample is carried out in multiple sampled points of contaminated site respectively, by carrying out pollutant concentration to the soil sample of acquisition
Test, obtains soil sample data.
7. a kind of computing system of contaminated site soil remediation amount characterized by comprising
Pollution risk evaluation module obtains the difference soil of the contaminated site for the soil sample data according to contaminated site
The pollution risk distribution map of layer;
Area to be repaired simplifies module, for selecting pollution risk grade to reach predetermined level according to the pollution risk distribution map
Region area's boundary profile to be repaired is simplified as area to be repaired, and according to preset simplified tolerance;
Transfer Depth evaluation module, for the pollution according to the soil sample data, using Thiessen polygon to different soil
Object Transfer Depth carries out piecemeal zoning, obtains contaminant transportation depth block diagram;
Reparation amount computing module, by the pollution risk distribution map of identical soil layer and the contaminant transportation depth block diagram into
Row superposition calculates contaminated site soil remediation amount according to the area to be repaired on stacking chart, and obtains the inflection point in the area to be repaired
Coordinate is constructed for subsequent soil remediation.
8. the computing system of contaminated site soil remediation amount according to claim 7, it is characterised in that:
The pollution risk evaluation module is further used for according to place geotechnical stratified situation, to the soil sample of contaminated site
Data are handled, and each soil layer sample data is obtained;According to each soil layer sample data, gone out using inverse distance weighting interpolation
The data of other unmeasured positions, to obtain the pollutant concentration of each position of each soil layer;According to each position of each soil layer
The pollutant concentration set is classified to pollution risk referring to soil background and soil pollutant concentration limit, is obtained
The pollution risk distribution map of corresponding soil layer.
9. the computing system of contaminated site soil remediation amount according to claim 7, it is characterised in that:
The area to be repaired simplifies module, is further used for based on Douglas-Pu Ke algorithm, according to presetting maximum allowable offset
Amount removes extra break to area's boundary profile to be repaired, obtains first and simplifies boundary profile;Based on Zhou-Jones algorithm, press
Simplify boundary profile to described first according to default minimum area and removes extra break.
10. the computing system of contaminated site soil remediation amount according to claim 7, which is characterized in that further include:
Sample data obtains module, for carrying out soil sample in multiple sampled points of contaminated site respectively, by acquisition
Soil sample carries out pollutant concentration test, obtains soil sample data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811648129.8A CN109754182B (en) | 2018-12-29 | 2018-12-29 | Method and system for calculating soil remediation quantity of contaminated site |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811648129.8A CN109754182B (en) | 2018-12-29 | 2018-12-29 | Method and system for calculating soil remediation quantity of contaminated site |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109754182A true CN109754182A (en) | 2019-05-14 |
CN109754182B CN109754182B (en) | 2020-09-11 |
Family
ID=66404952
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811648129.8A Active CN109754182B (en) | 2018-12-29 | 2018-12-29 | Method and system for calculating soil remediation quantity of contaminated site |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109754182B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110826236A (en) * | 2019-11-11 | 2020-02-21 | 北京建工环境修复股份有限公司 | Method for judging soil remediation pollution range |
CN111008920A (en) * | 2019-11-07 | 2020-04-14 | 上海勘察设计研究院(集团)有限公司 | Polluted site investigation method based on groundwater level fluctuation effect |
CN111400660A (en) * | 2020-03-05 | 2020-07-10 | 北京市环境保护科学研究院 | Method for calculating soil restoration target value of construction land based on trial algorithm |
CN112308290A (en) * | 2020-09-30 | 2021-02-02 | 湘潭大学 | Method for investigating and stationing polluted site |
CN113112057A (en) * | 2021-03-29 | 2021-07-13 | 广东省建筑工程监理有限公司 | Method for managing polluted site and repairing earth volume by combining unmanned aerial vehicle aerial survey and BIM technology |
CN114023042A (en) * | 2021-12-17 | 2022-02-08 | 中煤科工集团沈阳研究院有限公司 | Life safety monitoring and early warning system and method |
CN117196156A (en) * | 2023-10-19 | 2023-12-08 | 北京建工环境修复股份有限公司 | Intelligent planning method, system and medium for pollution in-situ treatment |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090094181A1 (en) * | 2007-09-28 | 2009-04-09 | Korean Electrical Power Corperation | Method of estimating salt contamination |
CN105117610A (en) * | 2015-09-21 | 2015-12-02 | 黄波 | Method for estimating fine particulate matter concentration in real time based on temporal and spatial characteristics |
CN105678454A (en) * | 2016-01-04 | 2016-06-15 | 北京师范大学 | Method for determining soil restoration amount of polluted field on the basis of condition simulation |
CN106446435A (en) * | 2016-10-09 | 2017-02-22 | 北京建工环境修复股份有限公司 | Method for optimizing contaminated soil restoration scope and boundary |
CN106528788A (en) * | 2016-11-09 | 2017-03-22 | 上海市政工程设计研究总院(集团)有限公司 | Method for analyzing space distribution feature of ground rainfall runoff pollution based on GIS (Geographic Information System) technology |
CN107705002A (en) * | 2017-09-21 | 2018-02-16 | 中国矿业大学(北京) | The determination method of mining soil content of beary metal sampled point exception high level coverage |
CN108268977A (en) * | 2018-01-12 | 2018-07-10 | 武汉智博创享科技股份有限公司 | Reduce the method and system that the extremely extra-high value of pollutant in soil remediation influences |
-
2018
- 2018-12-29 CN CN201811648129.8A patent/CN109754182B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090094181A1 (en) * | 2007-09-28 | 2009-04-09 | Korean Electrical Power Corperation | Method of estimating salt contamination |
CN105117610A (en) * | 2015-09-21 | 2015-12-02 | 黄波 | Method for estimating fine particulate matter concentration in real time based on temporal and spatial characteristics |
CN105678454A (en) * | 2016-01-04 | 2016-06-15 | 北京师范大学 | Method for determining soil restoration amount of polluted field on the basis of condition simulation |
CN106446435A (en) * | 2016-10-09 | 2017-02-22 | 北京建工环境修复股份有限公司 | Method for optimizing contaminated soil restoration scope and boundary |
CN106528788A (en) * | 2016-11-09 | 2017-03-22 | 上海市政工程设计研究总院(集团)有限公司 | Method for analyzing space distribution feature of ground rainfall runoff pollution based on GIS (Geographic Information System) technology |
CN107705002A (en) * | 2017-09-21 | 2018-02-16 | 中国矿业大学(北京) | The determination method of mining soil content of beary metal sampled point exception high level coverage |
CN108268977A (en) * | 2018-01-12 | 2018-07-10 | 武汉智博创享科技股份有限公司 | Reduce the method and system that the extremely extra-high value of pollutant in soil remediation influences |
Non-Patent Citations (3)
Title |
---|
刘庚 等: "某铅酸蓄电池污染场地表层土壤重金属Pb空间分布预测研究", 《环境科学》 * |
孙晓松 等: "典型铬污染场地健康风险评估及污染土壤的定量计算", 《2011污染场地修复产业国际论坛暨重庆市环境科学学会第九届学术年会》 * |
蒋世杰 等: "基于条件模拟的污染场地土壤修复量的确定研究", 《环境科学学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111008920A (en) * | 2019-11-07 | 2020-04-14 | 上海勘察设计研究院(集团)有限公司 | Polluted site investigation method based on groundwater level fluctuation effect |
CN111008920B (en) * | 2019-11-07 | 2023-05-16 | 上海勘察设计研究院(集团)有限公司 | Pollution site investigation method based on underground water level fluctuation effect |
CN110826236B (en) * | 2019-11-11 | 2023-11-21 | 北京建工环境修复股份有限公司 | Determination method for soil remediation pollution range |
CN110826236A (en) * | 2019-11-11 | 2020-02-21 | 北京建工环境修复股份有限公司 | Method for judging soil remediation pollution range |
CN111400660A (en) * | 2020-03-05 | 2020-07-10 | 北京市环境保护科学研究院 | Method for calculating soil restoration target value of construction land based on trial algorithm |
CN111400660B (en) * | 2020-03-05 | 2021-08-27 | 北京市环境保护科学研究院 | Method for calculating soil restoration target value of construction land based on trial algorithm |
CN112308290B (en) * | 2020-09-30 | 2023-07-28 | 湘潭大学 | Method for investigating and distributing points in polluted site |
CN112308290A (en) * | 2020-09-30 | 2021-02-02 | 湘潭大学 | Method for investigating and stationing polluted site |
CN113112057B (en) * | 2021-03-29 | 2023-08-29 | 广东省建筑工程监理有限公司 | Method for managing soil volume restoration of polluted site by combining unmanned aerial vehicle aerial survey and BIM technology |
CN113112057A (en) * | 2021-03-29 | 2021-07-13 | 广东省建筑工程监理有限公司 | Method for managing polluted site and repairing earth volume by combining unmanned aerial vehicle aerial survey and BIM technology |
CN114023042A (en) * | 2021-12-17 | 2022-02-08 | 中煤科工集团沈阳研究院有限公司 | Life safety monitoring and early warning system and method |
CN117196156A (en) * | 2023-10-19 | 2023-12-08 | 北京建工环境修复股份有限公司 | Intelligent planning method, system and medium for pollution in-situ treatment |
CN117196156B (en) * | 2023-10-19 | 2024-02-06 | 北京建工环境修复股份有限公司 | Intelligent planning method, system and medium for pollution in-situ treatment |
Also Published As
Publication number | Publication date |
---|---|
CN109754182B (en) | 2020-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109754182A (en) | A kind of calculation method and system of contaminated site soil remediation amount | |
CN110009226B (en) | Site pollution risk level evaluation and sensitive receptor influence identification system | |
Burrough et al. | Continuous classification in soil survey: spatial correlation, confusion and boundaries | |
CN112417573B (en) | GA-LSSVM and NSGA-II shield tunneling multi-objective optimization method based on existing tunnel construction | |
CN106446435B (en) | A kind of optimization method on contaminated soil remediation range and boundary | |
CN111260208A (en) | Regional gridding accumulative environment risk assessment system and method based on risk field | |
CN110852557A (en) | GIS-based refuse landfill site selection method and device | |
CN103984997A (en) | Transmission project site and line selecting method based on GIS space information | |
CN110597932B (en) | Environment comprehensive evaluation prediction method based on remote sensing image | |
CN106570334A (en) | Method for calculating the pollution load of urban rainfall runoff by geostatistical model | |
CN110570031A (en) | Method for predicting river channel split ratio of two-stage branched river network in plain | |
CN113720639B (en) | Method for investigating and supplementing sampling of polluted site | |
CN110059870B (en) | BIM and GIS based navigation channel renovation building maintenance analysis method | |
CN108647826A (en) | A kind of construction method and device of heavy metal-polluted soil environmental risk prediction model | |
CN111882243A (en) | Karst region territory space development suitability evaluation method and system | |
Li et al. | Wind tunnel study on the morphological parameterization of building non-uniformity | |
CN116401327A (en) | Storm flood calculation auxiliary system for small and medium-sized watershed design in non-data area | |
Mitasova et al. | Terrain modeling and soil erosion simulations for Fort Hood and Fort Polk test areas | |
CN117348092A (en) | Groundwater pollutant tracing method based on aggregate smoothing algorithm and geophysical prospecting data | |
CN110717233A (en) | Rainwater pipe network flow calculation method and system based on GIS underlying surface analysis | |
CN115526544A (en) | Green mine construction planning method and system based on multi-type data | |
CN115239097A (en) | Photo-thermal base planning method based on environmental impact evaluation key technology index system | |
Yao | Application of GIS remote sensing information integration in eco-environmental quality monitoring | |
Xiao et al. | Landscape architecture construction using landscape urbanism and dynamic nonlinear system thinking | |
Li et al. | Study on pollutant model construction and three-dimensional spatial interpolation in soil environmental survey |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |