王會(huì)，劉濱輝

(東北林業(yè)大學(xué)，黑龍江 哈爾濱 150040)

黑龍江省近47年生長(zhǎng)季淺層地溫變化特征

王會(huì)，劉濱輝

(東北林業(yè)大學(xué)，黑龍江 哈爾濱 150040)

利用黑龍江省15個(gè)國(guó)家氣象站1965—2011年生長(zhǎng)季(5—9月)0、5、10、15、20 cm處逐日平均地溫資料，計(jì)算黑龍江省各站點(diǎn)生長(zhǎng)季平均地溫的氣候傾向率和年際變化，分析黑龍江省47年來生長(zhǎng)季平均地溫隨緯度和土層深度的變化特征。結(jié)果表明：生長(zhǎng)季各月份平均地溫的垂直變化形式基本一致，且隨著土層深度增加地溫逐漸降低，在0～5 cm熱量衰減較快，而在5～10 cm、10～15 cm、15～20 cm的熱量衰減速率基本相同,衰減相對(duì)較慢；生長(zhǎng)季各月份各層平均地溫均呈升高趨勢(shì)，0 cm處平均地溫氣候傾向率在5個(gè)深度中最高，5 cm、10 cm、15 cm、20 cm平均地溫升溫率相差不大；5—7月平均地溫的升溫趨勢(shì)比8—9月明顯，8—9月的升溫對(duì)整個(gè)生長(zhǎng)季的升溫趨勢(shì)貢獻(xiàn)較??；地溫與空氣溫度呈正相關(guān)，與降水呈負(fù)相關(guān)。

黑龍江??；生長(zhǎng)季；淺層地溫；變化特征

近幾十年來全球氣候變化已成為世界各國(guó)研究的熱點(diǎn)，地表溫度與人類的生存息息相關(guān)，是氣候變化研究的重要內(nèi)容之一。IPCC報(bào)告指出在1906—2005年100年間，地表溫度升高了0.74 ℃，且這種升溫主要集中在近50年，在中高緯度地區(qū)尤為突出[1]。氣候和氣候變化與人類社會(huì)、經(jīng)濟(jì)生活息息相關(guān)，全球變暖對(duì)世界和我國(guó)的各類生態(tài)系統(tǒng)及社會(huì)經(jīng)濟(jì)產(chǎn)生了并將持續(xù)產(chǎn)生重大影響[2]。相對(duì)比其他的生態(tài)系統(tǒng)，森林生態(tài)系統(tǒng)受全球氣候變暖的影響更加明顯[3]，且其對(duì)土壤溫度和水分變化的響應(yīng)敏感[4]。 闊葉紅松林土壤呼吸對(duì)土壤溫度的反應(yīng)敏感, 不論何種土壤含水量, 土壤呼吸速率均隨溫度的增加而增高，土壤的呼吸速率升高促進(jìn)植物根部生長(zhǎng)加速[5]。土壤溫度主要通過下墊面反射率、土壤濕度及土壤熱儲(chǔ)量三個(gè)方面來影響森林生態(tài)系統(tǒng)的變化。因此,地溫變化規(guī)律的研究有助于全面了解氣候變化、生態(tài)系統(tǒng)演替、植物生長(zhǎng)、土壤碳吸收等多方面的影響。部分學(xué)者對(duì)氣候變化的分析主要集中于對(duì)氣溫和降水變化規(guī)律的研究，對(duì)于地溫的研究相對(duì)較少。Hu等分析了美國(guó)北部和西北部10 cm地溫的變化，發(fā)現(xiàn)大多數(shù)站點(diǎn)地溫均呈現(xiàn)升溫趨勢(shì)[6]。我國(guó)有關(guān)地溫方面的研究始于20世紀(jì)90年代, 杜軍等對(duì)我國(guó)淺層地溫進(jìn)行了研究，結(jié)果表明淺層地溫呈現(xiàn)明顯升高趨勢(shì)[7，8]。

