μ子能损超前探测技术研究

席振铢，刘骏华，何航，王亮，彭星亮，周胜, ，王鹤



μ子能损超前探测技术研究

席振铢1, 2，刘骏华1，何航1，王亮3，彭星亮1，周胜1, 3，王鹤1, 2

(1. 中南大学 地球科学与信息物理学院，湖南 长沙，410083；2. 教育部有色金属成矿预测重点实验室，湖南 长沙，410083；3. 湖南五维地质科技有限公司，湖南 长沙，410205)

依据宇宙射线μ子穿过高密度体物质时能量损失比穿过低密度体物质时的大的物理特性，推导电离损失、韧致辐射、电子偶衍生以及核相互作用的能损数学表达式，运用Geant4软件分别模拟μ子探测充填空气和水的溶洞穿透过程，记录其能量，构建μ子入射与出射能损成像，并与其散射成像图进行对比。研究结果表明：μ子能损成像有效解决了μ子散射对相近原子序数的元素组成物质分辨率低的问题，有望利用μ子低能损成像实现对充气和水等低密度空洞的超前探测，为超前探测提供一种新的探测方法。

μ子；能量损失；超前探测；溶洞

地质超前探测技术是一个国际前沿研究课题，人们对有关隧道地质超前预报技术研究应用已有40多年，但是预报的准确性和可靠性与隧道工程施工要求还有一定差距。早期地质超前预报主要是采用超前导洞、超前钻探方法勘探掌子面前方的地质情况，由于其效率较低、成本较高，因而人们逐步对无损地球物理超前探测技术进行研究[1]。目前，无损隧道超前预报技术主要有地震反射、地质雷达、直流电阻率、瞬变电磁、激发极化以及核磁共振等方法[2]，由于这些方法本身的分辨率较低、探测距离较短以及受隧道施工有限场地等因素制约，无损地球物理超前探测技术预报的准确率和探测距离亟需进一步提高。近20年来，宇宙μ子探测技术取得突破性进展。2003—2006年，LANL美国国家实验室科学家先后提出了最快径迹成像算法(point of closest approach, PoCA)、最大似然散射重建算法(the maximum likelihood scattering, MLS)以及最大似然散射与位移重建算法(the maximum likelihood scattering and displacement, MLSD)，实现了利用μ子对物体的散射成像计算[3−4]，奠定了μ子物理模拟基础。刘圆圆等[5−7]致力于μ子在高原子序数物质上的成像效果研究，分别对MLS和 MLSD 算法进行改进，提出了MLS-OSEM 和MLSD-OSEM 算法，提高了成像计算速度。LESPARRE等[8−9]针对标准岩石厚度与μ子穿透能量的理论计算及运用Soudan II探测器在地下700 m成功对月亮成像，证明μ子可以穿透地下几千米岩石层；DARIJANI等[10]应用μ子技术探测伊朗Pabdana煤矿厚度；SARACINO等[11]应用μ子技术探测那不勒斯Mt. Echia地下洞穴；MORISHIMA等[12]利用μ子对胡夫金字塔进行成像，利用宇宙射线μ子成像技术发现了埃及吉萨最大的金字塔——胡夫金字塔内的1个隐藏结构。可见，随着μ子成像技术以及探测器的进步，应用μ子成像技术开展超前地质预报成为可能。

1 能损成像原理

当高能量的初级宇宙射线进入大气层时，会与大气中的原子核连续发生作用，产生大量次级粒子。在这些次级粒子中，大部分次级粒子为π介子，分别带有3种不同的电性：带正电的π+介子、带负电的π−介子和中性介子π0，这3种粒子的数目都近乎相等。带正电的π+介子很快会衰变为μ子和μ子中微子，而带负电的π−将衰变成μ子和反μ子中微子，中性介子π0则很快会衰变为光子[13]。

通过测量μ子在通过物体时前后的角度偏差，可以确定物质的原子序数。散射角度偏差越小，则散射程度越小，物体的原子序数也越小，这是μ子在基于散射角偏差上成像的物理基础[17−18]。但对于原子序数差别不大的物质探测，散射成像相对困难，而μ子能损成像具有较大优势。

