CN112799054B - Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar - Google Patents

Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar Download PDF

Info

Publication number
CN112799054B
CN112799054B CN202011566874.5A CN202011566874A CN112799054B CN 112799054 B CN112799054 B CN 112799054B CN 202011566874 A CN202011566874 A CN 202011566874A CN 112799054 B CN112799054 B CN 112799054B
Authority
CN
China
Prior art keywords
crack
underground
film
ground penetrating
penetrating radar
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011566874.5A
Other languages
Chinese (zh)
Other versions
CN112799054A (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.)
Shaanxi Shanbei Mining Hanjiawan Coal Co ltd
China University of Mining and Technology Beijing CUMTB
Original Assignee
Shaanxi Shanbei Mining Hanjiawan Coal Co ltd
China University of Mining and Technology Beijing CUMTB
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 Shaanxi Shanbei Mining Hanjiawan Coal Co ltd, China University of Mining and Technology Beijing CUMTB filed Critical Shaanxi Shanbei Mining Hanjiawan Coal Co ltd
Priority to CN202011566874.5A priority Critical patent/CN112799054B/en
Publication of CN112799054A publication Critical patent/CN112799054A/en
Application granted granted Critical
Publication of CN112799054B publication Critical patent/CN112799054B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/885Radar or analogous systems specially adapted for specific applications for ground probing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures

Abstract

The invention discloses a method for acquiring a multi-period underground three-dimensional form of a dynamic crack based on a ground penetrating radar, which comprises the following steps: paving a film on the surface of the area where the crack is located, and ensuring that the surface of the crack is completely covered by the film; pouring liquid on the film; aligning a ground penetrating radar to the crack at the near-surface position of the crack, and performing interval scanning along the vertical direction of the row of cracks as a path to obtain electromagnetic wave data of a plurality of sections of the underground form of the crack at the present time; drawing the liquid out for storage and removing the film; carrying out data processing on all acquired electromagnetic wave data in the current period, and drawing and obtaining the outlines of cracks at each section in the current period by a visual interpretation method; constructing a three-dimensional model of the current-stage underground form of the crack by utilizing the contour of the crack to obtain the current-stage underground three-dimensional form of the crack and related data thereof; and continuously repeating the steps to obtain the multi-stage underground three-dimensional form and related data thereof in the dynamic change process of crack development. The method can obtain the multi-stage underground form of the dynamic fracture without damage.

