考虑交叉管线影响的城市供水管网地震响应分析及震害评估

钟紫蓝; 张亚波; 侯本伟; 韩俊艳; 杜修力

doi:10.11779/CJGE20220201

考虑交叉管线影响的城市供水管网地震响应分析及震害评估 English Version

北京工业大学城市与工程安全减灾教育部重点实验室, 北京 100124

基金项目:

国家自然科学基金面上项目 51978020

北京市教委科技计划一般项目 KM201910005022

国家重点研发计划项目 2018YFC1504305

详细信息

作者简介:
钟紫蓝(1986—)，男，博士，副研究员，主要从事生命线地震工程方面的研究工作。E-mail: zilanzhong@126.com

通讯作者:
侯本伟, E-mail: benweihou@bjut.edu.cn

中图分类号: TU990.3
计量
- 文章访问数: 225
- HTML全文浏览量: 34
- PDF下载量: 63
出版历程
- 收稿日期: 2022-02-24
- 网络出版日期: 2023-05-18
- 刊出日期: 2023-04-30

Seismic response analysis and damage assessment of urban water supply networks considering influences of crossing pipelines

Key Laboratory of Urban Security and Disaster Engineering of Ministry of Education, Beijing University of Technology, Beijing 100124, China

摘要

摘要: 以北京市某地区埋地供水管网为实例，通过直线型和交叉型典型管道二维有限元模型，系统分析了管道规格、接口形式、场地条件、地震动强度等级和地震动入射角等关键参数对管道接口轴向和转角变形的影响，得到了不同设防烈度地震动作用下该供水管网的损伤情况。基于国内外管道接口试验数据统计分析结果，确定了基于接口变形量的管道损伤判定准则，并用于不同管道接口类型地震损伤评价，建立了不同场地条件下典型管道接口地震损伤数据库。根据管道属性、场地信息和管道接口地震损伤数据库，基于GIS软件绘制管网地震损伤分布图。结果表明：相同烈度地震作用下，管线交叉处接口峰值变形量约为直管线接口的1.5倍~2.0倍，其中与法兰接口邻接的承插式接口存在峰值变形突变。罕遇烈度地震动作用下的管网地震损伤程度远高于设防烈度，且破坏多集中于Ⅳ类场地和管线交叉处。
- 供水管网 /
- 交叉管线 /
- 接口变形 /
- 有限元分析 /
- 震害评价
Abstract: Based on a buried water supply network in Beijing, the two-dimensional finite element models for the network are developed in this study. The influences of the critical parameters, such as the pipe diameter, joint type, site condition, intensity level of ground motion and incident angle of seismic wave, on the axial and bending deformations of pipe joints are systematically investigated, and the seismic damage status of the water supply network under different intensity levels of earthquakes is evaluated. Moreover, the criteria for damage assessment of the pipelines based on joint deformation are developed through the statistical analysis of the test results of the worldwide pipeline joints. These criteria are subsequently used for the seismic damage assessment of different types of pipeline joints. A seismic damage database of typical pipeline joints buried in different engineering sites is established. Finally, according to the pipeline properties, engineering site conditions and seismic damage database of typical pipeline joints, the seismic damage distribution maps of water supply networks are developed using the GIS. It is found that the peak deformations of the joints at the pipeline cross junctions are about 1.5 to 2.0 times those of the joints in a straight pipeline under the same intensity of earthquake ground motions. Besides, sudden changes of the peak seismic deformations occur at the push-on joints adjacent to the flange joints. The pipeline network suffers much more severe seismic damage under the considered maximum earthquake than under the design level of earthquake. The seismic damage mainly concentrates in the site class Ⅳ and the cross junction of the pipelines.
- water supply network /
- crossing pipeline /
- joint deformation /
- finite element analysis /
- seismic damage assessment

HTML全文

0. 引言

近年来,盾构施工技术已经在地铁隧道中得到了广泛的应用,如何正确预测城市地下工程施工对周围地层变形尤其是地表位移显得异常重要^[1-2]。纵观目前的地下施工技术,理论研究明显不足,实际施工大多数依赖于工程经验,因此,许多学者对盾构隧道施工引起的地层变形进行了研究。Peck^[3]于1969年提出了一种经验公式来预测地表沉降变形,许多工程实例已经证明了它的使用效果,但该公式中参数的选取往往取决于实际工程经验,未考虑土层条件,也不能预测深部土层位移情况；吴昌胜等^[4]研究了不同盾构直径对隧道地层损失率和地表沉降的影响；张冬梅等^[5]在研究隧道地表沉降时考虑了衬砌局部渗流作用的影响。以上研究均没有考虑土体的小应变特性,而大量工程实测数据表明,实际工程诸如地铁、基坑变形等中岩土体应变都很小,主要集中在0.01%～0.3%^[6-7],在此情况下,土体的剪切刚度是随着应变的增加而减小的,土体在小应变阶段的刚度要远大于较大应变阶段的刚度。众多的岩土试验^[8-9]也表明,土体的剪切模量与剪应变之间存在一种非线性关系,图1是相应的归一化刚度退化曲线。

