CN110308486A - A kind of Polluted area boundary demarcation method - Google Patents

A kind of Polluted area boundary demarcation method Download PDF

Info

Publication number
CN110308486A
CN110308486A CN201910742282.5A CN201910742282A CN110308486A CN 110308486 A CN110308486 A CN 110308486A CN 201910742282 A CN201910742282 A CN 201910742282A CN 110308486 A CN110308486 A CN 110308486A
Authority
CN
China
Prior art keywords
grid block
sensor
grid
rock mass
wave
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.)
Granted
Application number
CN201910742282.5A
Other languages
Chinese (zh)
Other versions
CN110308486B (en
Inventor
董陇军
张义涵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Central South University
Original Assignee
Central South University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Central South University filed Critical Central South University
Priority to CN201910742282.5A priority Critical patent/CN110308486B/en
Publication of CN110308486A publication Critical patent/CN110308486A/en
Application granted granted Critical
Publication of CN110308486B publication Critical patent/CN110308486B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/288—Event detection in seismic signals, e.g. microseismics
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G01V1/305—Travel times
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a kind of Polluted area boundary demarcation method, include the following steps: to obtain various types of rock mass in reference velocity of wave uncontaminated and when having polluted；Grid dividing is carried out to detection zone, grid extends certain depth to vertically apparent bearing and forms grid block；Monitoring grid block rock mass velocity；The measurement of grid block depth is polluted, grid block is divided into L layers by vertically apparent bearing, solves the rock mass velocity that grid block is respectively layered；Pollution range delimited, if certain grid block/reference velocity of wave of the grid block layering measurement velocity of wave relative to same type rock mass when uncontaminated, closer to reference velocity of wave of the same type rock mass when having polluted, then grid block/the grid block is layered as Polluted area, determines Polluted area according to all pollution grid blocks/grid block layering periphery.The present invention is based on rock mass medium wave speed difference is different, it is not necessary that detection zone sampling determination, the pollution range that efficiently can easily complete detection zone delimited.

Description

A kind of Polluted area boundary demarcation method
Technical field
The invention belongs to environment rock mass engineering project technical fields, and in particular to a kind of side delimited for Polluted area boundary Method.
Background technique
The heavy polluted industries such as mining, petrochemical industry, weaving can generate a large amount of pollutants, wherein industrial wastewater in development process The problems such as discharge of waste residue, heavy metal contaminants are stacked causes the pollution of rock mass environment.In recent years, some " dirt at random " enterprises care for And number one, waste is not handled according to the provisions of the relevant regulations issued by the State, direct discharge in violation of regulations, periphery rock mass environment is caused to dislike Change, has seriously threatened rock mass engineering project and the quality of living environment, caused great social safety hidden danger.With whole society's safety collar Consciousness enhancing is protected, the solution of this problem is extremely urgent.
The range for delimiting contaminated site simultaneously obtains more reliable boundary, be carry out the environmental impact assessment of rock mass engineering project stability estimate, field The committed step of the work such as ground reparation and contaminating enterprises' regulation.Geostatistics in Soil Environmental Pollution detection have been widely applied It is predicted in the spatial distribution of ground pollutant.But need to sample area to be tested in actual mechanical process, take back room Interior carry out sample test, often due to being influenced by the Spatial Variability of sampling density and pollutant, the pollution situation of sampling point cannot be complete Practical pollution situation in grid delimited in face reflection, moreover lack effective pollution situation in the environment based on rock mass in the past Analysis means.
Therefore, it is necessary to design a kind of side of quick judgement contaminated sites range that can be suitably used for the environment based on rock mass Method.
Summary of the invention
Technical problem solved by the invention is in view of the deficiencies of the prior art, to provide a kind of completely new Polluted area side Boundary's demarcation method efficiently can easily complete pollution range without to detection zone sampling determination and delimit.
For achieving the above object, the present invention adopts the following technical scheme:
A kind of Polluted area boundary demarcation method, comprising the following steps: firstly, obtaining various types of rock mass uncontaminated Reference velocity of wave when having polluted；Then, grid dividing is carried out to detection zone, each grid is extended one to vertically apparent bearing Depthkeeping degree forms grid block；Each grid block is divided into L layers by vertically apparent bearing；It measures in each grid block and its each classification Rock mass velocity；Finally, carrying out the delimitation of Polluted area boundary；If certain grid block/grid block layering measurement velocity of wave is relative to similar Reference velocity of wave of type rock mass when uncontaminated, closer to reference velocity of wave of the same type rock mass when having polluted, then grid block/the net Lattice block is layered as pollution grid block/grid block layering；One closed region is formed by the periphery of all pollution grid blocks, by this Horizontal boundary of the horizontal boundary in region as Polluted area；By it is all pollution grid blocks layering periphery formed one it is closed Region, using the vertical boundary in the region as the vertical boundary of Polluted area；Thus Polluted area boundary is completed to delimit.
A kind of Polluted area boundary demarcation method, includes the following steps:
Step 1 establishes the velocity of wave table of comparisons, for store various types of rock mass polluted with it is uncontaminated when reference wave Speed；
Step 2 carries out grid dividing to detection zone；
Region division be will test into multiple square nets；Each grid is extended into certain depth to the direction perpendicular to earth's surface, Form grid block；
Step 3 detects rock mass velocity in each grid block respectively；
Step 4, the horizontal boundary for determining pollution grid block and Polluted area；
For each grid block, according to its rock body types, the rock mass of same type is searched in the velocity of wave table of comparisons Reference velocity of wave when polluting and is uncontaminated；If the grid block solves the phase that obtained rock mass velocity is obtained with lookup by step 3 The difference of reference velocity of wave of the same type rock mass when having polluted is less than its reference velocity of wave with same type rock mass when uncontaminated Difference, then assert the grid block for pollution grid block；After having judged all grid blocks, by the outer of all pollution grid blocks It encloses to form a closed region, using the horizontal boundary in the region as the horizontal boundary of Polluted area；
Step 5, each pollution grid block are layered wave velocity measurement；
To each pollution grid block, it is divided into L layering by the direction perpendicular to earth's surface, is detected in its each layering Rock mass velocity, so that it is determined that grid block pollution depth；
Step 6 determines the pollution layering in each pollution grid block and the vertical boundary of Polluted area.
For each layering in pollution grid block, according to its rock body types, same type is searched in the velocity of wave table of comparisons Rock mass polluted with it is uncontaminated when reference velocity of wave, if pollution grid block layering solves obtained rock mass by step 5 Velocity of wave and the difference for searching reference velocity of wave of the obtained same type rock mass when having polluted are less than it and exist with same type rock mass The difference of reference velocity of wave when uncontaminated then assert that this is layered as pollution layering；All layerings of all pollution grid blocks are sentenced After having broken, a closed region is formed by the periphery of all pollution layerings, using the vertical boundary in the region as contaminated area The vertical boundary in domain；
The vertical edges for the Polluted area that step 7, the horizontal boundary of the Polluted area determined according to step 4 and step 6 determine The boundary of boundary's delimitation Polluted area.
Further, in the step 3, to each grid block, the rock mass velocity in it is determined by the following method respectively；
Step 3.1, different location arranges that N number of sensor, N are the integer more than or equal to 4 on the grid block；Each sensing Known to the position coordinates of device；
Step 3.2 arranges a man-made explosion on the grid block, records each sensor on the grid block and receives this The real time for the sound wave that man-made explosion generates；
Step 3.3, the real time recorded using step 3.2, are obtained between each sensor arranged on the grid block Actual measurement arrival time difference；It is listed according to most short propagation path, velocity of wave and the functional relation of time of sound wave and to be arranged on the grid block The function expression of theoretical arrival time difference between each sensor；Based on the reality between each sensor arranged on the grid block It measures the time difference and theoretical arrival time difference and solves the rock mass velocity of the grid block using least square method.
Further, in the step 3, in the step 3.3, the objective function of the rock mass velocity in grid block is solved Are as follows:
Δtij=ti-tj
Wherein, Q indicates target function value, x0, y0, z0, c is variable to be solved, (x0, y0, z0) it is cloth on the grid block The position coordinates for the man-made explosion set, c are the rock mass velocity in the grid block；For i-th of the biography arranged on the grid block The theoretical arrival time difference of sensor and j-th of sensor, Δ tijFor i-th of the sensor and j-th of sensor arranged on the grid block Actual measurement arrival time difference；WithThe theory of i-th of the sensor and j-th of sensor respectively arranged on the grid block is then；li And ljThe artificial shake arranged on i-th of the sensor and j-th of sensor and the grid block respectively arranged on the grid block I-th sensor being arranged in the man-made explosion that the shortest distance in source, i.e. sound wave are arranged from the grid block to the grid block and The most short propagation path of j-th of sensor；(xi, yi, zi) and (xj, yj, zj) it is respectively i-th of the sensing arranged on the grid block The position coordinates of device and j-th of sensor；t0Man-made explosion to arrange on the grid block generates the time of sound wave；tiAnd tjPoint I-th of the sensor and j-th of sensor that Wei not arrange on the grid block receive the man-made explosion arranged on the grid block and produce The real time of raw sound wave；I, j=1,2 ..., N.
Further, in the step 5, the rock mass wave in a pollution each layering of grid block is determined by following steps Speed；
Step 5.1, different location arranges that M sensor, M are the integer more than or equal to 4 on the grid block；Each sensing Known to the position coordinates of device；
Step 5.2 is respectively arranged a man-made explosion in each layering of the grid block, successively excites each artificial Focus, and record each sensor arranged on the grid block and receive the man-made explosion generation that the grid block is respectively layered interior arrangement Sound wave real time；
Step 5.3, the real time recorded using step 5.2, are obtained between each sensor arranged on the grid block Actual measurement arrival time difference；It is listed according to most short propagation path, velocity of wave and the functional relation of time of sound wave and to be arranged on the grid block The function expression of theoretical arrival time difference between each sensor；Based on the reality between each sensor arranged on the grid block It measures the time difference and theoretical arrival time difference and solves the rock mass velocity of each layering in the grid block using least square method.
Further, in the step 5, if the man-made explosion for being layered interior arrangement for k-th of a grid block is located at point DkPlace, is asked Solve the objective function of the rock mass velocity of each layering in the grid block are as follows:
Δtijk=tik-tjk
Wherein, P indicates target function value, c1, c2..., cLFor variable to be solved, i.e., the grid block the 1~L layering Interior rock mass velocity；Indicate that i-th of the sensor arranged on the grid block and j-th of sensor receive point DkThe people at place The theoretical time difference for the sound wave that work focus generates,WithI-th of the sensor and jth respectively arranged on the grid block A sensor receives point DkThe theoretical time for the sound wave that the man-made explosion at place generates；ΔtijkFor arranged on the grid block i-th A sensor and j-th of sensor receive point DkThe real time for the sound wave that the man-made explosion at place generates is poor, tikAnd tjkRespectively Point D is received for i-th of the sensor arranged on the grid block and j-th of sensorkThe sound wave that the man-made explosion at place generates Real time.
Further, to each grid block, by its L layering, from earth's surface, number consecutively is l=1,2 ..., L, l down The thickness of a layering is denoted as hl, hlIt is known；The left upper apex for taking grid block is coordinate origin, is down z-axis pros perpendicular to earth's surface To；
If i-th of the sensor arranged on the grid block is located at the grid block top edge, point D is receivedkThe people at place The theoretical time for the sound wave that work focus generates are as follows:
If i-th of the sensor arranged on the grid block is located at the grid block lower edge, point D is receivedkThe people at place The theoretical time for the sound wave that work focus generates are as follows:
Wherein,
Δhk1Indicate the man-made explosion and the longitudinal direction of the grid block kth -1 layering bottom surface in k-th of grid block layering Distance；Δhk2Indicate the man-made explosion and the fore-and-aft distance of kth+1 layering top surface in k-th of grid block layering；cosθik For i-th of the sensor arranged on the grid block and k-th of the grid block shortest path being layered between the man-made explosion of interior arrangement The angle in diameter direction and grid block longitudinal direction, 0 ° < θik<90°；(Xk, Yk, Zk) it is point DkLocate the position coordinates of man-made explosion；likFor The shortest distance of i-th of the sensor arranged on the grid block and the man-made explosion in k-th of grid block layering, i.e., k-th Most short propagation path of the sound wave that man-made explosion in layering generates to i-th of sensor；(xi, yi, zi) it is i-th of sensing The position coordinates of device.
Further, in the step 3.1,5 sensors are arranged on each grid block；The arrangement side of 5 sensors Formula and position are as follows:
Perpendicular to ground apparent bearing along the two neighboring rib drilling of grid block, two sensors are arranged in grid block bottom along drilling On two vertex of face side, excess-three sensor is located at the grid top surface (i.e. the grid block rock mass surface), two of them On two vertex of the grid block top surface other side, another sensor is located in the rib of grid block top surface side sensor arrangement Between on position.
The step 5.1 keeps the arrangement number and position that sense in step 3.1 constant.
The principle of the invention are as follows:
After trade waste is in contact with rock, rock with water gradually generates chemical reaction, physical reactions change rock Structure contact inside stone, be reflected in Rock Velocity (spread speed of the wave in rock mass) have in contaminated front and back it is apparent poor It is different, such as a large number of studies show that after being corroded by heavy metal rock mass internal structure and mineral constituent can change, wherein granulated State and grain spacing, the change of the microstructures such as hole, crack size cause rock mass porosity to increase, the rock mass of lithology of the same race Velocity of wave declines.The present invention is on this basis, micro- by the arrival time difference without measured in advance velocity of wave by microseismic detection equipment Shake seismic source location method, the rock mass velocity in each grid of exact inversion detection zone, by comparing detection zone rock mass velocity with Same type rock mass polluted with it is not comtaminated when reference velocity of wave difference, determine rock mass pollution range.The present invention can be In the case where to detection zone sampling determination, scene is measured rock mass velocity in detection zone, compares detection zone Velocity of wave difference determines rock mass pollution range, and operation sequence is simple, efficiently can easily complete pollution range delimitation.
The utility model has the advantages that
The present invention determines rock mass pollution range by Determination region velocity of wave difference and delimit boundary, and detection process is just Victory more can subtly mark off the pollution range of area to be tested, provide efficiently for rock mass engineering project stability entry evaluation Ground means.
Detailed description of the invention
Fig. 1 is the shortest path schematic diagram of man-made explosion and sensor in grid block layering
Fig. 2 is sensor and man-made explosion location arrangements schematic diagram
Fig. 3 is the schematic diagram that pollution range boundary delimited
Specific embodiment
The present invention is further described in detail below in conjunction with attached drawing.
The invention discloses a kind of Polluted area boundary demarcation methods, include the following steps:
Step 1 establishes the velocity of wave table of comparisons, for store various types of rock mass polluted with it is uncontaminated when reference wave Speed；
The rock body types of detection zone can be determined by consulting related geologic information；
It is sampled by having been polluted to detection zone periphery with rock mass uncontaminated and identical with detection zone type, and Measurement sampling rock mass velocity, so that it is determined that all types of rock mass polluted with it is uncontaminated when velocity of wave, as reference velocity of wave, storage In the velocity of wave table of comparisons；
Step 2 carries out grid dividing to detection zone；
Region division be will test into multiple square nets；Each grid is extended into certain depth to the direction perpendicular to earth's surface, Form grid block；
Step 3 detects rock mass velocity in each grid block respectively；
Specifically, to each grid block, the rock mass velocity in it is determined by the following method respectively；
Step 3.1, (grid block edge) different location arranges that N number of sensor, N are more than or equal to 4 on the grid block Integer；Known to the position coordinates of each sensor；As shown in Fig. 2, arranging 5 sensings on each grid block in the present embodiment Device；The arrangement of 5 sensors and position are as follows:
Perpendicular to ground apparent bearing along the two neighboring rib drilling of grid block, two sensors are arranged in grid block bottom along drilling On two vertex of face side, excess-three sensor is located at the grid block top surface (i.e. the grid block rock mass surface), wherein two For a sensor arrangement on two vertex of the grid block top surface other side, another sensor is located at the rib of grid block top surface side On middle position.Sensor arrangement makes the line between sensor run through entire grid as far as possible at unit grid body edge Body, to form cube monitoring network；
Step 3.2, (grid block surface or internal any position) one man-made explosion of arrangement on the grid block, record should Each sensor receives the real time of the sound wave of man-made explosion generation on grid block；Man-made explosion need to be by grid block Top surface drilling is arranged in inside grid block.
Step 3.3, the real time recorded using step 3.2 obtain actual measurement arrival time difference；Road is broadcast according to the most short pass of sound wave Diameter, velocity of wave and the functional relation of time list the function expression of theoretical arrival time difference；Based on actual measurement arrival time difference and theoretical arrival time difference, Using least square method, the rock mass velocity (i.e. spread speed of the sound wave in the grid block) of the grid block is solved；Solve grid The objective function of rock mass velocity in block are as follows:
Δtij=ti-tj
Wherein, Q indicates target function value, x0, y0, z0, c is variable to be solved, (x0, y0, z0) it is cloth on the grid block The position coordinates for the man-made explosion set, c are the rock mass velocity in the grid block；For i-th of the biography arranged on the grid block The theoretical arrival time difference of sensor and j-th of sensor, Δ tijFor i-th of the sensor and j-th of sensor arranged on the grid block Actual measurement arrival time difference；WithThe theory of i-th of the sensor and j-th of sensor respectively arranged on the grid block is then；li And ljThe artificial shake arranged on i-th of the sensor and j-th of sensor and the grid block respectively arranged on the grid block I-th sensor being arranged in the man-made explosion that the shortest distance in source, i.e. sound wave are arranged from the grid block to the grid block and The most short propagation path of j-th of sensor；(xi, yi, zi) and (xj, yj, zj) it is respectively i-th of the sensing arranged on the grid block The position coordinates of device and j-th of sensor；t0Man-made explosion to arrange on the grid block generates the time of sound wave；tiAnd tjPoint I-th of the sensor and j-th of sensor that Wei not arrange on the grid block receive the man-made explosion arranged on the grid block and produce The real time of raw sound wave；I, j=1,2 ..., N.
For (xi, yi, zi, xj, yj, zj), if theoretical valueWith measured value Δ tijDeviation it is smaller, then it is assumed that straight line and The degree of fitting of all testing sites is better.Then (x0, y0, z0, c) and Q (x should be made0, y0, z0, c) and reach minimum.The above objective function Differential nonlinear least square fitting is as sought, the rock mass velocity that can be obtained in man-made explosion position coordinates and grid block is solved.
Step 4, the horizontal boundary for determining pollution grid block and Polluted area；
For each grid block, according to its rock body types, the rock mass of same type is searched in the velocity of wave table of comparisons Reference velocity of wave when polluting and is uncontaminated；If the grid block solves the phase that obtained rock mass velocity is obtained with lookup by step 3 The difference of reference velocity of wave of the same type rock mass when having polluted is less than its reference velocity of wave with same type rock mass when uncontaminated Difference, then assert the grid block for pollution grid block；After having judged all grid blocks, by the outer of all pollution grid blocks It encloses to form a closed region (region is pollution depth range to be determined), using the horizontal boundary in the region as pollution The horizontal boundary in region；
Step 5, each pollution grid block layering wave velocity measurement (i.e. the pollution depth of measurement pollution grid block)；
To each pollution grid block, it is divided into L layering by the direction perpendicular to earth's surface, is detected in its each layering Rock mass velocity, so that it is determined that grid block pollution depth；Specifically, each point of grid block of pollution is determined by following steps Rock mass velocity in layer；
Step 5.1, (grid block edge) different location arranges that M sensor, M are more than or equal to 4 on the grid block Integer；Known to the position coordinates of each sensor；In this implementation, the step keep on the grid block sensor position and Number constant (M=N), i.e. sensor arrangement in holding step 3.1；
Step 5.2 is respectively arranged a man-made explosion in each layering of the grid block, successively excites each artificial Focus, and record each sensor arranged on the grid block and receive the man-made explosion generation that the grid block is respectively layered interior arrangement Sound wave real time；
Step 5.3, the real time recorded using step 5.2, obtain the actual measurement arrival time difference between each sensor；According to Most short propagation path, velocity of wave and the functional relation of time of sound wave list the function table of the theoretical arrival time difference between each sensor Up to formula；It is solved in the grid block with theoretical arrival time difference using least square method based on the actual measurement arrival time difference between each sensor The rock mass velocity of each layering.Specifically, if the man-made explosion for being layered interior arrangement for k-th of a grid block is located at point DkPlace solves The objective function of the rock mass velocity of each layering in the grid block are as follows:
Δtijk=tik-tjk
Wherein, P indicates target function value, c1, c2..., cLFor variable to be solved, i.e., the grid block the 1~L layering Interior rock mass velocity (i.e. spread speed of the sound wave in the grid block the 1~L layering)；The above objective function is to ask differential non- Linear least square fitting, according to principle of least square method, (c1, c2..., cL) target function value P should be made minimum；Table Show that i-th of the sensor arranged on the grid block and j-th of sensor receive point DkThe sound wave that the man-made explosion at place generates Theoretical time difference,WithI-th of the sensor and j-th of sensor respectively arranged on the grid block receive point Dk The theoretical time for the sound wave that the man-made explosion at place generates；ΔtijkIt is passed for i-th of the sensor arranged on the grid block and j-th Sensor receives point DkThe real time for the sound wave that the man-made explosion at place generates is poor, tikAnd tjkRespectively arranged on the grid block I-th of sensor and j-th of sensor receive point DkThe real time for the sound wave that the man-made explosion at place generates.
Further, in the step 5.1, to each grid block, by its L layering, from earth's surface, number consecutively is l down =1, the thickness of 2 ..., L, first of layering are denoted as hl, hlIt is known；The left upper apex for taking grid block is coordinate origin, perpendicular to ground Table is z-axis positive direction down；
If i-th of the sensor arranged on the grid block is located at the grid block top edge, point D is receivedkThe people at place The theoretical time for the sound wave that work focus generates are as follows:
If i-th of the sensor arranged on the grid block is located at the grid block lower edge, point D is receivedkThe people at place The theoretical time for the sound wave that work focus generates are as follows:
Wherein,
Δhk1Indicate the man-made explosion and the longitudinal direction of the grid block kth -1 layering bottom surface in k-th of grid block layering Distance；Δhk2Indicate the man-made explosion and the fore-and-aft distance of kth+1 layering top surface in k-th of grid block layering；cosθik For i-th of the sensor arranged on the grid block and k-th of the grid block shortest path being layered between the man-made explosion of interior arrangement The angle in diameter direction and grid block longitudinal direction, 0 ° < θik<90°；(Xk, Yk, Zk) it is point DkLocate the position coordinates of man-made explosion；likFor The shortest distance of i-th of the sensor arranged on the grid block and the man-made explosion in k-th of grid block layering, i.e., k-th Most short propagation path of the sound wave that man-made explosion in layering generates to i-th of sensor；(xi, yi, zi) it is i-th of sensing The position coordinates of device.
By taking two sensors of grid block top edge arrangement as an example, for each group observations (h1, h2..., hk,xi, yi, zi, xk, yk, zk, xj, yj, zj), it can determine a regressand value:
It is possible thereby to solve objective function.
Step 6 determines the pollution layering in each pollution grid block and the vertical boundary of Polluted area.
For each layering in pollution grid block, according to its rock body types, same type is searched in the velocity of wave table of comparisons Rock mass polluted with it is uncontaminated when reference velocity of wave, if pollution grid block layering solves obtained rock mass by step 5 Velocity of wave and the difference for searching reference velocity of wave of the obtained same type rock mass when having polluted are less than it and exist with same type rock mass The difference of reference velocity of wave when uncontaminated then assert that this is layered as pollution layering；All layerings of all pollution grid blocks are sentenced After having broken, a closed region is formed by the periphery of all pollution layerings, using the vertical boundary in the region as contaminated area The vertical boundary in domain；
The vertical edges for the Polluted area that step 7, the horizontal boundary of the Polluted area determined according to step 4 and step 6 determine The boundary of boundary's delimitation Polluted area.The Polluted area boundary that the present embodiment delimited is as shown in Figure 3.
The present invention is based on rock mass medium wave speed difference is different, without efficiently can easily complete to detect to detection zone sampling determination The pollution range in region delimited.

Claims (9)

1. a kind of Polluted area boundary demarcation method, which comprises the following steps: firstly, obtaining various types of rocks Body is in reference velocity of wave uncontaminated and when having polluted；Then, grid dividing is carried out to detection zone, by each grid to vertical earth's surface Direction extends certain depth and forms grid block；Each grid block is divided into L layers by vertically apparent bearing；Measure each grid block and its Rock mass velocity in each classification；Finally, carrying out the delimitation of Polluted area boundary；If certain grid block/grid block layering measurement velocity of wave Relative to reference velocity of wave of same type rock mass when uncontaminated, closer to reference velocity of wave of the same type rock mass when having polluted, then Grid block/the grid block is layered as pollution grid block/grid block layering；One closing is formed by the periphery of all pollution grid blocks Region, using the horizontal boundary in the region as the horizontal boundary of Polluted area；By the peripheral shape of all pollution grid block layerings At a closed region, using the vertical boundary in the region as the vertical boundary of Polluted area；Thus Polluted area side is completed Boundary delimit.
2. a kind of Polluted area boundary demarcation method, which comprises the steps of:
Step 1 establishes the velocity of wave table of comparisons, for store various types of rock mass polluted with it is uncontaminated when reference velocity of wave；
Step 2 carries out grid dividing to detection zone；
Region division be will test into multiple square nets；Each grid is extended into certain depth to the direction perpendicular to earth's surface, is formed Grid block；
Step 3 detects rock mass velocity in each grid block respectively；
Step 4, the horizontal boundary for determining pollution grid block and Polluted area；
For each grid block, according to its rock body types, the rock mass that same type is searched in the velocity of wave table of comparisons is polluting With it is uncontaminated when reference velocity of wave；If the grid block is obtained by the rock mass velocity that step 3 solution obtains with lookup mutually similar The difference of reference velocity of wave of the type rock mass when having polluted is less than the difference of its reference velocity of wave with same type rock mass when uncontaminated Value then assert the grid block for pollution grid block；After having judged all grid blocks, by the peripheral shape of all pollution grid blocks At a closed region, using the horizontal boundary in the region as the horizontal boundary of Polluted area；
Step 5, each pollution grid block are layered wave velocity measurement；
To each pollution grid block, it is divided into L layering by the direction perpendicular to earth's surface, detects the rock in its each layering Bulk wave speed, so that it is determined that grid block pollution depth；
Step 6 determines the pollution layering in each pollution grid block and the vertical boundary of Polluted area；
The rock of same type is searched in the velocity of wave table of comparisons according to its rock body types for each layering in pollution grid block Body polluted with it is uncontaminated when reference velocity of wave, if pollution grid block layering solves obtained rock mass velocity by step 5 It is less than it with the difference for searching reference velocity of wave of the obtained same type rock mass when having polluted with same type rock mass not dirty The difference of reference velocity of wave when dye then assert that this is layered as pollution layering；All layerings of all pollution grid blocks have been judged Later, a closed region is formed by the periphery of all pollution layerings, using the vertical boundary in the region as Polluted area Vertical boundary；
The vertical boundary for the Polluted area that step 7, the horizontal boundary of the Polluted area determined according to step 4 and step 6 determine is drawn Determine the boundary of Polluted area.
3. Polluted area boundary according to claim 2 demarcation method, which is characterized in that in the step 3, to each net Lattice block determines the rock mass velocity in it by the following method respectively；
Step 3.1, different location arranges that N number of sensor, N are the integer more than or equal to 4 on the grid block；Each sensor Known to position coordinates；
Step 3.2 arranges a man-made explosion on the grid block, records each sensor on the grid block and receives this manually The real time for the sound wave that focus generates；
Step 3.3, the real time recorded using step 3.2, obtain the reality between each sensor arranged on the grid block Measure the time difference；Listed according to most short propagation path, velocity of wave and the functional relation of time of sound wave arranged on the grid block it is each The function expression of theoretical arrival time difference between sensor；It is arrived based on the actual measurement between each sensor arranged on the grid block The time difference and theoretical arrival time difference solve the rock mass velocity of the grid block using least square method.
4. Polluted area boundary according to claim 3 demarcation method, which is characterized in that in the step 3, the step In 3.3, the objective function of the rock mass velocity in grid block is solved are as follows:
Δtij=ti-tj
Wherein, Q indicates target function value, x0, y0, z0, c is variable to be solved, (x0, y0, z0) it is to be arranged on the grid block The position coordinates of man-made explosion, c are the rock mass velocity in the grid block；For i-th sensor being arranged on the grid block and The theoretical arrival time difference of j-th of sensor, Δ tijFor the actual measurement of i-th of the sensor and j-th of sensor arranged on the grid block Arrival time difference；WithThe theory of i-th of the sensor and j-th of sensor respectively arranged on the grid block is then；liAnd ljPoint The man-made explosion arranged on i-th of the sensor and j-th of sensor and the grid block that Wei not arrange on the grid block is most I-th of the sensor arranged in the man-made explosion that short distance, i.e. sound wave are arranged from the grid block to the grid block and j-th of biography The most short propagation path of sensor；(xi, yi, zi) and (xj, yj, zj) it is respectively i-th of the sensor arranged on the grid block and jth The position coordinates of a sensor；t0Man-made explosion to arrange on the grid block generates the time of sound wave；tiAnd tjThe respectively net I-th of the sensor and j-th of sensor arranged on lattice block receive the sound wave that the man-made explosion arranged on the grid block generates Real time；I, j=1,2 ..., N.
5. Polluted area boundary according to claim 2 demarcation method, which is characterized in that in the step 5, by following Step determines the rock mass velocity in a pollution each layering of grid block；
Step 5.1, different location arranges that M sensor, M are the integer more than or equal to 4 on the grid block；Each sensor Known to position coordinates；
Step 5.2 is respectively arranged a man-made explosion in each layering of the grid block, successively excites each artificial shake Source, and record each sensor arranged on the grid block receive the grid block be respectively layered in arrange man-made explosion generate The real time of sound wave；
Step 5.3, the real time recorded using step 5.2, obtain the reality between each sensor arranged on the grid block Measure the time difference；Listed according to most short propagation path, velocity of wave and the functional relation of time of sound wave arranged on the grid block it is each The function expression of theoretical arrival time difference between sensor；It is arrived based on the actual measurement between each sensor arranged on the grid block The time difference and theoretical arrival time difference solve the rock mass velocity of each layering in the grid block using least square method.
6. Polluted area boundary according to claim 5 demarcation method, which is characterized in that in the step 5, if a grid K-th of the block man-made explosion for being layered interior arrangement is located at point DkPlace, solves the target of the rock mass velocity of each layering in the grid block Function are as follows:
Δtijk=tik-tjk
Wherein, P indicates target function value, c1, c2..., cLFor variable to be solved, i.e., in the grid block the 1~L layering Rock mass velocity；Indicate that i-th of the sensor arranged on the grid block and j-th of sensor receive point DkThe artificial shake at place The theoretical time difference for the sound wave that source generates,WithI-th of the sensor respectively arranged on the grid block and j-th of biography Sensor receives point DkThe theoretical time for the sound wave that the man-made explosion at place generates；ΔtijkFor i-th of the biography arranged on the grid block Sensor and j-th of sensor receive point DkThe real time for the sound wave that the man-made explosion at place generates is poor, tikAnd tjkRespectively should I-th of the sensor and j-th of sensor arranged on grid block receive point DkThe reality for the sound wave that the man-made explosion at place generates Time.
7. Polluted area boundary according to claim 6 demarcation method, which is characterized in that each grid block, by its L From earth's surface, number consecutively is l=1,2 ..., L down for layering, and the thickness of first of layering is denoted as hl, hlIt is known；Take a left side for grid block Upper vertex is coordinate origin, is down z-axis positive direction perpendicular to earth's surface；
If i-th of the sensor arranged on the grid block is located at the grid block top edge, point D is receivedkThe man-made explosion at place The theoretical time of the sound wave of generation are as follows:
If i-th of the sensor arranged on the grid block is located at the grid block lower edge, point D is receivedkThe man-made explosion at place The theoretical time of the sound wave of generation are as follows:
Wherein,
Δhk1Indicate the man-made explosion and the fore-and-aft distance of the grid block kth -1 layering bottom surface in k-th of grid block layering； Δhk2Indicate the man-made explosion and the fore-and-aft distance of kth+1 layering top surface in k-th of grid block layering；cosθikFor the net I-th of the sensor arranged on lattice block and k-th of the grid block shortest path direction being layered between the man-made explosion of interior arrangement With the angle of grid block longitudinal direction, 0 ° < θik<90°；(Xk, Yk, Zk) it is point DkLocate the position coordinates of man-made explosion；likFor the grid The shortest distance of i-th sensor arranged on block and the man-made explosion in k-th of grid block layering, i.e., in k-th layering Man-made explosion generate sound wave to i-th of sensor most short propagation path；(xi, yi, zi) be i-th of sensor position Set coordinate.
8. Polluted area boundary according to claim 6 demarcation method, which is characterized in that in the step 3.1, each 5 sensors are arranged on grid block；The arrangement of 5 sensors and position are as follows:
Perpendicular to ground apparent bearing along the two neighboring rib drilling of grid block, two sensors are arranged in grid block bottom surface one along drilling On two vertex of side, excess-three sensor is located at the grid block top surface, and two of them sensor arrangement is in grid block top surface On two vertex of the other side, another sensor is located on the rib middle position of grid block top surface side.
9. Polluted area boundary according to claim 6 demarcation method, which is characterized in that in the step 5.1, each 5 sensors are arranged on grid block；The arrangement of 5 sensors and position are as follows:
Perpendicular to ground apparent bearing along the two neighboring rib drilling of grid block, two sensors are arranged in grid block bottom surface one along drilling On two vertex of side, excess-three sensor is located at the grid block top surface, and two of them sensor arrangement is in grid block top surface On two vertex of the other side, another sensor is located on the rib middle position of grid block top surface side.
CN201910742282.5A 2019-08-13 2019-08-13 Method for defining boundary of polluted area Active CN110308486B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910742282.5A CN110308486B (en) 2019-08-13 2019-08-13 Method for defining boundary of polluted area

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910742282.5A CN110308486B (en) 2019-08-13 2019-08-13 Method for defining boundary of polluted area

Publications (2)

Publication Number Publication Date
CN110308486A true CN110308486A (en) 2019-10-08
CN110308486B CN110308486B (en) 2020-11-10

Family

ID=68083300

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910742282.5A Active CN110308486B (en) 2019-08-13 2019-08-13 Method for defining boundary of polluted area

Country Status (1)

Country Link
CN (1) CN110308486B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111781641A (en) * 2020-08-10 2020-10-16 中南大学 Passive seismic source positioning method, system and terminal based on active seismic source correction and readable storage medium
CN111855961A (en) * 2020-07-24 2020-10-30 中南大学 Rock mass drilling quality detection method, drilling machine, server and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105004372A (en) * 2015-07-09 2015-10-28 北京市地质工程设计研究院 Geological environment monitoring method
CN106094021A (en) * 2016-06-01 2016-11-09 北京科技大学 A kind of microseism focus method for rapidly positioning based on arrival time difference data base
CN107884822A (en) * 2017-11-13 2018-04-06 北京矿冶研究总院 Method for improving positioning precision of mining micro-seismic source

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105004372A (en) * 2015-07-09 2015-10-28 北京市地质工程设计研究院 Geological environment monitoring method
CN106094021A (en) * 2016-06-01 2016-11-09 北京科技大学 A kind of microseism focus method for rapidly positioning based on arrival time difference data base
CN107884822A (en) * 2017-11-13 2018-04-06 北京矿冶研究总院 Method for improving positioning precision of mining micro-seismic source

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LONGJUN DONG ET AL.: "Theoretical and experimental studies of localization methodology for AE and microseismic sources without pre-measured wave velocity in mines", 《IEEE ACCESS》 *
徐磊: "基于原位电阻率的膨胀土工程性能评价方法研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111855961A (en) * 2020-07-24 2020-10-30 中南大学 Rock mass drilling quality detection method, drilling machine, server and storage medium
CN111781641A (en) * 2020-08-10 2020-10-16 中南大学 Passive seismic source positioning method, system and terminal based on active seismic source correction and readable storage medium

Also Published As

Publication number Publication date
CN110308486B (en) 2020-11-10

Similar Documents

Publication Publication Date Title
Hess et al. Large‐scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts: 3. Hydraulic conductivity variability and calculated macrodispersivities
Sotoudeheian et al. Impact of Middle Eastern dust sources on PM10 in Iran: Highlighting the impact of Tigris‐Euphrates basin sources and Lake Urmia desiccation
Johnston et al. Quantitative methods for studying landscape boundaries
Tavakoli Carbonate reservoir heterogeneity: overcoming the challenges
Bauer et al. Modeling, parameterization and evaluation of monitoring methods for CO 2 storage in deep saline formations: the CO 2-MoPa project
Miller et al. Three‐dimensional mapping of Mt. Ruapehu volcano, New Zealand, from aeromagnetic data inversion and hyperspectral imaging
CN110308486A (en) A kind of Polluted area boundary demarcation method
Robinson et al. Multiscale characterisation of chimneys/pipes: Fluid escape structures within sedimentary basins
De Vita et al. Detailed rock failure susceptibility mapping in steep rocky coasts by means of non-contact geostructural surveys: the case study of the Tigullio Gulf (Eastern Liguria, Northern Italy)
RU2423306C1 (en) Method to assess impact of geodynamic factors at safety of underground gas storage operation in porous bed
Puglia et al. Extensive characterization of Italian accelerometric stations from single-station ambient-vibration measurements
CN114626583A (en) Prediction and evaluation method for rebound potential of groundwater in heterogeneous polluted site after remediation
Pierotti et al. Anomalous CO2 content in the Gallicano thermo-mineral spring (Serchio Valley, Italy) before the 21 June 2013, Alpi Apuane earthquake (M= 5.2)
Asnawi et al. Soil classification in a seismically active environment based on join analysis of seismic parameters
Benson et al. Spatial sampling considerations and their applications to characterizing fractured rock and karst systems
Lewis Outcrop-derived quantitative models of permeability heterogeneity for genetically different sand bodies
CN108051482A (en) A kind of monitoring and pre-alarming method of seawater invasion underground water
Skinner High-resolution hydrogeological characterization of a fractured dolostone municipal supply aquifer to create a refined 3-D conceptual site model with hydrogeologic units
Ward et al. Environmental monitoring: phase 4 final report (April 2018-March 2019)
Grimaz et al. Advancements from a posteriori studies on the damage to buildings caused by the 1976 Friuli earthquake (north-eastern Italy)
Wessel Automatic shear wave splitting measurements at Mt. Ruapehu volcano, New Zealand
Bruines et al. Development and Application of the GeoDFN and HydroDFN at the Mizunami Underground Research Laboratory
Cohn et al. Remote observations of aeolian saltation
Fabriol et al. Results of investigations to design a monitoring program for a CO2 storage project in the Paris Basin (France)
CN108896397A (en) Roof greening charge of surety evaluation method based on On Microseismic Monitoring Technique

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