吳良杰 薄景山 彭達(dá)

摘要：

黃土地震古/老滑坡原始地形恢復(fù)是開展黃土地震古/老滑坡研究的基礎(chǔ)。以寧夏回族自治區(qū)固原市西吉縣平峰鄉(xiāng)三合村潘溝組1號(hào)滑坡為例，利用衛(wèi)星遙感解譯和野外調(diào)查等技術(shù)和方法采集滑坡現(xiàn)今地形資料數(shù)據(jù);并基于滑坡滑動(dòng)前后滑坡體體積相等的原則，以MATLAB軟件為平臺(tái)，利用三維曲面插值技術(shù)恢復(fù)滑坡的原始地形，恢復(fù)后的滑坡原始地形斜坡坡面形態(tài)整體呈外凸型，坡度分布范圍在12°～15°間；最后，利用黃土地貌單元的坡度分布范圍和滑坡穩(wěn)定性計(jì)算論證恢復(fù)結(jié)果的合理性和可行性。文章研究成果可為缺乏滑動(dòng)前地形資料的黃土地震古/老滑坡原始地形的恢復(fù)提供新方法，對(duì)黃土地震古/老滑坡的研究具有重要的理論意義和工程應(yīng)用價(jià)值。

關(guān)鍵詞：

黃土古老滑坡; 地形恢復(fù); 三維曲面插值; 地形數(shù)據(jù); 穩(wěn)定性計(jì)算

中圖分類號(hào)： P694;P642????? 文獻(xiàn)標(biāo)志碼：A?? 文章編號(hào)： 1000-0844（2023）02-0329-09

DOI：10.20000/j.1000-0844.20211006002

Restoration method of the original topography of

ancient seismic landslides in loess areas

WU Liangjie1， BO Jingshan2，3， PENG Da3

（1. Geophysical Exploration Center， CEA， Zhengzhou? 450011， Henan， China;

2. Key Laboratory of Earthquake Engineering and Engineering Vibration， Institute of Engineering Mechanics，

CEA，Harbin 150080， Heilongjiang， China;

3. Institute of Disaster Prevention， Sanhe 065201， Hebei， China）

Abstract：

Restoration of the original topography is the basis for the study of ancient seismic landslides in loess areas. The No.1 landslide in the Pangou Group， Sanhe Village， Pingfeng Township （Xiji County， Ningxia Hui Autonomous Region） was taken as an example in this paper. Satellite remote sensing interpretation and field survey were adopted to collect the current topographic data of the landslide. Based on the principle of equal volumes of landslide bodies before and after sliding， MATLAB was used to restore the original topography of the No. 1 landslide using 3D surface interpolation technology. The surface shape of the restored landslide's original topography is generally convex， with a slope distribution range of 12°—15°. Finally， the rationality and feasibility of the restoration results were demonstrated by the slope distribution range of the loess geomorphic units and the landslide stability calculation. Results provide a novel approach for the restoration of the original topography of ancient seismic loess landslides that lacks the topography data before sliding， thus providing important theoretical significance and engineering application value for related studies.

Keywords：

ancient loess landslide; topography restoration; 3D surface interpolation; topography data; stability calculation

0 引言

黃土地震滑坡是黃土地區(qū)典型的地震地質(zhì)災(zāi)害［1］。地震災(zāi)害資料顯示，黃土地區(qū)歷次中強(qiáng)地震都誘發(fā)了大量的黃土滑坡，如1654年天水8.0級(jí)特大地震誘發(fā)了180多個(gè)規(guī)模不等的黃土地震滑坡，滑坡分布總面積達(dá)1 000 km2［2-3］;1718年通渭7.5級(jí)大地震共誘發(fā)了300多起大規(guī)模的黃土地震滑坡，分布總面積達(dá)3 500 km2［4］;1920年海原8.5級(jí)特大地震約誘發(fā)了600多起大規(guī)模的黃土地震滑坡，分布面積達(dá)4 000～5 000 km2［5-7］。因此，開展黃土地震滑坡研究對(duì)減輕黃土地區(qū)的地震災(zāi)害具有重要的理論意義和工程應(yīng)用價(jià)值。本文所稱古/老滑坡是工程地質(zhì)學(xué)中古滑坡和老滑坡的統(tǒng)稱［8］。由于形成在黃土地區(qū)的古/老滑坡數(shù)量多，種類齊全，且形態(tài)保留較完整，因此可通過(guò)黃土古/老滑坡的研究來(lái)探索黃土地震滑坡的形成機(jī)理和形成條件。開展黃土地震滑坡研究需要掌握滑坡滑動(dòng)前的地形資料，但是黃土古/老滑坡由于形成年代久遠(yuǎn)，無(wú)法獲取相關(guān)地形資料，因此開展黃土地震古/老滑坡原始地形恢復(fù)方法的研究已成為當(dāng)前滑坡研究領(lǐng)域一個(gè)亟待解決的問(wèn)題。古/老滑坡地形恢復(fù)研究工作始于21世紀(jì)初，李忠生［9］、李萍［10］、Conoscenti C等［11］、鄭旭輝［12］、周鑫等［13］針對(duì)不同的研究問(wèn)題提出了一些古/老滑坡地形恢復(fù)方法，取得了一定的研究成果，但現(xiàn)有的古/老滑坡原始地形恢復(fù)方法精度不高，遠(yuǎn)不能滿足實(shí)際工作的需要。

