吳俊杰，潘旭東，李少杰,李坤儒

(1.新疆水利水電勘測設計研究院有限責任公司,新疆 烏魯木齊 830000；2.新疆新華葉爾羌河流域水利水電開發(fā)有限公司,新疆 喀什 844000)

新疆阿爾塔什水利樞紐工程中的聯(lián)合進水口塔群在距壩軸線約1.03km上游右岸南—北向的溝內(nèi)，進水塔部位地形北高南低，基巖裸露。受右岸山體地形及高陡邊坡的限制，塔群布置場地較小且無法進行大開挖，因此，1—3號進水塔最終采用“緊湊型”布置，即1號、2號發(fā)電引水洞進水塔布置兩側(cè)，3號深孔放空排沙洞進水塔布置于中間。聯(lián)合進水口塔群段處于微風化巖體上，岸坡走向340°，自然邊坡為45°～65°，岸坡基巖裸露，層理發(fā)育，巖層產(chǎn)狀340°SW∠58°，傾向坡外，巖層走向與岸坡走向平行，在聯(lián)合進水口塔群后邊坡和基礎開始施工過程中，塔群背后高邊坡出現(xiàn)2～3小規(guī)模順層面滑坡。雖然后邊坡設置了一些100t級預應力錨索和錨梁?？紤]上述原因還是擔心后期該處巖層還會順層面滑坡的可能性，如圖1所示。

圖1 阿爾塔什水利樞紐工程聯(lián)合進水口塔群后邊坡實際情況

因此，要求將聯(lián)合進口進水閘塔塔群整體前移10m，受地形條件限制塔群右側(cè)的1號進水塔右邊墻回填混凝土高程有所下降且范圍較大。根據(jù)吳俊杰[1- 2]在2015年技施階段計算的阿爾塔什聯(lián)合進水口塔群靜動力計算分析成果，以及期間進行的塔背回填混凝土最優(yōu)高度敏感性分析成果可知，塔身應力在地震工況下受回填混凝土高度影響很大，該結(jié)論也在崔煒等[3]、楊光等[4]、李鋒[5]、徐東芝等[6]、苗宏恩等[7]、曹偉等[8]、孔科等[9]、張紫璇等[10]、唐青山等[11]的研究成果中被證實過。同時，根據(jù)趙曉紅等[12]對紫坪鋪水電站的2號進水塔進行了抗震復核，結(jié)果表明在原設計7度和遭遇9度以上地震時，位于回填混凝土一定高程高程范圍內(nèi)的進水塔邊墻結(jié)構(gòu)應力均超出允許范圍，證明該處在地震工況時最薄弱。

可見進水塔塔體的穩(wěn)定性受回填混凝土高度的影響很大，這勢必會影響1號、2號發(fā)電引水系統(tǒng)進水口進水塔和3號塔的受力狀態(tài)，所以亟需進行三維有限元靜、動力分析并與原設計方案的塔群應力進行比較，得出前移10m方案后塔體應力的變化范圍，自振特性、應力與變形影響分析，在塔群尚未澆筑前為前移方案提供依據(jù)。

1 概況

新疆阿爾塔什水利樞紐工程位于喀什地區(qū)莎車縣，工程在建設時面臨著“三高一深”等諸多設計難點，即最大壩高164.8m，600m級高邊坡，Ⅷ度高地震烈度，94.5m深厚覆蓋層。該樞紐工程大體右岸是本次研究對象1號、2號進水塔是發(fā)電引水洞前段控制段和3號塔深孔放空排沙洞控制段，其中1號、2號進水塔塔頂高程為1827.00m，底板高程為1750m，進水塔高度77m，底板厚4m。3號塔設置在1號、2號進水塔中間。

2017年現(xiàn)場聯(lián)合進水口塔群后邊坡、基礎開挖施工階段，前后經(jīng)過2～3次小規(guī)模的順層滑坡，擔心后期該處巖層還會有順層面滑坡的可能，為后期安全施工與運行，經(jīng)過論證后將聯(lián)合進口進水口塔群整體前移10m就可避免順層開挖，但是，受地形條件限制塔群右側(cè)的1號進水塔右邊墻回填混凝土高程會有所下降且范圍較大。右側(cè)1號進水塔后胸墻右邊墻回填混凝土從原來1795m降至1785m，右邊墻回填混凝土從原有的1785m以3級臺階形式逐級降至1770、1765、1760m，以下是原設計方案與前移10m方案對比平面圖，如圖2所示。

