DRIFT TUBE BOREHOLE MUON DETECTOR SYSTEM, APPARATUS, AND METHOD FOR MUON RADIOGRAPHY AND TOMOGRAPHY

Abstract

A borehole muon detector for muon radiography or geotomography is provided, the borehole muon detector including a substantially cylindrical housing, which defines a bore, a pair of end caps, each end cap sealing an end of the cylindrical housing and a plurality of sealed drift tubes which are longitudinally disposed in the bore of the housing to form a bundle of drift tubes, wherein each sealed drift tube comprises: a centrally located anode wire disposed on a longitudinal axis; an inner surface which is coated with a cathode coating, the cathode coating divided into a first cathode pad and a second cathode pad by a Vernier pattern; and a timer in electrical communication with the anode wire for measuring a drift time. A system and a method are also provided.

Claims

1.-10. (cancelled)

11. A muon detector for use in a borehole, the muon detector comprising: a plurality of elongated drift tubes forming a bundle of drift tubes, the bundle of drift tubes being arranged to be discretely circularly symmetric about a longitudinal axis; and data-acquisition electronics connected to receive electronic signals from each of the plurality of drift tubes and configured to determine trajectories of muons traversing the muon detector based on at least one of the signals received from the drift tubes.

12. The muon detector of claim 11 wherein the data-acquisition electronics are configured to determine an azimuthal angle about the longitudinal axis of a muon traversing the muon detector based on drift times of multiple drift tubes in the bundle.

13. The muon detector of claim 12 wherein each drift tube comprises a centrally located anode wire disposed on a longitudinal axis of the drift tube and the data-acquisition electronics are configured to determine the azimuthal angle about the longitudinal axis of the muon traversing the muon detector by: determining, for each of the multiple drift tubes, an isochrone comprising a radial distance of the muon from the anode wire; and, fitting a line through the isochrones of at least three drift tubes from among the multiple drift tubes to thereby determine the azimuthal angle.

14. The muon detector of claim 11 wherein the data acquisition electronics are configured to determine a zenith angle relative to the longitudinal axis of a muon traversing the muon detector based on longitudinal position estimates from a sub-plurality of drift tubes in the bundle.

15. The muon detector of claim 14 wherein the data-acquisition electronics are configured to determine the zenith angle relative to the longitudinal axis of the muon traversing the muon detector by fitting a line through the longitudinal position estimates of the sub-plurality of drift tubes.

16. The muon detector of claim 14 wherein each drift tube comprises a centrally located anode wire disposed on a longitudinal axis of the drift tube and a bore-defining cathode surface and the bore-defining cathode surface of each drift tube is divided into a plurality (e.g. 2) of cathode parts to form a Vernier pattern.

17. The muon detector of claim 16 wherein the data acquisition electronics are configured to determine, for each drift tube in the sub-plurality of drift tubes, a coarse longitudinal position estimate of the muon traversing the drift tube based on at least one of: charge division where the anode wire is resistive and relative timing of electrical pulses at longitudinal ends of the cathode parts.

18. The muon detector of claim 17 wherein the coarse longitudinal position estimate has an uncertainty that is less than a spatial period of the Vernier pattern.

19. The muon detector of claim 17 wherein the data acquisition electronics are configured to determine, for each drift tube in the sub-plurality of drift tubes, a precise longitudinal position estimate of the muon traversing the drift tube based on a ratio of the voltages on the plurality of cathode parts.

20. A method for detecting muons in a borehole, the method comprising: arranging a plurality of elongated drift tubes in a bundle, the bundle of drift tubes discretely circularly symmetric about a longitudinal axis; placing a muon detector comprising the plurality of drift tubes into a borehole; determining trajectories of muons traversing the muon detector based on signals generated by at least one of the drift tubes.

21. The method of claim 20 comprising determining, for a muon traversing the muon detector, an azimuthal angle about the longitudinal axis of the muon traversing the muon detector based on drift times associated with multiple drift tubes in the bundle.