当μ子穿透物质时，其能量损失主要表现为电离损失(I)、韧致辐射(B)、电子偶的产生(P)以及核相互作用(N)这几种形式。设为面密度，其表达式为

设为μ子入射能量，则能量损失可以写为

式中：=lg。除少数情况外，几乎所有岩石的原子序数和质量比例都相同，区别在于各自的密度不同。密度和能量关系为

当μ子穿过密度为1 g/cm2的物质时，能量损失约为2.2 MeV[19]，因此，可以观测μ子穿过不同密度的围岩和溶洞时能量的衰减程度圈出溶洞，实现对低密度地质灾害的超前探测。

将μ子探测器置于隧道、坑道以及平硐等地下空间时，通过不同的阵列布置，接收来自于各方向上的μ子，记录其能量分布，得到μ子在各方向穿透地质体时的能量损失Δ，如图1所示。

2 岩溶仿真模拟

由于Geant4软件对模型边界条件有一定要求，本次模拟围岩的边界是有限的。为了计算方便，围岩设计长×宽×高为16 cm×16 cm×4 cm的长方体，探测目标球体溶洞的半径为2 cm，球心与长方体的中心重合。其中长方体背景围岩设定为碳酸盐岩，主要化学成分为CaCO3，其元素组成(质量分数)为Ca 40%，C 12%，O 48%，密度为2.9 g/cm3；在溶洞球体内部分别填充空气和淡水，模拟μ子能损探测溶洞过程。

2.1 充填空气岩溶模拟实验

地下岩溶极易形成空洞，隧道施工容易造成地面塌方，对人员和设备安全具有重大危害。在一般情况下，空洞充填空气主要物质成分是氮气和氧气，且氮气的体积分数约为78%，氧气的体积分数约为21%，密度为1.29 mg/cm3。为了计算简便，假设空洞为球体，围岩为长方体，球心与正方体几何对称中心重合，并建立图2所示三维直角坐标系。以模型中心为坐标原点，粒子入射方向向量为(0，0，1)，平行轴入射，发射200束能量为4 GeV的μ子，穿透长方体灰岩围岩和球状溶洞，其中入射和出射探测器位于=−2和=2的16 cm×16 cm长方体表面上，分别记录μ子入射和出射能量。

图1 μ子探测溶洞原理示意图

图2 溶洞填充空气仿真试验示意图

如图2所示，当μ子穿透模型时，由于电离损失、韧致辐射、电子偶衍生以及核相互作用，射线方向发生细微的偏转并产生大量电子，造成能量损失。

2.2 充水岩溶模拟实验

若岩溶充水，则隧道施工时容易造成透水事故。超前探测岩溶充水是地下工程施工预报地质灾害最重要的工作之一。一般来说，溶洞充填淡水，其结构式为H2O，密度为1 g/cm3。围岩参数与上述实验参数相同，同样发射200束能量为4.000 GeV的μ子，观察其能量变化，所得结果见图4。

从图4 (a)和图4 (b)可见：尽管μ子在入射面能量分布随机，没有规律性，但出射能量在模型中心部位幅值较大，而且等值线稀疏均匀，能量差异较小，其他位置能量幅值减弱，等值线变化较大，说明能量差异较大。通过μ子在穿越背景场和目标体时的能量损失Δ的位置分布(图4(c))可见：粒子在穿过模型时，模型中心位置呈现出1个能量低于30.000 MeV的低能异常，跨度范围在−2＜＜3 cm，−2＜＜2 cm之间，异常边界形态虽与实际溶洞形态相比存在失真现象，但异常位置在设计模型中心，范围与正演模型比对相差不大，反映了充水岩溶低密度的特征。其他位置的能量损失大都大于0.030 GeV，是相对高密度灰岩的反映。

(a) 入射能量；(b) 出射能量；(c) 亏损能量

(a) 入射能量；(b) 出射能量；(c) 亏损能量

3 能损与偏角成像比对

为了验证能损成像效果，采用上述能损仿真模型和仿真过程，分别模拟充气和充水μ子偏角成像，模拟结果如图5所示。

