李冬冬 陳捷平 吳秋波 孫海清

收稿日期：2023-04-23；接受日期：2023-06-26

基金項目：國家自然科學基金項目（52079097）

作者簡介：李冬冬，女，高級工程師，博士，主要從事地下洞室工程研究。E-mail：lidongdong3@cjwsjy.com.cn

Editorial Office of Yangtze River. This is an open access article under the CC BY-NC-ND 4.0 license.

文章編號：1001-4179（2024） 06-0194-10

引用本文：李冬冬，陳捷平，吳秋波，等.旭龍電站尾水調壓室洞間巖柱圍巖宏細觀特性分析

［J］.人民長江，2024，55（6）：194-203.

摘要：為了將連續(xù)和非連續(xù)分析方法耦合用于地下廠房圍巖穩(wěn)定性分析，提出了一種基于軟硬混合平行黏結模型的細觀參數(shù)快速標定法，可根據(jù)給定巖石的宏觀力學參數(shù)快速給出其顆粒流模型細觀接觸參數(shù)；在此基礎上，構建了基于顆粒離散元PFC 3D與有限差分FLAC 3D耦合的旭龍地下電站尾水調壓室三維離散-連續(xù)耦合分析模型，研究了開挖期尾水調壓室洞間巖柱圍巖變形與破壞宏細觀特性，以及檢修期間側向水壓力對洞間巖柱穩(wěn)定性的影響。結果表明：檢修期間尾水調壓室水壓力引起的圍巖變形與應力增量變化是極小的，但洞間巖柱上半段22.7 m的高度存在塑形貫通區(qū)，局部穩(wěn)定性較差，有必要在工程建設期開展尾水調壓室及尾水支洞巖墩塑性貫通區(qū)專項支護方案研究，以確保旭龍地下電站工程的安全建設與運行。

關? 鍵? 詞：尾水調壓室； 洞間巖柱； FLAC 3D； 離散-連續(xù)耦合； 塑形貫通區(qū)

中圖法分類號： TV731

文獻標志碼： A

DOI：10.16232/j.cnki.1001-4179.2024.06.027

0? 引 言

水電站地下廠房正朝裝機規(guī)模大、廠房跨度大、高邊墻的超大型化方向發(fā)展。擬建在建的一大批水電站采用單機容量70萬kW級別的機組，地下廠房跨度多達20余m、甚至超過30 m，高度多達50余m、最大已接近90 m［1-2］。主廠房、主變洞、尾水調壓室等引水發(fā)電系統(tǒng)各洞室縱橫交錯布置，形成龐大、復雜的大規(guī)模密集地下洞室群［3-4］?？紤]工程地質條件、主要樞紐建筑物布置、工程造價等因素，地下洞室群常采用集約布置理念，各大洞室間巖柱間距較小，可能影響巖柱與洞室群的局部或整體穩(wěn)定性［5-6］。

實際工程中，洞室間距對圍巖穩(wěn)定性、支護設計和施工方案等具有重要影響［7-10］。以旭龍地下電站為例，尾水調壓室及尾水洞具有開挖體型復雜、岔管較多、山體挖空率高的特點，尾水調壓室洞間巖體厚度為16.0 m，高54.2 m，在開挖工況下存在巖柱破壞區(qū)貫穿的風險，危及圍巖穩(wěn)定性。檢修工況下同一水力單元的兩臺機組檢修、另一水力單元正常運行，調壓室洞間巖柱將承擔水壓力，需要明確該工況時洞間巖柱的穩(wěn)定性。

目前對于地下洞室群洞間巖柱的分析研究，常采用工程類比法和數(shù)值分析法。前者主要基于國內外已建地下廠房洞室群的洞室布置參數(shù)進行類比分析［11-12］，后者隨著計算機技術發(fā)展已經(jīng)被廣泛應用到地下巖石工程領域的研究中?？焖倮窭嗜沼邢薏罘殖绦騀LAC 3D和顆粒離散元程序PFC 3D均為美國ITASCA公司研發(fā)的大型仿真計算軟件，前者屬于連續(xù)介質力學法的范疇，可以模擬各類巖土體和結構單元（錨桿、錨索、襯砌等）的變形與受力狀態(tài)，得到了國內外的廣泛應用和認可，積累了豐富的工程經(jīng)驗［13～15］；后者則是從非連續(xù)介質力學法（離散元）的角度研究介質細觀力學特性和行為的工具，目前鮮少用于大型地下洞室群圍巖穩(wěn)定分析與研究。二者的程序可在計算過程中相互傳輸變量并擁有共同的二次開發(fā)語言FISH［16-18］。

