文 萍，張 慶，李 傳，鄧文安

(中國(guó)石油大學(xué)(華東)化學(xué)工程學(xué)院，山東 青島 266555)

瀝青質(zhì)對(duì)稠油黏度的影響

文 萍，張 慶，李 傳，鄧文安

(中國(guó)石油大學(xué)(華東)化學(xué)工程學(xué)院，山東 青島 266555)

選用新疆地區(qū)的6種高黏稠油為主要研究對(duì)象，測(cè)定各稠油的性質(zhì)并配制不同瀝青質(zhì)含量的同源調(diào)合稠油、異源調(diào)合稠油，利用黏度法考察瀝青質(zhì)含量對(duì)稠油黏度的影響。結(jié)果表明：無論是同源調(diào)合稠油還是異源調(diào)合稠油，隨著瀝青質(zhì)添加量的增加，稠油黏度均呈增大的趨勢(shì)，且在瀝青質(zhì)添加量達(dá)到一定值時(shí)黏度出現(xiàn)突變；相同溫度、相同瀝青質(zhì)含量下，四氫呋喃處理過的稠油黏度比甲苯處理過的稠油黏度要大，且增加幅度較甲苯處理過的要大。

稠油 黏度 瀝青質(zhì) 同源 異源

隨著經(jīng)濟(jì)的發(fā)展，石油消耗量逐年增加，而輕質(zhì)原油可采量逐年減少，這使得稠油開發(fā)變得十分必要[1]。同時(shí)隨著稠油開采技術(shù)的不斷發(fā)展與提高，全球稠油的開采量也不斷提高。稠油是一類雜原子含量高、沸點(diǎn)高、相對(duì)分子質(zhì)量較大的復(fù)雜混合物[2-4]，具有密度高、黏度大、膠質(zhì)瀝青質(zhì)含量高等特點(diǎn)，嚴(yán)重制約著稠油的開采、集輸及加工利用。因此，降低稠油黏度、改善其流動(dòng)性是解決這些問題的關(guān)鍵。國(guó)內(nèi)外學(xué)者在稠油降黏技術(shù)的研究上做了大量工作。由于開采及加工的需要，人們對(duì)稠油的研究主要集中在流變性模型的構(gòu)建以及添加劑對(duì)稠油黏度的改善方面，對(duì)稠油本身的化學(xué)組成結(jié)構(gòu)與黏度關(guān)系的研究還沒有統(tǒng)一的定論?；诩僬f和模型的關(guān)于常規(guī)組分含量對(duì)黏度影響[5-6]的研究結(jié)果表明，膠質(zhì)和瀝青質(zhì)含量是影響稠油黏度的主要因素[7]。前期實(shí)驗(yàn)[8-11]考察了稠油瀝青質(zhì)的基本化學(xué)組成結(jié)構(gòu)與結(jié)構(gòu)參數(shù)，本課題選用新疆地區(qū)的6種高黏稠油為主要研究對(duì)象，測(cè)定各稠油的性質(zhì)，并配制不同瀝青質(zhì)含量的同源調(diào)合稠油、異源調(diào)合稠油，利用黏度法考察瀝青質(zhì)含量對(duì)稠油黏度的影響。

1 實(shí) 驗(yàn)

1.1 稠油的基本性質(zhì)

實(shí)驗(yàn)選用6種新疆地區(qū)高黏稠油，編號(hào)分別為Z32、Z18、208P、103P、風(fēng)城一站、LG15-9C，6種稠油的主要性質(zhì)見表1。由表1可見：稠油中輕質(zhì)餾分少，膠質(zhì)、瀝青質(zhì)含量較高，且隨著膠質(zhì)、瀝青質(zhì)含量的增加，稠油的密度、黏度也相應(yīng)增加；6種稠油黏度大小的順序與瀝青質(zhì)含量、殘?zhí)恳约俺碛蛿?shù)均相對(duì)分子質(zhì)量的順序基本一致，說明瀝青質(zhì)含量與黏度有一定的相關(guān)性。

1.2 四組分分離[12]

先將稠油用正庚烷沉淀，分離出庚烷瀝青質(zhì)(C7瀝青質(zhì))和C7可溶質(zhì)，然后將C7可溶質(zhì)在中性Al2O3色譜柱上進(jìn)行四組分分離，得到C7飽和分、C7芳香分、C7膠質(zhì)。

1.3 稠油正戊烷瀝青質(zhì)及可溶質(zhì)的分離[13]

將稠油用正戊烷沉淀，分離出戊烷瀝青質(zhì)(C5瀝青質(zhì))和C5可溶質(zhì)。

1.4 同源調(diào)合稠油的配制

同源調(diào)合稠油是指同一油樣的瀝青質(zhì)與可溶質(zhì)按不同比例調(diào)合后的稠油。

1.4.1 甲苯體系 在稠油的C5可溶質(zhì)中，分別加入質(zhì)量分?jǐn)?shù)0，2.0%，4.0%，6.0%，8.0%，10.0%，12.0%各自稠油的C7瀝青質(zhì)，用甲苯加熱回流溶解后，將溶劑蒸干，真空干燥，得到甲苯體系下的同源調(diào)合稠油。

表1 6種稠油的主要性質(zhì)

1.4.2 四氫呋喃體系 按1.4.1節(jié)方法，用極性相對(duì)較大的四氫呋喃溶劑代替甲苯溶劑，得到四氫呋喃體系下的同源調(diào)合稠油。

1.5 異源調(diào)合稠油的配制

實(shí)驗(yàn)在甲苯體系下完成，具體操作過程與同源調(diào)合稠油的配制過程基本一致，不同的是將其中的同源瀝青質(zhì)換成結(jié)構(gòu)性質(zhì)差距較大的異源瀝青質(zhì)。

1.5.1 Z32 以Z32可溶質(zhì)為母體，分別添加LG15-9C瀝青質(zhì)、Z32瀝青質(zhì)和208P瀝青質(zhì)，配制成瀝青質(zhì)質(zhì)量分?jǐn)?shù)分別為6%，8%，10%，12%的異源調(diào)合稠油，記為：A系列(Z32可溶質(zhì)+LG15-9C瀝青質(zhì))；B系列(Z32可溶質(zhì)+Z32瀝青質(zhì))；C系列(Z32可溶質(zhì)+208P瀝青質(zhì))。

1.5.2 208P 以208P可溶質(zhì)為母體，分別添加LG15-9C瀝青質(zhì)、Z32瀝青質(zhì)和208P瀝青質(zhì)，配制成瀝青質(zhì)質(zhì)量分?jǐn)?shù)分別為2%， 4%，6%，8%，10%的異源調(diào)合稠油，記為：D系列(208P可溶質(zhì)+LG15-9C瀝青質(zhì))；E系列(208P可溶質(zhì)+Z32瀝青質(zhì))；F系列(208P可溶質(zhì)+208P瀝青質(zhì))。

1.5.3 103P 以103P可溶質(zhì)為母體，添加LG15-9C瀝青質(zhì)、208P瀝青質(zhì)和103P瀝青質(zhì)，配制成瀝青質(zhì)質(zhì)量分?jǐn)?shù)分別為2%，4%，6%，8%，10%的調(diào)合稠油，記為：G系列(103P可溶質(zhì)+LG15-9C瀝青質(zhì))；H系列(103P可溶質(zhì)+208P瀝青質(zhì))；I系列(103P可溶質(zhì)+103P瀝青質(zhì))。1.5.4 Z18 以Z18可溶質(zhì)為母體，添加103P瀝青質(zhì)、208P瀝青質(zhì)和Z18瀝青質(zhì)，配制成瀝青質(zhì)質(zhì)量分?jǐn)?shù)分別為 4%，6%，8%，10%的調(diào)合稠油，記為L(zhǎng)系列(Z18可溶質(zhì)+103P瀝青質(zhì))；J系列(Z18可溶質(zhì)+208P瀝青質(zhì))；K系列(Z18可溶質(zhì)+ Z18瀝青質(zhì))。

1.6 調(diào)合稠油黏度的測(cè)定

利用旋轉(zhuǎn)黏度計(jì)測(cè)定不同溫度下的同源、異源調(diào)合稠油的黏度，繪制黏溫曲線。

2 結(jié)果與討論

2.1 瀝青質(zhì)對(duì)同源調(diào)合稠油黏度的影響