图 1 刚度衰减曲线和应变范围

Figure 1. Curves of stiffness decreasing with strain and ranges of strain

下载: 全尺寸图片幻灯片

Lee等^[10-11]的研究表明,利用理想弹塑性本构模型（莫尔-库仑模型）来预测实际盾构施工变形,其结果偏于保守。因此,土体小应变特性会对预测盾构施工地层尤其是地表变形产生重要的影响。

鉴于此,本文从小应变塑性硬化本构（PH-SS模型）^[12]和莫尔-库仑本构（MC模型）出发,对盾构施工地层变形规律展开了研究。首先基于三轴模拟试验证实了PH-SS模型反映土体小应变特征的有效性；依托上海某一软土隧道工程,比较PH-SS模型和MC模型在模拟盾构施工中地表及深部土体变形规律的异同；最后结合现场监测数据,说明PH-SS模型在盾构施工地表变形预测中的先进性和合理性。

1. PH-SS模型的介绍

1.1 有限差分程序中的小应变本构模型

针对有限差分程序中的小应变特性研究,尹骥^[13]在室内试验的基础上,利用FLAC^3D中内嵌的fish语言将试验结果应用到有限差分程序的二次开发中,但这往往难度较大,同时自定义的本构模型往往受到试验结果较大的影响。本文利用FLAC^3D内置的Plastic Hardening-Small Strain model（简称PH-SS模型）^[12]来反映土体的小应变特性。