本文基于滑坡體滑動(dòng)前后體積相等的原則，結(jié)合MATLAB中三維曲面插值程序Gridfit［14］和三次樣條插值程序Spline［15］，提出一種黃土地震古/老滑坡原始地形的恢復(fù)方法。該方法的思路為：在采集滑坡現(xiàn)今地形數(shù)據(jù)的基礎(chǔ)上，估算黃土地震滑坡體的體積，并以剪出口最低點(diǎn)與滑坡后緣最高點(diǎn)連線的傾角作為滑坡滑動(dòng)前斜坡的原始坡度，通過(guò)三維曲面插值技術(shù)對(duì)滑坡周界的地形數(shù)據(jù)和滑坡剪出口與滑坡后緣連線上的點(diǎn)進(jìn)行插值和擬合，從而近似得到滑坡滑動(dòng)前斜坡的地表形狀；在此基礎(chǔ)上，利用滑坡的地形數(shù)據(jù)建立一個(gè)滑坡滑動(dòng)前的初始斜坡地形模型，在此模型的基礎(chǔ)上估算出滑坡體滑動(dòng)前的體積；最后，基于滑坡體滑動(dòng)前與滑動(dòng)后體積相等的原則，通過(guò)調(diào)整模型坡面的凹凸程度得到滑坡滑動(dòng)前的最終地形模型。并利用該方法對(duì)寧夏固原市西吉縣平峰鄉(xiāng)三合村潘溝組1號(hào)滑坡開展原始地形恢復(fù)工作，論證恢復(fù)結(jié)果的合理性。

1 滑坡概況

三合村1號(hào)滑坡位于寧夏回族自治區(qū)固原市西吉縣平峰鄉(xiāng)三合村潘溝組，系1920年海原8.5級(jí)特大地震引發(fā)的滑坡。其位于地震烈度IX度區(qū)內(nèi)，地理坐標(biāo)為105.49° E，35.79° N。由于西北黃土地區(qū)氣候干旱，降雨量少，滑坡的地貌形態(tài)保留較完整。滑坡現(xiàn)今地貌和等高線分別如圖1和圖2所示。我們通過(guò)野外調(diào)查和衛(wèi)星遙感解譯，獲得該滑坡的基本信息：滑坡主滑動(dòng)方向約為70°方向，主滑長(zhǎng)度約832 m，寬度約298 m，原始坡高為 154 m;周邊未滑動(dòng)區(qū)域的坡角為 16.3°，滑坡體坡角為 6°;滑坡后壁高差為58 m，后壁傾角為39°;滑坡左側(cè)壁高差為24.5 m，左側(cè)壁傾角為46.5°;滑坡右側(cè)壁高差為20 m，右側(cè)壁傾角為40°。通過(guò)現(xiàn)場(chǎng)估算，滑坡體最大滑距約為360 m。該滑坡在平面上大致呈矩形，地層結(jié)構(gòu)為泥巖-黃土二元結(jié)構(gòu)：下部為白堊紀(jì)的紅色泥巖，隔水性能較好，遇水強(qiáng)度降低;其上為第四紀(jì)馬蘭黃土，厚度約為60 m。在滑坡滑動(dòng)面上未見明顯的紅色泥巖出露，因此可判定其為黃土層內(nèi)滑坡。綜合衛(wèi)星遙感解譯和現(xiàn)場(chǎng)調(diào)查的結(jié)果，判定該滑坡形成于黃土梁地貌，其地質(zhì)剖面圖如圖3所示。

2 滑坡邊界確定和地形數(shù)據(jù)采集

2.1 滑坡周界、剪出口和原始周界的確定