2.1.1 甲苯體系 在甲苯體系下，Z32、Z18、208P、103P、風(fēng)城一站同源調(diào)合稠油在50～80 ℃以及LG15-9C在較高溫度(90 ℃)下的黏度見圖1。從圖1可以看出：前5種稠油，在相同溫度下瀝青質(zhì)含量低時(shí)各同源調(diào)合稠油黏度較小；瀝青質(zhì)含量相同時(shí)，低溫下稠油黏度比高溫時(shí)要大；隨著瀝青質(zhì)含量的增加，各溫度下稠油黏度均呈增大的趨勢(shì)；瀝青質(zhì)含量相同時(shí)，各同源調(diào)合稠油高溫下黏度增加緩慢，說明瀝青質(zhì)含量是影響體系黏度的重要因素。產(chǎn)生以上現(xiàn)象的原因可能是膠質(zhì)中極性雜原子能夠起到瀝青質(zhì)膠溶劑的作用[14]，天然膠質(zhì)可作膠溶劑，對(duì)瀝青質(zhì)起到“稀釋”作用，并且瀝青質(zhì)被膠溶化后也產(chǎn)生一定的解締作用從而表現(xiàn)出黏度較小。低溫下布朗運(yùn)動(dòng)不劇烈，網(wǎng)狀結(jié)構(gòu)緊密，而高溫則對(duì)瀝青質(zhì)締合體具有一定的解締作用，使黏度變小。稠油是一個(gè)膠體體系，膠體體系中膠質(zhì)是抑制瀝青質(zhì)聚沉的關(guān)鍵，當(dāng)瀝青質(zhì)含量相對(duì)較高而膠質(zhì)不足以包裹瀝青質(zhì)時(shí)，對(duì)瀝青質(zhì)的穩(wěn)定作用不足，體系的穩(wěn)定性遭到破壞，導(dǎo)致瀝青質(zhì)的大量聚沉，宏觀表現(xiàn)為稠油黏度的驟然增加；另外，稠油瀝青質(zhì)中所含的硫、氮、氧雜原子是氫鍵、電荷轉(zhuǎn)移、偶極等非共價(jià)鍵作用的主要來源，在這幾種相互作用下，瀝青質(zhì)分子產(chǎn)生很強(qiáng)的內(nèi)聚力，容易發(fā)生締合聚集，使瀝青質(zhì)分子聚集成層狀堆積狀態(tài)，當(dāng)稠油分子間發(fā)生相對(duì)位移時(shí)可以產(chǎn)生較大的內(nèi)摩擦力，使得稠油產(chǎn)生高黏度。因?yàn)轲ざ确从车氖橇黧w內(nèi)部分子之間的摩擦力，瀝青質(zhì)含量的增加增大了流體內(nèi)部分子之間的摩擦力，從而使黏度增大。從圖1還可以看出，在較高溫度(90 ℃)時(shí)，LG15-9C瀝青質(zhì)質(zhì)量分?jǐn)?shù)為30%處出現(xiàn)非常明顯的突躍點(diǎn)，主要原因是當(dāng)添加瀝青質(zhì)的量超過其原油自身所含瀝青質(zhì)的量時(shí)，膠質(zhì)含量相對(duì)減少，瀝青質(zhì)含量偏高，它與膠質(zhì)的比例發(fā)生了變化，膠質(zhì)可以抑制瀝青質(zhì)納米聚集體的形成，促進(jìn)瀝青質(zhì)分散穩(wěn)定，但此時(shí)膠質(zhì)比例減小破壞了體系的穩(wěn)定性，進(jìn)而使黏度驟然增加。

圖1 甲苯體系下瀝青質(zhì)添加量對(duì)同源調(diào)合稠油黏度的影響■—50 ℃； ●—60 ℃； ▲—70 ℃； —80 ℃； ◆—90 ℃

2.1.2 四氫呋喃體系 在四氫呋喃體系下，Z32、Z18、208P、103P、風(fēng)城一站同源調(diào)合稠油在50～70 ℃以及LG15-9C在較高溫度(90 ℃)下的黏度見圖2。從圖2可以看出：前5種稠油，隨瀝青質(zhì)含量的增加，同源調(diào)合稠油的黏度也在不斷增大，且對(duì)溫度也有較強(qiáng)的敏感性，當(dāng)瀝青質(zhì)質(zhì)量分?jǐn)?shù)小于8%時(shí)，調(diào)合稠油黏度增加緩慢，瀝青質(zhì)質(zhì)量分?jǐn)?shù)大于8%時(shí)，稠油黏度急劇增大，出現(xiàn)突躍點(diǎn)。從圖2還可以看出，LG15-9C調(diào)合稠油的黏度在瀝青質(zhì)質(zhì)量分?jǐn)?shù)為30%處出現(xiàn)突躍。溫度升高時(shí)稠油黏度降低，且低溫下隨瀝青質(zhì)含量的增加稠油黏度增長(zhǎng)的幅度較大，說明低溫、瀝青質(zhì)含量高時(shí)黏度對(duì)溫度最敏感。

