摘要：水文過程改變、物理形態(tài)收縮是制約長江中游天鵝洲與上車灣故道作為國家級保護物種重要遷地保護區(qū)高質(zhì)量發(fā)展的主要因子。通過遙感解譯提取兩條故道1990～2021年逐年枯水期水域形態(tài)，識別主要形狀參數(shù)，分析計算岸線發(fā)育系數(shù)（SDI）、形態(tài)發(fā)育系數(shù)（Φ）及收縮特征系數(shù)（Rs）等，然后利用LASSO回歸模型構(gòu)建消除水域面積影響的形態(tài)參數(shù)序列，再利用PELT算法和分段擬合方法分析形態(tài)參數(shù)年際間變化規(guī)律。分析結(jié)果表明：閘控通江故道可用形態(tài)發(fā)育系數(shù)表征，自然通江故道可用收縮系數(shù)表征；閘控后天鵝洲故道牛軛湖化進程基本停滯，岸線逐漸趨向于簡單及平滑；隨著下荊江河道持續(xù)沖刷，上車灣故道牛軛湖化進程可能會逐步加速；建議采用水道疏浚、岸線整治、涵閘節(jié)制與調(diào)度等措施對故道形態(tài)進行保護。研究成果對于辨析三峽水庫運行、閘控等對故道物理形態(tài)的影響具有指導(dǎo)意義。

關(guān) 鍵 詞：天鵝洲故道； 上車灣故道； 遙感解譯； 物理形態(tài)； LASSO模型； PELT算法

中圖法分類號： TV147

文獻標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2024.11.004

0 引 言

故道是彎曲河流發(fā)生自然裁彎后的遺留河道，原河道的進口或出口在自然演變與人類活動的共同影響下發(fā)生淤塞，經(jīng)歷若干年后變成封閉或半封閉的淺水湖泊。近年來，通過沉積學(xué)和年代學(xué)方法研究故道內(nèi)泥沙沉積物，以反演沖積河流演變過程、裁彎時間及歷史洪水過程已經(jīng)形成較為成熟的研究體系^［1］。故道獨特的水文地質(zhì)條件，使其成為了一類典型的濱河濕地生物棲息地^［2］，故道與原河道的水文聯(lián)系改變引起的水環(huán)境及水生生物系統(tǒng)變化已受到相關(guān)學(xué)者的關(guān)注^［3-5］。盡管河道和洪泛湖泊（故道）之間的依存關(guān)系主要是由季節(jié)性洪水泛濫和水位變化所決定的^［6］，但是故道的物理形態(tài)也是決定其物質(zhì)交換特征的關(guān)鍵因素^［3］。河流故道物理形態(tài)變化是河湖水文連通變化與人類活動影響的重要表征，也是反映故道生境特征的關(guān)鍵指標(biāo)^［7-8］。

長江荊江河道獨特的水文地質(zhì)條件形成了中國規(guī)模最大的河流故道分布區(qū)，其中以位于下荊江的天鵝洲與上車灣故道最為典型^［9］。天鵝洲故道通過天鵝洲閘與長江干流連通，是典型的閘控通江故道；上車灣故道與長江干流自然連通，是典型的自然通江故道。兩條故道是眾多物種的良好棲息場所，依托該故道建有2個國家級自然保護區(qū)和1個省級自然保護區(qū)，是長江江豚、中華鱘、麋鹿等國家級保護物種的主要遷地保護區(qū)，也是長江生物多樣性保護的一個天然基因庫，具有重要的生態(tài)保護價值。故道物理形態(tài)的變化主要受相鄰荊江河道水沙情勢變化、河道演變及涵閘、堤壩等人為控制工程的影響。此外，圍墾、水產(chǎn)養(yǎng)殖和圩垸建設(shè)等改變了故道的物理形態(tài)，對故道容積及其與長江水文連通性產(chǎn)生不利影響，可能威脅故道社會與自然服務(wù)功能^［10］。因此，本文以典型閘控及自然通江故道為研究對象，通過遙感解譯獲取近30 a枯水期故道水域形態(tài)，利用形態(tài)識別方法分析各形態(tài)參數(shù)變化，并基于統(tǒng)計學(xué)模型分析形態(tài)參數(shù)年際間變化特征及其主要影響因素，以辨析三峽水庫運行、閘控等對故道物理形態(tài)的影響，進而支撐不同類型故道生境演變特征分析并提出保護修復(fù)措施。

