數(shù)字化核能譜測(cè)量系統(tǒng)中的脈沖基線(xiàn)估計(jì)方法
【專(zhuān)利摘要】本發(fā)明公開(kāi)一種數(shù)字化核能譜測(cè)量系統(tǒng)中的脈沖基線(xiàn)估計(jì)方法，該方法以FPGA芯片為核心，由FPGA控制高速模數(shù)轉(zhuǎn)換芯片進(jìn)行核輻射信號(hào)采樣，在FPGA中引入循環(huán)緩沖區(qū)的方法來(lái)實(shí)現(xiàn)數(shù)據(jù)滑動(dòng)平均濾波，先對(duì)所有采樣數(shù)據(jù)進(jìn)行滑動(dòng)平均濾波求得粗基線(xiàn)，結(jié)合粗基線(xiàn)和一階數(shù)值微分法判斷并扣除采樣數(shù)據(jù)中的異常信號(hào)（包括脈沖或干擾），對(duì)扣除異常信號(hào)后的數(shù)據(jù)再次使用滑動(dòng)平均濾波實(shí)現(xiàn)高精度的脈沖基線(xiàn)估計(jì)。本發(fā)明所述去異常脈沖基線(xiàn)估計(jì)方法，降低了算法復(fù)雜度，提高了運(yùn)算速度，不受溫度等干擾因素的影響，能夠?qū)崟r(shí)準(zhǔn)確跟蹤基線(xiàn)的變化。
【專(zhuān)利說(shuō)明】
數(shù)字化核能譜測(cè)量系統(tǒng)中的脈沖基線(xiàn)估計(jì)方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明屬于核輻射測(cè)量領(lǐng)域，可應(yīng)用于數(shù)字化核能譜測(cè)量系統(tǒng)中的高精度脈沖基 線(xiàn)估計(jì)及穩(wěn)譜處理，特別適用于核測(cè)井等環(huán)境溫度變化較大情況下的能譜測(cè)量。
【背景技術(shù)】
[0002] 在核能譜測(cè)量系統(tǒng)中，由于探測(cè)器漏電流、極零相消處理的不完善、電子學(xué)元器件 溫漂、電源紋波、以及數(shù)字信號(hào)對(duì)模擬信號(hào)的高頻干擾等影響，使得在核脈沖幅度測(cè)量過(guò)程 中，有用信號(hào)總是疊加在一個(gè)不穩(wěn)定的直流基線(xiàn)上，從而導(dǎo)致測(cè)量結(jié)果的波動(dòng)。不僅降低了 系統(tǒng)的能量分辨率，而且會(huì)引起能量線(xiàn)性刻度曲線(xiàn)的斜率變化和直流基線(xiàn)漂移。特別是在 測(cè)井環(huán)境中，因井下溫度變化較大導(dǎo)致系統(tǒng)的基線(xiàn)漂移嚴(yán)重。如果能夠進(jìn)行準(zhǔn)確的基線(xiàn)估 計(jì)和扣除，可確保系統(tǒng)的能量刻度曲線(xiàn)經(jīng)過(guò)坐標(biāo)軸的零點(diǎn)，這必然會(huì)降低后續(xù)譜線(xiàn)處理的 復(fù)雜度。近年來(lái)，傳統(tǒng)的模擬核能譜測(cè)量系統(tǒng)逐漸被數(shù)字化核能譜測(cè)量系統(tǒng)所取代，相應(yīng)的 能譜測(cè)量過(guò)程中的脈沖基線(xiàn)處理方法也由傳統(tǒng)的基線(xiàn)恢復(fù)電路發(fā)展到數(shù)字基線(xiàn)估計(jì)，數(shù)字 基線(xiàn)估計(jì)與模擬的基線(xiàn)回復(fù)電路相比，可顯著提高系統(tǒng)穩(wěn)定性和可靠性，具有較高的靈活 性和適應(yīng)性，可進(jìn)一步改善系統(tǒng)能量分辨率，從而成為數(shù)字化核信號(hào)處理方法中特別是數(shù) 字多道分析中的一個(gè)重要部分。
[0003] 目前國(guó)內(nèi)外對(duì)數(shù)字基線(xiàn)估計(jì)和扣除的研究方法有很多，比如通過(guò)插值的方法求出 新的基線(xiàn)采樣值來(lái)代替原來(lái)的脈沖信號(hào)采樣值，采用分段濾波法中的一個(gè)信號(hào)濾波器和三 個(gè)基線(xiàn)濾波器實(shí)現(xiàn)基線(xiàn)估計(jì)值。國(guó)內(nèi)學(xué)者發(fā)表的文獻(xiàn)主要包括:《原子能科學(xué)技術(shù)》的《核輻 射脈沖幅度分析的基線(xiàn)卡爾曼濾波估計(jì)》，《核電子學(xué)與探測(cè)技術(shù)》的《多道脈沖幅度分析中 的數(shù)字基線(xiàn)估計(jì)方法》。此類(lèi)技術(shù)主要采用了FIR基線(xiàn)濾波器、自適應(yīng)Kalman濾波器實(shí)現(xiàn)脈 沖基線(xiàn)估計(jì)，該類(lèi)方法算法復(fù)雜，多是在計(jì)算機(jī)上完成，較難實(shí)現(xiàn)小型化和集成化。《核技 術(shù)》的《基于FPGA脈沖幅度分析器的數(shù)字化基線(xiàn)估計(jì)方法》提出使用直方圖統(tǒng)計(jì)法、平均值 法在FPGA上實(shí)現(xiàn)脈沖基線(xiàn)估計(jì)，該方法需要對(duì)基線(xiàn)進(jìn)行較長(zhǎng)時(shí)間的跟蹤，且判斷過(guò)程較為 復(fù)雜，而且在估計(jì)基線(xiàn)時(shí)沒(méi)有完全剔除正常核脈沖信號(hào)對(duì)基線(xiàn)的影響。
[0004] 利用在FPGA上引入去異?；瑒?dòng)平均濾波法來(lái)估算基線(xiàn)，降低了算法的復(fù)雜度，精 度好，不受溫度等干擾因素的影響，能夠?qū)崟r(shí)準(zhǔn)確跟蹤基線(xiàn)的變化。迄今為止，尚未見(jiàn)到將 該項(xiàng)技術(shù)直接應(yīng)用于數(shù)字化能譜測(cè)量的報(bào)道。