本文集連續(xù)-非連續(xù)分析方法之所長，基于PFC 3D-FLAC 3D離散-連續(xù)耦合分析方法，構建旭龍地下電站尾水調壓室三維離散-連續(xù)分析模型，重點分析開挖期和檢修期受側向水壓力時，尾水調壓室洞間巖柱的變形受力特征、破壞區(qū)分布、微裂紋發(fā)育規(guī)律等宏-細觀特性，以期為施工期尾水調壓室及尾水支洞巖墩塑性貫通區(qū)專項支護方案提供參考。

1? 旭龍尾水調壓室簡介

旭龍水電站位于云南省德欽縣與四川省得榮縣交界的金沙江干流上游河段，是金沙江上游河段“一庫十三級”梯級開發(fā)方案中的第12級，是西電東送骨干電源點之一。地下廠房集中布置在右岸，安裝4臺單機容量為600 MW的混流式水輪發(fā)電機組，裝機容量為2 400 MW。主廠房、主變室和尾水調壓室等三大洞室平行布置（圖1），主廠房與主變洞洞間巖體厚度為40.0 m，主變洞與尾水調壓室洞間巖體厚度為30.5 m。

尾水調壓室為阻抗式、長廊型，調壓室前尾水隧洞為一機一洞，調壓室后為兩機一洞布置，尾水調壓室開挖尺寸為161.5 m×19.1 m（20.0 m）×87.6 m。尾水調壓室外側端墻距岸邊240.0 m，洞室埋深290.0～410.0 m，調壓室軸線方向NE55°。圍巖以灰白色中細?；◢弾r為主，近岸側為灰黑色斜長角閃巖脈體（圖2）；花崗巖巖質新鮮、堅硬，構造不發(fā)育，巖體較完整—完整，單軸飽和抗壓強度90～110 MPa，巖體聲波測試Vp均值為5.29 km/s；斜長角閃巖脈體長英質含量較高，巖質較堅硬—堅硬狀，巖體完整性差—較完整，單軸飽和抗壓強度80～90 MPa，巖體聲波測試Vp均值為4.67 km/s。

根據(jù)尾水調壓室下游54.0 m處地應力測試孔的試驗成果，下游范圍內最大水平主應力范圍主要為4.1～32.5 MPa，最小水平主應力范圍主要為2.8～18.1 MPa，最大水平主應力為中等—高地應力水平。最大水平主應力方向為N27°～47°E（平均N38°E），與洞室軸線夾角8°～28°。

尾水調壓室圍巖花崗巖為塊狀構造，斜長角閃巖脈體以次塊狀構造為主，巖性分界面走向355°～5°，傾向W，傾角70°～80°，巖性界線與調壓室軸線夾角50°～60°。地表f3斷層斜穿尾水調壓室中部，斷層走向30°，傾向300°，傾角75°～80°，斷層走向與尾水調壓室軸線夾角25°，斷層帶寬0.2～0.3 m，構造巖主要為碎裂巖，原巖為花崗巖或斜長角閃巖，較破碎，膠結一般；主斷面為厚0.1～2.0 cm斷層泥，局部為鈣質及巖屑。尾水調壓室洞室圍巖Ⅱ類圍巖長度為125.0 m，Ⅲ類圍巖洞室長度為25.0 m，Ⅱ類圍巖占83%，Ⅲ類圍巖占17%。

2? 尾調室離散-連續(xù)耦合模型

2.1? PFC 3D-FLAC 3D耦合分析原理

PFC 3D-FLAC 3D離散-連續(xù)耦合分析是指將FLAC 3D連續(xù)區(qū)域與PFC離散區(qū)域的接觸面指定為PFC的墻體，采用等效力方法將顆粒作用于墻上的接觸力和接觸彎矩分配到墻面的頂點上，墻體頂點則附著于實體單元的網(wǎng)絡點或結構單元的節(jié)點上，由此實現(xiàn)了墻體頂點運動與實體單元網(wǎng)格點（或shell單元節(jié)點）的同步運動。PFC墻體傳遞來的力和彎矩，參與連續(xù)域的FLAC 3D模型計算分析進程，連續(xù)域節(jié)點的變形也帶動了墻體的運動，進而將位移和速度通過墻傳遞到離散域中的顆粒。通過上述各項值的不斷更新，改變模型幾何參數(shù)以及單元或結構單元的剛度。