22. The method of claim 21 wherein determining the azimuthal angle about the longitudinal axis of the muon traversing the at least one muon detector comprises determining, for each of the multiple drift tubes, an isochrone comprising a radial distance of the muon from the anode wire; and, fitting a line through the isochrones of at least three drift tubes from among the multiple drift tubes to thereby determine the azimuthal angle.

23. The method of claim 20 comprising determining, for a muon traversing the muon detector, a zenith angle relative to the longitudinal axis of the muon traversing the muon detector based on longitudinal position estimates from a sub-plurality of drift tubes in the bundle.

24. The method of claim 23 wherein determining the zenith angle relative to the longitudinal axis of the muon traversing the muon detector comprises fitting a line through the longitudinal position estimates of the sub-plurality of drift tubes.

25. The method of claim 23 comprising shaping a bore-defining cathode surface of each drift tube into a plurality (e.g. 2) of cathode parts to form a Vernier pattern.

26. The method of claim 25 comprising determining, for each drift tube in the sub-plurality of drift tubes, a coarse longitudinal position estimate of the muon traversing the drift tube based on at least one of: charge division where an anode wire disposed along a central longitudinal axis of the drift tube is resistive and relative timing of electrical pulses at longitudinal ends of the cathode parts.

27. The method of claim 26 wherein the coarse longitudinal position estimate has an uncertainty that is less than a spatial period of the Vernier pattern.

28. The method of claim 27 comprising determining, for each drift tube in the sub-plurality of drift tubes in the bundle, a precise longitudinal position estimate of the muon traversing the drift tube based on a ratio of the voltages on the plurality of cathode parts.

29. A system for estimating subsurface densities, the system comprising: one or more muon detectors located in one or more boreholes, each muon detector comprising: a plurality of elongated drift tubes forming a bundle of drift tubes, the bundle of drift tubes being arranged to be discretely circularly symmetric about a longitudinal axis; and data-acquisition electronics connected to receive electronic signals from each of the plurality of drift tubes and configured to determine trajectories of muons traversing the muon detector based on the signals received from the drift tubes; and, a processing unit in communication with each of the one or more muon detectors and configured to generate a three-dimensional density profile of a region of interest based on a plurality of trajectories of muons traversing the one or more muon detectors.

30. The system of claim 29 wherein, for each muon traversing a muon detector of the one or more muon detectors, the data-acquisition electronics are configured to determine an azimuthal angle about the longitudinal axis of the muon traversing the muon detector based on drift times of multiple drift tubes in the bundle.

Description

FIGURES

[0029] FIG. 1 is a schematic of a borehole muon detector system of the present technology.

[0030] FIG. 2A is a cross sectional view of a borehole muon detector of the system of FIG. 1; and FIG. 2B is a longitudinal sectional view of a borehole muon detector of the system of FIG. 1.

[0031] FIG. 3 is a perspective longitudinal view of a drift tube of the borehole muon detector of FIGS. 2A and 2B.

[0032] FIG. 4 shows an exemplary Vernier pattern cathode system on the inner wall of a drift tube.

[0033] FIG. 5 is a schematic of the electronics of the borehole muon detector of FIGS. 2A and 2B.

[0034] FIG. 6A is a cross sectional view showing how a muon is tracked through a single drift tube of the borehole muon detector of FIG. 2A and 2B; FIG. 6B is a longitudinal view showing how a muon is tracked through a number of drift tubes of the borehole muon detector of FIG. 2A and 2B.

[0035] FIG. 7A is a cross sectional view showing how a muon is tracked through the borehole muon detector of FIGS. 2A and 2B; and FIG. 7B is a longitudinal view showing how a muon is tracked through the borehole muon detector of FIGS. 2A and 2B.

[0036] FIG. 8A and 8B shows sensor images for a simulated cuboid, 100 meters above a muon sensor situated underground. FIG. 8A shows the image with a sliding window filtering algorithm; FIG. 8B shows with image without the sliding window filtering algorithm; and FIGS. 8C and 8D show that each pixel represents a region of solid angle.