滑坡周界是滑坡體與周邊未滑區(qū)域的分界線，它圈定了滑坡在平面上的范圍［16］。在衛(wèi)星影像中，通常取滑坡滑動(dòng)區(qū)域與周邊未滑區(qū)域之間地表輪廓不連續(xù)或銜接不和諧的分界線作為滑坡周界，如圖4所示?；录舫隹谑腔瑒?dòng)面最下端與原地面相交剪出的破裂帶，或滑動(dòng)面最下端與原始地面的交線［17］，通常作為滑坡滑動(dòng)前原始周界最下端的界限，但剪出口通常被掩埋于滑坡體之下，不易獲取。工程上獲取剪出口信息的方法有鉆探、物探等，但這些方法成本高、費(fèi)時(shí)費(fèi)力，對(duì)研究區(qū)域性滑坡不具備普適性。根據(jù)野外調(diào)查的經(jīng)驗(yàn)，滑坡剪出口的位置通常位于滑坡側(cè)壁最下端附近，因此，為了方便和快速地確定滑坡剪出口的位置，本研究根據(jù)野外調(diào)查的結(jié)果，取滑坡兩側(cè)壁邊界線最下端位置的連線作為滑坡剪出口的大致位置，并取剪出口的中點(diǎn)位置為其最低點(diǎn)，將滑坡前緣最低點(diǎn)的高程值作為剪出口最低點(diǎn)的高程值。同時(shí)，取滑坡后緣、滑坡側(cè)壁頂部邊界線和滑坡剪出口三者所圍成的閉合線作為滑坡滑動(dòng)前的原始周界，如圖5所示。

2.2 滑坡地形資料的采集

基于30 m精度的DEM高程數(shù)據(jù)，可根據(jù)不同地形要素的特點(diǎn)，按照不同的采樣間隔來(lái)采集滑坡不同要素的地形數(shù)據(jù)。利用91衛(wèi)圖助手軟件分別采集三合村1號(hào)滑坡的滑坡后壁、滑坡后緣、滑坡體、滑坡側(cè)壁頂部邊界、滑坡舌、滑坡剪出口和未滑區(qū)域不動(dòng)體等部位的地形數(shù)據(jù)，采樣點(diǎn)位置分布如圖6所示。本文采集地形數(shù)據(jù)坐標(biāo)時(shí)，采用國(guó)際通用的地理投影坐標(biāo)系——WGS84高斯投影坐標(biāo)系，地形數(shù)據(jù)單位統(tǒng)一采用m。

3 滑動(dòng)面確定和滑體體積的估算

滑動(dòng)面位置和滑坡體體積是滑坡原始地形恢復(fù)中兩個(gè)重要的幾何參數(shù)。在已采集的滑坡地形資料的基礎(chǔ)上，利用MATLAB軟件中的三次樣條插值程序進(jìn)行插值，可近似得到滑動(dòng)面的位置。在此基礎(chǔ)上，結(jié)合滑動(dòng)面與滑坡體表面之間的高差和滑坡體的表面積，近似估算出滑坡體的體積。

3.1 滑動(dòng)面的確定

結(jié)合野外調(diào)查經(jīng)驗(yàn)，提出一種利用三次樣條函數(shù)插值確定滑動(dòng)面的快捷方法。由于三次樣條函數(shù)插值在三維坐標(biāo)系中有一定的局限性，但在二維坐標(biāo)系中則比較成熟，因此在插值前需將三維坐標(biāo)系中的插值轉(zhuǎn)化為二維坐標(biāo)系的插值。在三維坐標(biāo)系的基礎(chǔ)上，以滑坡后緣最高點(diǎn)的水平投影點(diǎn)為原點(diǎn)O′，豎向高程為Z軸，滑坡后緣最高點(diǎn)與剪出口最低點(diǎn)的連線的水平向投影方向?yàn)镾軸，建立一個(gè)S-Z二維坐標(biāo)系進(jìn)行插值。具體步驟如下：

第一步，將滑坡后壁、滑坡體、剪出口最低點(diǎn)等位置的采樣點(diǎn)換算至S-Z坐標(biāo)系中;

第二步，確定一個(gè)調(diào)整系數(shù)，確定方法為取滑坡后壁高差、滑坡左側(cè)壁高差以及右側(cè)壁高差二分之一的平均值作為該系數(shù)值，取滑坡體中點(diǎn)位置采樣點(diǎn)的高程值減去調(diào)整系數(shù)所得到的點(diǎn)作為滑動(dòng)面的中點(diǎn)位置;