假設一個球與一個三角形墻面接觸，C為接觸點，CP為墻面上距離C點最近的點，從CP到三角形的頂點采用重心插值法外推。三角形頂點為Vik，（i=1，2，3，k=1，2，3），通過鏈接3個頂點與CP，得到3個三角形（圖3），面積分別為Ai，3個三角形的總面積為A=A1+A2+A3。

頂點加權因子定義為采取與頂點相對三角形面積除以三角形的總面積，ωi=Ai/A，等式ωi=1恒成立，以保證從CP到頂點外推值總和等于CP處的值。當完整計算模式關閉時，接觸力和平動剛度也通過這種方式外推到節(jié)點處，但由于頂點力和潛在的頂點力矩可能會不平衡，重心外推不能確保外推力和力矩與CP位置的瞬時接觸力和力矩一致。

假定Ri（i=1，2，3）為CP到3個頂點的向量，Ri=Vik－CP，施加在每一個網(wǎng)格點或節(jié)點的力為Fi（i=1，2，3），施加在接觸點C上的力為F，接觸點處由于黏結產(chǎn)生的力矩為Mb。由于接觸點C和接觸平面上點CP可能不共點，所以作用在接觸平面上的總力矩為

M=Mb+C－CP×F（1）

當全運算模式處于活動狀態(tài)時，耦合方案確定完全一致的等效力系統(tǒng)，即：3i=1Fi=F（2）

Ri×Fi=M（3）

n為三角形的單元法向向量，沿著三角面的剪切力矢量為Fs，即：Fs=F－F·n（4）

切向單位矢量為

S=FsFs（5）

局部坐標系中x軸的方向與法向n的方向一致，y軸的方向與剪切方向s方向一致，由于CP點在三角面上，Ri在z方向上均為0，即Ri，z=0。這種簡化可以直接確定局部坐標系統(tǒng)中頂點力和力矩的x、y、z方向上分量Fi，x、Fi，y、Fi，z，即：Fi，x=Fx

Fi，y=Fy

Fi，z=Fz=0（6）

Ri，y×Fi，z－Ri，z×Fi，y=Mx

Ri，z×Fi，x－Ri，x×Fi，z=My

Ri，x×Fi，y－Ri，y×Fi，x=Mz（7）

由于之前的重點加權項使得力的y分量分布Fi，y=ωiFy，因此在三角形平面中最大接觸力的方向上施加重點加權。這個假設的結果導致以下2個方程有3個未知數(shù)：Fi，z=0（8）

Ri，y×Fi，z=Mx+Ri，z×Fi，y（9）

對于這種欠正定方程組，在無約束條件下有無窮多解。控制方程式不能明確地確定約束的性質來找到一個等效的系統(tǒng)。因此為獲得一個特解，需要提供額外的約束條件，緊隨著實體單位或結構單元運動方程，以根據(jù)相應的網(wǎng)格點節(jié)點位置和速度設置墻面頂點的位置和速度。

Ri，y×Fi，z=0（10）

式（10）可解釋為：如果開始時在局部z方向上的每個網(wǎng)格點或節(jié)點施加力，則從CP指向局部坐標系中的三角形頂點向量的點積總和為0。所得到的合力隨后被轉換為全局坐標系統(tǒng)并應用到適當?shù)膮^(qū)域網(wǎng)格點或結構單元節(jié)點。此時觸發(fā)區(qū)域和結構單元更新。在區(qū)域和結構元件更新期間，對網(wǎng)格節(jié)點、結構單元節(jié)點添加剛度。當接觸處存在旋轉剛度時，這些也必須考慮到以確保穩(wěn)定性。通過將旋轉剛度的大小除以Ri來確定等效的平移剛度。

2.2? 尾調室整體模型與耦合模型

通過FLAC 3D建立的尾水調壓室數(shù)值分析模型其范圍為527.3 m×450.1 m×524.5 m（X×Y×Z），其中調壓室開挖尺寸為161.5 m×19.1 m（20.0 m）×87.6 m（長×寬×高）。計算模型及網(wǎng)格劃分見圖4：劃分單元262 403個、節(jié)點147 642個，x軸平行于尾水調壓室縱軸線，指向山體內部為正；y軸順水流向，向下游為正；z軸鉛直向上，底部高程為2 101.0 m。模型包括尾水調壓室以及部分尾水支洞與尾水主洞。同時考慮了f3斷層及不同巖層分布產(chǎn)狀。尾水調壓室自上而下共分22層開挖，為了保護洞間巖柱的穩(wěn)定性，開挖到洞室中下部時先開挖岔管部分（第5期）、后掏空2.5 m寬的巖柱（第6期），再按照2 m每層的厚度逐層向下進行開挖。

