CN114488346B - Underground space abnormal body detection method, device, equipment and medium - Google Patents
Underground space abnormal body detection method, device, equipment and medium Download PDFInfo
- Publication number
- CN114488346B CN114488346B CN202210108195.6A CN202210108195A CN114488346B CN 114488346 B CN114488346 B CN 114488346B CN 202210108195 A CN202210108195 A CN 202210108195A CN 114488346 B CN114488346 B CN 114488346B
- Authority
- CN
- China
- Prior art keywords
- data
- time period
- theoretical data
- underground space
- sampling time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
The invention provides a method, a device, equipment and a medium for detecting an underground space abnormal body, wherein the method comprises the following steps: acquiring theoretical data of a preset wave velocity model, and correcting parameters of the theoretical model according to the measured data to obtain an optimal value of the wave velocity model of the measured field; acquiring theoretical data in different time periods under a preset apparent resistivity model, and correcting the theoretical model according to actually measured data in a corresponding time period to obtain an optimal value of the apparent resistivity model of the actually measured field; correcting the conversion parameter model according to the multi-field underground optimal wave velocity and apparent resistivity parameter set, and solving the problem of multiple solutions of abnormal characteristics; by adopting the embodiment of the invention, the optimized seismic mapping method of double-channel acquisition is adopted to obtain the wave field signal of the underground space, and the array type transient electromagnetic method is utilized to obtain the vector data of the electrical structure of the stratum, thereby further improving the detection accuracy of the abnormal body in the underground space.
Description
Technical Field
The invention relates to the technical field of underground space exploration, in particular to a method, a device, equipment and a medium for detecting an underground space abnormal body.
Background
The underground space is a space developed, constructed and utilized below the ground surface to meet the requirements of production, life, traffic, environmental protection, energy, safety, disaster prevention, disaster reduction and the like of human society in urban and rural planning areas, and generally comprises basements, underground workshops, underground pipe galleries, underground passages, civil air defense projects and the like. The geological conditions under the earth surface or in deep structures and the like all affect underground space facilities and buildings, the factors such as artificial construction disturbance and movement of structural belts cause settlement and collapse of the ground, damage and deformation of the underground space and the like, and at present, a relatively perfect underground space detection technical system is formed.
However, the inventor of the present invention found in the research on the prior art that, currently, the prior art uses a single elastic wave field, and the detection accuracy of the underground space abnormal body is low due to the lack of electrical characteristics.
Disclosure of Invention
The invention provides a method, a device, equipment and a medium for detecting an abnormal body in an underground space, which can further improve the detection accuracy of the abnormal body in the underground space.
In order to achieve the above object, an embodiment of the present invention provides a method for detecting an underground spatial anomaly, including the following steps:
acquiring first theoretical data and first measured data which change along with the wave velocity in a preset sampling time period, and correcting the first theoretical data according to the first measured data to obtain first optimal theoretical data in the preset sampling time period; the first theoretical data are obtained according to an underground space wave impedance model established by wave impedance inversion, and the first measured data are high-frequency band signals and low-frequency band signals of an abnormal body acquired based on the improved dual-channel detector in the preset sampling time period;
acquiring second theoretical data and second measured data which change along with the apparent resistivity in the preset sampling time period, and correcting the second theoretical data according to the second measured data to obtain second optimal theoretical data in the preset sampling time period; the second theoretical data are obtained through an apparent resistivity model of the underground space established through smoke ring inversion, and the second measured data are induced electromotive force obtained based on the improved array type receiving coil in the preset sampling time period;
and substituting the first optimal theoretical data and the second optimal theoretical data into a preset wave velocity-apparent resistivity conversion model to obtain a detection result of the abnormal body.
As an alternative embodiment, the wave impedance model is established by:
based on the cavity structure of the underground space, the dual-channel detector is used for collecting high-frequency band signals and low-frequency band signals of the abnormal body, the impedance information of the reflection interface wave of the surrounding rock body, the cavity modal response characteristic and the structure frequency characteristic of the abnormal body are calculated according to the time-frequency characteristics of the signals, and the corresponding underground space wave impedance model of the collected signals is built according to the wave impedance inversion.
As an alternative embodiment, the apparent resistivity model is established by:
and acquiring abnormal signals of the reinforcing mesh structure of the peripheral cavity of the underground space based on the improved symmetrically distributed receiving coil system, and establishing an apparent resistivity model of the surrounding rock body and the underground space according to apparent resistivity information of the smoke ring inversion stratum structure.
As an optional embodiment, the modifying the first theoretical data according to the first measured data to obtain the first optimal theoretical data in the preset sampling time period includes:
correcting the first theoretical data of each sampling moment in the preset sampling time period according to the first measured data of the moment to obtain corrected data of the moment;
and acquiring a correction data set in the preset sampling time period, and taking the minimum value in the correction data set as the first optimal theoretical data in the preset sampling time period.
As an optional embodiment, the modifying the first theoretical data at the time according to the first measured data at the time to obtain modified data at the time includes:
and taking a value obtained by a two-norm difference between the first measured data and the first theoretical data at the moment as correction data at the moment.
As an optional embodiment, the modifying the second theoretical data according to the second measured data to obtain second optimal theoretical data in the preset sampling time period includes:
for each sampling moment in the preset sampling time period, correcting second theoretical data of the moment according to second measured data of the moment to obtain corrected data of the moment;
and acquiring a correction data set in the preset sampling time period, and taking the minimum value in the correction data set as second optimal theoretical data in the preset sampling time period.
As an alternative embodiment, the expression of the preset wave velocity-apparent resistivity conversion model is as follows:
wherein K, C, d are all constraint parameters, V is the first theoretical data, and R is t For the second theoretical data, H is depth.
Correspondingly, another embodiment of the present invention provides a device for detecting an underground space abnormal body, including:
the wave impedance model optimal value output module is used for acquiring first theoretical data and first measured data which change along with the wave speed in a preset sampling time period, and correcting the first theoretical data according to the first measured data to obtain first optimal theoretical data in the preset sampling time period; the first theoretical data are obtained according to an underground space wave impedance model established by wave impedance inversion, and the first measured data are high-frequency band signals and low-frequency band signals of an abnormal body acquired based on the improved dual-channel detector in the preset sampling time period;
the apparent resistivity model optimal value output module is used for acquiring second theoretical data and second measured data which change along with the apparent resistivity in the preset sampling time period, and correcting the second theoretical data according to the second measured data to obtain second optimal theoretical data in the preset sampling time period; the second theoretical data are obtained through an apparent resistivity model of the underground space established through smoke ring inversion, and the second measured data are induced electromotive force obtained based on the improved array type receiving coil in the preset sampling time period;
and the abnormal body detection result acquisition module is used for substituting the first optimal theoretical data and the second optimal theoretical data into a preset wave velocity-apparent resistivity conversion model to obtain a detection result of the abnormal body.
Correspondingly, another embodiment of the present invention provides a terminal device, which includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, where the processor executes the computer program to implement the underground space abnormal body detection method according to the above-described embodiment of the present invention.
Another embodiment of the present invention provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program, where when the computer program runs, a device in which the computer-readable storage medium is located is controlled to execute the method for detecting an underground space abnormal body according to the above-mentioned embodiment of the present invention.
Compared with the prior art, the method, the device, the equipment and the medium for detecting the underground space abnormal body provided by the embodiment of the invention have the advantages that the wave field signal of the underground space is obtained by adopting the double-channel detector excited by the artificial seismic source, the electrical structure of the stratum is analyzed by utilizing the array type transient electromagnetic method, the unique structure of the underground space which is formed by pouring reinforced concrete, is similar to a cavity and has a complete interface is combined, the wave velocity and the electrical characteristics reflected by the elastic wave field and the electromagnetic field are jointly analyzed, the azimuth determination and the scale detection of the underground space are facilitated, and the detection accuracy of the abnormal body in the underground space is further improved.
Drawings
FIG. 1 is a schematic flow chart of a method for detecting an underground spatial anomaly according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a dual channel detector according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of an array type transient electromagnetic acquisition device according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of an underground space abnormal body detecting device according to an embodiment of the present invention;
fig. 5 is a schematic structural diagram of a terminal device according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
Referring to fig. 1, a schematic flow chart of a method for detecting an underground spatial anomaly according to an embodiment of the present invention is shown, where the method includes steps S11 to S13:
s11, acquiring first theoretical data and first measured data which change along with wave velocity in a preset sampling time period, and correcting the first theoretical data according to the first measured data to obtain first optimal theoretical data in the preset sampling time period; the first theoretical data are obtained according to an underground space wave impedance model established by wave impedance inversion, and the first measured data are high-frequency band signals and low-frequency band signals of an abnormal body acquired based on the improved dual-channel detector in the preset sampling time period;
s12, second theoretical data and second measured data which change along with the apparent resistivity in the preset sampling time period are obtained, and the second theoretical data are corrected according to the second measured data to obtain second optimal theoretical data in the preset sampling time period; the second theoretical data are obtained through an apparent resistivity model of the underground space established through smoke ring inversion, and the second measured data are induced electromotive force obtained based on the improved array type receiving coil in the preset sampling time period;
and S13, substituting the first optimal theoretical data and the second optimal theoretical data into a preset wave velocity-apparent resistivity conversion model to obtain a detection result of the abnormal body.
It can be understood that, at present, the prior art adopts a single elastic wave field, which can effectively identify the interfaces of different structural layers of the stratum and has higher longitudinal resolution for integrity analysis, but due to the lack of electrical characteristics, the filling materials in the cavity structure can not be distinguished from sludge, water or cavities. The prior art techniques for detecting subterranean spaces may suffer from the following disadvantages:
(1) A single electromagnetic field is used. The ground penetrating radar cannot detect when meeting the conditions of ponding road surfaces, shallow steel plate and other good conductor coverage, the conventional transient electromagnetic method has insufficient transverse resolution, and when meeting shallow and small-size target bodies, the target bodies cannot be identified due to skin effect and volume effect;
(2) A single electric field is used. The conduction direct current and high-density resistivity methods are good at cavity structures such as metal pipes, water-rich karsts and the like, but the results of the distribution of the electrical characteristics of the stratum are obviously different due to different devices, so that false abnormity misjudgment is easily caused;
(3) A single elastic wave field is used. The interface of different structural layers of the stratum can be effectively identified, and the longitudinal resolution is high for integrity analysis; but the electric characteristics are lost, so that the filler in the cavity structure can not be distinguished to be sludge, water or a cavity, and a pipeline and a cavity can not be identified;
(4) And (6) micro-motion. The method is only suitable for underground space scale detection of surrounding spaciousness, interference-free regional performance and long-time scale observation, has higher requirements on stratum environment, rock stratum characteristics, surrounding environment and the like, and is difficult to meet the detection requirements of underground space between cities.
Compared with the prior art, the method for detecting the abnormal bodies in the underground space provided by the embodiment of the invention has the advantages that the wave field signals of the underground space are obtained by adopting the dual-channel detector excited by the artificial seismic source, the array type transient electromagnetic method is utilized for carrying out stratum electrical structure analysis, the unique structure of the underground space which is formed by pouring reinforced concrete and is shaped like a cavity and has a complete interface is combined, the wave velocity and the electrical characteristics reflected by the elastic wave field and the electromagnetic field are jointly analyzed, the azimuth determination and the scale detection of the underground space are facilitated, and the detection accuracy of the abnormal bodies in the underground space is further improved. By providing the detection method which integrates the elastic wave field and the time domain electromagnetic field and is suitable for the underground space, the technical method which can effectively improve the signal-to-noise ratio and the abnormal resolution of the underground space observation signal represented by the underground pipe gallery is provided, and the method is suitable for detecting the urban underground space with a complex ground surface structure and has the characteristics of high efficiency, no damage, effectiveness and the like.
As an alternative embodiment, the wave impedance model is established by:
based on the cavity structure of the underground space, the dual-channel detector is used for collecting high-frequency band signals and low-frequency band signals of the abnormal body, the impedance information of the reflection interface wave of the surrounding rock body, the cavity modal response characteristic and the structure frequency characteristic of the abnormal body are calculated according to the time-frequency characteristics of the signals, and the corresponding underground space wave impedance model of the collected signals is built according to the wave impedance inversion.
It can be understood that the embodiment of the invention provides a novel device, and the device improves the mode of point-to-station acquisition of a single seismic source and a single detector of a seismic image at equal offset distances by improving the detectors in the seismic image and utilizing the cavity structure characteristics of an underground pipe gallery. The method specifically comprises the following steps: the method comprises the steps of surrounding a unique cavity structure of an underground space, utilizing two channel detectors at equal intervals to be arranged on the ground surface in the direction of a vertical line, collecting high-frequency band signals and low-frequency band signals of surrounding rocks and a target body, calculating reflected interface wave impedance information of the surrounding rocks, cavity modal response characteristics and structural frequency characteristics of the target body according to time-frequency characteristics of the signals, and establishing a wave impedance model of the underground space.
Fig. 2 is a schematic structural diagram of a dual-channel detector according to an embodiment of the present invention. 21 is a hammer source, 22 is a hammer pad, 23 is an instrument, 24 is a transmitting source receiving point, 25 and 26 are receiving sources, 27 is a high-frequency band detector, 28 is a low-frequency band detector, and the specific process of collecting signals is as follows:
(1) Arranging a seismic source at a specified position, and respectively placing a high-frequency detector 27 and a low-frequency detector 28 at a position which is horizontally separated from the seismic source by L1;
(2) Arranging two detectors perpendicular to the ground surface and the measuring line at equal intervals L2 to form determinant arrangement, and receiving signals from the stratum in parallel;
(3) The method comprises the steps that a hammer source 21 is used for knocking a hammer pad 22 to excite a seismic source, a low-frequency band detector 25 obtains a wide-band signal with certain signal intensity by adjusting gain, a high-frequency band detector 26 obtains a high-sensitivity signal with certain bandwidth through a high-pass filter and a compensation circuit, a high-frequency receiving source 26 and a low-frequency receiving source 25 respectively receive seismic records from the two detectors, signals of the two detectors are input to an instrument at the same time and in different channels through cables, single-point excitation is achieved, seismic signals of high and low separation frequency bands are synchronously obtained, and two groups of vibration signals are obtained; before the seismic source is excited next time, the seismic source is horizontally moved to the right by L3 towards the direction of the detector, the relative positions of the seismic source and the detector are kept unchanged, and then a vibration signal set of the whole measuring line is acquired;
(4) The offset distance and the track distance are determined according to the size and the buried depth of the target structure. Synchronously moving and performing seismic record acquisition at equal offset distance to obtain two groups of signals;
(5) And after the two groups of signal sets are subjected to static correction, amplitude compensation and difference inversion and the like, analyzing the shape and the size of the structural mode set of the target abnormal body.
It is worth to be noted that, through a mode of synchronous acquisition of high and low frequency channels, two detectors with different bandwidths are used for acquisition, a broadband signal record with certain signal intensity and a high-sensitivity signal with certain bandwidth can be simultaneously acquired, and the modal response of the cavity is analyzed through a frequency domain, so that a new explanation thought and direction are provided for the structural form of an abnormal body, and more reliable data are provided for the refined analysis of the abnormality.
As an alternative embodiment, the apparent resistivity model is established by:
and acquiring abnormal signals of the reinforcing mesh structure of the peripheral cavity of the underground space based on the improved symmetrically distributed receiving coil system, and establishing an apparent resistivity model of the surrounding rock body and the underground space according to the apparent resistivity information of the stratal structure inversed by the smoke ring.
Illustratively, the embodiment of the invention provides a novel device, which improves the mode of a single coil transmitting signal and a single receiving coil receiving signal of the conventional transient electromagnetic method by improving the receiving coil of the transient electromagnetic method. Fig. 3 is a schematic structural diagram of an array transient electromagnetic acquisition device according to an embodiment of the present invention. Because the receiving and the transmitting are integrated, the device is vertically placed on the ground surface, 30 is a framework, 31 and 32 are receiving coils, 33 is a transmitting coil, 34 and 35 receiving points of the transmitting source are 34, 36, 37, 38 and 39 are all receiving points of the receiving coil, and 40 is an instrument host. The specific process of acquisition is as follows:
(1) Selecting a long cylindrical plastic pipe with moderate size (the size is 1m in radius and 2m in length) as a framework 40;
(2) Two receiving coils are distributed above the earth surface in parallel and in equal difference with the earth surface, a transmitting coil is arranged flatly on the earth surface, and three coils are fixed on a long cylindrical plastic pipe at equal intervals to form a new collecting device. 2 receiving coils 31 and 32,1 transmitting coils 33 which have the same number of turns, the same material and the same specification are wound on the direct head of the long cylindrical plastic pipe;
(3) The transmitter coil 33 is placed in close proximity to the earth's surface to form a central loop-like transmission pattern.
(4) The receiving coil systems are distributed at equal intervals in space to form the array type transient electromagnetic method receiving mode to obtain two groups of signals; the receiving coils 31 and 32 adopt a multi-coil receiving mode to ensure that the centers of the straight joints of the receiving coils are respectively placed at positions 1m and 2m away from the earth surface to form an array receiving mode and obtain two groups of signals;
(5) The signals of the two receiving coils 31 and 32 are respectively received by using the multi-channel receiving points 38 and 39, and 36 and 37 of the instrument 40, differential calculation is carried out to obtain gradient field signals and vector characteristics of mutual induction signals under different time sequences, direction judgment and range delineation are carried out on the resistance characteristic of a front target body according to proper characteristics of forward attenuation of the mutual induction signals and enhancement of the gradient signals of the target body in induction signals, and the junction electrical structure of the target body is analyzed.
It is worth to be noted that, by adopting the improved scheme of the embodiment, the transient electromagnetic response signal can have a gradient characteristic, the transmitting coil source has a better earth electric characteristic close to the earth surface, the receiving coil system and the transmitting coil present equal difference distribution in space, the primary field signal of the source coil in space and the secondary field signal of the target body are received, and differential calculation is performed to obtain gradient values of different spatial dimensions of the surrounding rock and the target body, so that the change rule of the spatial gradient of the target body is realized. In addition, the smaller the distance between the hidden trouble target and the hidden trouble target, the larger the distance between the hidden trouble target and the hidden trouble target, the smaller the gradient value of the secondary field, and the fixed distance of the transmitting and receiving coils is utilized to greatly reduce the early-stage interference of the primary field, thereby providing more accurate and fine analysis for the abnormity of the hidden trouble body.
As an optional embodiment, the correcting the first theoretical data according to the first measured data to obtain first optimal theoretical data in the preset sampling time period includes:
correcting the first theoretical data of each sampling moment in the preset sampling time period according to the first measured data of the moment to obtain corrected data of the moment;
and acquiring a correction data set in the preset sampling time period, and taking the minimum value in the correction data set as the first optimal theoretical data in the preset sampling time period.
As an optional embodiment, the modifying the first theoretical data at the time according to the first measured data at the time to obtain modified data at the time includes:
and taking a value obtained by a two-norm difference between the first measured data and the first theoretical data at the moment as correction data at the moment.
As an optional embodiment, the modifying the second theoretical data according to the second measured data to obtain second optimal theoretical data in the preset sampling time period includes:
correcting the second theoretical data of each sampling moment in the preset sampling time period according to the second measured data of the sampling moment to obtain corrected data of the sampling moment;
and acquiring a correction data set in the preset sampling time period, and taking the minimum value in the correction data set as second optimal theoretical data in the preset sampling time period.
Illustratively, the modified specific expression is:
A(v,t)＝‖B(v,t)-C(v,t)‖ 2 (1)
wherein, a (v, t) is a corresponding data set of first theoretical data or a corresponding data set of second theoretical data in the optimal solution process, B (v, t) is actually measured data of the first theoretical data or actually measured data of the second theoretical data actually measured in the field, and C (v, t) is theoretical data of the first theoretical data or theoretical data of the second theoretical data calculated through a model.
It should be noted that, in the process of obtaining the optimal solution, iteration is performed continuously, and after the optimal solution is obtained, iteration is not stopped, and then iteration is continued until an iteration threshold is reached, and calculation is stopped. If the optimal solution is not obtained after the iteration threshold is reached, iteration is continued, and whether the numerical result meets the requirement of the optimal solution or not is calculated after the iteration is not performed once until the optimal solution is obtained. Wherein the optimal solution may be a value that satisfies a predetermined error range.
As an alternative embodiment, the expression of the preset wave velocity-apparent resistivity conversion model is as follows:
wherein K, C, d are all constraint parameters, V is first theoretical data, and R t For the second theoretical data, H is depth.
Illustratively, after obtaining the optimal parameter of the first theoretical data or the second theoretical data in a period of time through formula (1), substituting the optimal parameter into formula (2) to obtain the real geoelectric structural model of the target abnormal body.
Specifically, the step S13 may be:
after the optimal solution of the second theoretical data within a period of time is obtained, the optimal solution is substituted into the formula (2), the first theoretical data is the theoretical data calculated according to the wave impedance model, and the optimized constraint parameter K, C, d is obtained through a single variable method so as to optimize the formula (2). In a popular way, one of the parameters is brought into the continuous optimization expression (2), namely the constraint coefficient is calibrated until the current optimal formula is obtained, and for an abnormal body, the abnormal body may be a water-containing cavity and a mud-containing cavity, the characteristics of the wave velocity representation may be the same, but the actual values are different), and various parameters are represented, so that the theoretical value needs to be corrected according to the measured value.
It can be understood that according to parameters under the optimal model, the first theoretical data changing with wave velocity and the second theoretical data changing with apparent resistivity are output through the formula (2), imaging analysis is carried out, a real earth electricity and structure model of the target body is obtained, after a plurality of groups of first theoretical data and second theoretical data are output, a multi-parameter set of the underground structure is obtained, and the depth and the object properties of the abnormal body, such as a low-resistance body or a high-resistance body, a water-containing cavity or a mud-containing cavity, are obtained. And the multi-solution of joint inversion is reduced through parameter sets, the inversion result of transient electromagnetism and the seismic mapping result are improved, and the problem of multi-solution of abnormal characteristics is solved.
It is worth to be noted that wave field signals of the underground space are obtained through a seismic mapping method of using an artificial seismic source to excite and using a double-channel detector to receive in a detection area, then stratum electrical structure analysis is carried out through an array type transient electromagnetic receiving mode, finally, the underground space is combined with a unique structure which is formed by pouring reinforced concrete and is shaped like a cavity and provided with a complete interface, and the elastic wave field and wave speed and electrical characteristics reflected by an electromagnetic field are subjected to combined analysis, so that azimuth determination and scale detection of the underground space are facilitated. By adopting the technical method suitable for underground space detection provided by the embodiment of the invention, the underground space detection precision represented by the underground pipe gallery can be improved.
Referring to fig. 4, a schematic structural diagram of an underground space abnormal body detecting device 4 provided in an embodiment of the present invention includes:
the wave impedance model optimal value output module 41 is configured to acquire first theoretical data and first actually measured data that change with a wave speed within a preset sampling time period, and correct the first theoretical data according to the first actually measured data to obtain first optimal theoretical data within the preset sampling time period; the first theoretical data are obtained according to an underground space wave impedance model established by wave impedance inversion, and the first measured data are high-frequency band signals and low-frequency band signals of an abnormal body collected based on the improved two-channel detector in the preset sampling time period;
the apparent resistivity model optimal value output module 42 is configured to obtain second theoretical data and second measured data that change with the apparent resistivity within the preset sampling time period, and correct the second theoretical data according to the second measured data to obtain second optimal theoretical data within the preset sampling time period; the second theoretical data are obtained through an apparent resistivity model of the underground space established through smoke ring inversion, and the second measured data are induced electromotive force obtained based on the improved array type receiving coil in the preset sampling time period;
and an abnormal body detection result obtaining module 43, configured to substitute the first optimal theoretical data and the second optimal theoretical data into a preset wave velocity-apparent resistivity conversion model to obtain a detection result of the abnormal body.
Compared with the prior art, the underground space abnormal body detection device provided by the embodiment of the invention obtains the wave field signal of the underground space by adopting the dual-channel wave detector excited by the artificial seismic source, utilizes the array type transient electromagnetic method to analyze the electrical structure of the stratum, combines the unique structure of the underground space which is formed by pouring reinforced concrete and is shaped like a cavity and has a complete interface, and performs combined analysis on the wave velocity and the electrical characteristics reflected by the elastic wave field and the electromagnetic field, thereby being beneficial to azimuth determination and scale detection of the underground space and further improving the detection accuracy of the abnormal body in the underground space
In addition, it should be noted that for the specific description and the beneficial effects related to each embodiment of the underground space abnormal body detection device in this embodiment, reference may be made to the specific description and the beneficial effects related to each embodiment of the underground space abnormal body detection method described above, and details are not described herein again.
Fig. 5 is a schematic structural diagram of a terminal device according to an embodiment of the present invention. The terminal device 5 of this embodiment includes: a processor 50, a memory 51 and a computer program stored in said memory 51 and executable on said processor 50. The processor 50, when executing the computer program, implements the steps in the above-described embodiments of the control method for the vehicle-mounted atmosphere lamp. Alternatively, the processor 50 implements the functions of the modules in the above device embodiments when executing the computer program.
Illustratively, the computer program may be divided into one or more modules, which are stored in the memory 51 and executed by the processor 50 to accomplish the present invention. The one or more modules may be a series of computer program instruction segments capable of performing specific functions, which are used to describe the execution process of the computer program in the terminal device 5.
The terminal device 5 may be a desktop computer, a notebook, a palm computer, a cloud server, or other computing devices. The terminal device 5 may include, but is not limited to, a processor 50 and a memory 51. It will be appreciated by those skilled in the art that the schematic diagram is merely an example of a terminal device, and does not constitute a limitation of the terminal device, and may include more or less components than those shown, or combine some components, or different components, for example, the terminal device 5 may further include an input-output device, a network access device, a bus, etc.
The Processor 50 may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic, discrete hardware components, etc. The general purpose processor may be a microprocessor or the processor may be any conventional processor or the like, and the processor 50 is the control center of the terminal device 5 and connects the various parts of the whole terminal device 5 by means of various interfaces and lines.
The memory 51 may be used for storing the computer programs and/or modules, and the processor 50 implements various functions of the terminal device 5 by running or executing the computer programs and/or modules stored in the memory 51 and calling data stored in the memory 51. The memory 51 may mainly include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application program (such as a sound playing function, an image playing function, etc.) required by at least one function, and the like; the storage data area may store data (such as audio data, a phonebook, etc.) created according to the use of the cellular phone, etc. Further, the memory 51 may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
Wherein, the module integrated by the terminal device 5 can be stored in a computer readable storage medium if it is implemented in the form of software functional unit and sold or used as a stand-alone product. Based on such understanding, all or part of the flow in the method according to the above embodiments may be implemented by a computer program, which may be stored in a computer readable storage medium and used by the processor 50 to implement the steps of the above embodiments of the method. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, read-Only Memory (ROM), random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, and the like. It should be noted that the computer readable medium may contain content that is subject to appropriate increase or decrease as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, computer readable media does not include electrical carrier signals and telecommunications signals as is required by legislation and patent practice.
It should be noted that the above-described embodiments of the apparatus are merely illustrative, where the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. In addition, in the drawings of the embodiment of the apparatus provided by the present invention, the connection relationship between the modules indicates that there is a communication connection between them, and may be specifically implemented as one or more communication buses or signal lines. One of ordinary skill in the art can understand and implement it without inventive effort.
The embodiment of the invention also provides a computer-readable storage medium, which includes a stored computer program, wherein when the computer program runs, the apparatus where the computer-readable storage medium is located is controlled to execute the method for detecting an underground space abnormal body as described above.
Those skilled in the art will appreciate that the modules in the devices in the embodiments may be adaptively changed and arranged in one or more devices different from the embodiments. The modules or units in the embodiments may be combined into one module or unit, and furthermore, they may be divided into a plurality of sub-modules or sub-units. All of the features disclosed in this specification (including any accompanying claims, abstract and drawings), and all of the processes or elements of any method or apparatus so disclosed, may be combined in any combination, except combinations where at least some of such features and/or processes or elements are mutually exclusive. Each feature disclosed in this specification (including any accompanying claims, abstract and drawings) may be replaced by alternative features serving the same, equivalent or similar purpose, unless expressly stated otherwise.
It should be noted that the above-mentioned embodiments illustrate rather than limit the invention, and that those skilled in the art will be able to design alternative embodiments without departing from the scope of the appended claims. In the claims, any reference signs placed between parentheses shall not be construed as limiting the claim. The word "comprising" does not exclude the presence of elements or steps not listed in a claim. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements.
Furthermore, those skilled in the art will appreciate that while some embodiments described herein include some features included in other embodiments, rather than other features, combinations of features of different embodiments are meant to be within the scope of the invention and form different embodiments. For example, in the claims of the present invention, any of the claimed embodiments may be used in any combination.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.
Claims (9)
1. A method for detecting an underground space abnormal body is characterized by comprising the following steps:
acquiring first theoretical data and first measured data which change along with wave velocity in a preset sampling time period, and correcting the first theoretical data according to the first measured data to obtain first optimal theoretical data in the preset sampling time period; the first theoretical data are obtained according to an underground space wave impedance model established by wave impedance inversion, and the first measured data are high-frequency band signals and low-frequency band signals of an abnormal body acquired based on the improved dual-channel detector in the preset sampling time period;
acquiring second theoretical data and second measured data which change along with the apparent resistivity in the preset sampling time period, and correcting the second theoretical data according to the second measured data to obtain second optimal theoretical data in the preset sampling time period; the second theoretical data are obtained through an apparent resistivity model of the underground space established through smoke ring inversion, and the second measured data are induced electromotive force obtained based on the improved array type receiving coil in the preset sampling time period;
substituting the first optimal theoretical data and the second optimal theoretical data into a preset wave velocity-apparent resistivity conversion model to obtain a detection result of the abnormal body;
wherein the wave impedance model is established by:
based on the cavity structure of the underground space, the dual-channel detector is used for collecting high-frequency band signals and low-frequency band signals of the abnormal body, the wave impedance information of the reflection interface of the surrounding rock body, the cavity modal response characteristics and the structure frequency characteristics of the abnormal body are calculated according to the time-frequency characteristics of the signals, and the corresponding underground space wave impedance model of the collected signals is established according to wave impedance inversion.
2. A method of detecting a subterranean spatial anomaly according to claim 1, wherein said apparent resistivity model is established by:
and acquiring abnormal signals of the reinforcing mesh structure of the peripheral cavity of the underground space based on the improved symmetrically distributed receiving coil system, and establishing an apparent resistivity model of the surrounding rock body and the underground space according to apparent resistivity information of the smoke ring inversion stratum structure.
3. The method for detecting an underground space abnormal body according to claim 1, wherein the step of correcting the first theoretical data according to the first measured data to obtain the first optimal theoretical data in the preset sampling time period comprises:
correcting the first theoretical data of each sampling moment in the preset sampling time period according to the first measured data of the moment to obtain corrected data of the moment;
and acquiring a correction data set in the preset sampling time period, and taking the minimum value in the correction data set as the first optimal theoretical data in the preset sampling time period.
4. A method according to claim 3, wherein the step of correcting the first theoretical data at the time based on the first measured data at the time to obtain corrected data at the time comprises:
and taking a value obtained by a two-norm difference between the first measured data and the first theoretical data at the moment as correction data at the moment.
5. The method for detecting an underground space abnormal body according to claim 1, wherein the step of correcting the second theoretical data according to the second measured data to obtain second optimal theoretical data in the preset sampling time period comprises:
correcting the second theoretical data of each sampling moment in the preset sampling time period according to the second measured data of the sampling moment to obtain corrected data of the sampling moment;
and acquiring a correction data set in the preset sampling time period, and taking the minimum value in the correction data set as second optimal theoretical data in the preset sampling time period.
6. The method according to claim 1, wherein the expression of the predetermined wave velocity-apparent resistivity conversion model is:
wherein, the first and the second end of the pipe are connected with each other,
are all constraint parameters->
For the first theoretical data->
For the second theoretical data->
Is the depth.
7. An underground space anomaly detection device, comprising:
the wave impedance model optimal value output module is used for acquiring first theoretical data and first measured data which change along with wave speed in a preset sampling time period, and correcting the first theoretical data according to the first measured data to obtain first optimal theoretical data in the preset sampling time period; the first theoretical data are obtained according to an underground space wave impedance model established by wave impedance inversion, and the first measured data are high-frequency band signals and low-frequency band signals of an abnormal body acquired based on the improved dual-channel detector in the preset sampling time period;
the apparent resistivity model optimal value output module is used for acquiring second theoretical data and second measured data which change along with apparent resistivity in the preset sampling time period, and correcting the second theoretical data according to the second measured data to obtain second optimal theoretical data in the preset sampling time period; the second theoretical data are obtained through an apparent resistivity model of the underground space established through smoke ring inversion, and the second measured data are induced electromotive force obtained based on the improved array receiving coils in the preset sampling time period;
the abnormal body detection result acquisition module is used for substituting the first optimal theoretical data and the second optimal theoretical data into a preset wave velocity-apparent resistivity conversion model to obtain a detection result of the abnormal body;
wherein the wave impedance model is established by:
based on the cavity structure of the underground space, the dual-channel detector is used for collecting high-frequency band signals and low-frequency band signals of the abnormal body, the impedance information of the reflection interface wave of the surrounding rock body, the cavity modal response characteristic and the structure frequency characteristic of the abnormal body are calculated according to the time-frequency characteristics of the signals, and the corresponding underground space wave impedance model of the collected signals is built according to the wave impedance inversion.
8. A terminal device comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, the processor implementing the underground space anomaly detection method according to any one of claims 1 to 6 when executing the computer program.
9. A computer-readable storage medium, comprising a stored computer program, wherein the computer program, when executed, controls an apparatus in which the computer-readable storage medium is located to perform the method for detecting a subterranean space abnormal body according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210108195.6A CN114488346B (en) | 2022-01-28 | 2022-01-28 | Underground space abnormal body detection method, device, equipment and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210108195.6A CN114488346B (en) | 2022-01-28 | 2022-01-28 | Underground space abnormal body detection method, device, equipment and medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114488346A CN114488346A (en) | 2022-05-13 |
| CN114488346B true CN114488346B (en) | 2023-04-07 |
Family
ID=81477526
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210108195.6A Active CN114488346B (en) | 2022-01-28 | 2022-01-28 | Underground space abnormal body detection method, device, equipment and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114488346B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105068118A (en) * | 2015-06-10 | 2015-11-18 | 中国矿业大学 | Method for interpreting range of magmatic rock intrusion into coal layer based on seismic information |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107272068A (en) * | 2017-06-09 | 2017-10-20 | 北京市地质调查研究院 | It is a kind of to utilize the method without anti-interference across the hole resistivity method Underground space of cable |
| CN107516191B (en) * | 2017-08-10 | 2020-09-04 | 广州市城市规划勘测设计研究院 | Multi-rule-in-one and one-graph project collaborative approval system, method and tool |
| US10920585B2 (en) * | 2017-12-26 | 2021-02-16 | Saudi Arabian Oil Company | Determining sand-dune velocity variations |
| CN112285802A (en) * | 2020-10-20 | 2021-01-29 | 重庆璀陆探测技术有限公司 | Tunnel earthquake and transient electromagnetism combined detection method |
-
2022
- 2022-01-28 CN CN202210108195.6A patent/CN114488346B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105068118A (en) * | 2015-06-10 | 2015-11-18 | 中国矿业大学 | Method for interpreting range of magmatic rock intrusion into coal layer based on seismic information |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114488346A (en) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102662195B (en) | Underground cavity detection system | |
| CN103869371B (en) | Manual field source frequency domain full-gradient electromagnetic measuring method | |
| CN110531422B (en) | Tensor artificial source electromagnetic signal data acquisition processing method and device | |
| Longoni et al. | Surface and subsurface non-invasive investigations to improve the characterization of a fractured rock mass | |
| CN104375191A (en) | Magnetotelluric data acquisition system and method | |
| CN109387834B (en) | Dam damage detection method for underground reservoir and electronic equipment | |
| CN111708094B (en) | Sandstone-type uranium ore sand body identification method and system based on wide-area electromagnetic method | |
| CN104020496A (en) | Ground controlled source magnetotelluric method based on axial collinear manner | |
| CN112327372B (en) | Comprehensive detection method for internal hidden danger of long-distance linear engineering canal dike | |
| CN106371142A (en) | Coil for transient electromagnetic explorations and the application method thereof | |
| CN107085240A (en) | A kind of side slope magnetic fluid detection system and method | |
| CN116931101A (en) | Method and device for detecting underground shallow disease body, electronic equipment and storage medium | |
| CN102662191A (en) | Underground pipeline tomography system for underground space | |
| CN103376443B (en) | Ground penetrating radar terrestrial interference detecting and fast eliminating method | |
| CN112180444B (en) | Detection method and device for stratum velocity structure and storage medium | |
| CN114488346B (en) | Underground space abnormal body detection method, device, equipment and medium | |
| Sørensen et al. | TDEM in groundwater mapping—a continuous approach | |
| CN112433253A (en) | Three-dimensional detection method for surrounding rock tunnel with weak broken flowing water | |
| CN114791633B (en) | Method, system and medium for monitoring shale gas fracturing | |
| CN113419288B (en) | Data detection and pretreatment method for inversion of underground shelter | |
| JP5562237B2 (en) | Method and apparatus for monitoring storm activity on the surface of the earth in real time | |
| CN112904348B (en) | Three-dimensional detection method, device, equipment and storage medium | |
| CN116736394A (en) | Exploration method and system based on ocean electromagnetic difference | |
| Sun et al. | Research on detection and visualization of underground pipelines | |
| WO2021098667A1 (en) | Full-coverage anti-interference controlled-source electromagnetic detection method |
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 |