MUON TOMOGRAPHY SYSTEM, APPARATUS, AND METHOD FOR TUNNEL DETECTION

Abstract

A system for passively monitoring territory proximate to or at restrictive boundaries for tunnels, the system comprising a plurality of muon sensors, a data network in communication with each muon sensor, a power network in electrical communication with each muon sensor, and a data analysis unit, the data analysis unit in communication with each muon sensor via the data network, the data analysis unit comprising a memory and a processor, the memory configured to instruct the processor to analyse data from the plurality of muon detectors to identify and locate a new or emerging tunnel. A method of locating tunnels is also provided.

Claims

1. A method for detecting elongated subterranean voids in a region of interest, the method comprising: locating a plurality of muon detectors in a plurality of spaced apart subterranean locations in, or in a vicinity of, the region of interest; detecting, by each muon detector, muon intensity data corresponding to a rate of muons that intersect the muon detector with trajectories in each of a plurality of solid angles; processing the muon intensity data from the plurality of muon detectors to identify one or more elongated subterranean voids in the region of interest.

2. The method of claim 1 wherein processing the muon intensity data comprises, for at least a first one of the plurality of muon detectors, identifying a first elongated feature comprising a first contiguous region of solid angles extending from the first one of the plurality of muon detectors and having muon intensities greater than threshold.

3. The method of claim 2 wherein processing the muon intensity data comprises, for at least a second one of the plurality of muon detectors, identifying a second elongated feature comprising a second contiguous region of solid angles extending from the second one of the plurality of muon detectors and having muon intensities greater than threshold.

4. The method of claim 3 wherein processing the muon intensity data comprises identifying an elongated subterranean void based on overlapping or connected regions of solid angle corresponding to the first and second elongated features.

5. The method of claim 4 wherein processing the muon intensity data comprises identifying an elongated subterranean void in three dimensions.

6. The method of claim 3 wherein processing the muon intensity data comprises identifying an elongated subterranean void based on shared fields of view of the first and second ones of the plurality of muon sensors.

7. The method of claim 1 wherein processing the muon intensity data comprises identifying an elongated subterranean void based on shared fields of view of a first muon sensor from among the plurality of muon sensors and a second muon sensor from among the plurality of muon sensors.

8. A method for monitoring a region of interest, the system comprising: locating a plurality of muon detectors in, or in a vicinity of, the region of interest, each muon detector configured to detect muons passing therethrough to generate muon intensity data; and, processing the muon intensity data from each of the plurality of muon detectors, wherein processing the muon intensity data comprises: detecting, based on the muon intensity data from at least some of the plurality of muon sensors, spatially connected regions in the muon intensity data to identify a feature of interest; and, calculating a statistical significance of the feature of interest.

9. The method of claim 8 wherein processing the muon intensity data from each of the plurality of muon detectors comprises determining whether the statistical significance exceeds a threshold level.

10. The method of claim 8 wherein processing the muon intensity data from each of the plurality of muon detectors comprises analyzing a first set of muon intensity data associated with a first muon sensor of the plurality of muon sensors and a second set of muon intensity data associated with a second muon sensor of the plurality of muon sensors to determine that the feature of interest is an elongated subterranean void.

11. The method of claim 8 wherein, in the event a common feature of interest is detected in muon intensity data of two or more muon sensors, processing the muon intensity data from each of the plurality of muon detectors comprises determine a three-dimensional region of interest based on shared fields of view of the two or more muon sensors.

12. The method of claim 8 wherein the muon detectors are configured to determine muon intensity data which comprises trajectories of the muons passing through each of the plurality of detectors.

13. The method of claim 12 wherein the muon intensity data comprises, for each of the plurality of muon detectors, a rate of muons that pass through the detector in each of a plurality of trajectories.

14. The method of claim 13 wherein, for each of the plurality of muon detectors, the plurality of trajectories comprise a plurality of solid angles originating at the detector.

15. The method of claim 12 wherein detecting spatially connected regions in the muon intensity data to identify the feature of interest is based at least in part on the trajectories.

