曹凈++張丙軍++劉海明

摘要：隨著地鐵、越江隧道等地下工程的日益增多，涌現(xiàn)出了大量的對(duì)稱式基坑?；诠残斡成淅碚摚茖?dǎo)了對(duì)稱式基坑涌水量計(jì)算表達(dá)式，并對(duì)其計(jì)算參數(shù)進(jìn)行分析。研究表明，對(duì)稱式基坑的單寬涌水量與滲透系數(shù)、基坑內(nèi)外水頭差成正比，與（相對(duì)）隔水層層頂?shù)娇拥椎木嚯x和基坑兩側(cè)止水帷幕間距的比值、止水帷幕嵌固深度和（相對(duì)）隔水層層頂?shù)娇拥拙嚯x的比值、基坑內(nèi)外水頭差和（相對(duì)）隔水層層頂?shù)娇拥拙嚯x的比值成反相關(guān)。

關(guān)鍵詞：對(duì)稱式基坑；共形映射；滲流；涌水量

中圖分類號(hào)：TU46.3文獻(xiàn)標(biāo)志碼：A文章編號(hào)：16744764（2017）01010606

收稿日期：20160712

基金項(xiàng)目：國家自然科學(xué)基金（51304088）

作者簡介：曹凈（1966），男，博士，教授，主要從事巖土工程研究，（Email） cjlb117@sina.com。

Received：20160712

Foundation item：National Natural Science Foundation of China （ No.51304088）

Author brief：Cao Jing （1966），PhD， Professor， main research interest： environmental geotechnical engineering，（Email） cjlb117@sina.com.Calculation Method of seepage discharge for symmetric foundation pits

Cao Jing1，Zhang Bingjun1，2 ，Liu Haiming1

（1. Faculty of Civil Engineering and Mechanics， Kunming University of Science and Technology， Kunming 650500，

P. R. China； 2.Weifang Architectural Design & Research Institute Co.， LTD. Weifang 261041， Shandong， P. R. China）

Abstract：In recent years， with the development of crossingriver tunnels and subways， there are more and more symmetric foundation pits. Based on conformal mapping， the mathematical modeling is built and the expression of seepage discharge in pensile curtain foundation pit is deduced， then calculation parameters of seepage discharge are analyzed. The results showed that the seepage discharge of symmetric foundation pit has direct ratio with permeability coefficient and the water head between the outside and inside of foundation pit. And the seepage discharge of symmetric foundation pit is in inverse correlation to the ratio of the distance between the top of relative impermeable layer to the bottom of foundation pit and the distance between pensile curtains， to the ratio of the embedded depth of pensile curtain and the distance between the top of relative impermeable layer to the bottom of foundation pit， and to the ratio of water head between the outside and inside of foundation pit and the distance between the top of relative impermeable layer to the bottom of foundation pit.

Keywords：symmetric foundation pit； conformal mapping； seepage； discharge

隨著中國城市基礎(chǔ)建設(shè)的快速發(fā)展，地下空間的開發(fā)利用正成為各大城市的新熱點(diǎn)，隨之涌現(xiàn)出大量的深基坑。在基坑施工過程中常常因地下水的存在會(huì)引起一系列的問題，其中基坑的涌水量成為基坑降水設(shè)計(jì)和施工過程中所要重點(diǎn)考慮的問題。當(dāng)基坑工程地質(zhì)條件、水文地質(zhì)條件變化不大或?qū)挾容^小時(shí)，可近似地認(rèn)為基坑兩側(cè)采取的止水帷幕、基坑內(nèi)外地下水位、基坑滲流區(qū)土層厚度及滲透系數(shù)成縱向?qū)ΨQ，本文針對(duì)此類對(duì)稱式基坑進(jìn)行涌水量計(jì)算和滲流分析。

共形映射又稱保角映射，是復(fù)變函數(shù)的一個(gè)重要分支，廣泛應(yīng)用于流體力學(xué)[15]、復(fù)合材料斷裂[68]、隧道圍巖應(yīng)力變形[911]、電磁場(chǎng)理論[1214]。本文利用共形映射將有止水帷幕的對(duì)稱式基坑滲流平面轉(zhuǎn)化為帶狀區(qū)域，推導(dǎo)其單寬涌水量的計(jì)算表達(dá)式，并對(duì)對(duì)稱式基坑滲流單寬涌水量計(jì)算參數(shù)進(jìn)行分析。

