干施工圍堰滲流穩(wěn)定數(shù)學模型試驗研究

黃靈武 楊宗宇 王 衎

(國核電力規(guī)劃設計研究院，北京 100095)

?

干施工圍堰滲流穩(wěn)定數(shù)學模型試驗研究

黃靈武 楊宗宇 王 衎

(國核電力規(guī)劃設計研究院，北京 100095)

以某電廠取水工程為例，對干施工圍堰滲流穩(wěn)定性進行了有限元數(shù)值仿真模擬計算，論證了各工況下圍堤典型斷面的滲流穩(wěn)定性，研究了帷幕深度對滲流控制效果的影響，得出了一些有價值的結論。

圍堰，滲流場，防滲墻，穩(wěn)定性

1 工程概況

某電廠擬采用海水直流循環(huán)供水系統(tǒng)，需要在廠區(qū)建設循環(huán)水取、排水設施，截滲圍護設施由部分防波堤、明渠施工圍堰、廠區(qū)護岸內截滲結構及天然巖盤順次相接，共同組成一個封閉空間。為滿足干施工條件，需設干施工圍堰止水。干施工圍堰采用拋石斜坡堤結構，斜坡堤內設柔性地連墻止水，并對底部巖基進行帷幕灌漿處理。針對截滲結構的重要性，在保證安全、經濟、可行的前提下，需選擇合適的截滲結構方案。

干施工圍堰滲流穩(wěn)定數(shù)學模型試驗的目的就是要根據(jù)其結構特點，采用二維滲流模擬計算分析，在保證結構滲流穩(wěn)定的前提下，開展取水工程滲流數(shù)值模擬計算，分析研究取水工程滲流穩(wěn)定性和防滲措施的可靠性，確定截滲結構的合理入巖深度及帷幕灌漿深度，對截流結構方案提出優(yōu)化建議，以滿足設計需要。

2 基本理論

采用二維非穩(wěn)定滲流數(shù)學模型：

對于符合達西定律的二向均質各向異性的土體的滲流，其水頭函數(shù)為以下偏微分方程式：

(1)

與之相應的定解條件為：

初始條件：

ht=0=h(x，z，0)

(2)

邊界條件：水頭邊界：

ht=t1=h(x，z，t)

(3)

流量邊界：

(4)

其中,z為坐標位置高程；q為自由面下降或上升時從自由邊界流入或流出滲流場的單寬流量。

式(1)加上相應的初始條件和邊界條件式(2)，式(3)就是描述地下水滲流的二維數(shù)學模型。采用有限單元法求解式(1)及相應定解條件，其原理和實現(xiàn)方法具體參閱文獻[1]。

（6）Who Owns Your Dinner?A Discussion of America's Patented Genetically Engineered Food Sources,and Why Reform is Necessary（誰將擁有你的餐桌？論美國基因工程食品專利及其改革的必要性）。

3 計算與分析

選擇防波堤典型堤段來進行不同工況下的二維穩(wěn)定和非穩(wěn)定滲流有限元法數(shù)值模擬計算，分析不同工況組合下滲流場和滲流特征值(滲流量、滲流水位、滲透坡降)，判別各工況下滲流穩(wěn)定性，并進行柔性防滲墻和帷幕灌漿的模擬計算，論證防滲措施的可靠性，提出合理可行的優(yōu)化防滲措施。

3.1 計算參數(shù)

1)水位。50年一遇高水位：1.69 m。

流場分析用平均潮位水位：1.23 m，滲流岸坡穩(wěn)定計算水位：1.69 m。

2)波浪要素：考慮25年一遇波浪。

3)堤身及各土層的滲透系數(shù)或透水率，見表1。

3.2 模型布置

有限元網格布置見圖1。

整個模型沿防波堤軸線垂向剖分30個斷面，每個斷面上按實際土層高程進行各土層的精確模擬，總計剖分結點26 244個，計算單元采用空間四面體單元，可適應復雜多變的地下結構與土層分布。單元布置原則是，滲流急變區(qū)(如防滲墻、帷幕等)附近單元密些，其他區(qū)相對疏些。

表1 滲流模型計算參數(shù)

3.3 計算工況