本文主要研究黑龍江省淺層地溫的變化特征，黑龍江省地處中高緯度，是氣候變化的敏感區(qū)域，有自己的氣候變化特點(diǎn)，季節(jié)性凍融更增加了本區(qū)地溫研究的復(fù)雜性[9]。黑龍江省森林覆蓋率達(dá)43.6%，地溫變化規(guī)律的研究對(duì)深入研究森林的土壤呼吸和植物初級(jí)生產(chǎn)力變化以及碳循環(huán)過程均有意義[10]。筆者基于黑龍江省1965—2011年5—9月0 cm、5 cm、10 cm、15 cm、20 cm處逐日平均地溫資料，分析生長(zhǎng)季各層地溫的空間分布和時(shí)間演變規(guī)律，以期揭示地溫對(duì)森林生態(tài)環(huán)境及植物生長(zhǎng)的影響。

1 材料與方法

1.1 數(shù)據(jù)來源與處理

研究數(shù)據(jù)選取黑龍江省15個(gè)國(guó)家基準(zhǔn)氣象觀測(cè)站1965—2011年0、5、10、15和20 cm處生長(zhǎng)季逐日平均地溫資料。大部分氣象站是從1950年初期開始觀測(cè)的，由于觀測(cè)初期檢測(cè)的數(shù)據(jù)缺失較嚴(yán)重，因此該研究?jī)H采用1965—2011年的氣象數(shù)據(jù)進(jìn)行研究，又由于黑龍江省冬半年為凍土，5 cm、10 cm、15 cm、20 cm地溫資料從10月到次年4月幾乎全部缺失，所以統(tǒng)一采用5—9月的地溫?cái)?shù)據(jù)。另外在本研究中，如果一個(gè)站點(diǎn)數(shù)據(jù)連續(xù)缺失超過30 d，則刪除該站點(diǎn)，如漠河站10 cm和15 cm的地溫?cái)?shù)據(jù)在1960—1961年全部缺失，因此本研究排除漠河站。剩下15個(gè)站點(diǎn)5—9月的地溫?cái)?shù)據(jù)從1965年開始比較完整，缺失數(shù)據(jù)占總觀測(cè)值的0.02%。所選取的15個(gè)氣象站點(diǎn)在黑龍江省的空間分布合理，有較好的代表性(圖1)。由于檢測(cè)數(shù)據(jù)較全，只是在連續(xù)缺失在7d以內(nèi)的情況則使用簡(jiǎn)單的線性插補(bǔ)法，且將插補(bǔ)過的數(shù)據(jù)與原始數(shù)據(jù)的計(jì)算結(jié)果進(jìn)行了比對(duì)，差異不顯著[11]。

1.2 研究方法

為了揭示黑龍江省生長(zhǎng)季淺層地溫的平均狀態(tài)，根據(jù)日平均地溫資料，得到黑龍江省各個(gè)站點(diǎn)各層深度1965—2011年的多年生長(zhǎng)季平均溫度，用ARCVIEW空間分析功能插入等值線，插入等值線有助于對(duì)地溫空間變化的分析與評(píng)價(jià)。

將黑龍江省所有站點(diǎn)的地溫月平均值進(jìn)行月份平均，X軸月平均溫度，Y軸不同深度土層用來述地溫在生長(zhǎng)季垂直變化趨勢(shì)。運(yùn)用多年淺層地溫?cái)?shù)據(jù)進(jìn)行距平，對(duì)時(shí)間序列進(jìn)行了九點(diǎn)二次平滑處理，以分析整個(gè)黑龍江省生長(zhǎng)季地溫的年際變化，揭示黑龍江省區(qū)域地溫隨時(shí)間的變化特征，然后應(yīng)用氣候傾向率和相關(guān)分析來尋找地溫與大氣溫度和降水的關(guān)系。

2 結(jié)果與分析

2.1 地溫的空間分布特征

2.1.1 地溫隨緯度的變化特征 黑龍江省位于我國(guó)最北部，南北跨越中溫帶和寒溫帶，地形復(fù)雜，溫差較大，緯度位置是造成溫度差異的主要原因[9]。由于排除漠河站，黑龍江省最北部的地溫分布趨勢(shì)并不準(zhǔn)確。因此本文重點(diǎn)討論呼瑪站以南黑龍江省大部分地區(qū)地溫的空間分布。之前的學(xué)者對(duì)黑龍江省1965—2011年地表溫度的分析得出，隨著緯度的升高，由南至北多年平均地溫逐漸降低，等溫線基本與緯線平行，南部地區(qū)等溫線產(chǎn)生沿經(jīng)線方向的分布特征[11]。本文通過對(duì)0、5、10、15、20 cm生長(zhǎng)季平均地溫的空間分布分析，得出相似的結(jié)果：黑龍江省各層生長(zhǎng)季平均地溫的空間分布模式基本一致，隨著緯度升高地溫降低，大部分地區(qū)等溫線與緯線平行，而東部等溫線也產(chǎn)生沿經(jīng)線方向分布的特點(diǎn)，東部地區(qū)地溫要低于中部和西部。0、5、10、15、20 cm地溫最低的站點(diǎn)都是黑河站，分別為18.29 ℃、16.55 ℃、15.84 ℃、15.11 ℃、14.47 ℃；溫度最高的站點(diǎn)都是哈爾濱站，分別為22.00 ℃、20.17 ℃、19.52 ℃、18.93 ℃、18.45 ℃；南北溫差分別為3.71 ℃、3.62 ℃、3.68 ℃、3.82 ℃和3.98 ℃?？梢姡麄€(gè)黑龍江各層生長(zhǎng)季平均地溫的南北溫度差異并不大。因此，緯度不是影響黑龍江省地溫空間差異的唯一因素。

2.1.2 地溫的垂直變化特征 由于土壤熱量衰減的性質(zhì)，在土壤淺層，隨著深度增加地溫逐漸降低。圖2給出0 cm、5 cm、10 cm、15 cm、20 cm的黑龍江省生長(zhǎng)季月平均地溫的垂直分布圖。各個(gè)站點(diǎn)生長(zhǎng)季月平均地溫的垂直變化形式基本一致，0 cm的溫度最高，20 cm的溫度最低。圖中相鄰兩層之間線段斜率的大小可以顯示出熱量衰減的快慢，在0～5 cm熱量衰減較快，因此0 cm地溫和5 cm地溫的溫差較大。這是因?yàn)橥寥辣韺咏咏髿?，土壤與大氣之間熱量交換頻繁，有一部分熱量散發(fā)到大氣中。5～10 cm、10～15 cm、15～20 cm的熱量衰減速率基本相同，說明5 cm以下地溫受大氣的影響較小。在5月熱量衰減的幅度較大，0 cm與20 cm相差6.12度，地表的熱量迅速向下傳遞，6—8月份每層溫度變化不大，9月份每層溫度幾乎沒有變化，甚至到20 cm時(shí)溫度還有上升的趨勢(shì)，主要是因?yàn)樘鞖廪D(zhuǎn)涼、日照時(shí)數(shù)減少，地表的熱量降低，深層的熱量有向上傳遞的趨勢(shì)。

2.2 生長(zhǎng)季地溫的年際變化特征

圖4為黑龍江省各層生長(zhǎng)季平均地溫的年際變化。為了更好地了解數(shù)據(jù)的規(guī)律性，將多年淺層地溫?cái)?shù)據(jù)進(jìn)行距平，對(duì)時(shí)間序列進(jìn)行了九點(diǎn)二次平滑處理。每層地溫的升溫幅度依次為0.048、0.040、0.041、0.044、0.043 ℃a-1。淺層地溫年際變化規(guī)律圖中可以看出，在1965—2011年期間都呈升溫趨勢(shì)，在1965—1995年之間地溫的上升幅度較平穩(wěn)，在1995—2011年之間溫度迅速上升，由此可以將淺層土壤溫度的變化劃分為兩個(gè)時(shí)期：平穩(wěn)期和上升期。在五個(gè)深度的地溫中0 cm地溫緊靠大氣層，受溫度影響最強(qiáng)，升溫幅度最大，變化最慢的是20 cm。各層地溫的年際變化均具有很好的一致性，這說明地溫變化主要受地表熱量條件控制，隨著深度加深并沒有其他熱量來源。因此任何深度的地溫都能用來印證及推測(cè)地表的氣候條件。

2.3 不同深度地溫與氣溫、降水的關(guān)系特征

表1 不同層次土壤溫度氣溫、降水相關(guān)分析

注：★表示通過90%的顯著性檢驗(yàn)

據(jù)對(duì)黑龍江省1965—2011年地表溫度年及其季節(jié)空間變化的分析，季風(fēng)環(huán)流對(duì)黑龍江省氣候有重要的影響。冬季風(fēng)加劇了空氣降溫，使冬季氣候嚴(yán)寒，降水稀少；夏季風(fēng)使黑龍江省氣溫普遍升高，降水主要來源于夏季風(fēng)[11]。由表1可知，不同層次的土壤溫度之間的相關(guān)性較好，所有土壤溫度相關(guān)性檢驗(yàn)均通過90%的顯著性檢驗(yàn)。5～20 cm各級(jí)土層的相關(guān)系數(shù)均超過0.96。氣溫與各層地溫之間均是正相關(guān)，在其中0cm與氣溫的相關(guān)性最強(qiáng)，是因?yàn)樵?—9月太陽(yáng)光直射到地表，地表折射到大氣中的熱量較多，土壤下層受太陽(yáng)輻射逐漸減少與氣溫的相關(guān)性也逐漸降低。表中所示降水與各層地溫均是負(fù)相關(guān)，而且相關(guān)系數(shù)較低，也未通過90%的顯著性檢驗(yàn)。黑龍江5—9月主要受夏季風(fēng)所影響。緯度和氣溫升高的雙重影響，是生長(zhǎng)季平均地溫隨緯度升高而降低但南北溫差不大的主要原因。至于東部地區(qū)地溫低于中部和西部，是因?yàn)橄募撅L(fēng)使得降水量自東向西逐漸減少，而降水和地溫呈負(fù)相關(guān)關(guān)系，東部地區(qū)較多的降水使得地溫較低。

2.4 各站點(diǎn)地溫變化特征

表2 按緯度排列的各站點(diǎn)1965～2011年生長(zhǎng)季平均地溫的氣候傾向率 ℃a-1

表2 按緯度排列的各站點(diǎn)1965～2011年生長(zhǎng)季平均地溫的氣候傾向率 ℃a-1

站點(diǎn)緯度0cm5cm10cm15cm20cm呼瑪51.430.069★0.044★0.041★0.037★0.040★黑河50.150.077★0.064★0.061★0.060★0.057★嫩江49.10.073★0.058★0.057★0.059★0.062★克山48.030.060★0.065★0.062★0.064★0.059★海倫47.260.061★0.052★0.055★0.059★0.057★齊齊哈爾47.230.065★0.067★0.063★0.073★0.072★富錦47.140.0490.031★0.036★0.039★0.038★佳木斯46.490.072★0.052★0.049★0.050★0.045★安達(dá)46.230.053★0.040★0.039★0.038★0.026通河45.580.049★0.054★0.049★0.052★0.056★哈爾濱45.450.050★0.044★0.043★0.043★0.040★雞西45.170.073★0.051★0.055★0.058★0.055★尚志45.130.066★0.046★0.050★0.053★0.049★牡丹江44.340.064★0.050★0.049★0.046★0.043★綏芬河44.230.047★0.0280.0240.0220.019

注：★表示通過90%的顯著性檢驗(yàn)

由表2可知， 15個(gè)站點(diǎn)5個(gè)深度的生長(zhǎng)季平均地溫都呈現(xiàn)升高趨勢(shì)，升幅為0.019～0.077 ℃a-1，大部分通過90%顯著性檢驗(yàn)，且北部站點(diǎn)的升溫率要高于南部站點(diǎn)。0 cm生長(zhǎng)季平均地溫變化趨勢(shì)，富錦、通河、哈爾濱及綏芬河的升溫趨勢(shì)不明顯，氣候傾向率在0.047～0.050 ℃a-1之間，其余11個(gè)站的升溫趨勢(shì)顯著，增溫幅度為0.053～0.077 ℃a-1；5 cm生長(zhǎng)季平均地溫以黑河、嫩江、克山、齊齊哈爾和通河的增幅最大，為0.054～0.067 ℃a-1，其余站點(diǎn)氣候傾向率在0.028～0.052 ℃a-1之間；10 cm生長(zhǎng)季平均地溫，除安達(dá)、富錦及綏芬河外，其余各站點(diǎn)均表現(xiàn)顯著的升溫趨勢(shì)，升溫率為0.041～0.063 ℃a-1；15 cm生長(zhǎng)季平均地溫在富錦、呼瑪及綏芬河升溫趨勢(shì)不明顯，其余站點(diǎn)升溫趨勢(shì)顯著，氣候傾向率為0.043～0.073 ℃a-1；20 cm生長(zhǎng)季平均地溫除富錦、安達(dá)和綏芬河升溫趨勢(shì)不明顯外，其他站點(diǎn)以0.040～0.072 ℃a-1的幅度升高。

可見，5個(gè)深度的大部分站點(diǎn)均出現(xiàn)顯著的升溫趨勢(shì)。其中黑河、嫩江、海倫、安達(dá)、通河、哈爾濱、雞西、尚志、牡丹江9個(gè)站，在5個(gè)深度均出現(xiàn)顯著升溫趨勢(shì)，綏芬河和富錦站5個(gè)深度升溫率都比較低，呼瑪和佳木斯站只有0 cm地溫升高趨勢(shì)明顯，而齊齊哈爾和克山站分別是0 cm和5 cm平均地溫的升高趨勢(shì)不顯著。

3 討論與結(jié)論

3.1 討論

本項(xiàng)研究表明，地溫的年際變化具有很強(qiáng)的一致性，其主要受地表熱量控制。在過去的47年中，黑龍江省生長(zhǎng)季5—9月地溫隨氣溫升高而下降，東部地區(qū)比中部和西部地區(qū)下降的快，但其上升較快，降低較慢，整體呈上升趨勢(shì)。王連喜等研究結(jié)果顯示地溫升高對(duì)氣候變暖有直接影響[12]。而肖輝的研究表明森林的進(jìn)化與衰退東部比西部快，特別是黑河地區(qū)，這可能是地溫東西部變化不一致導(dǎo)致的[13]。地溫上升會(huì)影響林地植物殘?bào)w的分解量、分解強(qiáng)度和泥炭物質(zhì)聚集，地溫越高，越有利于有機(jī)物質(zhì)分解，土壤中有機(jī)物質(zhì)含量增高，促進(jìn)樹木根系的生長(zhǎng)，提高產(chǎn)量[10]。

3.2 結(jié)論

利用黑龍江省1965—2011年15個(gè)國(guó)家氣象站淺層逐日地溫、大氣溫度及降水資料，分析黑龍江省近47年來淺層地溫的時(shí)空變化特征，初步得出如下結(jié)論：

(1)生長(zhǎng)季不同深度平均地溫均隨著緯度的升高而降低，大部分地區(qū)等溫線與緯線平行，東部地區(qū)等溫線產(chǎn)生沿經(jīng)線方向分布的特點(diǎn)，但整個(gè)黑龍江省生長(zhǎng)季平均地溫的南北差異不大。

(2)各個(gè)站點(diǎn)生長(zhǎng)季平均地溫的垂直變化形式基本一致，在0～5 cm深度熱量衰減較快，5～10 cm、10～15 cm、15～20 cm的熱量衰減速率基本相同。在47年的地溫年際變化趨勢(shì)中看出，整體呈上升趨勢(shì)，在1995年前地溫平穩(wěn)上升，在1995年后迅速上升。

(3)黑龍江省地溫在7月份達(dá)到最大，5月到7月升溫過程較快，而7月到9月的降溫過程緩慢。

(4)黑龍江省淺層地溫之間的相關(guān)性較好，與大氣溫度呈正相關(guān)，與降水呈負(fù)相關(guān)。

(5)在1965—2011年47年間，黑龍江省15個(gè)站生長(zhǎng)季5個(gè)深度的平均地溫都呈現(xiàn)升高趨勢(shì)，且大部分站點(diǎn)的升溫趨勢(shì)明顯，有9個(gè)站在5個(gè)深度均出現(xiàn)顯著的升溫。其中0 cm地溫氣候傾向率在5個(gè)深度中最高，為0.077 ℃y-1；5 cm、10 cm、15 cm、20 cm深度平均地溫升溫率相差不大。全省各層地溫月平均地溫都呈現(xiàn)升溫趨勢(shì)，整體升溫趨勢(shì)為5月>7月>6月>9月>8月， 8月、9月的升溫對(duì)整個(gè)生長(zhǎng)季升溫趨勢(shì)的貢獻(xiàn)較小。

[1] IPCC. Climate change 2007: Synthesis report[R]. Oslo: Intergovernmental Panel on Climate Change, 2007

[2] 秦大河.氣候變化:科學(xué)、影響和對(duì)策[M]//中國(guó)氣象年鑒.北京:氣象出版社,2004:119

