融合空間統(tǒng)計(jì)特征的長三角生態(tài)空間供需平衡研究

陳藝華 ， 陳振杰 ， 李飛雪 ， 黃秋昊

(南京大學(xué) 地理與海洋科學(xué)學(xué)院，南京 210023)

0 引言

生態(tài)空間是主要提供生態(tài)服務(wù)或生態(tài)產(chǎn)品的國土空間[1-2]。生態(tài)空間質(zhì)量直接影響其提供人類日常生活必要的食物、能源、服務(wù)的能力[3-4]。生態(tài)需求是人類社會(huì)希望獲得的生態(tài)系統(tǒng)服務(wù)以及通過一定方式(如金錢、時(shí)間、距離成本等)能夠或想要獲取的生態(tài)系統(tǒng)服務(wù)[5-8]。隨著經(jīng)濟(jì)全球化發(fā)展，快速的城市化對(duì)生態(tài)空間產(chǎn)生巨大影響，經(jīng)濟(jì)發(fā)展、城市擴(kuò)張、人口激增帶來的后果之一是日益增長的資源需求與不斷下降的生態(tài)空間質(zhì)量之間的嚴(yán)重沖突[9-10]，這給區(qū)域可持續(xù)發(fā)展和生態(tài)文明建設(shè)帶來了嚴(yán)峻挑戰(zhàn)。

隨著城市化的發(fā)展和生活水平的提高，人們對(duì)高質(zhì)量生態(tài)空間的需求日漸強(qiáng)烈[11]，生態(tài)空間供需關(guān)系受到廣泛關(guān)注。有關(guān)生態(tài)空間質(zhì)量和需求的研究方法主要分三類：(1)指標(biāo)構(gòu)建方法是根據(jù)一定范圍內(nèi)提供的生態(tài)服務(wù)與人口、經(jīng)濟(jì)的比例來衡量供需匹配程度。如Y.Deng等[12]通過構(gòu)建休閑服務(wù)供給指數(shù)和社會(huì)需求指數(shù)來分析城市休閑服務(wù)供給與需求的空間錯(cuò)配；壽飛云等[13]通過土地開發(fā)強(qiáng)度、人口密度、GDP測(cè)算生態(tài)系統(tǒng)服務(wù)需求，分析長三角生態(tài)系統(tǒng)供求變化趨勢(shì)。(2)基于地理信息方法分析供需匹配。如K.Gupta等[14]利用緩沖區(qū)分析和網(wǎng)絡(luò)分析評(píng)估綠地可達(dá)性及綠地需求程度。(3)基于調(diào)查問卷、統(tǒng)計(jì)數(shù)據(jù)或依賴專家經(jīng)驗(yàn)、知識(shí)的經(jīng)驗(yàn)調(diào)查法。如B.Burkhard等[5]提出將土地利用類型與生態(tài)系統(tǒng)服務(wù)潛力、生態(tài)系統(tǒng)服務(wù)流和生態(tài)系統(tǒng)服務(wù)需求聯(lián)系起來的專家生態(tài)系統(tǒng)服務(wù)供求矩陣；K.Vierikko等[15]使用訪談和統(tǒng)計(jì)數(shù)據(jù)，通過綜合生態(tài)系統(tǒng)服務(wù)方法，分析了芬蘭城市湖泊生態(tài)系統(tǒng)娛樂供需關(guān)系。研究表明，不同區(qū)域的自然和社會(huì)經(jīng)濟(jì)狀況差異會(huì)影響其生態(tài)服務(wù)供給和需求，存在空間不平衡的問題[16-17]。

上述方法能夠一定程度上反映生態(tài)供求關(guān)系，但仍存在以下問題：(1)對(duì)人們的生態(tài)需求研究往往停留在人口和社會(huì)經(jīng)濟(jì)層面，較少考慮到交通可達(dá)性對(duì)人們享受生態(tài)空間的影響。(2)生態(tài)供給與需求的平衡因地而異，并隨經(jīng)濟(jì)社會(huì)發(fā)展而變化，難以綜合量化。

