基于FLAC3D模擬庫水位改變對周圍傾斜場的影響★

李宏偉 王 霞 李文超

(1.山西省地震局，山西 太原 030021； 2.山西省地震局大同中心地震臺，山西 大同 037000)

?

·水利工程·

基于FLAC3D模擬庫水位改變對周圍傾斜場的影響★

李宏偉1王 霞1李文超2

(1.山西省地震局，山西 太原 030021； 2.山西省地震局大同中心地震臺，山西 大同 037000)

采用FLAC3D軟件，對郭莊水庫庫水位的改變引起周圍傾斜場的變化進行了分析，結果表明，2013年7月庫水位的上升對昔陽傾斜儀測點處的影響約為0.01角秒，遠遠小于實際觀測值的2角秒，因此該測點2013年7月的大幅度變化不是由庫水位的上升造成的。

水庫庫水位，FLAC3D，傾斜場，位移量

地震前兆異常跟蹤過程中，最重要的核心問題是科學判定異常。如何從定性到定量分析異常變化與干擾的影響規律，為地震預測研究提供堅實可靠的依據，是地震預測研究的探索和努力方向。

隨著基礎設施建設的增多，各個臺站的觀測環境受到越來越嚴重的影響，臺站周邊存在著越來越多的載荷影響[1，2]，其中包括水庫庫水位的改變，這些干擾在某種程度上可以等效為載荷的變化，許多學者也從理論上分析了載荷對周邊應變場觀測的影響[3-5]，也有不少學者利用數值分析技術探討水庫庫水位的改變對周圍巖體穩定性的影響[6-9]，水庫庫水位的改變確實可以引起周圍應力場的改變，但定量化的研究對周圍傾斜場影響的文獻還涉及不多，鑒于此，采用數值模擬的方法對郭莊水庫庫水位上升引起的傾斜場變化進行模擬計算。

1 水庫及測點概況

郭莊水庫位于昔陽縣城南的郭莊村附近松溪河上游，屬海河流域子牙河水系，設計水庫庫容為2 000萬m3左右的中型水庫[10]，2013年7月，8月，9月水庫進、出水量如表1，表2所示。

表1 郭莊水庫2013年進庫水量

表2 郭莊水庫2013年出庫水量

昔陽水平擺傾斜儀安放于昔陽地震臺山洞中，山洞開挖于奧陶系中統上馬家溝組深灰色厚層石灰巖中，山洞覆蓋厚度約12 m，每年雨季洞內濕度增大，從近幾年的觀測資料來看，每年雨季曲線均出現大幅臺階變化。2013年7月東西向最大傾斜角度達2 000毫角秒左右(見圖1)，遠超過正常的背景值，觀測點西南方向約4.63 km處的郭莊水庫在測點傾斜變化前有約300萬m3水量進庫，因此初步判斷測點大幅度傾斜變化與水庫庫容的增加有關，但具體原因仍需進一步論證。

2 水庫庫容改變影響的定量分析

2.1 三維地質模型的構建

將水庫庫水位上升變化等效為在水庫所在位置的上方額外增加載荷，將水庫庫水位的上升量轉化為豎直向下的法向應力，相當于均勻載荷的加載。由表1，表2可知，水庫進、出水量的最大變化為2013年7月份的進庫量，達310.91萬m3，為了獲取最大影響范圍及影響量級，以2013年7月份的變化量為計算依據來模擬獲取水庫荷載的變化量，經計算得到等效載荷約為3.1×1010N。

采用FLAC3D數值模擬軟件中的Brick六面塊體網格單元進行建模，該單元由8個節點組成；設定模型大小為8 km×8 km×4 km，單個網格長度為200 m，單元個數為3.2×104個，水庫面積簡化為(750×750)m2，建模時將水庫置于模型中央，在模型底部和側面施加位移約束，水庫表面上施加面載荷。

觀測點位于太行山斷塊隆起區，周圍出露的地層巖性以灰巖和砂巖為主；取彈性模量E=3.0×104MPa，泊松比μ=0.25，比重ρ=2.7 kg/m3。

2.2 位移量計算

圖2是加載計算平衡后的垂直向位移分布圖，可以看出，在水庫上方加載垂向應力后所引起的垂直向位移場是以水庫為中心，快速向水庫外圍衰減，最大垂向位移為1.87×10-4m。

圖3是加載計算平衡后沿X方向的水平位移圖，可以看出，面載荷模型加載所引起的水平向位移場最大位移點不是垂向最大位移點，而是存在一定的水平偏離，其中最大水平位移量為3.0×10-5m。

2.3 傾斜量計算

任意點M(x，y，z)的傾斜角度α可以通過質點的水平位移Wxy與垂直位移Wz的比值求反正切值計算(見圖4，式(1))。由以上計算得到的垂直位移場與水平位移場分布可知，水庫庫水位的改變引起的最大傾斜點位置不在水庫中心的位置，而應位于離開水庫中心一定的位置。

