CN118153279A - 一种焦化企业地块土壤污染范围划定的方法 - Google Patents
一种焦化企业地块土壤污染范围划定的方法 Download PDFInfo
- Publication number
- CN118153279A CN118153279A CN202410140339.5A CN202410140339A CN118153279A CN 118153279 A CN118153279 A CN 118153279A CN 202410140339 A CN202410140339 A CN 202410140339A CN 118153279 A CN118153279 A CN 118153279A
- Authority
- CN
- China
- Prior art keywords
- range
- point
- pollution
- partitioning
- sampling points
- 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.)
- Pending
Links
- 238000004939 coking Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000003900 soil pollution Methods 0.000 title claims abstract description 27
- 238000005070 sampling Methods 0.000 claims abstract description 60
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 9
- 231100000719 pollutant Toxicity 0.000 claims abstract description 9
- 238000013316 zoning Methods 0.000 claims abstract 3
- 238000000638 solvent extraction Methods 0.000 claims description 39
- 238000009826 distribution Methods 0.000 claims description 36
- 238000004088 simulation Methods 0.000 claims description 24
- 238000005192 partition Methods 0.000 claims description 20
- 238000011835 investigation Methods 0.000 claims description 19
- 238000004519 manufacturing process Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 12
- 239000003802 soil pollutant Substances 0.000 claims description 12
- 239000002689 soil Substances 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 8
- 239000000126 substance Substances 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 5
- 150000001555 benzenes Chemical class 0.000 claims description 4
- 239000000356 contaminant Substances 0.000 claims description 4
- 238000011109 contamination Methods 0.000 claims description 4
- 229910001385 heavy metal Inorganic materials 0.000 claims description 4
- 229930195733 hydrocarbon Natural products 0.000 claims description 4
- 150000002430 hydrocarbons Chemical class 0.000 claims description 4
- 238000011068 loading method Methods 0.000 claims description 4
- 239000003208 petroleum Substances 0.000 claims description 4
- 150000002989 phenols Chemical class 0.000 claims description 4
- 125000005575 polycyclic aromatic hydrocarbon group Chemical group 0.000 claims description 4
- 239000002910 solid waste Substances 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 4
- XFXPMWWXUTWYJX-UHFFFAOYSA-N Cyanide Chemical compound N#[C-] XFXPMWWXUTWYJX-UHFFFAOYSA-N 0.000 claims description 3
- 239000002253 acid Substances 0.000 claims description 3
- 150000007513 acids Chemical class 0.000 claims description 3
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000003912 environmental pollution Methods 0.000 claims description 3
- 239000000796 flavoring agent Substances 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 239000002920 hazardous waste Substances 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 239000000047 product Substances 0.000 claims description 3
- 231100000331 toxic Toxicity 0.000 claims description 3
- 230000002588 toxic effect Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 5
- 239000003245 coal Substances 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 201000004569 Blindness Diseases 0.000 description 2
- 239000000571 coke Substances 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000003895 groundwater pollution Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000000746 purification Methods 0.000 description 2
- 238000004065 wastewater treatment Methods 0.000 description 2
- 229910017052 cobalt Inorganic materials 0.000 description 1
- 239000010941 cobalt Substances 0.000 description 1
- GUTLYIVDDKVIGB-UHFFFAOYSA-N cobalt atom Chemical compound [Co] GUTLYIVDDKVIGB-UHFFFAOYSA-N 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 150000002825 nitriles Chemical class 0.000 description 1
- AXFBAIOSECPASO-UHFFFAOYSA-N pentacyclo[6.6.2.02,7.04,16.011,15]hexadeca-1(14),2(7),3,5,8(16),9,11(15),12-octaene Chemical compound C1=C(C=C23)C4=C5C3=CC=CC5=CC=C4C2=C1 AXFBAIOSECPASO-UHFFFAOYSA-N 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 229910052720 vanadium Inorganic materials 0.000 description 1
- GPPXJZIENCGNKB-UHFFFAOYSA-N vanadium Chemical compound [V]#[V] GPPXJZIENCGNKB-UHFFFAOYSA-N 0.000 description 1
- 239000002912 waste gas Substances 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种焦化企业地块土壤污染范围划定的方法，通过污染物识别获得疑似污染区域，针对疑似污染区域进行布设初步采样点位，初步点位布设完成后，利用ArcGIS软件及克里格插值模拟土壤污染浓度来布设的详细采样点位，将初步采样点位和详细采样点位合并，利用ArcGIS软件分析得到超标点位分区范围，然后利用克里格插值模拟空白区域浓度，得到模拟浓度超标范围，最后根据超标点位分区范围与模拟浓度超标范围对污染范围进行修订，得到目标焦化企业地块的污染范围。本发明解决了现有土壤污染范围划定方法中点位布设及污染范围划定不准确的问题。
Description
技术领域
本发明属于土壤污染分析领域，具体涉及一种焦化企业地块土壤污染范围划定的方法。
背景技术
依据土壤污染防治的要求，关闭搬迁的焦化企业地块需要进行土壤污染状况调查活动，确定土壤污染范围。建设用地土壤污染调查活动主要依据一系列建设用地土壤污染状况调查等技术导则，进行污染识别、点位布设，这一系列技术导则为企业用地调查提供了重要的指导作用。
当前地块土壤污染调查中，常用的布点方法有判断布点法、随机布点法、分区布点法和***布点法。地块内布点主要依据调查目标及污染物类型确定采样位置、采样数量、采样深度等；地块内污染范围的划定主要利用周边不超标点位连线法，但对超标点位周边范围点位布设的具体方法未作详细指导介绍，在调查活动的实际操作中，存在很大的主观性、盲目性；污染范围划定存在污染遗漏或范围扩大的情况。
发明内容
为了克服现有土壤污染范围划定方法中点位布设及污染范围划定存在不准确的缺陷，本发明提供了一种对地块调查布点可以进行定量指导，减少了布点的盲目性，并优化污染范围划定的一种焦化企业地块土壤污染范围划定的方法。
本发明为了实现上述目的所采用的技术方案是：
一种焦化企业地块土壤污染范围划定的方法，包括以下步骤：
步骤一、污染识别：通过收集目标焦化企业地块资料、人员访谈、现场踏勘识别地块主要污染物，并依据主要污染物所在区域划定疑似污染区域；
步骤二、初步布点采样：根据疑似污染区域，利用***随机与专业判断布点法对目标焦化企业地块进行布置初步采样点位，然后采集初步采样点位的土壤样品，并将土壤样品送至实验室进行检测分析，获取初步采样点位的土壤污染物浓度数据，并将初步采样点位划分超标点位和未超标点位；
步骤三、详细布点采样：依据初步采样点位空间分布，利用ArcGIS软件进行泰森多边形分区并计算分区面积，每个分区内有一个点，根据土壤污染物浓度数据，在ArcGIS软件中进行克里格插值模拟目标焦化企业地块内土壤污染物浓度分布，超标点位所在分区内布点按照每1600m2布设1个点位；未超标点位所在分区内按照每6400m2布设1个点位，同时选择在插值模拟浓度超标的范围内布设点位，得到目标焦化企业地块布设的详细采样点位；
步骤四、污染范围划定：先导入目标焦化企业地块布设的初步采样点位和详细采样点位，并利用ArcGIS软件对所有点位进行泰森多边形分区，得到目标焦化企业地块进行划定分区范围，即超标点位分区范围和未超标点位分区范围，然后利用克里格插值模拟空白区域浓度，得到模拟浓度超标范围，最后根据超标点位分区范围与模拟浓度超标范围对污染范围进行修订，得到目标焦化企业地块的污染范围。
进一步地，所述主要污染物包括：多环芳烃、重金属、氰化物、酚类、酸、碱、氨氮、苯系物、石油烃；疑似污染区域包括：(1)根据已有资料或前期调查表明可能存在污染的区域；(2)曾发生泄露或环境污染事故的区域；(3)各类地下罐槽、管线、集水井、检查井等所在的区域；(4)固体废物堆放或填埋的区域；(5)原辅材料、产品、化学品、有毒有害物质以及危险废物生产、贮存、装卸、使用和处置的区域；(6)其他存在明显污染痕迹或存在异味的区域。
优选的，所述***随机与专业判断布点法依据HJ 25.1-2019中规定的布点方法执行。
进一步地，步骤三中所述利用ArcGIS软件进行泰森多边形分区的步骤是：先将初步采样点位数据导入软件，包括初步采样点位的空间信息及初步采样点位的土壤污染物浓度数据；然后利用ArcGIS分析工具，创建泰森多边形，选择初步采样点位的空间信息为分区要素。
进一步地，所述根据超标点位分区范围与模拟浓度超标范围对污染范围进行修订的方法：①超标点位分区范围＞模拟浓度超标范围，按照超标点位分区范围划定污染范围；②超标点位分区范围＜模拟浓度超标范围，按照超标点位分区范围加模拟浓度超标范围划入污染范围：由于调查数据的有限性和克里格插值模拟的平滑效应，易将超标点位被平滑为未超标点而导致修复范围的缩小，此时需将平滑的超标点位分区划入污染范围。
进一步地，所述克里格插值的步骤如下：
先进行半变异函数的计算，计算公示如式(1)所示：
式中：r(h3)为半变异函数；
h3为两采样点位之间的空间距离；
N(h3)是以h3为间距所有样本对数；
Z(xi)和Z(xi+h3)分别是xi和xi+h3处的值；
然后进行克里格插值计算，计算公式如式(2)：
其中，Z0为区域化变量在(X0,Y0)位置的预测值；λi为第i个已知点的权重系数，Zi为第i个已知点的实测值；
克里格方差公式如式(3)：
引入拉格朗日乘子构造拉格朗日函数如式(4)：
对权值与朗格拉日乘子求一阶偏微分如式(5)：
偏微分求解结果如式(6)：
将上述求和部分展开如式(7)：
将上式进一步写为矩阵相乘形式，化为式(8)：
其中A代表在原有变异函数矩阵中额外添加全1行与全1列即交界处1换为0后的矩阵，λ代表个权重组成的列向量，代表拉格朗日乘子，B为各位置与待求解位置对应距离的变异函数值组成的列向量，且在列尾增加一个1；通过矩阵求逆，求解方程组即可得到待求解位置与其它已知点的权重；最后对每一个待插点进行同样操作，完成克里格插值。
步骤四将地块内所有点位数据导入ArcGIS，重复步骤三中泰森多变形分区和克里格插值步骤；划定的土壤污染范围＝分区范围+模拟浓度范围。
本发明通过污染物识别获得疑似污染区域，针对疑似污染区域进行布设初步采样点位，并将初步采样点位划分为超标点位及未超标点位，根据超标点位及未超标点位，利用ArcGIS软件及克里格插值模拟土壤污染浓度来布设目标焦化企业地块布设的详细采样点位，将初步采样点位和详细采样点位合并，利用ArcGIS软件分析得到超标点位分区范围，然后利用克里格插值模拟空白区域浓度，得到模拟浓度超标范围，最后根据超标点位分区范围与模拟浓度超标范围对污染范围进行修订，得到目标焦化企业地块的污染范围。本发明通过对地块调查布点进行更针对性指导，对地块详细调查布点及污染范围的划定提供了定量方法。
本发明运用泰森多边形法和克里格插值法对土壤污染范围进行了划定，结合了两个方法的优点，既从几何原理的角度较为准确地划定了土壤的修复范围，实现了不同深度超标点位修复范围的一致性，并通过克里格插值模拟空白区浓度，降低了布点采样的调查成本，同时最大可能避免污染范围遗漏的问题。
附图说明
下面附图对本发明作进一步描述，其中：
图1为本发明流程图；
图2为山西某焦化厂地块平面布置图；
图3为企业生产工艺流程及产污环节流程图；
图4为初步调查点位布设图；
图5为初步调查污染分区及插值模拟图；
图6为详细调查点位布设图；
图7至14为ArcGIS软件操作步骤图；
图15为地块调查点位泰森多边形分区与克里格插值模拟结果图；
图16为地块污染范围划定结果图。
附图标记说明：1、煤气净化车间，2、酚氰废水处理车间,3、炼焦车间,4、堆煤场。
具体实施方式
如图1所示，本实施例的焦化企业地块土壤污染范围划定的方法，包括以下步骤：
步骤一、污染识别：通过收集目标焦化企业地块资料、人员访谈、现场踏勘识别地块主要污染物，并依据主要污染物所在区域划定疑似污染区域；
步骤二、初步布点采样：根据疑似污染区域，利用***随机与专业判断布点法对目标焦化企业地块进行布置初步采样点位，然后采集初步采样点位的土壤样品，并将土壤样品送至实验室进行检测分析，获取初步采样点位的土壤污染物浓度数据，并将初步采样点位划分超标点位和未超标点位；
步骤三、详细布点采样：依据初步采样点位空间分布，利用ArcGIS软件进行泰森多边形分区并计算分区面积，每个分区内有一个点，根据土壤污染物浓度数据，在ArcGIS软件中进行克里格插值模拟目标焦化企业地块内土壤污染物浓度分布，超标点位所在分区内布点按照每1600m2布设1个点位；未超标点位所在分区内按照每6400m2布设1个点位，同时选择在插值模拟浓度超标的范围内，结合区域设施设备分布情况布设点位，得到目标焦化企业地块布设的详细采样点位。
步骤四、污染范围划定：先导入目标焦化企业地块布设的初步采样点位和详细采样点位，并利用ArcGIS软件对所有点位进行泰森多边形分区，得到目标焦化企业地块进行划定分区范围，即超标点位分区范围和未超标点位分区范围，然后利用克里格插值模拟空白区域浓度，得到模拟浓度超标范围，最后根据超标点位分区范围与模拟浓度超标范围对污染范围进行修订，得到目标焦化企业地块的污染范围。
进一步地，所述主要污染物包括：多环芳烃、重金属、氰化物、酚类、酸、碱、氨氮、苯系物、石油烃；疑似污染区域包括：(1)根据已有资料或前期调查表明可能存在污染的区域；(2)曾发生泄露或环境污染事故的区域；(3)各类地下罐槽、管线、集水井、检查井等所在的区域；(4)固体废物堆放或填埋的区域；(5)原辅材料、产品、化学品、有毒有害物质以及危险废物生产、贮存、装卸、使用和处置的区域；(6)其他存在明显污染痕迹或存在异味的区域。
优选的，所述***随机与专业判断布点法依据HJ 25.1-2019中规定的布点方法执行。
进一步地，步骤三中所述利用ArcGIS软件进行泰森多边形分区的步骤是：先将初步采样点位数据导入软件，包括初步采样点位的空间信息及初步采样点位的土壤污染物浓度数据；然后利用ArcGIS分析工具，创建泰森多边形，选择初步采样点位的空间信息为分区要素。
进一步地，所述根据超标点位分区范围与模拟浓度超标范围对污染范围进行修订的方法：①超标点位分区范围＞模拟浓度超标范围，按照超标点位分区范围划定污染范围；②超标点位分区范围＜模拟浓度超标范围，按照超标点位分区范围加模拟浓度超标范围划入污染范围：由于调查数据的有限性和克里格插值模拟的平滑效应，易将超标点位被平滑为未超标点而导致修复范围的缩小，此时需将平滑的超标点位分区划入污染范围。
进一步地，所述克里格插值的步骤如下：
先进行半变异函数的计算，计算公示如式(1)所示：
式中：r(h3)为半变异函数；
h3为两采样点位之间的空间距离；
N(h3)是以h3为间距所有样本对数；
Z(xi)和Z(xi+h3)分别是xi和xi+h3处的值；
然后进行克里格插值计算，计算公式如式(2)：
其中，Z0为区域化变量在(X0,Y0)位置的预测值；λi为第i个已知点的权重系数，Zi为第i个已知点的实测值；
克里格方差公式如式(3)：
引入拉格朗日乘子构造拉格朗日函数如式(4)：
对权值与朗格拉日乘子求一阶偏微分如式(5)：
偏微分求解结果如式(6)：
将上述求和部分展开如式(7)：
将上式进一步写为矩阵相乘形式，化为式(8)：
其中A代表在原有变异函数矩阵中额外添加全1行与全1列即交界处1换为0后的矩阵，λ代表个权重组成的列向量，代表拉格朗日乘子，B为各位置与待求解位置对应距离的变异函数值组成的列向量，且在列尾增加一个1；通过矩阵求逆，求解方程组即可得到待求解位置与其它已知点的权重；最后对每一个待插点进行同样操作，完成克里格插值。
具体案例：
本次以山西某焦化厂地块，表层土壤苯并(α)芘污染范围划定为例：
步骤一：
1.资料收集、现场踏勘、人员访谈
资料收集：焦化厂平面布置图如图2所示、焦化厂环评、厂区地勘报告、地块规划；现场踏勘：厂区范围，主要厂区生产工艺流程、主要生产、贮存，输送设施情况，废气、污水及固废处理设施情况。场地主要地下管线走向、危化品装卸点与转运路线、可疑污染痕迹、可疑堆弃物等分布情况；人员访谈：通过电话、座谈、邮件等方式，访谈与场地相关人员，核实地块内历史沿革、地块周边环境情况、是否发生过土壤、地下水污染事件等信息。
焦化厂设计生产能力为年产干焦93.6万吨，2001年12月23日开工建设，2004年4月30日投产，2012年8月22日关停，生产经营时间约为8年。焦化厂地块面积为16.87万m2(253.05亩)，主要生产车间分别是备煤车间、炼焦车间、净化车间等；生产经营期间未发生过土壤、地下水污染事件。
企业地块涉及主要构筑物或设备及物料情况见表1。
表1企业地块涉及主要设备及物料情况
企业涉及主要生产工艺及产污环节，见图3。
2.划定潜在污染区域
由上述资料分析可知，本地块内涉及潜在特征污染物为：多环芳烃、氰化物、重金属(铅、汞、钒、钴)、苯系物、石油烃、酚类；潜在污染区域：煤气净化车间1、酚氰废水处理车间2、炼焦车间3、堆煤场4。
步骤二：初步采样点位布设，如图4所示。
步骤三：详细布点采样
将初步调查点位数据导入ArcGIS中，利用ArcGIS软件进行泰森多边形分区，利用ArcGIS计算分区面积，按照超标点位所在分区内布点按照每1600m2布设1个点；未超标所在分区内按照6400m2布设1个点，同时选择在插值模拟浓度超标的范围布点。初步调查污染分区及插值模拟图如图5所示；详细调查点位布设图如图6所示；
ArcGIS软件的操作过程如图7至14所示。
步骤四：污染范围划定
将初步调查与详细调查点位合并，并重复步骤三操作，对地块内所有点位进行泰森多边形分区，并利用克里格插值模拟空白区浓度，结果如图15所示。
污染范围划定：先按照分区范围将超标分区划入污染范围，然后将克里格插值模拟浓度超标范围补充到划入范围，划定结果如图16所示。
Claims (6)
1.一种焦化企业地块土壤污染范围划定的方法，其特征在于，包括以下步骤：
步骤一、污染识别：通过收集目标焦化企业地块资料、人员访谈、现场踏勘识别地块主要污染物，并依据主要污染物所在区域划定疑似污染区域；
步骤二、初步布点采样：根据疑似污染区域，利用***随机与专业判断布点法对目标焦化企业地块进行布置初步采样点位，然后采集初步采样点位的土壤样品，并将土壤样品送至实验室进行检测分析，获取初步采样点位的土壤污染物浓度数据，并将初步采样点位划分超标点位和未超标点位；
步骤三、详细布点采样：依据初步采样点位空间分布，利用ArcGIS软件进行泰森多边形分区并计算分区面积，每个分区内有一个点，根据土壤污染物浓度数据，在ArcGIS软件中进行克里格插值模拟目标焦化企业地块内土壤污染物浓度分布，超标点位所在分区内布点按照每1600m2布设1个点位；未超标点位所在分区内按照每6400m2布设1个点位，同时选择在插值模拟浓度超标的范围内布设点位，得到目标焦化企业地块布设的详细采样点位；
步骤四、污染范围划定：先导入目标焦化企业地块布设的初步采样点位和详细采样点位，并利用ArcGIS软件对所有点位进行泰森多边形分区，得到目标焦化企业地块进行划定分区范围，即超标点位分区范围和未超标点位分区范围，然后利用克里格插值模拟空白区域浓度，得到模拟浓度超标范围，最后根据超标点位分区范围与模拟浓度超标范围对污染范围进行修订，得到目标焦化企业地块的污染范围。
2.根据权利要求1所述的焦化企业地块土壤污染范围划定方法，其特征在于，所述主要污染物包括：多环芳烃、重金属、氰化物、酚类、酸、碱、氨氮、苯系物、石油烃；疑似污染区域包括：(1)根据已有资料或前期调查表明可能存在污染的区域；(2)曾发生泄露或环境污染事故的区域；(3)各类地下罐槽、管线、集水井、检查井等所在的区域；(4)固体废物堆放或填埋的区域；(5)原辅材料、产品、化学品、有毒有害物质以及危险废物生产、贮存、装卸、使用和处置的区域；(6)其他存在明显污染痕迹或存在异味的区域。
3.根据权利要求1所述的焦化企业地块土壤污染范围划定方法，其特征在于，所述***随机与专业判断布点法依据HJ 25.1-2019中规定的布点方法执行。
4.根据权利要求1所述的焦化企业地块土壤污染范围划定方法，其特征在于，步骤三中所述利用ArcGIS软件进行泰森多边形分区的步骤是：先将初步采样点位数据导入软件，包括初步采样点位的空间信息及初步采样点位的土壤污染物浓度数据；然后利用ArcGIS分析工具，创建泰森多边形，选择初步采样点位的空间信息为分区要素。
5.根据权利要求1所述的焦化企业地块土壤污染范围划定方法，其特征在于，所述根据超标点位分区范围与模拟浓度超标范围对污染范围进行修订的方法：①超标点位分区范围＞模拟浓度超标范围，按照超标点位分区范围划定污染范围；②超标点位分区范围＜模拟浓度超标范围，按照超标点位分区范围加模拟浓度超标范围划入污染范围：由于调查数据的有限性和克里格插值模拟的平滑效应，易将超标点位被平滑为未超标点而导致修复范围的缩小，此时需将平滑的超标点位分区划入污染范围。
6.根据权利要求1或5所述的焦化企业地块土壤污染范围划定的方法，其特征在于，所述克里格插值的步骤如下：
先进行半变异函数的计算，计算公示如式(1)所示：
式中：r(h3)为半变异函数；
h3为两采样点位之间的空间距离；
N(h3)是以h3为间距所有样本对数；
Z(xi)和Z(xi+h3)分别是xi和xi+h3处的值；
然后进行克里格插值计算，计算公式如式(2)：
其中，Z0为区域化变量在(X0,Y0)位置的预测值；λi为第i
个已知点的权重系数，Zi为第i个已知点的实测值；
克里格方差公式如式(3)：
引入拉格朗日乘子构造拉格朗日函数如式(4)：
对权值与朗格拉日乘子求一阶偏微分如式(5)：
偏微分求解结果如式(6)：
将上述求和部分展开如式(7)：
将上式进一步写为矩阵相乘形式，化为式(8)：
其中A代表在原有变异函数矩阵中额外添加全1行与全1列即交界处1换为0后的矩阵，λ代表个权重组成的列向量，代表拉格朗日乘子，B为各位置与待求解位置对应距离的变异函数值组成的列向量，且在列尾增加一个1；通过矩阵求逆，求解方程组即可得到待求解位置与其它已知点的权重；最后对每一个待插点进行同样操作，完成克里格插值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410140339.5A CN118153279A (zh) | 2024-02-01 | 2024-02-01 | 一种焦化企业地块土壤污染范围划定的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410140339.5A CN118153279A (zh) | 2024-02-01 | 2024-02-01 | 一种焦化企业地块土壤污染范围划定的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN118153279A true CN118153279A (zh) | 2024-06-07 |
Family
ID=91293908
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410140339.5A Pending CN118153279A (zh) | 2024-02-01 | 2024-02-01 | 一种焦化企业地块土壤污染范围划定的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118153279A (zh) |
-
2024
- 2024-02-01 CN CN202410140339.5A patent/CN118153279A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Saffari et al. | Environmental impact assessment (EIA) by using the Fuzzy Delphi Folchi (FDF) method (case study: Shahrood cement plant, Iran) | |
| Dong et al. | Quantitative evaluation and case studies of cleaner mining with multiple indexes considering uncertainty factors for phosphorus mines | |
| CN114544913B (zh) | 一种土壤调查加密采样布点优化方法 | |
| Marchant et al. | The assessment of point‐source and diffuse soil metal pollution using robust geostatistical methods: a case study in Swansea (Wales, UK) | |
| CN110826236A (zh) | 一种土壤修复污染范围的判定方法 | |
| Tolmachev et al. | Experience in collapse risk assessment of building on covered karst landscapes in Russia | |
| Pashkevich et al. | Development of an operational environmental monitoring system for hazardous industrial facilities of Gazprom Dobycha Urengoy | |
| Doulati Ardejani et al. | Developing a conceptual framework of green mining strategy in coal mines: integrating socio-economic, health, and environmental factors | |
| Campbell | Sensor Systems for Environmental Monitoring: Volume Two: Environmental Monitoring | |
| Sun et al. | Contamination and source-specific health risk assessment of polycyclic aromatic hydrocarbons in soil from a mega iron and steel site in China | |
| CN118153279A (zh) | 一种焦化企业地块土壤污染范围划定的方法 | |
| CN109739898A (zh) | 基于使用场地的自下而上的工程机械排放清单编制方法 | |
| CN115166210A (zh) | 一种基于场地概念模型的工业场地土壤采样布点方法 | |
| Qian et al. | A prospective ecological risk assessment method based on exposure and ecological scenarios (ERA-EES) to determine soil ecological risks around metal mining areas | |
| Okmyanskaya et al. | Ecological and Radiation Assessment of the State of Drilling Waste and their Disposal Products on the Example of the Khanty-Mansiysk Autonomous Okrug–Yugra | |
| CN116956842A (zh) | 工业场地土壤污染物排放清单的构建方法及*** | |
| Fordyce et al. | Assessing threats to shallow groundwater quality from soil pollutants in Glasgow, UK: development of a new screening tool | |
| Yang et al. | Soil assessment | |
| Khaustov et al. | IT for the remediation of the geological environment poluted with the petroleum products: Experience of the Kazakh-Belarus Russian joint project | |
| Wang et al. | GIS applications in seismic loss estimation model for Portland, Oregon water and sewer systems | |
| Glinsky et al. | Organization of the SM–NI System | |
| Titov et al. | Geoenvironmetal assessment of different types of cryolithic soils in Western Yakutia under the conditions of diamond-mining operations | |
| Yang et al. | SPATIAL IDENTIFICATION OF HEAVY METAL POLLUTION SOURCES IN THE ZHUXIANZHUANG MINING REGION, CHINA. | |
| None | Streamlined Approach for Environmental Restoration (SAFER) Plan for Corrective Action Unit 572: Test Cell C Ancillary Building and Structures-Nevada National Security Site, Nevada | |
| Tseytlin | THE ROLE OF FORECASTING IN THE MANAGEMENT OF ENVIRONMENTAL ENGINEERING ISSUES AT ENTERPRISES OF THE MINERALS SECTOR IN THE CONTEXT OF NEW CHALLENGES AND THREATS |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |