CN111259543B - 一种石油类污染场地中sve工艺参数的设计方法 - Google Patents
一种石油类污染场地中sve工艺参数的设计方法 Download PDFInfo
- Publication number
- CN111259543B CN111259543B CN202010042282.7A CN202010042282A CN111259543B CN 111259543 B CN111259543 B CN 111259543B CN 202010042282 A CN202010042282 A CN 202010042282A CN 111259543 B CN111259543 B CN 111259543B
- Authority
- CN
- China
- Prior art keywords
- sve
- parameters
- site
- model
- petroleum
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 230000008569 process Effects 0.000 title claims abstract description 44
- 239000003208 petroleum Substances 0.000 title claims abstract description 43
- 238000013461 design Methods 0.000 title claims abstract description 32
- 230000008439 repair process Effects 0.000 claims abstract description 59
- 238000012417 linear regression Methods 0.000 claims abstract description 42
- 238000004088 simulation Methods 0.000 claims abstract description 28
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 17
- 231100000719 pollutant Toxicity 0.000 claims abstract description 17
- 238000012216 screening Methods 0.000 claims abstract description 13
- 238000013508 migration Methods 0.000 claims abstract description 12
- 230000005012 migration Effects 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims description 31
- 238000000605 extraction Methods 0.000 claims description 25
- 239000002689 soil Substances 0.000 claims description 24
- 230000000694 effects Effects 0.000 claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000035699 permeability Effects 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000011065 in-situ storage Methods 0.000 claims description 5
- 239000002957 persistent organic pollutant Substances 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 239000000356 contaminant Substances 0.000 claims description 4
- 238000001764 infiltration Methods 0.000 claims description 4
- 230000008595 infiltration Effects 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 4
- 239000012855 volatile organic compound Substances 0.000 claims description 4
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 3
- 230000006399 behavior Effects 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 239000001301 oxygen Substances 0.000 claims description 3
- 229910052760 oxygen Inorganic materials 0.000 claims description 3
- 108010014173 Factor X Proteins 0.000 claims description 2
- 239000003570 air Substances 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 239000002680 soil gas Substances 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000012795 verification Methods 0.000 abstract description 2
- 238000011109 contamination Methods 0.000 description 7
- 239000003673 groundwater Substances 0.000 description 7
- 239000012071 phase Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- UHOVQNZJYSORNB-UHFFFAOYSA-N Benzene Chemical compound C1=CC=CC=C1 UHOVQNZJYSORNB-UHFFFAOYSA-N 0.000 description 3
- YXFVVABEGXRONW-UHFFFAOYSA-N Toluene Chemical compound CC1=CC=CC=C1 YXFVVABEGXRONW-UHFFFAOYSA-N 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- IVSZLXZYQVIEFR-UHFFFAOYSA-N 1,3-Dimethylbenzene Natural products CC1=CC=CC(C)=C1 IVSZLXZYQVIEFR-UHFFFAOYSA-N 0.000 description 2
- YNQLUTRBYVCPMQ-UHFFFAOYSA-N Ethylbenzene Chemical compound CCC1=CC=CC=C1 YNQLUTRBYVCPMQ-UHFFFAOYSA-N 0.000 description 2
- 238000003912 environmental pollution Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- URLKBWYHVLBVBO-UHFFFAOYSA-N p-dimethylbenzene Natural products CC1=CC=C(C)C=C1 URLKBWYHVLBVBO-UHFFFAOYSA-N 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- WYLIRYQDDKDHLT-UHFFFAOYSA-N CC1=CC=CC=C1C.CC1=CC=CC=C1C Chemical group CC1=CC=CC=C1C.CC1=CC=CC=C1C WYLIRYQDDKDHLT-UHFFFAOYSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 239000008346 aqueous phase Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000013064 chemical raw material Substances 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- -1 monocyclic aromatic hydrocarbon Chemical class 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 239000000575 pesticide Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- VYMDGNCVAMGZFE-UHFFFAOYSA-N phenylbutazonum Chemical compound O=C1C(CCCC)C(=O)N(C=2C=CC=CC=2)N1C1=CC=CC=C1 VYMDGNCVAMGZFE-UHFFFAOYSA-N 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000005067 remediation Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000012209 synthetic fiber Substances 0.000 description 1
- 229920002994 synthetic fiber Polymers 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 150000003738 xylenes Chemical class 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种石油类污染场地中SVE工艺参数的设计方法，步骤如下：步骤一、明确场地和石油类污染状况；步骤二、结合场地参数、污染参数和SVE工艺参数，建立TOUGH软件关于SVE的修复模型并得出修复率y；步骤三、利用灰色关联度方法，筛选p个主控因子；步骤五、将根据多元线性回归方程的模拟结果进行拟合验证，并确定其模拟精度是否复合设计要求。步骤六、筛选最佳的SVE工艺参数组合并应用。本发明通过典型场地和多因素变量的设计并结合TOUGH模拟，可明晰污染物在场地中不同场景下的运移规律，利于对SVE工艺的参数设计选取；基于灰色关联度方法，找出主控因子，利于应对污染治理时的重点设计；通过建立多元线性回归方程，对应类似的场地更具有适用性。
Description
技术领域
本发明属于石油类污染场地治理技术领域，特别涉及一种石油类污染场地中SVE工艺参数的设计方法。
背景技术
随着工业化的进程，石油的使用也随之增加。在长期使用中，石油的跑冒滴漏时有发生，而石油类污染也日益成为环境污染中重要的组成部分。因石油类有机污染物具有挥发性和流动性，其在渗透性较大的土层中，会发生迁移有的甚至渗透到地下水中，在水土相互作用下产生更大范围的污染，加重污染程度。
对石油污染物进行治理，需明晰石油类污染物在土壤中迁移转化过程。TOUGH是非饱和地下水流及热流传输(Transport of Unsaturated Groundwater and Heat)的英文缩写,是一个模拟一维、二维和三维孔隙或裂隙介质中多相流(multi2phase)、多组分(multi2component)及非等温(non2isothermal)的水流及热量运移的数值模拟程序。应用TOUGH软件对石油类污染场地进行模拟，可准确的掌握石油类污染物的迁移转化过程。
在对石油类污染场地的治理技术中土壤气相抽提(soil vapor extraction，SVE)是对土壤挥发性有机污染进行原位修复的一种方法，用来处理包气带中地层介质的污染问题。目前有关SVE的研究大多停留在室内实验和现场试验上，其设计与运行大多根据经验公式或者在有限的场地实践开展，针对SVE数值模拟方面的研究存在不足。对于SVE修复中应用不同场地、不同污染物的影响因素以及影响比重没有针对性设计、对于SVE技术参数与数学模型的联用没有相对设计，以及缺少对于不同场地的验证适用性的设计。
发明内容
本发明提供了一种石油类污染场地中SVE工艺参数的设计方法，用以解决石油类污染场地中不同影响因素对于SVE修复率的影响、不同影响因素的筛选和数学模型的组合应用以及对于SVE工艺参数方程模型的验证等技术问题。
为实现上述目的，本发明采用如下技术方案：
一种石油类污染场地中SVE工艺参数的设计方法，具体步骤如下：
步骤一、根据场地实际踏勘、原位测试及土工试验结果，结合场地的地质资料，明确场地的地质类型、土壤类型和分布状况以及地下水分布等状况；并确定石油类污染的类型和位置；
步骤二、结合场地参数、污染参数和SVE工艺参数，利用TOUGH软件建立石油类污染场地的关于SVE(土壤气相抽提)修复模型，并得出不同影响因子条件下SVE修复率y；
其中，SVE修复率反映了SVE对场地中石油类污染的去除效果，SVE修复率yk计算公式如下所示：
式中：mk为SVE修复前模型中待去除污染物的总质量，kg；m′k为SVE修复后模型中污染物总质量，kg；其中k＝1,2,3，...，w，w为场地数量。
步骤三、利用灰色关联度分析方法，进行不同场地类型下SVE在不同影响因子修复率的关联度对比排名，通过排名前后筛选p个主控因子；
其中，不同场地的编号序列计为k(k＝1,2,3，......，n)，设Xi为SVE修复率影响因子，xi(k)为因素xi在场地k的观测数据；则{xi(k)│k＝1,2,3，......，n}是SVE效果行为序列，其中i＝1,2,3，...，m，m为影响因子的数量；设y(k)为场地k的SVE修复率；
计算关联度ri如下所示：
式中，ζi(k)为关联度系数。
步骤四、利用灰色关联度分析筛选出与SVE修复率相关的p个主控因子，建立因变量修复率y与自变量主控因子Xi(i＝1，2，……，p)之间的多元线性回归方程；
y与p个主控因子Xi之间的多元线性回归方程关系式如下所示：
y＝b0+b1X1+b2X2+…+bpXp
式中，b0，b1，b2，……，bp，b0是一个常量和其他为p个主控因子的待定系数；由最小二乘法求得多元线性回归方程的b0和其他待定系数。
步骤五、基于多元线性回归方程的模拟结果，采用拟合优度的方式进行检验并判断其模拟的准确性；再通过显著性检验判定多元线性回归方程的模型显著性和参数的显著性；最后再对多元线性回归方程模型的误差进行精度对比，确定其精度是否复合设计要求。
步骤六、利用建立的多元线性回归方程，将掌握的新场地特征参数代入多元线性回归方程，通过设定SVE修复率目标，筛选最佳的SVE工艺参数组合，为污染场地SVE修复技术参数设计提供技术参考。
进一步的，对于步骤一中场地的选择为具有代表性的典型地块或对同一地块根据地质进行区域划分，并对场地垂向的土质分层进行概化，概化土层包含石油类有机污染物迁移转化土层和SVE施用土层。
进一步的，对于步骤二中，通过多个典型污染场地进行TOUGH软件概化模拟，得出场地参数、污染参数和SVE工艺参数中不同影响因子对于SVE修复率的影响，并得出相应SVE修复率，且对于同一影响因子的量级变化幅度同样进行对比模拟。
进一步的，对于步骤二中选取SVE修复效率的影响因子包括场地参数中的入渗量、非饱和带厚度、孔隙度、渗透率、含氧量、温度和pH值；污染参数中包括污染物类型、深度、宽度和面积；SVE工艺参数中包括抽提井内流量、影响半径、抽提井深度和抽提井数量。
进一步的，对于TOUGH软件模拟过程根据污染物的不同选取不同的模块，模块包含T2VOC模块和TMVOC模块；T2VOC模块为三种组分三相流，包括水、空气及VOCs数值模拟，TMVOC模块为水、土壤气体和多组分混合挥发性有机化合物在多层面、非均质、多孔介质中的三相非等温流数值模拟；TOUGH软件模拟过程通过PetraSim软件进行可视化操作。
进一步的，在TOUGH软件模拟过程中对于不同的污染物，通过设置相同的泄露速率、泄漏点，泄露持续时间等初始参数，以及同一场地参数和SVE工艺参数，并通过已有实验或现场数据校正模型，得出具有可对比性的SVE修复率。
进一步的，对于步骤三中关联度ri的计算如下：
第1步，无量纲化如下所示：
第2步，求差序列如下所示：
Δi(k)＝|y(k)-x′i(k)|
i＝1,2,…,m；k＝1,2,…,n
第3步，求两级最大差与最小差如下所示：
第4步，求关联度系数如下所示：
第5步，计算关联度如下所示：
进一步的，对于步骤五中，关于多元线性回归方程的模型拟合程度通过拟合优度检验；
拟合优度检验公式如下所示：
式中：TSS是总离差平方和，ESS是回归平方和，RSS是剩余平方和。R2越接近于1，多元线性回归方程的模型拟合度越好。
进一步的，对于步骤五中，多元线性回归方程的显著性检验F如下所示：
式中，n是样本数、p是选取的变量，若F≥Fα(p，n-p-1)，回归模型显著成立；若F<Fα(p，n-p-1)，回归模型无显著意义，即回归模型不显著；
参数的显著性检验t如下所示：
式中，bi代表为回归系数，S(bi)代表为回归系数bi的标准差，若表明xi对y有显著性作用；若/>表明xi对y的作用不显著；此外，参数t值的检验可简化为概率检验，t值的概率若小于0.05则此自变量是显著的。
进一步的，对于步骤五中，多元线性回归方程的误差分析具体步骤如下：
①求原始数据均值如下所示：
②求原始数据方差S1如下所示：
③求残差ε均值如下所示：
ε(k)＝Y(k)-Y＇(k)
④求残差方差S2如下所示：
⑤计算方差比C与小概率误差P如下所示：
当后验差比值C小于0.5时，认为模型精度合格，C越小模型精度越高。小误差概率P大于0.8时，认为模型精度合格，P越大模型精度越高。
本发明的有益效果体现在：
1)本发明通过多场地多因素的设计，并结合TOUGH软件模拟不同场景下的污染状况，可明晰石油类污染物在场地中不同场景下的运移规律，确定污染范围，利于后期对于SVE工艺的参数设计；
2)本发明基于灰色关联度方法，对典型场地的不同影响因素进行排列，找出主控因子，利于应对同种污染在同种地块上的SVE工艺设计时的重点设计，提高施用效率；
3)本发明通过建立多元线性回归方程，将SVE修复率相关的主控因子进一步排序和进行修正，使其对应类似的典型场地更具有适用性；
4)本发明通过拟合优度、显著性检验和误差检验等方法，进一步确保多元线性回归方程于实际应用中的准确性和可靠性；
此外，本发明可根据不同的场地和对应污染物进行单独的模拟和设计，也可根据以往和现有设计预先确定相关影响参数，而后进行关联度的筛选，可适用性强；本发明的其它特征和优点将在随后的说明书中阐述，并且部分地从说明书中变得显而易见，或者通过实施本发明而了解；本发明的主要目的和其它优点可通过在说明书中所特别指出的方案来实现和获得。
附图说明
图1是石油类污染物中BTEX污染场地SVE工艺参数设计方法流程示意图；
图2是基于TOUGH软件模拟BTEX污染场地PetraSim软件可视化示意图；
图3是S1污染场地BTEX污染晕剖面示意图；
图4是S1污染场地BTEX污染应有抽提技术剖面示意图。
具体实施方式
石油类有机污染物中BTEX污染是具有代表性的一类。BTEX是石油中常见的苯(benzene)，甲苯(toluene)、乙基苯(ethylbenzene)、三种二甲基苯的异构体(o-xylene邻二甲苯，m-xylene间二甲苯，p-xylene对二甲苯)的合称，属于单环芳烃类质。BTEX主要存在于原油和石油产品中，作为化工原料，广泛应用于农药、塑料和合成纤维等制造业。BTEX在生产、储存和运输过程中，容易释放到环境内，造成环境污染，并对生态系统和人体健康构成危害。在BTEX污染处理中，因BTEX具有挥发性，在渗透性大的土壤中可通过SVE去除。
以BTEX为实施例，如图1所示，BTEX污染场地SVE工艺参数设计方法流程示意图，具体步骤如下：
步骤一、根据场地实际踏勘、原位测试及土工试验结果，结合场地的地质资料，明确场地的地质类型、土壤类型和分布状况以及地下水分布等状况；并确定石油类污染的类型和位置；
其中，对于步骤一中场地的选择为具有代表性的典型地块或对同一地块根据地质进行区域划分，并对场地垂向的土质分层进行概化，概化土层包含石油类有机污染物迁移转化土层和适应SVE施用土层。
根据BTEX的密度、蒸气压和可溶性等物理和化学参数的不同，本实施例中以收集的来自全国的15个典型地区(用S1，S2，S3，......S15表示)的数据进行不同场地下BTEX污染的模拟设计，15个典型地区的场地参数如表1所示。
表1场地参数
步骤二、步骤二、结合场地参数、污染参数和SVE工艺参数，利用TOUGH软件建立石油类污染场地的关于SVE(土壤气相抽提)修复模型，并得出不同影响因子条件下SVE修复率y；
其中，SVE修复率反映了SVE对场地中石油类污染的去除效果，SVE修复率yk计算公式如式(1)所示：
式中：mk为SVE修复前模型中待去除BTEX的总质量，kg；m′k为SVE修复后模型中BTEX总质量，kg；其中k＝1,2,3，...，w，w为场地数量。
选取SVE修复效率的影响因子包括场地参数中的入渗量、非饱和带厚度、孔隙度、渗透率、含氧量、温度和pH值；污染参数中包括污染物类型、深度、宽度和面积；SVE工艺参数中包括抽提井内流量、影响半径、抽提井深度和抽提井数量。
选取TOUGH软件进行过程模拟，根据BTEX选取TOUGH软件中TMVOC子模块进行建立模型，并通过PetraSim软件进行可视化操作。以S1场地为例，通过建立概念模型、基本的参数设置、边界条件和初始条件的设置、模型的运行和调试以及最后建立完整的模型，其中，模型建立界面如图2所示。对15个场地均进行BTEX泄露及迁移的模拟，设定BTEX以相同的速率及泄漏点，泄露持续时间1年。BTEX泄露及迁移后场地污染情况如表2所示，其中NAPL(Non-aqueous Phase Liquid)为非水相流体，为BTEX在场地中存在的一种相态。
表2 BTEX泄露迁移后场地污染模拟值
以场地S1为例，模型抽提前BTEX污染分布如图3所示，抽提过程为图4所示；图中h代表深度，b代表宽度，w代表BTEX的质量分数，三角形处为地下水水位面；在图4中深度为0m处间隔分布的竖向线段为抽提井。由此，对应各场地模拟BTEX污染情况，设计个场地SVE抽提参数如表3所示。
表3场地SVE参数
根据SVE模拟前后场地中BTEX的统计，计算SVE修复率如表4所示。
表4场地SVE去除率模拟值
步骤三、利用灰色关联度分析方法，进行不同场地类型下SVE在不同影响因子修复率的关联度对比排名，通过排名前后筛选p个主控因子；
其中，不同场地的编号序列计为k(k＝1,2,3，......，n)，设xi为SVE修复率影响因子，xi(k)为因素xi在场地k的观测数据；则{xi(k)│k＝1,2,3，......，n}是SVE效果行为序列，其中i＝1,2,3，...，m，m为影响因子的数量；设y(k)为场地k的SVE修复率；
第1步，无量纲化如式(2)所示：
第2步，求差序列如式(3)所示：
Δi(k)＝|y(k)-x′i(k)| (3)第3步，求两级最大差与最小差如式(4)所示：
第4步，求关联度系数如式(5)所示：
第5步，计算关联度如式(6)所示：
选择SVE修复效率的影响因子共11个，包括场地参中的入渗量(x1)、非饱和带厚度(x2)、孔隙度(x3)、渗透率(x4)；污染参数中的污染深度(x5)、宽度(x6)、面积(x7)以及抽提参数中的井内流量(x8)、影响半径(x9)、抽提井深度(x10)、抽提井数量(x11)，各参数数值如表5所示。
表5灰色关联度计算结果
根据上表格给出的数据信息，利用灰色关联度法分析与SVE修复率较高的主控因子。计算各个影响因子xi与SVE去除率模拟值y之间的相关度ri关联度结果见表。
比较得知：x4>x7>x2>x8>x10>x11>x5>x1>x6>x3>x9
将关联度计算结果和场地数据获取的难易程度两个方面作为回归模型参数选取的依据。采用场地参数中非饱和带厚度x2、非饱和带渗透率x4，污染参数中污染面积x7，抽提参数中井内流量x8，抽提井深度x10，抽提井数量x11作为建立多元线性回归模型的6个因子。
步骤四、利用灰色关联度分析筛选出与SVE修复率相关的p个主控因子，建立因变量SVE修复率y与自变量主控因子xi(i＝1，2，……，p)之间的多元线性回归方程；
多元线性回归方程的建立如式(7)所示，将6个主控因子代入到下式中得：
Y＇＝b0+b2x2+b4x4+b7x7+b8x8+b10x10+b11x11 (7)
利用SPSS软件计算得出方程关系式，如式(8)所示：
Y＇＝0.596+0.016x2-3.52×109×x4+0.001x7+3.71x8-0.096x10+0.034x11 (8)
步骤五、基于多元线性回归方程的模拟结果，采用拟合优度的方式进行检验并判断模拟的拟合准确性；再通过显著性检验判定多元线性回归方程的模型显著性和参数的显著性；最后再对多元线性回归方程的模型的误差进行精度对比，确定其精度是否复合设计要求。
拟合优度检验公式如式(9)所示：
式中：TSS是总离差平方和，ESS是回归平方和，RSS是剩余平方和。R2为0.819接近于1，多元线性回归方程的模型拟合度好。
多元线性回归方程的显著性检验F如式(10)所示：
式中，n是样本数、p是选取变量，若F≥Fα(p，n-p-1)，回归模型显著成立；若F<Fα(p，n-p-1)，回归模型无显著意义，即回归模型不显著。
在SPSS软件建立好回归方程后自动提供的模型检验参数中，拟合优度检验参数F为6.047，变量(p)的数量为6；样本数(n)为15。从标准F统计表查得，当显着性α＝0.05时，F(6，8)＝3.581。可知6.047>3.581，方程的显著性很高，具有统计学意义。
2)参数的显著性检验t如式(11)所示：
式中，bi代表为回归系数，S(bi)代表为回归系数bi的标准差，若表明xi对y有显著性作用；若/>表明xi对y的作用不显著；此外，参数t值的检验可简化为概率检验，t值的概率若小于0.05则此自变量是显著的，通过SPSS软件计算，得出表6。
表6 T检验结果
表中常量和自变量x7、x10的显著性sig值小于0.05，则这两个变量的系数和常量十分显著，而其他4个参数并不显著。由于已经先通过灰色关联度计算出其他4个参数与SVE修复效率密切相关，因此，这些参数被保留。
经过模型检验后，最终得到的回归方程关系式如式(12)所示：
Y＇＝0.596+0.016x2-3.52×109×x4+0.001x7+3.71x8-0.096x10+0.034x11 (12)
3)多元线性回归方程的误差分析具体步骤如下：
①求原始数据均值如式(13)所示：
②求原始数据方差S1如式(14)所示：
③求残差ε均值如式(15)和(16)所示：
ε(k)＝Y(k)-Y＇(k) (15)
④求残差方差S2如下所示：
⑤计算方差比C与小概率误差P，如式(18)和(19)所示：
当后验差比值C小于0.5时，认为模型精度合格，C越小模型精度越高。小误差概率P大于0.8时，认为模型精度合格，P越大模型精度越高。通过计算本次建立的方程的后验差比值C为0.43，该值小于0.5；方程的小误差概率P为0.93，该值大于0.8。因此，方程是合格的。
步骤六、利用建立的多元线性回归方程，将掌握的新场地特征参数代入多元线性回归方程，通过设定SVE修复率目标，筛选最佳的SVE工艺参数组合，为污染场地SVE修复技术参数设计提供技术参考。
以中国南方某工厂为例，对SVE修复率的多元线性回归模型进行验证。厂区位于中国华东地区，该厂区所在区域地表均为松散堆积物覆盖，根据场地实际踏勘、原位测试及土工试验结果，可知研究区为第四季地层所覆盖，第四纪地层可划分为下更新统(Q1)、中更新统(Q2)、上更新统(Q3)及全新统(Q4)，主要岩性为河湖相沉积的粘土、亚粘土夹粉细砂层。厂区内含水层主要由第四季上更新统粉土、粉砂、细砂组成，分布较连续、稳定。该区地下水主要为潜水，目前厂区地下水埋深3～4.5m。
参数赋值及SVE修复率的计算，模型相关参数设置见表7和表8。
表7场地参数
表8抽提参数
经过TOUGH软件的模拟后，计算SVE修复率的结果为：67％。结合厂区场地的实际情况给出6个参数的数值分别为：x2为3m，x4为3.07E-12m2，x7为220m2，x8为0.00814m3/s，x10为3.5m，x11为5个，将其建立的多元线性回归方程计算得y的值为72％。评价结果对比，TOUGH2的SVE修复率模拟值为67％，建立的多元线性回归方程计算得y的值为72％。结论差距在5％以内，结论相符。
SVE修复石油类污染场地是一个复杂的动力学过程，由于针对SVE修复率的影响因素较多，不同的影响因素对SVE修复率的贡献不同。尽管SVE技术在许多实地研究中应用广泛，但是目前的理论研究特别是对SVE过程中流体的运移机制、污染物传质机理和现场尺寸放大效应和综合数学模拟仍不充分。SVE法修复污染场地工艺设计的精确度关系到修复效果好坏及修复成本的高低，因此科学、快速、准确地设定SVE工艺参数，对于典型污染场地SVE的准确设计、缩短修复时间和节约修复成本等工作具有重要指导性意义。
以上所述仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内所想到的变化或替换，都应涵盖在本发明的保护范围之内。
Claims (7)
1.一种石油类污染场地中SVE工艺参数的设计方法，其特征在于，具体步骤如下：
步骤一、根据场地实际踏勘、原位测试及土工试验结果，结合场地的地质资料，明确场地的地质类型、土壤类型和分布状况以及地下水分布状况；并确定石油类污染的类型和位置；
对于步骤一中场地的选择为具有代表性的典型地块或对同一地块根据地质进行区域划分，并对场地垂向的土质分层进行概化，概化土层包含石油类有机污染物迁移转化土层和SVE施用土层；
步骤二、结合场地参数、污染参数和SVE工艺参数，利用TOUGH软件建立石油类污染场地的关于SVE修复模型，并得出不同影响因子条件下SVE修复率y；
其中，SVE修复率反映了SVE对场地中石油类污染的去除效果，SVE修复率yk计算公式如下所示：
式中：mk为SVE修复前模型中待去除污染物的总质量，kg；m′k为SVE修复后模型中污染物总质量，kg；其中k＝1,2,3，...，w，w为场地数量；
对于步骤二中，通过多个典型污染场地进行TOUGH软件概化模拟，得出场地参数、污染参数和SVE工艺参数中不同影响因子对于SVE修复率的影响，并得出相应SVE修复率，且对于同一影响因子的量级变化幅度同样进行对比模拟；
对于步骤二中选取SVE修复效率的影响因子包括场地参数中的入渗量、非饱和带厚度、孔隙度、渗透率、含氧量、温度和pH值；污染参数中包括污染物类型、深度、宽度和面积；SVE工艺参数中包括抽提井内流量、影响半径、抽提井深度和抽提井数量；
步骤三、利用灰色关联度分析方法，进行不同场地类型下SVE在不同影响因子修复率的关联度对比排名，通过排名前后筛选p个主控因子；
其中，不同场地的编号序列计为k，设Xi为SVE修复率影响因子，xi(k)为因素xi在场地k的观测数据；则{xi(k)│k＝1,2,3，......，n}是SVE效果行为序列，其中i＝1,2,3，...，m，m为影响因子的数量；设y(k)为场地k的SVE修复率；
计算关联度ri如下所示：
式中，ζi(k)为关联度系数；
步骤四、利用灰色关联度分析筛选出与SVE修复率相关的p个主控因子，建立因变量修复率y与自变量主控因子Xi之间的多元线性回归方程；
y与p个主控因子Xi之间的多元线性回归方程关系式如下所示：
y＝b0+b1X1+b2X2+…+bpXp
式中，b0，b1，b2，……，bp，b0是一个常量和其他为p个主控因子的待定系数；由最小二乘法求得多元线性回归方程的b0和其他待定系数；
步骤五、基于多元线性回归方程的模拟结果，采用拟合优度的方式进行检验并判断其模拟的准确性；再通过显著性检验判定多元线性回归方程的模型显著性和参数的显著性；最后再对多元线性回归方程模型的误差进行精度对比，确定其精度是否复合设计要求；
步骤六、利用建立的多元线性回归方程，将掌握的新场地特征参数代入多元线性回归方程，通过设定SVE修复率目标，筛选最佳的SVE工艺参数组合，为污染场地SVE修复技术参数设计提供技术参考。
2.如权利要求1所述的一种石油类污染场地中SVE工艺参数的设计方法，其特征在于，对于TOUGH软件模拟过程根据污染物的不同选取不同的模块，模块包含T2VOC模块和TMVOC模块；T2VOC模块为三种组分三相流，包括水、空气及VOCs数值模拟，TMVOC模块为水、土壤气体和多组分混合挥发性有机化合物在多层面、非均质、多孔介质中的三相非等温流数值模拟；TOUGH软件模拟过程通过PetraSim软件进行可视化操作。
3.如权利要求2所述的一种石油类污染场地中SVE工艺参数的设计方法，其特征在于，在TOUGH软件模拟过程中对于不同的污染物，通过设置相同的泄露速率、泄漏点，泄露持续时间初始参数，以及同一场地参数和SVE工艺参数，并通过已有实验或现场数据校正模型，得出具有可对比性的SVE修复率。
4.如权利要求1所述的一种石油类污染场地中SVE工艺参数的设计方法，其特征在于，对于步骤三中关联度ri的计算如下：
第1步，无量纲化如下所示：
第2步，求差序列如下所示：
Δi(k)＝|y(k)-x′i(k)|
i＝1,2,…,m；k＝1,2,…,n
第3步，求两级最大差与最小差如下所示：
第4步，求关联度系数如下所示：
第5步，计算关联度如下所示：
5.如权利要求1所述的一种石油类污染场地中SVE工艺参数的设计方法，其特征在于，对于步骤五中，关于多元线性回归方程的模型拟合程度通过拟合优度检验；
拟合优度检验公式如下所示：
式中：TSS是总离差平方和，ESS是回归平方和，RSS是剩余平方和；R2越接近于1，多元线性回归方程的模型拟合度越好。
6.如权利要求1所述的一种石油类污染场地中SVE工艺参数的设计方法，其特征在于，对于步骤五中，多元线性回归方程的显著性检验F如下所示：
式中，n是样本数、p是选取的变量，若F≥Fα(p，n-p-1)，回归模型显著成立；若F<Fα(p，n-p-1)，回归模型无显著意义，即回归模型不显著；
参数的显著性检验t如下所示：
式中，bi代表为回归系数，S(bi)代表为回归系数bi的标准差，若表明xi对y有显著性作用；若/>表明xi对y的作用不显著；此外，参数t值的检验简化为概率检验，t值的概率若小于0.05则此自变量是显著的。
7.如权利要求1所述的一种石油类污染场地中SVE工艺参数的设计方法，其特征在于，对于步骤五中，多元线性回归方程的误差分析具体步骤如下：
①求原始数据均值如下所示：
②求原始数据方差S1如下所示：
③求残差ε均值如下所示：
ε(k)＝Y(k)-Y＇(k)
④求残差方差S2如下所示：
⑤计算方差比C与小概率误差P如下所示：
当后验差比值C小于0.5时，认为模型精度合格，C越小模型精度越高；小误差概率P大于0.8时，认为模型精度合格，P越大模型精度越高。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010042282.7A CN111259543B (zh) | 2020-01-15 | 2020-01-15 | 一种石油类污染场地中sve工艺参数的设计方法 |
| PCT/CN2020/077574 WO2021142895A1 (zh) | 2020-01-15 | 2020-03-03 | 一种石油类污染场地中sve工艺参数的设计方法 |
| US17/060,937 US20210216681A1 (en) | 2020-01-15 | 2020-10-01 | Method for designing sve process parameters in petroleum-type polluted field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010042282.7A CN111259543B (zh) | 2020-01-15 | 2020-01-15 | 一种石油类污染场地中sve工艺参数的设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111259543A CN111259543A (zh) | 2020-06-09 |
| CN111259543B true CN111259543B (zh) | 2024-01-23 |
Family
ID=70948741
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010042282.7A Active CN111259543B (zh) | 2020-01-15 | 2020-01-15 | 一种石油类污染场地中sve工艺参数的设计方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN111259543B (zh) |
| WO (1) | WO2021142895A1 (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117789847A (zh) * | 2024-02-26 | 2024-03-29 | 中国石油大学(华东) | 一种基于化学需氧量的压裂返排液影响因素主次分析方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE4315550A1 (de) * | 1993-05-10 | 1994-11-17 | Mabo Fahrzeug Anlagenbau | Verfahren und Vorrichtung zum Dekontaminieren von mit Kohlenwasserstoffen verunreinigten Bodenmassen |
| CN102910712A (zh) * | 2012-11-01 | 2013-02-06 | 广西大学 | 电渗析处理酵母废水工艺条件的选择方法 |
| CN104759459A (zh) * | 2015-04-07 | 2015-07-08 | 山东师范大学 | 高浓度石油污染土壤异位-原位联合生物修复方法 |
| CN105260514A (zh) * | 2015-09-24 | 2016-01-20 | 中国环境科学研究院 | 地下水污染源强定量评价方法 |
| CN105787297A (zh) * | 2016-03-12 | 2016-07-20 | 云南圣清环境监测科技有限公司 | 一种微生物修复系统活性评价的方法 |
| JP2017154073A (ja) * | 2016-03-01 | 2017-09-07 | 株式会社竹中工務店 | 汚染対処方法、及び汚染対処プログラム |
-
2020
- 2020-01-15 CN CN202010042282.7A patent/CN111259543B/zh active Active
- 2020-03-03 WO PCT/CN2020/077574 patent/WO2021142895A1/zh active Application Filing
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE4315550A1 (de) * | 1993-05-10 | 1994-11-17 | Mabo Fahrzeug Anlagenbau | Verfahren und Vorrichtung zum Dekontaminieren von mit Kohlenwasserstoffen verunreinigten Bodenmassen |
| CN102910712A (zh) * | 2012-11-01 | 2013-02-06 | 广西大学 | 电渗析处理酵母废水工艺条件的选择方法 |
| CN104759459A (zh) * | 2015-04-07 | 2015-07-08 | 山东师范大学 | 高浓度石油污染土壤异位-原位联合生物修复方法 |
| CN105260514A (zh) * | 2015-09-24 | 2016-01-20 | 中国环境科学研究院 | 地下水污染源强定量评价方法 |
| JP2017154073A (ja) * | 2016-03-01 | 2017-09-07 | 株式会社竹中工務店 | 汚染対処方法、及び汚染対処プログラム |
| CN105787297A (zh) * | 2016-03-12 | 2016-07-20 | 云南圣清环境监测科技有限公司 | 一种微生物修复系统活性评价的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 李娟 等.灰色关联度分析在地下水污染源强定量评价中的应用.《安全与环境学报》.2015,第15卷(第4期),第342-347页第1-3节. * |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2021142895A1 (zh) | 2021-07-22 |
| CN111259543A (zh) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210216681A1 (en) | Method for designing sve process parameters in petroleum-type polluted field | |
| Johnson et al. | A practical approach to the design, operation, and monitoring of in situ soil‐venting systems | |
| AU2020102747A4 (en) | A decision-making method for in-situ remediation of petroleum-contaminated groundwater | |
| CN110826236A (zh) | 一种土壤修复污染范围的判定方法 | |
| CN111259543B (zh) | 一种石油类污染场地中sve工艺参数的设计方法 | |
| Chang et al. | Extended Kalman filtering to improve the accuracy of a subsurface contaminant transport model | |
| Schäfer et al. | Simulating transport and removal of xylene during remediation of a sandy aquifer | |
| Reilly et al. | Effects of physical and chemical heterogeneity on water‐quality samples obtained from wells | |
| Yang et al. | Investigation of soil and groundwater environment in urban area during post-industrial era: A case study of brownfield in Zhenjiang, Jiangsu Province, China | |
| Sepehr et al. | In situ soil remediation using vapor extraction wells, development and testing of a three‐dimensional finite‐difference model | |
| Guo et al. | The impact of well‐field configuration on contaminant mass removal and plume persistence for homogeneous versus layered systems | |
| Rosenbloom et al. | Application of VLEACH to vadose zone transport of VOCs at an Arizona superfund site | |
| CN114798706A (zh) | 一种用于污染土壤地下水原位修复药剂的高精度注入方法 | |
| CN109307690B (zh) | 一种调查石油烃土壤污染的方法、装置及电子设备 | |
| Petrula et al. | Evaluation of groundwater contamination in former military airport area | |
| Agah et al. | 2D simulation of dynamic transportation of volatile hydrocarbons in vadose zone of tehran oil refinery and industrial area of Ray, Tehran, Iran | |
| Alberti et al. | Fingerprinting and groundwater model application to evaluate hydraulic barrier efficiency (Italy) | |
| Chao | Scale dependence of transport parameters estimated from forced-gradient tracer tests in heterogeneous formations | |
| Mescher | Evaluating the Applicability of MODFLOW2005 Plus MT3D-USGS: Borden Case Study | |
| CN117314186A (zh) | 一种预测有机污染场地自然衰减修复效果的方法 | |
| Bear | Transport modeling | |
| Keshari | Recent trends in flow through porous media | |
| Lin et al. | Tracking chlorinated contaminants in the subsurface using analytical, numerical and geophysical methods | |
| CN118153279A (zh) | 一种焦化企业地块土壤污染范围划定的方法 | |
| Chakraborty et al. | Identification of Clandestine Groundwater Pollution Sources in Aquifer Polluted with BTEX |
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 |