趙 娟 高正明

?

支撐體對(duì)γ輻射物探測(cè)數(shù)據(jù)的影響

趙 娟1高正明2

1（荊楚理工學(xué)院電子信息工程學(xué)院 荊門 448000）2（寶雞高新技術(shù)研究所 寶雞 721013）

在對(duì)大質(zhì)量高安全性要求的γ源進(jìn)行輻射探測(cè)時(shí)，發(fā)現(xiàn)數(shù)據(jù)誤差很大，不能用于。為此，基于探測(cè)器端面注量率參數(shù)理論探析了地面或支撐臺(tái)對(duì)探測(cè)器數(shù)據(jù)影響的簡(jiǎn)化計(jì)算方法，采用蒙特卡羅方法分析了混凝土支撐體對(duì)粒子注量率的影響程度，研究了鉛作為次級(jí)粒子吸收物質(zhì)時(shí)，鉛厚度和半徑與注量率測(cè)量誤差之間的關(guān)系。研究結(jié)果表明，在放射性材料輻射探測(cè)過程中，為減少混凝土支撐體對(duì)輻射探測(cè)數(shù)據(jù)的影響，應(yīng)采用厚度不小于0.5 cm、半徑不小于源高3倍的鉛板作為次級(jí)粒子吸收體，可將探測(cè)數(shù)據(jù)最大誤差由22%減少至15%左右，最小可接近10%，實(shí)驗(yàn)結(jié)果證實(shí)了該方法的有效性。

支撐體，放射性材料，γ輻射探測(cè)，蒙特卡羅方法

由于實(shí)際實(shí)驗(yàn)條件限制，輻射探測(cè)實(shí)驗(yàn)數(shù)據(jù)與理論真值之間存在著偏差，偏差的數(shù)值大小決定著實(shí)驗(yàn)數(shù)據(jù)的可信程度，估算輻射探測(cè)數(shù)據(jù)偏差是數(shù)據(jù)處理與分析的重要組成部分，在核技術(shù)應(yīng)用研究中具有重要的地位。

輻射探測(cè)數(shù)據(jù)與真值之間的偏差來源于多個(gè)方面，如探測(cè)器自身性能、輻射粒子特性和實(shí)驗(yàn)環(huán)境因素等。在探測(cè)器制造過程中，對(duì)晶體進(jìn)行表面工藝處理有助于提升探測(cè)器性能，如晶體表面光滑與蝕刻工藝處理等[1?2]。在輻射探測(cè)數(shù)據(jù)處理過程中，實(shí)驗(yàn)儀器已經(jīng)選定，因此人們對(duì)實(shí)驗(yàn)環(huán)境條件的影響更為關(guān)注。Bartlett等[3]研究了鹽沉積對(duì)α放射性氣溶膠連續(xù)監(jiān)測(cè)儀的影響，指出由于沉積鹽層對(duì)α粒子的阻止作用，α放射性活度的探測(cè)精度將會(huì)降低；陳凌等[4]研究了大氣壓強(qiáng)或空氣密度對(duì)固體核徑跡氡探測(cè)器刻度系數(shù)的影響，得出大氣壓強(qiáng)降低或空氣密度增大可導(dǎo)致氡徑跡探測(cè)器刻度系數(shù)增大17%到31%不等的結(jié)論。此外，人們還研究了探測(cè)目標(biāo)的輸運(yùn)與運(yùn)動(dòng)狀態(tài)的影響程度，顧行發(fā)等[5]根據(jù)輻射在地表－大氣界面的反射和折射性質(zhì)，采用相對(duì)誤差因子定量分析了地表反射類型變化對(duì)輻射矢量的影響，認(rèn)為地表反射對(duì)光輻射的影響較大，數(shù)據(jù)處理需要考慮地表的方向反射特性和自然界的偏振效應(yīng)；稅舉等[6]用蒙特卡羅方法研究了放射性材料移動(dòng)對(duì)輻射探測(cè)陣列探測(cè)質(zhì)量下限的影響，推論在600 cm距離，探測(cè)勻速移動(dòng)時(shí)的放射性材料下限比靜止的下限大一倍。