[0037] FIG. 9 shows an array of boreholes with borehole muon detectors installed and the detection of a density anomaly using an array of borehole muon detectors in the borehole muon detector system of FIG. 1.

DESCRIPTION

[0038] 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 including are not limiting. Further, the terms comprising, having, including, and containing are to be construed as open-ended terms (i.e., meaning including, but not limited to,) unless otherwise noted.

[0039] 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.

[0040] 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.

[0041] A borehole muon detector system, generally referred to as 8 is shown in FIG. 1. It includes one or more borehole muon detectors 10, a computer processing unit 12 and a user interface 14. Support cables 16 and electrical connections 18 extend from the borehole muon detector 10 to the computer processing unit 12. The borehole muon detector system 8 can be used in resource exploration and monitoring, and in security applications. The system measures the trajectory of muons that pass through it, from which a differential intensity measurement (radiographic image) can be derived. The differential intensity measurements from multiple borehole muon detection systems 8 can be inverted to estimate the three-dimensional distribution of subsurface density (tomographic analysis).

[0042] As shown in FIG. 2A, each borehole muon detector 10 includes drift tubes 20 which are arranged in a bundle, generally referred to as 22, with a plurality of scintillator members 24 around the periphery 26 of the bundle 22. The drift tubes 20 and the scintillator members 24 are encased in a housing 30, which is tubular. As shown in FIG. 2B, the drift tubes 20 and the scintillator members 24 are longitudinally disposed in the bore 32 of the housing 30 and lie overtop the drift tube bundle 22. It can be seen that they are segmented from one another. There are no orthogonally disposed drift tubes 20. The housing 30 has an endcap 34 at each end 36, 38. The endcaps 34 include apertures 40 for the electrical connections 16 and support cables 18. At one end 36 of the housing 30 there is a power input 42 and data output connector 44. The data output connector 44 consists of a ruggedized Ethernet or coaxial connector. The power input 42 is a suitable ruggedized DC electrical connector. At the other end 38 the housing 30 there is a power output 46 and data input connector 48. The data input connector 48 consists of a ruggedized Ethernet or coaxial connector. The power output 46 is a suitable ruggedized DC electrical connector. A support matrix 50 locates and retains the drift tubes 20.

[0043] As shown in FIG. 3, each drift tube 20 includes a tubular wall 56 of rigid material such as plastic, aluminum or carbon fiber, or a Mylar film that is appropriately glued or supported to maintain rigidity. The drift tube 20 has a first end 58 and a second end 60. Each cylindrical drift tube 20 is a few centimeters in diameter and many centimeters (up to a few meters) in length. The inner surface 62 of each drift tube wall 56 has an electrically conductive surface coating 64, which is a cathode.

[0044] The tubular wall 56 defines a bore 72. The bore 72 is filled with an ionizing gas mixture such as argon/carbon dioxide or other appropriate mixture. Centrally located in the bore 72 is an anode wire 74. The central anode wire 74 is an appropriate tungsten/gold/steel wire, affixed in place by plugs 76 at both ends of the tubular wall 60, and maintained at a high potential of a few thousand volts with respect to the cathode coating 64. The anode wire 74 has a first end 78 and a second end 80. It is disposed on the longitudinal axis 82. The drift tube 20 is sealed.

[0045] As shown in FIG. 4, the cathode coating 64 is divided into two separate parts, forming a Vernier pattern 66 which extends from the first end 36 to the second end 38, to provide a first cathode pad 68 and a second cathode pad 70. There may be more than two cathode pads 68, 70.

[0046] As shown in FIG. 5, each borehole muon detector 8 includes a printed circuit board (PCB) 90. The power input 42, data output connector 44, the power output 46 and the data input connector 48 are in electrical and electronic communication with the PCB 90.