1 研究區(qū)域基本情況

天鵝洲與上車灣故道均位于長江干流下荊江河段，天鵝洲故道位于石首市，上車灣故道位于監(jiān)利市，上車灣故道位于天鵝洲故道下游約60 km（圖1）。

天鵝洲故道為長江中游典型閘控通江牛軛湖，1972年因河道發(fā)生自然裁彎而形成，其上口已經(jīng)完全淤塞，下口通過長約6.4 km的河道與長江干流連通，由天鵝洲閘（1999年建成）節(jié)制，與長江干流連通口門位于石首水文站與監(jiān)利水文站之間，上游距石首水文站約22 km，下游距監(jiān)利水文站約42 km。天鵝洲保護區(qū)是國際認(rèn)可的長江江豚唯一成功遷地保護區(qū)范例，依托該故道建有石首麋鹿國家級自然保護區(qū)和長江天鵝洲白鱀豚國家級自然保護區(qū)。

上車灣為長江中游典型自然通江牛軛湖，1971年因?qū)嵤┤斯げ脧澒こ绦纬?，其上口已?jīng)完全淤塞，下口通過長約4 km的連通水道與長江干流自然連通。其生境條件與長江干流有諸多相似性，是眾多物種良好的棲息場所，依托該牛軛湖建有長江江豚及中華鱘遷移保護區(qū)，是長江生物多樣性保護的一個天然基因庫，具有重要的現(xiàn)實及潛在生態(tài)保護價值。

2 數(shù)據(jù)來源與分析方法

本文遙感影像及解譯手段均基于Google Earth Engine（GEE），遙感影像數(shù)據(jù)源為Landsat 5、Landsat 7 ETM 和 Landsat 8。在研究時段（1990年1月至2021年12月）內(nèi)，應(yīng)用時間過濾器來選擇1 a內(nèi)所有可用的Landsat表面反射率圖像，并將有效的Landsat圖像整理到年度圖像集合中。將基于像素質(zhì)量評估的CFmask算法應(yīng)用于集合中的每個圖像，以掩蓋來自云和云陰影像素的障礙物，利用修正后的歸一化差異水體指數(shù)MNDWI對故道水域范圍進行識別。根據(jù)Boothroyd等^［11］的研究成果，本文故道水域范圍識別的MNDWI閾值可被定義為淹水區(qū)域（gt;-0.05）和非淹水區(qū)域（lt;-0.05）。在獲取故道長序列水域范圍后，利用ArcGIS中Spatial Analyst等工具獲取不同年份故道水邊線、面積、周長、中心軸線、連通水道長度等。

本文采用特征指標(biāo)法分析長江中游典型故道近30 a物理形態(tài)變化特征。故道形態(tài)可表征彎曲河流曾經(jīng)的演變痕跡，以指示彎曲河流演變過程。Weihaupt^［12］通過實地查勘全面統(tǒng)計了位于北美洲育空河（Yukon River）的河流故道，共計817條，并根據(jù)故道形態(tài)特征將其分為簡單型、復(fù)合型和復(fù)雜型，根據(jù)封閉程度將其分為開放型、正常型和封閉型。相對應(yīng)地，李志威等^［1］將河流故道分為“Ω”形、“U”形及月牙形，以呈現(xiàn)故道在全生命期內(nèi)不同階段的形態(tài)特征。Wang等^［13］提出使用收縮率和曲率來表征故道物理形態(tài)，以確定故道與河道之間的連通性，以及故道水質(zhì)之間的相關(guān)性。

根據(jù)Weihaupt^［12］和李志威^［1］等提出的分類方法，本文所研究的故道均呈“Ω”形，處于故道演化的初期，兼顧湖泊與河道的特點。基于圖形周界測度算法的彎曲指數(shù)可以很好地描述故道的形態(tài)變化特征，本節(jié)引入形態(tài)發(fā)育系數(shù)（Φ）來表征故道的彎曲程度變化特征：

Φ=L/4A/π

（1）

式中：L為故道的中軸線長度，m，見圖2；

