劉裕華，陳征宙，蔣 鑫

（南京大學 地球科學與工程學院，南京 210093）

1 引 言

滑坡是僅次于地震的第二大地質災害，對地震、火山、降雨等因素對滑坡的觸發(fā)機制研究以及滑坡發(fā)生時間的分布特征是當前滑坡研究的熱點之一，李生才等[1]對地質災害的統(tǒng)計分析結果表明滑坡呈現(xiàn)明顯的時間分布規(guī)律。在滑坡預報方面，除了采用一些有效方法外[2－3]，也可從滑坡的時間分布規(guī)律入手，樊曉一[4]對三峽水庫區(qū)滑坡發(fā)生時間進行了R/S分析研究，得出庫區(qū)滑坡發(fā)生時間的分布規(guī)律；李軍等[5]從滑坡統(tǒng)計數(shù)據(jù)的角度上詳細分析了香港地區(qū)滑坡災害在不同年份月份以及時刻的分布特征。目前的研究集中在以年份或月份等季節(jié)性方面[6]，對日周期方面的研究甚少。而滑坡發(fā)生時間的分布特征是由其內在規(guī)律所決定的，除了與降雨等季節(jié)性周期密切相關外，有一種觀點認為天體引力的作用是造成滑坡有日周期特征的原因，如山田剛二等[7]研究發(fā)現(xiàn)滑坡在地球自轉變化的年份明顯多于正常的年份；王運生等[8]也做過地球自轉及日月引潮力與滑坡災害發(fā)育的相關性研究。

固體潮相關地質問題的研究一直以來都集中在地震、火山、間歇泉活動等方面，是由于日月的引力，在地殼內產(chǎn)生一個附加應力Δσ，并定義其為附加潮汐應力，Δσ的量級與構造應力相比是很小的，在應變能積蓄的過程中不能產(chǎn)生地震觸發(fā)作用，但在活動構造帶的應力將要達到巖體強度極限的時候，斷裂帶鄰近破壞點，則Δσ對地震產(chǎn)生觸發(fā)作用。關于這方面的論述見諸于國內大量的文獻當中，尤其是淺源小震，不斷被人們通過固體潮監(jiān)測數(shù)據(jù)與地震統(tǒng)計數(shù)據(jù)對比分析而得到證實[9－19]。

對滑坡這種淺表生地質災害而言，附加應力將與坡體原有應力迭加，使得坡體總應力受到潮汐應力的影響。以往有關滑坡發(fā)生時間分布特征的研究中，大多通過對滑坡的統(tǒng)計分析結果建立相關關系及提出初步的看法。然而，從附加潮汐應力的角度深入研究其觸發(fā)機制及日周期特征的較少。本文從起潮力位和起潮力場著手，推導地表附加潮汐應力各分量，探討其對滑坡體的影響機制及其與時間的相關關系。

2 附加潮汐應力的推導

2.1 計算模型

如圖1所示，取地球為簡化地球模型，即假設地球為各向同性的均勻彈性球體，在地表任取一點P點，其與月心的距離為l，地心與月心的距離為R，地球半徑為 r，的夾角為θ，并在 P點建立沿南北方向為x，東西方向為y，半徑方向為z的坐標系。

圖1 計算模型Fig.1 Calculation model

2.2 附加潮汐應力推導

根據(jù)固體潮汐理論，P點的起潮力位Wp的表達式為

式中：G為引力常量；M為吸引地球的天體（如月亮或太陽）質量；，代入式（1）可得

半徑r0、南北方向x、東西方向y方向的起潮力的分別為

將式（2）、（3）代入下述平衡方程：

式中：u、v、w分別為x、y、z方向上的潮汐位移；λ、μ分別為拉梅常數(shù)；ρ為地球密度；θ為與的夾角，

采用二階起潮力位求解，可以得到球坐標系下附加潮汐應力的表達式：

式中：r0為地球半徑。

