摘" 要: 針對(duì)串聯(lián)鋰電池組均衡過程中出現(xiàn)的均衡速度慢、能量轉(zhuǎn)移效率低等問題,提出一種基于Buck?Boost電路的分層均衡拓?fù)浣Y(jié)構(gòu),以各單體電池的荷電狀態(tài)(SOC)作為均衡變量,采用相對(duì)極差法的均衡控制策略以實(shí)現(xiàn)各層均衡單元的同步均衡。首先闡述了改進(jìn)后均衡拓?fù)涞恼w結(jié)構(gòu)以及分層分組的方法、最小均衡單元的均衡原理以及均衡過程中相關(guān)參數(shù)的計(jì)算過程;之后根據(jù)所設(shè)計(jì)的分層拓?fù)浣Y(jié)構(gòu)設(shè)計(jì)了與其相適應(yīng)的均衡控制策略;最后在Matlab/Simulink中分別搭建了不同拓?fù)湎碌姆抡婺P筒⑦M(jìn)行了均衡仿真驗(yàn)證。仿真結(jié)果表明,在8節(jié)電池與9節(jié)電池的不同情況下與傳統(tǒng)單層均衡拓?fù)湎啾?,文中所提出的采用相?duì)極差控制策略的分層均衡拓?fù)浣Y(jié)構(gòu)在均衡速度上分別提高了86.06%和84.21%,在能量轉(zhuǎn)移效率上分別提高了3.3%和4.22%,驗(yàn)證了所提均衡拓?fù)湟约熬獠呗缘挠行浴?/p>
關(guān)鍵詞: Buck?Boost電路; 鋰離子電池; 分層主動(dòng)均衡; 相對(duì)極差; 能量轉(zhuǎn)移效率; 均衡速度
中圖分類號(hào): TN609?34; TM91" " " " " " " " " " " 文獻(xiàn)標(biāo)識(shí)碼: A" " " " " " " " " "文章編號(hào): 1004?373X(2024)19?0169?06
Lithium battery pack hierarchical equalization technology
based on buck?boost circuit and relative range strategy
WEI Yewen1, 2, LI Yibo1, LIU Jiou1, QIU Tianyi1
(1. College of Electrical Engineering and New Energy, China Three Gorges University, Yichang 443002, China;
2. Hubei Provincial Collaborative Innovation Center for New Energy Microgrid, China Three Gorges University, Yichang 443002, China)
Abstract: In view of the slow equalization speed and low energy transfer efficiency in the equalization process of series?connected lithium battery packs, a hierarchical equalization topology based on the buck?boost circuit is proposed. In this topology, the state of charge (SOC) of each single cell is taken as the equalization variable, and the equalization control strategy of the relative range method is adopted to achieve the simultaneous equalization of the equalization units in each layer. The overall structure of the improved equalization topology and the method of stratified grouping, the equalization principle of the smallest equalization unit, and the calculation process of the relevant parameters in the equalization process are described, and then the equalization control strategy is designed according to the designed stratified topology. The simulation models of different topologies are constructed in Matlab/Simulink and the equalization simulation verification is carried out. In comparison with the traditional single?layer equalization topology, the equalization speed of the proposed hierarchical equalization topology with relative range control strategy is improved by 86.06% and 84.21%, respectively, and its energy transfer efficiency is improved by 3.3% and 4.22%, respectively, in the different cases of 8 cells and 9 cells. To sum up, the effectiveness of the proposed equalization topology and the equalization strategy are verified.
Keywords: buck?boost circuit; lithium?ion battery; hierarchical active equalization; relative range; energy transfer efficiency; equalization speed
0" 引" 言
隨著國家“雙碳”目標(biāo)的提出,構(gòu)建清潔高效的能源架構(gòu),保證能源的可持續(xù)發(fā)展已成為當(dāng)今世界各國共同努力的方向,在此背景下,風(fēng)能、太陽能等清潔可再生能源迅速發(fā)展,并成為目前主要的發(fā)電方式,但由于其存在間歇性,若直接并入電網(wǎng)會(huì)影響電網(wǎng)的安全可靠運(yùn)行。而大規(guī)模的電池儲(chǔ)能技術(shù)則能很好地解決上述新能源發(fā)電的間歇性和隨機(jī)性問題,更好地消納新能源,保證電網(wǎng)的安全可靠運(yùn)行[1]。
鋰電池憑借其高功率、高能量密度、使用壽命長(zhǎng)的特性成為目前廣泛應(yīng)用的儲(chǔ)能電池,但由于電池制作工藝的差異以及不同使用環(huán)境對(duì)于電池性能的影響,導(dǎo)致電池組中的單體電池在容量、阻抗等方面出現(xiàn)差異,而這種差異會(huì)隨著使用周期的變長(zhǎng)而逐漸增大,長(zhǎng)期如此不僅會(huì)導(dǎo)致單體電池組使用的性能下降,減少電池壽命,甚至?xí)斐呻姵亟M起火爆炸。為減少電池組的不一致性,提高電池組壽命和能量利用率,引入電池均衡技術(shù)。
目前,電池均衡技術(shù)按照是否消耗能量可將均衡電路分為兩種類型,一種是能量耗散型(被動(dòng)均衡),另一種是能量非耗散型(主動(dòng)均衡)。被動(dòng)均衡是指將電池能量高的電池通過旁路電阻耗散掉從而達(dá)到均衡的目的,這種均衡模式結(jié)構(gòu)簡(jiǎn)單、成本低,但存在著能量利用率低、熱管理困難等缺陷[2]。其中文獻(xiàn)[3]提出利用MOS管作為旁路開關(guān),每個(gè)單體電池都有一個(gè)串聯(lián)開關(guān)和并聯(lián)開關(guān)的均衡拓?fù)浣Y(jié)構(gòu),不僅會(huì)造成嚴(yán)重的開關(guān)損耗,大量的開關(guān)器件與單體電池之間的串聯(lián)還會(huì)導(dǎo)致電池組的供電可靠性下降。
近年來,主動(dòng)均衡以其出色的性能成為均衡技術(shù)研究的熱點(diǎn)。它通過借助電感、電容等儲(chǔ)能元件,將電池中高能量的電量轉(zhuǎn)移到低能量的電池中。根據(jù)所使用的儲(chǔ)能元件的不同,可以分為電容式、電感式、變壓器式和電力變換器式均衡電路。電容型均衡電路體積小,易于控制,但其精度較低,均衡時(shí)間長(zhǎng)且均衡目標(biāo)只能是電壓,存在一定的局限性[4]。電感型均衡電路不受電壓差的影響,同時(shí)相比于電容型均衡電路有較高的均衡精度,但其往往伴隨著較大的開關(guān)陣列,結(jié)構(gòu)復(fù)雜,控制過程復(fù)雜,可靠性低[5]?;谧儔浩鞯木馔?fù)渚哂芯庑矢?、控制?jiǎn)單,且易于隔離的優(yōu)點(diǎn),以及均衡過程中變壓器的型號(hào)對(duì)均衡效果的影響較大,不同的電池組需要單獨(dú)設(shè)計(jì)不同的變壓器,導(dǎo)致變壓器型均衡拓?fù)涞耐ㄓ眯暂^差的缺陷[6]。文獻(xiàn)[7]中提出的基于開關(guān)電感的Buck?Boost型電路,可實(shí)現(xiàn)電池之間的任意組合均衡,但均衡過程中會(huì)造成較大的開關(guān)損耗,使得能量利用率低。文獻(xiàn)[8]所提出的分層式均衡拓?fù)浣Y(jié)構(gòu)采用了電感以及反激式變壓器,可實(shí)現(xiàn)任意電池間的能量傳輸,但是成本高且控制復(fù)雜。
為解決現(xiàn)有拓?fù)渲写嬖诘木馑俣嚷?、均衡開關(guān)數(shù)量多、結(jié)構(gòu)復(fù)雜等問題,本文以Buck?Boost電路為基礎(chǔ)均衡單元,以電池的SOC為均衡目標(biāo),提出了基于Buck?Boost電路的分層均衡拓?fù)洌捎孟鄬?duì)極差法的均衡控制策略對(duì)電池組進(jìn)行分層同步均衡,最終在Matlab中搭建了8節(jié)電池和9節(jié)電池的仿真模型,驗(yàn)證了所提均衡拓?fù)渑c均衡策略的可行性。
1" 均衡電路拓?fù)浼霸矸治?/p>
1.1" 均衡電路拓?fù)?/p>
傳統(tǒng)的單層均衡拓?fù)渲荒苁鼓芰吭谙噜彽膯误w電池或電池組之間傳遞,能量傳輸通道單一,均衡時(shí)間長(zhǎng)、均衡效率低。改進(jìn)后的雙層均衡拓?fù)潆m解決了能量傳輸通道單一的問題,但當(dāng)電池組中單體電池?cái)?shù)量過多時(shí),往往需要更多的均衡開關(guān),導(dǎo)致均衡過程中的開關(guān)損耗增大,均衡速度降低[9]。為進(jìn)一步提高均衡速度,減少開關(guān)損耗,本文在傳統(tǒng)單層均衡拓?fù)涞幕A(chǔ)上提出了可有效提高均衡速度、減少開關(guān)損耗的分層均衡拓?fù)?。均衡拓?fù)湔w結(jié)構(gòu)如圖1所示。
第一層由兩節(jié)單體電池構(gòu)成一個(gè)電路單元,每個(gè)電路單元內(nèi)部可實(shí)現(xiàn)兩節(jié)單體電池之間的均衡,同時(shí)兩個(gè)相鄰的電路單元由第二層的均衡模塊進(jìn)行電路單元間的均衡,以此類推。當(dāng)單體電池?cái)?shù)[n]為偶數(shù)時(shí),每個(gè)電池單元由兩節(jié)單體電池構(gòu)成;當(dāng)[n]為奇數(shù)時(shí),最末端單體電池與其之前的兩節(jié)單體電池共同構(gòu)成一個(gè)三電池單元,其余各電池依舊按照上述兩兩成組構(gòu)成電路單元,以此類推,各層也按照此方法進(jìn)行均衡分組。
1.2" 均衡原理
本文采用Buck?Boost變換器作為分層均衡電路的均衡器,以8節(jié)單體電池為例,分層式主動(dòng)均衡拓?fù)涞恼w結(jié)構(gòu)如圖2所示。
相較于普通的串聯(lián)式均衡拓?fù)?,圖2中的分層式均衡拓?fù)淇扇龑踊芈罚瑫r(shí)開啟均衡實(shí)現(xiàn)任意單體電池間的能量均衡,大大提高了能量均衡效率。最小均衡單元結(jié)構(gòu)如圖3所示。
最小均衡單元由兩個(gè)MOS管、一個(gè)電感和兩節(jié)單體電池組成,其中電感作為能量?jī)?chǔ)存的中間載體,兩個(gè)MOS管作為開關(guān)由PWM信號(hào)驅(qū)動(dòng),兩者相互配合實(shí)現(xiàn)電池能量從一個(gè)單體電池向另一個(gè)單體電池的傳輸。為避免電感的磁飽和,選用電感的斷續(xù)導(dǎo)通模式,以保證在一個(gè)脈沖調(diào)制開關(guān)周期內(nèi)電感所吸收的能量能被完全釋放。
假設(shè)B1能量高于B2,整個(gè)均衡過程分為如下三個(gè)階段。
第一階段:PWM信號(hào)控制開關(guān)管[SH]導(dǎo)通,此時(shí)電池B1、續(xù)流二集管D1、[SH]和電感[L]組成一個(gè)回路,B1將多余的能量釋放到電感[L]中以磁能的形式儲(chǔ)存,此時(shí)回路中的電流記作[IL1],此過程如圖4所示。
第二階段:PWM信號(hào)控制關(guān)斷MOS管[SH],此時(shí)[SL]、電池B2、續(xù)流二極管[D2]以及電感[L]組成一個(gè)回路,將儲(chǔ)存在電感中的能量通過此回路釋放給電池B2,直至電感電流降為0,此階段結(jié)束。整個(gè)過程如圖5所示,回路中的電流記作[IL2]。
第三階段:MOS管[SH]和[SL]均截止,一個(gè)周期內(nèi)的均衡完成,等待下一個(gè)PWM信號(hào)的驅(qū)動(dòng)。如此循環(huán)使得兩個(gè)單體電池的能量逐漸趨于一致,最終小于均衡閾值,完成均衡。
1.3" 參數(shù)計(jì)算
忽略開關(guān)管的正向?qū)妷阂约皟蓚€(gè)階段中回路的總電阻,設(shè)[UB1]、[UB2]分別為電池B1和電池B2兩端的電壓,[D]為PWM控制信號(hào)的占空比,[toff]為第二階段的截止時(shí)間,[T]為整個(gè)PWM控制信號(hào)的周期。
不同階段電路中的電流如下:
[iL=iL1=UB1tL," " " "0lt;tlt;DTiL2=UB1DtL-UB2(t-DT)L," " " "DTlt;tlt;ts0," " " "tslt;tlt;T] (1)
整個(gè)開關(guān)周期內(nèi)電池B1向B2轉(zhuǎn)移的能量為:
[E=12UB1L(DT)2] (2)
在第一階段高電量單體電池放電過程中,回路中的電流[IL1]逐漸增大,等[t=DT]時(shí),回路中的電流達(dá)到最大,此時(shí)最大電流為[UB1tL],之后進(jìn)入低電量單體電池吸收電量的過程,此過程中MOS管[SH]關(guān)斷,[SL]導(dǎo)通,回路中的電流[IL2]逐漸減小。
為使電路工作在電感的斷續(xù)導(dǎo)通模式(DCM),需保證均衡結(jié)束時(shí)間[ts]小于PWM控制周期[T],即當(dāng)[t]lt;[T]時(shí),回路二中電流[IL2]需降至0,代入公式可得,占空比[D]需滿足:
[D≤UB2 (UB1+UB2)] (3)
2" 均衡控制策略
2.1" 均衡變量選取
電池的均衡控制策略中第一步就是要選取合適的均衡控制指標(biāo),以均衡控制指標(biāo)作為均衡效果的判斷依據(jù),目前常用的均衡控制指標(biāo)有電池的端電壓、電池的剩余容量以及電池的荷電狀態(tài)三種[10]。電池電壓的獲取較為簡(jiǎn)單,可直接采集,也是目前最為常用的均衡控制指標(biāo),但電池電壓的測(cè)量易受干擾,同時(shí)電池在充放電以及靜置等不同情況下同一單體電池的電壓值也會(huì)不同,因此以電壓為均衡目標(biāo)會(huì)導(dǎo)致均衡誤差較大,無法達(dá)到準(zhǔn)確的均衡。
相比于用電壓作為均衡目標(biāo),以電池的SOC作為均衡指標(biāo)時(shí),均衡精度更高,但由于電池的SOC無法直接獲取,往往需要借助其他變量進(jìn)行估算,目前常用的估算方法是開路電壓法。
選用額定容量為42 A·h、額定電壓為3.6 V的鋰電池進(jìn)行橫流間隙充放電,測(cè)得充放電過程中電池的開路電壓和荷電狀態(tài)特性曲線如圖6所示。
由鋰電池OCV?SOC特性曲線可知,電池SOC處于20%~90%階段時(shí),電池的開路電壓變化小,此時(shí)若選用電壓作為均衡指標(biāo),會(huì)導(dǎo)致電池單體之間電壓差距很小,但SOC相差很大,不僅會(huì)嚴(yán)重降低均衡速度,而且還容易出現(xiàn)誤均衡操作,因此本文選用電池的SOC作為均衡電路的均衡指標(biāo)。
2.2" 均衡控制策略
目前常用的均衡控制策略主要有最值法、平均值法、差值法以及模糊控制法等[11]。根據(jù)本文所采用的分層結(jié)構(gòu),提出了“相對(duì)極差法”,以單體電池或電池組的極差作為均衡閾值條件,采用分層同步均衡的方法實(shí)現(xiàn)單體電池以及電池組的一致性均衡。
相對(duì)極差是指電池組中單體電池SOC最大值或最小值與該電池組所有單體電池SOC平均值之差。記電池組中第[i]節(jié)單體電池的SOC值為[SOCi],電池組中第[k]個(gè)電池單元中所包含單體電池的平均值為[SOCk],電池單元中最大的SOC值和最小的SOC值分別為[SOCmax]和[SOCmin],[α+]和[α-]分別為正向相對(duì)極差和負(fù)向相對(duì)極差。則:
[SOCk=1ni=1nSOCi," " "k=1,2,…,n] (4)
[α+=SOCmax-SOCk] (5)
[α-=SOCk-SOCmin] (6)
均衡控制策略流程圖如圖7所示。
首先利用電池檢測(cè)模塊獲取電池的電壓、電流、溫度等參數(shù),利用開路電壓法估算各單體電池的SOC值[12],之后計(jì)算各電池單元以及各電池組的平均SOC以及相對(duì)差異度[α],判斷相對(duì)差異度是否達(dá)到均衡開啟條件,若達(dá)到則控制相應(yīng)的MOS管導(dǎo)通,對(duì)第一層電池單元內(nèi)部的兩節(jié)單體電池通過Buck?Boost均衡模塊進(jìn)行均衡,均衡完成后判斷第二層整體電池組的SOC平均值以及最大、最小SOC單元[13],按照上述方法判斷是否達(dá)到均衡閾值,通過上層的均衡模塊對(duì)該層電池單元進(jìn)行均衡,如此循環(huán)直至電池單元內(nèi)以及電池組之間差異小于所設(shè)定的均衡閾值,均衡結(jié)束。
3" 仿真結(jié)果分析
為驗(yàn)證所提分層拓?fù)浣Y(jié)構(gòu)及均衡控制策略的可行性,選用8節(jié)鋰電池,在Matlab/Simulink仿真實(shí)驗(yàn)平臺(tái)搭建了仿真模型,并分別在傳統(tǒng)單層均衡以及本文所提分層均衡之間進(jìn)行了對(duì)比仿真驗(yàn)證,具體的仿真參數(shù)如表1所示。
8節(jié)電池的SOC值分別設(shè)置為90%、88%、86%、84%、82%、80%、78%、72%。均衡過程及單體電池SOC變化值如圖8所示。
采用8節(jié)電池做仿真實(shí)驗(yàn),初始狀態(tài)下8節(jié)電池的SOC均值為82.5%,采用傳統(tǒng)單層均衡以及本文分層均衡后均衡所用時(shí)間以及最終達(dá)到均衡時(shí)SOC均值對(duì)比如表2所示。
由表2數(shù)據(jù)可以看出,本文所提分層均衡拓?fù)渌镁鈺r(shí)間比傳統(tǒng)單層均衡縮短了269.88 s,均衡速度提高了86.06%,同時(shí)根據(jù)均衡后的SOC值可算出本文分層均衡結(jié)構(gòu)較傳統(tǒng)單層結(jié)構(gòu)能量轉(zhuǎn)移效率提高了3.3%。
采用9節(jié)初始SOC值分別為90%、88%、86%、84%、78%、76%、75%、73%、72%的鋰電池,使用傳統(tǒng)單層以及本文所提分層均衡結(jié)構(gòu)在Matlab中搭建仿真模型并設(shè)置相關(guān)仿真參數(shù)后進(jìn)行仿真,各單體電池SOC的變化情況如圖9所示。
采用9節(jié)鋰電池做仿真實(shí)驗(yàn),初始狀態(tài)下9節(jié)電池的SOC均值為80.22%,采用傳統(tǒng)單層拓?fù)湟约氨疚姆謱油負(fù)溥_(dá)到均衡所用時(shí)間以及最終達(dá)到均衡后的SOC均值如表3所示。
由表3數(shù)據(jù)可以看出,本文所提分層均衡拓?fù)渌镁鈺r(shí)間比傳統(tǒng)單層均衡縮短了319.01 s,均衡速度提高了84.21%,同時(shí)根據(jù)均衡后的SOC值可算出本文分層均衡結(jié)構(gòu)較傳統(tǒng)單層結(jié)構(gòu)能量轉(zhuǎn)移效率提高了4.22%。
4" 結(jié)" 論
本文提出一種基于Buck?Boost電路的分層式主動(dòng)均衡結(jié)構(gòu),采用相對(duì)極差的均衡控制策略,通過對(duì)所提拓?fù)涞脑磉M(jìn)行分析,利用仿真驗(yàn)證得出了以下結(jié)論:
1) 相對(duì)于傳統(tǒng)的單層均衡模式,所提出的基于Buck?Boost電路的分層式主動(dòng)均衡可有效減少均衡時(shí)間,提高均衡效率,并且結(jié)構(gòu)相對(duì)簡(jiǎn)單,控制難度較低,可用于大規(guī)模儲(chǔ)能系統(tǒng)中,有效防止單體電池的過充過放,提高了電池組的安全可靠性。
2) 采用電池的SOC作為均衡目標(biāo),采用相對(duì)極差法,對(duì)各層實(shí)現(xiàn)同步均衡,相比于傳統(tǒng)的分層分步均衡,提高了均衡效率。
參考文獻(xiàn)
[1] 劉紅銳,李海瑞,韋向陽,等.一種模塊化的高性能蓄電池均衡器研究[J].電工技術(shù)學(xué)報(bào),2023,38(17):4574?4585.
[2] 劉政,張向文,黃斌.基于LTC3300?1的電動(dòng)汽車鋰電池組主動(dòng)均衡系統(tǒng)研究[J].電源技術(shù),2018,42(7):993?996.
[3] TASHAKOR N, FARJAH E, GHANBARI T. A bidirectional battery charger with modular integrated charge equalization circuit [J]. IEEE transactions on power electronics, 2019, 32(3): 2133?2145.
[4] 劉俊峰,高鵬舉,廖武兵,等.基于儲(chǔ)能電容復(fù)用的超級(jí)電容陣列自均衡拓?fù)浣Y(jié)構(gòu)的分析[J].電工技術(shù)學(xué)報(bào),2021,36(z1):372?380.
[5] WEI Y W, LIU G T, XIONG S N, et al. Research on power equalization using a low?loss DC?DC chopper for lithium?ion batteries in electric vehicle [C]// IOP Conference Series: Earth and Environmental Science. [S.l.: s.n.], 2017: 012038.
[6] LI Y, XU J, MEI X, et al. A unitized multiwinding transformer?based equalization method for series?connected battery strings [J]. IEEE transactions on power electronics, 2019, 34(12): 11981?11989.
[7] 李婷,朱孔軍,熊克.鋰離子電池自由成組分層均衡技術(shù)研究[J].電力電子技術(shù),2020,54(11):34?37.
[8] 巫春玲,程琰清,劉智軒,等.一種改進(jìn)的儲(chǔ)能鋰電池主動(dòng)均衡拓?fù)溲芯縖J].電子測(cè)量技術(shù),2021,44(3):32?37.
[9] 徐元中,郭純,吳鐵洲,等.基于Buck?Boost的鋰離子電池組分層均衡研究[J].電子器件,2023,46(5):1422?1428.
[10] 蔡敏怡,張娥,林靖,等.串聯(lián)鋰離子電池組均衡拓?fù)渚C述[J].中國電機(jī)工程學(xué)報(bào),2021,41(15):5294?5311.
[11] WANG L Y, WANG C S, YIN G, et al. Balanced control strategies for interconnected heterogeneous battery systems [J]. IEEE transactions on sustainable energy, 2016, 7(1): 189?199.
[12] GHAEMINEZHAD N, OUYANG Q, HU X S, et al. Active cell equalization topologies analysis for battery packs: A systematic review [J]. IEEE transactions on power electronics, 2021, 36(8): 9119?9135.
[13] 郭向偉,司陽,李文彪,等.串聯(lián)電池組雙目標(biāo)混合均衡控制研究[J].電源學(xué)報(bào),2021,19(5):185?190.
[14] SHANG Y L, XIA B, YANG J F, et al. A delta?structured switched?capacitor equalizer for series?connected battery strings [C]// 2019 IEEE Energy Conversion Congress and Exposition (ECCE). New York: IEEE, 2019: 4493?4496.
[15] UNO M, KUKITA A. String?to?battery voltage equalizer based on a half?bridge converter with multi?stacked current doublers for series?connected batteries [J]. IEEE transactions on power electronics, 2018, 34(2): 1286?1298.
作者簡(jiǎn)介:魏業(yè)文(1987—),男,湖北宜昌人,博士研究生,講師,碩士生導(dǎo)師,主要研究方向?yàn)樾履茉床⒕W(wǎng)與消納。
李藝博(2000—),男,河南平頂山人,碩士研究生,主要研究方向?yàn)殡姵亟M均衡控制拓?fù)?、電池組均衡控制策略。
劉霽歐(2000—),女,湖北武漢人,碩士研究生,主要研究方向?yàn)榫C合能源系統(tǒng)優(yōu)化調(diào)度。
邱天一(1999—),男,湖北黃岡人,碩士研究生,主要研究方向?yàn)殡娏﹄娮印?/p>
收稿日期:2024?03?15" " " " " "修回日期:2024?04?03
基金項(xiàng)目:國家自然科學(xué)基金面上項(xiàng)目資助(52377191)