A為故道水域面積，m2。Φ（Φgt;1）可以用于表征牛軛湖發(fā)育情況，Φ越小表明故道湖泊化程度越高，Φ越大表明故道河道化程度越高。Φ由低到高，表明故道由“Ω”形向“U”形轉(zhuǎn)變。

為表征故道周邊生產(chǎn)生活對其物理形態(tài)的影響，分析幾何形態(tài)及衡量場地空間數(shù)據(jù)圖形的完整性和聚集性，本節(jié)引入岸線發(fā)育系數(shù)（SDI）以表征湖泊幾何形態(tài)及衡量場地空間數(shù)據(jù)圖形的完整性和聚集性^［14］：

SDI=P/2πA

（2）

式中：P為故道岸線長度，m。SDI越大表示故道岸線越曲折，由于故道存在河道屬性，其岸線越曲折表明受人類活動影響程度越大，自然岸線保有率越低。

在故道演化過程中，發(fā)生變化最早的一般為口門和與河道連通的水道，收縮指數(shù)可以很好地描述這一過程。這里引入收縮系數(shù)（Rs）來表征故道湖泊化程度：

Rs=W/（W+Ld+Lu）

（3）

式中：W為故道平均寬度，m；

Ld為下通江水道長度，m；

Lu為上通江水道長度，m。

Rs取值區(qū)間為（0，1），趨向于0表示故道湖泊化程度越高，與河道的聯(lián)系更微弱，趨向于1表示故道河道化程度越高，受河道的影響越劇烈。

枯水期故道形態(tài)變化特征可以較好地反映故道演變過程。由于衛(wèi)星影像質(zhì)量和數(shù)量的限制，很難尋找到同一水位或水域面積條件下的影像。而本文所選取的衛(wèi)星影像均是在枯水期（1～3月）提取的，每個年度形狀系數(shù)對應(yīng)的水域面積并不一致，因此，本文采用LASSO回歸模型^［15］對時間序列數(shù)據(jù)進行處理，以獲取消除水域面積變化影響的形狀系數(shù)，進而分析形狀系數(shù)的年際間變化規(guī)律。在此基礎(chǔ)上，采用一種高效、準(zhǔn)確的序列數(shù)據(jù)斷點檢測算法——PELT（pruned exact linear time）算法^［16］分析形狀系數(shù)時間序列的變化趨勢，提取變化趨勢的分段點，利用分段回歸模型分析形狀系數(shù)時間序列的變化規(guī)律。

3 結(jié)果分析

基于遙感解譯及形狀參數(shù)提取獲取了各年度形狀系數(shù)，采用LASSO回歸模型，以年份及水域面積變化為解釋變量，以3類形狀系數(shù)為響應(yīng)變量，通過正則化參數(shù)Lambda的自適應(yīng)調(diào)整，解析形狀系數(shù)年際間變化趨勢。圖3、4分別給出了天鵝洲和上車灣故道形狀系數(shù)預(yù)測值與實測值的對比情況，預(yù)測值為消除了不同季節(jié)水域面積不一致的影響，可以真實反映故道的形態(tài)特征。

基于PELT算法和分段回歸模型識別了各形狀參數(shù)年際間變化趨勢。圖5給出了天鵝洲及上車灣故道形態(tài)發(fā)育系數(shù)（Φ）年際間變化趨勢?？梢钥闯鎏禊Z洲故道Φ值可分為3個區(qū)間：1990～1999年間呈下降趨勢，下降率為-0.010 12；1999～2003年呈上升趨勢，上升率為0.014 13；2004～2021年無顯著變化。上車灣故道Φ值可分為兩個區(qū)間：1990～2012年呈顯著上升趨勢，上升率為0.020 83；2013～2021年無顯著變化。上車灣故道Φ值總體上顯著高于天鵝洲故道，表明上車灣故道河道化程度要顯著高于天鵝洲故道；同時，上車灣故道河道化程度在2012年前顯著增加，該區(qū)間內(nèi)故道由“Ω”形向“U”形轉(zhuǎn)變。天鵝洲故道Φ值總體穩(wěn)定，在天鵝洲閘建成前，呈小幅度下降趨勢，至三峽水庫運行，呈小幅上升趨勢，此后基本穩(wěn)定。形態(tài)發(fā)育系數(shù)主要受故道與長江干流水體交換強度的影響，長江干流水位是決定Φ值的關(guān)鍵因素，并以中高水位為主。上車灣故道與長江干流連通口門位于監(jiān)利站下游約20 km，受地形調(diào)整、綜合糙率及洞庭湖頂托的影響，上車灣故道口門所在干流中流量條件下水位變化不顯著^［17］，高流量條件下水位略有抬升^［18］，上車灣故道與長江干流水力聯(lián)系呈增長趨勢。