明晰區(qū)域的生態(tài)空間供需特征對(duì)生態(tài)空間科學(xué)管理、城鄉(xiāng)高質(zhì)量發(fā)展、人居環(huán)境改善意義重大。本研究選擇長三角為研究區(qū)域，構(gòu)建生態(tài)空間質(zhì)量評(píng)價(jià)體系，建立生態(tài)空間供需平衡指數(shù)，分析長三角生態(tài)空間質(zhì)量、供需平衡的時(shí)空演變特征，為長三角區(qū)域生態(tài)空間的合理管控提供支持。

1 研究區(qū)域與數(shù)據(jù)來源

1.1 研究區(qū)域

研究區(qū)域?yàn)椤堕L江三角洲區(qū)域一體化發(fā)展規(guī)劃綱要》中的長三角區(qū)域(圖1)。長三角是中國經(jīng)濟(jì)社會(huì)發(fā)展中最活躍、最具創(chuàng)新力、能源資源需求最高的地區(qū)之一。由于工業(yè)化、現(xiàn)代化建設(shè)的飛速發(fā)展，經(jīng)歷了“先污染后治理”之路，導(dǎo)致生態(tài)空間污染嚴(yán)重，生態(tài)環(huán)境問題逐漸加劇。隨著城鄉(xiāng)一體化發(fā)展等政策的相繼出臺(tái)，堅(jiān)持生態(tài)保護(hù)優(yōu)先、加強(qiáng)生態(tài)空間共保、推動(dòng)環(huán)境協(xié)同治理成為長三角未來生態(tài)環(huán)境管治的主要任務(wù)。

圖1 研究區(qū)域

1.2 數(shù)據(jù)來源

1995年和2015年土地利用、降水量、土壤侵蝕、GDP、人口和道路網(wǎng)數(shù)據(jù)均來源于資源環(huán)境科學(xué)與數(shù)據(jù)中心(http://www.resdc.cn/)；Landsat影像數(shù)據(jù)來源于地理空間數(shù)據(jù)云(http://www.gscloud.cn/)；社會(huì)經(jīng)濟(jì)數(shù)據(jù)來源于1996年和2016年上海市、江蘇省、浙江省、安徽省統(tǒng)計(jì)年鑒。

2 研究方法

2.1 生態(tài)空間質(zhì)量評(píng)價(jià)指標(biāo)

居民生活普遍關(guān)注生物多樣性、植被覆蓋率、供水等方面，有關(guān)領(lǐng)域?qū)＜姨貏e關(guān)注景觀格局和城市生態(tài)空間對(duì)人類的益處[18]。因此，本研究選擇以下指標(biāo)來評(píng)價(jià)生態(tài)空間質(zhì)量。(1)生物豐度指數(shù)(OI)。基于土地利用數(shù)據(jù)，由生物多樣性和生境質(zhì)量的平均值計(jì)算得到[18]。參照謝高地等[19]的研究結(jié)果，生物多樣性由土地利用類型對(duì)應(yīng)的生物多樣性服務(wù)價(jià)值計(jì)算得到。生境質(zhì)量以中國環(huán)境保護(hù)部的生境質(zhì)量指數(shù)為標(biāo)準(zhǔn)計(jì)算[20]。通過土地利用類型劃分，得到生物豐度指數(shù)(表1)。(2)植被指數(shù)(VI)。反映植被覆蓋程度，由Landsat影像計(jì)算一年中6—8月NDVI均值。(3)水網(wǎng)密度指數(shù)(WI)。由區(qū)域河流長度(RL)、區(qū)域水域面積(WA)和區(qū)域水資源(WR)三者均值計(jì)算得到[18]，RL和WA來自土地利用數(shù)據(jù)，WR來自降水量數(shù)據(jù)。(4)土地脅迫指數(shù)(LSI)。該指數(shù)參照中國環(huán)境保護(hù)部標(biāo)準(zhǔn)[20]，由區(qū)域重度侵蝕面積(SE)、區(qū)域中度侵蝕面積(ME)和區(qū)域建設(shè)用地面積(CA)加權(quán)計(jì)算得到，權(quán)重分別為0.5，0.3，0.2，SE和ME來自土壤侵蝕數(shù)據(jù)，CA來自土地利用數(shù)據(jù)。(5)景觀連通性指數(shù)(LCI)。該指數(shù)是對(duì)城市群生態(tài)空間中影響多種生境物種多樣性和生態(tài)安全格局的資源斑塊之間運(yùn)動(dòng)的促進(jìn)或阻礙程度的評(píng)價(jià)。采用Fragstats 4.2中的移動(dòng)窗口法,基于土地利用數(shù)據(jù)計(jì)算景觀指數(shù)PCI(patch cohesion index)。比較不同窗口大小的計(jì)算結(jié)果效果，選擇合適的窗口大小。通過比較,選擇9 000為移動(dòng)窗口大小。