由于邊坡為地球表生地質災害，即r≈r0，代入式（5）可得潮汐應力：

由式（6）可見，南北向的潮汐應力σx為一常量，而東西向的的潮汐應力σy為以θ值為變量的函數(shù)，可寫成如下表達式：

式中：θ為隨P點的位置、地球自轉、天體M的運動位置而變化的變量，它們的位置隨時間而不同，即可由矢量的乘積可知：

式中：α、δ分別為P點的經(jīng)緯度。

由坐標轉換有

式中：Rz(-SG)為矢量′繞 z軸旋轉-SG度時的轉換，

SG為格林尼治恒星時，

在計算精度要求不高時可將SG計算式中的 2次方和3次方項略去。

同理，有

式中：Rx(-ε)為矢量繞 z軸旋轉-ε度時的轉換，

ε為黃赤交角，ε=23?26′ 21′.4 48-46′.8150T 。

關于月球在赤道坐標系下的(X Y Z)坐標可以參考David的求解公式[19]：

上述的T均為儒略日，指由公元前4713年1月1日，協(xié)調世界時中午12：00開始所經(jīng)過的天數(shù)（day），多為天文學家采用，用以作為天文學的單一歷法，把不同歷法的年表統(tǒng)一起來。也上述式中須將公歷年（year）月（mouth）日（day）T轉化為儒略日JDN，轉換公式為

式中：a、y、m的取值為

將式(10)、（15）～（17）代入式（8），即可得θ與時間T的函數(shù)式，再將式（8）代入式（7），最終將東西向的潮汐應力σy轉化為以時間及經(jīng)緯度為變量的函數(shù)，即

式中：α、δ分別為P點的經(jīng)緯度。

3 潮汐應力觸發(fā)滑坡的作用方式分析

（1）瑞典條分法

將滑動土體豎直分成若干條，并把土條當成剛體，分別求作用于各土條上的力對圓心的滑動力矩和抗滑力矩，安全系數(shù)Fs被定義為兩者之比，即

若令剛體土條在在滑帶處所受的附加潮汐應力為Δσ，則可通過將此附加應力分解到平行于滑面和垂直于滑面兩個方向，從而將其對坡體安全系數(shù)的影響考慮進去，得到考慮附加潮汐應力影響的安全系數(shù)的公式（22），受力分析如圖2(a)所示。

由式（22）可知，附加潮汐應力會改變坡體安全系數(shù)，由于潮汐應力與時間有關，因此安全系數(shù)也是隨時間變化的函數(shù)，當安全系數(shù)減小至接近1.0時，將促使邊坡滑動。

圖2 潮汐應力影響下的邊坡各破壞模式的作用力分析Fig.2 The forces analysis with slope failure modes by taking tidal stress into account

（2）折線滑動法

通過勘查定出滑坡的滑動面是折線形，根據(jù)折線將滑體條分后進行計算的一種方法。由圖2（b），加入附加潮汐應力后的穩(wěn)定性系數(shù)Fs公式為

孔隙壓力比：

（3）平面滑動法

巖質邊坡中常見的破壞模式，將被單結構面切割的滑體視為剛體，因此該滑塊的穩(wěn)定性系數(shù)Fs被定義為抗滑力與下滑力之比，即

將附加潮汐應力Δσ加入到滑體的受力分析中，見圖 2(c)，得到考慮潮汐應力后的滑塊穩(wěn)定性安全系數(shù)：

（4）楔形滑動

巖質邊坡中另一常見的破壞模式，滑塊被 AB 2組結構面切割，在認為滑塊將沿著兩結構面組合交線滑動的情況下，該滑塊的穩(wěn)定性系數(shù)Fs被定義為抗滑力與下滑力之比，即

將附加潮汐應力Δσ加入到楔形體的受力分析中（見圖2(d)），得到考慮潮汐應力的楔形體穩(wěn)定性安全系數(shù)，式中的參數(shù)如圖2所示。

4 算 例