實(shí)際放射性材料輻射探測(cè)過程中，出于安全性考慮，質(zhì)量較大、輻射較強(qiáng)或不能移動(dòng)等原因，有時(shí)不能將源吊于空中[7]進(jìn)行輻射探測(cè)，而只能通過限定的操作將源放置于地面或工作臺(tái)上[8]。此時(shí)，地面或工作臺(tái)物質(zhì)與輻射粒子相互作用，散射或產(chǎn)生的次級(jí)粒子對(duì)探測(cè)數(shù)據(jù)有一定的影響。為此，本文采用蒙特卡羅方法分析放射性材料放置于混凝土地面上時(shí)的γ輻射探測(cè)數(shù)據(jù)誤差，并采用常用的鉛屏蔽材料作為次級(jí)粒子吸收體，計(jì)算鉛層厚度與覆蓋面尺寸與探測(cè)數(shù)據(jù)誤差之間的關(guān)系。通過改變輻射探測(cè)環(huán)境減少數(shù)據(jù)偏差，為特定條件下放射性材料輻射特性參數(shù)測(cè)量提供技術(shù)支撐。

1 影響機(jī)理分析

1.1 原理分析

γ粒子與物質(zhì)相互作用的主要形式是光電效應(yīng)、康普頓效應(yīng)和電子對(duì)生成效應(yīng)等。當(dāng)γ粒子能量較低時(shí)，γ粒子主要與原子殼層轉(zhuǎn)道上的內(nèi)層電子碰撞產(chǎn)生光電子；當(dāng)γ粒子能量較高時(shí)，一部分γ粒子與自由電子發(fā)生碰撞，將一部分能量轉(zhuǎn)移給電子，而自身能量和運(yùn)動(dòng)方向發(fā)生改變；當(dāng)γ粒子能量大于1.022 MeV時(shí)，一部分γ粒子經(jīng)過原子核附近時(shí)在庫(kù)侖電場(chǎng)的作用下轉(zhuǎn)換成一對(duì)正、負(fù)電子。次級(jí)電子通過電離和激發(fā)、軔致輻射將能量傳給原子或轉(zhuǎn)換為光子，直至能量完全被物質(zhì)吸收或逃逸。

在放射性材料的γ輻射探測(cè)實(shí)驗(yàn)中，到達(dá)探測(cè)器的粒子來源有：(1) 源發(fā)射出的未與物質(zhì)發(fā)生相互作用的粒子；(2) 源粒子與其它粒子碰撞后改變方向到達(dá)的源輻射粒子；(3) 源粒子與物質(zhì)相互作用產(chǎn)生的各級(jí)電子在空氣或支撐體中軔致輻射產(chǎn)生的光子；(4) 源粒子與物質(zhì)相互作用產(chǎn)生的各級(jí)正電子湮滅光子；(5) 源粒子與物質(zhì)相互作用產(chǎn)生的其它如瑞利散射效應(yīng)產(chǎn)生的光子；(6) 環(huán)境本底。這些粒子中只有(1)是期望測(cè)量的。

1.2 地面或支撐臺(tái)影響程度的數(shù)值分析

地面或支撐臺(tái)與源粒子相互作用，產(chǎn)生的次級(jí)產(chǎn)物隨著遠(yuǎn)離源距離的增大而減弱。假定柱源源強(qiáng)，半徑，高，各向同性。以柱源與地面或支撐臺(tái)接觸的圓面中心為原點(diǎn)，中心垂線為軸，建立直角坐標(biāo)系。忽略圓柱截面粒子作用，則該柱源放置于地面或支撐臺(tái)上時(shí)，軸上任一點(diǎn)(0,0,0)能夠接收來自于源圓柱面部分的粒子：

曲面上的任一點(diǎn)(,,)鄰域d在點(diǎn)(0,0,0)的粒子注量率為：

根據(jù)軸對(duì)稱性和式，對(duì)式進(jìn)行曲面積分，可得：

式(3)可根據(jù)積分的數(shù)值求解法進(jìn)行求解，如圖1所示。

若忽略次級(jí)粒子角分布和能量分布規(guī)律的影響，認(rèn)為地面或支撐臺(tái)與源粒子相互作用產(chǎn)生的次級(jí)光子注量率正比于源粒子在該點(diǎn)的注量率，并以各向同性分布被探測(cè)器探測(cè)，此時(shí)探測(cè)器的幾何因子可采用蒙特卡羅方法進(jìn)行計(jì)算，抽樣點(diǎn)(,0,)滿足：