分別在尾水調壓室洞間巖柱的頂部、中部和底板處設置耦合區(qū)域（圖5）。耦合區(qū)域1和2沿尾水調壓室縱軸線方向貫穿巖柱，長16.0 m，橫截面積為2.0 m×2.0 m，分別生成耦合墻體256個、顆粒22 145個，顆粒半徑為7.5～10.0 cm，主要用于研究洞間巖柱在開挖過程中的變形與破壞宏細觀特性。

耦合區(qū)域3設置在底板附近，主要用于研究檢修工況下洞間巖柱承擔單方向水壓力時，洞間巖柱邊墻中下部的圍巖宏細觀特性，耦合區(qū)域尺寸為4.0 m×3.0 m×4.0 m（X×Y×Z），生成耦合墻體192個、顆粒12 060個，顆粒半徑為7.9～11.5 cm。顆粒間為平行黏結接觸模型，顆粒與墻體之間為線性接觸模型。

2.3? 巖體宏-細觀力學參數(shù)快速標定

在室內、現(xiàn)場試驗成果的基礎上，結合具體工程地質條件，確定地下廠房區(qū)主要巖石（體）和結構面的物理性質和力學參數(shù)建議值見表1和表2。地下廠房洞室群區(qū)域地層巖性主要為花崗巖、混合巖。綜合原位試驗成果、洞室群開挖卸荷及可能存在的裂隙等不利地質構造等因素，圍巖穩(wěn)定分析時，各巖層和結構面所采用的物理力學參數(shù)建議取低值。

在采用PFC 3D-FLAC 3D離散-連續(xù)耦合分析方法研究洞室群巖體宏細觀力學特性之前，首先需要解決的、也是關鍵問題之一的是PFC顆粒流模型細觀力學參數(shù)的選取。實際巖體中都是由剛度顯著不同的顆粒構成的，其變形過程中不同顆粒之間可以協(xié)調變形，而采用均質的、相同黏結和剛度參數(shù)時，則無法反映這種現(xiàn)象。

為此，可以通過混合黏結模型參數(shù)來更真實地逼近實際情況。假定巖石試樣由軟硬兩種顆粒接觸構成，其中40%的接觸為軟接觸。軟接觸模量為硬接觸模量的0.1倍，黏結強度為硬接觸黏結強度的0.08倍。FISH語言實現(xiàn)方法為：遍歷模型所有接觸，對于任一接觸，定義一個介于0和1的隨機數(shù)，當隨機數(shù)的值小于0.4時，則賦予該接觸為軟接觸的參數(shù)。通過上述基于混合平行黏結模型的細觀參數(shù)快速標定法，確定旭龍地下廠房區(qū)巖體宏觀參數(shù)對應的細觀接觸參數(shù)取值。下面以廠區(qū)某花崗巖為例闡述顆粒流模型細觀參數(shù)快速標定的過程。

待標定的花崗巖宏觀試驗參數(shù)為：拉伸彈性模量40 GPa，壓縮彈性模量40 GPa，泊松比0.23，單軸抗壓強度90 MPa，單軸抗拉強度9 MPa，拉壓比0.1。

用于標定細觀參數(shù)的PFC顆粒流模型如下：寬2.0 m、高4.0 m，粒徑在1.0～2.0 cm間隨機分布，初始平衡得到模型后，修改參數(shù)并激活接觸，開展各種數(shù)值模型試驗，獲取試樣的宏觀力學參數(shù)，并不斷調整參數(shù)直至完成標定。

巖石試樣的宏觀力學試驗特性與顆粒流模型數(shù)值試驗結果的對應關系如表3所列。

（1） 將線性有效接觸模量emod保持為相對較小的數(shù)值（0.1 MPa），依次設置平行黏結有效模量pb_emod為60，80，100，120 GPa，分別獲取直接拉伸條件下的應力-應變曲線（圖6），試樣的宏觀拉伸彈性模量即為峰值應力與峰值應變的比值。擬合得到拉伸彈性模量與平行黏結模量pb_emod的關系，將待標定巖石試樣的拉彈性模量Et=40 GPa，代入趨勢線公式（圖7）Et=0.385 8x+ 2.172 2，求出平行黏結模量pb_emod的值為98.50 GPa。