以湖北省宜昌市秭歸縣新灘鎮(zhèn)大型滑坡作為算例，東經(jīng) 110°49′26.58″，北緯為 30°55′53.06″，堆積斜坡近南北向展布，長約2 km，北窄南寬，總面積近1.1 km2。后緣高程約900 m，向南伸展至高程約 65 m的江邊，平均坡度為 23°?；瑤劣杀婪e物、崩坡積物、沖積與崩坡積混合堆積物組成，一般厚 30～40 m，最厚達 86 m，總體積約為 3×107m3。下臥基巖，從東向西依次為志留系砂頁巖、泥盆系砂巖、石炭系灰?guī)r及二疊系灰?guī)r與煤層（見圖 3）?；禄瑒用鏋橥翈r分界面，滑帶土計算參數(shù)中黏聚力 c=10 kPa，內摩擦角φ=17.7°。該滑坡的平面形態(tài)示意圖見圖 4，主滑段工程地質剖面圖見圖3。采用條分法的計算模型見圖5（1985年 5月20日～6月20日），滑坡主滑段條分計算參數(shù)見表1。

圖3 滑坡主滑段工程地質剖面圖Fig.3 The engineering geological profile of main section

圖4 滑坡平面形態(tài)圖（單位：m）Fig.4 The sketch of morphology of the research area (unit: m)

圖5 潮汐應力與時間的關系Fig.5 Tidal stress vs. time

表1 滑坡主滑段條分計算參數(shù)Table1 The parameters for calculation of the main sliding section of the landslide

滑坡于1985年6月12日3:45發(fā)生滑動，原因是后緣荷重和斜坡物質的積累，根據(jù)剩余推力公式計算出的整個堆積物的剩余水平推力為8.0×106kN。主要誘因是降雨，根據(jù)薛果夫等的研究，新灘斜坡的位移滯后于雨季起止時間，且不同的變形階段，位移與降雨的關系也不同[20]。本文考慮潮汐應力的因素，取計算起始時間為1985年5月20日，終止時間為1985年6月20日，計算單位為d。集中計算1985年6月11日12:00至1985年6月12日12:00之間潮汐應力，計算單位為 h。按照上述公式編制程序進行計算，附加潮汐應力的計算結果見圖 6。1985年6月11日12:00至1985年6月12日12:00的計算結果見圖7。

由圖5、6可知，潮汐應力呈明顯周期性變化，從30 d的計算結果來看，最大的周期為15 d，但上半6708與下半月的各峰值仍有一定差異。另外，日周期以及半日周期也十分明顯，但每天的峰值大小有差異。這與地球自轉以及月球繞地球運動的位置有密切關系。另外值得注意的是，由于潮汐應力南北方向σx為一常量，而東西分量σy的大小及方向均變化明顯，并影響潮汐應力總量的變化。

由圖6可知，滑坡發(fā)生的時間（1985年6月12日3:45）與潮汐應力處于峰值的時間較為接近（1985年6月12日4:00），說明除了降雨這個主要觸滑因素之外，潮汐應力也可能是觸發(fā)因素之一。

根據(jù)折線法將主滑段坡體分為 27個條塊并采用式（23）進行計算，條分圖見圖7。參數(shù)見表1。計算結果為：在不考慮潮汐應力的作用下安全系數(shù)的計算結果為1.03。在同等條件下，考慮潮汐應力作用后（計算時間，每1 h取一個計算點），計算結果如圖8所示，由圖可知，邊坡安全系數(shù)隨時間呈明顯日周期和半日周期性變化，最大與最小安全系數(shù)的差值與不考慮潮汐應力影響的安全系數(shù)的比值達3.4%。在2008年9月23日18:00和24日的3:00，安全系數(shù)達到極小值，與新灘滑坡發(fā)生的時間較為接近。由此推測，潮汐應力很可能在降雨等其他主要滑動因素的前提下起到了觸發(fā)滑動的作用。