【發(fā)明內(nèi)容】

[0005] 本發(fā)明的目的是為了解決數(shù)字化核能譜測(cè)量系統(tǒng)中的能譜測(cè)量精度和穩(wěn)定性問(wèn) 題，而提出一種脈沖基線(xiàn)估計(jì)方法。
[0006] 本發(fā)明的技術(shù)方案為:本發(fā)明的目的是通過(guò)下述技術(shù)方案實(shí)現(xiàn)的。
[0007] 數(shù)字化核能譜測(cè)量系統(tǒng)中的脈沖基線(xiàn)估計(jì)方法，首先在FPGA芯片內(nèi)求得粗基線(xiàn)：
[0008] (1)在FPGA芯片內(nèi)開(kāi)辟一組大小為N=2k個(gè)單元的循環(huán)緩沖區(qū)A，在執(zhí)行平均濾波 之前，將緩沖區(qū)的內(nèi)容全部初始化為〇;
[0009] (2)每次采樣得到的數(shù)據(jù)x(n)存入循環(huán)緩沖區(qū)A頭指針指向的位置，剔除循環(huán)緩沖 區(qū)A尾指針指向的數(shù)據(jù)x(n-N-l)，并利用如下公式求得粗基線(xiàn)? :
[0010] Sum = SumPre+x(n)-x(n-N-l)
[0011] ^ = Sum / N = Sum / 2/' = Sum ? k
[0012] 式中Sum為循環(huán)緩沖區(qū)A中N個(gè)數(shù)據(jù)的累加和，SumPre為Sum延遲一個(gè)時(shí)鐘序列的 值，SumPr e的初始值為0，Sum >> k表示將Sum右移k個(gè)二進(jìn)制位。
[0013] 其次根據(jù)粗基線(xiàn)扣除采樣數(shù)據(jù)中的異常信號(hào)：
[0014] (3)將采樣數(shù)據(jù)與粗基線(xiàn)進(jìn)行比較，如果大于粗基線(xiàn)，則認(rèn)為該數(shù)據(jù)是疊加在基線(xiàn) 上的異常信號(hào)。如果小于或等于粗基線(xiàn)（即，則對(duì)其進(jìn)行一階數(shù)值微分，微分結(jié)果 用dx(n)表示，公式為：
[0015] dx(n) =χ(η)-χ(η-1)
[0016] 式中χ(η)表示當(dāng)次采樣數(shù)據(jù)，x(n-l)表示前一次采樣數(shù)據(jù)；
[0017] (4)選取p個(gè)數(shù)據(jù)來(lái)判斷輸入數(shù)據(jù)進(jìn)行一階數(shù)值微分后，其數(shù)據(jù)變化幅度大小及正 負(fù)符號(hào)是否連續(xù)，如果微分后的數(shù)據(jù)變化幅度較小且正負(fù)符號(hào)不連續(xù)，則認(rèn)為該數(shù)據(jù)為基 線(xiàn)，通過(guò)下列公式求得扣除異常信號(hào)后的基線(xiàn)數(shù)據(jù)V (η):
[0018]