Description

Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar
Technical Field
The invention relates to the field of land reclamation of coal mining subsidence lands, in particular to a method for acquiring a multi-period underground three-dimensional form of a dynamic crack based on a ground penetrating radar.
Background
In western regions of China, particularly in western aeolian sand regions and loess regions, ground cracks are one of geological environment problems caused by coal exploitation, and cause the problems of building deformation, underground pipeline damage, cultivated land damage, accelerated soil moisture evaporation, vegetation damage, water and soil loss and the like, so that great difficulty is brought to mining area management workers, and the method is also an important link for mining area land reclamation. In order to study the effect of the ground fissures on the ecological environment, obtaining the underground morphology of the ground fissures is essential for assessing the risk thereof and for studying the development law of the ground fissures. At present, a learner obtains the underground form of the crack through a ground penetrating radar technology by injecting gypsum slurry into the crack, after gypsum is solidified, but the method can damage the ground crack structure, influence the subsequent development process of the ground crack, and insufficient filling of the crack can be caused by uneven gypsum slurry filling, so that the method cannot be used for researching the underground form of the dynamic crack. Meanwhile, a learner digs a working section on one side of the ground fracture, and then obtains the underground form of the fracture through a ground penetrating radar technology transverse scanner, and the method is nondestructive to the fracture, but the required engineering amount is extremely large and time and labor are consumed when the depth of the fracture is large, and when the width of the fracture is small, the result of the ground penetrating radar can generate larger error, and meanwhile, the surrounding ground surface environment can be damaged greatly.
Disclosure of Invention
In order to solve the problems, the invention provides a method for acquiring the multi-period underground three-dimensional form of the dynamic crack based on the ground penetrating radar, which has the characteristics of simple flow, easiness in realization, no influence on the surface environment, extremely small required engineering amount, time saving and labor saving, no influence on the form of the crack and the subsequent development process, capability of acquiring the multi-period underground three-dimensional form in the development process of the dynamic crack in a lossless manner, and great significance for researching the development process and the rule of the underground form of the dynamic crack along with the change of time.
In order to achieve the above purpose, the invention adopts the following technical scheme:
the invention provides a method for acquiring a multi-period underground three-dimensional form of a dynamic crack based on a ground penetrating radar, which is characterized by comprising the following steps:
step 1), paving a film on the surface of the area where the crack is located, and ensuring that the surface of the crack is completely covered by the film; pouring liquid on the film in the crack, filling the liquid in the area A formed by enclosing the film in the whole crack, and isolating the liquid from the crack through the film; aligning a ground penetrating radar to a crack at the near surface of the area A, and scanning at intervals from one end to the other end along the vertical direction of the row of cracks to obtain electromagnetic wave data of a plurality of sections of the underground form of the crack at the present time; drawing out, storing and removing the film from the liquid in the area A;
step 2) carrying out data processing on all electromagnetic wave data in the current period obtained in the step 1), and drawing by a visual interpretation method to obtain the outlines of cracks at each section in the current period;
step 3) drawing the contour of the crack at each section in the current period obtained in the step 2) in three-dimensional image processing software, and constructing a three-dimensional model of the underground form of the crack in the current period in an interpolation mode so as to obtain the three-dimensional form of the underground form of the crack in the current period and relevant data thereof;
and 4) continuously repeating the steps 1) to 3) according to the set frequency to obtain the multi-stage underground three-dimensional form and related data thereof in the dynamic change process of the crack development.
Further, the liquid in step 1) is selected to avoid interference with radar wave energy.
Further, in step 1), the ground penetrating radar is not in contact with the liquid. Optionally, the ground penetrating radar is kept out of contact with the liquid by laying a sheet on the membrane.
Further, the sheet selects a sheet that avoids interference with radar wave energy.
In the step 1), when interval scanning is performed from one end to the other end along the line seam trend, a marking line is made at intervals of 10cm in the vertical direction of the seam trend, and the ground penetrating radar is used for scanning from the position of 30cm on one side of the seam to the position of 30cm on the other side of the seam along the marking line.
Further, in step 2), the data processing includes background de-drying, one-dimensional filtering, wavelet transformation, and gain operation.
Further, in step 3), the correlation data includes an average width, an average depth, a surface area, and a volume of the fracture.
The invention has the characteristics and beneficial effects that:
the method for acquiring the multi-stage underground three-dimensional form of the dynamic crack based on the ground penetrating radar solves the defect that the underground form of the crack can damage the ground crack structure by injecting gypsum slurry into the crack, simultaneously avoids the defect that the required engineering quantity is extremely large and time and labor are consumed when a ground penetrating radar technology is adopted for excavating a working section at one side of the ground crack, ensures that the crack cannot be penetrated into soil by injecting liquid which can be extracted and cannot damage the underground form into the crack by using a film, ensures that the ground penetrating radar cannot be contacted with the liquid, enables the ground penetrating radar to scan, has the characteristics of simple flow, easiness in realization, no influence on the ground surface environment, extremely small required engineering quantity, time saving and labor saving, does not influence the form of the crack and the subsequent development process, can be used for nondestructively modeling the underground form of the dynamic crack, acquires the form parameters such as the average width, the average depth, the surface area and the volume of the crack, and the like, and has important significance on the development process and rule of the change of the three-dimensional form of the dynamic crack along with time.
Drawings
Fig. 1 (a) - (d) are three-dimensional, front, side and top views, respectively, of a method embodiment of the present invention in step 1) of laying a film over a ground fracture and filling with liquid;
fig. 2 (a) - (d) are a three-dimensional view, a front view, a side view and a top view of the ground penetrating radar scanning mode in step 2) of the method embodiment of the present invention, respectively;
FIG. 3 is an electromagnetic wave image of a ground penetrating radar scan in an embodiment of the present invention;
FIG. 4 is a cross-sectional view of a fracture subsurface morphology depicted using visual interpretation in an embodiment of the present invention;
FIG. 5 is a cross-sectional view of a fracture subsurface morphology of different scan profiles obtained in an embodiment of the present invention;
fig. 6 is a schematic diagram of a multi-stage three-dimensional model in the development process of a dynamic crack constructed according to an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the detailed description is presented by way of example only and is not intended to limit the scope of the invention.
In order to better understand the present invention, an application example of the method for acquiring the multi-period underground three-dimensional form of the dynamic fracture based on the ground penetrating radar is described in detail below.
The embodiment of the invention discloses a method for acquiring a multi-period underground three-dimensional form of a dynamic crack based on a ground penetrating radar, which specifically comprises the following steps:
step 1) acquiring electromagnetic wave data of an area where a crack is located by utilizing ground penetrating radar scanning
Referring to fig. 1, a film large enough to cover the surface of the crack is laid on the surface of the area where the crack is located, and the crack surface is ensured to be completely covered by the film; pouring liquid on the film in the crack to enable the liquid to fill the area A formed by enclosing the film in the whole crack, and isolating the liquid from the crack through the film; the method comprises the steps that liquid with small interference degree to radar wave energy, such as water, oil, carbohydrate and the like, is selected, a solvent-resistant, friction-resistant, tear-resistant and liquid-sealed film is selected as the film, a ground penetrating radar is aligned to a crack at the near surface of an area A, interval scanning is carried out from one end to the other end along the trend of a row of cracks, and electromagnetic wave data of a plurality of sections of the underground form of the crack at the present period are obtained by taking the vertical direction of the trend of the crack as a scanning path; the liquid in zone a is then pumped out for storage and the membrane is removed.
In this embodiment, the specific scanning process is as follows:
referring to fig. 2, in order to ensure that the liquid in the crack is not contacted with the ground penetrating radar so as to avoid damaging the ground penetrating radar, a thin plate (other ways of isolating the liquid and the ground penetrating radar can be used besides isolating the liquid and the ground penetrating radar by using the thin plate) is paved on the plastic film in a scanning range, the thin plate can be a plastic plate or other polymeric thin plates which do not interfere with the scanning result of the ground penetrating radar, a mark line is made at intervals of 10cm in the vertical direction of the crack trend, and the high-frequency ground penetrating radar is used for scanning from a position of 30cm on one side of the crack to a position of 30cm on the other side along the mark line, so that electromagnetic wave data of the current crack subsurface morphology trend to a plurality of sections is obtained.
In the embodiment, the ground penetrating radar adopts 1.6GHz high-frequency ground penetrating radar of coal resources of China mining university (Beijing) and key laboratories of safe exploitation countries.
Step 2) data processing is carried out on all electromagnetic wave data in the current period obtained in the step 1), and the outline of the crack at each section in the current period is obtained through drawing by a visual interpretation method
And (3) carrying out data processing on the image data of the area where the crack obtained in the step (1) is located by utilizing data processing software. In this embodiment, the data processing software adopts a GR radar processing analysis system developed by the university of chinese mining (beijing). And (3) performing conventional background drying, one-dimensional filtering, wavelet transformation, gain and other operations on each image data obtained in the step 1) to obtain an electromagnetic wave image with an interference area as shown in fig. 3.
Referring to fig. 4, the profile of the cross-sectional crack is determined from the interference area in the electromagnetic wave image, and the profile of the crack at each cross-section at the present time is depicted by visual interpretation.
Step 3) SolidWorks-based three-dimensional morphological modeling
Referring to fig. 5, the contour of the crack at each section in the current period obtained in the step 2) is drawn in the SolidWorks software of the three-dimensional image processing software according to the ratio of 1:1, and a three-dimensional model of the underground form of the crack in the current period is constructed in an interpolation mode, so that the underground three-dimensional form of the crack in the current period and relevant data thereof are obtained.
And 4) continuously repeating the steps 1) to 3) according to the set frequency to obtain the multi-stage underground three-dimensional form and related data thereof in the dynamic change process of the crack development.
Referring to fig. 6, three-phase underground three-dimensional morphology is obtained in this example, and the morphology parameters are shown in table 1, and the volume accuracy is the volume of the fracture and the volume of the liquid injected into the fracture.
TABLE 1 crack morphology parameters
As can be seen from Table 1, the method of the invention can obtain the three-dimensional morphological modeling volume of the crack without damage, the error is about 20%, and the main reasons for the error are the error of the scanned image of the ground penetrating radar and the error in the three-dimensional modeling process, if more accurate depth information is needed, the number of layers can be scanned by encryption or more accurate ground penetrating radar equipment and image processing software can be used.
The invention and its embodiments have been described above schematically, without limitation, and the invention is illustrated in the drawings as one of its embodiments and is not limited to practice. Therefore, those skilled in the art, having benefit of this disclosure, will appreciate that other embodiments can be devised which do not depart from the gist of the invention.