（2） 固定平行黏結模量pb_emod的值為98.5 GPa，采用無側限壓縮試驗標定線性接觸模量，研究線性有效接觸模量emod對壓縮模量的貢獻。依次設置線性有效接觸模量emod為10，15，20，25，30，40 GPa，分別獲取單軸壓縮時試樣的應力-應變曲線（圖8），壓縮彈性模量即為峰值應力與峰值應變的比值。擬合得到壓縮模量與線性有效接觸模量emod的線性取值關系Ec=0.425 1x+27.723。根據(jù)待標定試樣的壓彈性模量Ec=40 GPa，求出線性有效接觸模量emod的值為23.52 GPa。

（3） 固定線性有效接觸模量emod為23.52 GPa和平行黏結模量pb_emod為98.50 GPa，假定平行黏結分量與線性接觸分量的剛度比kratio相同，研究剛度比與宏觀泊松比的對應關系。依次設置剛度比為1.5，2.0，2.5，3.0，3.5的情況進行單軸壓縮試驗，取名義應變（峰值強度一半對應的應變）所處位置計算泊松比（圖9）。擬合得到巖石試樣泊松比與剛度比kratio的對應關系μ=0.074 1x+0.024 7。將巖石實際泊松比μ=0.23代入，即可求得剛度比kratio的值為2.77。

（4） 當法向與切向黏結強度比值越大，試樣越容易出現(xiàn)剪切破壞；當比值越小，試樣越容易出現(xiàn)脆性破壞（圖10）。應根據(jù)實際試驗中的巖石破壞情況，近似在0.5～2之間取值。綜合數(shù)值模型試驗與試驗機結果，此處黏結強度比取為1.2，確定后固定該比值不變。

假定切向黏結強度的取值為10 MPa，根據(jù)法向與切向黏結強度比1.2，求得法向黏結強度為12 MPa，二者稱為基準黏結強度。在基準黏結強度的基礎之上，依次同時乘以黏結強度變化系數(shù)3，6，8，10，12.5，15，18，分別得到不同的單軸抗壓強度。擬合得到巖石試樣單軸抗壓強度σc與黏結強度變化系數(shù)間的對應關系為σc=6.936 4x+0.336 7。巖石試樣實際抗壓強度為90 MPa，求得放大系數(shù)為12.926，再乘以基準黏結強度10 MPa，即可得到法向黏結強度為15.512 MPa，切向黏結強度為12.926 MPa。

（5） 根據(jù)上述步驟得到了該花崗巖顆粒流模型細觀接觸參數(shù)如下：平行黏結模量pb_emod為98.50 GPa，線性有效接觸模量emod為23.52 GPa，剛度比kratio為2.77，法向黏結強度pb_ten為15.512 MPa，切向黏結強度pb_coh為12.926 MPa。其他參數(shù)：摩擦角pb_fa為45°，半徑乘子pb_rmul為1.0，彎矩貢獻系數(shù)pb_mcf為1.0，法向應力更新模式lin_mode為1，摩擦系數(shù)fric為0.5。

（6） 由于各參數(shù)間會相互交叉影響，按照上述參數(shù)試算，巖石試樣拉彈性模量為31.398 GPa，壓彈性模量為33.90 GPa，單軸抗拉強度為8.086 MPa，單軸抗壓強度為90.257 MPa，拉壓比為0.089 6。模量比待標定值略偏小，此時可以按照實際數(shù)據(jù)進行微調，提高平行黏結模量與線性有效接觸模量，并適當降低黏結強度（表4）。通過若干次試算，最終標定得到的參數(shù)（表5），已經(jīng)非常接近巖石真實的宏觀力學參數(shù)。

3? 開挖期洞間巖柱宏細觀特性

第6期開挖后洞室圍巖變形最大值22.0 mm出現(xiàn)在兩側邊墻，方向指向洞內。洞間巖柱頂端以豎直

向上的變形為主，第6期開挖后位移最大值9.5 mm出現(xiàn)在巖柱頂端轉角處，巖柱兩側位移較小。隨后在尾調室下半段開挖過程中，以每層2.0 m左右的高度向下開挖，巖柱頂端位移持續(xù)增大，巖柱兩側的變形隨高度降低而呈線性減小趨勢。到開挖完畢后，洞間巖柱頂端最大位移為14.6 mm，兩側位移矢量相對巖柱中心軸基本呈對稱分布（圖11）。

