太湖區(qū)小微水體氮素內(nèi)源釋放-沉降動態(tài)規(guī)律

摘要：為探究小微水體氮素內(nèi)源釋放—沉降的動態(tài)規(guī)律，本研究于2023年1月至11月期間，在太湖平原水網(wǎng)區(qū)選取了5個具有代表性的溝、塘、浜等小微水體進行連續(xù)監(jiān)測。每月定期對上覆水體、沉積物以及沉降顆粒物樣品進行采集，并同步測定氮濃度及理化因子?；诳紫端?dāng)U散模型，計算了沉積物—水界面NH+4、NO3擴散通量，并同步計算了NH+4、NO3沉降通量和NH+4、NO3內(nèi)源釋放通量。研究結(jié)果表明，太湖平原水網(wǎng)區(qū)小微水體的氮素內(nèi)源釋放以NH+4為主，NH+4內(nèi)源釋放總量在224-9 631 kg·km-2·a-1之間，擴散通量介于50-312 kg·km-2·a-1，對水體氮濃度的貢獻率為2.2%-15.1%。相較于河流、湖泊等大型水體，小微水體的內(nèi)源貢獻占比相對較低。進一步分析表明，上覆水與沉積物之間的濃度差以及流速等因素，是影響小微水體內(nèi)源釋放的關(guān)鍵。小微水體氮素內(nèi)源污染較為嚴重，為有效降低氮素內(nèi)源污染，應(yīng)強化沉積物管理，嚴格把控外源污染物的輸入，積極提升水體的自然凈化能力，合理調(diào)控水體流速，并構(gòu)建完善的監(jiān)測與預(yù)警體系。

關(guān)鍵詞：小微水體；內(nèi)源污染；沉積物—水界面；內(nèi)源氮釋放

中圖分類號：X52 文獻標(biāo)志碼：A 文章編號：1672-2043（2024）08-1859-10 doi：10.11654/jaes.2024-0024

隨著我國社會經(jīng)濟的高速發(fā)展以及人類活動的加劇，陸地系統(tǒng)中的氮、磷等營養(yǎng)物質(zhì)逐漸飽和，大量氮、磷進入水體，其中很大比例會通過沉降、擴散等形式最終蓄積在沉積物中，使得沉積物的氮素內(nèi)源釋放潛力增大。在外源性污染得到有效控制后，沉積物內(nèi)源釋放已經(jīng)成為水體富營養(yǎng)化的關(guān)鍵因素。有研究表明，富營養(yǎng)化嚴重的水體即使停止了外源性輸入，水體富營養(yǎng)化依舊會持續(xù)數(shù)十年，這也是部分水體富營養(yǎng)化治理成效不顯著的主要原因。

大量研究表明，許多湖泊、河流的沉積物內(nèi)源性氮污染較強。逄勇等在研究風(fēng)浪對營養(yǎng)鹽內(nèi)源釋放的影響中發(fā)現(xiàn)，太湖年均TN內(nèi)源釋放量高達7 773.0 t。張?zhí)鹉鹊仍趯Π籽蟮沓练e物—水界面的氮、磷營養(yǎng)鹽交換通量研究中發(fā)現(xiàn)，白洋淀沉積物—水界面的NH+4擴散通量介于1.71-7.43 mg·m-2·d-1之間，總?cè)芙庑缘獢U散通量的均值高達9.11 mg·m-2·d-1。張淑珍等在對釜陽河沉積物的氮素分布特征及交換通量的研究中發(fā)現(xiàn)，釜陽河部分河段的NH+4擴散通量介于48.9-1471.0 μmol·m-2·d-1。同時，也有許多學(xué)者對氮素的內(nèi)源釋放規(guī)律進行了深入研究，沉積物內(nèi)源釋放能力的強弱主要取決于沉積物及其上覆水體的物理化學(xué)性質(zhì)。首先，沉積物中的氮含量決定了其釋放能力，當(dāng)沉積物中的氮含量越高，其釋放到上覆水中的潛力也就越大。崔會芳等在橫江水庫的研究中發(fā)現(xiàn)，沉積物氮釋放能力和沉積物氮含量呈顯著正相關(guān)關(guān)系。其次，水動力因素對沉積物內(nèi)源釋放也有著顯著的影響。胡開明等、童亞莉等在研究水體擾動對沉積物氮素內(nèi)源釋放的研究中發(fā)現(xiàn)，隨著流速的增大，沉積物氮素內(nèi)源釋放強度逐漸增大。此外，水溫、pH、DO等因素對沉積物氮素內(nèi)源釋放也有一定的影響。

上述研究主要聚焦于氮素的內(nèi)源釋放潛力，多局限于實驗室單一因素，而在野外實際環(huán)境中，多因素綜合作用下的內(nèi)源釋放規(guī)律及其對水質(zhì)的具體貢獻，尚缺乏深入的研究。另外，內(nèi)源釋放通量主要基于Fick定律的孔隙水濃度擴散模型法。在計算過程中，大多研究只考慮了界面擴散通量總量，并未明確區(qū)分釋放和沉降這兩個過程。此外，現(xiàn)有的研究多集中于湖泊、河流等大型水體，對于農(nóng)村小微水體沉積物的內(nèi)源釋放研究相對匱乏。小微水體是接納污染物的第一道水域，由于受高強度人類活動的影響，大量生活污水直接排放，加之小微水體較淺、水力條件差等因素，氮素更容易在沉積物中積聚，小微水體內(nèi)源釋放潛力巨大。基于此，本研究以太湖地區(qū)典型小微水體溝、塘、浜為研究對象，基于Fick擴散定律模型，估算沉積物—水界面氮素釋放—沉降過程，并進一步分析內(nèi)源釋放、沉降與各環(huán)境因子之間的影響關(guān)系。

1 材料與方法

1.1 研究區(qū)概況

本研究選定江蘇省常熟市辛莊鎮(zhèn)為研究區(qū)域，辛莊鎮(zhèn)地處東經(jīng)120°33＇-121°03＇、北緯31°31＇-31°50＇，屬太湖流域下游，是平原水網(wǎng)地帶的典型代表。該區(qū)域水網(wǎng)密集，河道縱橫交錯，湖泊星羅棋布，水網(wǎng)面積占比達15.3%。區(qū)域內(nèi)溝渠、池塘、河浜等小微水體分布廣泛，農(nóng)業(yè)灌溉和養(yǎng)殖活動頻繁。此外，區(qū)域內(nèi)生活和農(nóng)業(yè)區(qū)塊分布明顯，河浜貫穿其中，水體面源污染嚴重。為深入研究小微水體氮素內(nèi)源釋放-沉降的動態(tài)規(guī)律，綜合考慮研究區(qū)域內(nèi)不同類型的小微水體的分布、占比及上覆水體、沉積物氮素濃度情況，本研究在辛莊鎮(zhèn)精心選取了5個采樣點：1個魚塘（Q1）、2條生活區(qū)河浜（Q2、Q3）、1條農(nóng)業(yè)溝渠（Q4）及1條農(nóng)業(yè)區(qū)河浜（Q5），采樣點分布見圖1。在這些采樣點進行了上覆水、沉積物和沉降顆粒物的采集工作。

1.2 樣品采集與測定

于2023年1-11月每月1次在上述采樣點進行采樣工作，在每個采樣點，使用200 mL無菌水樣袋收集3袋水樣，使用柱狀樣采樣器于距岸邊1-2 m處采集3柱沉積物，每個沉積物柱狀樣的采集深度為10cm。沉降顆粒物通過改進后的沉降收集裝置進行收集，裝置由收集桶和濾網(wǎng)兩部分組成，濾網(wǎng)防止魚蝦、碎石等進入，收集桶則收集沉降顆粒物。沉降收集裝置收集一個月之后進行打撈，將所收集到的沉降顆粒物使用自封袋保存，并重新投放繼續(xù)收集。在上述采樣的同時，現(xiàn)場使用便攜式多參數(shù)水質(zhì)監(jiān)測儀（YSI EXO Multiparameter probe，美國）測定水樣pH、DO、鹽度（S）和水溫（WT），使用濁度懸浮物檢測儀測定水樣濁度、懸浮物（SS），使用便攜式流速儀測定水體流速。樣品采集完成后，在無擾動情況下，將所有柱狀樣品運至實驗室內(nèi)，確保沉積物—水界面處于接近原位環(huán)境的狀態(tài)，將采集到的柱狀樣在與原位溫度相近的溫度條件下穩(wěn)定24 h，以確保柱狀樣沉積物和上覆水達到動態(tài)平衡。平衡后在柱狀樣水-沉積物界面下方1、3、5、7、9 cm處插入Rhizon孔隙水采樣管，使用30 mL注射器進行孔隙水的采集。

上覆水樣帶回實驗室經(jīng)過0.45 μm過濾后，使用流動分析儀（荷蘭SAKLAR公司San++System）測定濾液的氨態(tài)氮（NH+4）、硝態(tài)氮（NO-3），總?cè)芙庑缘═DN）。因為Rhizon孔徑為0.25 μm，所以孔隙水無需過濾，直接使用流動分析儀測定其氨態(tài)氮（NH+4）、硝態(tài)氮（NO-3）。沉積物樣和沉降顆粒物樣取部分使用2 mol·L-1 KCl溶液按1：5的比例浸提，振蕩2h后過濾，使用流動分析儀測定濾液的氨態(tài)氮（NH+4）、硝態(tài)氮（NO-3）。

1.3 沉積物氮素交換速率及其對水質(zhì)的貢獻計算方法

沉積物—水界面營養(yǎng)鹽進行交換的過程主要包括對流通量和擴散通量，一般在流速較低的情況下，對流通量遠小于擴散通量，因此營養(yǎng)鹽的遷移過程主要由界面上下的濃度梯度所引起的擴散轉(zhuǎn)移過程所主導(dǎo)叫。營養(yǎng)鹽的分布主要受一維垂向傳輸和反應(yīng)過程控制，因此，營養(yǎng)鹽在沉積物-水界面的擴散服從Fick第一擴散定律，由濃度梯度所引起的內(nèi)源釋放是沉積物與上覆水之間進行物質(zhì)交換的主要形式。其通量可以通過Fick第一定律進行計算。根據(jù)Fick第一定律，可以用下式估算溶質(zhì)的擴散通量（F）：

F=φ×Ds×?c/?x（1）

式中：F為通過沉積物—水界面的擴散通量，mg·m-2·d-1；φ為沉積物孔隙度，可通過表層沉積物含水率推算：?c/?x為沉積物—水界面間隙水中溶質(zhì)的濃度梯度，通過表層沉積物孔隙水濃度c與深度x指數(shù)擬合在沉積物-水界交界面處求導(dǎo)可得，mg·L-1·cm-1；Ds為實際溶液中溶質(zhì)的分子擴散系數(shù)，m2·s-1。通常當(dāng)φ≤0.7時，Ds=φ×D0；當(dāng)φ＞0.7時，Ds=φ2×D0；式中D0為理想溶液中溶質(zhì)的分子擴散系數(shù)，m2·s-1。φ的計算公式如下：

φ＝（Ww - Wd）×100%/[（Ww-wd）+ Wd／ρ]（2）

式中：Ww為沉積物鮮質(zhì)量，g；Wd為沉積物干質(zhì)量，g;ρ為表層沉積物平均密度與水密度的比值，一般取2.5。根據(jù)φ和D0，結(jié)合間隙水剖面分析，由公式（1）可以計算出各沉積物-水界面氮鹽的擴散通量。

沉降顆粒物通過投放收集裝置進行收集，沉降通量（則可按下式進行估算：

R=-c×M/（A×T）（3）

式中：及為沉降通量，mg·m-2·d-1；c為收集的顆粒物氮素濃度，mg·kg-1；M為收集的顆粒物質(zhì)量；A為收集裝置底面積，m2；T為收集時間，d。

沉積物-水界面的擴散通量（F）是沉積物內(nèi)源釋放和上覆水沉降的綜合作用結(jié)果，F(xiàn)＞0時，內(nèi)源釋放強度大于沉降強度，此時沉積物表現(xiàn)為“源”；F＜0時，內(nèi)源釋放強度小于沉降強度，沉積物表現(xiàn)為“匯”。由此，基于沉積物—水界面的擴散通量（F），結(jié)合收集沉降顆粒物估算的沉降通量（R），可對內(nèi)源釋放通量（J）進行估算：

J=F-R（4）

式中：J為1d時間下，1 m2沉積物向上覆水釋放的營養(yǎng)液質(zhì)量。引入水深（H），J為1d時間下，沉積物向底面積為1 m2，高為H的上覆水柱釋放的營養(yǎng)液質(zhì)量。內(nèi)源釋放對水質(zhì)的負荷可用下式進行估算：

P=J×T/H（5）

式中：P為內(nèi)源釋放對水質(zhì)的負荷，mg·L-1；T為釋放時間，d；H為采樣點水深，m。

沉積物—水界面擴散過程（釋放凈值）對水質(zhì)的負荷將公式（5）中的J替換為擴散通量（F）計算。

2 結(jié)果和分析

2.1 沉積物-水界面理化性質(zhì)與NH+4、NO-3濃度剖面分布

研究區(qū)水體理化性質(zhì)如表1所示，各采樣點D0含量為2.59-14.41 mg·L-1；pH值在7.56-8.38之間，水體呈弱堿性；流速范圍為0.010-0.027 m·s-1；在采樣期間，研究區(qū)域上覆水體NH+4濃度變化范圍為0.03-2.53 mg·L-1，平均值為0.31 mg·L-1；NO-3濃度變化范圍為0.36-1.48 mg·L-1，平均值為0.89 mg·L-1。

各采樣點孔隙水的NH+4、NO3-濃度剖面分布如圖2所示，孔隙水NH+4濃度變化范圍為0.219-16.058 mg·L-1，平均值為3.096 mg·L-1，孔隙水與上覆水NH+4濃度比值范圍為24.14-52.33；孔隙水NO-3濃度變化范圍為0.005-1.967 mg·L-1，平均值為0.034 mg·L-1，孔隙水與上覆水NO-3濃度比值范圍為0.19-0.61。各采樣點中孔隙水NH+4的濃度都遠高于NO-3，且孔隙水NH+4、NO-3濃度都隨深度增加而增大，表明各點都具有很強的內(nèi)源釋放潛力，且NH+4內(nèi)源釋放潛力遠高于NO-3。各采樣點孔隙水NH+4的濃度Q1最大，其次是Q2、Q3，而Q4、Q5最小，孔隙水作為氮素從沉積物釋放到上覆水的過渡區(qū)域，其氮素濃度分布與沉積物中氮素濃度分布一致。

2.2 沉積物-水界面NH+4、NO-3擴散通量

營養(yǎng)鹽在沉積物—水界面的擴散服從Fick第一擴散定律，一般情況下，孔隙水中營養(yǎng)鹽濃度隨深度變化服從指數(shù)分布規(guī)律。因此，將沉積物孔隙水（0-10 cm）處的NH+4、NO-3濃度對深度進行指數(shù)擬合，進一步按照擬合的表達式對深度。進行求導(dǎo)，求出在沉積物—水界面處的?c/?x，并利用Fick定律計算出NH+1、NO-3在沉積物—水界面的擴散通量。

從孔隙水與上覆水體NH+4、NO-3濃度對深度的指數(shù)擬合結(jié)果來看，各采樣點的相關(guān)系數(shù)都較高，R2介于0.501-0.954之間（表2），說明對于該區(qū)域各點沉積物—水界面NH+4、NO-3擴散過程可以通過一級反應(yīng)動力學(xué)來描述。根據(jù)各采樣點的實測數(shù)據(jù)，通過Fick擴散定律計算得到的沉積物—水界面NH+4、NO-3擴散通量見表2，NH+4、NO-3擴散通量分別在0.068-0.688 mg·m-2·d-1和-0.013-0.035 mg·m-2·d-1之間。氮素界面擴散以NH+4為主，所在區(qū)域的擴散通量（F）基本為正值，沉積物向上覆水?dāng)U散較為明顯，整個區(qū)域的沉積物表現(xiàn)為“源”。

各小微水體沉積物—水界面NH+4、NO-3擴散通量的季節(jié)特征如圖3所示，各點之間沉積物-水界面的NH+4、NO-3擴散通量在不同季節(jié)之間存在顯著差異，總體呈現(xiàn)出夏季高、冬季低的趨勢。這一趨勢與袁軼君、盧昶雨等在太湖、鄱陽湖等大型水體中的研究結(jié)果一致。這種季節(jié)性變化主要的影響因素是溫度，在夏季，水溫相對較高，沉積物中的NH+4、NO-3分子運動加快，增強了沉積物—水界面的NH+4、NO-3擴散速率。同時，溫度的上升還有利于微生物將有機氮轉(zhuǎn)化為無機氮，進一步促進了沉積物—水界面NH+4、NO-3的擴散。此外，由于小微水體自身的水位較淺，夏季頻繁的降雨對沉積物的擾動較強，也可能促進沉積物—水界面NH+4、NO-3擴散。

各采樣點之間的沉積物—水界面NH+4、NO-3擴散通量也存在顯著差異，其主要原因是各點水體—沉積物的氮素本底濃度存在差異。小微水體沉積物—水界面的氮素擴散以NH+4為主，沉積物中NH+4的濃度越高，其在沉積物—水界面的擴散通量通常也越大。這是因為高濃度的NH+4為擴散過程提供了更多的物質(zhì)來源，從而擴散作用更為強烈。

2.3 沉積物-水界面NH+4、NO-3沉降通量與內(nèi)源釋放通量

通過沉降收集裝置，基于公式（3）估算出各采樣點的NH+4、NO-3沉降通量。如圖4所示，沉降通量受上覆水NH+4濃度影響，無顯著季節(jié)性差異。但是5個采樣點之間的NH+4沉降通量具有顯著性差異，其中Q1點的NH+4沉降通量最大，高達25.82 mg·m-2·d-1。這主要是由于Q1點為養(yǎng)殖魚塘點，魚苗會在每年3月左右開始投放，4-10月是魚苗的生長發(fā)育期，在這期間魚類食量增大、活動頻繁，漁民會定時投喂大量魚食餌料，日常投喂的魚食餌料中含有大量的有機氮，這些有機氮沒有被魚群及時攝入便沉降到了塘底，并在微生物的作用下礦化為NH+4。此外魚群每日也會排泄出大量的糞便沉降到底部，糞便中也含有大量的NH+4，這也是導(dǎo)致Q1點NH+4沉降通量大的原因之一。

基于公式（4），通過上述沉降通量與擴散通量進一步計算了沉積物—水界面NH+4、NO-3內(nèi)源釋放通量。如圖5所示，各點的氮素內(nèi)源釋放以NH+4為主。各點NH+4內(nèi)源釋放通量呈現(xiàn)夏季高、冬季低的趨勢，但總體差異并不顯著。而各采樣點NH+4內(nèi)源釋放通量之間存在顯著性差異，各點NH+4內(nèi)源釋放通量主要受其沉積物中NH+4濃度的影響。

由于沉積物-水界面NH+4擴散、沉降、內(nèi)源釋放通量都遠大于NO-3，因此本研究主要對NH+4的內(nèi)源釋放量和擴散進行水質(zhì)負荷分析。由表3可知，各個采樣點沉積物NH+4內(nèi)源釋放對上覆水體的氮負荷為223.99-9 630.62 kg·km-2·a-1，沉積物釋放對水質(zhì)的貢獻占比介于25.40%-415.17%。其中Q1點的貢獻超過100%，主要是由于該點受高強度人類活動影響，長期的投喂魚食餌料導(dǎo)致其沉降量很大，同時由于沉積物NH+4濃度較高，高密度魚群活動促進了沉積物再懸浮，導(dǎo)致其釋放量也很大。由于氮沉降量與釋放量都遠大于上覆水體中的氮，導(dǎo)致其貢獻值大于100%。但是，沉積物向上覆水?dāng)U散作用（即沉降過程與釋放過程的綜合作用）對上覆水體的負荷為50.75-311.52kg·km-2·a-1，其對水質(zhì)的貢獻占比介于2.07%-15.10%。

對小微水體沉積物—水界面NH+4擴散通量、沉降通量、內(nèi)源釋放通量及各水理化性質(zhì)指標(biāo)進行相關(guān)性分析。如圖6所示，沉積物—水界面NH+4擴散通量與沉積物與上覆水NH+4濃度差顯著正相關(guān)，與pH、D0顯著負相關(guān)，這與太湖、洞庭湖、西湖等水體的研究結(jié)果一致。NH+4沉降通量與水體流速、NH+4擴散通量顯著正相關(guān)，這與胡開明等、逄勇等的研究一致，水體流速增大能夠促進沉積物再懸浮，NH+4隨之?dāng)U散和沉降。NH+4內(nèi)源釋放通量與沉積物和上覆水NH+4濃度差、水體流速顯著正相關(guān)。

3 討論