从图5(a)可以看出：μ子在穿越空气溶洞时偏转角度均未超过4 mrad，溶洞所在位置出现低异常，但异常形态呈现出不规则椭圆形状，与实际模型比对发现存在失真现象。从图5(b)可见：μ子在穿越充水溶洞时，异常的偏转角度为1~5 mrad，由于其物质组成元素皆为低原子序数元素，因此，溶洞弧度偏转角与灰岩围岩的偏转角度相差较小，成像边界不够清晰。

为了突出异常，设定μ子能损20 MeV和偏角 3 mrad为异常边界，构建岩溶μ子能损与偏角成像异常比对分布图，如图6所示。当溶洞内部填充物质为空气时，能损和偏角成像二者均能准确反映异常所在位置，但能损成像反映异常边界和大小的效果明显优于弧度偏转成像；而当溶洞内部填充物质为水时，由于水与碳酸盐岩的组成元素原子序数相近，从弧度偏转图中无法判定出异常位置，而能损成像效果则更明显，从而有效解决了μ子散射对相近原子序数元素组成物质分辨率低的问题。

(a) 充气岩溶μ子偏转；(b) 充水岩溶μ子偏转

(a) 充气岩溶μ子偏转；(b) 充水岩溶μ子偏转

4 结论

1)低密度的空气和水充填的岩溶相对高密度的灰岩围岩μ子能量亏损较小，利用相对围岩μ子低能损成像异常反映超前预报前方的溶洞。

2)μ子能损成像有效解决了μ子散射对相近原子序数元素组成物质分辨率低的问题。

3) 采用μ子能损可探测不同介质的充填物，所反映的异常形态不尽相同，充填水介质相对空气介质岩溶形态失真较大。

4) 利用现有隧道、巷道以及平硐等地下空间布置探测器，有望实现μ子的地下超前探测，为地质超前探测提供一种新的探测方法。

[1] 赵永贵, 刘浩, 孙宇, 等. 隧道地质超前预报研究进展[J]. 地球物理学进展, 2003, 18(3): 460−464. ZHAO Yonggui, LIU Hao, SUN Yu, et al. Advances in geological advance prediction of tunnels[J]. Advances in Geophysics, 2003, 18(3): 460−464.

[2] 李术才, 刘斌, 孙怀凤, 等. 隧道施工超前地质预报研究现状及发展趋势[J]. 岩石力学与工程学报, 2014, 33(6): 1090−1113. LI Shucai, LIU Bin, SUN Huaifeng, et al. Research status and development trend of geological prediction in tunnel construction[J]. Journal of Rock Mechanics and Engineering, 2014, 33(6): 1090−1113.

[3] BOROZDIN K N, HOGAN G E, MORRIS C, et al. Surveillance: radiographic imaging with cosmic-ray muons[J]. Nature, 2003, 422(6929): 277−278.

[4] SCHULTZ L, BLANPIED G, BOROZDIN K, et al. ML/EM reconstruction algorithm for cosmic ray muon tomography[C]// Nuclear Science Symposium Conference. Piscataway, NJ, USA: IEEE, 2006: 2574−2577.

[5] 刘圆圆, 赵自然, 陈志强, 等. 用于宇宙射线μ子成像的MLS-EM重建算法加速研究[J]. CT 理论与应用研究, 2007, 16(3): 1−5. LIU Yuanyuan, ZHAO Ziran, CHEN Zhiqiang, et al. Research on accelerating MLS-EM reconstruction algorithm for cosmic ray muon radiography[J]. CT Theory and Applications, 2007, 16(3): 1−5.

[6] 叶瑾. 利用宇宙线对高物质进行成像的技术研究[D]. 北京: 清华大学工程物理系, 2009: 126−129. YE Jin. Radiography of highmaterial with cosmic ray muon[D]. Beijing: Tsinghua University. Department of Engineering Physics, 2009: 126−129.

[7] 罗志飞. 缪子能量测量及其在缪子散射成像中的应用研究[D]. 北京: 清华大学工程物理系, 2016: 93−95. LUO Zhifei. Research on muon energy measurement and application in muon tomography[D]. Beijing: Tsinghua University. Department of Engineering Physics, 2016: 93−95.