16. The method of claim 15 wherein detecting spatially connected regions in the muon intensity data comprises identifying contiguous regions of solid angle.

17. A system comprising: a plurality of muon detectors located in a plurality of spaced apart subterranean locations in, or in a vicinity of, the region of interest; a processor configured to: detect, for each muon detector, muon intensity data corresponding to a rate of muons that intersect the muon detector with trajectories in each of a plurality of solid angles; process the muon intensity data from the plurality of muon detectors to identify one or more elongated subterranean voids in the region of interest.

Description

FIGURES

[0034] FIG. 1 is a prior art schematic of the relationship between a pixel in a radiographic image and a section of solid angle extending from the underground muon sensor to the surface of the earth.

[0035] FIG. 2 is a partial cutaway schematic of the device of the present technology.

[0036] FIG. 3 is a perspective view of the scintillator bars of FIG. 2.

[0037] FIG. 4 is a schematic longitudinal view through two adjacent scintillator bars of the muon sensor of FIG. 3.

[0038] FIGS. 5A and 5B are schematics showing a muon trajectory through a scintillator bar array 32. FIG. 5A shows the (x.sub.1,x.sub.3) projection, and FIG. 5B shows the (x.sub.2,x.sub.3) projection.

[0039] FIG. 6 is a schematic showing an installation of the present technology.

[0040] FIGS. 7A, B and C are schematics showing the identification of linear features of higher muon intensity in images and the relationship of these features to planar sections in three dimensions, and the method of locating a tunnel target in three dimensions by combining measurements from a series of muon sensors. FIG. 7A shows muon radiographs for two muon sensors; FIG. 7B shows muon radiographs for two other two muon sensors; and FIG. 7C shows the reconstruction of the radiographs to provide the location and size of a tunnel.

[0041] FIG. 8 is a block diagram of the method of the present technology.

DESCRIPTION

[0042] Except as otherwise expressly provided, the following rules of interpretation apply to this specification (written description and claims): (a) all words used herein shall be construed to be of such gender or number (singular or plural) as the circumstances require; (b) the singular terms a, an, and the, as used in the specification and the appended claims include plural references unless the context clearly dictates otherwise; (c) the antecedent term about applied to a recited range or value denotes an approximation within the deviation in the range or value known or expected in the art from the measurements method; (d) the words herein, hereby, hereof, hereto, hereinbefore, and hereinafter, and words of similar import, refer to this specification in its entirety and not to any particular paragraph, claim or other subdivision, unless otherwise specified; (e) descriptive headings are for convenience only and shall not control or affect the meaning or construction of any part of the specification; and (f) or and any are not exclusive and include and comprising are not limiting. Further, the terms comprising, having, comprising, and containing are to be construed as open-ended terms (i.e., meaning comprising, but not limited to,) unless otherwise noted.

[0043] Recitation of ranges of values herein are merely intended to serve as a shorthand method of referring individually to each separate value falling within the range, unless otherwise indicated herein, and each separate value is incorporated into the specification as if it were individually recited herein. Where a specific range of values is provided, it is understood that each intervening value, to the tenth of the unit of the lower limit unless the context clearly dictates otherwise, between the upper and lower limit of that range and any other stated or intervening value in that stated range, is included therein. All smaller sub ranges are also included. The upper and lower limits of these smaller ranges are also included therein, subject to any specifically excluded limit in the stated range.

[0044] Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the relevant art. Although any methods and materials similar or equivalent to those described herein can also be used, the acceptable methods and materials are now described.

Definitions:

[0045] Muon sensorin the context of the present technology, a muon sensor is any one of: one or more drift tubes longitudinally disposed within a housing, each tube consisting of a bore filled with an ionizing gas containing an axial anode wire; or one or more drift chambers each consisting of a volume of ionizing gas with a plurality of anode wires, the drift chambers either in a cylindrical form factor or in a plane transverse to the enclosure; or a micropattern ionizing gas detector in a cylindrical form factor or in the transverse plane to the enclosure; or a resistive plate chamber in a planar form factor transverse to the enclosure; or segmented plastic scintillators (bars) in the transverse plane to the enclosure, disposed in two orthogonal directions in multiple planes and each bar coupled to photodetectors at one or both ends; or elongated scintillators arranged in a helical bundle and coupled to photodetectors at one or both ends.

