CN109614654A - 一种精确获取基坑降水井影响半径的方法 - Google Patents
一种精确获取基坑降水井影响半径的方法 Download PDFInfo
- Publication number
- CN109614654A CN109614654A CN201811349919.6A CN201811349919A CN109614654A CN 109614654 A CN109614654 A CN 109614654A CN 201811349919 A CN201811349919 A CN 201811349919A CN 109614654 A CN109614654 A CN 109614654A
- Authority
- CN
- China
- Prior art keywords
- radius
- point
- influence
- base pit
- pit dewatering
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 68
- 239000002349 well water Substances 0.000 claims abstract description 16
- 235000020681 well water Nutrition 0.000 claims abstract description 16
- 238000001556 precipitation Methods 0.000 claims abstract description 14
- 238000011144 upstream manufacturing Methods 0.000 claims abstract description 11
- 238000005381 potential energy Methods 0.000 claims description 19
- 239000003673 groundwater Substances 0.000 claims description 7
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 6
- 230000008595 infiltration Effects 0.000 claims description 6
- 238000001764 infiltration Methods 0.000 claims description 6
- 239000004567 concrete Substances 0.000 claims description 5
- 239000002689 soil Substances 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 229910052742 iron Inorganic materials 0.000 claims description 3
- 230000033228 biological regulation Effects 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 230000000052 comparative effect Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 239000004927 clay Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000010998 test method Methods 0.000 description 3
- 241000209202 Bromus secalinus Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical group O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003204 osmotic effect Effects 0.000 description 1
- 239000011150 reinforced concrete Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种精确获取基坑降水井影响半径的方法。本发明将地下水原始稳定水位看作上游边界水头，降水井水位看作下游边界水头，控制水位与降水井水位之间的出渗段为逸出边界，降水影响半径为渗流场的水平尺寸，将求解降水井影响半径简化成已知上、下游边界水头及逸出边界求解模型水平尺寸。相比现有技术，本发明的方法误差更小、精度更高，适用于工程中的多种水文地质条件，能够精确、高效地计算出降水影响半径。
Description
技术领域
本发明涉及一种获取基坑降水井影响半径的方法，属于基坑工程技术领域。
背景技术
基坑降水影响半径是指在基坑降水过程中，基坑周围形成降水漏斗，当地下水水位降幅为零时，降水漏斗边缘与降水井中心之间的距离，如图1所示。降水影响半径是基坑涌水量的一个重要参数，它的取值直接影响到基坑总涌水量的准确性，进而影响到基坑降水方案设计的合理性。基坑降水方案设计是基坑地下水控制的关键环节，合理的基坑降水方案可以有效地控制地下水。这样不仅可以减少管涌、流土等工程灾害的发生，还能够减轻由于不合理降水而引发的基坑周围不均匀沉降、建筑倾斜等破坏。
对于潜水层基坑降水影响半径的获得，《建筑基坑支护技术规程》(JGJ120-2012)(以下简称“规程”)推荐由现场抽水试验结合《基坑降水手册》中的公式确定。
但是，规程采用的是经验公式，仅依赖于渗透系数和井水位降深的大小，对水跃现象等影响因素考虑不足；经验值法是一种估算法，计算结果存在较大误差；图解法多采用s-lgr直线交汇法，观测孔水位降的测量及绘图微小误差都可能导致计算结果几米到几十米的偏差。
发明内容
为解决现有技术的上述不足，获得更准确的基坑降水影响半径，本发明提供了一种新的获取基坑降水井影响半径的方法。相比现有技术，本发明的方法误差更小、精度更高，适用于工程中的多种水文地质条件，能够精确、高效地计算出降水影响半径。
一种获取基坑降水井影响半径的方法，包括下述步骤：
(1)计算水跃值Δh，确定基坑降水控制水位的高度sd，
其中，单井涌水量q＝120π0.5dlk/3，
过滤器工作部分表面积F＝πdl，
d、l分别为过滤器的直径和长度，
α为过滤器构造影响系数，
k为含水层的渗透系数，根据具体工程案例的地质勘查报告确定，
sw为降水井水位降深，根据具体工程案例的需要确定，
然后，确定基坑降水控制水位的高度sd：
sd＝b+Δh (2)
b为降水井的井中水位到不透水层之间的距离；
(2)根据图2所示的均质潜水含水层基坑的二维稳定渗流场，建立图3所示的基坑降水影响半径计算模型，
2.1)设A1为上游水位点即地下水水位降幅为零的点、B1点为A1点在不透水边界的垂直投影，D1为下游水位点即降水井水位点、C1为D1点在不透水边界的垂直投影点，F1为降水影响半径计算模型所对应的逸出点，
假设a为模型水平尺寸即基坑降水影响半径R，则a＝R，
B1A1为上游边界水头h′，D1C1为下游边界水头h″，
并且，h′＝H，h″＝b，H为潜水含水层厚度；
2.2)然后，以B1为坐标原点，即以上游水位点在不透水边界的垂直投影点为坐标原点，B1C1为x轴，B1A1为z轴建立直角坐标系，即基坑降水影响半径计算模型，其中，所述模型高度为H，
记p为基坑降水影响半径模型水平尺寸时需要调整的次数，初始时设p＝0；
(3)所述基坑降水影响半径计算模型的有限单元划分，
对所述基坑降水影响半径计算模型进行n行、u列剖分，得到n×u个四边形的有限单元，再对各四边形的有限单元进行同方向对角线连接，进一步剖分为2×n×u个三角形的有限单元，各节点按照由上向下、自左向右的顺序依次编号为1，2，3……(n+1)×(u+1)，其中，n、u为非零自然数，所述各节点是指各有限单元在顶角处相互联接的联接点；
(4)求解基坑降水影响半径计算模型的渗流场的逸出点的高度sp
4.1)设h为步骤3)划分的各有限单元节点的水头值，则h满足边界条件，即边界上的h等于该有限单元节点的纵坐标，
4.2)假定任意逸出点Sm(a,zm)的水头值为hm＝zm，其中，m为假定逸出点的选取序号，依次取值为1，2，3……m,
然后，将步骤2.1)的所述上游边界水头h′、下游边界水头h″和假定逸出点的水头值hm，且h″≤hm≤h′，代入线性方程组：
[K]{h}＝{f} (3)
其中，[K]为总渗透矩阵，
{h}为待求各有限单元节点的水头值，其中：
{f}为自由项，根据各单元渗透矩阵及边界条件确定，
然后，求解线性方程组(3)即可得到所述基坑降水影响半径计算模型内各节点的水头值h；
4.3)计算步骤3)划分的各有限单元节点的水头值h与各节点位置纵坐标z的差值ΔL并划分虚区和实区，
当有限单元内各节点的差值ΔL≤0时，则该有限单元位于虚区，
否则，该有限单元位于实区；
4.4)将虚区内的有限单元丢弃，由下列公式(5)计算当前假定逸出点Sm(a,zm)位置时的实域总势能ESm，
其中，E为实域总势能，为土的孔隙率，ρ为水的密度，g为重力加速度，hi、hj、hm分别为三角形单元三个节点的水头值，A为三角形单元的面积，
4.5)将假定逸出点Sm(a,zm)位置升高一个节点Sm+1(a，zm+1)，其水头值为hm+1＝zm+1，h″＜hm+1＜h′，按照步骤4.2)、4.3)、4.4)相同操作求解出假定的逸出点Sm+1(a，zm+1)的实域总势能
判断ESm与的大小，若说明真实逸出点在假定逸出点Sm(a,zm)、Sm+1(a，zm+1)之间的区段内，执行下一步骤4.6)，
若说明真实逸出点不在假定逸出点Sm(a,zm)、Sm+1(a，zm+1)之间的区段内，返回步骤4.5)继续升高一个节点直至步骤4.2)中假定的逸出点高度为H时停止，即假定逸出点Sm(a,zm)已经超出了边界条件；
4.6)在假设的逸出点Sm(a,zm)、Sm+1(a，zm+1)的区段内，在Sm(a,zm)点的基础上升高β，重新开始假定逸出点Yt(a，z′t)，其中，z′t＝zm+β，β为假设逸出点最小误差值，t依次取值为1，2，3……t，
然后，依次按照步骤4.2)、4.3)、4.4)相同操作求解实域总势能
4.7)若说明真实逸出点不在假定逸出点Sm(a,zm)、Yt(a，z′t)之间的区段内，令t＝t+1，执行步骤4.8)继续升高一个β值，
若此时所对应的逸出点Sm(a,zm)的纵坐标即为真实逸出点高度sp，即sp＝zm，并执行步骤5)，
4.8)在Yt(a，z′t)点的基础上升高一个β值，得到新的假定逸出点Yt+1(a，z′t+1)，其中，z′t+1＝z′t+β，然后，依次按照步骤4.2)、4.3)、4.4)相同操作求解实域总势能
4.9)若说明真实逸出点不在假定逸出点Yt(a，z′t)、Yt+1(a，z′t+1)之间的区段内，令t＝t+1，返回到步骤4.8)继续升高一个β值，
若此时所对应的逸出点Yt(a，z′t)的纵坐标即为真实逸出点高度sp，即sp＝z′t，并执行步骤5)；
(5)比较sp与sd的大小，若其中，Δe为误差限值，则步骤2)中所述基坑降水影响半径模型假设的水平尺寸a即为真实的基坑降水影响半径R，
若或执行步骤6)；
(6)若说明步骤2)中所述基坑降水影响半径模型假设的水平尺寸a偏小，则令p＝p+1，重新假设基坑降水影响半径计算模型水平尺寸为(1+2-p)a，此处a为前一个模型假设的水平尺寸，并保持有限单元总数不变，增加各有限单元的水平尺寸，执行步骤3)相同操作对重新假设的基坑降水影响半径计算模型进行有限单元划分，
若说明步骤2)中所述基坑降水影响半径模型假设的水平尺寸a偏大，则令p＝p+1，重新假设基坑降水影响半径模型水平尺寸为(1-2-p)a，此处a为前一个模型假设的水平尺寸，并保持有限单元总数不变，减小各有限单元的水平尺寸，执行步骤3)相同操作对重新假设的基坑降水影响半径计算模型进行有限单元划分。
本发明方法的流程图如图4所示。
优选，所述假设逸出点最小误差值β≤0.1，所述误差限值Δe≤5％，设置偏于保守的误差取值限制条件。
优选，步骤(3)的基坑降水影响半径计算模型的有限单元划分为四边形的有限单元，具体如下所述：
对所述基坑降水影响半径计算模型进行n行、u列剖分，得到n×u个四边形的有限单元，各节点按照由上向下、自左向右的顺序依次编号为1，2，3……(n+1)×(u+1)，其中，n、u为非零自然数，降低计算复杂度。
优选，网状和砾石过滤器α取值为20，对缝隙和铁丝过滤器α取值为7。
与现有技术相比，本发明的技术效果是：可以精确、高效地计算影响半径，提高基坑涌水量的计算精度，使基坑降水设计方案更加合理。
并且，
(1)当基坑降水达到控制水位时，将地下水原始稳定水位看作上游边界水头，降水井水位看作下游边界水头，控制水位与降水井水位之间的出渗段(水跃区段)看作逸出边界，降水影响半径看作渗流场的水平尺寸，求解降水影响半径便可简化成已知上、下游边界水头及逸出边界，求解模型水平尺寸的问题。
(2)给出了一种新的降水影响半径的计算方法，同时将Δe取为+5％，即所求影响半径计算结果总是偏小5％左右，从而使基坑涌水量计算结果总是偏大，本专利的方法是一种偏于保守的取值方法。
附图说明
图1为现有技术中均质潜水含水层基坑降水影响半径的示意图，
图2为现有技术中均质潜水含水层基坑的二维稳定渗流场剖面图，
图3为本发明一种获取基坑降水井影响半径的方法建立的基坑降水影响半径计算模型示意图，
图4为本发明一种获取基坑降水井影响半径的方法的流程图，
图5为实施例1的基坑降水井点布置图，
图6(a)为实施例1的基坑二维稳定渗流场剖面图，
图6(b)为实施例1的基坑建立的基坑降水影响半径计算模型示意图，
图7实施例1建立的基坑降水影响半径计算模型有限元划分及实域总势能示意图。
其中，图1-7中的数字或字母标示的指代说明如下：
100为基坑，200为降水井，R为基坑降水影响半径，r0为基坑的引用半径，sw为降水井水位降深，H为潜水含水层厚度，Δh为井壁水位与降水井水位之差即水跃值，k为含水层的渗透系数，b为降水井井中水位到不透水层之间的距离，sd为基坑降水控制水位的高度，sp为计算模型真实逸出点的高度，z为坐标轴点，A表示从A点往左地下水水位降幅为零，G点为地下水原始水位线与降水井井壁的交点，O为与A在同一水平面上的井孔中心点，B点为A点在不透水层上的垂直投影，C点为降水井与不透水层的交点，D点为降水井水位点，F点为距坑底0.5m的基坑降水控制水位点，A1为上游水位点，G1为不透水边界交点，B1点为A1点在不透水边界的垂直投影，D1为下游水位点，C1为D1点在不透水边界的垂直投影，F1'、F1″为假定逸出点，F1为此模型所对应的真实逸出点，E为实域总势能，o为坐标原点，a为基坑降水影响半径计算模型水平尺寸，1′、2′、3′………200′表示单元序号，1、2……121表示节点序号，4.0m静止水位深，13.3m为上游水位高度，5.94m为下游水位高度，7.3m为基坑降水控制水位的高度，1.36m为水跃值，7.36m为降水深度，12.07为真实逸出点高度，33491.12为实域总势能值，为实域总势能单位。
具体实施方式
以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
以具体实例为例，该例具体情况如下：
某钢铁厂氧化铁皮坑工程为直径18m圆形钢筋混凝土结构，深度15.8m。
设计要求降水深度为-10m(自然地面设为0.00m)，静止水位深4m，由于存在水跃值，实际降水深度sw＝7.36m，图6(a)为该例的工程的基本概况，降水面积为F0＝1900m2，井点布置见图5。
在深度20m范围内的土层主要由粉质黏土所组成，但在不同深度的粉质黏土层内夹有薄层黏质粉土和砂层以及黏土层；地面以下4.00m为静止地下水位，按照规程规定可将以上不同含水层概化成一层含水层，其渗透系数按各含水层厚度加权平均，计算得含水层的渗透系数k＝5.7m/d，含水层厚度H＝13.3m。
井点过滤器半径r＝0.054m，d＝2r，长度l＝1.5m。
单个井点抽水量：q＝58m3/d。
实施例1：
根据上述基坑100，本发明的方法计算降水井200的影响半径具体步骤如下所述：
1)计算水跃值Δh，含水层厚度H＝13.3m，考虑到水跃值的存在，水位降深sw＝7.36m，井点过滤器半径r＝0.054m，长度l＝1.5m，实测每个井点抽水量q＝58m3/d。F＝πdl＝π×0.054×2×1.5＝0.509，α＝20，水跃值Δh由公式(1)确定：
b＝5.94m，由公式(2)确定基坑降水控制水位的高度sd：
sd＝b+Δh＝5.94+1.36＝7.3m
得，计算水跃值Δh为1.36m，基坑降水控制水位高度sd＝7.3m。
2)建立如图6(b)所示的基坑降水影响半径计算模型，
设A1为上游水位点即地下水水位降幅为零的点、B1点为A1点在不透水边界的垂直投影、D1为下游水位点即降水井水位点，
C1为D1点在不透水边界的垂直投影点、F1为模型所对应的真实逸出点，a为模型水平尺寸即基坑降水影响半径R。B1A1为上游边界水头，B1C1为下游边界水头。
本例中第一次假设时，假设a＝50m，为模型水平尺寸即基坑降水影响半径R，B1A1为上游边界水头h′＝13.3m，D1C1为下游边界水头h″＝5.94m，
然后，以B1为坐标原点，B1C1为x轴，B1A1为z轴建立直角坐标系，设模型高度为H即潜水含水层厚度。第一次对a值假设时，记p＝0，p为计算降水影响半径模型水平尺寸时需要调整的次数。
3)基坑降水影响半径计算模型的有限单元划分，
对所述基坑降水影响半径计算模型进行n行、u列剖分，得到n×u个四边形的有限单元，再对各四边形的有限单元进行同方向对角线连接，进一步剖分为2×n×u个三角形的有限单元，各节点按照由上向下、自左向右的顺序依次编号为1，2，3……(n+1)×(u+1)，采用直角边边长为5m×1.33m的三角形单元将计算模型划分为图7所示的2×10×10＝200个单元。节点按照由上向下，由左向右的顺序依次编号为1，2，3，……，121。
4)求解基坑降水影响半径计算模型的渗流场逸出点的高度sp，
设h为步骤3)划分的各有限单元节点的水头值，h满足边界条件，即边界上的h等于该节点的纵坐标，
假定任意逸出点Sm(a,zm)的水头值为hm＝zm，即5.94m≤hm≤13.3m，其中，m为假定逸出点的选取序号，依次取值为1，2，3……121，由图7可知，本例中假定逸出节点有7个，纵坐标分别为5.94m、6.65m、7.98m、9.31m、10.64m、11.97m、13.3m。
不管先取哪个点为逸出点，每个假定的逸出点所对应的实域总势能不会变。随机取一个为假定逸出点，以纵坐标是7.98m的逸出点为例，取S1(50，7.98)。
计算步骤3)划分的各有限单元节点的水头值h与各节点位置纵坐标z的差值ΔL并划分虚区和实区，
当有限单元内各节点的差值ΔL≤0时，则该有限单元位于虚区，通过计算机计算得到本实施例中位于虚区的有限单元，在图7中的其代号分别是1′，22′，23′，41′，42′，61′，62′，63′，81′，82′，83′，84′，101′，102′，103′，104′，105′，121′，122′，124′，125′，126′，141′，142′，143′，144′，145′，146′，161′，162′，163′，164′，165′，166′，167′，181′，182′，183′，184′，185′，186′，187′，188′；
否则，该有限单元位于实区，本实施例中去除虚区单元剩余单元代号位于实区；
去掉虚区有限单元，用公式(5)计算当前假定逸出点位置时的实域总势能ESm，
本例中，为土的孔隙率，ρ为水的密度1g/cm3，g为重力加速度10m/s2，hi、hj、hm分别为三角形单元三个节点的水头值，例如有限单元3′的三个节点水头值分别是88.44、80.922386、80.922394，A为三角形单元的面积为3.325㎡，
用公式(5)计算出实域总势能为
将假定逸出点Sm(a,zm)位置升高一个节点S2(50，9.31)，用公式(3)计算出各节点水头值h，去掉虚区单元，由公式(5)计算出实域总势能为
由于ES2＜ES1，说明真实逸出点不在假定逸出点S1(50，7.98)、S2(50，9.31)之间的区段内，继续升高一个节点值重新假定逸出点，直至假定的逸出点水头值为hm+1大于13.3m时停止，即重新假定逸出点Sm(a,zm)已经超出了边界条件，计算的各逸出点对应的实域总势能如表1所示，
表1实域总势能
由表1可知说明真实逸出点在假定逸出点S4(50,11.97)、S5(50，13.3)之间的区段内，
在假设点逸出点S4(50,11.97)、S5(50，13.3)的区段内，从S4(50,11.97)点开始升高β，重新开始假定逸出点Y1(50，z′1)，其中，z′1＝11.97+0.1＝12.07，β为假设逸出点最小误差值，取β＝0.1，
然后，按照上述相同操作求解实域总势能
由于EY1＜ES4，说明真实逸出点不在假定逸出点S4(50,11.97)、Y1(50，12.07)之间的区段内，令t＝t+1，在Y1(50，12.07)点的基础上升高β，得到新的假定逸出点Y2(50，z′2)，其中，z′2＝z′1+β＝12.07+0.1＝12.17，然后，求解实域总势能
由于EY2＞EY1，此时EY1所对应的逸出点Y1(50，12.07)的纵坐标即为真实逸出点高度sp，即sp＝12.07m；
5)此时模型水平尺寸为50m，将真实逸出点高度12.07m与基坑降水控制水位高7.3m进行比较：
(12.07-7.3)/7.3＝0.653×100％＝65.3％＞Δe＝5％
说明此时基坑降水影响半径模型假设的水平尺寸50m偏小，重新假设计算模型水平尺寸为(1+2-1)×50＝75m，并保持有限单元总数不变，增加各有限单元的水平尺寸，执行步骤3)，
计算得真实逸出点高度10.0m与基坑降水控制水位高7.3m进行比较，误差为36.99％＞Δe＝5％，说明此时基坑降水影响半径模型假设的水平尺寸75m偏小，重新假设计算模型水平尺寸为(1+2-2)×75＝93.75m，并保持有限单元总数不变，增加各有限单元的水平尺寸，执行步骤3)，直至基坑降水影响半径模型的真实逸出点高度与基坑降水控制水位高度之间的误差满足要求，计算结果如表2所示，
表2模型逸出点高度sp与控制水位高度sd对比
6)当模型水平尺寸为104.17m时，计算的逸出点为7.4m。与基坑降水控制水位高度7.3m进行比较：
(7.4-7.3)/7.3＝0.0137×100％＝1.37％＜Δe＝5％，
则，满足误差限值Δe，水平尺寸104.17m为真实的基坑降水影响半径R，故求解出影响半径R＝104.17m。
此时，误差仅为1.37％。
实施例2：
与实施例1同一基坑，采用与实施例1相同操作，不同仅在于首次假设模型水平尺寸从a＝60m开始，按照上述步骤2)、3)、4)、5)、6)进行操作，结果如表3所示，
表3模型逸出点与控制水位高度sd对比
当模型水平尺寸为105m时，计算的逸出点为7.39m。与基坑降水控制水位高度7.3m进行比较：
(7.39-7.3)/7.3＝0.012×100％＝1.2％＜Δe＝5％
同样满足误差限值Δe，求解出影响半径R＝105m，误差仅为1.2％。
可见，实施例1和2对同一基坑得出的结果相差不大，说明，在满足误差要求Δe的前提下本发明的方法求解影响半径R都是有效的，首次假设模型水平尺寸a值的大小不对本发明方法精度有影响，即首次假设模型水平尺寸a可以任意取值。Δe取为+5％，即所求影响半径计算结果总是偏小5％左右，由于涌水量计算与影响半径呈反比例关系，从而使基坑涌水量计算结果总是偏大，是一种偏于保守的误差取值限制条件。
对比例1：
与实施例1同一基坑，经验公式求解降水影响半径R
降水影响半径R按照经验公式确定：
式中，R为影响半径(m)；sw为井水位降深(m)；当井水位降深小于10m时，取sw＝10m；k为含水层的渗透系数(m/d)；H为潜水含水层厚度(m)。
故：
对比例2：
与实施例1同一基坑，现场抽水试验法求解降水影响半径R
1)计算引用半径r0：
2)
其中，Q为基坑总出水量(m3/d)，本工程总抽水量为1750m3/d；
r0为引用半径(m)。
故由抽水试验法计算的降水影响半径：
R＝102.01m≈103m。
与对比例1和2相比，本发明(以实施例1为例)获得的影响半径与现场抽水试验结果和规程公式解对比情况如表4所示。
表4实施例1精度对比
由表4可知，经验公式计算的影响半径与现场抽水试验计算的影响半径相比，经验公式的绝对误差为68.67m，相对误差为66.7％；本发明的计算结果与现场抽水试验法的计算结果相比，本发明绝对误差为1.17m，相对误差为1.1％，更接近现场抽水试验结果。故本发明可提高基坑总涌水量的计算精度，使基坑降水设计方案更加合理，具有实际的工程应用价值。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种获取基坑降水井影响半径的方法，包括下述步骤：
(1)计算水跃值Δh，确定基坑降水控制水位的高度sd，
其中，单井涌水量q＝120π0.5dlk/3，
过滤器工作部分表面积F＝πdl，
d、l分别为过滤器的直径和长度，
α为过滤器构造影响系数，
k为含水层的渗透系数，根据具体工程案例的地质勘查报告确定，
sw为抽水井水位降深，根据具体工程案例的需要确定，
然后，确定基坑降水控制水位的高度sd：
sd＝b+Δh (2)
b为降水井的井中水位到不透水层之间的距离；
其特征在于，以下步骤：
(2)建立基坑降水影响半径计算模型，
2.1)设A1为上游水位点即地下水水位降幅为零的点、B1点为A1点在不透水边界的垂直投影，D1为下游水位点即降水井水位点、C1为D1点在不透水边界的垂直投影点，F1为降水影响半径计算模型所对应的逸出点，
先假设a为模型水平尺寸即基坑降水影响半径R，
B1A1为上游边界水头h′，D1C1为下游边界水头h″，
并且，h′＝H，h″＝b，H为潜水含水层厚度；
2.2)然后，以B1为坐标原点，即以上游水位点在不透水边界的垂直投影点为坐标原点，B1C1为x轴，B1A1为z轴建立直角坐标系，即建立基坑降水影响半径计算模型，其中，所述模型高度为H，
记p为基坑降水影响半径模型水平尺寸时需要调整的次数，初始时设p＝0；
(3)所述基坑降水影响半径计算模型的有限单元划分，
对所述基坑降水影响半径计算模型进行n行、u列剖分，得到n×u个四边形的有限单元，再对各四边形的有限单元进行同方向对角线连接，进一步剖分为2×n×u个三角形的有限单元，各节点按照由上向下、自左向右的顺序依次编号为1，2，3……(n+1)×(u+1)，其中，所述各节点是指各有限单元在顶角处相互联接的联接点；
(4)求解基坑降水影响半径计算模型的渗流场逸出点的高度sp，
4.1)设h为步骤3)划分的各有限单元节点的水头值，则h满足边界条件，即边界上的h等于该节点的纵坐标，
4.2)假定任意逸出点Sm(a,zm)的水头值为hm＝zm，其中，m为假定逸出点的选取序号，依次取值为1，2，3……m,
然后，将步骤2.1)的所述上游边界水头h′、下游边界水头h″和假定逸出点的水头值hm，且h″≤hm≤h′，代入线性方程组：
[K]{h}＝{f} (3)
其中，[K]为总渗透矩阵，
{h}为待求节点的水头值，其中：
{f}为自由项，根据各单元渗透矩阵及边界条件确定，
然后，求解线性方程组(3)即可得到所述基坑降水影响半径计算模型内各节点的水头值h；
4.3)计算步骤3)划分的各有限单元节点的水头值h与各节点位置纵坐标z的差值ΔL并划分虚区和实区，
当有限单元内各节点的差值ΔL≤0时，则该有限单元位于虚区，
否则，该有限单元位于实区；
4.4)将虚区内的有限单元丢弃，计算当前假定逸出点位置时的实域总势能ESm，
其中，E为实域总势能，为土的孔隙率，ρ为水的密度，g为重力加速度，hi、hj、hm分别为三角形单元三个节点的水头值，A为三角形单元的面积，
4.5)将假定逸出点Sm(a,zm)位置升高一个节点Sm+1(a，zm+1)，其水头值为hm+1＝zm+1，h″＜hm+1＜h′，按照步骤4.2)、4.3)、4.4)相同操作求解出假定的逸出点Sm+1(a，zm+1)的实域总势能
判断ESm与的大小，若说明真实逸出点在假定逸出点Sm(a,zm)、Sm+1(a，zm+1)之间的区段内，执行步骤4.6)，
若说明真实逸出点不在假定逸出点Sm(a,zm)、Sm+1(a，zm+1)之间的区段内，返回步骤4.5)继续升高一个节点直至步骤4.2)中假定的逸出点水头值为hm+1大于H时停止，即假定逸出点Sm(a,zm)已经超出了边界条件；
4.6)在假设的逸出点Sm(a,zm)、Sm+1(a，zm+1)的区段内，在Sm(a,zm)点的基础上升高β，重新开始假定逸出点Yt(a，z′t)，其中，z′t＝zm+β，β为假设逸出点最小误差值，t依次取值为1，2，3……t，
然后，依次按照步骤4.2)、4.3)、4.4)相同操作求解实域总势能
4.7)若说明真实逸出点不在假定逸出点Sm(a,zm)、Yt(a，z′t)之间的区段内，令t＝t+1，执行步骤4.8)，
若此时所对应的逸出点Sm(a,zm)的纵坐标即为真实逸出点高度sp，即sp＝zm，并执行步骤5)，
4.8)在Yt(a，z′t)点的基础上升高β，得到新的假定逸出点Yt+1(a，z′t+1)，其中，z′t+1＝z′t+β，然后，依次按照步骤4.2)、4.3)、4.4)相同操作求解实域总势能
4.9)若说明真实逸出点不在假定逸出点Yt(a，z′t)、Yt+1(a，z′t+1)之间的区段内，令t＝t+1，返回到步骤4.8)，
若此时所对应的逸出点Yt(a，z′t)的纵坐标即为真实逸出点高度sp，即sp＝z′t，并执行步骤5)；
(5)比较sp与sd的大小，若其中，Δe为误差限值，则步骤2)中所述基坑降水影响半径模型假设的水平尺寸a即为真实的基坑降水影响半径R，
若或执行步骤6)；
(6)若说明步骤2)中所述基坑降水影响半径模型假设的水平尺寸a偏小，则令p＝p+1，重新假设基坑降水影响半径模型水平尺寸为(1+2-p)a，此处a为前一个模型假设的水平尺寸，并保持有限单元总数不变，增加各有限单元的水平尺寸，执行步骤3)对新假设的基坑降水影响半径模型进行有限单元的划分，
若说明步骤2)中所述基坑降水影响半径模型假设的水平尺寸a偏大，则令p＝p+1，重新假设基坑降水影响半径模型水平尺寸为(1-2-p)a，此处a为前一个模型假设的水平尺寸，并保持单元总数不变，减小各单元的水平尺寸，执行步骤3)对新假设的基坑降水影响半径模型进行有限单元的划分。
2.根据权利要求1所述的一种获取基坑降水影响半径的方法，其特征在于，所述假设逸出点最小误差值β≤0.1，所述误差限值Δe≤5％。
3.根据权利要求1所述的一种获取基坑降水影响半径的方法，其特征在于，步骤(3)的基坑降水影响半径计算模型的有限单元划分为四边形的有限单元，具体如下所述：
对所述基坑降水影响半径计算模型进行n行、u列剖分，得到n×u个四边形的有限单元，各节点按照由上向下、自左向右的顺序依次编号为1，2，3……(n+1)×(u+1)，其中，n、u为非零自然数。
4.根据权利要求1所述的一种获取基坑降水影响半径的方法，其特征在于，网状和砾石过滤器α取值为20，对缝隙和铁丝过滤器α取值为7。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811349919.6A CN109614654B (zh) | 2018-11-13 | 2018-11-13 | 一种精确获取基坑降水井影响半径的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811349919.6A CN109614654B (zh) | 2018-11-13 | 2018-11-13 | 一种精确获取基坑降水井影响半径的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109614654A true CN109614654A (zh) | 2019-04-12 |
| CN109614654B CN109614654B (zh) | 2022-11-01 |
Family
ID=66003137
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811349919.6A Active CN109614654B (zh) | 2018-11-13 | 2018-11-13 | 一种精确获取基坑降水井影响半径的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109614654B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110387896A (zh) * | 2019-06-28 | 2019-10-29 | 中交路桥北方工程有限公司 | 基坑降水井设计方法 |
| CN110863504A (zh) * | 2019-12-02 | 2020-03-06 | 河南科技大学 | 一种管井降水施工方法 |
| CN111898283A (zh) * | 2020-08-14 | 2020-11-06 | 中铁十六局集团有限公司 | 隧道形变量的确定方法、装置、计算机设备及存储介质 |
| US20210181375A1 (en) * | 2019-12-12 | 2021-06-17 | China Institute Of Water Resources And Hydropower Research | Numerical method for simulating a karez well in association with a groundwater model |
| CN114067095A (zh) * | 2021-11-29 | 2022-02-18 | 黄河勘测规划设计研究院有限公司 | 基于水尺字符检测识别的水位识别方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005250599A (ja) * | 2004-03-01 | 2005-09-15 | Fuji Heavy Ind Ltd | 仮想境界法を用いた流れ解析装置、及び、仮想境界法を用いた流れ解析方法 |
| JP2006249764A (ja) * | 2005-03-10 | 2006-09-21 | Shimizu Corp | 地下水揚水時の許容揚水量の算出方法 |
| CN101092818A (zh) * | 2007-06-07 | 2007-12-26 | 上海交通大学 | 基坑潜水与承压含水层一井分层联合降水的方法 |
| CN105239611A (zh) * | 2015-10-22 | 2016-01-13 | 上海交通大学 | 确定基坑开挖面下方止水帷幕渗漏对周边环境影响的方法 |
| CN105386474A (zh) * | 2015-10-22 | 2016-03-09 | 上海交通大学 | 确定基坑开挖面上方止水帷幕渗漏对周边环境影响的方法 |
| CN107451399A (zh) * | 2017-07-10 | 2017-12-08 | 烟台大学 | 一种基于实域总势能计算渗流溢出点的方法 |
| CN107491624A (zh) * | 2017-09-20 | 2017-12-19 | 西安建筑科技大学 | 一种利用虚拟材料等效栓接结合部动态特性的方法 |
| CN108090291A (zh) * | 2017-12-23 | 2018-05-29 | 北京地矿工程建设有限责任公司 | 一种预测承压水悬挂式帷幕基坑涌水量的计算方法 |
| CN108549774A (zh) * | 2018-04-17 | 2018-09-18 | 中南大学 | 非均质承压水层悬挂式帷幕基坑定水位抽水量的确定方法 |
-
2018
- 2018-11-13 CN CN201811349919.6A patent/CN109614654B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005250599A (ja) * | 2004-03-01 | 2005-09-15 | Fuji Heavy Ind Ltd | 仮想境界法を用いた流れ解析装置、及び、仮想境界法を用いた流れ解析方法 |
| JP2006249764A (ja) * | 2005-03-10 | 2006-09-21 | Shimizu Corp | 地下水揚水時の許容揚水量の算出方法 |
| CN101092818A (zh) * | 2007-06-07 | 2007-12-26 | 上海交通大学 | 基坑潜水与承压含水层一井分层联合降水的方法 |
| CN105239611A (zh) * | 2015-10-22 | 2016-01-13 | 上海交通大学 | 确定基坑开挖面下方止水帷幕渗漏对周边环境影响的方法 |
| CN105386474A (zh) * | 2015-10-22 | 2016-03-09 | 上海交通大学 | 确定基坑开挖面上方止水帷幕渗漏对周边环境影响的方法 |
| CN107451399A (zh) * | 2017-07-10 | 2017-12-08 | 烟台大学 | 一种基于实域总势能计算渗流溢出点的方法 |
| CN107491624A (zh) * | 2017-09-20 | 2017-12-19 | 西安建筑科技大学 | 一种利用虚拟材料等效栓接结合部动态特性的方法 |
| CN108090291A (zh) * | 2017-12-23 | 2018-05-29 | 北京地矿工程建设有限责任公司 | 一种预测承压水悬挂式帷幕基坑涌水量的计算方法 |
| CN108549774A (zh) * | 2018-04-17 | 2018-09-18 | 中南大学 | 非均质承压水层悬挂式帷幕基坑定水位抽水量的确定方法 |
Non-Patent Citations (5)
| Title |
|---|
| K R RAKESHASHARMA: "A literature survey on silt erosion and cavitation in hydro turbine", 《2016 INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING, COMMUNICATION, POWER AND EMBEDDED SYSTEM (SCOPES)》 * |
| SHEN CY: "DLVO Interaction Energies between Hollow Spherical Particles and Collector Surfaces", 《LANGMUIR》 * |
| 张晓欣: "基坑降水数值模拟研究", 《科技风》 * |
| 张莲花: "基坑降水引起的沉降变形时空规律及降水控制研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
| 李远东: "渗流有限元数值计算及其在基坑工程中的应用", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110387896A (zh) * | 2019-06-28 | 2019-10-29 | 中交路桥北方工程有限公司 | 基坑降水井设计方法 |
| CN110387896B (zh) * | 2019-06-28 | 2021-04-09 | 中交路桥北方工程有限公司 | 基坑降水井设计方法 |
| CN110863504A (zh) * | 2019-12-02 | 2020-03-06 | 河南科技大学 | 一种管井降水施工方法 |
| US20210181375A1 (en) * | 2019-12-12 | 2021-06-17 | China Institute Of Water Resources And Hydropower Research | Numerical method for simulating a karez well in association with a groundwater model |
| CN111898283A (zh) * | 2020-08-14 | 2020-11-06 | 中铁十六局集团有限公司 | 隧道形变量的确定方法、装置、计算机设备及存储介质 |
| CN111898283B (zh) * | 2020-08-14 | 2024-05-28 | 中铁十六局集团有限公司 | 隧道形变量的确定方法、装置、计算机设备及存储介质 |
| CN114067095A (zh) * | 2021-11-29 | 2022-02-18 | 黄河勘测规划设计研究院有限公司 | 基于水尺字符检测识别的水位识别方法 |
| CN114067095B (zh) * | 2021-11-29 | 2023-11-10 | 黄河勘测规划设计研究院有限公司 | 基于水尺字符检测识别的水位识别方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109614654B (zh) | 2022-11-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109614654A (zh) | 一种精确获取基坑降水井影响半径的方法 | |
| CN110046469B (zh) | 多约束条件下水电站坝前河床冲淤变形的计算方法 | |
| CN107587516B (zh) | 一种跟进式高效止水阻砂降水复合系统及其施工方法 | |
| CN108824505B (zh) | 一种近距离沉井群下沉施工方法 | |
| CN106018739A (zh) | 一种潮滩-潮沟系统地貌演变物理模型试验系统及方法 | |
| KR101817560B1 (ko) | 지하수 인공함양 실험 시스템 | |
| CN102902844A (zh) | 基于大数据量dem数据的子流域划分方法 | |
| CN109253717B (zh) | 一种矿区地表沉降三维激光扫描地表沉降监测设站方法 | |
| CN108446489A (zh) | 海岸线变迁对湿地地下水影响的测量方法及其处理装置 | |
| CN112033641A (zh) | 一种海岸地貌演变物理模型试验方法 | |
| CN106320256B (zh) | 一种多沙河流抽水蓄能电站过机泥沙通量的确定方法 | |
| CN102087677A (zh) | 一种弹性地基梁的求解和m值反分析方法 | |
| CN109063230A (zh) | 不同突水或者透水条件下矿井淹没过程的评价方法 | |
| Zeelie | OMDEL dam and recharge ponds to enhance recharge in the Namib Desert | |
| CN106547974A (zh) | 一种采煤沉陷区蓄洪滞洪效果的定量方法 | |
| CN105045963A (zh) | 基于局部网格自适应技术的岸滩侧蚀数值模拟方法 | |
| CN117273263A (zh) | 一种基于dem数据工程抽水蓄能选址的方法 | |
| CN116151421A (zh) | 抽水蓄能电站选址方法、数据处理终端及可读存储介质 | |
| CN111877257B (zh) | 一种基于大数据处理的水坝建造方法 | |
| CN115630507A (zh) | 侵蚀基面下降后复式河道河床下切突变响应临界条件预测方法 | |
| KR101650251B1 (ko) | 3차원 수치해석을 통한 부유물 차단시설 입지 선정 방법 및 시스템 | |
| Corti et al. | Evaluation of the applicability of sediment transport models to dam filling prediction in different Italian geological contexts | |
| CN115994396A (zh) | 一种水系连通数学模型的构建方法 | |
| Mitew-Czajewska | Numerical analysis of a 36 m deep diaphragm wall–parametric study | |
| CN106934156A (zh) | 一种沿河道建筑危险点识别方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240118 Address after: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province Patentee after: Dragon totem Technology (Hefei) Co.,Ltd. Address before: 264005, Qingquan Road, Laishan District, Shandong, Yantai, 32 Patentee before: Yantai University |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240515 Address after: 116000, No. 5 Shazhongtun, Jinfang Village, Chengshan Town, Zhuanghe City, Dalian City, Liaoning Province Patentee after: Zhu Qingtao Country or region after: China Address before: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province Patentee before: Dragon totem Technology (Hefei) Co.,Ltd. Country or region before: China |
|
| TR01 | Transfer of patent right |