[8] LESPARRE N, GIBERT D, MARTEAU J, et al. Geophysical muon imaging: feasibility and limits[J]. Geophysical Journal International, 2010, 183(3): 1348−1361.

[9] COBB J H, MARSHAK M L, ALLISON W W M, et al. Observation of a shadow of the Moon in the underground muon flux in the Soudan 2 detector[J]. Physical Review D, 2000, 61(9): 375−389.

[10] DARIJANI R, NEGARESTANI A, REZAIE M R, et al. A new approach in coal mine exploration using cosmic ray muons[J]. Acta Geophysica, 2016, 64(4): 1034−1050.

[11] SARACINO G, AMATO L, AMBROSINO F, et al. Imaging of underground cavities with cosmic-ray muons from observations at Mt. Echia (Naples)[J]. Scientific Reports, 2017, 7(1): 1−12.

[12] MORISHIMA K, KUNO M, NISHIO A, et al. Discovery of a big void in Khufu’s Pyramid by observation of cosmic-ray muons[J]. Nature, 2017, 552(7685): 386−390.

[13] MOLLERACH S, ROULET E. Progress in high-energy cosmic ray physics[J]. Progress in Particle and Nuclear Physics, 2018, 98: 85−118.

[14] BETHE H A. Moliere's theory of multiple scattering[J]. Physical Review, 1953, 89(6): 1256−1266.

[15] MOTZ J W, OLSEN H, KOCH H W. Electron scattering without atomic or nuclear excitation[J]. Reviews of Modern Physics, 1964, 36(4): 881−928.

[16] SCOTT W T. The theory of small-angle multiple scattering of fast charged particles[J]. Reviews of Modern Physics, 1963, 35(2): 231−313.

[17] SCHULTZ L J. Cosmic ray muon radiography[D]. Portland, USA: Portland State University, 2003: 104−145.

[18] XIAO S, HE W B, LAN M C, et al. A modified multi-group model of angular and momentum distribution of cosmic ray muons for thickness measurement and material discrimination of slabs[J]. Nuclear Science and Techniques, 2018, 29(2): 1−7.

[19] NAGAMINE K. Geo-tomographic observation of inner-structure of volcano with cosmic-ray muons[J]. Journal of Geography (Chigaku Zasshi), 1995, 104(7): 998−1007.

(编辑 陈灿华)

Research on muon advanced detection technology

XI Zhenzhu1, 2, LIU Junhua1, HE Hang1, WANG Liang3, PENG Xingliang1, ZHOU Sheng1, 3, WANG He1, 2

(1. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China; 2. Key Laboratory of Nonferrous Metal Metallogenic Prognosis of Ministry of Education, Changsha 410083, China; 3. Hunan 5D Geophyson Co. Ltd., Changsha 410205, China)

Based on that when cosmic ray muon penetrates high density material, the energy loss of cosmic ray muon is greater than that when it penetrates the low density material, the mathematical formula of ionization loss, bremsstrahlung, electron pair and nuclear interaction energy loss were deduced. Then, the Geant4 software was used to simulate the process that muon penetrated the karst cave filled air and water, record incident energy and ejection energy and construct the image of energy loss of muon, and muon low energy imaging was compared with its scattergram. The results show that muon energy loss imaging can effectively solve the problem of low resolution of muon scattering for elements with similar atomic number components, and muon low energy imaging is expected to be used in advanced detection with low density cavity filled with air and water, and it will provide a new way for advanced detection.

muon; energy loss; advanced detection; karst cave

10.11817/j.issn.1672-7207.2018.11.016

P631

A

1672−7207(2018)11−2753−06

2018−03−10；

2018−05−11

国家重点研发计划专项(2016YFC0303100)；国际海域资源调查与开发“十三五”计划项目(DY135-S1-1-07) (Project(2016YFC0303100) supported by the National Key R & D Plan; Project(DY135-S1-1-07) supported by the “13th Five-Year Plan” of International Sea Area Resources Survey and Development)

席振铢，博士，教授，从事电磁法勘探理论与应用研究；E-mail: xizhenzhu@163.com