1解析解推導(dǎo)

在城市基坑降水施工過程中，為保證基坑周邊已有建構(gòu)筑物、設(shè)施、管線的正常使用，《建筑基坑工程檢測(cè)技術(shù)規(guī)范》要求基坑坑外地下水位變化累計(jì)值不超過1 000 mm，變化速率不超過500 mm/d。故假定基坑坑外地下水位保持不變，如圖1所示。

圖1對(duì)稱式基坑滲流示意圖

Fig.1 schematic diagram of seepage field

in symmetric foundation pits圖1為對(duì)稱式基坑示意圖，為了簡化計(jì)算，假定滲流區(qū)土層為均質(zhì)、各向同性，地下水運(yùn)動(dòng)符合達(dá)西定律，基坑開挖深度為H，基坑坑外地下水位與地面齊平，坑內(nèi)水位降到坑底，即H也為基坑內(nèi)外水頭差，圖中邊agfg′a′為隔水邊界或相對(duì)隔水邊界，止水帷幕嵌固深度為h，滲流區(qū)土層厚度為T，基坑兩側(cè)止水帷幕間距為2S，則S為基坑兩側(cè)止水帷幕間距的一半。

止水帷幕厚度相對(duì)于基坑開挖深度、嵌固深度等可以忽略不計(jì)，含水層為雙連通區(qū)域，可沿虛線nn′將基坑滲流區(qū)劃分為兩對(duì)稱區(qū)域，且虛線為一條流線，因而可利用圖1中的左邊平面進(jìn)行基坑滲流求解。

對(duì)圖1的左側(cè)平面的單連通區(qū)域（如圖2所示）進(jìn)行分析。對(duì)稱式基坑滲流問題的解析解可以表達(dá)為w（z）=φ（x，y）+iψ（x，y），（z=x+iy）（1）式中：w為復(fù)勢(shì)平面；φ為勢(shì)函數(shù)；ψ為流函數(shù)；z為復(fù)數(shù)變量。

圖2z平面滲流區(qū)示意圖

Fig.2z plane schematic diagram of

seepage field in symmetric foundation pits 在ab邊界上，φ=-kH；

沿agfe邊界上，ψ=q；

沿bcd邊界上，ψ=0；

在de邊界上，φ=0。

為了將圖2中的z平面映射到圖3中的t平面，由施瓦茲克里斯托費(fèi)爾映射公式可得z=S-iP∫t0x-βx+αdxx（x-ε）（x-1）（2）式中：α>0；0<ε<1<β；P>0。

圖3t平面中的滲流區(qū)

Fig.3seepage field in t planef→e：T-HP=∫κ20β-tt+αdtt（ε-t）（1-t）（3a）e→d：SP=∫1κ2β-tt+αdtt（t-ε）（1-t）（3b）d→c：hP=∫β1β-tt+αdtt（t-ε）（t-1）（3c）c→b：H+hP=∫∞βt-βt+αdtt（t-ε）（t-1）（3d）令：ε=κ2，得： κ′=1-κ2（4）式中：κ為模數(shù)；κ′為補(bǔ)模數(shù)。

根據(jù)文[15]可得T-H2P=1+βα∏0-κ2α，κ-K（κ）（5a）

S2P=α+β1+αΠ0κ′2α+1，κ′-K（κ′） （5b）

h2P=β-κ2α+κ2Fβ-1β-κ2，κ-

κ′2（α+β）（α+κ2）（α+1）Πβ-1β-κ2，α+κ21+α，κ（5c）

H+h2P=1+βαΠ1β，-α，κ-

βαF1β，κ（5d）式中：F（z，κ）=∫z0dx（1-x2）（1-κ2x2）為第一類非完全橢圓積分；K（κ）=F（1，κ）為第一類完全橢圓積分；κ為模數(shù)；Π（z，n，κ）=∫z0dx（1+nx2）（1-x2）（1-κ2x2）為第三類非完全橢圓積分；Π0（n，κ）=Π0（1，n，κ）為第三類完全橢圓積分。簡化式（5），可得T-HS=1+βαΠ0-κ2α，κ-Kκα+β1+αΠ0κ′2α+1，κ′-K（κ′）（6a）