從圖1和圖2對(duì)比可以看出，在相同溫度、相同瀝青質(zhì)含量的條件下，四氫呋喃處理過的稠油黏度比甲苯處理過的稠油黏度要大。出現(xiàn)上述現(xiàn)象的原因可能是因?yàn)闉r青質(zhì)是稠油四組分中相對(duì)分子質(zhì)量最大、極性最強(qiáng)的組分，其富集的雜原子導(dǎo)致瀝青質(zhì)內(nèi)部存在氫鍵、電荷轉(zhuǎn)移、偶極等非共價(jià)鍵作用。甲苯是極性溶劑，改變?yōu)r青質(zhì)的極性導(dǎo)致黏度增大，四氫呋喃溶劑的極性比甲苯的大，且存在O原子容易與瀝青質(zhì)形成氫鍵，從而增加瀝青質(zhì)間的相互作用。四氫呋喃溶劑溶解瀝青質(zhì)后又將其蒸干除去，這一過程較甲苯處理時(shí)可能改變了瀝青質(zhì)的分布狀態(tài)，在完全除去溶劑后瀝青質(zhì)的非共價(jià)鍵作用特別是氫鍵作用的改變導(dǎo)致瀝青質(zhì)容易富集、聚沉，從而使稠油黏度增加，且增加幅度較甲苯(不含雜原子的極性溶劑)處理過的大。這主要與瀝青質(zhì)的締合性相關(guān)，而瀝青質(zhì)的締合性又與膠質(zhì)、瀝青質(zhì)在數(shù)量、組成結(jié)構(gòu)性質(zhì)的相容匹配性有關(guān)，同時(shí)也與其分子中所含雜原子的數(shù)目及形態(tài)有關(guān)[15]。

圖2 四氫呋喃體系下瀝青質(zhì)添加量對(duì)同源調(diào)合稠油黏度的影響■—50 ℃； ●—60 ℃； ▲—70 ℃； ◆—90 ℃

2.2 瀝青質(zhì)對(duì)異源調(diào)合稠油黏度的影響

2.2.1 Z32異源調(diào)合稠油 在50 ℃下，瀝青質(zhì)添加量對(duì)Z32異源調(diào)合稠油黏度的影響見圖3。從圖3可以看出，隨著瀝青質(zhì)含量的增大，A，B，C系列異源調(diào)合稠油的黏度均增大，且在瀝青質(zhì)質(zhì)量分?jǐn)?shù)為8%處出現(xiàn)突變，但增加的幅度各有差別，溫度及瀝青質(zhì)含量相同時(shí)，調(diào)合稠油黏度由大到小的順序?yàn)锳系列>B系列>C系列。

圖3 瀝青質(zhì)添加量對(duì)Z32異源調(diào)合稠油黏度的影響▲—A系列； ■—B系列； ●—C系列

2.2.2 208P異源調(diào)合稠油 在50 ℃下，瀝青質(zhì)添加量對(duì)208P異源調(diào)合稠油黏度的影響見圖4。從圖4可以看出，隨著瀝青質(zhì)含量的增大，D，E，F(xiàn)系列異源調(diào)合稠油的黏度均增大，且在瀝青質(zhì)質(zhì)量分?jǐn)?shù)為6%處出現(xiàn)突變，溫度及瀝青質(zhì)含量相同時(shí)，調(diào)合稠油黏度由大到小的順序?yàn)镈系列>E系列>F系列。

圖4 瀝青質(zhì)添加量對(duì)208P異源調(diào)合稠油黏度的影響▲—D系列； ●—E系列； ■—F系列

