單顯明,劉業(yè)峰,那崇正,靳 新
(1.沈陽工學院 信息與控制學院,遼寧 撫順 113122;2.遼寧省數(shù)控機床信息物理融合與智能制造重點實驗室,遼寧 撫順 113122;3.沈陽工學院 機械工程與自動化學院,遼寧 撫順 113122)
球磨機是磨礦工作環(huán)節(jié)中的重要設(shè)備,它的作用是通過轉(zhuǎn)動來推動球體對礦石進行沖擊和粉碎,從而將粉碎的礦石磨成顆粒大小適宜的礦漿,在水泥、硅酸鹽制品、化肥、黑色和有色金屬的選礦、玻璃、陶瓷等行業(yè)中都有著廣泛的應用。球磨機的符合參數(shù)與磨礦質(zhì)量、磨礦效率等因素密切相關(guān)。然而球磨機在實際磨礦工作過程中,由于磨機筒體內(nèi)部環(huán)境惡劣且工況復雜,其負荷參數(shù)難以準確檢測和控制。
為了提高磨礦工業(yè)過程的生產(chǎn)率和生產(chǎn)質(zhì)量,有必要對球磨機在不同工況下的負荷參數(shù)進行測量,為球磨機的控制提供參考數(shù)據(jù)。軟測量是球磨機負荷參數(shù)常見的測量方式,其測量原理是將自動化技術(shù)引入到工業(yè)生產(chǎn)中,利用狀態(tài)估算技術(shù)在線估算難以直接檢測的參數(shù),從而達到軟件替代硬件的效果。一般情況下,對于不易測量的變量,或無法直接測量的變量,可以通過獲取易于測量的變量,并建立易測變量和難測變量的關(guān)系,從而實現(xiàn)在線測量。
現(xiàn)階段支恩瑋[1]等設(shè)計了基于域適應支持向量回歸的負荷軟測量系統(tǒng),對多工況下的數(shù)據(jù)進行預處理,并對其頻譜特征進行提取。在此基礎(chǔ)上,利用目標區(qū)域內(nèi)的小樣本所包含的特征和知識結(jié)構(gòu),提高源域數(shù)據(jù)建立模式與目標領(lǐng)域數(shù)據(jù)相匹配的能力。采用回歸方法預測了磨機的載荷參數(shù)。該系統(tǒng)具有一定的有效性,然而該系統(tǒng)在實際應用過程,其測量功能和運行性能較低。
針對上述問題,設(shè)計了基于ASOS-ELM的球磨機負荷參數(shù)軟測量系統(tǒng)。ASOS算法指的是共生生物搜索算法,通過模仿生物行為,引導個體逐漸進化。而ELM算法是指極限學習機算法,適用于監(jiān)督學習和非監(jiān)督學習問題。通過ASOS-ELM算法的聯(lián)合應用,以期能夠提升球磨機負荷參數(shù)軟測量系統(tǒng)的測量功能與運行性能,以此提高磨礦質(zhì)量與生產(chǎn)效率。
由于球磨機在負載工況下會產(chǎn)生振動信號,根據(jù)振動信號的特征可以推導出球磨機的工作負荷參數(shù),因此需要在硬件系統(tǒng)環(huán)境中安裝球磨機工作數(shù)據(jù)采集器,主要采集對象為球磨機的振動工作數(shù)據(jù)。球磨機工作數(shù)據(jù)采集器的基本結(jié)構(gòu)如圖1所示。
圖1 球磨機工作數(shù)據(jù)采集器結(jié)構(gòu)圖
圖1中,DH5922N型號動態(tài)數(shù)據(jù)采集芯片能夠?qū)Α?、振動、沖擊等各種物理參數(shù)進行檢測和分析[2]。DH131型加速度傳感器的輸出為電荷信號,它必須用相應的DH5857-1型電荷適配器來進行數(shù)據(jù)類型轉(zhuǎn)換[3]。將電荷轉(zhuǎn)換成電壓信號后,將其送入動態(tài)數(shù)據(jù)采集芯片中。安裝的球磨機工作數(shù)據(jù)采集器的頻率頻響范圍為[1 Hz,10 000 Hz],而球磨機的頻率在8 686 Hz以下,所以設(shè)計的數(shù)據(jù)采集器能夠滿足對球磨機信號的采集要求。除了振動數(shù)據(jù)采集裝置外,為保證系統(tǒng)處理器設(shè)備能夠直接處理采集信號,在采集裝置輸出前端分別嵌入前置放大器和抗混疊濾波器,前置放大器的作用是對采集器輸出結(jié)果中的阻抗進行表換處理,選擇了一種用于增益FET的電壓互補源跟蹤電路,它可以在FET中增加相同的源間電壓,從而提高了輸入阻抗,降低輸出阻抗的同時降低噪聲[4]。在采集信號時,根據(jù)取樣原理,若有一種干擾信號,其頻率超過取樣頻率的一半,則會出現(xiàn)一種無法用數(shù)字濾波方法消除的混疊信號。因此,在取樣前連接MAX293抗混雜濾波器,用于抗混疊濾波。
系統(tǒng)處理器的處理任務包括球磨機初始振動數(shù)據(jù)的處理、負荷參數(shù)軟測量程序的運行以及ASOS-ELM算法的運行,安裝的系統(tǒng)處理器型號為OMAP-L138。該處理器為達芬奇架構(gòu)的嵌入式處理器,在處理器內(nèi)部同時嵌入DSP內(nèi)核和ARM內(nèi)核[5]。OMAP-L138處理器既具備了復雜的指令集,又具備了指令集的簡化技術(shù)。因此,OMAP-L138相對于傳統(tǒng)系統(tǒng)中使用的單核結(jié)構(gòu)處理器有更好的性能。OMAP-L138包含一個ARM9內(nèi)核和一個C674x的DSP內(nèi)核,這兩種內(nèi)核都可以在456 MHz的頻率下工作,所以O(shè)MAP-L138不僅具有很好的控制性能,而且還具有很好的計算能力[6]。另外,OMAP-L138還整合了外部設(shè)備,包括增強型直接數(shù)據(jù)接入、多通道緩沖串口、通用串行通信接口等。
為了給球磨機負荷數(shù)據(jù)以及軟測量程序提供充足的存儲空間,在傳統(tǒng)硬件系統(tǒng)存儲器的基礎(chǔ)上進行存儲擴展,擴展結(jié)果如圖2所示。
圖2 擴展存儲器結(jié)構(gòu)圖
設(shè)計的存儲器充分利用RAM的片內(nèi)資源,同時保證所有的地址不發(fā)生沖突。512KBRAM芯片的接線將A15僅被用作74138的控制信號,與P10直接連接[7]。由于A15連接到P10,故芯片的物理空間充分利用,所以在CPU運行時,所有的邏輯地址都是8000H-FFFFFH。
顯示屏的安裝為球磨機負荷參數(shù)軟測量結(jié)果的輸出提供了可視化硬件支持,所選用的顯示屏設(shè)備是DC48270N043_0lTW/NW,該設(shè)備自帶ASCII、GBK和GB2312字庫,速度快、成本低,適合圖片少畫面不復雜的場合,其能最多存儲27張全屏圖片,需要的供電電壓為5 V,其通訊接口為RS232,易于與系統(tǒng)電路主板相連接,顯示和控制符合系統(tǒng)需求。
在基于ASOS-ELM的球磨機負荷參數(shù)軟測量系統(tǒng)中,供電質(zhì)量將直接影響到系統(tǒng)的功能與性能的運行效果。硬件系統(tǒng)中各個組成元件所需電壓規(guī)格包括5 V、3.3 V、1.8 V、1.0 V,結(jié)合系統(tǒng)各個部分功耗以及對電流的要求,外部供電電壓為5 V,XC7Z010的COREPOWER為18 V和10 V,IOPOWER為33 V[8]。所以,系統(tǒng)選擇了TPS70151為電源芯片。電源電路的供電次序與處理器上的斷電次序一致,并能滿足系統(tǒng)所需的電流和功率[9]。除了上述系統(tǒng)元件外,通信網(wǎng)絡(luò)、程序主控芯片等部分均沿用傳統(tǒng)軟測量系統(tǒng),利用設(shè)計的系統(tǒng)電源電路,實現(xiàn)所有硬件設(shè)備的連接。
在硬件設(shè)備的支持下,通過ASOS-ELM算法的應用,實現(xiàn)球磨機負荷參數(shù)軟測量功能。根據(jù)不同負荷下球磨機振動信號的變化特征,實現(xiàn)實時負荷參數(shù)的軟測量,并將測量結(jié)果以可視化的形式輸出。
從結(jié)構(gòu)方面來看,球磨機包括進料部、出料部、回轉(zhuǎn)部、傳動部等多個部分[10]。中空軸以鑄鋼作為材料,筒體內(nèi)嵌耐磨襯墊,用來提高設(shè)備的耐磨性。筒體是球磨機的主要加工位置,筒體上裝有支撐圓筒和支撐圓筒轉(zhuǎn)動的支座,并利用電動機進行設(shè)備驅(qū)動。物料由進料機構(gòu)經(jīng)料中空心軸,以螺旋狀的方式均勻地送入磨機的第一倉,在料倉中裝有臺階墊片或波浪形墊片,里面裝有各種尺寸的鋼球,利用滾筒旋轉(zhuǎn)產(chǎn)生的離心力完成物料的研磨操作[11]。物料完成粗磨工作后,進入第二倉作進一步研磨。滾筒在旋轉(zhuǎn)的過程中,研磨體自身也會發(fā)生滑動,在滾筒的滑動過程中,物料受到磨料的影響,為了充分利用研磨作用,當物料力度較大時,將磨料圓筒分成兩個部分,當材料進入到第一倉時在鋼球作用下被擊碎,而當材料進入第二倉時,鋼段會對物料進行進一步碾壓,將合格的物料從出料端的中空軸排出[12]。將上述球磨機的工作原理代入到球磨機結(jié)構(gòu)中,得出球磨機數(shù)學模型的構(gòu)建結(jié)果,如圖3所示。
圖3 球磨機數(shù)學模型示意圖
球磨機數(shù)學模型中破碎環(huán)節(jié)的數(shù)學表達式如下:
gdischarge=(B·S+I-S)nfracture·ggive
(1)
公式(1)中,ggive和gdischarge分別為球磨機給進礦和排出礦的粒度,B和S分別對應的是球磨機的破碎函數(shù)和選擇函數(shù),I表示單元矩陣,nfracture為物料的破碎段數(shù)[13]。同理可以得出球磨機多個加工環(huán)節(jié)的數(shù)學表達式,將其與球磨機結(jié)構(gòu)融合在一起,得出球磨機數(shù)學模型的最終構(gòu)建結(jié)果。
由于球磨機工作負荷與振動信號之間存在相關(guān)關(guān)系,因此可以利用振動信號來反映球磨機的負荷參數(shù),為此利用硬件系統(tǒng)中的采集器設(shè)備對球磨機工作狀態(tài)下的振動信號進行采集[14]。在忽略筒體轉(zhuǎn)動效應的情況下,依據(jù)動量定理,可以得出球載下落對筒體的沖擊力為:
(2)
公式(2)中,Lmill為球磨機筒體長度,ω和θ分別為筒體轉(zhuǎn)動角速度和轉(zhuǎn)動速度,θ1和θ2為轉(zhuǎn)動速度分量,其計算公式如下:
(3)
公式(3)中,ν為磨機轉(zhuǎn)速,r1和r2分別為球載回轉(zhuǎn)半徑的最大值和最小值。另外,公式(2)中變量ρ表示的是球載松散密度。在沖擊力作用下,球磨機發(fā)生振動。由于球磨機筒體呈對稱式分布,采用不連續(xù)方式進行進出料,所以在檢測圓筒的振動時,應選擇一側(cè)軸承為測量點在筒體旋轉(zhuǎn)時,鋼珠和物料的拋出對筒體的沖擊會從筒體進出料端蓋傳遞到軸承座,同時,由于內(nèi)部負荷的變化,物料與磨介之間的相互作用使其產(chǎn)生振動[15-16]?;诹ζ胶庠恚梢源_定球磨機中振動幅度最大的位置,以此作為測點,并將采集器設(shè)備安裝在該測點上,通過采集器的運行與輸出,得出球磨機實時振動信號的采集結(jié)果[17]。為了提高球磨機振動信號的采集質(zhì)量,同時去除信號中的正、負兩種噪聲,利用內(nèi)置的濾波器進行信號預處理,濾波器的數(shù)學表達式如下:
(4)
公式(4)中,f(x)和g(x)分別為初始采集振動信號和轉(zhuǎn)換振動信號,HOC(f(x))和HCO(f(x))對應的是開-閉濾波器和閉-開濾波器,最終得出的球磨機振動信號去噪處理結(jié)果為:
(5)
分別將實時振動采集信號和兩個濾波函數(shù)代入到公式(5)中,便可得出滿足質(zhì)量要求的振動信號采集結(jié)果[18]。
采用ELM算法分別建立振動信號頻譜特征和磨機負荷參數(shù)之間的非線性模型,ELM結(jié)構(gòu)如圖4所示。
圖4 極限學習機結(jié)構(gòu)圖
假設(shè)極限學習機的輸入層、隱含層和輸出層的節(jié)點數(shù)量分別為n、p和q,對于輸入的球磨機振動信號,若極限學習機能以極低的誤差逼近振動信號樣本,則公式(6)成立:
(6)
公式(6)中,yj為預處理完成的球磨機振動信號,?i為隱含層第i個節(jié)點到輸出層的權(quán)重值,H(ai,bi,yj)為激活函數(shù),其函數(shù)表達式如下:
(7)
公式(7)中,ai和bi分別表示的是連接第i個隱含層結(jié)點的輸入權(quán)值和輸出權(quán)值[19]。為了提高系統(tǒng)的測量精度,利用ASOS算法對ELM工作系數(shù)進行優(yōu)化,優(yōu)化過程如圖5所示。
圖5 ASOS算法優(yōu)化ELM的流程圖
在ASOS算法運行過程中,利用公式(8)對“生物”種群進行初始化。
Xi=Wb+rand(1,D)(Ub-Wb)
(8)
公式(8)中,Ub和Wb分別為共生生物搜索空間的上界和下界,rand(1,D)為縮放因子向量,其中D為初始解的維數(shù)[20]。搜索種群中的最優(yōu)值Xbest,通過互利和共棲階段完成種群更新,更新過程可以表示為:
(9)
公式(9)中,Mv和κ1分別為互利向量和隨機因子。在寄生階段產(chǎn)生“寄生向量”,當達到終止條件時,終止優(yōu)化操作,否則重新執(zhí)行上述操作,直到滿足終止條件為止[21]。將采集并處理完成的球磨機振動信號代入到ASOS優(yōu)化的ELM中,得出特征提取結(jié)果,其中峰值特征的提取結(jié)果可以表示為:
τm-m=xmax(t)-xmin(t)
(10)
公式(10)中,xmin(t)和xmax(t)分別為球磨機振動信號的最小值和最大值。同理可以得出球磨機所有振動特征的提取結(jié)果。
根據(jù)球磨機振動信號特征提取結(jié)果,選擇待軟測量的球磨機負荷參數(shù)為球負荷、料負荷、水負荷,分別記為φb、φm和φw,上述負荷參數(shù)的計算公式如下:
(11)
式(11)中,變量φbmw、φmw、φmb以及φbf分別為充填率、磨礦濃度、介質(zhì)充填率以及料球比,Vmill表示球磨機容積,ρm、ρw和ρb對應的是物料、水和研磨介質(zhì)的密度值,δ為介質(zhì)空隙率[22]。因此可以通過測量料球比、磨礦濃度和充填率計算得出球磨機的負荷參數(shù),上述參數(shù)即為軟測量的輔助變量。最終在顯示屏的支持下,得出球磨機負荷參數(shù)軟測量結(jié)果的可視化輸出結(jié)果。
為了驗證設(shè)計的基于ASOS-ELM的球磨機負荷參數(shù)軟測量系統(tǒng)的有效性,融合硬件設(shè)備和軟件程序,模擬系統(tǒng)的實際應用需要和環(huán)境,檢驗設(shè)計的基于ASOS-ELM的球磨機負荷參數(shù)軟測量系統(tǒng)是否能夠達到預期效果。此次系統(tǒng)測試實驗分為兩個部分,分別為軟測量功能測試和運行性能測試,并通過與基于域適應支持向量回歸的負荷軟測量系統(tǒng)的對比,驗證設(shè)計的基于ASOS-ELM的球磨機負荷參數(shù)軟測量系統(tǒng)的應用優(yōu)勢。
此次實驗以某大型采礦工程作為實驗背景,選擇礦石的加工環(huán)境作為系統(tǒng)的開發(fā)與測試環(huán)境,在實驗環(huán)境中安裝主測計算機作為系統(tǒng)的開發(fā)設(shè)備。按照硬件系統(tǒng)的設(shè)計結(jié)果,分別安裝采集器、處理器等設(shè)備,并利用供電電源電路完成硬件設(shè)備的連接。在此基礎(chǔ)上,設(shè)置通信網(wǎng)絡(luò)路由器和網(wǎng)關(guān)的工作參數(shù),保證配置的網(wǎng)絡(luò)環(huán)境能夠覆蓋整個實驗環(huán)境,確保球磨機的實時工作參數(shù)能夠被正常采集與傳輸。由于球磨機屬于特種設(shè)備,為完成負荷參數(shù)的軟測量工作的安全性,需在實驗環(huán)境中加設(shè)防護網(wǎng)。系統(tǒng)測試環(huán)境中所有設(shè)備安裝完成后,需要對其進行聯(lián)合調(diào)試,若存在設(shè)備運行故障,則需要重新選擇系統(tǒng)運行設(shè)備,直至系統(tǒng)測試環(huán)境調(diào)試成功。
為滿足礦石的處理需求,選擇型號為Φ450 mm×450 mm的Bond指數(shù)球磨機作為系統(tǒng)測試實驗的研究對象,該設(shè)備的電機功率為0.8 kW。球磨機的圓筒在中線的兩端設(shè)有端蓋,用于填料、排料,并通過軸承的整體放置在機架上,在軸承的一端通過傳動裝置與蝸輪蝸桿減速裝置相連接,并由電機驅(qū)動。重復執(zhí)行“啟動-停止”操作,判斷準備的球磨機是否能夠在實驗環(huán)境中正常運行。
為了保證設(shè)計的基于ASOS-ELM的球磨機負荷參數(shù)軟測量系統(tǒng)適用于球磨機的多種工作狀況,此次系統(tǒng)測試實驗分別設(shè)置空磨、正常以及飽磨三種工況,空磨指的是球磨機圓筒內(nèi)不存在需要加工的礦石,正常工況是指球磨機圓筒內(nèi)的加工礦石含量占圓筒總量的10%~75%,當球磨機圓筒中加入的礦石量高于圓筒總?cè)萘康?5%時,認為當前球磨機處于飽磨工況。實驗準備的加工礦石物料為某礦山的鎢礦石,其普式硬度系數(shù)約為16,密度約為1 950 kg/m3。為了避免由于入料的粒度分布不均勻,影響球磨機運行工況的控制效果,要求在對球磨機進行入料前,對碎磨前的礦石物料進行初級破碎,并采用電動篩分機對礦石進行篩分,通過控制物料的輸入量,確定料球比、磨礦濃度和充填率的真實值,并求解出球磨機負荷參數(shù)的實際值,實際負荷參數(shù)的設(shè)置情況如表1所示。
表1 球磨機負荷參數(shù)設(shè)置表
以表1中的設(shè)置數(shù)據(jù)作為判斷系統(tǒng)軟測量功能的對比標準數(shù)據(jù)。
由于此次實驗應用了ASOS和ELM算法,需要對相關(guān)工作系數(shù)進行設(shè)置,ASOS算法中種群個數(shù)設(shè)置為50,最大迭代次數(shù)為100次。而ELM算法的隱藏層節(jié)點數(shù)量設(shè)置為20,學習率為0.001%。
利用編碼工具將基于ASOS-ELM的球磨機負荷參數(shù)軟測量系統(tǒng)的軟件功能設(shè)計結(jié)果轉(zhuǎn)換為程序代碼,分別在三種不同的球磨機運行工況下,利用采集器設(shè)備獲得振動信號數(shù)據(jù)作為實驗數(shù)據(jù),以振動信號作為輸入數(shù)據(jù),球磨機負荷參數(shù)作為輸出參數(shù),通過硬件設(shè)備與軟件程序的協(xié)同運行,得出系統(tǒng)運行結(jié)果。三種工況下球磨機負荷參數(shù)的軟測試輸出結(jié)果如圖6所示。
圖6 球磨機負荷參數(shù)軟測量系統(tǒng)運行結(jié)果
在三種工況下進行多次實驗,采用求解平均值的方式得出更為精準的系統(tǒng)測試結(jié)果。此外,為了體現(xiàn)出設(shè)計系統(tǒng)的運行優(yōu)勢,設(shè)置基于域適應支持向量回歸的負荷軟測量系統(tǒng)作為實驗的對比系統(tǒng),對比系統(tǒng)的運行環(huán)境以及測量對象均相同。
此次實驗分為系統(tǒng)功能測試和系統(tǒng)性能測試兩個部分,其中功能量化測試指標設(shè)置為負荷參數(shù)軟測量的均方誤差,其數(shù)值結(jié)果為:
(12)
Δt=tout-tin
(13)
公式(13)中,tin和tout分別表示球磨機振動數(shù)據(jù)的輸入時間以及球磨機負荷參數(shù)軟測量結(jié)果的輸出時間,計算得出系統(tǒng)程序的運行時間更短,說明對應系統(tǒng)的運行性能越優(yōu)。
為了驗證設(shè)計系統(tǒng)的整體有效性,利用ASOS-ELM算法和域適應支持向量回歸方法對球磨機負荷參數(shù)分別進行測試。選取球負荷、料負荷和水負荷為試驗參數(shù),在Bond指數(shù)球磨機筒體外側(cè)安裝傳感器,保證球磨機轉(zhuǎn)速穩(wěn)定后,依次加入一定量的鋼球、物料和水,在球磨機運行一段時間后,記錄測試結(jié)果,空磨工況下系統(tǒng)軟測量功能的測試結(jié)果如表2所示。
表2 空磨工況下軟測量系統(tǒng)功能測試結(jié)果
根據(jù)表2結(jié)果可知,在5次實驗中,對比系統(tǒng)的球負荷在393.8~399.2 kg之間,料負荷在0.2~0.9 kg之間,水負荷在44.1~53.2 kg之間;而設(shè)計系統(tǒng)的球負荷在399.4~400 kg之間,料負荷在0~0.5 kg之間,水負荷在45~54.7 kg之間,與對比系統(tǒng)相比,設(shè)計系統(tǒng)的整體球負荷和水負荷較高,料負荷較低。將表2與表1中的數(shù)據(jù)代入到公式(12)中,得出對比系統(tǒng)的球負荷、料負荷、水負荷參數(shù)的平均軟測量誤差分別為0.42 kg、0.007 kg和0.12 kg。而設(shè)計系統(tǒng)的球負荷、料負荷、水負荷參數(shù)軟測量誤差的平均值分別為0.002 kg、0.001 kg和0.001 kg。同理在正常工況和飽磨工況下,通過公式(12)的計算得出正常和飽磨工況下系統(tǒng)的軟測試功能測試對比結(jié)果,如圖7所示。
圖7 正常和飽磨工況下系統(tǒng)的軟測試功能測試對比結(jié)果
綜上所述,在球磨機的三種工況下,基于ASOS-ELM的球磨機負荷參數(shù)軟測量系統(tǒng)的測量誤差始終低于對比系統(tǒng),即設(shè)計系統(tǒng)在測量功能方面具有更高優(yōu)勢。
通過主測計算機后臺數(shù)據(jù)的提取以及公式(13)的計算,得出系統(tǒng)運行時間的測試結(jié)果,如表3所示。
表3 球磨機負荷參數(shù)軟測量系統(tǒng)運行性能測試數(shù)據(jù)表
經(jīng)過表3對比可以看出,球磨機在空磨工況下的負荷軟測量速度更快。在三種工況下,對比系統(tǒng)的運行時間始終高于設(shè)計系統(tǒng),通過平均值計算,兩種系統(tǒng)的平均運行時間分別為7.78 s和3.65 s,即設(shè)計的基于ASOS-ELM的球磨機負荷參數(shù)軟測量系統(tǒng)的運行性能更優(yōu)。
設(shè)計的基于ASOS-ELM的球磨機負荷參數(shù)軟測量系統(tǒng),通過ASOS-ELM算法的應用和軟測量系統(tǒng)的設(shè)計與開發(fā),保證了球磨機負荷參數(shù)測量的準確度,有效提升了系統(tǒng)運行性能。然而系統(tǒng)測試實驗中使用的球磨機型號較為單一,因此設(shè)計的軟測量系統(tǒng)可能會出現(xiàn)與部分球磨機設(shè)備不適配的問題,針對這一問題還需要在今后的研究工作中進一步解決。