Claims (8)

1. The method for acquiring the multi-period underground three-dimensional form of the dynamic crack based on the ground penetrating radar is characterized by comprising the following steps of:
step 1), paving a film on the surface of the area where the crack is located, and ensuring that the surface of the crack is completely covered by the film; pouring liquid on the film in the crack, filling the liquid in the area A formed by enclosing the film in the whole crack, and isolating the liquid from the crack through the film; aligning a ground penetrating radar to a crack at the near surface of the area A, and scanning at intervals from one end to the other end along the trend of the row of cracks to obtain electromagnetic wave data of a plurality of sections of the underground form of the crack at the present time; drawing out, storing and removing the film from the liquid in the area A;
step 2) carrying out data processing on all electromagnetic wave data in the current period obtained in the step 1), and drawing by a visual interpretation method to obtain the outlines of cracks at each section in the current period;
step 3) drawing the contour of the crack at each section in the current period obtained in the step 2) in three-dimensional image processing software, and constructing a three-dimensional model of the underground form of the crack in the current period in an interpolation mode so as to obtain the three-dimensional form of the underground form of the crack in the current period and relevant data thereof;
and 4) continuously repeating the steps 1) to 3) according to the set frequency to obtain the multi-stage underground three-dimensional form and related data thereof in the dynamic change process of the crack development.
2. The method of claim 1, wherein the liquid in step 1) is selected to avoid interference with radar wave energy.
3. The method according to claim 1, wherein in step 1), the ground penetrating radar is not in contact with the liquid.
4. A method according to claim 3, wherein the ground penetrating radar is kept out of contact with the liquid by laying a sheet on the film.
5. The method of claim 4, wherein the sheet selects a sheet that avoids interference with radar wave energy.
6. The method according to claim 1, wherein in the step 1), when scanning at intervals from one end to the other end along the line direction of the slit, a mark line is made at intervals of 10cm in the vertical direction of the slit direction, and the ground penetrating radar is used for scanning from a position of 30cm on one side of the slit to a position of 30cm on the other side along the mark line.
7. The method of claim 1, wherein in step 2) the data processing comprises background de-drying, one-dimensional filtering, wavelet transformation, and gain manipulation.
8. The method according to any one of claims 1 to 7, wherein in step 3) the relevant data comprises an average width, an average depth, a surface area and a volume of the crack.
CN202011566874.5A 2020-12-25 2020-12-25 Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar Active CN112799054B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011566874.5A CN112799054B (en) 2020-12-25 2020-12-25 Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011566874.5A CN112799054B (en) 2020-12-25 2020-12-25 Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar

Publications (2)

Publication Number Publication Date
CN112799054A CN112799054A (en) 2021-05-14
CN112799054B true CN112799054B (en) 2024-03-22

Family

ID=75804948

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011566874.5A Active CN112799054B (en) 2020-12-25 2020-12-25 Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar

Country Status (1)

Country Link
CN (1) CN112799054B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113433547A (en) * 2021-07-05 2021-09-24 长沙理工大学 Ground penetrating radar hidden crack offset imaging method, system, terminal and medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102590874A (en) * 2012-01-16 2012-07-18 中国矿业大学(北京) Method for detecting ground surface crack of upland coal-mining subsidence paddy field
CN102944195A (en) * 2012-11-28 2013-02-27 水利部交通运输部国家能源局南京水利科学研究院 Method for detecting depth of crack
CN104267449A (en) * 2014-09-15 2015-01-07 中国神华能源股份有限公司 Method for detecting coal mining earth surface crack patterns
CN108519596A (en) * 2018-03-05 2018-09-11 长沙理工大学 A kind of section of jurisdiction Hidden Fracture recognition methods based on match tracing and wavelet transformation
CN109696674A (en) * 2019-02-14 2019-04-30 中国矿业大学(北京) A method of ground fissure form is obtained based on Ground Penetrating Radar

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10392929B1 (en) * 2018-02-08 2019-08-27 Saudi Arabian Oil Company Mapping fracture length using downhole ground penetrating radar

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102590874A (en) * 2012-01-16 2012-07-18 中国矿业大学(北京) Method for detecting ground surface crack of upland coal-mining subsidence paddy field
CN102944195A (en) * 2012-11-28 2013-02-27 水利部交通运输部国家能源局南京水利科学研究院 Method for detecting depth of crack
CN104267449A (en) * 2014-09-15 2015-01-07 中国神华能源股份有限公司 Method for detecting coal mining earth surface crack patterns
CN108519596A (en) * 2018-03-05 2018-09-11 长沙理工大学 A kind of section of jurisdiction Hidden Fracture recognition methods based on match tracing and wavelet transformation
CN109696674A (en) * 2019-02-14 2019-04-30 中国矿业大学(北京) A method of ground fissure form is obtained based on Ground Penetrating Radar

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
定向井压裂裂缝三维扩展形态的可视化仿真;程远方;杨柳;吴百烈;董丙响;;计算机仿真;20121215(第12期);全文 *
形态参数下地表采动裂缝三维建模及可视化方法;戴华阳;王祥;李军;郭俊廷;阎跃观;;测绘通报;20180225(第02期);全文 *
砂煤岩互层水力裂缝扩展规律的数值模拟研究;兰天庆;胡泊洲;董文楠;张昕;;能源与环保;20181102(第10期);全文 *
采动地表裂缝三维形态探测方法及精度评价研究;张健;毕银丽;彭苏萍;;煤炭科学技术;20200915;48(09);全文 *