(1)

式中:α——傾斜角度；Wxy——M點相對于變形前的水平方向的位移；Wz——M點相對于變形前的垂向位移。

圖5是經過計算得到的傾斜場，結果顯示，傾斜量的最大值位置與水平位移量最大值位置一致，傾斜量自最大處向載荷中心與外圍逐漸縮小，經計算臺站所在處的傾斜量約為0.01角秒。

3 結語

地表傾斜監測作為重要的地震前兆監測手段，在地震監測預

報研究中發揮著重要的作用，尤其在地震短臨預報研究中取得了許多研究成果，為地震預報研究積累了豐富的觀測資料。在各類干擾因素中，載荷的變化是其主要的影響因素之一，如何定量分析載荷的變化對觀測值的影響規律，是進行分析預報的首要工作。

根據構建的FLAC3D模型計算結果得出，郭莊水庫2013年7月庫水 位改變對傾斜監測點的影響遠遠小于實際監測值，因此認為昔陽傾斜儀測點在2013年7月的大幅度變化與郭莊水庫庫水位的上升關系不大，可能同時段還有其他干擾因素的存在。據文獻[11]可知，觀測值大幅變化前測區附近有降雨，根據數值模擬計算結果得出的傾斜值與實際觀測較為吻合，據此推測此次變化可能受降雨的影響。

水庫庫水位改變對周圍傾斜場的影響在實際分析中，不能僅從時間進程的角度去給出判定，還應通過構建理論模型定量分析，給出影響范圍及量級，但模型結果的可靠性與地質模型構建的精細程度、參數的選取以及荷載變化量的多少有關。

[1] 朱 航,蘇 琴,楊 濤，等.耿達短水準觀測資料在汶川MS8.0地震前后異常的辯別[J].地震學報，2010，32(6)：649-658.

[2] 馬 棟,陳建國,鄭 爽，等.荷載對鉆孔應變觀測影響的實例分析[J].地震研究，2014，37(1)：79-81.

[3] 呂 芳，張淑亮，胡玉良，等.山西靜樂井水位異常調查與分析[J].山西地震，2013(3)：28-31.

[4] 張惠蘭.建筑物載荷位置及大小對地表及隧道的影響研究[J].交通與建筑科學，2013，5(10)：149-151.

[5] 胡衛建,張俊山,謝 智，等.載荷對鉆孔應變測值影響的實驗動力學解析[J].地震，2002，22(3)：95-104.

[6] 邱澤華.鉆孔應變觀測點到地面載荷干擾源最小“安靜”距離的理論分析[J].巖石力學與工程學報，2004，24(4)：80-86.

[7] 鄭爭鋒,龔 輝.水位漲落對庫岸滑坡穩定性影響研究[J].山西建筑,2016,42(12):56-57.

[8] 杜飛翔,張德成.某水庫右岸斜坡體穩定性分析研究[J].山西建筑,2016,42(16):75-77.

[9] 宋丹青.水庫蓄水對庫岸邊坡穩定性的影響[D].蘭州：蘭州大學,2015.

[10] 宋丹青,王 豐,梅明星,等.水庫蓄水對庫岸邊坡穩定性的影響[J].鄭州大學學報(工學版),2016(1):1.

[11] 劉黎明.針梁式鋼模臺車在郭莊水庫泄洪洞的應用[J].山西水利科技,2014(4):27-28.

[12] 李宏偉,劉瑞春,王 霞，等.山西昔陽地震臺水平擺觀測受降雨影響的數值分析[J].山西地震,2016(2):1-4.

The simulation of the change in the water-level of reservoir influencing tilting field based on FLAC3D★

Li Hongwei1Wang Xia1Li Wenchao2

(1.EarthquakeAdministrationofShanxiProvince,Taiyuan030021,China; 2.DatongSeismologicalofEarthquakeAdministrationofShanxiProvince,Datong037000,China)

The paper analyzes the change of the water-level of Guozhuang reservoir causing the variation of tilting field through FLAC3D. The result shows that the influence on the tiltmeter of Xiyang station caused by the rising of reservoir water-level in July, 2013 is about 0.01 arcsecond, which is far less than the observation value of 2 arcsecond. Therefore, the significant change of the tiltmeter of Xiyang station is not caused by the rising of reservoir water-level.

water-level of reservoir, FLAC3D, tilting field, displacement

1009-6825(2016)31-0221-02

2016-08-24 ★:山西省地震局科研項目(項目編號：SBK-1509)；2016年度震情跟蹤青年課題(項目編號：2016010219)

李宏偉(1986- )，男，碩士，工程師； 王 霞(1987- )，女，碩士，工程師； 李文超(1990- )，男，助理工程師

TV697.25

A