2. 中国科学院高能物理研究所, 北京 100049;
3. 内蒙古自治区机电控制重点实验室, 内蒙古 呼和浩特 010051
2. Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China;
3. Inner Mongolia Key Laboratory of Electromechanical Control, Hohhot 010051, China
LHAASO-WCDA实验的主要物理目标是实现甚高能伽马射线的全天候观测，其中包括甚高能伽马源的探测与监测、能谱测量以及扩展形态的研究[1]。水切伦科夫探测器阵列的总面积为78 000 m2，有效水深4 m，由3个水池构成，共需要3.5 × 108 kg净化水。每个水池被黑色隔光帘分隔成多个探测器单元，共3 120个，每个单元大小为5 m × 5 m。在1号水池的每个单元底部中心位置安装一支光阴极朝上的8英寸光电倍增管(Photo Multiplier Tube, PMT)，用于接收广延大气簇射的次级粒子在水中产生的切伦科夫光。根据每个光电倍增管上接收的光子数目和光子到达时间，重建原初粒子的簇射芯位、能量和方向[2]。在每个8英寸光电倍增管旁边放置了一个1.5英寸的光电倍增管，用于扩大簇射芯区可测量的动态范围，从而实现高能宇宙线的高精度测量[3]。
水作为水切伦科夫探测器阵列的关键探测介质，对长期运行的探测器阵列实验至关重要，水的衰减长度直接影响探测器的探测效率，水衰减长度越大，表示水质越好，光探测效率越高。不同实验对水质的要求也不相同。切伦科夫探测器(CHerenkov Detectors In PitS, CHIPS)实验中测量切伦科夫光的衰减长度随滤光时间的变化，研究了405 nm激光的透射特性，衰减长度高达100 m[4]。中国大亚湾中微子实验中水的衰减长度要求30 m以上[5]。液氩时间投影室(Liquid Argon Time Projection Chamber, LArTPC)探测器实验中测量出液氩中掺杂的氮气分子与闪烁光之间的反应截面为(7.14 ± 0.74)×10-21cm2·molecule-1。百亿分之二(ppm水平)的掺杂浓度对应的衰减长度为(30 ± 3) m[6]。平方千米阵列(Kilometer-square Array, KM2A)缪子探测器[7]组自行研制了一种测量装置，在8 m长的水箱中安装多光源，达到快速测量水衰减长度的目的。通过对山上超纯水水样以及不同水质水样进行测量，验证该装置在100 m衰减长度的测量精度可以达到20%[8]。
水衰减长度测量方法中，采用一种简便易行的直筒式测量方法，通过改变测量装置中水位的高低来改变光在不同光程下的衰减，即可得到水样的衰减长度[9]。前期水切伦科夫探测器阵列工程样机实验中，通过数据研究分析发现了一套利用单路计数能谱的第2个峰的峰位变化监测水质的方法[10]。此方法需要结合直接测量装置的测量结果进行对比分析。GEANT4模拟时需要给定不同波长下的水衰减长度参数，即需要给出一个和实验接近的水质模型，但通过调研得出Smith和Querry两种水质模型[11]，需要通过测量确定一个合理的水质模型，为模拟提供可靠的参量。
实验采用直筒式测量方法得出一定波长下的水衰减长度，并同紫外可见分光光度计的测量结果进行对比，标定探测介质的衰减长度，并进一步通过不同波长发光二极管的测量结果分析，为GEANT4模拟提供一组可靠的水质模型参数。
1 水衰减长度测量装置 1.1 测量装置的工作原理水衰减长度测量装置由发光二极管、光筒、主水管、光电倍增管组成，在某一波长下测量不同水位h的光电倍增管接收的信号幅度A，A0为发光二极管的发光强度，根据水衰减长度定义λ：
| $ A = {A_0}{{\rm{e}}^{ - h/\lambda }}, $ | (1) |
| $ \ln A = \ln {A_0} + ( - h/\lambda ), $ | (2) |
得到光电倍增管接收到的信号幅值lnA与水位高度h的线性关系，进而得到水样的λ值，即斜率绝对值的倒数。
1.2 水切伦科夫探测器阵列水衰减长度测量装置图 1为水切伦科夫探测器阵列利用该测量原理设计的水衰减长度测量装置。测量步骤：(1)用待测水样自动清洗测量装置，以降低实验误差；(2)开启对应的电水阀和水泵开始注水，当水位达到1.1 m时，关闭电水阀和水泵，打开发光二极管和光电倍增管电源，光电倍增管采集发光二极管的发光信号；(3)控制排水阀，改变测量装置中水的高度h，测量光电倍增管接收到的信号幅值为A，记录保存实验数据；(4)根据水衰减长度计算公式拟合得到水样的水衰减长度。
|
| 图 1 水衰减长度测量装置结构图 Fig. 1 Water attenuation length measuring device structure diagram |
水衰减长度测量可以分为自动和手动测量两种方式。在慢控制的系统下，实现自动测量，单次水样测量每2小时完成一次，达到监测水质的目的。
1.3 水切伦科夫探测器阵列水衰减长度测量装置稳定性分析和测量结果在正式测量之前，对该装置的稳定性进行了验证分析。利用该装置对水池内的水质进行了27日的监测，通过分析每日的数据误差，验证系统测量的稳定性。图 2为27日的测量误差，可以看出误差波动在0.22 m范围内，说明本装置测量稳定。
|
| 图 2 水衰减长度测量装置的测量误差 Fig. 2 Measurement error of water attenuation length measuring device |
图 3利用此装置测得在400 nm波长发光二极管下水池水样的结果，图 3 (a)为单次测量拟合得到的水衰减长度，测量误差小于3%；图 3 (b)为长期测量结果，目前水池内的水质基本维持在(7.53 ± 0.21) m，总体测量误差小于5%。
|
| 图 3 测量装置的测量结果。(a)单次测量结果；(b)总体测量结果 Fig. 3 Measuring device measuring results. (a) Single measurement result; (b) Overall measurement result |
为了和直接测量装置对比研究，实验中还利用紫外可见分光光度计对同一水样进行测量，得出吸光度和测量装置给出的水衰减长度的对应关系。
2 紫外可见分光光度计工作原理 2.1 紫外可见分光光度计工作原理朗伯-比尔定律又称为光吸收基本定律。它表明物质吸收入射光时，吸收强度与光程及吸收物质浓度之间的关系，是紫外吸收光谱和红外光谱定量分析的理论依据。
| $ A = \lg \left( {{I_0}/I} \right) = \lg (1/T) = - \lg T = KCL, $ | (3) |
其中，A为吸光度；T为透过率；K为吸光系数；C为样品浓度；L为光程；I0/I为入射光强度/透过光强度。
紫外可见分光光度计由光源、单色器、样品室、检测器、放大控制系统和结果显示系统构成。实验采用双光束，有两条光路和一个检测器，两条光路中分别放置参比和样品，光从单色器发出后经斩波器，依次测定参比和样品。检测器选用光电倍增管，当它受到辐射时，吸收光子能量，并转变为可测量的物理量，将样品的测试结果显示保存在计算机上。该方法测量水样的吸光度是相对吸光度。
2.2 水切伦科夫探测器阵列紫外可见分光光度计本次水切伦科夫探测器阵列水质测量使用安捷伦科技有限公司生产的Cary300型紫外可见分光光度计。首先打开电脑Scan软件和Cary300设备，预热1小时，再设置测量波长300~700 nm，采样时间为0.2 s，采样间隔为1 nm，样品室放入纯水样品进行基线校正，再对水样进行测量。各个水样经紫外可见分光光度计测量的结果是吸光度(A)，即：
| $ A = \lg \left( {{I_0}/I} \right). $ | (4) |
依据光的衰减遵从比尔(Beer)定律，假定水的衰减长度为λ，那么强度为I0的入射光通过长度为L的水样后的强度为
| $ I = {I_0}{{\rm{e}}^{ - L/\lambda }}, $ | (5) |
分光计测量时L取0.1 m，可计算出各个水样在每个波长下的衰减长度为
| $ \lambda = 1/(10A\ln 10). $ | (6) |
在正式测量之前，利用该装置对水池内的水质进行了27日的监测，分析每日的数据误差，图 4为27日的吸光度转化为水衰减长度的测量误差，对该装置的稳定性进行了验证分析。可以看出，误差波动在0.18 m范围内，说明本装置测量同样稳定。
|
| 图 4 紫外可见分光光度计测量误差 Fig. 4 Measurement error with uv-visible spectrophotometer |
紫外可见分光光度计在波长400 nm时测量各个水样吸光度如图 5，纯水的吸光度基本为0，吸光度越小，表明水样的透明度越高，水衰减长度数值越大。开始注水时，吸光度快速下降，是不断增加的水量对水池内污染物的稀释作用。循环站处理后的水的吸光度小于水池内水的吸光度，表明循环站处理有一定的净化水质的能力，但尚未达到设计要求。
|
| 图 5 紫外设备在400 nm波长下对不同水样的长期监测结果 Fig. 5 Long-term monitoring results of different water samples at 400nm wavelength by ultraviolet equipment |
由上述两种测量装置的测试误差分析可以看出，两种测量装置均存在一定误差，分析是由于目前水切伦科夫探测器阵列水净化循环设备不稳定，影响测量。此外，紫外可见分光光度计测量时，由于设备本身误差和每次取样时比色皿上会有不同程度不可避免的污染导致测量结果略微不同。但两种装置测量的数据误差都很小，均在可接受范围内。随后利用两种装置对相同水样进行了测量及对比分析。
图 6是用水衰减长度测量装置在18日内对水切伦科夫探测器阵列水池内水的监测结果，水衰减长度是(7.41 ± 0.21) m。同时利用紫外可见分光光度计监测同一水样的结果经转化后，水衰减长度(6.79 ± 0.37) m。两种测量方法得到的水衰减长度结果在误差范围内是符合的。
|
| 图 6 水池内的水衰减长度 Fig. 6 The attenuation length of water in a pool |
对两种装置测量的数据进行了水质模型计算。表 1为通过水衰减长度测量装置和紫外可见分光光度计测量得到的水池内的水样在发光二极管波长365 nm、385 nm和400 nm下的水衰减长度，通过3个点得出水衰减长度与波长关系的斜率，分别为0.08 ± 0.04和0.10 ± 0.02。图 7是水的吸收长度随波长的变化关系图，根据不同的实验结果，给出了两种水吸收长度模型。Smith和Querry水质模型在相同条件下的斜率分别为0.79 ± 0.05和0.16 ± 0.01。两种装置测量的数据计算出的斜率，均明确得出水池内水样更接近Querry水质模型，因此，GEANT4模拟确定使用Querry水质模型。
|
| 图 7 两种水质模型 Fig. 7 Two water quality models |
Wavelength /nm |
λ/m | Gradient | |
| Water attenuation length measuring device for measuring water samples in pools | 365 385 400 |
4.53 5.07 7.53 |
0.08 |
| UV-visible spectrophotometer for measuring water samples in pools | 365 385 400 |
4.23 5.62 7.7 |
0.10 |
| Smith model | 360 370 380 390 400 |
26.24 33.01 44.20 51.10 56.84 |
0.79 |
| Querry model | 350 375 400 |
12.11 16.15 20 |
0.16 |
本文对水衰减长度测量装置和紫外可见分光光度计两种装置的稳定性进行了验证分析，得出了两种装置的测量误差分别为0.22 m和0.18 m。利用两种装置对相同水样进行了对比测量分析，证明两种装置在误差范围内，得到的监测结果相符。
利用水衰减长度测量装置和紫外可见分光光度计两种装置的测量结果，分析不同发光二极管波长下监测得到的水池水样的水衰减长度，通过计算及对比分析，确定了GEANT4模拟中应采用Querry水质模型。
| [1] | 中国科学院高能物理研究所, 中国科学院基础研究局. 西藏羊八井宇宙线国家野外观测研究站的建设、运行与发展[J]. 中国科学院院刊, 2010(4): 461–464 |
| [2] | 曹臻, 刘加丽, 白云翔. 物理学中的世纪难题:高能宇宙线的起源之"谜"[J]. 自然杂志, 2009, 31(6): 342–347, 363 |
| [3] | 曹臻, 陈明君, 陈松战, 等. 高海拔宇宙线观测站LHAASO概况[J]. 天文学报, 2019, 60(3): 1–16 |
| [4] | AMAT F, BIZOUARD P, BRYANT J, et al. Measuring the attenuation length of water in the CHIPS-M water Cherenkov detector[J]. Nuclear Instruments and Methods in Physics Research Section A:Accelerators Sepctrometers Detectors and Associated Equipment, 2017, 844: 108–115. |
| [5] | 路浩奇.大亚湾反符合水切伦科夫探测器及其模型研究[D].北京: 中国科学院研究生院, 2009. |
| [6] | JONES B J P, CHIU C S, CONRAD J M, et al. A measurement of the absorption of liquid argon scintillation light by dissolved nitrogen at the part-per-million level[J]. Journal of Instrumentation, 2013, 8(7): 1–19. |
| [7] | ZUO X, XIAO G, FENG S H, et al. Design and performances of prototype muon detectors of LHAASO-KM2A[J]. Nuclear Instruments and Methods in Physics Research Section A:Accelerators, Spectrometers, Detectors and Associated Equipment, 2015, 789: 143–149. |
| [8] | LI C, XIAO G, FENG S H, et al. An apparatus to measure water optical attenuation length for LHAASO-MD[J]. Nuclear Instruments and Methods in Physics Research Section A:Accelerators, Spectrometers, Detectors and Associated Equipment, 2018, 892: 122–126. |
| [9] | LI H C, YAO Z G, YU C X, et al. A method to monitor and measure the water transparency in LHAASO-WCDA using cosmic muon signals[J]. Chinese Physics C, 2017, 41(2): 145–152. |
| [10] | LI H C, YAO Z G, CHEN M J, et al. Study on single-channel signals of water Cherenkov detector array for the LHAASO project[J]. Nuclear Instruments and Methods in Physics Research Section A:Accelerators, Spectrometers, Detectors and Associated Equipment, 2017, 854: 107–112. |
| [11] | STEVE J, SCOTT P. Optical absorption of water compendium[DB/OL].[2019-11-04]. http://omlc.org/spectra/water/abs/index.html. |