除洞間巖柱頂端臨近上部開挖臨空面的區(qū)域外，尾水調壓室頂拱、上部邊墻以及岔管的開挖未在洞間巖柱產(chǎn)生圍巖破壞區(qū)。下半段洞室邊墻開挖過程中洞間巖柱產(chǎn)生了破壞區(qū)，并在后續(xù)開挖時產(chǎn)生了塑性區(qū)貫通（圖12）。到開挖完畢時巖柱上半段22.7 m的高度存在塑形貫通區(qū)，下半段31.5 m高度范圍內塑性區(qū)深度呈線性逐漸降低，中間仍保留有一定范圍的彈性區(qū)，到巖柱底板處圍巖塑性區(qū)幾乎為0。本文中圍巖破壞區(qū)類型的釋義如下：shear為剪切破壞單元；tension為張拉破壞單元；n為now，指應力正位于屈服面；p為previous，指曾進入屈服狀態(tài)但現(xiàn)已退出；下同。

耦合區(qū)域1位于尾水調壓室洞間巖柱頂端，臨近上部洞室開挖臨空面，圍巖破壞區(qū)發(fā)展較快且深度較大，在第6期開挖后即出現(xiàn)了貫通塑性區(qū)，相應在PFC 3D耦合顆粒模型內部發(fā)育了分布均勻的微裂紋（圖13）；后續(xù)開挖過程中持續(xù)在耦合區(qū)域1內產(chǎn)生了一定程度的破壞，到開挖完畢時微裂紋單元數(shù)達到了6 190個。

耦合區(qū)域2位于尾水調壓室洞間巖柱中部，與尾水調壓室最后一期開挖的巖體相鄰，因此受上部洞室開挖影響相對較小，在尾調室第18期開挖以前，雖有少量微裂紋發(fā)育，但未在該區(qū)域內產(chǎn)生破壞單元。開挖完畢后，洞間巖柱兩側對稱產(chǎn)生了各5.3 m左右的塑性區(qū)，中間5.3 m為彈性區(qū)。PFC 3D耦合顆粒模型兩側的微裂紋集中發(fā)育、密度較大，中間微裂紋密度相對較?。▓D14），整個開挖過程共出現(xiàn)微裂紋單元2 536個。

4? 側向水壓力對巖柱穩(wěn)定影響

尾水調壓室檢修工況下同一水力單元的兩臺機組檢修，另一水力單元正常運行，調壓室間洞間巖柱將承擔水壓力。經(jīng)水力過渡過程計算，調壓室最高涌浪為2 170.35 m。計算出底板處最大水壓力為1 595 106 Pa，選取承受水壓力的尾調室?guī)r壁表面，利用梯度應力在該表面施加隨Z方向線性遞減的法向應力，來模擬呈三角形分布的水壓力（圖15）。

選取1-1′和2-2′作為典型剖面研究水壓力對洞室變形和受力的影響（圖16），其中2-2′剖面為理論上的水壓力集中作用面。

水壓力引起的調壓室典型剖面增量位移分布見圖17。尾調室一側施加水壓荷載后，隨著高程的降低，尾調室水壓引起的圍巖位移增量逐步增大。從山體外到山體內，隨著遠離荷載作用面，尾調室水壓引起的巖柱位移增量逐步減小。此外，尾調室?guī)r柱變形主要位于施加荷載一側的邊墻中下部，增量位移一般為0.1～0.7 mm，最大值為0.75 mm，位于高程2 126 m附近。

尾水調壓室水壓力作用前后，無水壓力時洞間巖柱最大主應力出現(xiàn)在山內側的巖柱與底板交匯處，局部拉應力達到0.195 MPa，出現(xiàn)在邊墻中上部（表6）。有水壓力時，洞柱巖體最大主應力分布云圖基本保持不變，最大壓應力降低0.2 MPa左右。承受水壓力一側的洞室圍巖，拉應力最大值（0.189 MPa）以及一般量值拉應力的范圍均有所降低。調壓室岔管部位圍巖最大壓應力大約降低0.2 MPa，最大拉應力降低0.01 MPa（圖18）。因此，總體來看，調壓室水壓力引起的圍巖應力狀態(tài)變化是極小的。