第三步，以滑坡后壁、滑動(dòng)面中點(diǎn)位置和剪出口最低點(diǎn)位置的采樣點(diǎn)為被插值點(diǎn)，利用三次樣條插值程序?qū)ζ溥M(jìn)行插值來(lái)近似確定滑動(dòng)面的位置，如圖7所示；

第四步，將已得到的滑動(dòng)面位置的插值點(diǎn)和滑坡地形數(shù)據(jù)換算至三維坐標(biāo)系中，利用MATLAB軟件的三維曲面插值程序Gridfit對(duì)滑動(dòng)面位置的插值點(diǎn)和滑坡其他地形采樣點(diǎn)進(jìn)行插值，得到滑動(dòng)面的近似三維地形數(shù)據(jù);

第五步，利用MATLAB中的meshz程序?qū)瑒?dòng)面的三維地形數(shù)據(jù)進(jìn)行擬合，得到滑動(dòng)面的三維地形圖（圖8）。

該方法目前僅適用于一般性的黃土層內(nèi)滑坡，不適用于其他非黃土層內(nèi)的滑坡。

3.2 滑坡滑動(dòng)后體積的估算

利用滑坡體表面積與滑坡體平均厚度的乘積來(lái)估算滑坡體的體積：

Va=Sa×Ha （1）

式中：Sa為滑動(dòng)后滑坡體的表面積，近似取自滑坡體表面積的水平投影面積，利用91衛(wèi)圖助手軟件的面積測(cè)量功能測(cè)得三合村1號(hào)滑坡的Sa值約為246 854.2 m2（圖9）;Ha為滑動(dòng)后滑坡體的平均厚度，取滑坡體現(xiàn)今地表與滑動(dòng)面以及滑坡舌以下原地面之間若干點(diǎn)高差的平均值（圖10），通過(guò)計(jì)算得到三合村1號(hào)滑坡的Ha值約為13 m;Va為滑動(dòng)后滑坡體的體積，估算得出三合村1號(hào)滑坡的Va值約為3 209 105.0 m3。

4 滑坡原始地形的恢復(fù)

野外調(diào)查和有關(guān)研究工作表明，自然形成的黃土斜坡坡面形態(tài)大致分為平直型、內(nèi)凹型、外凸型三種［18］。因此，在滑坡原始地形恢復(fù)中，可首先假定滑坡滑動(dòng)前的斜坡坡面形態(tài)為平直型，再利用MATLAB的三維曲面插值程序Gridfit和三次樣條插值程序Spline對(duì)所采集的滑坡地形數(shù)據(jù)進(jìn)行插值和擬合，近似得到其初始的三維地表形態(tài)。在此基礎(chǔ)上，基于滑坡體滑動(dòng)前后體積相等的原則，對(duì)斜坡的坡面形態(tài)進(jìn)行適當(dāng)調(diào)整，以此來(lái)確定滑坡滑動(dòng)前的三維地形。

4.1 滑動(dòng)前斜坡初始地形模型的建立

假定滑坡滑動(dòng)前斜坡的坡面形態(tài)為平直型，在S-Z坐標(biāo)系中繪制滑坡滑動(dòng)前的斜坡剖面圖（圖11）。取滑坡后緣最高點(diǎn)與剪出口最低點(diǎn)的連線為滑動(dòng)前斜坡的坡面，取該線的傾角為滑動(dòng)前斜坡的坡度α，通過(guò)計(jì)算得到三合村1號(hào)滑坡的α為14°。通過(guò)擬合得到該直線在Z-S坐標(biāo)中的數(shù)學(xué)表達(dá)式：

Z=-0.246 4S+1 963 （2）

利用式（2），采用等間隔的方式求得滑動(dòng)前斜坡坡面上若干點(diǎn)的高程值，將這些點(diǎn)換算至三維坐標(biāo)系中。在三維坐標(biāo)系中，利用MATLAB軟件中的三維曲面插值程序Gridfit對(duì)斜坡坡面上的若干點(diǎn)和滑坡其余地形采樣點(diǎn)進(jìn)行插值，再利用meshz程序?qū)λ玫降牟逯迭c(diǎn)數(shù)據(jù)進(jìn)行擬合，得到三合村1號(hào)滑坡滑動(dòng)前斜坡的三維地形圖（圖12）。

4.2 滑坡滑動(dòng)前滑體體積的估算

采用滑坡體表面積與平均厚度的乘積來(lái)估算滑動(dòng)前滑坡體的體積Vb值：