hT-H=β-κ2α+κ2Fβ-1β-κ2，κ1+βαΠ0-κ2α，κ-K（κ）

-κ′2（α+β）（α+κ2）（α+1）Πβ-1β-κ2，α+κ21+α，κ1+βαΠ0-κ2α，κ-K（κ）（6b）

H+hT-H=1+βαΠ1β，-α，κ1+βαΠ0-κ2α，κ-K（κ）-

βαF1β，κ1+βαΠ0-κ2α，κ-K（κ）（6c）求解等式（6），可求得參數(shù)α、κ、β的值。

為了使上方程組簡化，令T1=T-H（7）式中：T1即為基坑（相對(duì)）隔水層層頂?shù)娇拥椎木嚯x。

下面將圖2中的滲流區(qū)域映射到圖4中，在ω平面中，ω=φ+iψ，參數(shù)ψc、ψf、ψg、q為未知變量，為了方便計(jì)算，用式（8）（輔助函數(shù)）進(jìn)行變化。ω=ωqK（m′）-K（m）（8）式中：m為模數(shù)；0

K（m）=∫10dx（1-m2x2）（1-x2）（10）通過式（8）可將圖4中的滲流區(qū)域映射到圖5中。

圖4ω平面中的滲流區(qū)

Fig.4seepage field in w plane圖5ω平面中的滲流區(qū)

Fig.5seepage field in ω plane由圖4和圖5中b點(diǎn)對(duì)應(yīng)關(guān)系可得K（m）=kHqK（m′）-K（m）（11）整理得q=12kHK（m′）K（m）（12）為了計(jì)算方便，將圖5平面滲流區(qū)域共形映射到圖6的上半平面。

圖6ξ平面中的滲流區(qū)

Fig.6seepage field in ξ plane最終，可將問題簡化為圖3和圖6之間共形映射關(guān)系的求解，由文[16]知，tξ之間的關(guān)系可通過a、b、d三點(diǎn)的對(duì)應(yīng)關(guān)系進(jìn)行確定，故關(guān)系式可以表達(dá)為t=A+B1-ξ（13）由點(diǎn)a、點(diǎn)d的對(duì)應(yīng)關(guān)系得A+B1-1m=-α

A+B1-（-1）=1 （14）

t=2m+α（m-1）1+m+2（m-1）（α+1）（m+1）（ξ-1） （15）由點(diǎn)e在圖3和圖6中的對(duì)應(yīng)關(guān)系，代入式（15）可得κ2=ε=2m+α（m-1）1+m+2（m-1）（α+1）（m+1）-1m-1（16）整理式（16）得（α+ε）m2+2（ε-α-2）m+ε+α=0（17）求解式（17）得m=α-ε+2-2（1-ε）（1+α）α+ε（18）通過m值，可用式（12）求解對(duì)稱式基坑滲流的單寬涌水量。

前面公式推導(dǎo)過程知，計(jì)算對(duì)稱式基坑滲流單寬涌水量的步驟為

1）先聯(lián)立等式（6）三式，求得參數(shù)α、κ、β的值；

2）再利用式（18），求得m的值；

3）最后利用式（12）求得基坑的單寬涌水量。

通過式（6）、式（18）知，m值與T1S、hT1、HT1有關(guān)。

2影響參數(shù)分析

在上述模型中，影響對(duì)稱式基坑單寬涌水量q計(jì)算的參數(shù)有：基坑內(nèi)外水頭差H、基坑（相對(duì)）隔水層層頂?shù)娇拥椎木嚯x與基坑兩側(cè)止水帷幕間距一半的比值T1S、止水帷幕嵌固深度與（相對(duì)）隔水層層頂?shù)娇拥拙嚯x的比值hT1 、基坑內(nèi)外水頭差H與（相對(duì)）隔水層層頂?shù)娇拥拙嚯x的比值HT1 以及滲流區(qū)土體的滲透系數(shù)k，本文僅對(duì)上述參數(shù)的影響進(jìn)行分析。

圖7對(duì)稱式基坑示意圖

Fig.7schematic diagram of seepage field2.1H的影響

當(dāng)T1S、hT1 、HT1 、k一定時(shí)，由式（12）知：對(duì)稱式基坑的單寬涌水量q與基坑內(nèi)外水頭差H成正比。