[0046] Underground elongated voidin the context of the present technology, an underground elongated void is an underground passageway, tunnel or channel. An elongated void may be substantially linear. It is identified using data from two or more muon detectors that detect the same elongated feature and is confirmed through analytics.

[0047] Restrictive boundariesin the context of the present technology, a restrictive boundary is a wall, a fence, a gate, a political border and the like for stopping free passage of people, goods, animals and the like.

[0048] Elongated featurein the context of the present technology, an elongated feature is a putative or potential elongated void that is detected by at least one muon sensor and is yet to be processed and analysed in the determination of it being an elongated void.

[0049] Seriesin the context of the present technology, a series is three or more, in a substantially linear arrangement. A series may be, for example, between 3 and 10 sensors, between 10 and 100 sensors or greater than 100 sensors.

DETAILED DESCRIPTION

[0050] The present technology relates generally to systems, apparatus, and methods of tracking cosmic ray muons through an underground sensor in order to develop an image of subsurface density above the sensor (muon radiography), to identify features in a radiographic image consistent with the presence of an underground elongated void in a region above and beside the sensor, and to use a series of sensors to locate the underground elongated void in three dimensions. This is useful for searching for unknown and potentially hazardous tunnels. Muon tracking sensors can reconstruct the trajectory (with some finite precision) of almost all muons that pass through them within a large acceptance window, typically around +60 from vertical with respect to the sensor for a variety of sensor configurations. By recording the rate of muons that pass through the sensor from all angles within the acceptance window, a muon intensity image (radiograph) can be determined as shown in FIG. 1. This intensity image is a measure of the average density of the ground above the sensor.

[0051] By identifying lines of higher than expected muon intensity in muon radiographic images, an underground tunnel may be identified or targeted for further investigation. For a single muon sensor, such a line corresponds to a planar section of contiguous solid angle emanating from the sensor to the surface.

[0052] As shown in FIG. 2, the muon sensor, generally referred to as 10 has an enclosure 20, which is a rectangular shell with a first endcap 22 containing feedthroughs 24 for electrical connections and support cables and a second endcap 26. In an embodiment, the enclosure may be cylindrically shaped. The enclosure 20 is made of suitable plastic, steel, aluminum, or other metal alloy. The enclosure interior 28 has support structures 30 to hold a scintillator bar array 32 firmly in place in a fixed orientation, and it also has a shelf 34 to hold the electronic components (accelerometers 36, gyroscopes 38, humidity sensor 40, temperature sensor 42, printed circuit boards (PCB) 44, field programable gate array (FPGA), 46, magnetometer 48, processor 50 and memory 52, which are referred to collectively as the digital data acquisition system (DAQ) 60). The accelerometers 36 and gyroscopes 38 are for identifying the exact position and orientation once installed, and the humidity sensor 40 and temperature sensor 42 are for monitoring the operational state of the optical sensors and electronic data acquisition system. In the embodiment shown there are three scintillator bar arrays 32. An exemplary photodetector 62 is shown. The power input feedthrough is a suitable ruggedized DC electrical connector 64. The data input connector consists of a ruggedized shielded data cable connector 66.

[0053] As shown in FIG. 3, in a scintillator bar array 32, the scintillator bars 70 are arranged in a first plane 72 of scintillator bars 70 that are packed in a pattern that minimizes any space between the scintillator bars 70 and a second plane 74 in which the scintillator bars 70 are orthogonal to the scintillator bars 70 of the first plane 72 and are packed in a pattern that minimizes any space between the scintillator bars. The scintillator bars 70 have less than a few square centimeters in cross-sectional area, and many centimeters in length. The scintillator bars 70 have a central bore 76 through which a wave-length shifting (WLS) optical fiber 78 extends. The fibers 78 have a high internal reflection coefficient to transmit photons collected from the scintillation pulse to the fiber ends. The central bore 76 has a cross sectional area of less than a few square millimeters.