[0019] 式中f為已經(jīng)求得的粗基線(xiàn)值
3連續(xù)P個(gè)數(shù)據(jù)一階微分后的數(shù)值求 和，I dx(n) I表示數(shù)值微分dx(n)的絕對(duì)值，ε表示輸入數(shù)據(jù)的紋波范圍。
[0020] 最后根據(jù)扣除異常信號(hào)后的數(shù)據(jù)求得精密基線(xiàn)：
[0021] (5)在FPGA芯片內(nèi)開(kāi)辟一組大小為M=2m個(gè)單元的循環(huán)緩沖區(qū)Β，在執(zhí)行平均濾波 之前，將緩沖區(qū)的內(nèi)容全部初始化為〇;
[0022] (6)將扣除異常信號(hào)后的基線(xiàn)數(shù)據(jù)Y (η)存入循環(huán)緩沖區(qū)B頭指針指向的位置，剔 除循環(huán)緩沖區(qū)B尾指針指向的數(shù)據(jù)V (η-Ν-1)，并利用如下公式求得精密基線(xiàn)玄：
[0023] Sum7 =SumPre^x7 (Ii)-X7 (η-Ν-1)
[0024] X - Sum ? M = Sum ? 2"' = Suin' ? m
[0025] 式中Sun/為循環(huán)緩沖區(qū)B中M個(gè)數(shù)據(jù)的累加和，SumPrV為Sun/延遲一個(gè)時(shí)鐘序列 的值，SumPre 7的初始值為0,Sun/ >>m表示將Sun/右移m個(gè)二進(jìn)制位。
[0026] 本發(fā)明的優(yōu)點(diǎn)在于：（1)在FPGA中引入循環(huán)緩沖區(qū)來(lái)實(shí)現(xiàn)滑動(dòng)平均濾波算法后，算 法復(fù)雜度顯著降低，提高了運(yùn)算速度；（2)結(jié)合粗基線(xiàn)和一階數(shù)值微分方法，扣除疊加在直 流基線(xiàn)上的異常信號(hào)(包括脈沖和干擾），可以求得更加精密的基線(xiàn)，不受溫度等干擾因素 的影響，能夠?qū)崟r(shí)準(zhǔn)確跟蹤基線(xiàn)的變化。
【附圖說(shuō)明】
[0027] 圖1為本發(fā)明實(shí)施例1的系統(tǒng)結(jié)構(gòu)框圖；
[0028] 圖2為本發(fā)明實(shí)施例1的不同含量的鈾礦模型井中的實(shí)測(cè)基線(xiàn)漂移對(duì)比圖；
[0029] 圖3為本發(fā)明實(shí)施例1的不同溫度條件下的實(shí)測(cè)基線(xiàn)漂移對(duì)比圖。
[0030] 圖中：核輻射信號(hào)探測(cè)器-I，信號(hào)調(diào)理電路-2，高速模數(shù)轉(zhuǎn)換器-3，F(xiàn)PGA芯片-4， FPGA控制系統(tǒng)進(jìn)行高速模數(shù)采樣-5,求得粗基線(xiàn)-6,對(duì)采樣數(shù)據(jù)進(jìn)行一階數(shù)值微分-7,求得 精密基線(xiàn)-8,基線(xiàn)扣除_9,求得能譜曲線(xiàn)-10。
【具體實(shí)施方式】
[0031] 下面通過(guò)附圖和【具體實(shí)施方式】對(duì)本發(fā)明作更為詳細(xì)的描述。
[0032] 本發(fā)明的基本思想是在以FPGA為核心的數(shù)字化核能譜測(cè)量系統(tǒng)中（如圖1所示）， 利用FPGA芯片控制高速模數(shù)轉(zhuǎn)換器進(jìn)行數(shù)據(jù)連續(xù)采樣，先對(duì)采樣數(shù)據(jù)使用滑動(dòng)平均濾波法 求得粗基線(xiàn)，然后結(jié)合粗基線(xiàn)和一階數(shù)值微分法判斷并扣除采樣數(shù)據(jù)中的異常信號(hào)（異常 信號(hào)是指疊加在直流基線(xiàn)上的脈沖或干擾信號(hào)），再對(duì)扣除異常信號(hào)后的數(shù)據(jù)使用滑動(dòng)平 均濾波法求得精密基線(xiàn)。從而在FPGA芯片上實(shí)現(xiàn)數(shù)字化核能譜測(cè)量系統(tǒng)中的高精度脈沖基 線(xiàn)估計(jì)。
[0033] 本發(fā)明所述的在FPGA中進(jìn)行高精度脈沖基線(xiàn)估計(jì)的具體步驟如下：
[0034] a、初始化，在FPGA芯片內(nèi)開(kāi)辟兩組大小分別為N = 2$PM=2m個(gè)單元的循環(huán)緩沖區(qū)A 和B，在執(zhí)行平均濾波之前，將緩沖區(qū)的內(nèi)容全部初始化為0，同時(shí)將SumPre和SumPr^初始 化為〇;
[0035] b、由FPGA控制系統(tǒng)按一定的時(shí)鐘頻率進(jìn)行模數(shù)采樣，采樣數(shù)據(jù)存入循環(huán)緩沖區(qū)A 的頭指針處；
[0036] c、對(duì)循環(huán)緩沖區(qū)A中的數(shù)據(jù)使用滑動(dòng)平均濾波法求得粗基線(xiàn)；
[0037] d、利用一階數(shù)值微分法判斷采樣數(shù)據(jù)是否為脈沖或干擾信號(hào)；
[0038] e、如采樣數(shù)據(jù)不是脈沖或干擾信號(hào)，進(jìn)入步驟f，否則進(jìn)入步驟g;
[0039] f、將該采樣數(shù)據(jù)存入循環(huán)緩沖區(qū)B的頭指針處；
[0040] g、對(duì)循環(huán)緩沖區(qū)B中的數(shù)據(jù)使用滑動(dòng)平均濾波法求得精密基線(xiàn)；
[0041 ] h、返回步驟b。
[0042] 通常情況下，若要在FPGA中求取N個(gè)采樣數(shù)據(jù)的平均值J，可通過(guò)如下的基本原理 公式實(shí)
[0043]
[0044] 其中x(n-i)表示當(dāng)前采樣數(shù)據(jù)x(n)向后延遲i個(gè)單位(i個(gè)時(shí)鐘序列）后得到的值。 該方法需要在FPGA中占用N個(gè)寄存器，且需要進(jìn)行N-I次加法運(yùn)算和一次除法運(yùn)算(如果取N 為2的整數(shù)次冪即N=2k，則可以將二進(jìn)制除法轉(zhuǎn)化為移位來(lái)實(shí)現(xiàn)），當(dāng)N取值較大時(shí)，運(yùn)算較 為復(fù)雜，該方法難以在FPGA上予以實(shí)現(xiàn)。
[0045] 本發(fā)明提出的在FPGA中引入循環(huán)緩沖區(qū)來(lái)滑動(dòng)平均濾波，可最大程度降低運(yùn)算復(fù) 雜度，且運(yùn)算次數(shù)不受N的影響。在FPGA中開(kāi)辟一組N個(gè)單元的循環(huán)緩沖區(qū)，在執(zhí)行平均濾波 之前，將緩沖區(qū)的內(nèi)容全部初始化為〇,每次新采樣得到的數(shù)據(jù)x(n)需要存入緩沖區(qū)頭指針 BUF_IN指向的位置，尾指針BUF_0UT指向的x(n-N-l)為循環(huán)緩沖區(qū)中最老的一個(gè)數(shù)據(jù)，需要 予以剔除。利用如下公式求得N個(gè)采樣數(shù)據(jù)的平均值無(wú)：
[0046]
[0047]
[0048] 式中Sum為循環(huán)緩沖區(qū)中N個(gè)數(shù)據(jù)的累加和，SumPre為Sum延遲一個(gè)時(shí)鐘序列的值， SumPr e的初始值為0，Sum >> k表示將Sum右移k個(gè)二進(jìn)制位。該方法只需要進(jìn)行一次加法、 一次減法和一次移位操作即可實(shí)現(xiàn)N個(gè)數(shù)求平均的效果?？梢钥闯隼帽景l(fā)明提出的方法 實(shí)現(xiàn)基線(xiàn)估計(jì)的算法復(fù)雜度大大降低。
[0049] 本發(fā)明所述的精密基線(xiàn)的準(zhǔn)確估計(jì)時(shí)間不僅與粗基線(xiàn)計(jì)算的緩沖區(qū)長(zhǎng)度N、精密 基線(xiàn)計(jì)算的緩沖區(qū)長(zhǎng)度M及系統(tǒng)采樣周期Tsample(或采樣頻率f sample)有關(guān)，同時(shí)還與系統(tǒng)的 脈沖通過(guò)率有關(guān)，通常情況下系統(tǒng)默認(rèn)的可靠工作脈沖通過(guò)率是系統(tǒng)最高極限脈沖計(jì)數(shù)率 的0.1倍。若想求得準(zhǔn)確的粗基線(xiàn)，需要在循環(huán)緩沖區(qū)N中填充的所有數(shù)據(jù)均為有效的采樣 數(shù)據(jù);若想求得準(zhǔn)確的精密基線(xiàn)，需要在循環(huán)緩沖區(qū)M中填充的所有數(shù)據(jù)均為扣除異常(脈 沖和干擾)后的采樣數(shù)據(jù)，因?yàn)橥ǔＧ闆r下每M個(gè)采樣數(shù)據(jù)中最多包含0.1 XM個(gè)脈沖數(shù)據(jù)， 即需要M/(l-0.1)個(gè)采樣數(shù)據(jù)才能填滿(mǎn)循環(huán)緩沖區(qū)M?；€(xiàn)估計(jì)的最長(zhǎng)時(shí)間可通過(guò)下式求 得：
[0050] M
[0051] 式中NXTsample3表示粗基線(xiàn)的最長(zhǎng)準(zhǔn)確計(jì)算時(shí)間，.表示精密基線(xiàn)的最 1U, I 長(zhǎng)準(zhǔn)確計(jì)算時(shí)間，若~(yú) = 256、]\1=1024、1^_16=16.7118(即;1^_16 = 601抱），則1:/~23118。可以 看出利用本發(fā)明提出的方法實(shí)現(xiàn)基線(xiàn)估計(jì)的實(shí)時(shí)性非常好。
[0052]將本方法嵌入到數(shù)字化核能譜測(cè)量系統(tǒng)中，在常溫下于不同時(shí)間段對(duì)東華理工大 學(xué)的7個(gè)標(biāo)準(zhǔn)鈾礦模型井進(jìn)行了測(cè)試，產(chǎn)生的脈沖基線(xiàn)估計(jì)效果如圖2所示。從圖2可看出， 使用本方法前，儀器輸入脈沖的基線(xiàn)在_8mV左右波動(dòng);使用本方法后，基線(xiàn)在OmV左右波動(dòng)， 相應(yīng)的能譜測(cè)量過(guò)程中的能量刻度曲線(xiàn)均可看成經(jīng)過(guò)坐標(biāo)零點(diǎn)，基線(xiàn)漂移情況基本得到校 正。
[0053]將本方法嵌入到數(shù)字化核能譜測(cè)量系統(tǒng)中，在實(shí)驗(yàn)室的恒溫箱中對(duì)鈾礦石粉末進(jìn) 行了測(cè)試，在(TC~60°C的溫度變化條件下，產(chǎn)生的脈沖基線(xiàn)估計(jì)效果如圖3所示。從圖3可 看出，使用本方法前，儀器的脈沖直流基線(xiàn)隨溫度變化而變化;使用本方法后，基線(xiàn)穩(wěn)定在 2.5mV左右，在誤差允許范圍內(nèi)，相應(yīng)的能譜測(cè)量過(guò)程中的能量刻度曲線(xiàn)均可看成經(jīng)過(guò)坐標(biāo) 零點(diǎn)，溫度引起的基線(xiàn)漂移情況基本得到校正。
【主權(quán)項(xiàng)】
1. 一種數(shù)字化核能譜測(cè)量系統(tǒng)中的脈沖基線(xiàn)估計(jì)方法，其步驟： 首先在FPGA忍片內(nèi)求得粗基線(xiàn)， (1) 在FPGA忍片內(nèi)開(kāi)辟一組大小為N = 2k個(gè)單元的循環(huán)緩沖區(qū)A，在執(zhí)行平均濾波之前， 將緩沖區(qū)的內(nèi)容全部初始化為0; (2) 每次采樣得到的數(shù)據(jù)x(n)存入循環(huán)緩沖區(qū)A頭指針指向的位置，剔除循環(huán)緩沖區(qū)A 尾指針指向的數(shù)據(jù)x(n-N-l)，并利用如下公式求得粗基線(xiàn)定: Sum = SumPre+x(η)-X(n-N-1) 5. Sum ! Ν 二 Sum / 二 Slim ? k 式中Sum為循環(huán)緩沖區(qū)A中N個(gè)數(shù)據(jù)的累加和，SumPre為Sum延遲一個(gè)時(shí)鐘序列的值， SumPre的初始值為0，S皿> >k表示將S皿右移k個(gè)二進(jìn)制位； 其次根據(jù)粗基線(xiàn)扣除采樣數(shù)據(jù)中的異常信號(hào)： (3) 將采樣數(shù)據(jù)與粗基線(xiàn)進(jìn)行比較，如果大于粗基線(xiàn)，則認(rèn)為該數(shù)據(jù)是疊加在基線(xiàn)上的 異常信號(hào)；如果小于或等于粗基線(xiàn)（即.命0心叫，則對(duì)其進(jìn)行一階數(shù)值微分，微分結(jié)果用dx (η)表示，公式為： dx(n) = χ(η)-χ(η-1) 式中χ(η)表示當(dāng)次采樣數(shù)據(jù)，x(n-l)表示前一次采樣數(shù)據(jù)； (4) 選取P個(gè)數(shù)據(jù)來(lái)判斷輸入數(shù)據(jù)進(jìn)行一階數(shù)值微分后，其數(shù)據(jù)變化幅度大小及正負(fù)符 號(hào)是否連續(xù)，如果微分后的數(shù)據(jù)變化幅度較小且正負(fù)符號(hào)不連續(xù)，則認(rèn)為該數(shù)據(jù)為基線(xiàn)，通 過(guò)下列公式求得扣除異常信號(hào)后的基線(xiàn)數(shù)據(jù)χ/ (η):式中玄為已經(jīng)求得的粗基線(xiàn)值為連續(xù)Ρ個(gè)數(shù)據(jù)一階微分后的數(shù)值求和，|dx (η) I表示數(shù)值微分dx(n)的絕對(duì)值，ε表示輸入數(shù)據(jù)的紋波范圍； 最后根據(jù)扣除異常信號(hào)后的數(shù)據(jù)求得精密基線(xiàn)： (5) 在FPGA忍片內(nèi)開(kāi)辟一組大小為Μ = 2"個(gè)單元的循環(huán)緩沖區(qū)Β，在執(zhí)行平均濾波之前， 將緩沖區(qū)的內(nèi)容全部初始化為0; (6) 將扣除異常信號(hào)后的基線(xiàn)數(shù)據(jù)χ/(η)存入循環(huán)緩沖區(qū)Β頭指針指向的位置，剔除循 環(huán)緩沖區(qū)Β尾指針指向的數(shù)據(jù)χ/ (η-Ν-1)，并利用如下公式求得精密基線(xiàn)呆'：式中Sun/為循環(huán)緩沖區(qū)Β中Μ個(gè)數(shù)據(jù)的累加和，SumPre^為Sun/延遲一個(gè)時(shí)鐘序列的值， SumPre^的初始值為0，Sun/ >>m表示將Sun/右移m個(gè)二進(jìn)制位。
【文檔編號(hào)】G01T1/36GK105842728SQ201610171727
【公開(kāi)日】2016年8月10日
【申請(qǐng)日】2016年3月24日
【發(fā)明人】王海濤, 湯彬, 王仁波, 瞿金輝, 劉志鋒, 張雄杰, 周書(shū)民, 陳銳
【申請(qǐng)人】東華理工大學(xué)