[0047] The PCB 90 consists of preamplifiers 92, a high voltage generator 94 to supply potential to the drift tube anode wires 74, application specific integrated circuits (ASICs) 96, analog to digital converters (ADCs) 98 for measuring the voltage on the cathode pads 68, 70 of each drift tube 20, the current on the anode wires 74, and time to digital converters (TDCs) 100 for measuring the time of a pulses. Alternatively, charge division can be used with a resistive anode wire 74 instead of relative timing on the anode wires 74. The TDCs 100 determine both the drift time of the ionization and the relative timing along the anode wire 74 of the pulse arrival at each end 36, 38 of the drift tube 20. The latter is used to determine which of the repeated Vernier segments (a pair of cathode pads 68, 70) is aligned with the ionization event. The PCB 90 also includes a field programmable gate array (FPGA) 102, an accelerometer 104, a magnetometer 106, and a gyroscope 108. The TDCs 100 and ADCs 96 are connected to field programmable gate arrays (FPGAs) which process the data from all drift tubes 20 to identify coincidences consistent with the passage of a muon, and export this data to an online backend computer 18 for further processing, storage, and pipelining to a surface data collector 1 10. Additional electronics for monitoring temperature, gas pressure, and other slow controls are embedded within the system.

[0048] As shown in FIG. 6A, the electric field 120 is perpendicular to the longitudinal axis 82 and therefore the anode wire 72, as it traverses the bore 72 from the cathode coating 64 to the anode wire 72. The muons 122 create ionization 124 that drifts towards the anode wire 72 on drift paths 126. The radial position within the drift tube 20 (how close the muon passes by the anode wire 74) is determined by the time it takes for the closest electrons to drift to the anode wire 74. By measuring the difference in time or current between the current pulses measured on each end 78, 80 of the anode wire 74, a rough estimate of the position along the anode wire 74 where the muon passed by is achievable. Each drift tube 20 measures position along the anode wire 74 (Z measurement), and a drift time relative to a reference clock (may be provided by a scintillator pulse from the scintillator members 24).

[0049] As shown in FIG. 6B, the Vernier pattern 66 of the cathode coating 64 then allows for a much finer measurement of longitudinal position. Within a single Vernier pattern the fractional mirror charge measured on the respective cathodes is linearly related to longitudinal position. If a muon 122 passes through at least three drift tubes 20 within the housing 30, the azimuthal angle of the muon trajectory can be accurately determined from the radial measurements. The zenith angle of the muon trajectory 128 with respect to the longitudinal axis 82 of the borehole muon detector 10 is determined from regression on the longitudinal positions measured in each of the drift tubes 20. Thus, a borehole muon detector 10 comprised of an array of drift tubes 20 instrumented with a Vernier pattern 66 on the cathode coating 64 of the drift tubes 20, allows for precise determination of muons 122 that pass through the borehole muon detector 10. Difference in drift times between multiple drift tubes 20 in a bundle 22 provides X-Y measurements and with X-Y-Z measurements from multiple drift tubes 20, the full trajectory can be reconstructed.

[0050] FIG. 7A shows how a muon is tracked through the borehole muon detector 10 in the transverse plane. A line is fit through all the drift tubes 20 in which a pulse is observed, using the isochrones (white circles) to perform a best fit. FIG. 7B shows how a muon is tracked through the borehole muon detector 10 in the longitudinal plane. A line is fit using the longitudinal position estimates (white points) coming from time differences between first end 58 and second end 60 measurements of each drift tube 20 and from the voltage division among the Vernier pattern cathodes 64.

[0051] FIGS. 8A and 8B show two dimensional radiographic images of muon intensity underneath a simulated 50 m50 m20 m cuboid (1 g/cc contrast), 100 meters above a borehole muon detector that is situated 400 meters underground. The image is shown with (FIG. 8A) and without (FIG. 8B) a sliding window filtering algorithm applied. Each pixel represents a region of solid angle as shown in FIGS. 8C and 8D. The muon intensity shadow shows the location of a density anomaly above the borehole muon detector 10. Three dimensional images can be obtained from the two dimensional images by applying inversion algorithms to a set or plurality of two dimensional images.