圖6給出了天鵝洲及上車灣故道SDI值年際間變化趨勢。天鵝洲故道SDI指數(shù)介于3.0～3.4之間，上車灣故道SDI指數(shù)介于3.3～3.6之間，均高于一般天然湖泊^［19］。上車灣故道SDI值變化趨勢與Φ值相似，可分為兩個區(qū)間：1990～2012年呈顯著上升趨勢，上升率為0.011 80；2013～2021年無顯著變化。天鵝洲故道SDI值變化趨勢可分為3個區(qū)間：1990～1999年間呈上升趨勢，上升率為0.030 85；1999～2003年呈下降趨勢，下降率為-0.043 08；2004～2021年呈小幅下降趨勢，下降率為-0.008 86。1999年以前天鵝洲故道水位周期性漲落，岸線水力侵蝕強度較大，枯水期水域面積小，岸線不規(guī)則程度逐步加劇；天鵝洲閘建成后，故道年內(nèi)水位變化趨于平坦，枯水期水位大幅度抬升，加之對故道水域空間的管護力度逐步增加，故道岸線趨于簡單及平滑。上車灣故道SDI值總體上顯著高于天鵝洲故道，表明上車灣故道岸線受人類活動的影響程度要顯著高于天鵝洲故道，人類活動主要為人工魚塘、灘涂圍墾等。

圖7給出了天鵝洲及上車灣故道Rs值年際間變化趨勢，Rs值變化趨勢與Φ值基本相反。天鵝洲故道Rs值可分為3個區(qū)間：1990～1999年間呈上升趨勢，上升率為0.002 84；1999～2003年呈下降趨勢，下降率為-0.003 96；2004～2021年無顯著變化。上車灣故道Φ值可分為兩個區(qū)間：1990～2012年呈顯著下降趨勢，下降率為-0.007 09；2013～2021年無顯著變化。上車灣故道Rs值顯著高于天鵝洲故道，表明上車灣故道與長江干流水文聯(lián)系更為緊密，河道化程度更高，這一規(guī)律與Φ值是一致的。天鵝洲故道與長江水體交換強度在建閘后大幅降低，閘門在主汛期一般不開啟，汛期前后根據(jù)干流水位情況進行引水，枯水期用于蓄水，一般不開啟。因此，長江干流枯水同流量下水位下降對天鵝洲故道基本無影響，中高流量下水位變化的影響也不顯著。天鵝洲閘建成前，故道河道化程度呈上升趨勢，之后呈下降趨勢，故道趨向于封閉湖泊，三峽水庫運用后基本趨于穩(wěn)定。

對于閘控連通的天鵝洲故道，Rs與SDI指數(shù)的變化趨勢相同，與Φ的變化趨勢相逆。對于自然連通的上車灣故道，Φ與SDI的變化趨勢相同，與Rs的變化趨勢相反。Φ與Rs均可表征故道的牛軛湖化進程，但是側(cè)重點有所不同，Φ主要關(guān)注縱向發(fā)育，Rs側(cè)重于橫向發(fā)育與連通水道變化。對于閘控連通故道，與長江干流水流交換產(chǎn)生的縱向沖刷作用基本消失，通過風(fēng)生流、水位變化等產(chǎn)生的橫向侵蝕作用發(fā)揮主要作用，因此，Rs與SDI的變化趨勢相同。對于自然連通故道，與長江干流水流交換產(chǎn)生的縱向沖刷作用占主導(dǎo)，因此，Rs與SDI的變化趨勢相反。因而，在評估江河故道物理形態(tài)演變時，對于自然通江的故道，Rs更適合作為指示參數(shù)；對于閘控通江故道，Φ更適合作為指示參數(shù)。

4 結(jié)論及建議

本文利用遙感解譯、形態(tài)參數(shù)與統(tǒng)計學(xué)模型，對1990～2021年閘控連通的天鵝洲故道與自然連通的上車灣故道的岸線發(fā)育、形態(tài)收縮與形態(tài)發(fā)育等物理形態(tài)演變特征進行了分析，主要結(jié)論如下：