[3] Lindroth A, Grelle A, Moren A S. Long-term measurements of boreal forest carbon balance reveal large temperature sensitivity. Global Change Biology, 1998, 4(4):443-450

[4] Boone RD, Nadelhoffer KJ, Canary JD, KayeJP. Roots exert a strong influence on the temperature sensitivity of soil respiration. Nature, 1998, 396(6711):570-572

[5] 王淼,姬蘭柱,李秋榮，等.土壤溫度和水分對(duì)長(zhǎng)白山不同森林類型土壤呼吸的影響[J].應(yīng)用生態(tài)學(xué)報(bào),2003,14(8):1234-1238

[6] Hu Q，F(xiàn)eng S． A dalily soil temperature dataset and soil temperature cilmatology of the contiguous united states[J]． Journal of applied meteorology，2003，42(8) : 1139-1156

[7] Qi Hu, S. Feng How have soil temperatures been affected by the surface temperature and precipitation in the Eurasian continent? American Geophysical Union 2005

[8] 杜軍,李春,廖健,等.拉薩近45年淺層地溫的變化特征[J].干旱區(qū)地理,2007,30(6):826-831

[7] 代玉田.近49年德州淺層地溫的變化特征分析[J].中國(guó)農(nóng)業(yè)氣象,2011，32(S1):1-4

[8] 胡軍,杜軍,邊多,等.西藏地溫的年際和年代際變化[J].地理學(xué)報(bào),2007(9):925-934

[9] 李帥，王萍，陳莉，等．黑龍江省春季淺層(0～20 cm)地溫變化特征及預(yù)報(bào)[J]．冰川凍土，2014，36(1):55-62

[10] 高舉明,張一平,于貴瑞，等.西雙版納熱帶季節(jié)雨林地溫特征[J].生態(tài)學(xué)雜志,2008,27(6):880-887

[11] 袁文濤,劉濱輝,劉燕玲,等.黑龍江省1960-2000年0cm地溫的變化趨勢(shì)[J].東北林業(yè)大學(xué)學(xué)報(bào),2011，39(2):67-71

[12] 王連喜,陳懷亮,李琪,等.植物物候與氣候研究進(jìn)展[J].生態(tài)學(xué)報(bào),2010,30(2):0447-0454

[13] 肖輝林.土壤溫度上升與森林衰退[J].熱帶亞熱帶土壤科學(xué)，1995,4(4):246-249

Change Characteristics of Surface Temperature in Growing Season in Heilongjiang Province in Recent 47 Years

Wang Hui , Liu Binhui

(Northeast Forestry University, Harbin 150040,China)

Inclination rate & interannual variation of the climate of the average earth temperature in the growing season of the weather stations in Heilongjiang Province were calculated by using data of daily mean earth temperature at 0, 5, 10, 15 & 20 cm in growing season (from May to September) from 1965 to 2011 in 15 national weather stations of Heilongjiang Province;the change characteristics of the mean earth temperature in the growing season with latitude and depth in the past 47 years in Heilongjiang Province were analyzed.Result shows that the form of various levels for the mean earth temperature in the growing season in every month is basically identical, and the earth temperature decrease with the increase of the soil depth;heat decay faster at 0-5 cm layer;heat decay rate are basically the same at 5-10, 10-15,15-20 cm. The decay is relatively slow;the mean earth temperature shows rising trend in every month ;the inclination rate of climate at 0 cm is the highest among the five depth; temperature rise rate of the mean earth temperature is little at 5,10,15 &20 cm layer;temperature rise trend in May, June & July is significant than that in August & September;temperature rise in August and September contributes little to the temperature rise trend in the whole growing season;the earth temperature is positively correlated with the air temperature and negatively correlated with precipitation.

Heilongjiang Province;growth season; shallow earth temperature; change characteristics

1005-5215(2016)12-0013-05

2016-11-01

黑龍江省科技攻關(guān)項(xiàng)目(GA06B302-7)；北京林業(yè)大學(xué)985平臺(tái)建設(shè)項(xiàng)目(2011BAD32B05)

王會(huì)(1987- )，女，碩士研究生，現(xiàn)從事水土保持與荒漠化防治及氣候變化研究.

劉濱輝(1970- )，男，副教授，現(xiàn)從事水土保持與荒漠化防治及氣候變化研究，Email:947441282@qq.com

P468.021

A

10.13601/j.issn.1005-5215.2016.12.004