[0052] FIG. 9 is a schematic of a borehole muon detector system 8 with multiple borehole muon detectors 10 arranged in series in each of two boreholes.

[0053] In an embodiment, the borehole muon detector system 8 or the borehole muon detector 10 may determine an occupancy or muon intensity map (radiographic image) of a projected surface above the borehole muon detector system 8 or the borehole muon detector 10 using the muon trajectory for muons that pass through the borehole muon detector system 8 or the borehole muon detector 10. The muon intensity in each section of solid angle emanating from the borehole muon detector system 8 or the borehole muon detector 10 to the surface is a measure of the average density of the earth within that volume. Regions of lower or higher density will have correspondingly higher or lower muon occupancy within the corresponding pixels of the radiographic images. The radiographic images from the borehole muon detector system 8 may facilitate tomographic analysis in order to estimate the three-dimensional distribution of subsurface density above the borehole systems.

[0054] The invention comprises a number of non-limiting aspects. Non-limiting aspects of the invention provide: [0055] 1. A borehole muon detector for muon radiography or geotomography, the borehole muon detector including a substantially cylindrical housing, which defines a bore, a pair of end caps, each end cap sealing an end of the cylindrical housing and a plurality of sealed drift tubes which are longitudinally disposed in the bore of the housing to form a bundle of drift tubes, wherein each sealed drift tube comprises: a centrally located anode wire disposed on a longitudinal axis; an inner surface which is coated with a cathode coating, the cathode coating divided into at least first cathode pad and a second cathode pad by a Vernier pattern; and a printed circuit board in electrical communication with the anode wire for measuring a current. [0056] 2. The borehole muon detector of aspect 1, wherein there are at least three drift tubes in the bundle. [0057] 3. The borehole muon detector of aspect 2, further including a plurality of scintillator members disposed on the bundle. [0058] 4. The borehole muon detector of aspect 2 or 3, further including at least one position sensor. [0059] 5. The borehole muon detector of aspect 4, wherein the position sensor is one or more of an accelerometer, a magnetometer and a gyroscope. [0060] 6. A system for borehole muon geotomography, the system including at least one borehole muon detector for placement in at least one boreholes and a processor in electronic communication with each of the borehole muon detectors, each borehole muon detector including a substantially cylindrical housing, which defines a bore, a pair of end caps, each end cap sealing an end of the cylindrical housing and a plurality of sealed drift tubes which are longitudinally disposed in the bore of the housing to form a bundle of drift tubes, wherein each sealed drift tube comprises: a centrally located anode wire disposed on a longitudinal axis; an inner surface which is coated with a cathode coating; and printed circuit board in electrical communication with the anode wire and in electronic communication with the processor. [0061] 7. The system of aspect 6, further comprising a plurality of scintillator members disposed on at least one bundle. [0062] 8. The system of aspect 6 or 7, wherein each borehole muon detector further includes at least one position sensor. [0063] 9. The system of aspect 8, wherein the position sensor is one or more of an accelerometer, a magnetometer and a gyroscope. [0064] 10. The system of any one of aspects 6 to 9, wherein the cathode coating is divided into at least a first cathode pad and a second cathode pad by a Vernier pattern. [0065] 11. A method for borehole muon radiography or geotomography, the method comprising: positioning a plurality of borehole muon detectors in one or more boreholes, each borehole muon detector including a plurality of sealed drift tubes, each drift tube including a centrally located anode wire disposed on a longitudinal axis, an inner surface which is coated with a cathode coating and timer; generating a radial electric field between the anode wire and the cathode coating; and measuring a plurality of muon trajectories. [0066] 12. The method of aspect 11, wherein the plurality of muon trajectories are processed to provide a two dimensional image. [0067] 13. The method of aspect 11, wherein the plurality of muon trajectories are processed to provide a plurality of two dimensional images. [0068] 14. The method of aspect 13, wherein the plurality of two-dimensional images are further processed by applying inversion algorithms to provide at least one three dimensional image.

[0069] 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.