取水工程渠底開挖到-6.00 m時與外側潮位的滲流水頭差最大，因此考慮施工竣工期為研究典型工況?？紤]水位組合工況包括：1)施工竣工期外側設計最高潮水位下防護堤及取水渠穩(wěn)定滲流場計算；2)施工竣工期外側設計潮水位變化過程和設計波浪過程防護堤及取水渠非穩(wěn)定滲流場計算。

考慮的防滲措施工況組合包括：

1)無防滲措施；

2)采用柔性防滲墻。

3.4 計算成果及分析

防滲墻方案滲流場計算。對防滲墻成墻后進行了入巖深度和滲透性2種工況的數(shù)值模擬及敏感性計算分析，結果如下：

1)滲透量分析。

表2 典型斷面防滲墻

從表2計算結果來看，該斷面堤身由開山石構成，滲透系數(shù)較大，防滲墻入巖深度的變化對下滲的滲透量影響比較明顯，而防滲墻滲透系數(shù)大小對滲漏量具有一定影響，但影響滲漏量的程度不如入巖深度敏感。因此選用防滲墻的幾何參數(shù)應優(yōu)先考慮幾何尺寸，再考慮滲透性能。從防滲措施的施工工藝和安全角度，宜選擇防滲墻進入強風化層深度1 m～1.5 m，防滲墻滲透系數(shù)在10×10-7cm/s 以下可以滿足設計要求。

2)滲流穩(wěn)定性分析。

從表3可知，防滲墻滲透系數(shù)越低，防滲墻最大坡降越大，1號斷面防滲墻坡降在8～10，2號斷面防滲墻坡降在7～8，不同

滲透性的防滲墻均起到較好的防滲作用。

表3 上游水位1.23 m防滲墻不同滲透系數(shù)值典型斷面計算結果

表4為防滲墻不同入巖深度的滲流計算結果。從計算的11個斷面中可以看出，防滲墻坡降變化范圍比較大。最小防滲墻坡降為2號斷面，為3.8,該斷面堤身由原有土層構成，且滲透性較小，該斷面防滲墻防滲效果相對不明顯。而對于堤身由開山石筑成的斷面，防滲墻防滲效果較為明顯，防滲墻最大坡降在7～12。對于同一個斷面，防滲墻進入強風化層深度對防滲效果敏感度也不一樣，對于1號、6號、9號斷面防滲墻坡降隨入巖深度變化較為敏感，入巖深度每增加0.5 m，坡降變化范圍在0.2～0.9內。

表4 上游水位1.23 m防滲墻入巖不同深度計算結果

4 結論及建議

綜合有限單元法滲流模擬計算分析成果，有如下結論：

1)采用二維地下水滲流數(shù)值模型進行了模擬計算，結果表明：無防滲的干法施工方案滲漏總量較大，無法保證干法施工的有效實施，因此必須進行防滲處理，采用防滲墻截斷透水層結合基巖帷幕灌漿方案截斷基巖滲透的方案是合適的。

2)截滲圍護設施中柔性地下防滲墻防滲作用較為明顯，其下切高程必須截斷透水層，因此地下連續(xù)墻下切底高程需在強風化層內1 m～1.5 m，此時滲流量和出逸坡降均較小。

3)防滲墻滲透系數(shù)大小對滲漏量具有一定影響，但影響滲漏量的程度不如入巖深度敏感。因此選用防滲墻的幾何參數(shù)應優(yōu)先考慮幾何尺寸，再考慮滲透性能,施工時需根據(jù)實際地質情況調整防滲墻深度。

[1] Wibbeler,H.,Meissner,U..FEM Model for seepage flow through rubble mound breakwaters[J].Mathematical Modeling in Water Resources,1992(10):609.

Numerical simulation research on seepage stability of cofferdam under dry construction

Huang Lingwu Yang Zongyu Wang Kan

(StateNuclearElectricPowerPlanningDesignandResearchInstitute,Beijing100095,China)

Taking the water intake engineering of the power plant as an example, the paper carries out finite element simulation calculation of the dry-construction cofferdam seepage stability, discusses seepage stability of different cofferdam sections under various working conditions, studies the impact of curtain depth upon seepage control effect, and finally draws some valuable conclusions.

cofferdam, seepage field, diaphragm wall, stability

1009-6825(2016)28-0090-02

2016-07-28

黃靈武(1984- )，男，工程師

TU413.62

A