2.2.3 103P異源調(diào)合稠油 在50 ℃下，瀝青質(zhì)添加量對(duì)103P異源調(diào)合稠油黏度的影響見圖5。從圖5可以看出，隨著瀝青質(zhì)含量的增大，調(diào)合稠油的黏度均在瀝青質(zhì)質(zhì)量分?jǐn)?shù)為8%處發(fā)生突躍。可見當(dāng)瀝青質(zhì)含量在一定范圍(原油自身瀝青質(zhì)含量)內(nèi)，一定量的膠質(zhì)可以較好地膠溶化定量的瀝青質(zhì)，穩(wěn)定性較好，黏度隨不同來源瀝青質(zhì)含量的增加均緩慢增大，瀝青質(zhì)的量對(duì)黏度的影響占據(jù)主導(dǎo)地位；另一方面，當(dāng)瀝青質(zhì)的量超過一定值，膠質(zhì)不足以包裹瀝青質(zhì)后黏度會(huì)發(fā)生突變，但從異源調(diào)合稠油黏度突變后不同的增長(zhǎng)趨勢(shì)可以看出，此時(shí)瀝青質(zhì)的結(jié)構(gòu)與締合性對(duì)稠油黏度的影響占主導(dǎo)地位?？傊?，稠油的宏觀黏度是以上幾個(gè)方面綜合作用的結(jié)果。

圖5 瀝青質(zhì)添加量對(duì)103P異源調(diào)合稠油黏度的影響■—G系列； ●—L系列； ▲—H系列

2.2.4 Z18異源調(diào)合稠油 在50 ℃下，測(cè)定以Z18可溶質(zhì)為母體的異源調(diào)合稠油的黏度，結(jié)果見圖6。從圖6可以看出，在瀝青質(zhì)質(zhì)量分?jǐn)?shù)小于8%時(shí)，隨異源瀝青質(zhì)加入量的增多，黏度增長(zhǎng)趨勢(shì)沒有太大差別，但當(dāng)瀝青質(zhì)質(zhì)量分?jǐn)?shù)大于8%時(shí)，黏度發(fā)生突變，并且其增長(zhǎng)幅度也發(fā)生變化，Z18調(diào)合稠油黏度與103P稠油瀝青質(zhì)所配異源調(diào)合稠油黏度差距較小，而與208P稠油瀝青質(zhì)所配異源調(diào)合稠油黏度差距較大。

圖6 瀝青質(zhì)添加量對(duì)Z18異源調(diào)合稠油黏度的影響■—L系列； ●—K系列； ▲—J系列

3 結(jié) 論

(1) 隨著瀝青質(zhì)含量的增加，不同溫度下同源調(diào)合稠油黏度均呈增大的趨勢(shì)；溫度高時(shí)，黏度增幅減小。瀝青質(zhì)含量較低時(shí)，隨著瀝青質(zhì)含量增加黏度增加緩慢，對(duì)溫度變化也不敏感，近似成線性關(guān)系；當(dāng)瀝青質(zhì)含量較高時(shí)，曲線開始偏離線性關(guān)系，增加幅度變大，黏度對(duì)溫度變化也較為敏感；瀝青質(zhì)含量達(dá)到一定值后，黏度會(huì)出現(xiàn)突變點(diǎn)。