表1 生物豐度指數(shù)

將所有的評(píng)價(jià)指標(biāo)歸一化為0～1，參考中國環(huán)境保護(hù)部標(biāo)準(zhǔn)[20]，生態(tài)空間質(zhì)量ESQ計(jì)算公式如下：

ESQ=0.30×OI+0.25×VI+0.15×WI+0.15×(1-LSI)+0.15×LCI 。

2.2 可達(dá)生態(tài)空間質(zhì)量

可達(dá)性表示居民克服距離、費(fèi)用等阻力到達(dá)一個(gè)目的地的愿望或能力[21]。根據(jù)不同等級(jí)道路、下墊面的通行能力計(jì)算相應(yīng)通行距離下的通行成本，以最小成本作為路徑選擇標(biāo)準(zhǔn)來表達(dá)個(gè)體通過不同出行方式克服空間阻隔到達(dá)目的地的難易程度?；诼肪W(wǎng)數(shù)據(jù)，根據(jù)交通方式的不同，對(duì)矢量路網(wǎng)賦予不同的速度值，即各道路車輛的平均行駛速度，計(jì)算通過每個(gè)格網(wǎng)單元的時(shí)間成本，格網(wǎng)單元大小為3 km×3 km(表2)。

表2 時(shí)間成本值統(tǒng)計(jì)

對(duì)得到的時(shí)間成本柵格，計(jì)算每個(gè)格網(wǎng)單元單位時(shí)間內(nèi)可到達(dá)的最遠(yuǎn)范圍，將該范圍內(nèi)的生態(tài)空間質(zhì)量和生態(tài)空間面積等權(quán)相加作為生態(tài)空間供給值。

2.3 生態(tài)需求

生態(tài)需求主要用居民日常生活(休閑游憩、短途出行等)中希望從生態(tài)用地獲得的相應(yīng)服務(wù)或產(chǎn)品的多少來衡量[22]。生態(tài)需求較高的地區(qū)通常集中在經(jīng)濟(jì)發(fā)展水平較高的區(qū)域[22]。本研究選取土地開發(fā)強(qiáng)度(區(qū)域建設(shè)用地面積占比)、GDP、人口密度3個(gè)指標(biāo)，建立3 km×3 km的格網(wǎng)，分別計(jì)算每個(gè)格網(wǎng)單元的生態(tài)環(huán)境壓力以反映不同的生態(tài)需求。由于各個(gè)指標(biāo)的數(shù)值差異很大，對(duì)土地開發(fā)強(qiáng)度進(jìn)行歸一化處理，對(duì)人口密度和GDP取對(duì)數(shù)后進(jìn)行歸一化處理，公式如下[13,22]：

ESD=X1×lgX2×lgX3。

式中：ESD為生態(tài)需求；X1，X2，X3分別為土地開發(fā)強(qiáng)度、GDP、人口密度。

2.4 生態(tài)空間供需平衡指數(shù)

由于不同地區(qū)的自然和社會(huì)經(jīng)濟(jì)條件的差異，其生態(tài)供需平衡狀態(tài)也會(huì)有所不同[23]。而且，隨著經(jīng)濟(jì)社會(huì)發(fā)展水平的提高，人們對(duì)自然環(huán)境、城市品質(zhì)的要求也會(huì)提高，相應(yīng)的生態(tài)空間供需平衡狀態(tài)也會(huì)發(fā)生變化[24]。因此，生態(tài)供給與需求的平衡因地而異、因時(shí)而異，難以定量化分析。從區(qū)域生態(tài)文明建設(shè)、實(shí)施成本的角度出發(fā)，最迫切的應(yīng)該是補(bǔ)“短板”，即優(yōu)先在生態(tài)空間供需缺口大的區(qū)域增加生態(tài)要素。基于上述分析，將區(qū)域內(nèi)各格網(wǎng)單元的可達(dá)生態(tài)空間質(zhì)量均值、生態(tài)需求均值作為供需平衡分析的基準(zhǔn)，借鑒統(tǒng)計(jì)分析方法，分析生態(tài)空間供需平衡狀況。