PH-SS模型以硬化土模型（PH模型）为基础,充分考虑了BENZ^[8]提出的卸载—再加载应变相关刚度关系,将土体的小应变特性考虑到有限差分程序中,最终用小应变参数G₀（初始剪切模量）和γ₇₀（剪应变阈值）来反应土体的小应变特性。除此之外,PH-SS还包含了11个硬化土模型参数：有效黏聚力c′、膨胀角ψ、有效内摩擦角 $φ^{'}$ 、三轴固结排水剪切试验的参考割线模量 $E_{50}^{ref}$ 、固结试验的参考切线模量 $E_{oed}^{ref}$ 、三轴固结排水卸载再加载试验的参考加卸载模量 $E_{ur}^{ref}$ 、与模量应力水平相关的幂指数m、参考应力p^ref、破坏比R_f、泊松比 $ν$ 、正常固结时的静止侧压力系数K₀。

PH模型以增量弹性定律来描述弹性行为；采用的是剪切屈服流动准则,确定剪切硬化开始和发展的剪切屈服函数可以表示为

$f^{s} = \frac{E_{ur}}{E_{i}} \cdot \frac{q_{a} q}{(q_{a} - q)} - q - \frac{E_{ur} γ^{p}}{2} = 0$ 。

(1)

式中 E_i=2E₅₀/(2-R_f)是初始变形模量； $q_{a}$ = $q_{f}$ /R_f,破坏比R_f是一个比1还小的值（通常取0.9）；而最终偏应力q_f的定义和莫尔-库仑准则是一致的,q是偏应力,PH模型还引入一种截止因子f_cut（Cut-off）以防止偏应力过大； $γ^{p}$ 是剪切硬化参数（内部变量之一）, $Δ γ^{p} = - (Δ ε_{1}^{p} + Δ ε_{2}^{p} + Δ ε_{3}^{p})$ ,由此可见, $γ^{p}$ 与3个塑性应变有关。

PH模型中的体积塑性应变和剪切塑性应变之间的流动规律如下：

$Δ ε_{v}^{p} = \sin ψ_{m} Δ γ^{p}$ 。

(2)

式中, $\sin ψ_{m}$ 为动力膨胀角,应该小于或等于定义的最终膨胀角ψ。

PH模型的剪切势函数为

$g^{s} = (- 1 + sin ψ_{m}) + (1 + \sin ψ_{m}) σ_{3} / 2$ 。

(3)

式中, $σ_{1}$ 为最大主应力, $σ_{3}$ 为最小主应力。

为了避免土体达到临界孔隙状态e_max时发生过度膨胀,对膨胀角进行了稍微修改。主要是采用Schanz等^[14]提出的一种设定分割线方法：

$\begin{array}{l} \sin ψ_{m}^{cut} = 0 & (e \geq) e_{max} \end{array}$ 。

(4)

这时为避免膨胀角的突然变化,对膨胀角进行平滑处理：

$\sin ψ_{m}^{smooth} = 100 (1 - e / e_{\max}) \sin ψ_{m} (e \geq 0.99 e_{\max})$ 。

(5)

对有截止因子f_cut（Cut-off）,光滑处理（Smoothing）和没有截止因子f_cut（No Cut-off）3种情况下的土体膨胀规律进行了如图2所示的对比,可以看出,光滑处理（Smoothing）情况下的数值计算最为稳定（渐变过程）,与实际情况也最为贴近。

图 2 标准三轴压缩试验的体积应变曲线

Figure 2. Curves volumetric strain of standard triaxial compression tests

下载: 全尺寸图片幻灯片

PH模型的体积屈服函数定义为

$f^{v} = g^{v} = {\tilde{q}}^{2} / a^{2} + p^{2} - p_{c}^{2} = 0$ 。

(6)

式中 $α$ 为一个常量,是基于数值固结仪测试推导得到的； $\tilde{q}$ 为一个剪应力量度；p为硬化参数,而硬化参数初值 $p_{c}^{}$ （代表是先前固结压力）可以由初始应力状态和超固结比OCR来确定。此处的体积硬化采用相关联流动法则。

PH-SS模型中的小应变特性是由小应变参数 $G_{0}$ 和 $γ_{70}$ 来反映的,其相应的关系式为

$G / G_{0} = {(1 + 0.385 γ / γ_{70})}^{^{- 2}}$ ,

(7)

式中 $γ_{70}$ 为剪切模量衰减到初始剪切模量的70%所对应的剪应变；初始剪切模量G₀可以由参考压力p^ref下的模量 $G_{0}^{ref}$ 在极小的应变水平下得

$G_{0} = G_{0}^{ref} {(\frac{c^{'} \cos ϕ^{'} + {σ^{'}}_{3} \sin ϕ^{″}}{c^{'} \cos ϕ^{'} + p^{ref} \sin ϕ^{'}})}^{m}$ 。

(8)

考虑到土体的实际情况,有限差分程序在PH-SS模型中充分考虑了卸载剪切模量 $G_{u r}$ 和卸载弹性模量 $E_{u r}$ 之间的关系,并设置了一个下限值：

$G_{u r} = E_{u r} / [3 (1 - 2 v)]$ 。

(9)

1.2 基于三轴模拟试验的PH-SS模型分析

利用FLAC^3D构建三轴压缩数值计算模型,模型尺寸为1 m×1 m×1 m,围压为100 Pa。为说明PH-SS模型在反映软土加载和卸载过程的优越性,该三轴模拟试验利用不同的本构模型（MC模型,PH模型和PH-SS模型）进行计算,且每种情况下都采用相同的循环荷载过程。

图3为三轴试样偏应力随轴向应变的变化曲线。可以看出,MC模型的卸载模量和加载模量是相同的；PH模型能够反映加载和卸载模量的差异性,但是不区分加载路径和卸载路径；PH-SS模型不仅能够区分加载和卸载模量的差异,而且能较好地区分加载路径和卸载路径,这种情况最符合相关试验^[15]的结果,同时也和谢东武等^[16]的模拟结果保持一致。PH-SS模型中出现回滞环进一步说明,随着剪应变的增大,土体的剪切模量迅速衰减,在加载—卸载的过程中产生塑性积累。

图 3 三轴试验中加卸载的应力应变曲线

Figure 3. Stress-strain curves of loading and unloading in triaxial tests

下载: 全尺寸图片幻灯片

图4为三轴试样的剪切模量随剪应变的变化曲线。由此可知：MC模型和PH模型中的模量始终保持恒定；而PH-SS模型能够较好地反应土体模量随剪应变增大逐渐衰减的特性,并达到一个下限值,这与式（9）是保持一致的,说明了PH-SS模型能够很好地表现土体的小应变特性。

图 4 剪切模量随剪应变的关系曲线

Figure 4. Curve of shear modulus with shear strain

下载: 全尺寸图片幻灯片

2. 基于不同模型的盾构隧道地层变形分析

2.1 计算模型的建立

根据上海某隧道^[17]的实际工程情况,建立二维数值模型。其中,隧道埋深为11.9 m,隧道直径为D为6.2 m。为避免边界的影响,X两侧方向取6倍的直径尺寸,下底边界取3倍的直径尺寸,即模型尺寸为76 m×40 m（宽×高）。隧道周围取1倍的隧道半径进行加密处理,单元尺寸为0.5 m,其他单元尺寸为1 m,数值模型共计有11942个节点、5880个单元。

本文分别利用PH-SS模型与MC模型对土体的力学行为进行模拟。参照文献[17～20],确定上海地区典型土层的小应变模型参数见表1。为保证数值计算结果更具有对比性,MC模型参数选取与PH-SS模型中相同的数值（其中,弹性模量取对应土层的加卸载模量）。衬砌材料利用shell单元来模拟,衬砌结构为线弹性材料, $γ$ =24.5 kN·m^-3, $ν$ =0.2；E= 24.44×10³ MPa,衬砌厚度d=0.35m。

表 1 PH-SS模型的土体参数^[17-20]

Table 1. Parameters of soil of PH-SS model

介质	$γ$ /(kN·m^-3)	$ν$	φ/(°)	c/kPa	/MPa	$E_{ur}^{ref}$ /MPa	/MPa	$p_{}^{ref}$ /kPa	m	$γ_{70}$ /10^-4
1素填土	18.0	0.35	18.0	10.0	4.0	23.2	92.8	100	0.9	2
2粉质黏土	18.7	0.35	15.0	18.0	6.0	34.8	139.2	100	0.9	2
3淤泥质粉质黏土	18.0	0.35	24.0	10.0	5.0	29.0	116.0	100	0.9	2
4淤泥质黏土	17.0	0.35	19.5	10.0	4.0	23.2	92.8	100	0.9	2
5-1粉质黏土	18.3	0.35	18.5	12.0	5.0	29.0	116.0	100	0.9	2
5-2粉质黏土	18.3	0.35	21.0	15.0	6.5	37.7	150.8	100	0.9	2
6粉质黏土	19.9	0.35	17.0	30.0	6.4	37.1	148.4	100	0.9	2
7-粉质黏土夹粉砂	20.1	0.30	26.3	6.0	10.2	59.2	236.6	100	0.9	2

下载: 导出CSV

| 显示表格

2.2 计算方案

采用应力释放法模拟盾构隧道的开挖和支护过程^[21]。利用不同的应力释放系数（ $λ$ 为0.1,0.2,0.3,0.4）反映隧道周围的地应力因开挖所引起的扰动程度。本文分别采用PH-SS和MC两种本构模型,总计八组计算工况,如表2。

表 2 计算工况设计

Table 2. Design of numerical cases

工况名称	是否考虑小应变特性	本构模型	应力释放系数λ
MC-0.1	否	MC	0.1
MC-0.2			0.2
MC-0.3			0.3
MC-0.4			0.4
PH-SS-0.1	是	PH-SS	0.1
PH-SS-0.2			0.2
PH-SS-0.3			0.3
PH-SS-0.4			0.4

下载: 导出CSV

| 显示表格

为研究深部土体的水平位移情况,取5个剖面（编号为^#1,^#2,^#3,^#4和^#5,相对位置分别如图5所示）为研究对象。对于地层竖向位移情况,本文重点关注地表、地表以下5 m处、地表以下10 m处、以及隧道以下5 m处的地层隆沉规律。

图 5 深层位移分析剖面位置示意图

Figure 5. Profile locations for displacement analysis of deep soils layers

下载: 全尺寸图片幻灯片

2.3 地表位移分析

图6给出了利用PH-SS模型计算得到的地表位移曲线,图7是利用MC模型计算得到的地表位移曲线。

图 6 利用PH-SS模型计算得到的地表位移曲线

Figure 6. Curves of surface displacement based on PH-SS model

下载: 全尺寸图片幻灯片

图 7 利用MC模型计算得到的地表位移曲线

Figure 7. Curves of surface displacement based on MC model

下载: 全尺寸图片幻灯片

由图6,7可知,两种地表沉降规律大致相似,但是在较大应力释放系数情况下（如 $λ$ =0.4）,PH-SS模型下的最大地表位移明显比MC模型下的地表位移大很多；在同一种本构模型计算工况下,应力释放系数越大,地表变形（竖向沉降、水平位移）越大,这是因为应力释放系数较大时,地层的应力扰动较大,造成地表变形较大。由图8可以得到,当应力释放系数小于0.3时,采用MC模型和PH-SS模型计算得到的地层损失率和最大地表沉降差别相对较小；但是当应力释放系数大于0.3时,采用PH-SS模型比MC模型计算的结果大得多,此时考虑模型之间的差异是必要的,这是因为盾构开挖应力释放系数较小（地层扰动较小）时,地层变形的小应变行为表现得不明显,当应力释放系数较大（地层扰动较大）时,小应变特征则十分明显,此时因剪应变而引起的刚度变化则十分显著。

图 8 地表最大沉降值、地层损失率与应力释放系数之间的关系

Figure 8. Relationship among maximum surface settlement, volume loss rate and stress release coefficient

下载: 全尺寸图片幻灯片

2.4 深部土层位移分析

以应力释放系数为0.3为例进行两种本构模型下的深部土层位移分析。图9分别为地表以下5 m、地表以下10 m的土层沉降的对比曲线,可以得到,深部土层位移场中的考虑小应变特性的沉降槽宽度比不考虑小应变特性的情况要小一些。图10是剖面^#1,^#2,^#4,^#5剖面的水平位移曲线,鉴于^#3剖面处于中央的位置,水平位移的差异性不明显,所以不单独列出^#3剖面的曲线。对图10的结果进行分析,沿着隧道对称的剖面（^#1和^#4,^#2和^#5）的水平位移变化规律沿隧道轴线对称,当埋深较浅时,由于地层沉降槽的形成产生了趋向隧道轴线方向的水平位移,而当埋深达到开挖面附近时,由于盾构管片的挤压作用,水平位移变成正值（土体沿着坐标轴正向移动）,随着分析剖面与隧道轴向距离的增加,岩土体受开挖扰动的影响逐渐减小。由此可见,盾构推进引起隧道周围土体向开挖面临空移动,最终水平位移大致关于中轴线对称。同时,考虑到隧道开挖造成周围地层在不同方向都产生了卸载作用,因而考虑土体加载模量和卸载模量的差异性时,土体模量变化较大,PH-SS模型相较MC模型在隧道轴线埋深处的地层水平位移更大。

图 9 地表以下地层剖面的沉降曲线

Figure 9. Settlement curves of soil layers below surface

下载: 全尺寸图片幻灯片

图 10 深部土体剖面的水平位移曲线

Figure 10. Curves of horizontal displacement of deep soil profile

下载: 全尺寸图片幻灯片

2.5 隧道以下地层隆起分析

盾构开挖会引起隧道以下地层发生卸载,PH-SS模型能区分土体加载和卸载模量的差异性。图11给出了不同应力释放系数情况下的隧道以下地层的位移规律。可以看出,两种本构模型下的地层位移规律差异较大。图12给出了对应地层处的最大隆起量（隆起为正、沉降为负）变化规律。可以看出,PH-SS模型计算隆起量比MC模型计算值较小,并且当应力释放系数增大时,其差异性逐步增大,说明应力释放系数越大（地层扰动越大）,应力相关模量的力学响应越明显,这也与前述保持一致,说明了PH-SS模型在地层隆起研究中的适宜性。

图 11 不同应力释放系数下的隧道以下5 m地层处的隆沉规律

Figure 11. Heave or settlement curves of layer 5 m below tunnel under different stress release coefficients

下载: 全尺寸图片幻灯片

图 12 地层隆起量随应力释放系数的变化曲线

Figure 12. Curves of soil layer uplift with stress release coefficient

下载: 全尺寸图片幻灯片

3. 与实测数据对比

根据上述研究,地层变形受到应力释放系数的影响较大。而实际工程情况受很多因素的影响,本文从简单问题出发,以地层损失为衡量标准^[22],通过三维数值计算结果反演二维模型的应力释放系数。三维模拟按照现场工况进行,每步推进1个管片环宽即1.0 m,利用强度不同的等代层来模拟注浆层及超挖层部分,最终确定二维数值计算的应力释放系数为0.38。因此,将 $λ$ =0.38情况下基于PH-SS模型和MC模型的数值计算结果与文献[17]中的监测数据进行对比,如图13所示。

图 13 数值计算结果与实测数据对比

Figure 13. Comparison between numerical results and measured data

下载: 全尺寸图片幻灯片

由此可知,利用PH-SS模型计算的地表沉降值更大,其地表沉降曲线相较MC模型更加符合实际监测结果,因此验证了预测盾构隧道地层变形中考虑小应变特性的必要性,同时也说明了PH-SS模型在预测地表沉降变形中的适用性。

4. 结论

软土地区盾构施工中的大部分区域处于小应变的状态,本文基于FLAC^3D内置的MC模型和PH-SS模型对盾构施工地层变形规律展开研究。首先基于三轴模拟试验证明了PH-SS在表现土体小应变特性方面的合理性和优越性；依托上海某软土隧道工程,分别利用MC模型和PH-SS模型进行二维盾构施工数值计算研究,着重研究了地层变形尤其是地表沉降随应力释放系数的演化规律,以及利用两种本构模型计算所得结果之间的差异性；最后将数值结果与现场数据进行了对比。

（1）相较MC模型,PH-SS模型不仅可以区分软土加载和卸载模量的差异性,而且可以较好地反映土体剪切模量随剪应变增大而逐渐衰减的特性,在表现软土行为时更为合理。

（2）当应力释放系数小于0.3时,预测地表沉降时考虑小应变特性与不考虑时的情况区别不大；但是当应力释放系数大于0.3时,考虑土体小应变特性是十分必要的。

（3）应力释放系数越大,土体应力相关模量的力学响应越明显,地层变形越大。

（4）数值计算结果与实测数据进行对比的结果表明,考虑小应变特性的PH-SS模型地表沉降更大,与实测数据更为贴近,说明了PH-SS模型在预测盾构隧道地表变形规律更为合理。

图 1 基于物理模型的城市供水管网系统地震响应分析方法

Figure 1. Seismic response analysis method for urban water supply networks based on physical model

下载: 全尺寸图片幻灯片

图 2 交叉管线数值模型

Figure 2. Numerical model for crossing pipelines

下载: 全尺寸图片幻灯片

图 3 人工合成地震动加速度时程

Figure 3. Time histories of artificial ground motions

下载: 全尺寸图片幻灯片

图 4 地震动输入点坐标转换

Figure 4. Coordinate transformation of input ground motions

下载: 全尺寸图片幻灯片

图 5 DN200支干线接口地震峰值响应（Ⅲ类场地、E2地震）

Figure 5. Peak seismic responses of DN200 joints in main and branch pipelines (Site class Ⅲ, intensity E2)

下载: 全尺寸图片幻灯片

图 6 不同平面入射角

Figure 6. Schematic diagram of different plane incident angles

下载: 全尺寸图片幻灯片

图 7 不同入射角下L型管线接口响应（Ⅲ类场地，E2地震）

Figure 7. Responses of L-type pipeline joints at different incident angles (Site class Ⅲ, intensity E2)

下载: 全尺寸图片幻灯片

图 8 不同入射角下十字型管线接口响应（Ⅲ类场地，E2地震）

Figure 8. Responses of cross-type pipeline joints at different incident angles (Site class Ⅲ, intensity E2)

下载: 全尺寸图片幻灯片

图 9 两类场地下管线接口峰值变形

Figure 9. Peak deformations of pipeline joints at two engineering sites

下载: 全尺寸图片幻灯片

图 10 不同管径下的管道接口响应影响因子

Figure 10. Influence factors of response of pipeline joints under different pipe diameters

下载: 全尺寸图片幻灯片

图 11 供水管网地震损伤分布

Figure 11. Damage distribution map for water supply networks

下载: 全尺寸图片幻灯片

图 12 管网受损长度分布

Figure 12. Distribution of totally damaged length for pipeline networks

下载: 全尺寸图片幻灯片

表 1 研究区域代表性地层参数信息

Table 1 Representational parameters of soil layers in target region

场地编号	土层种类	厚度 H/m	密度 ρ/(kg·m^-3)	弹性模量 E/MPa	剪切波速 V_s/(m·s^-1)	场地类别
①	杂填土	3	1700	4	100	Ⅲ类
	粉土	4	1850	12	75
	中砂、细砂	5	1900	32	200
	粉质黏土	6	1900	70	375
②	碎石填土	4	1800	10	80	Ⅳ类
	细砂、粉砂	7	2300	56	180
	细砂、中砂	4	2400	80	200
	砂质粉土	7	1400	20	360

下载: 导出CSV

表 2 球墨铸铁管道接口模型参数^[15-24]

Table 2 Parameters of DI pipe joint model

管径	峰值轴力F_u/kN	接触转角θ₁/(°)	接触弯矩M₁/(kN·m)	极限转角θ₂/(°)	极限弯矩M₂/(kN·m)
150	9.2	7.0	1.46	13.1	14.83
200	11.6	6.6	2.86	12.5	7.60
300	17.4	5.8	3.57	—	—
400	21.0	3.9	6.00	10.5	10.00

下载: 导出CSV

表 3 球墨铸铁管道接口拉伸试验结果统计

Table 3 Statistics of tensile test results of DI pipe joints

样本编号	直径D/mm	管内水压P/MPa	接口渗漏张开量Δ_u/mm	安装深度d_p/mm	随机变量X (Δ_u/d_p)	参考文献
1	100~250	—	30.5	—	—	Singhal等^[15]
2	200	0.4	41~49.16	100	0.41~0.492	刘为民等^[16]
3	150	0.2	70~80	94	0.745	周静海等^[17]
4	200	0.2	65	100	0.65	周静海等^[17]
5	200	0.6	48	100	0.48	同济大学^[18-19]
6	200	0.2	51.7~66.98	100	0.517~0.67	同济大学^[18-19]
7	150	0.38	56	85.8	0.653	Wham等^[20]
8	200	1.0	60	100	0.6	王颂翔^[21]
9	150	0.2	52.6~53.1	100	0.526~0.531	李冠潮^[22]
10	200		55.3~56.2	100	0.553~0.562
11	300		51~54	100	0.51~0.54
12	150	0~0.3	54.15~55.77	94	0.576~0.593	李晓晓等^[23]，钟紫蓝等^[24]
13	200	0，0.1	60，58	100	0.6，0.58
14	400	0	50	110	0.455

下载: 导出CSV

表 4 研究区域管道规格

Table 4 Pipeline types in target region

管道编号	规格/mm
1	100×9×6000
2	200×12×6000
3	400×15×6000
4	600×15.8×6000
5	1000×18×6000
注：规格为内径×壁厚×长度。

下载: 导出CSV

表 5 交叉管线接口响应结果

Table 5 Seismic responses of cross-type pipeline joints

公称直径	承插式接口		法兰接口
公称直径	峰值张开量d_max/mm	峰值转角R_max/(°)	峰值压缩量 P_max/mm	峰值转角R_max/(°)
100	6.33	1.9×10^-2	0.872	2.20×10^-3
200	6.37	2.0×10^-2	0.611	0.90×10^-3
400	6.40	2.1×10^-2	0.518	0.26×10^-3
600	6.46	2.3×10^-2	0.494	0.18×10^-3
1000	6.54	3.8×10^-2	0.287	0.10×10^-3

下载: 导出CSV

表 6 管线接口损伤表

Table 6 Statistics of damage of pipeline joints

接口形式	场地类型	管径	多通形式	地震动等级	损伤分类
承插式接口	Ⅲ类	All	直线/L型	E1/E2/E3	基本完好
			T型/十字型/斜交型	E2	基本完好
			T型/十字型/斜交型	E3	中等破坏
	Ⅳ类	All	All	E1	基本完好
			直线/L型	E2	基本完好
			T型/十字型/斜交型	E2	中等破坏
			直线/L型	E3	中等破坏
		＜DN600	T型/十字型/斜交型		中等破坏
		≥DN600	T型/十字型/斜交型		严重破坏
法兰接口	Ⅲ类	All	L型	E1/E2/E3	基本完好
		DN100	T型/十字型/斜交型	E1	基本完好
		DN100		E2/E3	中等破坏
		DN200		E1/E2	基本完好
		DN200		E3	中等破坏
		≥DN400		E1/E2/E3	基本完好
	Ⅳ类	All	L型	E1/E2/E3	基本完好
		DN100	T型/十字型/斜交型	E1/E2/E3	中等破坏
		DN200		E1	基本完好
		DN200		E2/E3	中等破坏
		≥DN400		E1/E2/E3	基本完好

下载: 导出CSV

参考文献(31)

[1]	TOSHIMA T, IWAMOTO T, NAKAJIMA T. 2000. Study on behavior of buried pipes in liquefed ground[C]//12th World Conference on Earthquake Engineering, Auckland, 2000.
[2]	李乔, 赵世春. 汶川大地震工程震害分析[M]. 成都: 西南交通大学出版社, 2008. LI Qiao, ZHAO Shichun. Analysis of Seismic Damage of Engineering Structures in Wenchuan Earthquake[M]. Chengdu: Southwest Jiaotong University Press, 2008. (in Chinese)
[3]	WHAM B P, DASHTI S, FRANKE K, et al. Water supply damage caused by the 2016 Kumamoto Earthquake[J]. Lowland Technology International: The Official Journal of the International Association of Lowland Technology. 2017, 19(3): 151-160.
[4]	李杰. 生命线工程抗震: 基础理论与应用[M]. 北京: 科学出版社, 2005. LI Jie. Structural Response Network Disaster Simulation Lifeline Ground Motion Reliability[M]. Beijing: Science Press, 2005. (in Chinese)
[5]	O'ROURKE T D, JEON S S, TOPRAK S, et al. Earthquake response of underground pipeline networks in Christchurch, NZ[J]. Earthquake Spectra, 2014, 30(1): 183-204. doi: 10.1193/030413EQS062M
[6]	O'ROURKE M J, LIU X. Response of Buried Pipelines Subjected to Earthquake Eﬀects[R]. La Mirada: University of Buﬀalo, 1999.
[7]	SHIROZU T, YUNE S, ISOYAMA R, et al. Report on damage to water distribution pipes caused by the 1995 Hyogoken-Nanbu (Kobe) earthquake[J]. Terremotos, 1996: 93-110.
[8]	冯启民, 高惠瑛, 俞虹桥. 供水系统震害预测专家系统[J]. 地震工程与工程振动, 2000, 20(3): 67-75. doi: 10.3969/j.issn.1000-1301.2000.03.010 FENG Qimin, GAO Huiying, YU Hongqiao. An expert system for predicting damage to water delivery systems[J]. Earthquake Engineering and Engineering Vibration, 2000, 20(3): 67-75. (in Chinese) doi: 10.3969/j.issn.1000-1301.2000.03.010
[9]	冯启民, 高惠瑛, 张伟林. 天津开发区供水系统地震反应分析[J]. 自然灾害学报, 2000, 9(4): 39-44. doi: 10.3969/j.issn.1004-4574.2000.04.007 FENG Qimin, GAO Huiying, ZHANG Weilin. Seismic performance analysis of water delivery systems in Tianjin Economic Developed Area[J]. Journal of Natural Disasters, 2000, 9(4): 39-44. (in Chinese) doi: 10.3969/j.issn.1004-4574.2000.04.007
[10]	NOURZADEH D, MORTAZAVI P, GHALANDARZADEH A, et al. Numerical, experimental and fragility analysis of urban lifelines under seismic wave propagation: study on gas distribution pipelines in the greater Tehran area[J]. Tunnelling and Underground Space Technology. 2020, 106: 103607. doi: 10.1016/j.tust.2020.103607
[11]	MAKHOUL N, NAVARRO C, LEE J S, et al. A comparative study of buried pipeline fragilities using the seismic damage to the Byblos wastewater network[J]. International Journal of Disaster Risk Reduction. 2020, 51: 101775. doi: 10.1016/j.ijdrr.2020.101775
[12]	中华人民共和国建设部. 室外给水排水和燃气热力工程抗震设计规范: GB 50032—2003[S]. 北京: 中国标准出版社, 2003. Ministry of Construction of the People's Republic of China. Code for Seismic Design of Outdoor Water Supply Sewerage Gas and Heating Engineering: GB 50032—2003[S]. Beijing: Standards Press of China, 2003. (in Chinese)
[13]	Ductile Iron Pipes, Fittings, Accessories and Their Joints for Water Applications: ISO 2531: 2009[S]. 2009.
[14]	Japan Water Works Association JWWA. Seismic Design Specifications for Water Facilities – in Japanese[S]. 1997.
[15]	国家质量监督检验检疫总局, 中国国家标准化管理委员会. 压力容器第3部分: 设计: GB/T 150.3—2011[S]. 北京: 中国标准出版社, 2012. General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, Standardization Administration of the People's Republic of China. Pressure Vessels—Part 3: Design: GB/T 150.3—2011[S]. Beijing: Standards Press of China, 2012. (in Chinese)
[16]	American Lifelines Alliance (ALA). Seismic Guidelines for Water Pipelines[M]. ASCE, 2005.
[17]	O'ROURKE M, VARGAS-LONDONO T. Analytical model for segmented pipe response to tensile ground strain[J]. Earthquake Spectra, 2016, 32(4): 2533-2548. doi: 10.1193/050415EQS064M
[18]	洪华生. 工程中的概率概念[M]. 北京: 中国建筑工业出版社, 2017. HONG Huasheng. Probability Concepts in Engineering[M]. Beijing: China Architecture & Building Press, 2017. (in Chinese)
[19]	SINGHAL A C, BENAVIDES J C. Axial and bending behavior of buried pipelin e joints[C]// Proceedings of the Fourth National Congress on Pressure Vessel and Piping Technology, American Society of Mechanical Engineers. 1983.
[20]	刘为民, 孙绍平. 管道接口的抗震试验研究[C]//第五届全国地震工程会议. 北京, 1998. LIU Wei-min, SUN Shao-ping. Shaking table test of pipe with junctions[C]// Proceeding of Fifth China Symposia on Earthquake Engineering. Beijing, 1998. (in Chinese)
[21]	周静海, 赵海艳, 魏立群. 球墨铸铁供水管线在地震作用下功能性实验分析[J]. 沈阳建筑大学学报(自然科学版), 2008, 24(2): 196-199. https://www.cnki.com.cn/Article/CJFDTOTAL-SYJZ200802005.htm ZHOU Jinghai, ZHAO Haiyan, WEI Liqun. Experimental research on functionality of ductile cast iron pipelines under the earthquake[J]. Journal of Shenyang Jianzhu University (Natural Science), 2008, 24(2): 196-199. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SYJZ200802005.htm
[22]	刘威, 黄鹭娜, 李杰. 供水管线渗漏试验研究[J]. 地震工程与工程振动, 2011, 31(4): 167-173. LIU Wei, HUANG Luna, LI Jie. Experiment on leakage of water pipelines[J]. Journal of Earthquake Engineering and Engineering Vibration, 2011, 31(4): 167-173. (in Chinese)
[23]	傅俊. 供水管线接口变形试验与抗震可靠度研究[D]. 上海: 同济大学, 2013. FU Jun. Displacement Experiment and Seismic Reliability Study of Water Supply Pipeline Joint[D]. Shanghai: Tongji University, 2013. (in Chinese)
[24]	WHAM B P, DENIS O'ROURKE T. Jointed pipeline response to large ground deformation[J]. Journal of Pipeline Systems Engineering and Practice, 2016, 7(1): 4015001- 4015009.
[25]	王颂翔. 承插式给水管道安全评价研究[D]. 大连: 大连理工大学, 2015. WANG Songxiang. Study on Safety Evaluation of Water Supply Pipelines with the Socket and Spigot Joint[D]. Dalian: Dalian University of Technology, 2015. (in Chinese)
[26]	李冠潮. 管道柔性接口轴向力学性能及功能性试验研究[D]. 郑州: 河南工业大学, 2020. LI Guanchao. Experimental Study on The Axial Mechanical Properties and Functional Properties of Flexible Joints in Pipelines[D]. Zhengzhou: Henan University of Technology, 2020. (in Chinese)
[27]	李晓晓, 钟紫蓝, 侯本伟, 等. 大型球墨铸铁管承插式接口力学性能研究[J]. 特种结构, 2020, 37(4): 47-55. https://www.cnki.com.cn/Article/CJFDTOTAL-TZJG202004010.htm LI Xiaoxiao, ZHONG Zilan, HOU Benwei, et al. Study on mechanical behaviors of push-on joints of large-diameter ductile iron pipelines[J]. Special Structures, 2020, 37(4): 47-55. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TZJG202004010.htm
[28]	钟紫蓝, 王书锐, 杜修力, 等. 管道承插式接口轴向力学性能试验研究与数值模拟[J]. 工程力学, 2019, 36(3): 224-230, 239. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201903024.htm ZHONG Zilan, WANG Shurui, DU Xiuli, et al. Experimental and numerical study on axial mechanical properties of pipeline under pseudo-static loading[J]. Engineering Mechanics, 2019, 36(3): 224-230, 239. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201903024.htm
[29]	球墨铸铁给排水管道工程施工及验收规范技术要求: ZXB/T 0202—2013[S]. 2013. The Code for Construction and Acceptance of Water and Sewerage Ductile Iron Pipeline Works Specification Requirement: ZXB/T 0202—2013[S]. 2013. (in Chinese)
[30]	Rubber Seals-Joint Rings for Water: ISO 4633—2015[S]. 2015.
[31]	范重, 张康伟, 张郁山, 等. 视波速确定方法与行波效应研究[J]. 工程力学, 2021, 38(6): 47-61. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX202106005.htm FAN Zhong, ZHANG Kangwei, ZHANG Yushan, et al. Study on apparent wave velocity calculation method and on travelling wave effect[J]. Engineering Mechanics, 2021, 38(6): 47-61. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX202106005.htm