太湖地區(qū)小微水體氮素內(nèi)源釋放主要以NH+4釋放為主，這與河流、湖泊一致。外源進入水體的氮大多為無機氮，這些無機氮隨時間在沉積物中不斷積蓄，發(fā)生同化作用，轉(zhuǎn)化成有機氮。沉積物中大量的有機氮在微生物的作用下又被礦化為NH+4，而NH+4帶正電荷，易被吸附在沉積物中，部分NH+4通過硝化作用被氧化為NO-3，NO-3易通過反硝化和淋溶等途徑損失。因此，小微水體氮素內(nèi)源釋放主要以NH+4為主。

氮素在水體中的沉降和釋放通量主要受氮素在水體中的運移過程影響。濃度差決定了氮素在水體與沉積物界面的運移方向。當(dāng)水體中的氮素濃度高于沉積物中的濃度時，氮素會向沉積物中遷移，導(dǎo)致沉積物中的氮素濃度增加，而水體中的濃度減小，相反，如果沉積物中的氮素濃度高于水體中的濃度，沉積物中的氮素會向水體中釋放，導(dǎo)致水體中的氮素濃度增加，而沉積物中的濃度減小。此外，流速決定了氮素在水體中的運移速度。水體流速越快，氮素在水體與沉積物界面中的運移速度也越快，因為快速流動的水體會促進氮素的遷移和擴散，從而加快氮素在水體和沉積物之間的交換速率。

小微水體沉積物—水界面NH+4擴散通量為50.75-311.52 kg·km-2·a-1，其對水質(zhì)貢獻為2.07%- 15.10%。與大型水體如太湖相比，小微水體的擴散強度和水質(zhì)貢獻均相對較小。例如，太湖擴散通量為1149-7 773t·km-2·a-1，其對水質(zhì)的貢獻占比高達25.7%。這種差異主要源于小微水體沉積物氮濃度以及沉積物與水體氮濃度差顯著小于太湖。

與河流相比，太湖地區(qū)小微水體的內(nèi)源釋放通量也相對較低。以太倉市的釜陽河為例，其NH+4擴散通量高達321-9 665 t·km-2·a-1。這種差異主要源于流速和水動力條件差異。河流通常具有較大的流速，這種高流速對沉積物產(chǎn)生了強烈的擾動作用。這種擾動不僅使得沉積物中的營養(yǎng)鹽更容易釋放到水體中，還促進了水體與沉積物之間的物質(zhì)交換，從而加劇了內(nèi)源釋放過程。相比之下，小微水體往往流速較慢甚至靜止不動。這種緩慢的水流條件限制了沉積物與水體之間的物質(zhì)交換，使得沉積物中的內(nèi)源釋放過程相對較弱。

雖然小微水體內(nèi)源釋放的擴散強度和水質(zhì)貢獻較小，但由于小微水體數(shù)量眾多、分布廣泛，其累積效應(yīng)和對整體水環(huán)境質(zhì)量的潛在影響不容忽視。特別是在農(nóng)村地區(qū)和城市化進程中，小微水體往往成為生活污水、農(nóng)業(yè)廢水等污染物的集中排放地，加劇了水體富營養(yǎng)化的風(fēng)險，容易形成黑臭水體，阻礙了鄉(xiāng)村振興戰(zhàn)略的實施。因此，在水體富營養(yǎng)化治理中，除了關(guān)注大型水體，還應(yīng)將小微水體納入重點管理和控制范疇。通過加強小微水體的水質(zhì)監(jiān)測、污染源控制、生態(tài)修復(fù)等措施，可以有效減少其內(nèi)源釋放，降低水生態(tài)環(huán)境惡化和富營養(yǎng)化的風(fēng)險。

4 結(jié)論

（1）太湖流域小微水體氮素內(nèi)源釋放以NH+4-N為主，其內(nèi)源釋放通量范圍在223.99-9 630.62 kg·km-2·a-1之間。

（2）太湖流域小微水體內(nèi)源釋放現(xiàn)象顯著，其沉積物—水界面NH+4擴散通量范圍為50.75-311.52 kg·km-2·a-1，NH+4擴散對水質(zhì)的貢獻介于2.07% - 15.10%。

（3）水體流速和沉積物與上覆水之間的NH+4濃度差是影響小微水體內(nèi)源釋放的主要因素。

（責(zé)任編輯：葉飛）

基金項目：國家自然科學(xué)基金項目（42177401）