Vb=Sb×Hb （3）

式中：Sb為滑動(dòng)前滑坡體的表面積，主要根據(jù)滑坡原始周界所圈定的區(qū)域范圍內(nèi)滑坡體表面積的水平投影面積確定，圖13給出三合村1號(hào)滑坡Sb的測(cè)量結(jié)果約為115 002.9 m2;Hb為滑動(dòng)前滑坡體的平均厚度，取自滑動(dòng)前斜坡坡面與滑動(dòng)面之間若干點(diǎn)的高差的平均值，通過(guò)計(jì)算得到三合村1號(hào)滑坡的Hb值約為25 m。滑坡滑動(dòng)前斜坡坡面與滑動(dòng)面高程對(duì)比如圖14所示。經(jīng)過(guò)估算，三合村1號(hào)滑坡的Vb值約為2 875 073.8 m3。

4.3 基于滑體滑動(dòng)前后體積相等的原則調(diào)整微地形

理論上Vb=Va，若Vb與Va不相等，則可通過(guò)調(diào)整滑動(dòng)前斜坡的坡面形態(tài)，使得Vb=Va。具體調(diào)整方法為：當(dāng)Vb>Va時(shí)，將滑動(dòng)前斜坡的坡面形態(tài)調(diào)整為外凸型;當(dāng)VbVa時(shí)，增大調(diào)整點(diǎn)的高程值;Vb

由上述計(jì)算可知，三合村1號(hào)滑坡的Vb>Va，因此需按上述調(diào)整方法將滑動(dòng)前斜坡的坡面形態(tài)調(diào)整成凸坡，以滿足Vb=Va的要求。滑動(dòng)前斜坡的坡面形態(tài)調(diào)整示意圖如圖15所示，調(diào)整后滑動(dòng)前斜坡三維地形圖如圖16所示。

4.4 滑坡滑動(dòng)前斜坡坡度的驗(yàn)證

圖17是S-Z坐標(biāo)系下建立的調(diào)整后的斜坡剖面示意圖。通過(guò)對(duì)斜坡坡面地形數(shù)據(jù)的多項(xiàng)式擬合，可得到調(diào)整后斜坡坡面的數(shù)學(xué)表達(dá)式為：

Z=-5.247×10-5S2-0.218 4S+1 963 （4）

式中：S的定義域?yàn)椋?，533.7）。

通過(guò)對(duì)（4）進(jìn)行求導(dǎo)，得到斜坡坡面的斜率表達(dá)式：

Z′=-1.049 4×10-4S-0.218 4 （5）

利用反三角函數(shù)，可得到斜坡坡度的表達(dá)式為：

α=arctan（Z′） （6）

利用式（6），求得斜坡坡度的分布范圍在12°～15°之間。已知該滑坡形成于黃土梁之上，已有對(duì)黃土地貌單元斜坡坡度的統(tǒng)計(jì)分析表明黃土梁的斜坡坡度分布范圍在0.9°～44.1°之間［19］?？梢?，所確定的滑動(dòng)前斜坡的坡度分布范圍位于黃土梁統(tǒng)計(jì)的坡度范圍內(nèi)，因此本文給出的三合村1號(hào)滑坡滑動(dòng)前的斜坡坡角在合理區(qū)間內(nèi)。

5 三合村1號(hào)滑坡穩(wěn)定性驗(yàn)算

基于恢復(fù)后的三合村1號(hào)滑坡地形，利用AutoCAD、Geostudio軟件建立相應(yīng)的斜坡計(jì)算模型，并采用擬靜力法對(duì)滑動(dòng)前斜坡的計(jì)算模型進(jìn)行穩(wěn)定性分析。

5.1 計(jì)算模型的建立

根據(jù)三合村1號(hào)滑坡滑動(dòng)前的斜坡剖面圖和滑坡外擴(kuò)一定范圍內(nèi)的采樣點(diǎn)數(shù)據(jù)，利用AutoCAD、Geostudio軟件建立如圖18所示的滑坡計(jì)算模型。該滑坡屬于黃土層內(nèi)滑坡，因此計(jì)算模型的地層結(jié)構(gòu)可概化為均質(zhì)黃土結(jié)構(gòu)，黃土為第四系馬蘭黃土，斜坡的坡面形狀為外凸形。模型的水平向長(zhǎng)度約為900 m，豎直向高度約為200 m，單元網(wǎng)格的尺寸設(shè)置為3 m，共劃分了15 558節(jié)點(diǎn)，15 284個(gè)單元。依據(jù)黃土樣品的室內(nèi)實(shí)驗(yàn)結(jié)果，并參考黃土參數(shù)的經(jīng)驗(yàn)數(shù)據(jù)，將黃土層的物理力學(xué)參數(shù)列于表1。