圖2 原方案與前移方案對比平面圖

2 動力計算

結(jié)構(gòu)地震反應分析的方法一般有振型疊加反應譜法、時程分析法[13- 15]。其中振型分解反應譜法是目前結(jié)構(gòu)抗震設計規(guī)范規(guī)定的一種常規(guī)方法，本次與原設計方案計算方法一致，按規(guī)范要求采用振型疊加反應譜法對塔群結(jié)構(gòu)進行動力分析，該方法被廣泛應用。動力計算振動方程如下式[16- 18]。

(1)

阻尼采用瑞利阻尼，假定[C]阻尼矩陣，[K]為剛度矩陣，[M]為質(zhì)量矩陣的線性組合，即滿足下列公式[1- 2]：

[C]=α[M]+b[K]

(2)

式中，α、b—比例常數(shù)。

3 動力計算反應譜

壩址區(qū)基巖地震基本烈度為Ⅷ度，Ⅰ類場地，基巖峰值加速度為221.0g。發(fā)電引水建筑物為1級非壅水建筑物，根據(jù)抗震設計規(guī)范要求，發(fā)電引水建筑物抗震設防類別為乙類，按50年超越概率5%進行設防。特征周期按照規(guī)范要求以及地震安評報告給出Tg取0.2s，出于安全考慮取0.3s，按照水工建筑物抗震設計標準規(guī)定，設計反應譜因采用標準加速度反應譜根據(jù)特征參數(shù)擬合出本工程動力計算加速度反應譜如圖3所示。

圖3 場地基巖設計反應譜

4 三維有限元計算分析

通過大型有限元軟件ABAQUS[19- 23]建立本次前移10m方案的實體模型如圖4(b)所示，為了便于比較原設計方案也列于此如圖4(a)所示。由于時間較緊，因此1—3號進水塔全部采用C3D20四面體20結(jié)點單元，整體模型結(jié)點數(shù)量共計102686個，單元101218個。整體實體模型與網(wǎng)格剖分圖如圖5所示。

圖4 計算實體模型

圖5 聯(lián)合進水口塔群

4.1 計算工況

本次計算所有參數(shù)均有原方案三維有限元靜動力分析及配筋設計[1]相同，同時，為了能夠真實反映塔體與回填混凝土、塔體與基巖基礎之間的相互關系，計算利用ABAQUS外掛子程序?qū)崿F(xiàn)接觸之間的靜、動力相互作用，然后對進水塔群-閘房-整體地基-回填混凝土-塔體內(nèi)外水體進行靜、動力聯(lián)動耦合計算(與原方案[1]一致)。計算塔體結(jié)構(gòu)在前移10m方案后，完建期、正常運行期、正常運行期+地震、檢修期等工況下的應力應變情況。

在各個工況計算過程中，塔體荷載主要有水體重量，進水塔自重(含金結(jié)、水體壓力、浪壓力、底部揚壓力、地震力、動水附加質(zhì)量力)，按荷載規(guī)范要求，在不同工況時荷載需要進行組合見表1。

表1 荷載組合表

4.2 材料計算參數(shù)

工程區(qū)聯(lián)合進水口塔群在弱風化至微風化基巖基礎上，根據(jù)地質(zhì)提供地基承載力試驗可知，文獻[1]也有所交代，基礎承載力標準值在2.0～2.5MPa，具體數(shù)據(jù)見表2。

表2 材料的力學參數(shù)

4.3 計算成果分析

4.3.1靜力成果分析

通過計算可知聯(lián)合進水塔前移10m后，1號進水塔在完建期、運行期、檢修期邊墻外側(cè)回填混凝土高度降低后側(cè)向的摩擦力有減少，導致豎向荷載增大，從而使得地基基礎沉降增加，前移10m方案在靜力工況下對1號進水塔有影響相對原設計方案較小，2號進水塔和3號塔受力狀態(tài)基本沒有。