圖6 潮汐應力與時間的關系Fig.6 Tidal stress vs. Time

圖7 滑坡主滑段條分圖Fig.7 The slice computational module of the landslide

圖8 潮汐應力作用下邊坡穩(wěn)定性系數(shù)隨時間的變化Fig.8 FOS vs. time by taking the tidal stress into account

5 結 論

（1）附加潮汐應力是與經(jīng)度、緯度以及時間為變量的函數(shù)，通過對某給定經(jīng)緯度的地點進行附加潮汐應力的計算，可以了解其變化規(guī)律以及變化周期。

（2）潮汐應力對臨界狀態(tài)下邊坡的觸發(fā)機制可根據(jù)從安全系數(shù)的表達式的改進來解釋，也可以從巖體破壞判據(jù)來進行解釋，即附加潮汐應力以4種不同的情況使坡體滑帶應力摩爾圓與破壞線相切，使坡體破壞。

（3）通過對某滑坡進行30 d和3 d計算，得到附加潮汐應力的周期性變化情況，其中半月周期、日周期以及半日周期最為明顯。該滑坡發(fā)生的時間與潮汐應力出現(xiàn)峰值的時間較為接近。邊坡穩(wěn)定性安全系數(shù)隨時間的關系同樣呈現(xiàn)明顯日周期和半日周期變化，滑坡發(fā)生的時間與安全系數(shù)處于最低值的時間較接近。

（4）從固體潮理論上對月亮潮汐應力與時間的相關性以及對滑坡的觸發(fā)機理做出分析，用以闡述除降雨等明顯季節(jié)因素外，滑坡發(fā)生具有時間分布特征的另一重要影響因子。邊坡的日周期因子研究是一個復雜的問題，需要結合足夠多的監(jiān)測資料做進一步詳細地研究。

[1]李生才,王亞軍,黃平. 2005年7－8月國內環(huán)境事件數(shù)據(jù)[J]. 安全與環(huán)境學報,2005,4(5): 120－124.LI Sheng-cai,WANG Ya-jun,HUANG Ping. Statistics of environmental events in China during period from July to August in 2005[J]. Journal of Safety and Environment,2005,4(5): 120－124.

[2]徐金明,孫昆侖. 流體包裹體法在巖體滑坡周界預測中的應用[J]. 土木工程學報,2007,40(1): 69－73.XU Jin-ming,SUN Kun-lun. Determination of landslide boundary by fluid inclusions[J]. China Civil Engineering Journal,2007,40(1): 69－73.

[3]蘭恒星,伍法權,周成虎,等. GIS支持下的降雨型滑坡危險性空間分析預測[J]. 科學通報,2003,48(5): 507－512.LAN Heng-xing,WU Fa-quan,ZHOU Cheng-hu,et al.Spatial analysis and prediction of landslide hazard in rainfall using GIS technology[J]. Science Bulletin,2003,48(5): 507－512.

[4]樊曉一,喬建平. 三峽水庫區(qū)滑坡時間記錄的 R/S分析[J]. 中國地質災害與防治學報,2006,17(2): 99－101.FAN Xiao-yi,QIAO Jian-ping. R/S analysis of landslide time record in the Three Gorges reservior area[J].Chinese Journal of Geological Hazard and Control,2006,17(2): 99－101.

[5]李軍,周成虎,許增旺. 香港島地區(qū)滑坡災害的時空分布模式[J]. 山地學報,2001,19(3): 248－252.LI Jun,ZHOU Cheng-hu,XU Zeng-wang. Temporal and spatial distribution patterns of landslides in Hong Kong island[J]. Journal of Mountain Science,2001,19(3):248－252.

[6]楊媛媛,黃宏偉. 四川省滑坡災害記錄分析[J]. 巖石力學與工程學報,2005,24(增刊2): 5366－5370.YANG Yuan-yuan,HUANG Hong-wei. Landslide hazard probability analysis for Sichuan province using landslide inventory[J]. Chinese Journal of Rock Mechanics and Enginveering,2005,24(Supp.2): 5366－5370.