探測(cè)器截面半徑，截面圓表征曲線：

以概率=·()（其中為比例系數(shù)，用于表征次級(jí)粒子注量率與源粒子注量率的比值）模擬抽樣個(gè)粒子，飛行至=，(?)2+(?)2≤2的粒子數(shù)，則探測(cè)器幾何因子=/。

位于=處的探測(cè)器，地面或支撐臺(tái)次級(jí)產(chǎn)物對(duì)其探測(cè)數(shù)據(jù)的貢獻(xiàn)可簡(jiǎn)單表征為：

即地面或支撐臺(tái)的存在使探測(cè)器注量率測(cè)量值增大了D。

實(shí)際上，源粒子與地面或支撐臺(tái)物質(zhì)相互作用形式復(fù)雜，次級(jí)產(chǎn)物與物質(zhì)的作用形式多樣，且數(shù)學(xué)解析難度大，因此地面或支撐臺(tái)對(duì)源探測(cè)真值的影響一般采用數(shù)值方法進(jìn)行分析。

2 影響程度與對(duì)策分析

2.1 仿真環(huán)境與模型

MCNP程序是常用的粒子輸運(yùn)計(jì)算軟件之一，可用于處理多種粒子的輸運(yùn)計(jì)算，本文基于MCNP5對(duì)γ粒子輸運(yùn)過程進(jìn)行數(shù)值仿真分析。仿真分析過程中采用兩種計(jì)算模型：一是放射性材料自由位于計(jì)算空間，周邊環(huán)境為真空或空氣(圖2)，可用于理想條件下粒子通過真空或空氣環(huán)境后到達(dá)探測(cè)器真值的數(shù)值分析；二是無限厚混凝土基礎(chǔ)支撐，放射性材料位于混凝土或支撐板上(圖3)，用于計(jì)算空曠區(qū)域放射性材料放置于地面或地面支撐板上時(shí)的探測(cè)器計(jì)數(shù)。測(cè)量對(duì)象為60Co圓柱源，該源發(fā)射1.17MeV (99.85%)、1.332 MeV (99.98%)兩種能量的γ粒子。

2.2 混凝土地面支撐時(shí)γ輻射探測(cè)數(shù)據(jù)誤差分析

采用探測(cè)器前端面γ粒子注量率進(jìn)行評(píng)價(jià)[9]，由圖2模型仿真得到探測(cè)器前端面注量率真值和空氣環(huán)境的影響程度，采用圖2模擬放射性材料放置于地面或支撐的情況，結(jié)果如圖4、5所示。

可見，理想條件下，探測(cè)器前端面粒子注量率受空氣的影響較弱，空氣介質(zhì)的存在對(duì)粒子注量率的影響不大于5%。當(dāng)放射性材料放置于地面時(shí)，源粒子與混凝土物質(zhì)成分相互作用，產(chǎn)生了大量的次級(jí)粒子，使探測(cè)器前端面注量率受到顯著影響，仿真10?400 cm探測(cè)器前端面粒子注量率，平均偏差16.6%，最大偏差可達(dá)22%。隨著探測(cè)器逐漸遠(yuǎn)離放射源，計(jì)數(shù)誤差先增大后減小(圖4)，這是因?yàn)楫?dāng)兩者之間距離較小時(shí)，隨著探測(cè)器遠(yuǎn)離源物質(zhì)，探測(cè)器對(duì)混凝土地面所張立體角逐漸增大，放射源發(fā)出的γ粒子與混凝土相互作用產(chǎn)生的次級(jí)γ粒子和散射粒子到達(dá)探測(cè)器前端面的量增多；當(dāng)探測(cè)器與放射源足夠遠(yuǎn)時(shí)，能量較小的次級(jí)γ粒子被吸收或散射，不能到達(dá)探測(cè)器，同時(shí)探測(cè)器對(duì)次級(jí)粒子較多的近源區(qū)域所張的立體角減小，因此到達(dá)探測(cè)器前端面的粒子越來越多地來自于放射源，使探測(cè)器前端面粒子注量率更接近真實(shí)值。

圖4 前端面γ粒子注量率與測(cè)量距離的關(guān)系 Fig.4 The γ fluence rate on the front plate of the detectors vs. distances. 圖5 實(shí)際探測(cè)環(huán)境中注量率誤差與距離的關(guān)系 Fig.5 Error of the fluence rate vs. distances.