開挖完畢后，尾調室錨桿最大應力為203 MPa；調壓室水壓力作用前后，各特征點處錨桿應力變化極小。與開挖完畢時相比，水壓力作用下洞間巖柱的破壞區(qū)范圍基本沒有變化（圖19），頂端耦合區(qū)域1內微裂紋單元數(shù)增加11個，巖柱中部耦合區(qū)域2內微裂紋數(shù)量增加47個，說明調壓室水壓力引起的圍巖塑性區(qū)分布范圍及深度變化極小。

由圍巖宏觀受力特性可以看出，尾水調壓室洞間巖柱底板附近是巖柱最大主應力集中的部位，在此處設置耦合區(qū)域3以觀察底板圍巖在檢修工況下巖柱承擔單方向水壓力時的細觀變形與破壞特性。耦合區(qū)域3上半部分位于洞間巖柱內，下半部分位于巖柱底板以下圍巖內，左側臨近岔管開挖臨空面。在水壓力作用下洞間巖柱有以底部轉角為中心發(fā)生側向位移的趨勢，上半部分的側向位移明顯大于下半部分，位移最大值出現(xiàn)在施加荷載一側的左頂角，至右下角越來越小。說明洞間巖柱的最大側向位移不在耦合區(qū)域中間（即底板附近），結合連續(xù)單元的位移云圖（圖20）可以看出，巖柱1/3高度處側向位移最大。

開挖過程以及有無水壓力工況下的耦合區(qū)域3微裂紋分布形態(tài)見圖21?？梢钥闯?，底板圍巖內的細觀損傷主要出現(xiàn)在開挖期間，即第5期岔管開挖時，且靠近開挖面附近的微裂紋密度較大，遠離開挖面時微裂紋密度較小。水壓力作用下與開挖完畢時相比，微裂紋數(shù)量僅增加4個，說明水壓力幾乎沒有加深底板附近圍巖的破壞程度。

5? 結 論

本文提出了一種軟硬混合平行黏結模型的細觀參數(shù)快速標定法，其根據(jù)給定的巖石（體）力學參數(shù)，快速將宏觀-細觀參數(shù)對應起來，然后構建了基于顆粒離散元與有限差分耦合的旭龍地下電站尾水調壓室PFC 3D-FLAC 3D三維離散-連續(xù)耦合分析模型，用于開挖期和檢修期側向水壓力下尾調室洞間巖柱圍巖穩(wěn)定宏-細觀特性分析，主要結論如下：（1） 尾調室洞間巖柱頂端最大位移為14.6 mm，巖柱兩側的變形隨高度降低而呈線性減小趨勢，兩側位移矢量相對巖柱中心軸基本呈對稱分布。洞間巖柱上半段22.7 m的高度存在塑形貫通區(qū)，下半段31.5 m高度范圍內塑性區(qū)深度呈線性降低，中間仍保留有一定范圍的彈性區(qū)，到巖柱底板處圍巖塑性區(qū)幾乎為0。

（2） 尾水調壓室檢修工況下洞間巖柱承擔單方向水壓力時，最大位移0.75 mm位于施加荷載一側的邊墻中下部，從山體外到山體內隨著遠離荷載作用面，水壓引起的巖柱位移增量逐步減小。承受水壓力一側的洞室圍巖，最大壓應力降低0.2 MPa左右，最大拉應力降低0.006 MPa，因此調壓室水壓力引起的圍巖變形與應力增量變化是極小的。

（3） 研究結果表明，旭龍地下電站尾水調壓室洞間巖柱存在明顯的塑性貫通區(qū)，局部穩(wěn)定性較差，有必要在工程建設期開展尾水調壓室及尾水支洞巖墩塑性貫通區(qū)專項支護方案研究，確保旭龍地下電站工程的安全建設與運行。

參考文獻：［1］? 樊啟祥，林鵬，蔣樹，等.金沙江下游大型水電站巖石力學與工程綜述［J］.清華大學學報（自然科學版），2020，60（7）：537-556.

［2］? 樊啟祥，王義鋒，裴建良，等.大型水電工程建設巖石力學工程實踐［J］.人民長江，2018，49（16）：76-86.

［3］? 朱維申，李曉靜，郭彥雙，等.地下大型洞室群穩(wěn)定性的系統(tǒng)性研究［J］.巖石力學與工程學報，2004，23（10）：1689-1693.

［4］? 邵國建，卓家壽，章青.巖體穩(wěn)定性分析與評判準則研究［J］.巖石力學與工程學報，2003，22（5）：691-696.