（1） 對于閘控連通的天鵝洲故道，收縮系數(shù)（Rs）與岸線發(fā)育系數(shù)（SDI）的變化趨勢相同，與形態(tài)發(fā)育系數(shù)（Φ）的變化趨勢相反。對于自然連通的上車灣故道，Φ與SDI的變化趨勢相同，與Rs的變化趨勢相反。在評估江河故道物理形態(tài)演變時，對于自然通江的故道，Rs更適合作為指示參數(shù)；對于閘控通江故道，Φ適合作為指示參數(shù)。

（2） 近30 a天鵝洲故道形態(tài)參數(shù)變化可分為3個區(qū)間，節(jié)點分別為天鵝洲閘建成、三峽水庫運用。Rs與SDI在故道被天鵝洲閘節(jié)制前呈上升趨勢，至三峽水庫運用前呈下降趨勢，三峽水庫運用后基本穩(wěn)定。Φ在故道被天鵝洲閘節(jié)制前呈下降趨勢，至三峽水庫運用前呈上升趨勢，三峽水庫運用后基本穩(wěn)定。

上車灣故道形態(tài)參數(shù)變化可分為兩個區(qū)間，節(jié)點為2012年、三峽水庫滿蓄運用。Φ與SDI在2012年前呈顯著上升趨勢，之后則基本保持穩(wěn)定。Rs在2012年前呈顯著下降趨勢，之后則基本保持穩(wěn)定。

（3） 在閘控后，天鵝洲故道牛軛湖化進程基本停滯，岸線逐漸趨于簡單與平滑。隨著下荊江河道持續(xù)沖刷，中高流量條件下水位降低，上車灣故道牛軛湖化進程可能會逐步加速，岸線逐步趨于簡單與平滑。

為了進一步保護天鵝洲與上車灣故道，維持其作為基因庫及重要物種保護地的功能，建議開展以下工作：① 研究物理形態(tài)變化對故道水環(huán)境變化及重要物種生境的影響，提出關(guān)鍵閾值區(qū)間，以指導(dǎo)后續(xù)保護工作；② 進一步優(yōu)化閘控措施，研究提出故道形態(tài)調(diào)節(jié)、水質(zhì)改善等措施；③ 針對干流水位持續(xù)下降、故道與長江水沙聯(lián)系減弱的問題，可以適當(dāng)考慮疏浚、涵閘改造等工程措施，以加強水沙聯(lián)系，恢復(fù)故道水文節(jié)律。

參考文獻：

［1］ 李志威，王兆印，潘保柱.牛軛湖形成機理與長期演變規(guī)律［J］.泥沙研究，2012（5）：16-25.

［2］ STELLA J C，HAYDEN M K，BATTLES J J，et al.The role of abandoned channels as refugia for sustaining pioneer riparian forest ecosystems［J］.Ecosystems，2011，14：776-790.

［3］ OBOLEWSKI K，GLIńSKA-LEWCZUK K J.Effects of oxbow reconnection based on the distribution and structure of benthic macroinvertebrates［J］.Clean-Soil Water，2011，39（9）：853-862.

［4］ 劉心愿，渠庚，姚仕明，等.三峽工程運用后石首彎道段整治工程累積影響和演變趨勢研究［J］.水利水電快報，2020，41（1）：22-27.

［5］ GUO X，GAO P，LI Z.Hydrologic connectivity and morphologic variation of oxbow lakes in an alpine pristine fluvial system［J］.Journal of Hydrology，2023，623：129768.

［6］ ASPETSBERGER F，HUBER F，KARGL S，et al.Particulate organic matter dynamics in a river floodplain system：impact of hydrological connectivity［J］.Archiv Fur Hydrobiologie，2002，156（1）：23-42.

［7］ ANDERSEN T K，BOLDING K，NIELSEN A，et al.How morphology shapes the parameter sensitivity of lake ecosystem models［J］.Environmental Modelling amp; Software，2021（2）：136-159.

［8］ CHUKWUKA A V，EGWARE T U，OKALI K D，et al.The influence of lake morphology，landscape structure，and urbanization factors on bird community composition in wetlands of four tropical lakes［J］.Wetlands，2022，42（7）：1-15.