可見，當(dāng)特定條件下放射性材料只能放置于地面（或工作臺(tái)）上進(jìn)行輻射探測(cè)時(shí)，混凝土地面對(duì)γ輻射探測(cè)數(shù)據(jù)的影響很大。此時(shí)必須考慮地面或支撐體對(duì)輻射探測(cè)數(shù)據(jù)的影響，采取措施降低數(shù)據(jù)偏差。

2.3 分析與對(duì)策

由前述分析可知，測(cè)量放射性材料時(shí)，空氣介質(zhì)對(duì)探測(cè)數(shù)據(jù)的影響較小，測(cè)量偏差主要來源于源粒子與地面支撐體的相互作用產(chǎn)生的次級(jí)粒子的貢獻(xiàn)。因此，若采用一層γ粒子吸收層，盡可能減少到達(dá)地面或支撐體的γ粒子散射產(chǎn)物或其與支撐體相互作用產(chǎn)生的次級(jí)產(chǎn)物到達(dá)探測(cè)器，就可有效減少探測(cè)誤差。

鉛是常用的γ粒子屏蔽材料，可通過加工制作一定厚度的鉛板放置于放射性材料與地面之間，吸收放射性材料發(fā)射向地面的粒子，減少非期望粒子對(duì)探測(cè)數(shù)據(jù)的貢獻(xiàn)，降低輻射探測(cè)誤差[10]。仿真計(jì)算時(shí)在放射性材料與地面間放置一層鉛板，其厚度分別為0.2?1.0 cm不等，隨著鉛層厚度的不同，探測(cè)器前端面注量率相對(duì)于真空狀態(tài)數(shù)據(jù)的偏差也不同，結(jié)果如圖6、7所示。

隨鉛層厚度的增大，入射向地面的源粒子更多地被鉛層吸收，因此相應(yīng)計(jì)數(shù)比無鉛層時(shí)的數(shù)據(jù)更接近真值，但由于鉛本身也和γ粒子發(fā)生相互作用，產(chǎn)生散射光子或次級(jí)光子，因此添加鉛吸收層并不能完全消除誤差。需要注意的是，當(dāng)鉛層能夠完全吸收混凝土地面產(chǎn)生的次級(jí)光子時(shí)，再增大鉛層厚度已經(jīng)無明顯效果了。

鉛的密度很大，會(huì)導(dǎo)致中毒，因此實(shí)際使用時(shí)應(yīng)盡量減少鉛的用量，據(jù)圖6、7結(jié)果可選擇0.5 cm厚度的鉛層作為吸收層。因200 cm處探測(cè)數(shù)據(jù)誤差最大，以該點(diǎn)為參考，分析鉛層覆蓋面與200 cm處探測(cè)器前端面注量率誤差之間的關(guān)系，如圖8所示。

由圖8可知，隨著鉛板半徑的增大，200 cm處探測(cè)器前端面注量率誤差逐漸下漸，下降速度先增大后減小。因此，在條件許可時(shí)，可將源和探測(cè)器全部采用鉛板支撐，以減少次級(jí)粒子的影響，此時(shí)探測(cè)器前端面粒子注量率誤差可減少于10.1%；鉛量不足時(shí)，至少應(yīng)確保源高3倍以上吸收層半徑，此時(shí)前端面粒子注量率誤差仍然較大，但再小幅度增加覆蓋面尺寸（鉛板半徑在源高3?6倍之間）時(shí)的性價(jià)比較差。

2.4 實(shí)驗(yàn)驗(yàn)證

為驗(yàn)證仿真分析結(jié)果，采用澳大利亞Canberra公司生產(chǎn)的Inspector 1000 γ探測(cè)器對(duì)放置于室外開闊區(qū)的137Cs源進(jìn)行實(shí)驗(yàn)測(cè)量。加工了3個(gè)

100 cm×100 cm×0.3 cm的鉛板，放置于源和混凝土地面之間，分別測(cè)量源放置于混凝土地面、0.3 cm、0.6cm、0.9 cm鉛板時(shí)，離源100?300 cm之間間隔20 cm不同位置處的γ劑量率數(shù)值。測(cè)量結(jié)果如圖9所示。