2.2T1S的影響

當(dāng)hT1 、HT1 、k、H一定時(shí)，隨T1S的增大，q先保持不變后減小，如圖8所示。當(dāng)T1S≤0.5時(shí)，q基本保持不變，T1S>0.5時(shí)，q大體呈線性降低。

圖8當(dāng)hT1=0.5、HT1=1時(shí)，q值隨T1S的變化曲線

Fig.8The variation curve of q value with the change

of T1S when hT1=0.5 and HT1=12.3hT1 的影響

從圖9可以看出，當(dāng)T1S、HT1 、k、H一定時(shí)，隨著hT1的不斷增大，q逐漸減小，整體趨勢(shì)呈先快后慢再快。值得注意的是，當(dāng)hT1 = 1，即止水帷幕嵌固到（相對(duì)）隔水層中，對(duì)稱式基坑的單寬涌水量q為零。

圖9當(dāng)HT1=0.5、T1S=0.2時(shí)，q值隨hT1的變化曲線

Fig.9The variation curve of q value with the change

of hT1 when HT1=0.5 and T1S=0.22.4HT1 的影響

如圖10所示，當(dāng)hT1、T1S、k、H一定時(shí)，隨著HT1 的不斷增大，q逐漸減小，整體趨勢(shì)呈先快后慢。

圖10當(dāng)hT1=0.5、T1S=0.2時(shí)，q值隨HT1的變化曲線

Fig.10The variation curve of q value with the change

of HT1 when hT1=0.5 and T1S=0.22.5k的影響

由式（12）可知，對(duì)稱式基坑的單寬涌水量q與基坑土體的滲透系數(shù)k成正比，不受H、T1S、hT1、HT1 的影響。

3結(jié)論

假定滲流區(qū)土層為均質(zhì)、各向同性，且地下水運(yùn)動(dòng)符合達(dá)西定律，基于共形映射理論，推導(dǎo)了對(duì)稱式基坑單寬涌水量的計(jì)算表達(dá)式，通過對(duì)其參數(shù)進(jìn)行分析，可以得到以下結(jié)論：

1）當(dāng)T1S、hT1 、HT1 、k一定時(shí)，對(duì)稱式基坑的單寬涌水量q與基坑內(nèi)外水頭差H成正比。

2）當(dāng)hT1 、HT1 、k、H一定時(shí)，隨著T1S的增大，對(duì)稱式基坑的單寬涌水量q先保持不變后減小。

3）當(dāng)T1S、HT1 、k、H一定時(shí)，隨著hT1的不斷增大，對(duì)稱式基坑的單寬涌水量q逐漸減小，整體趨勢(shì)呈先快后慢再快。當(dāng)hT1 = 1，即止水帷幕嵌固到（相對(duì)）隔水層中，對(duì)稱式基坑的單寬涌水量q為零。

4）當(dāng)T1S、hT1、k、H一定時(shí)，隨著HT1 的不斷增大，對(duì)稱式基坑的單寬涌水量q逐漸減小，整體趨勢(shì)呈先快后慢。對(duì)稱式基坑的單寬涌水量q與基坑土體的滲透系數(shù)k成正比，不受H、T1S、hT1、HT1 的影響。

參考文獻(xiàn)：

[1] 童磊，謝康和，盧萌盟，等. 半無限含水層中帶襯砌隧洞滲流解析研究[J]. 巖土力學(xué)， 2011， 32（1）： 304308.

TONG L， XIE K H， LU M M， et al. Analytical study of seepage flow into a lined tunnel in a semiinfinite aquifer [J]. Rock and Soil Mechanics， 2011， 32（1）： 304308. （in Chinese）

[2] 毛海濤，侍克斌，李玉建.無限深透水地基上土石壩壩基垂直防滲的保角變換滲流計(jì)算[J].水利水運(yùn)工程學(xué)報(bào)， 2008（4）： 7177.

Mao H T， SHI K B， LI Y J. Seepage calculation of earthrockfill dam swith vertical cutoffwall on infinite deep pervious foundation by using transformation method [J]. Hydroscience and Engineering， 2008（4）： 7177. （in Chinese）

[3] CHOUDHARY M， CHAHAR B R. Recharge/seepage from an array of rectangular channels[J]. Journal of Hydrology， 2007， 343（1）： 7179.