[0054] In an alternative embodiment, the scintillator bar array 32 within each muon sensor consists of multiple pairs of orthogonally oriented planes 72, 74 of long scintillator members.

[0055] As shown in FIG. 4, photodetectors 62 are affixed to each end 90, 92 (or, in an alternative embodiment, one end) of the optical fiber 78 for each scintillator bar 70. The photodetectors 62 are connected to the PCB 44 through cables or soldered mounts. The PCB 44 consists of preamplifiers to amplify the current output from the photodetectors 62, a high voltage generator to supply potential to the photodetectors 62, application specific integrated circuits (ASIC's), analog to digital converters, a field programmable gate array (FPGA), an accelerometer, magnetometer, and gyroscope, and ethernet and/or coaxial connectors for digital data and power.

[0056] In an alternative embodiment, a clear optical fiber replaces the wave-length shifting optical fiber 78, if the wavelength of the emitted light from the scintillator bars 70 matches the optimal sensitivity of the photodetector. In another embodiment of the device, the scintillator bars 70 are very narrow and do not have a central bore 76 and no optical fiber 78. The photodetectors 62 are attached to each end of the scintillator bars 70.

[0057] In another embodiment, there are multiple scintillator bar arrays 32 within each muon sensor 10, each array 32 consisting of multiple pairs of orthogonally oriented planes of scintillator bars 70 with either a rectangular or triangular, or hexagonal cross-sectional shape.

[0058] In all embodiments, the scintillator bars 70 are organic polymer scintillators containing an appropriate dopant so that scintillation light is produced by the passage of charged particles. Each scintillator bar 70 is coated with a reflective layer on the sides and ends so that light produced inside the scintillator bar 70 is reflected back into it.

[0059] In all embodiments, the photodetectors 62 are multi-pixel photodiode arrays enabled by a bias voltage of a few tens of volts, supplied by a voltage regulator within the system. The electrical signals from the photodetectors 62 are amplified and passed to a PCB 44 within which is embedded an integrated electronic circuit that precisely times the arrival of the electrical signal. The integrated circuit propagates the timing and light yield information for all channels to a device such as an FPGA 46 for coincidence detection, when a user-configured number of photodetectors emit an electrical signal (detect light) within a user-configured period of time. The positions of the scintillator bars 70 from which light was detected are used to reconstruct the path of the muon through the device 10.

[0060] As shown in FIG. 5A and B, using a single muon as an example, the muon travels in straight-line trajectory 94 through a single pair of orthogonal planes 72, 74 of scintillator bars 70. FIG. 5A shows the (x.sub.1,x.sub.3) projection, and FIG. 5B shows the (x.sub.2,x.sub.3) projection, of the same plane pair. As the muon passes through each pair of orthogonal planes 72, 74, one or more scintillator bars 70 in each of the orthogonal planes emits scintillation light, which is captured by the WLS fiber 78 and transmitted to the photodetector 62, creating an electrical signal that is analyzed with the digital logic of the DAQ 60. Since the plane pairs consist of orthogonally oriented arrays 32 of scintillator bars 70, sets of (x.sub.1,x.sub.3) and (x.sub.2,x.sub.3) measurements of the muon position are performed by using the knowledge of the position of the scintillator bars 70, where x.sub.1, x.sub.2, and x.sub.3 are the three coordinates required to identify a point in three dimensions. Multiple pairs of orthogonal planes separated by a distance in x.sub.3 are used to determine the position of the muon at different points along its trajectory 94. These multiple position measurements are exploited to determine a straight-line muon trajectory 94 in three dimensions for each muon that passes through the sensor. The rate at which muons by a single muon sensor 10 are detected within certain ranges of angle (pixels) is a muon radiographic image. From one plane pair, a three-dimensional muon position measurement is ascertained. Criteria that distinguish muon trajectories from background noise (such as arising from latent radioactivity in the surrounding environment, or from electronic noise) are applied to the reconstructed tracks, and the path trajectory 94 is saved to the sensor radiograph data for further offline analysis. Data are sent to a data aggregation network that collects data from the data interfaces of the series of sensors.

[0061] As shown in FIG. 6, a muon sensor system or installation, generally referred to as 98, includes a series of muon detectors, generally referred to as 100, is installed underground below a region of interest in boreholes 204. A power distribution network 102 provides electrical power through cabling to the enclosure feedthroughs of each muon sensor 10, and the muon sensors 10 are interconnected on a data network 104 over which radiograph data for each muon sensor 10 can be queried. The series 100 of muon sensors 10 is suitably arranged to maximize three-dimensional reconstruction capability of any observed linear anomalies to target the location of tunnels and voids. The muon intensity measurements (radiographs) from the muon sensors 10 in the series 100 are collected over the data network 104 and sent to data aggregator 108 for further offline analysis and interpretation. Multiple trajectories are shown passing through the tunnel 206 and intersecting one muon sensor. For three dimensional images other multiple trajectories pass through the tunnel 206 and intersect at least one other muon sensor 10.

[0062] Over any period of exposure time the radiographic images from each muon sensor 10 in the series 100 can be analyzed offline using a specialized algorithm. Elongated voids of higher muon intensity are identified. Corresponding sections of solid angle are constructed and the combination thereof through specialized geophysical inversion algorithms or 3-dimensional surface intersection algorithms are used to identify the target location and extent in three dimensions.

[0063] As shown in FIGS. 7A, B and C by combining the images shown in FIG. 7A from a first series 100 of muon sensors 10 with the images shown in FIG. 7B from a second series 100 of muon sensors 10, the intersection of these planar sections for each muon sensor 10 can identify the location of candidate linear features, for example tunnels, in three dimensions. Elongated features (not necessarily aligned with a coordinate) can be identified in the images of FIGS. 7A and 7B. These features correspond to planar-like sections 200, 202 of solid angle emanating from each muon sensor 10 location as depicted in FIG. 7C. The boreholes 204 in which the muon sensors 10 are deployed is indicated by vertical lines, and the target tunnel location 206 is indicated by the horizontal line. The intersection 208 between the two planar sections 200, 202 shows the reconstructed position of the tunnel-like anomaly 206 identified by both of the muon sensors 10.

[0064] In the method, as shown in FIG. 8, a series of muon sensors is installed 400 in boreholes, other underground facilities, or adjacent to topographic features such as hills or mountains for monitoring for pre-existing and emerging tunnels and passageways in territory proximate or at political borders and other restrictive boundaries. Although the drawings show two or three muon sensors, it is envisioned that there will be tens to thousands of muon sensors in a series. The muon sensors each continuously and passively monitor 402 the incoming cosmic ray muon intensity from all angles within each sensor's field of view. The field of view of a sensor is an inverted cone or pyramidal shape (depending on the coordinate system chosen) emanating from the sensor. Light is emitted 404 from the scintillator bars in response to the incoming muons and is detected 406 by the photodetector, which then sends 408 a signal to the DAQ. The DAQ analyzes 410 the signals and determines 412 the trajectories of the incident cosmic ray muons.

[0065] The trajectories of muons are described in a straight-line sense by two discrete three-dimensional points along the trajectory to a fixed point in the sensor or equivalently by a single discrete three-dimensional point and two angles. The DAQ creates 414 a dataset of radiographic images defined in some coordinate system (such as polar angles q, f or rectilinear coordinates such as tan(0x), tan(0y)) which is then time-stamped 416 to provide 418 time-stamped intensity maps. The acquired muon flux data within the field of view is uploaded 420 over a common network at periodic intervals to the central data analysis unit. Alternatively, the raw data are time-stamped 422 and sent 424 to a computing device in the data analysis unit where the raw, time-stamped data are stored 426 for later analysis.

[0066] The computing device includes a processor and a memory, the memory comprising instructions thereon for instructing the processor. The processor, under control of the memory detects 428 spatially connected pixels in the muon intensity map which have either higher or lower (than the surrounding pixels) cosmic ray muon intensity, and identifies 430 where the connectedness is linear in the sense of being significantly elongated along a major axis that may be parallel to or at some angle with respect to either axis of the radiographic image. The processor, under control of the memory, identifies 432 elongated features on this basis.

[0067] The statistical significance of a given elongated feature is calculated 434 by the processor, under control of the memory. This can be applied to each muon sensor intensity map separately, or to a combined intensity map from multiple muon sensors. If data from a single sensor is analyzed, it will only provide two-dimensional information and will not discriminate between a surface anomaly and a subsurface anomaly, for example. Thus it cannot determine if an elongated feature is an elongated void, nor can it determine its position in three dimensions.

[0068] Either an operator or the processor under control of the memory sets 436 the alpha level to determine whether an elongated feature is significant or not. The alpha level takes into account the statistical fluctuations within the cosmic ray muon data as well as any inherent systematic errors of the sensor or data acquisition thereof. The processor under control of the memory reports 438 the significant elongated features and the memory stores 440 the processed data. Each elongated feature in the radiographic image derived from any muon sensor corresponds to an approximately planar section of solid angle emanating from the sensor.

[0069] In the case where an elongated feature is detected in two or more muon sensors, further processing may consist of identifying 442 three-dimensional regions of interest within shared fields of view of multiple muon sensors wherein there are overlapping or connected sections of solid angle corresponding to identified elongated features from multiple different radiographic images. The statistical significance of the three-dimensional regions of interest is calculated 444 by the processor, under control of the memory. Either an operator or the processor under control of the memory sets 446 the alpha level to determine whether an elongated feature is significant or not. The statistically significant three-dimensional elongated features are tagged 448 as elongated voids. The memory is configured to instruct the processor to utilize a feature extraction technique to analyse the images, utilize computer vision and conduct digital image processing in order to identify the elongated voids.

[0070] The processor optionally compares 450 tagged elongated voids over time to allow for four-dimensional interrogation (three spatial dimensions and time). Data from geographic locating devices comprising, but not limited to a Global Navigation Satellite System (GNSS), accelerometers and gyroscopes is sent 452 to the processor, which, under control of the memory, provides 454 a geographic location for each tagged elongated void. The processor under control of the memory then generates 456 a map of the area and the elongated void or voids which is stored 458 in the memory and is displayed 460 on a user interface.

[0071] With this method, an underground elongated void such as a tunnel, passageway, tunnel like void, tunnel-like region of reduced or enhanced density can be established for further targeted investigation, such as with underground drilling, seismic interrogation, or other techniques. The system may be an installation at or proximate a restrictive boundary such as a political border. The installation is a series of boreholes, each borehole retaining at least one muon sensor, the power and data networks in communication with each muon detector and the data analysis unit which is also in communication with the power and data networks.

[0072] While example embodiments have been described in connection with what is presently considered to be an example of a possible most practical and/or suitable embodiment, it is to be understood that the descriptions are not to be limited to the disclosed embodiments, but on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the example embodiment. Those skilled in the art will recognize or be able to ascertain using no more than routine experimentation, many equivalents to the specific example embodiments specifically described herein. Such equivalents are intended to be encompassed in the scope of the claims, if appended hereto or subsequently filed.

[0073] The invention has a number of non-limiting aspects. Non-limiting aspects of the invention comprise: [0074] 1. A system for passively monitoring territory proximate to or at restrictive boundaries, the system comprising a plurality of muon sensors, a data network in communication with each muon sensor, a power network in electrical communication with each muon sensor, and a data analysis unit, the data analysis unit in communication with each muon sensor via the data network, the data analysis unit comprising a memory and a processor, the memory configured to instruct the processor to analyse data from the plurality of muon detectors to identify and locate an elongated feature. [0075] 2. The system of aspect 1, wherein the memory is further configured to instruct the processor to analyze a data set from at least two muon sensors to determine if the elongated feature is a tunnel. [0076] 3. The system of aspect 2, wherein the memory is further configured to instruct the processor to create a three-dimensional map of the tunnel. [0077] 4. The system of aspect 3, wherein the memory is further configured to instruct the processor to identify a location of the tunnel. [0078] 5. The system of aspect 4, wherein each muon sensor includes at least two arrays, each array comprising elongated scintillator bars in a first plane and elongated scintillator bars in a second plane, the elongated scintillator bars in the first plane orthogonal to the elongated scintillator bars in the second plane. [0079] 6. The system of aspect 5, wherein each scintillator bar is individually coupled to a photodetector at a first end. [0080] 7. The system of aspect 6 wherein each scintillator bar is individually coupled to a photodetector at a second end. [0081] 8. The system of aspect 7, wherein each scintillator bar includes a bore and an optical fiber housed in the bore, the optical fiber comprising a first end and a second end, each of which are individually attached to the photodetector at each of the first end and the second end of the scintillator bar. [0082] 9. The system of aspect 8, wherein the muon detector includes a housing which houses the arrays and a plurality of electronic components. [0083] 10. The system of aspect 9, wherein the plurality of muon sensors is a series of muon sensors. [0084] 11. A method for monitoring a region of interest for elongated voids, the method comprising: [0085] providing a muon detection system comprising a plurality of muon sensors and a data analysis unit in communication with each muon sensor, the data analysis unit comprising a processor and a memory, the memory configured to instruct the processor; [0086] inserting each muon sensor of the plurality of muon sensors in a borehole; [0087] the plurality of muon sensors sensing muons; [0088] at least two muon sensors sending at least a first data set and a second data set to the processor, which under control of the memory: [0089] 1) identifies an elongated feature in both the first data set and the second data set; [0090] 2) calculates the statistical significance of the elongated feature; [0091] 3) analyzes the first data set and the second data set to determine whether or not the elongated feature is an elongated void; [0092] 4) calculates the statistical significance of the elongated void; [0093] 5) and creates a three-dimensional map of the elongated void. [0094] 12. The method of aspect 11, wherein the plurality of muon sensors is a series of muon sensors. [0095] 13. The method of aspect 12, further comprising the processor receiving a geographic position of the at least two muon sensors. [0096] 14. The method of aspect 13, further comprising mapping the geographic position of the elongated void. [0097] 15. The method of any one of aspects 1 1 to 14, wherein the region of interest is a restrictive boundary. [0098] 16. The method of aspect 15, wherein the restrictive boundary is a political border. [0099] 17. The method of any one of aspects 1 1 to 16, wherein the elongated void is a tunnel. [0100] 18. A method of reducing or eliminating tunneling through a restrictive boundary, the method comprising: [0101] providing a muon detection system comprising a plurality of muon sensors and a data analysis unit in communication with each muon sensor, the data analysis unit comprising a processor and a memory, the memory configured to instruct the processor; [0102] inserting each muon sensor of the plurality of muon sensors in a borehole along or proximate the restrictive boundary to provide a series of muon detectors; [0103] the series of muon sensors sensing muons; [0104] at least two muon sensors sending at least a first data set and a second data set to the processor, which under control of the memory: [0105] 1) identifies an elongated feature in both the first data set and the second data set; [0106] 2) calculates the statistical significance of the elongated feature; [0107] 3) analyzes the first data set and the second data set to determine whether or not the elongated feature is a tunnel; [0108] 4) calculates the statistical significance of the tunnel; [0109] 5) creates a three-dimensional map of the tunnel; [0110] 6) receives a geographic position of the at least two muon sensors; [0111] 7) maps the geographic position of the tunnel; and [0112] 8) reports the geographic position to an operator.