該137Cs源標(biāo)定于1991年，長(zhǎng)期使用后活度變化情況未知，因此實(shí)驗(yàn)中未計(jì)算相對(duì)誤差，但實(shí)驗(yàn)數(shù)據(jù)分布和仿真數(shù)據(jù)具有較強(qiáng)的相似性，0.6 cm厚度鉛屏蔽后的輻射探測(cè)數(shù)據(jù)已顯著降低了誤差，證實(shí)了本文分析結(jié)果。

圖9 實(shí)驗(yàn)探測(cè)數(shù)據(jù) Fig.9 Experimental results.

3 結(jié)語

據(jù)實(shí)際放射性材料輻射探測(cè)需求，本文以60Co圓柱源為分析對(duì)象，采用探測(cè)器朝向源的前端面的粒子注量率為特征參數(shù)，基于蒙特卡羅方法分析了實(shí)際探測(cè)環(huán)境對(duì)源探測(cè)數(shù)據(jù)的影響。研究結(jié)果表明：(1) 在不考慮支撐的情況下，空氣環(huán)境對(duì)γ粒子注量率的影響較少，與真空狀態(tài)數(shù)據(jù)相比偏差不大于5%；(2) 實(shí)際測(cè)量中，位于混凝土地面的源的γ輻射探測(cè)誤差較大，與真空狀態(tài)數(shù)據(jù)相比，平均偏差16.6%，最大誤差達(dá)22%，數(shù)據(jù)可信度低；(3) 采用0.5 cm厚的鉛板可有效降低混凝土地面的影響，為獲得更好的效果，板半徑應(yīng)盡量包圍源和探測(cè)器，至少應(yīng)大于源高3倍；(4) 測(cè)量誤差先增大后減少，以本文仿真設(shè)定，增大鉛板覆蓋面，可將探測(cè)數(shù)據(jù)最大誤差降低至15%左右，隨鉛板覆蓋面的進(jìn)一步增大，理論上最大誤差可減少至10.1%。

數(shù)值分析結(jié)果表明，采用一定厚度的鉛板可顯著減少源粒子與地面混凝土成分相互作用產(chǎn)生的次級(jí)粒子的影響。因此，在特定條件下，當(dāng)放射性材料只能放置于地面或工作臺(tái)上進(jìn)行γ輻射探測(cè)時(shí)，在源與地面或工作臺(tái)之間添加一層厚度不小于0.5cm的鉛板，可顯著提高測(cè)量準(zhǔn)確度。

1 Oliveira I B, Costa F E, Brazil S P,. Influence of crystalline surfaces on the TlBr radiation detector performance[J]. IEEE Transactions on Nuclear Science, 2005, 52(5): 2058?2062. DOI: 10.1109/TNS.2005.56788

2 Zheng Q, Dierre F, Crocco J,. Influence of surface preparation on CdZnTe nuclear radiation detectors[J]. Applied Surface Science, 2011, 257(20): 8742?8746. DOI: 10.1016/j.apsusc.2011.05.098

3 Bartlett W T, Walker B A. The influence of salt aerosol on alpha radiation detection by WIPP continuous air monitors[A]. 24thDOE/NRC nuclear air cleaning and treatment conference[C]. Oregon, Portland: The Harvard Air Cleaning Laboratory, 1997: 188?199. DOI: 10.2172/ 179208

4 陳凌, 劉森林, 武奕華, 等. 固體核徑跡氡探測(cè)器刻度系數(shù)及其受大氣壓強(qiáng)影響的蒙特卡羅模擬計(jì)算[J]. 原子能科學(xué)技術(shù), 2004, 38(1): 74?78. DOI: 10.3969/j.issn. 1000-6931.2004.01.015

CHEN Ling, LIU Senlin, WU Yihua,. Monte Carlo computation of the calibration factor of solid state nuclear track detector for radon and the effect by atomsphericpressure[J]. Atomic Energy Science and Technology, 2004, 38(1): 74?78. DOI: 10.3969/j.issn. 1000-6931.2004.01.015

5 顧行發(fā), 程天海, 余濤, 等. 地表雙向反射對(duì)天基矢量輻射探測(cè)的影響分析[J]. 物理學(xué)報(bào), 2009, 58(10): 7368?7375. DOI: 10.7498/aps.58.7368

