CN116773606B - 一种河湖水下沉积物污染状况获取方法 - Google Patents
一种河湖水下沉积物污染状况获取方法 Download PDFInfo
- Publication number
- CN116773606B CN116773606B CN202310932169.XA CN202310932169A CN116773606B CN 116773606 B CN116773606 B CN 116773606B CN 202310932169 A CN202310932169 A CN 202310932169A CN 116773606 B CN116773606 B CN 116773606B
- Authority
- CN
- China
- Prior art keywords
- resistivity
- sediment
- pollution
- river
- lake
- 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
- 239000013049 sediment Substances 0.000 title claims abstract description 151
- 238000000034 method Methods 0.000 title claims abstract description 48
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 93
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 12
- 231100000719 pollutant Toxicity 0.000 claims abstract description 12
- 238000011109 contamination Methods 0.000 claims description 47
- 238000005259 measurement Methods 0.000 claims description 28
- 238000010291 electrical method Methods 0.000 claims description 15
- 239000000356 contaminant Substances 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 9
- 238000007596 consolidation process Methods 0.000 claims description 8
- 238000005553 drilling Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 abstract 2
- 238000012544 monitoring process Methods 0.000 description 6
- 238000011835 investigation Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000003911 water pollution Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000149 chemical water pollutant Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005183 environmental health Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 239000010842 industrial wastewater Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/06—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a liquid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
- G01C13/008—Surveying specially adapted to open water, e.g. sea, lake, river or canal measuring depth of open water
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Remote Sensing (AREA)
- General Health & Medical Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Development Economics (AREA)
- Hydrology & Water Resources (AREA)
- Immunology (AREA)
- Radar, Positioning & Navigation (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Electrochemistry (AREA)
- Pathology (AREA)
- Educational Administration (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本发明涉及一种河湖水下沉积物污染状况获取方法,包括以下步骤:通过点测河湖水体的电阻率值和水深,初步推测河湖沉积物污染区域电阻率值的范围;利用水上连续高密度电法在圈定的疑似沉积物污染区域内快速地获取连续的沉积物电阻率数据,在疑似污染区域内采集沉积物原状样品,建立电阻率与沉积物污染物浓度之间的定量关系,根据所述定量关系,进而确定精确的沉积物污染电阻率阈值,对疑似污染区进行识别,从而精准定位了沉积物污染区的范围。采用本发明的方法可以对河湖水下沉积物污染实现快速、准确地圈定,可以更科学地制定清淤治理规划,缩短施工周期,节省施工费用。
Description
技术领域
本发明涉及水下污染场地调查技术领域,更具体的,涉及一种河湖水下沉积物污染状况获取方法。
背景技术
水资源作为基础性自然资源和战略性经济资源,支撑着经济社会的发展,是生态环境的重要控制性要素。近年来,绿色环保理念深入人心,人们对水污染带来的生态安全问题和环境健康风险愈发重视。
水污染的来源主要有工业废水、生活污水、农业面源污染、垃圾填埋渗滤液等,被污染的水跟随地下或地表径流进入河流、湖泊,导致河湖水质恶化,破坏生态环境,危害人类健康。污染物通过扩散、沉降等形式在河湖沉积物中聚集,沉积物中的污染物又向河湖中未被污染的水体释放,并随着水的流动不断扩大污染范围。且由于水层的遮挡,水下沉积物污染具有不易察觉、难以定位等特点,若对整个河湖进行全面的清淤治理,则花费巨大,因此亟需对河湖中的水下沉积物污染进行调查,圈定污染分布,评估其对水环境造成的风险,以制定科学的清淤治理规划。
目前河湖水下沉积物污染的调查主要以传统钻探方法为主,该方法在水上操作困难,且需要密集地布点钻探,否则将无法准确圈定污染区域,由于钻探数据具有在空间上非连续的特点,无法整体地掌握河湖水下沉积物污染的分布状况。
发明内容
为了解决上述技术问题,本发明提出了一种河湖水下沉积物污染状况获取方法。
本发明提供了一种河湖水下沉积物污染状况获取方法,包括:
S01:获取河湖水体的电阻率值ρw的最小值ρwmin和最大值ρwmax,推测河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax;
S02:获取河湖水深d的最大值dmax和平均值dave,确定水上连续高密度电法电缆的电极间距j;
S03:利用水上连续高密度电法进行快速测量,测量过程中搭载定位装置和测深装置,同时获取视电阻率数据、位置数据和水深数据,测量船只根据定位装置按照设定航线行驶;
S04:通过所述视电阻率数据对视电阻率进行反演,获得测量区域内沉积物的电阻率ρ,将电阻率ρ与推测河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax进行对比,圈定疑似沉积物污染区域Ω0;
S05:在疑似污染区域Ω0内的设定位置采样沉积物原状样品,测定沉积物样品中的最低污染物浓度cmin和最高污染物浓度cmax,根据对应位置的电阻率ρ,建立电阻率ρ与沉积物污染物浓度c之间的定量关系;
S06:根据所述定量关系确定精确的沉积物污染的电阻率阈值ρmin、ρmax,并根据电阻率阈值对疑似污染区Ω0进行识别,确定沉积物污染区范围Ω。
本方案中,圈定疑似沉积物污染区域Ω0,具体为:
采用水上连续高密度电阻率法进行水体和水下沉积物视电阻率的快速测量,获取河湖探测区域内的视电阻率数据;
对水上连续高密度电法获得的视电阻率数据进行数据处理和反演解译,得到探测区域的电阻率ρ;
对比反演得到的电阻率ρ与推测的河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax,确定疑似沉积物污染区Ω0。
本方案中,根据河湖水体的电阻率的最小值ρwmin和最大值ρwmax,推测河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax,推测公式为:
其中,0.2≤λ≤0.4,λ的取值大小由沉积物的孔隙度和固结程度决定,通常取λ=0.3进行估算。
本方案中,根据河湖水深d的最大值dmax和平均值dave,确定水上连续高密度电法电缆的电极间距j,所述电极间距j的确定公式为:
本方案中,将电阻率ρ与推测河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax进行对比,圈定疑似沉积物污染区域Ω0,具体为:
在圈定疑似沉积物污染区域时,若污染将导致沉积物电阻率的增加,则疑似污染区域的电阻率值ρ0应满足αρsmin≤ρ0≤ρsmax;
若污染将导致沉积物电阻率的减小,则疑似污染区域的电阻率值ρ0应满足ρsmin≤ρ0≤βρsmax;
其中:α,β为控制疑似污染区域圈定范围的系数,α的取值范围为β的取值范围为/>α的取值越大,疑似污染区域圈定的范围越小,β的取值越小,疑似污染区域圈定的范围越小。
本方案中,在疑似污染区域Ω0内的不同位置进行不同深度的钻探取样,以获取沉积物样品中的最低污染物浓度cmin和最高污染物浓度cmax。
本方案中,水上连续高密度电法测量的视电阻率数据的反演过程包含位置信息、水深和水体电阻率的约束。
本方案中,电阻率ρ与沉积物污染物浓度c之间的定量关系由线性拟合公式表述,沉积物污染的电阻率阈值ρmin、ρmax由下述公式确定:
若污染将导致沉积物电阻率的增加,则沉积物污染区域的电阻率阈值为:
若污染将导致沉积物电阻率的减小,则沉积物污染区域的电阻率阈值为:
其中,a为线性拟合公式中的斜率,若污染将导致沉积物电阻率的增加,则a为正值,若污染将导致沉积物电阻率的减小,则a为负值,污染导致沉积物电阻率的变化越大,则|a|越大;b为线性拟合公式中的截距,代表无污染状态下的沉积物电阻率。
本方案中,在疑似污染区Ω0中,若沉积物的电阻率ρ1满足ρmin≤ρ1≤ρmax,则确定为沉积物污染区域Ω。
本方案中,水体电阻率ρw的测量点数量、测量河湖水深d的测量点数量、水上连续高密度电阻率法测量航线的数量根据河湖的长度L、宽度W、面积S确定。
本发明解决了背景技术中存在的缺陷,本发明具备以下有益效果:
1.本发明通过简单的点测河湖水体的电阻率值和水深,便可初步推测河湖沉积物污染区域电阻率值的范围,并确定水上连续高密度电法电缆的电极间距,减少了繁琐地前期勘察准备工作。
2.本发明利用水上连续高密度电法在圈定的疑似沉积物污染区域内快速地获取连续的沉积物电阻率数据,大大提升了沉积物污染调查的效率。
3.本发明通过电阻率剖面初步圈定疑似污染区域,仅在疑似污染区域内采集沉积物原状样品,节省了钻孔和化验的费用。
4.本发明建立了电阻率与沉积物污染物浓度之间的定量关系,根据所述定量关系,进而确定精确的沉积物污染电阻率阈值,对疑似污染区进行识别,从而精准定位了沉积物污染区的范围,使得清淤治理规划的制定更具科学行、针对性,改善了污染区域定位不精确的缺陷,提高了河湖水下沉积物污染探测的效率和精度,能在整体上掌握河湖水下沉积物的污染状况。
附图说明
图1示出了本发明一种河湖水下沉积物污染状况获取方法的流程图;
图2示出了本发明一种河湖水下沉积物污染状况获取方法的示意图;
图3示出了本发明实施例提供的测点与测线分布示意图;
图4示出了本发明实施例中水上连续高密度电法测量示意图;
图5示出了本发明实施例中水下沉积物污染区域圈定示意图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述的其他方式来实施,因此,本发明的保护范围并不受下面公开的具体实施例的限制。
图1、2示出了本发明一种河湖水下沉积物污染状况获取方法的流程图及示意图。
本发明第一方面提供了一种河湖水下沉积物污染状况获取方法,包括:
S01:获取河湖水体的电阻率值ρw的最小值ρwmin和最大值ρwmax,推测河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax;
S02:获取河湖水深d的最大值dmax和平均值dave,确定水上连续高密度电法电缆的电极间距j;
S03:利用水上连续高密度电法进行快速测量,测量过程中搭载定位装置和测深装置,同时获取视电阻率数据、位置数据和水深数据,测量船只根据定位装置按照设定航线行驶;
S04:通过所述视电阻率数据对视电阻率进行反演,获得测量区域内沉积物的电阻率ρ,将电阻率ρ与推测河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax进行对比,圈定疑似沉积物污染区域Ω0;
S05:在疑似污染区域Ω0内的设定位置采样沉积物原状样品,测定沉积物样品中的最低污染物浓度cmin和最高污染物浓度cmax,根据对应位置的电阻率ρ,建立电阻率ρ与沉积物污染物浓度c之间的定量关系;
S06:根据所述定量关系确定精确的沉积物污染的电阻率阈值ρmin、ρmax,并根据电阻率阈值对疑似污染区Ω0进行识别,确定沉积物污染区范围Ω。
在该步骤S01中,河湖水体电阻率ρw的测量根据电阻率水质自动分析仪的行业标准HJ T97-2003执行,河湖水体电阻率的测量点的数量遵循如下原则:
其中:NRσ为河流中水体电阻率测量点的数量,LR为河流的长度,WR为河流的宽度,NLσ为湖泊中水体电阻率测量点的数量,SL为湖泊的面积,式中长度单位为km,面积单位为km2。即,河流中水体电阻率的测量点沿河流宽度方向每0.25km一个,沿河流长度方向每0.1km一个,湖泊中水体电阻率的测量点每0.1km2一个。
在本发明一个较佳实施例中,河流的长度为2km,宽度为0.2km,湖泊的面积为2km2,则河流中的水体电阻率测量点为8个,湖泊中的水体电阻率测量点为20个,如图3所示。
测得河湖水体的电阻率的最小值ρwmin和最大值ρwmax后,根据下述公式可估算河湖沉积物电阻率值的区间范围,求得沉积物电阻率的最小值ρsmin和最大值ρsmax:
其中:0.2≤λ≤0.4,对于孔隙度较大的沉积物,λ取值稍大,对于孔隙度较小的沉积物,λ取值稍小,对于固结良好的沉积物,λ取值稍小,对于未固结的沉积物,λ取值稍大。若沉积物的孔隙度和固结程度未知,则可取λ=0.3进行估算。
若污染将导致沉积物电阻率的增加,则最小值ρsmin表征水下沉积物污染浓度最低或无污染的位置,最大值ρsmax表征水下沉积物污染浓度最高的位置;反之,若污染将导致沉积物电阻率的减小,则最大值ρsmax表征水下沉积物污染浓度最低或无污染的位置,最小值ρsmin表征水下沉积物污染浓度最高的位置。
上述步骤中水体电阻率的测量点的数量和位置可根据实际情况进行调整,涉及河湖水下沉积物电阻率计算的参数可以通过相关水文地质资料获得,或按照经验给出沉积物的参数值。
在该步骤S02中,河湖水深的初步测量可通过声呐测深装置进行,测量点的数量遵循如下原则:
式中:NRd为河流中水深测量点的数量,SR为河流的面积,NLd为湖泊中水深测量点的数量,SL为湖泊的面积,式中长度单位为km,面积单位为km2。即,河流中河流中水深测量点每0.02km2一个,湖泊中水深测量点每0.1km2一个。
在本发明一个较佳实施例中,河流的面积为0.2km2,湖泊的面积为2km2,则河流中水深测量点为10个,湖泊中水深测量点为20个。
具体的,测得水深数据的最大值为dmax,平均值为dave,则水上连续高密度电法电极间距j的选用按照下述公式确定:
在本发明一个较佳实施例中,测得水深数据的最大值dmax为9m,平均值dave为7.5m,则选用的水上连续高密度电法电缆的电极间距j为5m。
上述步骤中水深测量点的数量和位置可根据实际情况进行调整,若水深变化较小,则可适当减少测量点的数量,若水深变化较大,则可适当增加测量点的数量,或在水深变化较大的位置,增加测量点的密度。
在该步骤S03中,利用水上连续高密度电法进行快速测量,测量过程中搭载定位装置和测深装置,同时获取视电阻率数据、位置数据和水深数据,测量船只根据定位装置按照设定航线行驶。
在本发明一个较佳实施例中的水上连续高密度电阻率法测量航线的布设按下述公式确定:
设河湖的长轴方向为i,长度为L,短轴方向为j,宽度为W,
其中:Ni为i方向上的测线数量,Nj为j方向上的测线数量,L为河湖的长度,W为河湖的宽度,x,y为控制测线密度的控制系数,式中长度单位为km。
在本发明一个较佳实施例中,河流的长度为2km,宽度为0.2km,湖泊的长度为2km,宽度为1km。则河流中的水上连续高密度电阻率法测线沿长轴方向每0.05km一条,共4条,沿短轴方向每0.2km一条,共5条,湖泊中的水上连续高密度电阻率法测线沿长轴方向每0.15km一条,共6条,沿宽度方向每0.3km一条,共6条,如图3所示;通过水上连续高密度电法进行快速测量,测量过程中搭载定位装置和测深装置,同时获取视电阻率数据、位置数据和水深数据,水上连续高密度电法测量示意图如图4所示。
在S03中,需要确定测量航线是否经过S01、S02中的水体电阻率测量点和水深测量点附近,若相距较远,需要在本步骤中对测量航线进行适当调整,以便后续圈定疑似沉积物污染区域。此外,在出现电阻率信号异常的区域,需适当增加该区域的测量航线。
具体的,进行水上高连续密度电法探测时,由船只搭载电法仪、定位装置和测深装置,同时获取视电阻率数据、位置数据和水深数据,三种数据的采集位置相对固定,因此可以根据位置坐标进行对应。测量船只根据定位装置按照设定航线行驶,偏移航线时,应及时进行修正,实际行驶轨迹与设计测量航线之间的距离应保持在20m以内。
在该步骤S04中,仪器采集的原始视电阻率数据经过处理,生成计算机可识别处理的文件格式,数据处理时需要对不良数据点进行剔除,并且需要将测得的水深数据d和水体电阻率数据ρw按照一定的格式输入到待反演的文件中,处理后的数据文件通过专业反演软件进行解译,位置信息、水深和水体电阻率数据将会约束视电阻率的反演,以获得更加接近真实情况的电阻率剖面。
将电阻率剖面上的电阻率ρ和S01中初步估算的沉积物电阻率最小值ρsmin和最大值ρsmax进行对比,圈定疑似污染区域Ω0。疑似污染区域Ω0的圈定遵循如下原则:
若污染将导致沉积物电阻率的增加,则疑似污染区域Ω0的电阻率值ρ0应满足如下条件:
αρsmin≤ρ0≤ρsmax
若污染将导致沉积物电阻率的减小,则疑似污染区域Ω0的电阻率值ρ0应满足如下条件:
ρsmin≤ρ0≤βρsmax
其中:α,β为控制疑似污染区域Ω0圈定范围的系数,α的取值范围为β的取值范围为/>α的取值越大,疑似污染区域Ω0圈定的范围越小,β的取值越小,疑似污染区域Ω0圈定的范围也越小。
具体的,S01中的水体电阻率的测量和沉积物电阻率的估算、S03中视电阻率数据的采集、S04中反演得到电阻率数据,均可能存在测量误差和计算误差,因此使用参数α,β调整本步骤中疑似污染区域Ω0的圈定范围,结合实际的水文地质背景和前期资料,使疑似污染区域Ω0的圈定更加真实合理。
在该步骤S05中,疑似污染区域Ω0内钻取沉积物原状样品位置的设定遵循如下原则:
Nip≥Ni,Njp≥Nj
其中,Nip为i方向测线上的沉积物原状样品采集点数量,Njp为j方向测线上的沉积物原状样品采集点数量。
沉积物原状样品采集点应大于等于测线的数量,每条测线上至少有一个采集点,且采集点的位置应位于测线上。若测线的电阻率剖面显示有多个疑似沉积物污染的电阻率异常位置,则需在异常位置加密采集点。
具体的,测定沉积物样品中的污染物浓度c后,电阻率剖面中的沉积物电阻率ρ与污染物浓度c之间的定量关系可以用线性拟合公式y=ax+b进行定量化表述,即:
ρ=ac+b
其中:a为线性拟合公式中的斜率,若污染将导致沉积物电阻率的增加,则a为正值,若污染将导致沉积物电阻率的减小,则a为负值,污染导致沉积物电阻率的变化越大,则|a|越大;b为线性拟合公式中的截距,代表无污染状态下的沉积物电阻率。
根据不同位置和不同深度处的沉积物电阻率和污染物浓度的关系,可以对上述拟合公式的参数a,b进行确定。
在该步骤S06中,测得沉积物样品中的最低污染物浓度cmin和最高污染物浓度cmax,则沉积物污染电阻率阈值ρmin、ρmax的确定遵循如下原则:
若污染将导致沉积物电阻率的增加,则沉积物污染区域的电阻率阈值为:
若污染将导致沉积物电阻率的减小,则沉积物污染区域的电阻率阈值为:
在疑似污染区Ω0中,若沉积物的电阻率ρ1满足ρmin≤ρ1≤ρmax,则确定为沉积物污染区域Ω。所述的沉积物污染区域Ω将在电阻率剖面中圈定,以更加直观地显示沉积物污染区域的分布信息,如图5所示。
构建河湖水下沉积物污染信息数据集,将目标区域内水下沉积物污染状况的历史监测信息存入数据库中,根据预设时间内的历史监测信息分析水下沉积物污染对应的污染趋势,获取污染平均变化速率作为基准速率,获取目标区域水下沉积物污染状况的当前监测信息与上一监测时间点的历史监测信息的偏差,并计算当前监测时间的污染变化速率,与基准速率进行对比,若速率偏差大于预设阈值,则生成水下沉积物污染预警,根据污染预警进行相关污染源的异常排放检测或者新增污染源的检测。
采用本发明的方法,设计了一套完整的河湖水下沉积物污染状况获取方法,为后续制定科学的清淤治理规划提供详细的参考依据,提高了河湖水下沉积物污染探测的效率和精度,只需在疑似污染区域内采集少量沉积物样品的污染信息,结合水上高密度电法的快速测量,就能在整体上掌握河湖水下沉积物的污染状况。
在本申请所提供的几个实施例中,应该理解到,所揭露的设备和方法,可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,如:多个单元或组件可以结合,或可以集成到另一个***,或一些特征可以忽略,或不执行。另外,所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口,设备或单元的间接耦合或通信连接,可以是电性的、机械的或其它形式的。
上述作为分离部件说明的单元可以是、或也可以不是物理上分开的,作为单元显示的部件可以是、或也可以不是物理单元;既可以位于一个地方,也可以分布到多个网络单元上;可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。
另外,在本发明各实施例中的各功能单元可以全部集成在一个处理单元中,也可以是各单元分别单独作为一个单元,也可以两个或两个以上单元集成在一个单元中;上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:移动存储设备、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
或者,本发明上述集成的单元如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括:移动存储设备、ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (8)
1.一种河湖水下沉积物污染状况获取方法,其特征在于,包括以下步骤:
S01:获取河湖水体的电阻率值ρw的最小值ρwmin和最大值ρwmax,推测河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax;
S02:获取河湖水深d的最大值dmax和平均值dave,确定水上连续高密度电法电缆的电极间距j;
S03:利用水上连续高密度电法进行快速测量,测量过程中获取视电阻率数据、位置数据和水深数据;
S04:通过所述视电阻率数据对视电阻率进行反演,获得测量区域内沉积物的电阻率ρ,将电阻率ρ与推测河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax进行对比,圈定疑似沉积物污染区域Ω0;
S05:在疑似污染区域Ω0内的设定位置采样沉积物原状样品,测定沉积物样品中的最低污染物浓度cmin和最高污染物浓度cmax,根据对应位置的电阻率ρ,建立电阻率ρ与沉积物污染物浓度c之间的定量关系;
S06:根据所述定量关系确定精确的沉积物污染的电阻率阈值ρmin、ρmax,并根据电阻率阈值对疑似污染区Ω0进行识别,确定沉积物污染区范围Ω;
根据河湖水体的电阻率的最小值ρwmin和最大值ρwmax,推测河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax,推测公式为:
其中,0.2≤λ≤0.4,λ的取值大小由沉积物的孔隙度和固结程度决定,对于孔隙度较大的沉积物,λ取值稍大,对于孔隙度较小的沉积物,λ取值稍小,对于固结良好的沉积物,λ取值稍小,对于未固结的沉积物,λ取值稍大。若沉积物的孔隙度和固结程度未知,则可取λ=0.3进行估算;
将电阻率ρ与推测河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax进行对比,圈定疑似沉积物污染区域Ω0,具体为:
在圈定疑似沉积物污染区域时,若污染将导致沉积物电阻率的增加,则疑似污染区域的电阻率值ρ0应满足αρsmin≤ρ0≤ρsmax;
若污染将导致沉积物电阻率的减小,则疑似污染区域的电阻率值ρ0应满足ρsmin≤ρ0≤βρsmax;
其中:α,β为控制疑似污染区域圈定范围的系数,α的取值范围为β的取值范围为/>α的取值越大,疑似污染区域圈定的范围越小,β的取值越小,疑似污染区域圈定的范围越小。
2.根据权利要求1所述的一种河湖水下沉积物污染状况获取方法,其特征在于,圈定疑似沉积物污染区域Ω0,具体为:
采用水上连续高密度电法进行水体和水下沉积物视电阻率的快速测量,获取河湖探测区域内的视电阻率数据;
对水上连续高密度电法获得的视电阻率数据进行数据处理和反演解译,得到探测区域的电阻率ρ;
对比反演得到的电阻率ρ与推测的河湖沉积物污染区域电阻率值的最小值ρsmin和最大值ρsmax,确定疑似沉积物污染区Ω0。
3.根据权利要求1所述的一种河湖水下沉积物污染状况获取方法,其特征在于,根据河湖水深d的最大值dmax和平均值dave,确定水上连续高密度电法电缆的电极间距j,所述电极间距j的确定公式为:
4.根据权利要求1所述的一种河湖水下沉积物污染状况获取方法,其特征在于,在疑似污染区域Ω0内的不同位置进行不同深度的钻探取样,以获取沉积物样品中的最低污染物浓度cmin和最高污染物浓度cmax。
5.根据权利要求1所述的一种河湖水下沉积物污染状况获取方法,其特征在于,水上连续高密度电法测量的视电阻率数据的反演过程包含位置信息、水深和水体电阻率的约束。
6.根据权利要求1所述的一种河湖水下沉积物污染状况获取方法,其特征在于,电阻率ρ与沉积物污染物浓度c之间的定量关系由线性拟合公式表述,沉积物污染的电阻率阈值ρmin、ρmax由下述公式确定:
若污染将导致沉积物电阻率的增加,则沉积物污染区域的电阻率阈值为:
若污染将导致沉积物电阻率的减小,则沉积物污染区域的电阻率阈值为:
其中,a为线性拟合公式中的斜率,若污染将导致沉积物电阻率的增加,则a为正值,若污染将导致沉积物电阻率的减小,则a为负值,污染导致沉积物电阻率的变化越大,则|a|越大;b为线性拟合公式中的截距,代表无污染状态下的沉积物电阻率。
7.根据权利要求1所述的一种河湖水下沉积物污染状况获取方法,其特征在于,在疑似污染区Ω0中,若沉积物的电阻率ρ1满足ρmin≤ρ1≤ρmax,则确定为沉积物污染区域Ω。
8.根据权利要求1所述的一种河湖水下沉积物污染状况获取方法,其特征在于,水体电阻率ρw的测量点数量、测量河湖水深d的测量点数量、水上连续高密度电法测量航线的数量根据河湖的长度L、宽度W、面积S确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310932169.XA CN116773606B (zh) | 2023-07-27 | 2023-07-27 | 一种河湖水下沉积物污染状况获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310932169.XA CN116773606B (zh) | 2023-07-27 | 2023-07-27 | 一种河湖水下沉积物污染状况获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116773606A CN116773606A (zh) | 2023-09-19 |
CN116773606B true CN116773606B (zh) | 2024-03-08 |
Family
ID=88008240
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310932169.XA Active CN116773606B (zh) | 2023-07-27 | 2023-07-27 | 一种河湖水下沉积物污染状况获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116773606B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117890553B (zh) * | 2024-03-14 | 2024-06-11 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种水文水资源监控*** |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR676075A (fr) * | 1928-08-03 | 1930-02-18 | Méthode électrique et appareil pour la détermination de la concentration du sang,du lait et d'autres suspensions | |
US5652519A (en) * | 1994-01-13 | 1997-07-29 | Jesse G. Robison | Method and apparatus for measuring pollutants contaminating earth formation |
KR20040029828A (ko) * | 2002-10-02 | 2004-04-08 | 김필선 | 다중전극탐침봉 및 이를 이용한 시추공 방위 비저항탐사방법 |
CN111664887A (zh) * | 2020-05-05 | 2020-09-15 | 中国海洋大学 | 基于电阻率探杆的海底浮泥层动态变化原位观测方法 |
CN114415242A (zh) * | 2022-01-26 | 2022-04-29 | 河北地质大学 | 一种用于土壤及地下水污染调查的地球物理探测设备 |
CN115096947A (zh) * | 2022-07-18 | 2022-09-23 | 东南大学 | 一种基于电阻率法的NAPLs污染物浓度反演方法 |
CN115290702A (zh) * | 2022-08-22 | 2022-11-04 | 合肥工业大学 | 一种升尺度裂隙污染物运移监测装置及监测方法 |
CN115933003A (zh) * | 2022-11-23 | 2023-04-07 | 山东省生态环境规划研究院 | 一种平原地形垃圾填埋场地下水污染状况获取方法 |
CN116310771A (zh) * | 2023-03-21 | 2023-06-23 | 中国电力工程顾问集团华北电力设计院有限公司 | 一种耦合遥感的面源污染源识别方法 |
-
2023
- 2023-07-27 CN CN202310932169.XA patent/CN116773606B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR676075A (fr) * | 1928-08-03 | 1930-02-18 | Méthode électrique et appareil pour la détermination de la concentration du sang,du lait et d'autres suspensions | |
US5652519A (en) * | 1994-01-13 | 1997-07-29 | Jesse G. Robison | Method and apparatus for measuring pollutants contaminating earth formation |
KR20040029828A (ko) * | 2002-10-02 | 2004-04-08 | 김필선 | 다중전극탐침봉 및 이를 이용한 시추공 방위 비저항탐사방법 |
CN111664887A (zh) * | 2020-05-05 | 2020-09-15 | 中国海洋大学 | 基于电阻率探杆的海底浮泥层动态变化原位观测方法 |
CN114415242A (zh) * | 2022-01-26 | 2022-04-29 | 河北地质大学 | 一种用于土壤及地下水污染调查的地球物理探测设备 |
CN115096947A (zh) * | 2022-07-18 | 2022-09-23 | 东南大学 | 一种基于电阻率法的NAPLs污染物浓度反演方法 |
CN115290702A (zh) * | 2022-08-22 | 2022-11-04 | 合肥工业大学 | 一种升尺度裂隙污染物运移监测装置及监测方法 |
CN115933003A (zh) * | 2022-11-23 | 2023-04-07 | 山东省生态环境规划研究院 | 一种平原地形垃圾填埋场地下水污染状况获取方法 |
CN116310771A (zh) * | 2023-03-21 | 2023-06-23 | 中国电力工程顾问集团华北电力设计院有限公司 | 一种耦合遥感的面源污染源识别方法 |
Non-Patent Citations (5)
Title |
---|
Monitoring underground water quality based on high density resistivity method;Yanyan Xu;Internatuinal Conference on Intelligent Earth Observing and Applications;第9808卷;98081T1-98081T12 * |
刘国栋 ; 罗强 ; 张沛垚 ; .高密度电阻率法在城市环境监测中的应用.地质装备.2009,(02),1246-1254. * |
刘汉乐 ; 周启友 ; 吴华桥 ; .轻非水相液体污染过程的高密度电阻率成像法室内监测.地球物理学报.2008,(04),31-36. * |
轻非水相液体污染过程的高密度电阻率成像法室内监测;刘汉乐;周启友;吴华桥;;地球物理学报(04);31-36 * |
高密度电阻率法在城市环境监测中的应用;刘国栋;罗强;张沛垚;;地质装备(02);1246-1254 * |
Also Published As
Publication number | Publication date |
---|---|
CN116773606A (zh) | 2023-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116773606B (zh) | 一种河湖水下沉积物污染状况获取方法 | |
CN111539867A (zh) | 一种针对地下水氯代烃污染的环境调查方法 | |
CN106484997B (zh) | 一种基于克里金插值的水岸带淤泥厚度计算及出图方法 | |
CN105203599A (zh) | 一种场地污染土的快速诊断方法 | |
CN105181758A (zh) | 一种基于电阻率物探技术的污染土快速诊断方法 | |
CN113203691A (zh) | 一种基于小波分析的水质污染物溯源方法 | |
CN107909192A (zh) | 土壤重金属含量的估测方法及装置 | |
CN114184760B (zh) | 基于水位在线观测的污染物通量监测方法、终端及存储介质 | |
CN114002381A (zh) | 一种大气污染溯源扩散分析方法及装置 | |
KR101570518B1 (ko) | 3차원 지도를 이용한 매립 폐기물 조사 방법 | |
CN115933003A (zh) | 一种平原地形垃圾填埋场地下水污染状况获取方法 | |
CN115420778B (zh) | 一种浅层土壤污染快速调查定位方法及*** | |
CN117113641A (zh) | 一种平原水库库盆渗漏状况的获取方法及*** | |
CN115541834A (zh) | 地下水污染源位置识别方法及装置、电子设备、存储介质 | |
CN116165718A (zh) | 土石坝渗流安全的测定方法、装置、电子设备及介质 | |
CN114117757A (zh) | 地下水污染物溯源方法、装置、设备及存储介质 | |
CN109060858A (zh) | 一种土壤重金属浓度空间分布预测准确性的定量评价方法 | |
CN115078190A (zh) | 一种悬浮体现场激光粒度数据处理方法及其装置 | |
CN109839395B (zh) | 土壤污染物分析方法、装置、设备和*** | |
CN105301660A (zh) | 一种用于污染土辨识的点式物探方法 | |
RU2319139C2 (ru) | Способ определения дефектов гидроизоляционного покрытия и коррозионных повреждений наружных поверхностей подземных и подводных трубопроводов | |
CN115983526A (zh) | 一种土壤有机物污染监测方法 | |
CN108229849A (zh) | 确定小型河道污染物降解系数不确定性及其风险程度的方法 | |
CN116953186A (zh) | 河道污染溯源方法、装置、设备及存储介质 | |
JP2005300260A (ja) | 汚染土壌の調査方法 |
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 |