CN111259558B - 一种非饱和土壤中胶体运动分布三维可视化的方法 - Google Patents
一种非饱和土壤中胶体运动分布三维可视化的方法 Download PDFInfo
- Publication number
- CN111259558B CN111259558B CN202010071548.0A CN202010071548A CN111259558B CN 111259558 B CN111259558 B CN 111259558B CN 202010071548 A CN202010071548 A CN 202010071548A CN 111259558 B CN111259558 B CN 111259558B
- Authority
- CN
- China
- Prior art keywords
- colloid
- soil
- distribution
- dimensional visualization
- potential
- 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.)
- Expired - Fee Related
Links
- 239000000084 colloidal system Substances 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000002689 soil Substances 0.000 title claims abstract description 38
- 238000009826 distribution Methods 0.000 title claims abstract description 28
- 238000012800 visualization Methods 0.000 title claims abstract description 15
- SZVJSHCCFOBDDC-UHFFFAOYSA-N iron(II,III) oxide Inorganic materials O=[Fe]O[Fe]O[Fe]=O SZVJSHCCFOBDDC-UHFFFAOYSA-N 0.000 claims abstract description 9
- 238000013528 artificial neural network Methods 0.000 claims abstract description 5
- 238000002347 injection Methods 0.000 claims abstract description 4
- 239000007924 injection Substances 0.000 claims abstract description 4
- 238000009738 saturating Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000002159 abnormal effect Effects 0.000 claims description 2
- 239000000700 radioactive tracer Substances 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims 1
- 239000000243 solution Substances 0.000 abstract description 16
- 239000002245 particle Substances 0.000 abstract description 6
- 238000003900 soil pollution Methods 0.000 abstract description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 6
- 239000012466 permeate Substances 0.000 abstract description 3
- 239000002681 soil colloid Substances 0.000 abstract description 3
- 238000002156 mixing Methods 0.000 abstract description 2
- 239000003344 environmental pollutant Substances 0.000 description 5
- 231100000719 pollutant Toxicity 0.000 description 5
- 230000005684 electric field Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005672 electromagnetic field Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 229910001385 heavy metal Inorganic materials 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000010865 sewage Substances 0.000 description 2
- 238000012271 agricultural production Methods 0.000 description 1
- 238000003915 air pollution Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000002957 persistent organic pollutant Substances 0.000 description 1
- 239000000575 pesticide Substances 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000003802 soil pollutant Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001988 toxicity Effects 0.000 description 1
- 231100000419 toxicity Toxicity 0.000 description 1
- 239000002351 wastewater Substances 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本发明公开了一种非饱和土壤中胶体运动分布三维可视化的方法,装取原状土于圆柱形透明容器;充分饱和土柱,并确保无漏水、渗水现象发生;在土柱表面插上测点,建立测量电路,等待至电路稳定;调配纳米Fe3O4胶体,将胶体浇于土柱表面,浇注前需晃动胶体溶液,使纳米Fe3O4颗粒均匀分布于溶液中,待胶体完全渗入土壤后,静置并保证Fe3O4胶体充分下渗;读取设备测得的电位数据,根据胶体注入时间,将数据分为两部分;使用基于中心差分公式的有限差分法进行电势分布的正演;利用正演结果,使用BP神将网络进行电导率分布的反演,并将反演结果进行三维可视化。本发明操作性强、方法清晰,可得到可靠的土壤胶体运动特点,为提高土壤污染治理效率提供了新方法。
Description
技术领域
本发明属于农业技术和环境生态技术领域,尤其涉及一种非饱和土壤中胶体运动分布三维可视化的方法。
背景技术
随着国家战略目标的改变,近几年来环境问题收到越来越多的关注,与饱受关注的空气污染、水污染相比,土壤污染的研究工作相对较少。
土壤污染物主要是由人类活动带入土壤的,垃圾填埋渗漏、大气重金属沉降、污水排放、农业生产、工业制造等活动最终都会以污水的形式将污染物携带入土壤,而废水中的部分污染物是以胶体的形式存在的,且在运动过程中,胶体会吸附重金属、有机农药以及有机污染物等物质,将这些污染物带离浅层土壤,进入深层土壤及地下水,因此,土壤中胶体的运动在很大程度上增加了污染物的危害范围,并给土壤污染的治理带来难度。研究胶体在土壤中的运动特点可以为污染物的治理提供指导,提高治理效率,降低治理成本。
电法勘探发展历史悠久,现广泛应用于地质勘探工程中,其原理是通过研究与总结天然电场或人工电场的电磁场关于时间与空间的一般分布规律来解决相关问题,但目前的电法装置多应用于大尺度研究,无法在实验室环境中进行研究,实验成本及试验周期较大。
发明内容
本发明所要解决的技术问题是,提供一种非饱和土壤中胶体运动分布三维可视化的方法,操作性强,可得到可靠的土壤胶体运动特点,为提高土壤污染治理效率提供了新方法。
本发明解决其技术问题所采用的技术方案是:提供一种非饱和土壤中胶体运动分布三维可视化的方法,该方法包括以下步骤,步骤一,装取原状土于圆柱形透明容器,充分饱和土柱;步骤二,在土柱表面布置测点电极,建立测量电路,设备上电,等待各测点电压稳定;步骤三,调配纳米Fe3O4胶体,将胶体浇于土柱表面,浇注前需晃动胶体溶液,使纳米Fe3O4颗粒均匀分布于溶液中,待液体完全渗入土壤后,静置并保证Fe3O4胶体充分下渗;步骤四,读取设备测得的电量数据,根据胶体注入时间,将数据分为两部分;步骤五,使用基于中心差分公式的有限差分法对两部分数据进行电势分布正演;步骤六,利用步骤五的正演结果,对两部分数据进行电导率分布反演,将反演结果进行三维可视化。
按上述技术方案,以Fe3O4胶体做为示踪剂。Fe3O4具有良好的导电性,胶体运动过程中滞留在土壤中的Fe3O4胶体颗粒可以显著改变运动轨迹上土壤的电导率,此外Fe3O4无毒,不会造成土壤污染。
按上述技术方案,采用单机供电法作为数据采集方法,电极采用二极排列,供电电极布置在直径两端。
按上述技术方案,所述步骤五中,电势正演方法的具体计算过程如下:
麦克斯韦方程,
将差分公式代入偏微分方程,化简后得:
按上述技术方案,所述步骤四中,读取设备测得的电量数据后转化为.csv文件。
按上述技术方案,所述步骤六中,使用BP神将网络对两部分数据进行电导率分布反演。
本发明产生的有益效果是:操作性强,可得到可靠的土壤胶体运动特点,为提高土壤污染治理效率提供了新方法。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1为本发明实施例中二极排列图;
图2为本发明实施例中电极布置示例图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1、图2所示,本实施例中,选择粒径在100nm以下的纳米Fe3O4,溶于水后形成胶体,用来模拟胶体在土壤中的运动,胶体中的纳米Fe3O4颗粒可以显著改变胶体运动区域的土壤电导率,引起电场分布的改变,以麦克斯韦方程为基础,使用数值模拟计算测点外的其他电压值,反演出胶体运动轨迹,该方法包括测量电路的设计、差分方程的推导以及胶体轨迹反演,包括以下步骤:步骤1:装取原状土于圆柱形透明容器;充分饱和土柱,并确保无漏水、渗水现象发生;步骤2:在土柱表面布置测点电极,建立测量电路,设备上电,等待各测点电压稳定;步骤3:调配纳米Fe3O4胶体,将胶体浇于土柱表面,浇注前需晃动胶体溶液,使纳米Fe3O4颗粒均匀分布于溶液中,待液体完全渗入土壤后,静置并保证Fe3O4胶体充分下渗;步骤4:读取设备测得的电位数据并转化为.csv文件,根据胶体注入时间,将数据分为两部分;步骤5:使用基于中心差分公式的有限差分法对两部分数据进行电势分布正演;步骤6:利用步骤5的正演结果,使用BP神将网络对两部分数据进行电导率分布反演,并将反演结果进行三维可视化。
为了充分使用土罐内的空间,采用单点供电法作为数据采集方法,为充分利用空间采集较多的电压信息,电极采用二极排列,为减小无穷远极的影响,供电两极的距离尽可能的大于正极与测量电极之间的距离,因此供电电极布置在直径两端,该布极方式在众多排列中观测信号较强。
电势正演方法:电势正演的目的是通过测得的电压电流,借助电磁场基本方程计算出电压在样本空间内的分布,其最终要解决的是稳定电流场的边值问题。本发明使用中心差分公式,得到一种更加简洁的计算式,计算过程如下。
基础公式:麦克斯韦方程
差分公式如下:
差分公式代入下式(4),可得:
化简得:
导电率反演则是使用正演过程算得的电势分布,计算对应的导电率分布,但由于投影数据空间小于要反演的图像数据空间,该问题的解必然依赖于正演或其他先验信息才能实际成像,因此数据采集应尽可能地密集、全方位、精确可靠,才能使成像更真实可靠。同时,该问题也是高度非线性问题,若使用数值解法,由于线性化时忽略高于一阶的导数项,会使成像有一定的误差,并且增加了迭代运算的复杂性。BP神经网络在解决非线性问题中具有良好的表现,通过使用“模拟退火”算法或随机梯度下降可避免其陷入局部最优解,因此选用该方法作为反演算法。
分布可视化:经过正演计算后,可得到土样的电导率分布,但直接用所得的电导率进行做图,不能直接反映胶体运动对计算单位内土壤电导率的影响,与之相比,电导率的变化更能直接表现胶体运动对土样的影响。
Δσ=σs-σ0 (10)
步骤四进行数据拆分的目的即在于此,σ0为第一部分计算结果,作为背景电导率,σs为第二部分计算结果,为胶体运动稳定后的电导率,对Δσ进行可视化,即可反映出胶体对土样的影响。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (5)
1.一种非饱和土壤中胶体运动分布三维可视化的方法,其特征在于,该方法包括以下步骤,步骤一,装取原状土于容器,充分饱和土柱;步骤二,在土柱表面布置测点电极,建立测量电路;步骤三,调配纳米Fe3O4胶体,将胶体浇于土柱表面;步骤四,读取设备测得的电量数据,根据胶体注入时间,将数据分为两部分,由背景电位与异常电位组成,
步骤五,使用基于中心差分公式的有限差分法对两部分数据进行电势分布正演;步骤六,利用步骤五的正演结果,对两部分数据进行电导率分布反演,将反演结果进行三维可视化,电势正演方法的具体计算过程如下:
麦克斯韦方程,
将差分公式代入偏微分方程,化简后得:
2.根据权利要求1所述的非饱和土壤中胶体运动分布三维可视化的方法,其特征在于,以Fe3O4胶体做为示踪剂。
3.根据权利要求1或2所述的非饱和土壤中胶体运动分布三维可视化的方法,其特征在于,采用单机供电法作为数据采集方法,电极采用二极排列,供电电极布置在直径两端。
4.根据权利要求1或2所述的非饱和土壤中胶体运动分布三维可视化的方法,其特征在于,所述步骤四中,读取设备测得的电量数据后转化为.csv文件。
5.根据权利要求1或2所述的非饱和土壤中胶体运动分布三维可视化的方法,其特征在于,所述步骤六中,使用BP神将网络对两部分数据进行电导率分布反演。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010071548.0A CN111259558B (zh) | 2020-01-21 | 2020-01-21 | 一种非饱和土壤中胶体运动分布三维可视化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010071548.0A CN111259558B (zh) | 2020-01-21 | 2020-01-21 | 一种非饱和土壤中胶体运动分布三维可视化的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111259558A CN111259558A (zh) | 2020-06-09 |
CN111259558B true CN111259558B (zh) | 2022-06-10 |
Family
ID=70952593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010071548.0A Expired - Fee Related CN111259558B (zh) | 2020-01-21 | 2020-01-21 | 一种非饱和土壤中胶体运动分布三维可视化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111259558B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106683745A (zh) * | 2016-12-16 | 2017-05-17 | 北京先研科技有限责任公司 | 用于改变土壤电阻的改阻剂 |
CN107064645A (zh) * | 2017-04-24 | 2017-08-18 | 南方电网科学研究院有限责任公司 | 一种土壤电阻率和温度对应关系的测量系统及方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013039750A1 (en) * | 2011-09-15 | 2013-03-21 | Planar Solutions, Llc | Homogeneous blending |
CN106443189A (zh) * | 2016-12-05 | 2017-02-22 | 国家电网公司 | 一种接地极极址及其周边土壤电阻率三维探测方法和系统 |
CN108646308B (zh) * | 2018-06-15 | 2019-02-05 | 山东大学 | 一种基于四维电阻率反演的不良地质体灌浆实时监测方法 |
CN110208262A (zh) * | 2018-12-25 | 2019-09-06 | 武汉大学 | 一种测量非饱和土壤中胶体运动分布的方法 |
CN110496854A (zh) * | 2019-08-22 | 2019-11-26 | 上海市环境科学研究院 | 一种基于跨孔高密度电阻率法监控原位注入修复土壤与地下水过程的方法 |
-
2020
- 2020-01-21 CN CN202010071548.0A patent/CN111259558B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106683745A (zh) * | 2016-12-16 | 2017-05-17 | 北京先研科技有限责任公司 | 用于改变土壤电阻的改阻剂 |
CN107064645A (zh) * | 2017-04-24 | 2017-08-18 | 南方电网科学研究院有限责任公司 | 一种土壤电阻率和温度对应关系的测量系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111259558A (zh) | 2020-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103454333B (zh) | 基于PPy/RGO修饰玻碳电极对农药吡虫啉的电化学检测方法 | |
CN107766666B (zh) | 一种基于分数阶差分法的三维时域电磁反常扩散模拟方法 | |
CN115238550B (zh) | 自适应非结构网格的滑坡降雨的地电场数值模拟计算方法 | |
CN204086568U (zh) | 一种超高密度电法数据采集装置 | |
CN111259558B (zh) | 一种非饱和土壤中胶体运动分布三维可视化的方法 | |
CN106969750A (zh) | 磁性液体全方位水平倾角传感器 | |
Touzani et al. | Mapping the pollution plume using the self-potential geophysical method: Case of Oum Azza Landfill, Rabat, Morocco | |
González-Gamboa et al. | Effect of operating parameters on the performance evaluation of benthic microbial fuel cells using sediments from the Bay of Campeche, Mexico | |
Yulong et al. | Fast recognition on shallow groundwater and anomaly analysis using frequency selection sounding method | |
Yang et al. | Output Characteristics of an Electromagnetic–Triboelectric Hybrid Energy Harvester Based on Magnetic Liquid | |
CN106180159A (zh) | 一种基于数据采集卡的污染土壤电动力修复系统电流实时采集装置 | |
Lu et al. | Migration and diffusion of heavy metal Cu from the interior of sediment during wave-induced sediment liquefaction process | |
Shao et al. | Dynamic accumulation of the Quaternary shale biogas in Sanhu Area of the Qaidam Basin, China | |
Lovecchio et al. | Customized Multichannel Measurement System for Microbial Fuel Cell Characterization | |
CN104360246A (zh) | 非均匀场下液体电介质耐电强度的测试装置及方法 | |
Chang‐Ying et al. | A global weak form element free method for direct current resistivity forward simulation | |
Dade-Robertson et al. | Thinking soils: a synthetic biology approach to material-based design computation | |
La Licata et al. | Hydrogeological study of the glacial—Fluvioglacial territory of grandate (Como, Italy) and stochastical modeling of groundwater rising | |
López-Vizcaíno et al. | Two-dimensional modelling approach for electrokinetic water transport in unsaturated kaolinite | |
Anăstăsoaie et al. | Surface Characterization of New Azulene-Based CMEs for Sensing | |
Fang et al. | Mechanics of bio-sediment transport | |
Lu et al. | Detection of Landfill Leachate Leakage Based on ERT and OCTEM | |
Wang et al. | Simulation of an ionic rare earth leaching process based on the darcy law-chemical reaction engineering-transfer of dilute substance coupling | |
Huang et al. | Current pulse characteristic analysis of typical void and point partial discharge based on PIC‐MCC method | |
Freire et al. | Rio Madeira HVDC System: Commissioning of the Ground Electrode for Bipole 2 at Porto Velho |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220610 |