GU Xingfa, CHENG Tianhai, YU Tao,. Effect of surface reflectances on the space-based vector radiative detection[J]. Acta Physica Sinica, 2009, 58(10): 7368?7375. DOI: 10.7498/aps.58.7368

6 稅舉, 何彬, 付廣智, 等. 核材料移動(dòng)對(duì)輻射探測(cè)器陣列探測(cè)下限影響的蒙特卡羅模擬[J]. 原子能科學(xué)技術(shù), 2007, 41(4): 497?500. DOI: 10.3969/j.issn.1000-6931. 2007.04.024

SHUI Ju, HE Bin, FU Guangzhi,. Monte-Carlo simulation of movement of nuclear matter influencing on minimum detectable limit of radiation detector array[J]. Atomic Energy Science and Technology, 2007, 41(4): 497?500. DOI: 10.3969/j.issn.1000-6931.2007.04.024

7 蔣詩(shī)平.241Am-Be中子源快中子成像研究[J]. 同位素, 2007, 20(2): 87?90. DOI: 10.7538/tws.2007.20.02.0087

JIANG Shiping. Fast neutron imaging using241Am-Be neutron source[J]. Journal of Isotopes, 2007, 20(2): 87?90. DOI: 10.7538/tws.2007.20.02.0087

8 Lu Z Y, Gao L, Peng G X. Density isomer of nuclear matter in an equivalent mass approach[J]. Nuclear Science and Techniques, 2013, 24(5): 050505. DOI: 10.13538/j.1001-8042/nst.24.050505

9 胡廣春, 郝樊華, 劉素萍, 等. 钚體源樣品γ能譜計(jì)算的蒙特卡羅方法[J]. 物理學(xué)報(bào), 2005, 54(8): 3523?3529. DOI: 10.7498/aps.54.3523

HU Guangchun, HAO Fanhua, LIU Suping,. Monte-Carlo method in calculating the γ spectrum of plutonium volume source[J]. Acta Physica Sinica, 2005, 54(8): 3523?3529. DOI: 10.7498/aps.54.3523

10 Ochbelagh D R, Azimkhani S, Mosavinejad H G. Gamma-ray shielding properties of concrete with different percentages of lead[J]. Nuclear Science and Techniques, 2012, 23(4): 215?218. DOI: 10.13538/j.1001-8042/nst. 23.215?218

國(guó)防預(yù)先研究項(xiàng)目基金(No.51331040501)資助

Supported by National Defense Pre-research Projects (No.51331040501)

Influence of the support body on radiation detection for gamma radiator

ZHAO Juan1GAO Zhengming2

1(School of Electronics and Information Engineering, Jingchu University of Technology, Jingmen 448000, China)2(Baoji Research Institute of High Technology, Baoji 721013, China)

Background: Theof the gamma detection on large-mass and high security controlled radioactive materials is too large for quantitative analysis. Purpose: This study aims to accurately measure the radiation of gamma radioactive materials placed upon support body. Methods & Results: Based on the theory of particle flux parameters on the cross section of detection, the method to calculate the influence of the support body on radiation detection was simplified. Then, simulations were done to find the error distributions with Monte Carlo method when the materials were placed on the concrete ground. In order to reduce the detection error, a plate of lead was recommended to be inserted between the materials and the concrete ground. The best thickness of the lead layer was simulated and experimental verification was performed. Results: Simulation results show that if a lead layer thicker than 0.5 cm exists between the radioactive materials and the concrete ground. The relative detection error could be reduced from 22% to 15% or even the least of 10% in some particular situations. Conclusion: If a large-mass or high security controlled gamma radiator would be detected, it should be placed on a lead layer thicker than 0.5 cm in case of reducing the datum errors.

Support body, Radioactive material, Gamma radiation detection, Monte Carlo method

ZHAO Juan, female, born in 1979, graduated from Huazhong University of Science and Technology with a master’s degree in 2009, focusing on the blind equalization, the non-destructive testing and monitoring technologies

TL81，TL67

10.11889/j.0253-3219.2016.hjs.39.020401

趙娟，女，1979年出生，2009年于華中科技大學(xué)獲碩士學(xué)位，研究領(lǐng)域?yàn)橥ㄐ畔到y(tǒng)盲均衡技術(shù)和無損檢測(cè)與監(jiān)控技術(shù)

2015-10-12，

2015-12-16