5.2 斜坡初始應(yīng)力分析

將上述計(jì)算模型和土層的物理力學(xué)參數(shù)導(dǎo)入Geostudio軟件，利用該軟件的SIGMA/W模塊快速完成斜坡的初始應(yīng)力分析計(jì)算（自重作用下的應(yīng)力分布），計(jì)算結(jié)果如圖19、20所示，其中σY代表豎直向的土壓應(yīng)力，σX代表水平向土的側(cè)向壓應(yīng)力。該結(jié)果可作為三合村1號(hào)滑坡穩(wěn)定性分析的初始條件。

5.3 利用擬靜力法計(jì)算穩(wěn)定性系數(shù)

擬靜力法是將地震作用等效為作用于滑坡體質(zhì)心的對(duì)斜坡穩(wěn)定最不利的慣性力，根據(jù)極限平衡理論可計(jì)算出沿指定滑動(dòng)面的穩(wěn)定性系數(shù)［20-21］。該方法極大地簡(jiǎn)化了地震動(dòng)對(duì)斜坡體的作用過(guò)程，由于計(jì)算過(guò)程比較簡(jiǎn)單，因此被廣泛應(yīng)用于巖土工程設(shè)計(jì)領(lǐng)域。對(duì)于作用在每個(gè)土條質(zhì)心位置的水平向地震荷載，可按《水利水電工程邊坡設(shè)計(jì)規(guī)范》［22］推薦的公式進(jìn)行計(jì)算：

Fhi=αhξWiαi/g （7）

式中：Fhi為作用在第i條土條上的水平向地震荷載;ξ為地震作用的動(dòng)態(tài)分布系數(shù);αi為地震作用的折減系數(shù);Wi為第i個(gè)土條的重量;αh為地震動(dòng)峰值加速度。經(jīng)換算，施加在每個(gè)土條質(zhì)心位置的水平向地震加速度α為：

α=αhξαi （8）

按照規(guī)范要求，ξ，αi分別取值為1.0和0.25。已知三合村1號(hào)滑坡地處1920年海原8.5級(jí)特大地震IX度區(qū)，可取其水平地震動(dòng)峰值加速度為0.4g。利用Geostudio軟件的SLOPE/W模塊，分別采用瑞典條分法［23］、畢肖普條分法［24］和簡(jiǎn)布條分法［25］計(jì)算了三合村1號(hào)滑坡的穩(wěn)定性系數(shù)，計(jì)算結(jié)果列于表2。

由表2可知，三合村1號(hào)滑坡在0.4g的水平向地震動(dòng)峰值加速度作用下，采用不同方法計(jì)算的穩(wěn)定性系數(shù)均小于1，說(shuō)明在該地震動(dòng)作用下會(huì)發(fā)生失穩(wěn)破壞。這與實(shí)際情況吻合，在一定程度上說(shuō)明了本文地震古/老滑坡原始地形恢復(fù)方法的合理性。

6 結(jié)論

黃土地震古/老滑坡原始地形恢復(fù)是開展地震古/老滑坡研究的一項(xiàng)基礎(chǔ)性和必不可少的工作，開展這項(xiàng)工作的難度較大。本文在已有研究工作的基礎(chǔ)上［26］，利用滑坡地形資料，基于滑坡體滑動(dòng)前后體積相等的原則，以MATLAB軟件為操作平臺(tái)，初步提出了一種黃土地震古/老滑坡原始地形恢復(fù)方法，并利用該方法恢復(fù)了寧夏回族自治區(qū)西吉縣平峰鄉(xiāng)三合村1號(hào)滑坡的原始地形，驗(yàn)證了地形恢復(fù)結(jié)果的合理性和可行性。

本文提出的恢復(fù)方法具有方便、快捷的優(yōu)點(diǎn)，且考慮到了較多的因素，恢復(fù)精度較以往的方法有所提高，適用于區(qū)域性地震黃土滑坡研究，對(duì)黃土古/老滑坡的研究有一定的理論意義和工程應(yīng)用價(jià)值。但該方法在滑動(dòng)面位置、滑坡體體積、斜坡坡面具體凹凸形態(tài)的確定等方面還有待深入開展研究工作，在精度上仍需進(jìn)一步提高。

