CN116774295A - Method for determining solid waste landfill position based on high-density resistivity method - Google Patents
Method for determining solid waste landfill position based on high-density resistivity method Download PDFInfo
- Publication number
- CN116774295A CN116774295A CN202310576239.2A CN202310576239A CN116774295A CN 116774295 A CN116774295 A CN 116774295A CN 202310576239 A CN202310576239 A CN 202310576239A CN 116774295 A CN116774295 A CN 116774295A
- Authority
- CN
- China
- Prior art keywords
- resistivity
- landfill
- inversion
- solid waste
- area
- 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
- 238000000034 method Methods 0.000 title claims abstract description 94
- 239000002910 solid waste Substances 0.000 title claims abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims description 18
- 230000002159 abnormal effect Effects 0.000 claims description 15
- 238000001514 detection method Methods 0.000 claims description 14
- 239000002689 soil Substances 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 11
- 238000002955 isolation Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 4
- 238000005553 drilling Methods 0.000 claims description 4
- 239000012267 brine Substances 0.000 claims description 3
- 238000010291 electrical method Methods 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims description 3
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 claims description 3
- 238000012800 visualization Methods 0.000 claims description 3
- 238000013102 re-test Methods 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 8
- 238000011160 research Methods 0.000 description 7
- 238000005755 formation reaction Methods 0.000 description 5
- 239000011780 sodium chloride Substances 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 238000004080 punching Methods 0.000 description 4
- 239000004927 clay Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000011835 investigation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000005684 electric field Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000002893 slag Substances 0.000 description 1
Classifications
-
- 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
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
Abstract
The invention discloses a method for determining a solid waste landfill position based on a high-density resistivity method, which comprises the following steps of: the method has the advantages of high observation precision, large data acquisition amount, rich geological information, high production efficiency and the like, solves the problems of inaccurate parameter model and the like in the data processing in the inversion process through the pre-judgment and the data processing, improves the reliability of the data, and further accurately determines the solid waste landfill area.
Description
Technical Field
The invention relates to the technical field of environmental protection, in particular to a method for determining a solid waste landfill position based on a high-density resistivity method.
Background
The high-density resistivity method belongs to the category of direct current resistivity methods, is an exploration method developed on the basis of conventional electrical exploration, is still based on the electrical difference of rock and soil bodies, and is used for researching the distribution rule of underground conduction current under the action of an applied electric field. The resistivity method can quickly and efficiently realize the detection of the polluted site, and the detection result can effectively reflect the pollution distribution information of the polluted site, thereby realizing the preliminary judgment of the pollution degree and the pollution range of the polluted site. Compared with the traditional electrical method, the high-density resistivity method has the advantages of high observation precision, large data acquisition amount, rich geological information, high production efficiency and the like. The high-density resistivity method can complete the longitudinal and transverse two-dimensional exploration process, not only can reflect the electrical change condition of a rock-soil body at a certain depth in the underground along the horizontal direction, but also can provide the electrical change information of the formation lithology along the longitudinal direction, so that the method has the comprehensive detection capability of an electrical profile method and an electrical sounding method.
The high-density resistivity method is a direct current electric exploration method which is based on resistivity difference between an underground research target body and surrounding mediums, artificially establishes an underground stable direct current electric field, arranges a plurality of electrodes in advance, adopts a preset device arrangement for detection, and researches the distribution rule of conductive current in a certain underground range, thereby researching and solving the related geological problems. The method is characterized in that: when the device works in the field, a plurality of electrodes are controlled by the connection of the electrode converters and are connected in series on the same arrangement, and the free switching of various different devices in the resistivity method is realized according to the arrangement of preset devices, so that various apparent resistivity parameters under the condition of various devices can be measured by one-time electrode distribution.
Disclosure of Invention
Accordingly, the present invention provides a method for determining a location of a solid waste landfill based on a high density resistivity method, comprising the steps of:
s1, resistivity pre-judgment
Adopting two modes of ground sampling and drilling sampling, testing resistivity values of known landfill and stratum around the landfill in soil of a region to be confirmed, judging whether the known landfill and the stratum have obvious electrical difference according to the difference range of the resistivity values of the known landfill and the stratum around the known landfill, if so, further exploring the landfill by using a high-density resistivity method, otherwise, not researching the landfill;
s2, arranging the measuring lines:
arranging a plurality of high-density resistivity measuring lines in a solid waste landfill area to be confirmed;
s3, arranging electrodes:
selecting a plurality of measuring points with equal intervals on each high-density resistivity measuring line, arranging electrodes on the measuring points, and then connecting the electrodes to a program-controlled multi-channel electrode conversion switch by using a multi-core cable;
s4, collecting data:
sending the measurement signal generated by the electrode into a micro-electromechanical measuring instrument through an electrode change-over switch, and storing the measurement result into a random access memory; the data of the random access memory is firstly transmitted to a computer through a network, and the WDAFC.EXE software of the computer is utilized to convert the original data format into a format which can be read by Res2dinv software;
s5, processing data:
pre-processing the raw data using Res2dinv processing software, the pre-processing comprising: splicing and removing mutation points; then, establishing a model, and carrying out forward calculation by adopting Res2dmod, wherein the forward calculation is carried out to obtain a potential value and a resistivity profile at the grid node; the grid is defined by Res2dmod software, and then a least square method is adopted to perform two-dimensional inversion calculation to obtain an inversion result of the forward model;
in the inversion process, all models select 3-5 times of inversion iteration results; obtaining an RMS error value in the multi-iteration inversion result; when the obtained RMS error values are smaller than 0.2, calculating variance values of a plurality of groups of RMS error values, and repeating the inversion process when the RMS error values are larger than or equal to 0.2;
judging whether the variance value is more than 15%, if so, judging that the inversion result is reasonable, and if so, judging that the inversion result is unreasonable; when the inversion result is reasonable, the inversion result is derived, a multilayer inversion resistivity section diagram is drawn by using Res2dinvx64 software, and meanwhile, coordinate conversion is carried out; repeating the inversion process when the inversion result is unreasonable;
s6, confirming area
Merging the inversion resistivity section data of each layer in the inversion result to obtain full-region inversion resistivity three-dimensional space data, and carrying out three-dimensional visualization on the full-region inversion resistivity three-dimensional space data according to RockWorks processing software to obtain an inversion resistivity plane graph and a three-dimensional stereogram; and searching a common area with the largest resistance difference of the inversion resistivity section view, the inversion resistivity plane view and the three-dimensional stereogram to obtain the position coordinates of the abnormal resistance area, and determining the position coordinates as a solid waste landfill area.
Description: through the arrangement of the method, the solid waste landfill area can be accurately determined, and the high-density resistivity method has the advantages of high observation precision, large data acquisition amount, rich geological information, high production efficiency and the like. The high-density resistivity method can complete a longitudinal and transverse two-dimensional exploration process, not only can reflect the electrical change condition of a rock-soil body at a certain depth in the underground along the horizontal direction, but also can provide the electrical change information of stratum lithology along the longitudinal direction, and has the comprehensive detection capability of an electrical profile method and an electrical sounding method; through pre-judgment and data processing, the problems that a parameter model is not accurate enough and the like exist in the data processing of the inversion process are solved, the reliability of data is improved, and therefore the solid waste landfill area is accurately determined; the reliability of the data is comprehensively judged through the combination of the RMS error value and the variation stability of the RMS error value, so that the phenomenon that an RMS model can sometimes show large unrealistic variation on the model resistivity value is avoided.
Further, in step S1, the method for determining whether the difference between the resistivity values of the known landfill and the surrounding stratum has a significant electrical difference is:
according to the formulaWhether the calculated value is larger than 0.2 or not is judged to have obvious electrical property difference, otherwise, the calculated value does not have obvious electrical property difference; where R0 is the resistivity value of the surrounding formation and R1 is the resistivity value of the known landfill.
Description: the reliability of the high-density resistivity method for targeted research can be rapidly and effectively judged by the judging method, and the conditions of redundant test and inaccurate results are avoided.
Further, in step S4, a Chongqing Pentium WGMD-9 high density electrical method system is adopted to collect measurement signals, and the collection parameters are as follows: the number of the electrodes is 60-180, the electrode distance is 1m, the maximum isolation coefficient is 19-30, and the power supply voltage is 450V.
Description: through the use of the system, the acquisition process can be completed well, and meanwhile, the obtained data is accurate and comprehensive.
Further, the measured voltage in step S4 is 5 to 10mV.
Description: the observation precision of the WDA-1 host can reach +/-0.2% when the measurement voltage is greater than 5 mv; in order to ensure that the observed data has an optimal signal to noise ratio, the measured voltage is kept in a range above 5mv as much as possible.
Further, in step S5, after the calculation is completed, the quality of the profile curve or the resistivity profile is checked, whether there are bad points or not and the bad point area are judged, and when there are a few bad points or the bad point area and the number of the bad points is less than or equal to 5 and the area of the area is less than or equal to 1%, the single-point measurement function retest is started; and retesting the whole section when the number of the bad points is more and the number of the bad points is more than 5 and the area of the area is more than 1%.
Description: through the setting of above-mentioned inspection process, can make numerical value and section value that obtains more accurate, have good practicality.
Further, before the electrode is arranged in the step S3, holes are punched in the ground, soil is filled in the holes, and then saline water is poured until the grounding resistance of the electrode is reduced to be less than 2KΩ.
Description: the field grounding resistance measurement results show that the grounding resistance can be effectively controlled below 2KΩ through punching, mud plugging and water pouring, the grounding condition can be good, and the electrode grounding condition is obviously improved through punching, mud plugging, water pouring and other modes, so that a good foundation is laid for guaranteeing the quality of original data.
Further, before the data acquisition in step S4, a parameter selection test is performed, and data of the temperature receiving device, the temperature applying device and the dipole device are acquired simultaneously for each section, and parameters of the device, the device form, the isolation coefficient and the power supply voltage are determined according to inversion results of the three devices and position coordinates of the abnormal resistance areas.
Description: through a parameter selection test, determining parameters such as device form, isolation coefficient, power supply voltage and the like according to inversion results of three devices and position coordinates of a resistance abnormal region; the detection results can be mutually supplemented and mutually verified, and the data quality and the exploration effect are fully ensured.
Further, in step S5, the forward model is: the number of the electrodes is 81, the electrode distance is 1m, the low-resistance different area is 0-3 omega, and the surrounding stratum is 17-20 omega; the number of layers of resistivity data acquisition was 20.
Description: the data acquisition effect of the parameters on solid waste is good.
Further, after confirming the area in step S6, the planar area of the solid waste landfill is estimated, and the square quantity is estimated in the detection depth range.
Description: through the calculation, the position of the solid waste landfill area can be clarified, and the subsequent evaluation of the land and the corresponding measures can be conveniently carried out.
The beneficial effects of the invention are as follows:
(1) By the arrangement of the method, the solid waste landfill area can be accurately determined, and the high-density resistivity method has the advantages of high observation precision, large data acquisition amount, rich geological information, high production efficiency and the like. The high-density resistivity method can complete a longitudinal and transverse two-dimensional exploration process, not only can reflect the electrical change condition of a rock-soil body at a certain depth in the underground along the horizontal direction, but also can provide the electrical change information of stratum lithology along the longitudinal direction, and has the comprehensive detection capability of an electrical profile method and an electrical sounding method; through pre-judgment and data processing, the problems that a parameter model is not accurate enough and the like exist in the data processing of the inversion process are solved, the reliability of data is improved, and therefore the solid waste landfill area is accurately determined; the reliability of the data is comprehensively judged through the combination of the RMS error value and the variation stability of the RMS error value, so that the phenomenon that an RMS model can sometimes show large unrealistic variation on the model resistivity value is avoided.
(2) The reliability of the high-density resistivity method for targeted research can be rapidly and effectively judged through the pre-judging method, and the conditions of redundant test and inaccurate results are avoided; through a parameter selection test, determining parameters such as device form, isolation coefficient, power supply voltage and the like according to inversion results of three devices and position coordinates of a resistance abnormal region; the detection results can be mutually supplemented and mutually verified, so that the data quality and the exploration effect are fully ensured; the grounding resistance can be effectively controlled below 2KΩ through punching, mud plugging and water casting, so that the grounding condition is good, and a good foundation is laid for guaranteeing the quality of the original data.
Drawings
FIG. 1 is a survey line arrangement of the present invention;
FIG. 2 is a flow chart illustrating the high density resistivity process of the present invention;
FIGS. 3-4 are cross-sectional views of resistivity of three devices along line T1 in the present invention;
FIG. 5 is a plot of inversion resistivity and geological interpretation of the high density resistivity method along line T1 in the region of the invention, wherein (a) is the inversion resistivity profile and (b) is the geological interpretation profile;
FIG. 6 is a plot of inversion resistivity and geological interpretation of the high density resistivity method along line T2 of the region of the invention, wherein (a) is an inversion resistivity profile and (b) is a geological interpretation profile;
FIG. 7 is a plot of inversion resistivity and geological interpretation of the high density resistivity method along line T3 in the region of the invention, wherein (a) is an inversion resistivity profile and (b) is a geological interpretation profile;
FIG. 8 is a plot of inversion resistivity and geological interpretation of the high density resistivity method along line T4 in the region of the invention, wherein (a) is an inversion resistivity profile and (b) is a geological interpretation profile;
FIG. 9 is a plot of inversion resistivity and geological interpretation of the high density resistivity method along line T5 in the region of the invention, wherein (a) is an inversion resistivity profile and (b) is a geological interpretation profile;
FIG. 10 is a three-dimensional plot of the inversion resistivity of the high density resistivity method of the invention.
Detailed Description
The invention will be described in further detail with reference to the following embodiments to better embody the advantages of the invention.
Example 1:
a method for determining a location of a solid waste landfill based on a high density resistivity method, comprising the steps of:
s1, resistivity pre-judgment
Adopting two modes of ground sampling and drilling sampling, testing resistivity values of known landfill and stratum around the landfill in soil of a region to be confirmed, judging whether the known landfill and the stratum have obvious electrical difference according to the difference range of the resistivity values of the known landfill and the stratum around the known landfill, if so, further exploring the landfill by using a high-density resistivity method, otherwise, not researching the landfill;
the method for judging whether the resistivity value difference between the known landfill and the stratum around the known landfill has obvious electrical property difference is as follows:
according to the formulaWhether the calculated value is larger than 0.2 or not is judged to have obvious electrical property difference, otherwise, the calculated value does not have obvious electrical property difference; wherein R0 is the resistivity value of the surrounding formation and R1 is the resistivity value of the known landfill;
s2, arranging the measuring lines:
arranging a plurality of high-density resistivity measuring lines in a solid waste landfill area to be confirmed;
s3, arranging electrodes:
selecting a plurality of measuring points with equal intervals on each high-density resistivity measuring line, punching holes on the ground where the measuring points are located, filling local soil in the holes, and pouring 25g/L sodium chloride saline until the electrode grounding resistance is reduced to 1.8KΩ; arranging electrodes on the measuring points, and then connecting the electrodes to a program-controlled multi-path electrode conversion switch by using a multi-core cable;
s4, collecting data:
sending the measurement signal generated by the electrode into a micro-electromechanical measuring instrument through an electrode change-over switch, and storing the measurement result into a random access memory; the data of the random access memory is firstly transmitted to a computer through a network, and the WDAFC.EXE software of the computer is utilized to convert the original data format into a format which can be read by Res2dinv software;
wherein, a Chongqing Pentium WGMD-9 high-density electrical system is adopted to collect measurement signals and collect data of a temperature device; the acquisition parameters are as follows: the number of the electrodes is 81, the electrode distance is 1m, the maximum isolation coefficient is 30, and the power supply voltage is 450V; the measurement voltage was 7mV; collecting the layer number 20;
s5, processing data:
pre-processing the raw data using Res2dinv processing software, the pre-processing comprising: splicing and removing mutation points; then, a model is built, and Res2dmod is adopted to carry out forward calculation, wherein the forward model is as follows: the number of the electrodes is 81, the electrode distance is 1m, the low-resistance different area is 0-3 omega, and the surrounding stratum is 17-20 omega; the number of the resistivity data acquisition layers is 20; forward calculation is carried out to obtain a potential value and a resistivity profile at the grid node; the grid is defined by Res2dmod software, and then a least square method is adopted to perform two-dimensional inversion calculation to obtain an inversion result of the forward model;
in the inversion process, selecting 5 inversion iteration results for all models; obtaining an RMS error value in the inversion result of 5 iterations; when the obtained RMS error values are smaller than 0.2, calculating variance values of a plurality of groups of RMS error values, and repeating the inversion process when the RMS error values are larger than or equal to 0.2;
judging whether the variance value is more than 15%, if so, judging that the inversion result is reasonable, and if so, judging that the inversion result is unreasonable; when the inversion result is reasonable, the inversion result is derived, a multilayer inversion resistivity section diagram is drawn by using Res2dinvx64 software, and meanwhile, coordinate conversion is carried out; repeating the inversion process when the inversion result is unreasonable;
s6, confirming area
Merging the inversion resistivity section data of each layer in the inversion result to obtain full-region inversion resistivity three-dimensional space data, and carrying out three-dimensional visualization on the full-region inversion resistivity three-dimensional space data according to RockWorks processing software to obtain an inversion resistivity plane graph and a three-dimensional stereogram; searching a common area with the largest resistance difference of the inversion resistivity section view, the inversion resistivity plane view and the three-dimensional stereogram, obtaining the position coordinates of the abnormal resistance area, determining the position coordinates as a solid waste landfill area, estimating the plane area of the solid waste landfill area, and estimating the square quantity of the solid waste landfill area in the detection depth range.
Example 2
This example differs from example 1 in that in step S2, 30g/L sodium chloride brine is poured until the electrode ground resistance is reduced to 1.8KΩ; in the step S4, data of a temperature application device are adopted, and acquisition parameters are as follows: 180 electrodes, 1m electrode distance, 19 maximum isolation coefficient and 450V supply voltage; the measured voltage in step S3 was 10mV.
Example 3
This example differs from example 1 in that in step S2, 20g/L sodium chloride brine is poured until the electrode ground resistance is reduced to 1.8KΩ; in the step S4, data of a dipole device are adopted, and acquisition parameters are as follows: the number of the electrodes is 60, the electrode distance is 1m, the maximum isolation coefficient is 24, and the power supply voltage is 450V; the measured voltage in step S3 was 5mV.
Experimental example
1. Testing resistivity values of known landfill and formations around the landfill, and adopting two modes of ground sampling and drilling sampling; after field testing, the following results were obtained:
the known landfill is in the form of white slag soil, and the average value change range of the resistivity is 7.16-9.04 omega.m under the condition that the water content is more and the soil is more open; the known landfill has a resistivity change range of 9.76-11.6Ω·m when the soil is dry and the water content is low; the main soil in the factory is clay layer, and the resistivity change range is 17.42-20.5Ω.m after measurement; therefore, the overall resistivity of the target landfill in the research area is considered to be smaller than 12Ω & m due to different water content conditions, but the resistivity of the target landfill in the research area is considered to be lower due to the fact that the water level in the research area is lower and the water content is rich; the resistivity of the clay layer of the main covering stratum is wholly larger than 17 omega-m and reaches about 20 omega-m on average,
from this, it can be determined that the formulaA value greater than 0.2, i.e., a significant electrical difference between the known landfill and the surrounding formation is considered, which can be distinguished using resistivity detection methods;
2. determining a measuring area, wherein the measuring lines are T1-T5, calculating the parameters and the model according to the embodiments 1-3 to obtain the following results,
the elements of each measuring line are shown in table 1, and the positions of the measuring lines are shown in fig. 1;
table 1 survey area comprehensive physical survey line element list
If fig. 3-4 are three device resistivity sections of the line T1, fig. 3-4 show a measured apparent resistivity section of the temperature receiving device, fig. b shows an inverted resistivity section of the temperature receiving device, fig. c shows an measured apparent resistivity section of the temperature applying device, fig. d shows an inverted resistivity section of the temperature applying device, fig. e shows an measured apparent resistivity section of the dipole device, and fig. f shows an inverted resistivity section of the dipole device;
as can be seen from the comparison of the graphs a, c and e in fig. 3-4, the overall morphology of the section with apparent resistivity measured by the three devices is similar, the contour line of resistivity is smoother, no obvious abrupt change point exists, and the quality of the original data is reflected to be higher; the vertical change characteristics of the contour line of the measured apparent resistivity of the three devices are basically consistent, the shallow resistivity is higher and is reflected by hardened ground and miscellaneous fill, and the deep resistivity is relatively lower and is reflected by the fourth-series loose-bed silty clay, silty sand and local landfill. The characteristic of the change of the transverse resistivity contour line has certain difference, specifically, the transverse change of the resistivity contour line of the dipole device is most obvious, the temperature application device is inferior, and the temperature application device is extremely weak, which reflects the difference of the transverse resolutions of the three devices, namely, the transverse resolution of the dipole device is highest, the temperature application device is inferior, and the temperature application device is relatively low. The three devices can be used for inverting the overall shape of the resistivity section to better reflect the burial condition of the landfill as can be seen from the inversion resistivity sections of the figures b, d and f. In contrast, the temperature receiving device and the temperature applying device are basically consistent in position of the landfill on the inverted resistivity section, the dipole device is different from the temperature receiving device and the temperature applying device to a certain extent, and the reason for consideration is probably that the dipole device is an asymmetric device, and the abnormal shape is often complicated;
the comprehensive analysis can be achieved, and the temperature applying device has the characteristics of relatively high transverse resolution, simple abnormal shape and large data collection amount, and can be used as a first-choice device in the landfill investigation; the temperature and humidity device has the characteristics of strong anti-interference capability and simple abnormal shape; the dipole device has the highest transverse resolution and has advantages in layering capability, so that the temperature application device is used as a main device, the temperature application device and the dipole device are used as auxiliary devices, the detection results are mutually complemented and mutually verified, and the advantages of different devices are fully exerted;
the high-density resistivity method in the investigation region has the length of the survey line below 75m under the limitation of the field condition, and in order to obtain the maximum exploration depth under the condition of limited length of the survey line, the maximum isolation coefficient 30 allowed by the device is acquired by the determining device, and the power supply voltage is 450 and is better;
5-9 are inversion resistivity and geological interpretation sections of the high density resistivity method of the measuring lines T1-T5, wherein a diagram a is the inversion resistivity section, b is the geological interpretation section, the abscissa is the pile number, and the ordinate is the depth; deducing the position of a suspected landfill in a low-resistance abnormal area encircled in the graph;
as can be seen from fig. 4 to 8, the T1 coil defines the suspected landfill location 2, numbered (1) and (2); the T2 coil is used for determining the position 1 of the suspected landfill object, and the number (2) is given; the T3 coil is used for determining the position 1 of the suspected landfill object, and the number (2) is given; no obvious signs of suspected landfill were found on line T4; the T5 coil is positioned at the suspected landfill position 2, and is numbered (3) and (4); the low-resistance abnormal region is found in the T1-T5 lines except for T4;
3. three-dimensional spatial inference interpretation
As shown in fig. 10, which is a three-dimensional map of inversion resistivity of the high-density resistivity method of the measurement region, as can be seen from fig. 10, there is a low-resistance anomaly region 4 in the measurement region, and numbers (1) to (4), wherein anomalies (2) and (3) are located in the middle of the measurement region, and the anomaly range and intensity are relatively large; (1) The abnormal numbers (4) and (4) are respectively positioned in the south and north of the measuring area, and the abnormal range and the abnormal strength are relatively small; deducing that the low-resistance anomalies are the suspected landfill positions, and taking the anomaly areas of (2) and (3) as main landfill positions;
based on the analysis result and the obtained coordinate data, the plane area and square quantity are estimated as follows:
landfill information list
The estimated planar area at the point 4 of the total delineating suspected landfill in the investigation region is 372.2m 2 The method comprises the steps of carrying out a first treatment on the surface of the Estimated square 1914.2m in effective detection depth range 3 。
Claims (9)
1. A method for determining a location of a solid waste landfill based on a high density resistivity process, comprising the steps of:
s1, resistivity pre-judgment
Adopting two modes of ground sampling and drilling sampling, testing resistivity values of known landfill and stratum around the landfill in soil of a region to be confirmed, judging whether the known landfill and the stratum have obvious electrical difference according to the difference range of the resistivity values of the known landfill and the stratum around the known landfill, if so, further exploring the landfill by using a high-density resistivity method, otherwise, not researching the landfill;
s2, arranging the measuring lines:
arranging a plurality of high-density resistivity measuring lines in a solid waste landfill area to be confirmed;
s3, arranging electrodes:
selecting a plurality of measuring points with equal intervals on each high-density resistivity measuring line, arranging electrodes on the measuring points, and then connecting the electrodes to a program-controlled multi-channel electrode conversion switch by using a multi-core cable;
s4, collecting data:
sending the measurement signal generated by the electrode into a micro-electromechanical measuring instrument through an electrode change-over switch, and storing the measurement result into a random access memory; the data of the random access memory is firstly transmitted to a computer through a network, and the WDAFC.EXE software of the computer is utilized to convert the original data format into a format which can be read by Res2dinv software;
s5, processing data:
pre-processing the raw data using Res2dinv processing software, the pre-processing comprising: splicing and removing mutation points; then, a forward model is established, forward calculation is carried out by adopting Res2dmod, and potential values and resistivity profiles at grid nodes are obtained through the forward calculation; the grid is defined by Res2dmod software, and then a least square method is adopted to perform two-dimensional inversion calculation to obtain an inversion result of the forward model;
in the inversion process, all models select 3-5 times of inversion iteration results; obtaining an RMS error value in the multi-iteration inversion result; when the obtained RMS error values are smaller than 0.2, calculating variance values of a plurality of groups of RMS error values, and repeating the inversion process when the RMS error values are larger than or equal to 0.2;
judging whether the variance value is more than 15%, if so, judging that the inversion result is reasonable, and if so, judging that the inversion result is unreasonable; when the inversion result is reasonable, the inversion result is derived, a multilayer inversion resistivity section diagram is drawn by using Res2dinvx64 software, and meanwhile, coordinate conversion is carried out; repeating the inversion process when the inversion result is unreasonable;
s6, confirming area
Merging the inversion resistivity section data of each layer in the inversion result to obtain full-region inversion resistivity three-dimensional space data, and carrying out three-dimensional visualization on the full-region inversion resistivity three-dimensional space data according to RockWorks processing software to obtain an inversion resistivity plane graph and a three-dimensional stereogram; and searching a common area with the largest resistance difference of the inversion resistivity section view, the inversion resistivity plane view and the three-dimensional stereogram to obtain the position coordinates of the abnormal resistance area, and determining the position coordinates as a solid waste landfill area.
2. The method for determining a landfill location of solid waste based on high density resistivity method according to claim 1, wherein the determining in the step S1 whether the difference between the resistivity values of the known landfill and the surrounding stratum has a significant electrical difference is:
according to the formulaWhether the calculated value is larger than 0.2 or not is judged to have obvious electrical property difference, otherwise, the calculated value does not have obvious electrical property difference; where R0 is the resistivity value of the surrounding formation and R1 is the resistivity value of the known landfill.
3. The method for determining a landfill site of solid waste based on a high density resistivity method as claimed in claim 1, wherein the step S4 is to collect measurement signals by using a Chongqing Pentium WGMD-9 high density electrical method system, and the collection parameters are as follows: the number of the electrodes is 60-180, the electrode distance is 1m, the maximum isolation coefficient is 19-30, and the power supply voltage is 450V.
4. A method for determining a location of a solid waste landfill based on a high density resistivity method according to claim 3, wherein the measured voltage is 5 to 10mV.
5. The method for determining a landfill position of solid waste based on a high-density resistivity method according to claim 1, wherein in the step S5, after the calculation is completed, the quality of a section curve or a resistivity section is checked, whether there are dead spots or not is judged, and when there are few dead spots or dead spots, the dead spots are 5 or less and when the area of the area is 1 or less, a single-point measurement function retest is started; and retesting the whole section when the number of the bad points is more and the number of the bad points is more than 5 and the area of the area is more than 1%.
6. A method for determining a landfill location of solid waste based on a high density resistivity method according to claim 1, wherein before the electrodes are laid in step S3, holes are punched in the ground and filled with soil, and then brine is poured until the electrode ground resistance is reduced to less than 2kΩ.
7. The method for determining the landfill position of the solid waste based on the high-density resistivity method according to claim 1, wherein before the data acquisition in the step S4, a parameter selection test is performed, data of a temperature receiving device, a temperature applying device and a dipole device are simultaneously acquired by each section, and parameters adopting the device, the device form, the isolation coefficient and the power supply voltage are determined according to inversion results of the three devices and position coordinates of abnormal areas of resistance values.
8. The method for determining a solid waste landfill site based on a high density resistivity method as claimed in claim 1, wherein the forward model in the step S5 is: the number of the electrodes is 81, the electrode distance is 1m, the low-resistance different area is 0-3 omega, and the surrounding stratum is 17-20 omega; the number of layers of resistivity data acquisition was 20.
9. The method for determining a landfill location of solid waste based on a high density resistivity method according to claim 1, wherein after the area is confirmed in step S6, a planar area of the solid waste landfill area is estimated, and a square quantity thereof is estimated in a detection depth range.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310576239.2A CN116774295A (en) | 2023-05-22 | 2023-05-22 | Method for determining solid waste landfill position based on high-density resistivity method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310576239.2A CN116774295A (en) | 2023-05-22 | 2023-05-22 | Method for determining solid waste landfill position based on high-density resistivity method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116774295A true CN116774295A (en) | 2023-09-19 |
Family
ID=87990503
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310576239.2A Pending CN116774295A (en) | 2023-05-22 | 2023-05-22 | Method for determining solid waste landfill position based on high-density resistivity method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116774295A (en) |
-
2023
- 2023-05-22 CN CN202310576239.2A patent/CN116774295A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhou et al. | Three‐dimensional spatial and temporal monitoring of soil water content using electrical resistivity tomography | |
| CN107829453B (en) | Method and device for detecting leakage of vertical plastic-spreading impermeable curtain | |
| CN110702587B (en) | Earth-rock dam leakage diagnosis method based on Wennan joint inversion | |
| CN106405250B (en) | High-density ground resistivity measuring system and method suitable for complex terrain condition | |
| Ungureanu et al. | Use of electric resistivity tomography (ERT) for detecting underground voids on highly anthropized urban construction sites | |
| US4296379A (en) | Ground prospecting method utilizing electrical resistivity measurements for measuring the resistivity of unit blocks of the ground | |
| US5914603A (en) | Electrical resistance tomography using steel cased boreholes as electrodes | |
| CN110988999A (en) | Detection method and system for analyzing pile foundation based on cross-hole resistivity method CT inversion imaging | |
| Acworth | Investigation of dryland salinity using the electrical image method | |
| CN206378448U (en) | Polluted Soil detection arrangement of measuring-line structure based on comprehensive geophysical prospecting methods | |
| Jouen et al. | Evaluation and localization of an artificial drainage network by 3D time-lapse electrical resistivity tomography | |
| CN110673216B (en) | Single-hole resistivity karst cave detection method | |
| CN111721831A (en) | Three-dimensional tomography dam hidden leakage channel scanning method based on electric stimulation | |
| CN112485829B (en) | Method for calculating dynamic water content by resistivity method | |
| CN116774295A (en) | Method for determining solid waste landfill position based on high-density resistivity method | |
| CN116165718A (en) | Method and device for measuring seepage safety of earth and rockfill dam, electronic equipment and medium | |
| CN110967767A (en) | Method for detecting karst cave by cross-hole induced polarization | |
| Nobahar et al. | Moisture variation monitoring of failed and not-failed highway slope through resistivity imaging in Mississippi | |
| CN115421206A (en) | Measuring method for fast survey of land parcel based on electrical prospecting | |
| CN113406707A (en) | Magnetotelluric multi-scale and multi-time-period detection method | |
| RU2678535C1 (en) | Method of monitoring state of diaphragm from boron-cutting alumina-cement concrete piles in earth dam by electrical tomography | |
| CN113447991A (en) | Method and device for rebuilding underground electrical anomaly | |
| CN111721830A (en) | Three-dimensional tomography-based channel improvement dam body internal erosion detection, monitoring and early warning method | |
| CN116819647B (en) | Hydrologic geophysical data fusion method based on cross gradient structure constraint | |
| CN113447990B (en) | Method and device for observing abnormal electrical property of well site |
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 |