［9］ 蔡曉斌，燕然然，王學(xué)雷.下荊江故道通江特性及其演變趨勢分析［J］.長江流域資源與環(huán)境，2013，22（1）：53-58.

［10］劉小光，姚仕明，柴朝暉.湖泊與河流水文連通量化計算模型及應(yīng)用［J］.湖泊科學(xué)，2024，36（5）：1516-1524.

［11］BOOTHROYD R J，WILLIAMS R D，HOEY T B，et al.Applications of Google Earth Engine in fluvial geomorphology for detecting river channel change［J］.Wiley Interdisciplinary Reviews：Water，2021，8（1）：e21496.

［12］WEIHAUPT J G.Morphometric definitions and classifications of oxbow lakes，Yukon River Basin，Alaska［J］.Water Resources Research，1977，13（1）：195-196.

［13］WANG D，LI Z，LI Z，et al.Environmental gradient relative to oxbow lake-meandering river connectivity in Zoige Basin of the Tibetan Plateau［J］.Ecological Engineering，2020，156：105983.

［14］王哲，劉凱，詹鵬飛，等.近三十年青藏高原內(nèi)流區(qū)湖泊岸線形態(tài)的時空演變［J］.地理研究，2022，41（4）：980-996.

［15］韓曉育，郭潁奎.耦合 LASSO 回歸的 HHO-LSVR 中長期徑流預(yù)報模型［J］.水文，2021，41（3）：69-74.

［16］KILLICK R，F(xiàn)EARNHEAD P，ECKLEY I A.Optimal detection of changepoints with a linear computational cost［J］.Journal of the American Statistical Association，2012，107（500）：1590-1598.

［17］尹祖迎，柴元方.長江中游宜昌—螺山段同流量下中水位變化特征及機制研究［J］.水電能源科學(xué)，2022，40（10）：31-34，86.

［18］柴元方，鄧金運，楊云平，等.長江中游荊江河段同流量—水位演化特征及驅(qū)動成因［J］.地理學(xué)報，2021，76（1）：101-113.

［19］張鳳太，王臘春，冷輝，等.典型天然與人工湖泊形態(tài)特征比較分析［J］.中國農(nóng)村水利水電，2012（7）：38-41.

（編輯：胡旭東）

Morphological evolution characteristics of typical oxbow connecting middle reaches of Yangtze River

LIU Xiaoguang¹，SUN Guizhou¹，YAO Rui²，JIN Zhongwu¹，CHAI Zhaohui¹

（1.Key Laboratory of River and Lake Regulation and Flood Control in the Middle and Lower Reaches of the Yangtze River of Ministry of Water Resources，Changjiang River Scientific Research Institute，Wuhan 430010，China； 2.Hefei University of Technology Gongda Engineering Testing and Experiment Co.，Ltd.，Hefei 230041，China）

Abstract：

Hydrological connectivity barriers and physical morphology shrinkage are the main factors restricting Tian′ezhou and Shangchewan oxbow′s role as important translocation protection areas for nationally protected species.Through remote sensing interpretation，the water morphology during dry seasons from 1990 to 2021 was extracted annually，and the main shape parameters were identified.Analysis of physical morphology parameters such as shoreline development coefficient （SDI），morphology development coefficient（Φ），and shrinkage coefficient （R_s） were conducted.A morphological parameter sequence eliminating the influence of water area was constructed using the LASSO regression model.The annual variation in morphological parameters were analyzed using the PELT algorithm and segmented fitting method.The gate regulated communication channel can be characterized by the morphological development coefficient，while the natural communication channel can be characterized by the shrinkage coefficient.The process of transformation of Tian′ezhou oxbow lake has basically stagnated after regulation，and the shoreline has gradually tended to become simpler and smoother.With the continuous scouring of the Shangchewan river channel，the process of oxbow lake transformation may gradually accelerate.It is recommended to adopt measures such as waterway dredging，shoreline improvement，sluice gate controlling and regulation to protect the morphology of the oxbow.The research results have guiding significance for the analysis on influence of the Three Gorges Reservoir operation and gate regulation on the physical form of oxbows.

Key words：

Tian′ezhou oxbow lake； Shangchewan oxbow lake； remote sensing interpretation； physical morphology； LASSO model； PELT algorithm