參考文獻(xiàn)（References）

［1］ 王蘭民.黃土動(dòng)力學(xué)［M］.北京：地震出版社，2003：144-150.

WANG Lanmin.Loess dynamics［M］.Beijing：Seismological Press，2003：144-150.

［2］ 袁道陽(yáng)，雷中生，王愛國(guó).1654年甘肅天水南8級(jí)地震補(bǔ)充考證［J］.地震工程學(xué)報(bào)，2017，39（3）：509-520.

YUAN Daoyang，LEI Zhongsheng，WANG Aiguo.Additional textual criticism of southern Tianshui M8 earthquake in Gansu Province in 1654［J］.China Earthquake Engineering Journal，2017，39（3）：509-520.

［3］ 趙云.1654年甘肅天水地震烈度和震級(jí)考析［J］.高原地震，2008，20（2）：30-33，29.

ZHAO Yun.Textual research on the intensity and magnitude of Gansu Tianshui earthquake in 1654［J］.Plateau Earthquake Research，2008，20（2）：30-33，29.

［4］ 彭達(dá)，楊順，李孝波.甘肅通渭縣黃土地震滑坡分布特征及發(fā)育機(jī)理［J］.中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào)，2017，28（3）：31-38.

PENG Da，YANG Shun，LI Xiaobo.Distribution and mechanism of loess landslides induced by earthquake in Tongwei，Gansu Province［J］.The Chinese Journal of Geological Hazard and Control，2017，28（3）：31-38.

［5］ 郭增建，謝原定，李孟鑾，等.1920年12月16日的海原大地震［J］.地球物理學(xué)報(bào)，1976，19（1）：42-49.

GUO Zengjian，XIE Yuanding，LI Mengluan.Great Haiyuan （海原） earthquake on December 16，1920［J］.Chinese Journal of Sinica，1976，19（1）：42-49.

［6］ 鄧龍勝，范文.寧夏海原8.5級(jí)地震誘發(fā)黃土滑坡的變形破壞特征及發(fā)育機(jī)理［J］.災(zāi)害學(xué)，2013，28（3）：30-37.

DENG Longsheng，F(xiàn)AN Wen.Deformation breakage characteristics and development mechanism of loess landslide triggered by Haiyuan M8.5 earthquake in Ningxia［J］.Journal of Cata-strophology，2013，28（3）：30-37.

［7］ ZHUANG J，PENG J，XU C，et al.Distribution and characteristics of loess landslides triggered by the 1920 Haiyuan earthquake，Northwest of China［J］.Geomorphology，2018，314：1-12.

［8］ 《工程地質(zhì)手冊(cè)》編委會(huì).工程地質(zhì)手冊(cè)［M］.5版.北京：中國(guó)建筑工業(yè)出版社，2018：651-657.

Editorial Committee of Geological Engineering Handbook.Geological engineering handbook［M］.5th ed.Beijing：China Architecture & Building Press，2018：651-657.

［9］ 李忠生.黃土滑坡的滑前地貌重建［J］.自然災(zāi)害學(xué)報(bào)，2004，13（2）：96-100.

LI Zhongsheng.Pre-failure topography reconstruction for loess landslide［J］.Journal of Natural Disasters，2004，13（2）：96-100.

［10］ 李萍.黃土邊坡可靠性研究［D］.西安：長(zhǎng)安大學(xué)，2006.

LI Ping.Reliability applied to stability analysis of loess slope［D］.Xi'an：Chang'an University，2006.

［11］ CONOSCENTI C，CIACCIO M，ALMARU CARABALLO-ARIAS N，et al.Assessment of susceptibility to earth-flow landslide using logistic regression and multivariate adaptive regression splines：a case of the Belice River Basin （western Sicily，Italy）［J］.Geomorphology，2015，242：49-64.

［12］ 鄭旭輝.邊坡地震動(dòng)力失穩(wěn)概率分析方法研究：以寧夏西吉震湖斜坡為例［D］.北京：中國(guó)地質(zhì)大學(xué)（北京），2018.

ZHENG Xuhui.Study on probability analysis method of slope seismic dynamic instability：a case study of Xiji Zhenhu slope in Ningxia［D］.Beijing：China University of Geosciences，2018.

［13］ 周鑫，周慶，高帥坡，等.重慶小南?；略嫉匦位謴?fù)及滑坡體體積計(jì)算［J］.地震地質(zhì)，2020，42（4）：936-954.