[4] KOLYMBAS D， WAGNER P. Groundwater ingress to tunnelsthe exact analytical solution [J]. Tunnelling and Underground Space Technology， 2007， 22（1）： 2327.

[5] 曹凈， 張丙軍， 劉海明. 基坑懸掛式止水帷幕滲流的共形映射[J]. 工程勘察，2014（11）： 3438.

CAO J， ZHANG B J， LIU H M. Conformal mapping of unclosed cutoff curtain seepage in the foundation pit[J]. Geotechnical Investigation & Surveying， 2014（11）： 3438. （in Chinese）

[6] 陳永剛，張少琴.保角映射方法在復(fù)合材料斷裂分析中的應(yīng)用[J].太原科技大學(xué)學(xué)報(bào)， 2001， 22（2）： 97101.

CHEN Y G ， ZHANG S Q. Fracture analysis of orthotropic composite plate by conformal mapping method [J]. Journal of Taiyuan University of Science and Technology， 2001， 22（2）： 97101. （in Chinese）

[7] 程海霞， 李俊林.含曲線裂紋的各向異性板斷裂問題[J].中北大學(xué)學(xué)報(bào)（自然科學(xué)版），2012（5）：503508， 529.

CHENG H X， LI J L. Fracture problem in anisotropic composite plate including curve crack[J]. Journal of North University of China（Natural Science Edition）， 2012（5）：503508，529. （in Chinese）

[8] 吳華.L（o）wner理論對(duì)動(dòng)態(tài)斷裂問題的應(yīng)用初探[J]. 重慶電子工程職業(yè)學(xué)院學(xué)報(bào)，2010，19（5）：140143， 165.

WU H. Preliminary exploration on the application of L（o）wner theory in developmental fracture [J]. Journal of Chongqing Vocational & Technical Institute. 2010， 19（5）： 140143，165. （in Chinese）

[9] 王志良，申林方，姚激，高成杰.淺埋隧道圍巖應(yīng)力場(chǎng)的計(jì)算復(fù)變函數(shù)求解法[J].巖土力學(xué)，2010，31（Sup1）： 8690.

WANG Z L， SHEN L F， YAO J， et al. Calculation of stress field in surrounding rocks of shallow tunnel using computational function of complex variable method [J]. Rock and Soil Mechanics， 2010， 31（Sup1）： 8690. （in Chinese）

[10] 王立忠，呂學(xué)金.復(fù)變函數(shù)分析盾構(gòu)隧道施工引起的地基變形[J].巖土工程學(xué)報(bào)，2007， 29（3）：319327.

Wang L Z， LV X J. A complex variable solution for different kinds of oval deformation around circular tunnel in an elastic half plane [J]. Chinese Journal of Geotechnical Engineering， 2007， 29（3）： 319327. （in Chinese）

[11] 蔚立元，陳曉鵬，韓立軍，等.基于復(fù)變函數(shù)方法的水下隧道圍巖彈性分析[J].巖土力學(xué)，2012（Sup2）： 345351.

YU L Y， CHEN X P， HAN L J， et al. Elastic analysis of surrounding rock for underwater tunnels based on functions of complex variables method [J]. Rock and Soil Mechanics， 2012（Sup2）： 345351. （in Chinese）

[12] 任保文.對(duì)稱平板傳輸線分布電容的計(jì)算[J].物理通報(bào)， 2013（1）： 2426.

REN B W. Calculation on distributed capacitance of symmetrical slab transmission line[J]. Bulletin of Physics. 2013（1）： 2426. （in Chinese）

[13] COSTAMAGN E， FANNI A. Needle electrode fields with spacecharge effects[J]. Microwave and Optical Technology Letters， 2000， 26（5）： 321325.

[14] HSU J M， WU M L， LEE C T. Hybrid approach for quasistatic analysis of shielded strip lines[J]. Microwave and Optical Technology Letters， 1997， 14（2）： 111115.

[15] FRIEDLANDER S， SERRE D. Handbook of elliptic integrals for engineers and scientists [M]. Berlin： Springer， 2004.

[16] ELGANAINY M A. Flow underneath a pair of structures with intermediate filters on a drained stratum[J]. Applied Mathematical Modelling， 1986， 10（6）： 394400.