［5］? 李術才，李樹忱，朱維申，等.三峽右岸地下電站廠房圍巖穩(wěn)定性斷裂損傷分析［J］.巖土力學，2000，21（3）：193-197.

［6］? 李寧，陳蘊生，陳方方，等.地下洞室圍巖穩(wěn)定性評判方法新探討［J］.巖石力學與工程學報，2006，25（9）：1941-1944.

［7］? 晏莉，陽軍生，張學民，等.水平互層巖體并行隧道中間巖柱穩(wěn)定分析［J］.巖石力學與工程學報，2009，28（增1）：2898-2904.

［8］? 郭佳奇，喬春生，曹茜.側部高壓富水溶腔與隧道間巖柱安全厚度的研究［J］.現(xiàn)代隧道技術，2010，47（6）：10-16.

［9］? 武愛玲，宋寶生.南水北調中線吳莊隧洞洞間巖體厚度的分析［J］.水科學與工程技術，2009（6）：56-57.

［10］王登科，駱建軍，高立平，等.基于爆破等效荷載的大型地下洞室群合理間距分析［J］.中南大學學報（自然科學版），2022，53（6）：2224-2233.

［11］段亞剛.地下大型儲備庫洞室斷面形狀優(yōu)化及合理間距研究［D］.北京：北京交通大學，2007.

［12］崔思嚴.地下廠房洞室群巖柱厚度、洞跨與洞高的統(tǒng)計相關分析［J］.四川水力發(fā)電，2006，25（5）：81-84.

［13］李仲奎，戴榮，姜逸明.FLAC3D分析中的初始應力場生成及在大型地下洞室群計算中的應用［J］.巖石力學與工程學報，2002，21（增2）：2387-2392.

［14］曾靜，盛謙，廖紅建，等.佛子嶺抽水蓄能水電站地下廠房施工開挖過程的FLAC3D數(shù)值模擬［J］.巖土力學，2006（4）：637-642.

［15］劉建華，朱維申，李術才，等.小浪底水利樞紐地下廠房巖體流變與穩(wěn)定性FLAC3D數(shù)值分析［J］.巖石力學與工程學報，2005，24（14）：2484-2489.

［16］李冬冬，胡浩然，肖明，等.基于PFC-FLAC離散-連續(xù)耦合的隧洞開挖與錨桿支護顆粒流模型［J］.武漢大學學報（工學版），2024，57（3）：286-294.

［17］李冬冬，胡浩然，肖明，等.全長粘結錨桿支護隧洞圍巖的顆粒離散元模擬與細觀作用機理研究［J］.武漢大學學報（工學版），2024，57（1）：37-44.

［18］李冬冬，盛謙，肖明，等.基于改進顆粒流聲發(fā)射片的地下廠房洞室圍巖局部損傷細觀機制研究［J］.巖土力學，2022，43（增2）：117-129.

（編輯：郭甜甜）

Macro and micro characteristics of surrounding rock of rock pillar between tail-race surge chamber caverns in Xulong Power Station

LI Dongdong，CHEN Jieping，WU Qiubo，SUN Haiqing

（Changjiang Survey，Planning，Design and Research Co.，Ltd.，Wuhan 430010，China）

Abstract：

To couple continuous and discontinuous analysis methods for the stability analysis of surrounding rock in underground power stations，a rapid calibration method for mesoscopic parameters based on a soft-hard mixed parallel bond model was implemented.This method can quickly correspond to the given rock macroscopic mechanical parameters to PFC model mesoscopic contact parameters.A three-dimensional discrete continuous coupling analysis model for the tail-race surge chamber of Xulong Power Station based on the coupling of particle discrete element PFC 3D and finite difference FLAC 3D was constructed.The macro-micro deformation and failure characteristics of rock pillars between tail-race surge chamber caverns during excavation were studied，as well as the impact of lateral water pressure on the rock pillar stability.The results showed that the deformation and stress increment of surrounding rock caused by the water pressure in the surge chamber under maintenance condition were minimal，but there was a plastic penetration zone at the height of 22.7 m in the upper half of the rock pillar between the caverns，with poor local stability.Thus it is necessary to carry out a special support scheme study on the plastic penetration zone of the tail-race surge chamber and tail-race branch tunnel rock piers during the construction period to ensure the safe construction and operation of the Xulong Power Station project.

Key words：

tail-race surge chamber; rock pillar between caverns; FLAC 3D; discrete-continuous coupling model; plastic penetration zone