4.3.2動力成果分析

(1)前移10m方案與原方案結(jié)構(gòu)整體自振特性比較

兩組方案前10階振型頻率及周期見表3。

表3 進水口塔群自振特性

當1號進水塔邊墻回填混凝土高度下降后，整個聯(lián)合進水口塔群的邊界條件有所改變，導致1號進水塔、2號進水塔、3號塔(附帶各自閘房及金屬結(jié)構(gòu)質(zhì)量)的震動頻率有所降低，自振周期均有所提高，但幅度不大，這也表明建筑物回填混凝土比值越小，建筑物自振周期越長的特性。

(2)前移10m方案與原方案靜動疊加成果比較

在不改變外部邊界條件的情況下，原方案邊墻為1785m時，由于1號進水塔剛度比2號閘井剛度小，1號進水塔橫向+豎向地震中右震位移量最大為14.12mm如圖6(a)所示。前移10m方案1號進水塔邊右側(cè)邊墻降低后，右震位移量最大值15.58mm如圖6(b)所示，從增加的變形量也可說明，前移10m后塔體邊墻降低對1號進水塔橫向位移很大，水平位移增大了9.37%，從而導致邊墻變形增大。

圖6 運行期遇地震聯(lián)合進水口塔群變形矢量圖(單位:m)

前移10m方案在地震工況時邊墻應力如圖8(b)所示較原方案如圖7(a)所示有所增大，表明前移10方案1號進水塔右側(cè)地形變低后確實影響到了1號塔體應力，而2號塔體、3號塔體回填混凝土高度沒有改變，邊墻應力基本沒有改變，限于篇幅限制這里不在論述。

圖7 1號進水塔地震工況拉應力大于1.43MPa部分云圖(單位：Pa)

原設計1號進水塔地震工況計算成果見[1]，現(xiàn)將前移10方案1號進水塔在地震工況下的結(jié)構(gòu)應力(去除集中應力點后)表格的形式列于下表，同時進行配筋見表4。

表4 1#進水口塔群動力計算拉應力表

由上表計算成果可以得出前移10方案后，邊墻在受橫向地震作用后，1號進水塔左右邊墻配筋量均有所增加，左邊墻1780～1800m段配筋最大提高了29.88%，右邊墻1780～1800m段配筋最大提高了43.48%，并且左右邊墻配筋有下移加深趨勢，從配筋成果來看左右邊坡7000～8000mm2的配筋量已經(jīng)不合理，因此，該方案在地震工況時對1號塔配筋影響比較大。

5 結(jié)語

(1)由于本工程聯(lián)合進水口塔群在施工時遭遇2～3次小規(guī)模的滑坡，各方出于安全考慮提出將進水塔前移10m方案，以便后期施工安全。然后，受本工程地形條件限制，前移10m方案會造成1號進水塔右側(cè)限制高度降低，靜力工況塔體受力影響較小，而地震工況時右岸1號進水塔頂部最大水平位移增加了9.37%，同時自振周期也有所增加，導致邊墻配筋量較原方案有所增加同時配筋深度有所加深。對2號進水塔和3號塔底板稍有影響，主要是1號進水塔豎向荷載變大導致地基沉降量增大，從而影響到了2號進水塔和3號塔。對2號進水塔和3號塔邊墻影響很小，主要因為回填混您土沒有變化。最終，因前移10方案導致1號進水塔在地震工況時配筋過大且不合理，可能會影響1號塔體整體受力而否定該方案，現(xiàn)場最終對聯(lián)合進水口塔群后邊坡進行加強支護，增加多排預應力錨索與錨梁，經(jīng)復核邊坡整體穩(wěn)定性，為后續(xù)施工與安全運行提供保證。

(2)2021年8月工程全部機組并網(wǎng)發(fā)電進行涉網(wǎng)調(diào)試，到目前為止塔體后邊坡變形監(jiān)測儀器數(shù)據(jù)、塔群運行狀態(tài)一切正常，通過對比計算也證明了塔體外側(cè)約束高度對結(jié)構(gòu)在地震工況時的重要性，也為塔體結(jié)構(gòu)受力優(yōu)化提供借鑒。