[7]山田剛二,渡正亮. 滑坡和斜坡崩塌及防治[M]. 《滑坡和斜坡崩塌及其防治》翻譯組譯. 北京: 科學技術出版社,1980: 128－186.

[8]王運生,王士天. 地球自轉、日月引潮力與滑坡災害發(fā)育的相關性研究[J]. 大自然探索,1997,16(4): 55－58.WANG Yun-sheng,WANG Shi-tian. A study of the relation between tidal stress,rotation of the earth and landslides[J]. Exploration of Nature,1997,16(4): 55－58.

[9]吳小平,冒蔚,黃雍,等. 潮汐應力對發(fā)震斷層作用的統(tǒng)計分析[J]. 地球物理學報,1999,42(增刊): 65－74.WU Xiao-ping,MAO Wei,HUANG Yong,et al.Statistical analysis of tidal stress effect on seismic faults[J]. Chinese Journal of Geophysics,1999,42(Supp.): 65－74.

[10]楊學祥,韓延本,陳震,等. 強潮汐激發(fā)地震火山活動的新證據(jù)[J]. 地球物理學報,2004,47(4): 69－74.YANG Xue-xiang,HAN Yan-ben,CHEN Zhen,et al.New evidence of earthquake and volcano triggering by strong tides[J]. Chinese Journal of Geophysics,2004,47(4): 69－74.

[11]吳小平,黃雍,冒蔚,等. 云南地震的潮汐應力觸發(fā)機制及相關天體位置圖像[J]. 地球物理學報,2005,48(3):574－583.HUANG Xiao-ping,HUANG Yong,MAO Wei,et al.Tidal stress mechanism of earthquakes in Yunnan and related of celestial body positions[J]. Chinese Journal of Geophysics,2005,48(3): 574－583. .

[12]張晶,郗欽文,楊林章,等. 引潮力與潮汐應力對強震觸發(fā)的研究[J]. 地球物理學報,2007,50(2): 126－132.ZHANG Jing,XI Qin-wen,YANG Lin-zhang,et al. A study on tidal force/stress triggering of strong earthquakes[J]. Chinese Journal of Geophysics,2007,50(2): 126－132.

[13]周擎,山秀明,張立,等. 基于HHT提取昆明、下關重力固體潮的地震前兆信息[J]. 地球物理學報,2008,51(3): 836－844.ZHOU Zhi,SHAN Xiu-ming,ZHANG Li,et al. The gravity tide of Kunming and Xiaguan based on the HHT[J]. Chinese Journal of Geophysics,2008,51(3)836－844.

[14]ELI Z S,JOHN E V,SACHTKO T. Earth tides can trigger shallow thrust fault earthquakes[J]. Science,2004,306:1164－1166.

[15]LAURENT M,OLIVIER D V,CLINTON R,et al.Evidence of earthquake triggering by the solid earth tides[J]. Earth and Planetary Science Letters,2009,278: 370－375.

[16]CAI Z X. A continuous crustal stress monitoring method for earthquake prediction[J]. Pure and Applied Geophysics,165(9-10): 1879－1889.

[17]STROUP D F,BOHNENSTIEHL D R,TOLSTOY M,et al. Pulse of the seafloor: Tidal triggering of microearthquakes at 90 degrees 50’N East Pacific Rise[J].Geophysical Research Letters,34(15): 6.

[18]TANAKA S,SATO H,MATSUMURA S,OHTAKE M.Tidal triggering of earthquakes in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region,Japan[J]. Tectonophysics,417(1-2): 69－80.

[19]DAVID G S. An alternative lunar ephemeris model for on-board flight software use[R]. Maryland: Flight Mechanics Symposium,NASA GSFC,1999.

[20]薛果夫,呂貴芳,任江. 新灘滑坡研究[M]. 北京: 科學出版社,1988: 48－54.