按照統(tǒng)計(jì)學(xué)的思想，各格網(wǎng)單元的可達(dá)生態(tài)空間質(zhì)量、生態(tài)需求將形成各自的統(tǒng)計(jì)分布，且二者的均值和方差均有差異。為了便于比較，將兩個(gè)分布標(biāo)準(zhǔn)化，并根據(jù)二者標(biāo)準(zhǔn)化值的差異，建立區(qū)域生態(tài)空間供需平衡指數(shù)，公式如下：

3 結(jié)果分析

3.1 長三角生態(tài)空間質(zhì)量格局

土地利用與生態(tài)空間質(zhì)量具有較強(qiáng)的空間相關(guān)性。按自然斷裂法將生態(tài)空間質(zhì)量評(píng)價(jià)結(jié)果分為5個(gè)等級(jí)(圖2)。長三角地區(qū)的生態(tài)空間質(zhì)量存在明顯的區(qū)域異質(zhì)性，上海生態(tài)空間數(shù)量和質(zhì)量均較低，江蘇處于中下水平，浙江和安徽相對(duì)較高。

圖2 生態(tài)空間質(zhì)量分布

從空間分布來看，長三角生態(tài)空間質(zhì)量呈現(xiàn)由北向南逐漸升高的趨勢(shì)。生態(tài)空間質(zhì)量高值區(qū)域主要位于浙南、皖西南丘陵區(qū)域，該區(qū)域自然環(huán)境良好，氣候適宜，降水充足，植被覆蓋度高，生態(tài)用地集中連片，能夠提供較高質(zhì)量的生態(tài)服務(wù)，且山區(qū)高程起伏大，與城區(qū)距離較遠(yuǎn)，受城市建設(shè)用地?cái)U(kuò)張等人類活動(dòng)的干擾較??；低值區(qū)主要包括上海、蘇南、皖中、浙北城鎮(zhèn)及周邊地區(qū)，該區(qū)域發(fā)展水平高，受人類活動(dòng)的干擾較大。由于南北經(jīng)濟(jì)發(fā)展的差異，蘇北和皖北地區(qū)的生態(tài)空間質(zhì)量較低，主要是該區(qū)域通過加強(qiáng)土地的開發(fā)利用和刺激經(jīng)濟(jì)增長來縮小區(qū)域間經(jīng)濟(jì)差異所導(dǎo)致。

從時(shí)空演變來看，1995—2015年生態(tài)空間質(zhì)量有較大的變化，等級(jí)為差和較差的區(qū)域占比均有上升，分別上升了3.52%和10.36%，等級(jí)為好、較好和一般的區(qū)域占比均下降，分別下降了3.94%，4.51%和5.42%，其中，較差區(qū)域有一定程度的增加。1995年生態(tài)空間質(zhì)量較差的區(qū)域主要分布在各大城市中心部分，2015年擴(kuò)大到城市周邊地區(qū)，并逐漸成連片狀。生態(tài)空間質(zhì)量下降的主要原因是社會(huì)經(jīng)濟(jì)的快速發(fā)展和城鎮(zhèn)化水平的提高以及建設(shè)用地大規(guī)模的擴(kuò)張。1995—2015年長三角耕地面積減少1.16萬km2，建設(shè)用地面積增加約1.33萬km2。其中，大部分耕地、部分林地和草地變?yōu)榻ㄔO(shè)用地，少部分轉(zhuǎn)換為水域。

生態(tài)空間質(zhì)量受城市化發(fā)展的影響顯著[18]?？紤]數(shù)據(jù)的可獲取性和統(tǒng)計(jì)尺度，選擇城市化進(jìn)程中社會(huì)經(jīng)濟(jì)相關(guān)因素。數(shù)據(jù)經(jīng)標(biāo)準(zhǔn)化處理，通過普通最小二乘法(OLS)模型剔除共線性變量，配合使用地理加權(quán)回歸(GWR)模型對(duì)城市化相關(guān)影響因素進(jìn)行分析(表3)。總體而言，不同城市化相關(guān)因素對(duì)不同城市生態(tài)空間質(zhì)量的影響差異較大。例如，教育水平和第三產(chǎn)業(yè)比重是影響上海和南京生態(tài)空間質(zhì)量的主要因素，第一產(chǎn)業(yè)比重是影響宿州和六安的主要因素。廢棄物處理水平對(duì)生態(tài)空間質(zhì)量的影響主要是正向影響，交通水平對(duì)生態(tài)空間質(zhì)量的影響主要是負(fù)向影響。