(2) 相同溫度、相同瀝青質(zhì)含量下，四氫呋喃處理過的稠油黏度比甲苯處理過的稠油黏度要大。且增加幅度較甲苯處理過的要大?？赡苁怯盟臍溥秽渲普{(diào)合稠油時(shí)改變了瀝青質(zhì)的分布狀態(tài)，在完全除去溶劑后瀝青質(zhì)的非共價(jià)鍵作用特別是氫鍵作用的增強(qiáng)導(dǎo)致瀝青質(zhì)容易富集、聚沉，從而使稠油黏度增加，且增加幅度較甲苯(不含雜原子的極性溶劑)處理過的大。

(3) 隨著瀝青質(zhì)添加量的增大，各異源調(diào)合稠油黏度均增大，且在瀝青質(zhì)含量達(dá)到一定值時(shí)出現(xiàn)突變，突變點(diǎn)與同源調(diào)合稠油黏度突變點(diǎn)基本一致，但異源調(diào)合稠油黏度增加的幅度各有差別。說明當(dāng)瀝青質(zhì)含量在一定范圍(原油自身瀝青質(zhì)含量)內(nèi)，瀝青質(zhì)的量對(duì)黏度的影響占據(jù)主導(dǎo)地位。

[1] Hinkle A，Batzle M.Heavy oils：A worldwide overview[J].The Leading Edge，2006，25(6)：742-749

[2] 黃娟，李本高，秦冰，等.塔河稠油催化降黏機(jī)理研究[J].石油煉制與化工，2014，45(5)：16-19

[3] 李庶峰，文萍，鄧文安，等.輪古稠油與煤焦油混合原料懸浮床加氫裂化研究[J].石油煉制與化工，2007，38(10)：25-28

[4] 沈明歡，王振宇，秦冰，等.化學(xué)降黏劑對(duì)稠油采出液破乳的影響[J].石油煉制與化工，2013，44(1)：90-94

[5] 汪雙清，沈斌，林壬子.稠油黏度與化學(xué)組成的關(guān)系[J].石油學(xué)報(bào)(石油加工)，2010，26(5)：795-798

[6] 朱戰(zhàn)軍，林壬子，汪雙清.稠油主要族組分對(duì)其黏度的影響[J].新疆石油地質(zhì)，2004，25(10)：512-513

[7] 周風(fēng)山，吳瑾光.稠油化學(xué)降黏技術(shù)研究進(jìn)展[J].油田化學(xué)，2001，18(3)：268-27

[8] 張慶，鄧文安，李傳，等.稠油瀝青質(zhì)的基本化學(xué)組成結(jié)構(gòu)與締合性研究[J].石油煉制與化工，2014，45(6)：20-24

[9] 文萍，崔敏，李傳，等.稠油及其組分結(jié)構(gòu)參數(shù)和紅外光譜的考察[J].石油化工高等學(xué)校學(xué)報(bào)，2013，26(4)：11-15

[10]文萍，王曉丹，鞏磊，等.塔里木稠油及其組分結(jié)構(gòu)參數(shù)的研究[J].石油化工高等學(xué)校學(xué)報(bào)，2013，26(5)：6-9

[11]文萍，崔敏，李傳，等.塔里木稠油及其組分中元素的分布和紅外光譜分析[J].石油化工高等學(xué)校學(xué)報(bào)，2013，26(6)：16-21

[13]文萍，石斌，于道永，等.大慶減壓渣油熱轉(zhuǎn)化前后鎳分布的研究[J].燃料化學(xué)學(xué)報(bào)，2002，30(5)：477-480

[14]張春明，趙紅靜，肖乾華，等.稠油黏度預(yù)測(cè)模型[J].長(zhǎng)江大學(xué)學(xué)報(bào)(自科版)，2005，2(7)：216-217

[15]Mojelsky T W，Ignasiak T M，F(xiàn)rakman Z，et al.Structural features of Alberta oil sand bitumen and heavy oil asphaltenes[J].Energy & Fuels，1992，6(1)：83-96

EFFECT OF ASPHALTENE ON VISCOSITY OF VISCOUS OILS

Wen Ping， Zhang Qing， Li Chuan， Deng Wen’an

(CollegeofChemicalEngineering，UniversityofPetroleum(EastChina)，Qingdao，Shandong266555)

Six kinds of Xinjiang heavy oils were used as the main research objects. The basic compositions and properties of the heavy oils were determined. The homologous and heterologous blended heavy oils were prepared with different contents of asphaltenes to study the influence of content of asphaltene on viscosity of blended heavy oils by viscosity method. The results show that the viscosity of both of homologous and heterologous heavy oils increases as asphaltene content increases； the viscosity increases sharply when the asphaltene content reaches a certain value. At the same temperature and asphaltene content， the viscosity of heavy oil treated by tetrahydrofuran is higher than that by methylbenzene. The extent of viscosity increase of the oils treated by tetrahydrofuran is larger than that by methylbenzene.

viscous oil； viscosity； asphaltene； homology； heterology

2015-12-08； 修改稿收到日期： 2016-03-10。

文萍，碩士，高級(jí)工程師，主要從事應(yīng)用化學(xué)與油田化學(xué)的實(shí)驗(yàn)教學(xué)及重油化學(xué)組成與結(jié)構(gòu)、劣質(zhì)重油輕質(zhì)化等方面的研究工作，發(fā)表論文30余篇。

文萍，E-mail：wenping@upc.edu.cn。

國(guó)家自然科學(xué)基金青年基金項(xiàng)目(21106186)。