ZHOU Xin，ZHOU Qing，GAO Shuaipo，et al.Restoration of the original topography of the Xiaonanhai landslide in Chongqing and calculation of its volume［J］.Seismology and Geology，2020，42（4）：936-954.

［14］ 王廣，邢林芳.MATLAB GUI程序設(shè)計(jì)［M］.北京：清華大學(xué)出版社，2018：115-128.

WANG Guang，XING Linfang.MATLAB GUI programming［M］.Beijing：Tsinghua University Press，2018：115-128.

［15］ 邢麗.三次樣條插值端點(diǎn)約束條件的構(gòu)造與Matlab實(shí)現(xiàn)［J］.上海第二工業(yè)大學(xué)學(xué)報(bào)，2012，29（4）：319-323.

XING Li.Structure and Matlab implementation of cubic spline interpolation endpoint constraints［J］.Journal of Shanghai Second Polytechnic University，2012，29（4）：319-323.

［16］ 鄭穎人，陳祖煜，王恭先.邊坡與滑坡工程治理［M］.2版.北京：人民交通出版社，2010：16-20.

ZHENG Yingren，CHEN Zuyu，WANG Gongxian.Engineering treatment of slope & landslide［M］.2nd ed.Beijing：China Communications Press，2010：16-20.

［17］ 成永剛.滑坡區(qū)域性分布規(guī)律與防治［M］.北京：人民交通出版社，2014：24-28.

CHENG Yonggang.Regional distribution rules and prevention of landslides［M］.Beijing：China Communications Press，2014：24-28.

［18］ 王振宇，張常亮，周瀟朗，等.土質(zhì)滑坡滑前三維地形恢復(fù)研究［J］.甘肅科學(xué)學(xué)報(bào)，2020，32（1）：72-78.

WANG Zhenyu，ZHANG Changliang，ZHOU Xiaolang，et al.Study on the restoration of three-dimensional terrain before the landslide［J］.Journal of Gansu Sciences，2020，32（1）：72-78.

［19］ 吳良杰，薄景山，牛潔.黃土地區(qū)地形特征參數(shù)分布特征的統(tǒng)計(jì)分析［J］.科學(xué)技術(shù)與工程，2021，21（21）：8797-8806.

WU Liangjie，BO Jingshan，NIU Jie.Statistical analysis on the distribution characteristics of topographical parameters in loess area［J］.Science Technology and Engineering，2021，21（21）：8797-8806.

［20］ 薄景山，段玉石，常晁瑜，等.斜坡地震穩(wěn)定性研究的若干問(wèn)題［J］.自然災(zāi)害學(xué)報(bào)，2019，28（1）：1-8.

BO Jingshan，DUAN Yushi，CHANG Chaoyu，et al.Some problems of study on slope stability under earthquake［J］.Journal of Natural Disasters，2019，28（1）：1-8.

［21］ 唐世雄，鄧廣輝，劉衡秋.斜坡抗震穩(wěn)定性分析中擬靜力法的應(yīng)用研究［J］.城市地質(zhì)，2020，15（2）：181-186.

TANG Shixiong，DENG Guanghui，LIU Hengqiu.Study on the application of pseudo-static method in the analysis of seismic stability of slope［J］.Urban Geology，2020，15（2）：181-186.

［22］ 中華人民共和國(guó)水利部.水利水電工程邊坡設(shè)計(jì)規(guī)范：SL 386—2007［S］.Beijing：China Water & Power Press.

Ministry of Water Resources of the People's Republic of China.Design code for engineered slopes in water resources and hydropower projects：SL 386—2007［S］.Beijing：China Water & Power Press.

［23］ 東南大學(xué)，浙江大學(xué)，湖南大學(xué).土力學(xué)［M］.4版.北京：中國(guó)建筑工業(yè)出版社，2016.

Southeast University，Zhejiang University，Hunan University.Soil mechanics［M］.4th ed.Beijing：China Architecture & Building Press，2016.

［24］ BISHOP A W.The use of the slip circle in the stability analysis of slopes［J］.Géotechnique，1955，5（1）：7-17.

［25］ JANBU N，BJERRUM L，KJAERNSLI B，et al.Soil mechanics applied to some engineering problems［R］.Osio：Norwegian Geotechnical Institute Publication，1956：25-31.

［26］ 吳良杰.黃土古老滑坡原始地形恢復(fù)方法及其應(yīng)用研究［D］.廊坊：防災(zāi)科技學(xué)院，2021.

WU Liangjie.Research on the original terrain restoration method of loess ancient landslide and its application［D］.Langfang：Institute of Disaster Prevention，2021.