表3 城市化相關(guān)生態(tài)空間質(zhì)量影響因素系數(shù)

3.2 長三角可達(dá)生態(tài)空間質(zhì)量

研究表明，單次出行時(shí)間超過30 min時(shí)，人們普遍會(huì)感到焦慮和不耐煩[25]。本研究選擇計(jì)算30 min可達(dá)范圍內(nèi)生態(tài)空間質(zhì)量情況(圖3)?；诳蛇_(dá)性的生態(tài)空間質(zhì)量與原生態(tài)空間質(zhì)量評(píng)價(jià)結(jié)果相比，各等級(jí)連片分布，與路網(wǎng)相類似呈放射狀。1995—2015年，長三角地區(qū)路網(wǎng)逐漸完善，道路通達(dá)性提高。一般來說，道路通達(dá)性高，其可達(dá)的生態(tài)空間相對(duì)更多，但滬寧沿線城市、皖中、蘇北及杭州灣附近城市，由于城市化的快速發(fā)展，其本身生態(tài)空間質(zhì)量偏低，即使可達(dá)程度高，其基于可達(dá)性的生態(tài)空間質(zhì)量仍然偏低。而皖南和浙南地區(qū)道路通達(dá)性差，自然環(huán)境良好，生態(tài)空間質(zhì)量很高，其基于可達(dá)性的生態(tài)空間質(zhì)量較高。

圖3 基于可達(dá)性的生態(tài)空間質(zhì)量分布

3.3 長三角生態(tài)需求

按自然斷裂法將生態(tài)需求評(píng)價(jià)結(jié)果分為5個(gè)等級(jí)(圖4)。由于社會(huì)經(jīng)濟(jì)的發(fā)展，長三角人口和社會(huì)經(jīng)濟(jì)差異明顯，生態(tài)需求的時(shí)空差異顯著。在空間分布上，生態(tài)需求高值區(qū)集中在土地開發(fā)強(qiáng)度高、社會(huì)經(jīng)濟(jì)比較發(fā)達(dá)、人口密度較高的區(qū)域。其中，上海土地開發(fā)強(qiáng)度已接近飽和，皖西南和浙南丘陵地帶土地利用開發(fā)程度較低；上海人口密度最高，南京、蘇錫常等城市人口密度均處于前列，浙江和安徽除省會(huì)和臨海城市之外的人口密度較低；長三角社會(huì)經(jīng)濟(jì)總體比較發(fā)達(dá)，但各城市間的差別比較明顯，安徽各市GDP多處于下游水平。在時(shí)間演變上，1995—2015年，生態(tài)需求呈現(xiàn)擴(kuò)張趨勢(shì)，這與城市化發(fā)展和土地利用類型的轉(zhuǎn)變密切相關(guān)。考慮到實(shí)際居民的主觀因素，生態(tài)需求可能還與居民的收入和消費(fèi)水平、人口年齡層次、城市環(huán)境污染情況以及教育水平等因素相關(guān)。由于相關(guān)統(tǒng)計(jì)年鑒數(shù)據(jù)是以行政單元為主體進(jìn)行粗口徑的統(tǒng)計(jì)，本研究僅考慮人口、經(jīng)濟(jì)、土地利用等客觀因素對(duì)生態(tài)需求的影響。

圖4 生態(tài)需求分布

3.4 長三角生態(tài)空間供需平衡情況

長三角生態(tài)空間供需平衡空間分布具有明顯的時(shí)空差異(圖5)。從空間上看，供給不足的區(qū)域主要為上海、蘇州等經(jīng)濟(jì)發(fā)達(dá)城市，生態(tài)空間的需求很大，且同時(shí)缺少高質(zhì)量的生態(tài)空間，生態(tài)空間供給不足；供大于需的區(qū)域主要分布在浙南以及皖西南黃山、大別山、太湖等區(qū)域，這些區(qū)域自然條件好，生態(tài)空間質(zhì)量高，生態(tài)空間供大于需。從時(shí)間上看，1995—2015年，城市周邊生態(tài)空間供給不足，且呈擴(kuò)張趨勢(shì)。由于交通的發(fā)展使得人們出行更加便利，可以去到相對(duì)更遠(yuǎn)的地方，因此，生態(tài)空間供需平衡的區(qū)域有一定的增加，緩解了一些供給不足區(qū)域的生態(tài)環(huán)境壓力，甚至道路沿線的區(qū)域從供需平衡轉(zhuǎn)變?yōu)楣┐笥谛琛?/p>