Also Published As

Publication number Publication date
CN112799054A (en) 2021-05-14

Similar Documents

Publication Publication Date Title
WO2020248817A1 (en) Water supply pipeline leakage detection method based on three-dimensional image attribute analysis using ground penetrating radar
CN105239611B (en) Determine the method that water stopping curtain seepage influences on surrounding enviroment below excavation of foundation pit face
CN103290836B (en) Vacuum decompression compacting and consolidating construction method for soft foundation treatment
CN101906771B (en) Three-dimensional soil texture analyzing method for assisting survey, design and construction of dredging engineering
Wu et al. The pixel crack reconstruction method: From fracture image to crack geological model for fracture evolution simulation
CN112799054B (en) Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar
CN110275204B (en) Quantitative analysis method for horizontal slip distance of underground slip fracture
CN106501090A (en) Crack characterizing method for hydraulic fracturing simulation experiment
CN102590874A (en) Method for detecting ground surface crack of upland coal-mining subsidence paddy field
CN110987765A (en) Rock mass microscopic fracture testing method based on three-dimensional digital speckle
CN108959681A (en) A kind of barrier lake sluice way design method based on point cloud data
CN112817057A (en) Method for economically, quickly and accurately detecting underground space distribution characteristics of landfill pond
CN109696674B (en) Method for acquiring ground fracture form based on ground penetrating radar
CN110019594B (en) Geological knowledge base building method based on numerical simulation
CN110363855B (en) Rock-fill dam transparentization modeling method
CN103091471A (en) Model for promoting solidification of hydraulic fill and test method for model
CN106680454A (en) Method for calculating soil erosion modulus of governed collapse mound with sand-blocking dam
CN105297752A (en) Method for determining water pumping quantity of dewatering well when fixed-water-level water pumping is carried out under action of water-proof curtain
CN106338257A (en) Method for quantitatively measuring embedding depth of proppant
CN109143381B (en) A kind of loess high roadbed aquitard detection method
CN114609682A (en) Method and system for extracting parameters of electrical imaging seam hole and computer readable storage medium
CN113722663B (en) Method for evaluating underwater dredging and mud throwing of limited space
Gillespie et al. A 3d model of complex A, la Venta, Mexico
CN113592823B (en) Rock stratum porosity distribution calculation method
Meng et al. Quantifying vein attributes in massive mudstones (Triassic, SW England): implications for progressive evolution of opening-mode fracture networks

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