圖5 生態(tài)空間供需平衡指數(shù)

利用空間統(tǒng)計(jì)分析軟件Geoda計(jì)算得到1995年和2015年的生態(tài)空間供需平衡指數(shù)與部分生態(tài)空間質(zhì)量評(píng)價(jià)因子以及生態(tài)需求因子之間的雙變量局部自相關(guān)LISA圖(圖6，圖7)，分析生態(tài)空間供需平衡指數(shù)與自然因素、社會(huì)因素之間的關(guān)系。浙南區(qū)域生態(tài)空間供需平衡指數(shù)與生態(tài)空間質(zhì)量因子呈“高-高”集聚，但同時(shí)生態(tài)空間供需平衡指數(shù)與生態(tài)需求因子呈“高-低”集聚，說明該區(qū)域資源條件良好，人口經(jīng)濟(jì)發(fā)展壓力較小，生態(tài)空間供大于需；蘇北和皖北以及省會(huì)和發(fā)達(dá)城市區(qū)域生態(tài)空間供需平衡指數(shù)與生態(tài)空間質(zhì)量因子呈“低-低”集聚，生態(tài)空間供需平衡指數(shù)與生態(tài)需求因子呈“低-高”集聚，說明該區(qū)域資源條件一般，且人口密集，生態(tài)空間供給不足。此外，這些因子之間還存在一定的趨同性。

圖6 1995年生態(tài)空間供需平衡指數(shù)與評(píng)價(jià)因子LISA圖

圖7 2015年生態(tài)空間供需平衡指數(shù)與評(píng)價(jià)因子LISA圖

4 結(jié)論與建議

4.1 結(jié)論

(1)基于融合可達(dá)性分析、統(tǒng)計(jì)分析方法建立的生態(tài)空間供需平衡指數(shù)，能夠很好地識(shí)別生態(tài)空間供需“短板”，可為區(qū)域生態(tài)空間的合理管控提供支持以提高生態(tài)空間的服務(wù)效率和價(jià)值。

(2)長三角生態(tài)空間質(zhì)量與土地利用類型密切相關(guān)，總體呈現(xiàn)由北向南逐漸升高的趨勢(shì)。生態(tài)質(zhì)量高值區(qū)分布在浙南、皖西南丘陵等受城市土地?cái)U(kuò)張建設(shè)及人類活動(dòng)干擾影響較小的區(qū)域，低值區(qū)分布在上海、蘇南、皖中、浙北城鎮(zhèn)及周邊。

(3)1995—2015年生態(tài)空間質(zhì)量較差的區(qū)域從各大城市中心部分?jǐn)U大到城市周邊地區(qū)，并逐漸成連片狀。各大城市城鎮(zhèn)周邊地區(qū)生態(tài)需求明顯增加，長三角生態(tài)空間供需平衡情況空間分布具有明顯的時(shí)空差異，供給不足的區(qū)域主要是經(jīng)濟(jì)和城市化發(fā)展較快的城市，供大于需的區(qū)域主要分布在生態(tài)高質(zhì)量區(qū)域，城市周邊生態(tài)空間供給不足，且呈擴(kuò)張趨勢(shì)。

4.2 建議

根據(jù)居民的休閑需求，精細(xì)化管理生態(tài)服務(wù)，構(gòu)建優(yōu)質(zhì)的城市內(nèi)部生態(tài)綠色空間，促進(jìn)宜居城市建設(shè)，深化綠色生活理念。因地制宜發(fā)展適合本地的生態(tài)服務(wù)功能，重點(diǎn)保護(hù)生態(tài)敏感區(qū)公益林，提升河湖水系包括太湖、長江、運(yùn)河等的水生態(tài)環(huán)境質(zhì)量，推廣綠色生態(tài)農(nóng)業(yè)。