3.5 溃坝后果计算

溃坝后果是指水库溃坝造成的下游可能淹没范围内生命损失、经济损失、社会与环境影响等。在进行溃坝后果计算之前，应进行溃坝洪水分析，制作溃坝洪水风险图。

3.5.1 溃坝洪水分析

溃坝洪水分析包括溃口洪水分析和溃坝洪水演进计算。

3.5.1.1 溃口洪水分析

进行溃口洪水分析时，首先确定溃坝模式，然后进行溃口流量变化过程计算。土石坝和面板堆石坝宜采用逐步溃坝模式，重力坝和拱坝宜采用瞬时全溃或瞬时局部溃决模式。溃口洪水计算可参考文献［33］，下面对文献［33］进行补充完善。

1.逐步溃坝溃口洪水计算

（1）BREACH模型。BREACH模型是基于Fread（1984）预报土坝溃坝洪水过程线而开发的一个数学模型。该模型建立在水力学、泥沙输移、土力学、大坝几何尺寸与数学特征、库容特征、溢洪道特征及入库流量随时间变化的基础上。模型有7个主要部分：①溃口形成；②溃口宽度；③库水位；④溃口泄槽水力学；⑤泥沙输移；⑥突然坍塌引起溃口的扩大；⑦计算方法。模型可以模拟因漫顶或管涌引起的溃坝，大坝可以是均质的，也可以是由两组不同特性材料组成的坝壳和心墙。

（2）面板堆石坝溃坝模型。面板堆石坝与一般土石坝溃决过程的区别在于钢筋混凝土面板在未被冲毁的下游坝体支撑下仍起挡水作用，而随着下游坝体冲刷的积累，面板悬空长度不足以承受面板自重和水荷载的共同作用便折断，此时水头突然增加，溃口处流量突增，溃决过程突然加速。其后，随着水头逐渐减小，溃口流速、流量及冲刷也逐渐减小，面板又起到挡水作用。如此往复，直至最后稳定在某一平衡位置。

1）李雷模型。为简化分析，设面板为单向板，取宽度b=1.0m，当下游支撑体被冲毁后，可视为一在自重和水荷载共同作用下的悬臂板，其计算简图见图3.6。

则自重荷载产生的弯矩M₁可按式（3.14）计算：

式中 m₁——上游坡比；

ρ_h——面板重度，kN/m³；

h——面板厚度，m。

水荷载产生的弯矩M₂可按式（3.15）计算：

式中 h₀——库水位，m；

Z_f——面板顶端高程，m；

L_d——面板折断长度，m；

ρ_w——水的重度，kN/m³。

图3.6 面板折断长度L_d推导简图

面板折断点总弯矩M按式（3.16）计算：

根据混凝土标号、面板配筋率和面板截面积可确定钢筋面积，再从钢筋混凝土结构设计手册查出面板所能承受的弯矩M。再考虑安全系数k，则面板所能承受的极限弯矩M_f可按式（3.17）计算：

故面板折断条件为

即

式（3.19）是一个一元三次方程，采用牛顿法可很快求得临界折断长度。

面板折断形成溃口后，如果库水无别的出处，也没有别的水量补给，则水库减少的水量就是溃口的出流量，则有

式中 A_s（h）——水库水面面积，m²；

dh/dt——水面高程的时间变化率，m/s；

Q——溃口流量，m³/s。

在Δt很小时，式（3.20）可改写为

溃口的出流量与平均流速u、过流面积A的乘积成正比，有

通常溃口出流量作为宽顶堰来计算的，即

对照式（3.22）、式（3.23），作为矩形断面，可知

式中 m——流量系数；

σ_c——侧收缩系数；

b——溃口平均宽度，m；

h ₀和Z_f如前所定义。

V.P.Singh等（1988）提出如式（3.25）所示的下游坝体冲刷速度公式：

式中 Z——下游坝体冲刷高程，m；

u——水流平均流速，m/s；

α ₂和β₂——经验系数，β₂大约为2，误差可用冲刷系数α₂修正，α₂=0.0008～0.0090。

当Δt取得很小时，式（3.25）可表示为

当溃口开始出流时，在起始条件明确时，即可根据式（3.24）、式（3.26）及式（3.18）确定在Δt₁时段内面板会否折断，从而确定Δt₂时段中面板顶部高程，再根据式（3.21）和式（3.22）求出剩余库容。根据库容曲线求得新的水面高程，再继续Δt₂时段的计算，循环往复，便可分析各时段水力要素（溃口流量、流速、下泄水量等）、上游库水位和面板顶高程，最后得到上述各要素的溃坝过程线。

2）陈生水模型。陈生水认为^［68］，要正确模拟面板砂砾石坝溃决过程，应重点解决以下两个问题：一是针对砂砾石材料粒径范围宽的特点，提出能较为合理计算其冲蚀率的经验表达式；二是正确分析下游坝体冲蚀量对混凝土面板受力状态的影响，合理确定面板的折断时刻。

a.砂砾石坝料临界起动流速。忽略砂砾石颗粒间的黏聚力，对如图3.7所示的代表颗粒1而言，所受的力一般有浮重度W，水流拖曳力F_d、上举力F_l，其表达式分别为

式中 γ_s——土颗粒的重度，kN/m³；

γ_w——水的重度，kN/m³；

d ₅₀——土体颗粒平均粒径，m。

式中 F_d——水流对土体颗粒的拖曳力，kN；

C_d——拖曳力系数，一般取0.4；

v——水流流速，m/s。

式中 F_l——水流对土颗粒的上举力，kN；

C_l——上举力系数，一般取0.1。

考虑到砂砾石料级配范围宽，最大颗粒与最小颗粒粒径相差大，为反映在漫坝水流作用下粗颗粒对细颗粒的阻拦、遮蔽作用，细颗粒对粗颗粒的包围、填实等，除上述作用力外，再引入一个与水流方向垂直的附加作用力R（见图3.7），并近似假定R与颗粒间的平均剪力成比例，即

图3.7 土体颗粒在坝坡上的受力示意图

式中 τ_s——不均匀颗粒的平均剪力，kPa；

K_m——无因次系数；

φ——比例系数与颗粒面积系数的乘积。

由此可得

于是式（3.31-1）可以写为

式中 K——可根据不均匀颗粒起动流速实验资料确定，约在0.785～1.727范围内变化，此处取1.3；

M——紧密系数，代表颗粒组成的密实程度，与不均匀系数C_u有关，可由式（3.32）确定。

如图3.7所示，土体颗粒起动时受到的摩擦力可表示为

式中 F_f——土颗粒受到的摩擦力，kN；

φ——土体颗粒间的内摩擦角；

θ——坝坡坡角；

c——土体的黏聚力，kPa。

通过受力分析可知，土体颗粒1起动的临界条件为

将式（3.27）、式（3.28）、式（3.29）、式（3.33）、式（3.34）代入式（3.35）可得到土体颗粒在坝坡上的临界起动流速为

b.砂砾石坝体冲蚀率。砂砾石坝体在溃坝水流的作用下，坝顶溃口和下游坝坡将发生冲蚀。针对砂砾石料级配范围宽，最大颗粒与最小颗粒粒径相差大，在分析不同土体陡水槽冲蚀试验结果的基础上，选择d₉₀与d₃₀作为代表粒径，建议计算砂砾石料单宽冲蚀率的经验表达式如下

式中 q_s——单宽冲蚀率，m²/s；

d ₉₀和d₃₀——小于某粒径的颗粒含量分别为90%和30%所对应的颗粒粒径，mm；

B——溃口宽度，m；

v _*——摩阻流速，m/s；

v_b——溃口底流速，m/s；

v——水流平均流速，m/s；

Q_b——溃口流量，m³/s；

J——水力梯度；

H——水库水位高程，m；

H_c——溃口底部高程，m；

N——溃口糙率。

溃口流量Q_b分别采用下述方法计算。

当漫顶溃坝发生后，水流沿着初始溃口冲蚀下游坝坡，可采用式（3.38）的宽顶堰公式计算：

式中 m——流量系数，此处取0.5。

随着下游坝料冲蚀，溃口向上游发展，溃口部位坝顶宽度逐渐减少。当坝顶冲蚀完毕后，由于面板还在发挥挡水作用，因此过流断面为薄壁堰，此时的溃口流量按式（3.39-1）计算：

式中 m₀——流量系数；

b——薄壁堰过流宽度，m；

Δh——堰顶水头，m；

p ₁——堰高，m。

c.溃口发展和面板折断时刻确定。由于溃坝水流的作用，下游坝坡发生冲蚀，某一时间段增量Δt_i内，溃口纵向冲蚀深度增量Δy_ci为

式中 Δy_c——溃口纵向冲蚀深度，m；

L ₂——下游坝坡长度，m；

n——坝壳料孔隙率。

假设初始溃口为梯形，在水流冲蚀情况下，溃口边坡保持极限稳定边坡，坡角为砂砾石材料的内摩擦角φ；溃口底部的冲蚀速率与溃口边坡的冲蚀速率相等，见图3.8，则水流对坝体溃口两侧的直接冲刷形成的溃口宽度增量ΔB可表达为

时间段Δt_i内水库水位变化量为

图3.8 面板砂砾石坝下游坝坡冲蚀示意图

式中 Q_in——入库流量，m³/s；

S_a——库水位为H时的水库面积，m²。

由于面板的挡水作用，溃口的发展主要以溢流水流对下游坝体的冲刷为主。随着坝体不断经受过坝水流冲刷，面板下方残存的堆石体下塌至极限稳定边坡（见图3.9），上游边坡坡比记为1∶m₁。取单位宽度面板进行分析，当下游支撑体被冲蚀后，面板可视为在自重和水荷载共同作用下的悬臂板，悬臂板承受自重荷载产生的弯矩M₁及水荷载产生的弯矩M₂可通过式（3.44）、式（3.45）计算。

式中 M₁——自重荷载产生的弯矩，kN·m；

ρ_m——面板的密度，kg/m³；

δ——面板的厚度，m；

L_d——折断面板的长度，m。

式中 M₂——水荷载产生的弯矩，kN·m；

ρ_w——水的密度，kg/m³；

H——库水位高程，m；

Z_f——面板顶部高程，m。

图3.9 面板砂砾石坝漫顶破坏溃口发展示意图

在自重荷载和水荷载的共同作用下面板承受的总弯矩M为

根据面板所用的混凝土标号、配筋率和面板的横截面积，查出面板所能承受的弯矩M₀，安全系数k，则面板的极限弯矩M_f为

其中k=1.4。

因此面板折断的条件为

则面板折断的临界条件为

则面板折断的临界长度可通过求解关于L_d的一元三次方程得到

坝体冲蚀量决定了面板折断时刻。当面板悬空长度为L_d时，残余坝体的单宽体积为

式中 L₁——面板的总长度，m。

坝体的单宽总体积可表示为

式中 W₁——坝顶宽度，m；

W ₂——坝底部宽度，m；

H ₀——坝高，m。

则面板折断的时刻在溃坝发生后的时间ΔT可按式（3.53）表达为

此时可通过式（3.42）、式（3.53）求出此时溃口的宽度和面板发生折断的时间。当面板折断后，溃口宽度和溃口流量将突然增大。如果面板砂砾石坝较高，面板可能不止发生一次折断，可继续采用上述方法计算面板发生折断的时间和溃口的发展。

（3）简化计算法。坝址处的最大溃坝流量Q_m可按式（3.54）计算：

式中 Q_m——最大溃坝流量，m³/s；

F——简化评估特征参数，F=1.3；

V——下泄水的总体积，10³m³；

H——最大水深，m。

式（3.54）为保守地利用地形与水文过程线资料和凭经验确定的溃决流量表达式，适用于坝高小于等于12m的均质土坝。对于坝高大于12m的均质土坝或坝高小于等于12m的非均质土石坝或非黏性材料坝，其溃口特性不同，最大溃坝流量应做相应调整。

2.瞬时溃坝溃口洪水计算

（1）混凝土坝瞬时全溃。假定坝下游无水，上下游河槽断面为矩形，槽底坡降i=0，并设溃坝时水流惯性力为主导，忽略水流阻力，则根据圣维南方程和特征线理论，溃坝波的波形为式（3.55）、式（3.56）所示的二次抛物线（如图3.10所示）方程：

图3.10 溃坝波的二次抛物线

式中 H₀——坝址上游水深，m。

当x=0时，坝址处的水深和流速即为常数，即

若矩形断面的宽度为B，则坝址处的最大流量可按式（3.59）计算：

若上下游河道断面不为矩形，设断面面积可表示为式（3.60）：

式中 K——常系数；

B——水面宽；

m——河槽断面形状系数。

则坝址溃口处流速和最大流量可分别按式（3.61）、式（3.62）计算：

此即所谓的圣维南公式解，也称A.Ritter解。

（2）混凝土坝横向局部一溃到底。如图3.11所示，再考虑溃口影响因子——（B／b）^α后，便可用瞬时全溃的公式来计算横向局部一溃到底时的坝址溃口处水深、流速及最大流量。若取α=0.25，则式（3.57）、式（3.58）及式（3.59）分别修正为

（3）混凝土坝瞬时垂向局部溃坝。与横向局部一溃到底不同的是，在坝高方向残留了一段坝体，如图3.12所示。

图3.11 混凝土坝横向局部一溃到底示意图

（箭头表示负波传播方向）

图3.12 混凝土坝垂向局部溃坝示意图

设坝体残留部分高度为h′，则此时溃口处的最大流量可按式（3.66）计算：

如果沿大坝横向和垂向同时局部溃决，则最大流量公式修改为式（3.67）：

美国水道实验站对式（3.67）进行了修正：

黄河水利委员会水利科学研究院根据试验提出式（3.69）：

（4）瞬时溃坝时坝址流量过程线估算。通过详细算法成果和模型实验资料分析，发现瞬时溃坝流量过程线的最大流量与溃坝前下泄流量及溃坝前可泄库容（溃坝库容）有关，可概化为4次或2.5次抛物线，见表3.19、表3.20，其中Q₀为溃坝前下泄流量，t为溃坝时刻，Q为t时刻溃坝流量。实际多用4次抛物线。

库水泄完时间T可按式（3.70）确定：

式中 W——库容；

——流量过程线纵坐标的平均值；

Q_m——溃口处最大流量；

Q ₀——溃坝前下泄流量。

3.5.1.2 溃坝洪水演进计算

溃坝洪水演进计算包括洪水向下游演进时的沿程洪水到达时间、流速、水深、历时等洪水要素的计算，山区、丘陵区的小型水库溃坝洪水演进计算可采用一维数学模型，其他水库应采用平面二维数学模型。目前用于溃坝洪水演进计算的软件很多，常用的如美国国家气象局（NWS）开发的DAMBRK、FLADWAV，美国陆军工程师团（USACE）开发的HEC-RAS，英国HRWallingford软件公司开发的InfoWorks软件，以及丹麦水利科学研究所（DHI）开发的MIKE系列。对于一些小型水库，若无计算软件，可采用简化分析法和经验公式法进行溃坝洪水演进计算，具体参见文献［33］。

3.5.2 溃坝洪水风险图制作

溃坝洪水风险图是融合洪水特征信息、地理信息、社会经济信息，通过洪水计算、风险判别、社会调查，反映溃坝发生后潜在风险区域洪水要素特性的专题地图。

溃坝洪水风险图制作的一般流程为：收集整编资料、确定计算范围和溃坝洪水风险分析方法、溃坝洪水风险分析、溃坝洪水风险图制作。

溃坝洪水风险图包括纸质溃坝洪水风险图、电子溃坝洪水风险图两种。纸质溃坝洪水风险图是在电子溃坝洪水风险图基础上，按照信息显示要求进行编辑加工后的打印输出，基本内容应与电子版溃坝洪水风险图保持一致。

溃坝洪水风险图可包括如下信息：工作底图信息、风险要素（洪水水深、流速、淹没历时、到达时间、严重性等）信息、防洪工程信息、防洪非工程信息、社会经济信息等。根据不同要求，信息可有所侧重。一般水库必需制作溃坝洪水淹没范围图，对大中型水库，还应制作溃坝洪水严重性分布图、溃坝洪水到达时间分布图。

（1）溃坝洪水淹没范围图制作。洪水淹没范围图制作及成果要求可参见水利行业标准《洪水风险图编制导则》（SL 483—2010）的规定。

（2）溃坝洪水严重性分布图制作。依据洪水演进计算成果，提取某一时刻计算单元或计算断面的水深和流速，生成淹没区的洪水严重性，并将其绘制到工作底图上，形成洪水严重性分布图。

（3）溃坝洪水到达时间分布图制作。依据洪水演进计算成果，提取淹没区的洪水到达时间信息，并将其绘制到工作底图上，形成洪水到达时间分布图。

3.5.3 溃坝生命损失计算

溃坝生命损失计算方法很多，主要有Brown&Graham法（1988）^［69］、DeKay&McClelland法（1993）^［70］、Assaf法（1997）^［71］、Graham法（1999）^{［72，73］}、McClelland&Bowles法（2001）^［74］、RESCDAM法（2001）^［75］和李—周法（2006）^［76］。

目前在我国推荐使用Dekay&McClelland法或李—周法。

（1）Dekay&McClelland法。采用Dekay&McClelland法时，生命损失L_OL可按式（3.71）计算：

式中 L_OL——溃坝洪水淹没范围内的生命损失，人；

P_AR——溃坝洪水淹没范围内的风险人口，人；

W_T——警报时间，h；

F——溃坝洪水严重性S_d的函数符号，取值范围为0～1。对于高度严重溃坝洪水，取F=1；对于低度严重溃坝洪水，取F=0；对中度严重溃坝洪水，取F=0.5。

风险人口计算方法有静态统计法和动态统计法，静态统计法又有人口密度估算法、居民点或居住单元数目累计估算法，以及基于土地利用类型的风险人口估算法。

采用人口密度估算法时，风险人口P_AR可按式（3.72）计算：

式中 ρ——溃坝洪水淹没范围内人口密度，假设均匀分布，人/km²；

A——溃坝洪水淹没面积，km²。

采用居民点或居住单元数目累计估算法时，风险人口P_AR可分别按式（3.73）和式（3.74）计算：

式中 P_ARi——第i个居民点的人口数量，通过调查统计和人口登记数据确定，人；

——单个居民点人口均值，通过调查统计和人口登记数据确定，人/个；

n——居民点数目，个。

采用基于土地利用类型的风险人口估算法时，风险人口P_AR可按式（3.75）计算：

式中 A_i——土地利用类型i的面积，km²；

ρ_i——土地利用类型i的人口密度，人/km²；

n——土地利用类型总数，个。

采用动态统计法时，风险人口P_AR可按式（3.76）计算：

式中 P_ARi（t）——第i个居民点在时刻t的人口，人；

n——居民点数目，个；

P_ARCi——第i个居民点人口常量，通过调查统计和人口登记数据确定，人；

f_i（t）——第i个居民点在时刻t的人口分布系数。

溃坝洪水严重性S_d可按式（3.77）计算：

式中 h——溃坝洪水淹没范围内某点的水深，m；

v——相应某点的流速，m/s。

溃坝洪水严重性S_d标准划分如下：当S_d≤3.0m²/s时，为低度严重；当3.0m²/s<S_d≤7.0m²/s时，为中度严重；当S_d>7.0m²/s时，为高度严重。

（2）李—周法。采用李—周法时，溃坝生命损失L_OL可按式（3.78）计算：

式中 L_OL——溃坝洪水淹没范围内的生命损失，人；

P_AR——溃坝洪水淹没范围内的风险人口，人；

f——风险人口死亡率，可按表3.21确定。夏天、晴天、白天宜取此表的下限值，冬天、雨天、夜间宜取此表的上限值。

表3.21中，当风险人口接到溃坝警报后，对溃坝洪水可能淹没范围和严重程度缺乏足够了解，对逃生的必要性、措施、路径没有正确的理解和反应时，可认为风险人口对溃坝洪水严重性的理解程度U_d是模糊的；反之则认为风险人口对溃坝洪水严重性的理解程度U_d是明确的。

警报时间W_T是指大坝下游风险人口接受到撤退警报到溃坝洪水到达之间的逃脱时间，是影响和确定溃坝生命损失的一个极重要参数。一般认为当W_T<15min时，无警报；当15min≤W_T≤60min时，部分警报；当W_T>60min时，充分警报。警报时间长短受客观环境和人为因素的双重制约，若溃坝发生在白天、现场有管理人员或其他人员、有仪器直接监测（控）的水库，容易及时发布警报，警报时间长；若溃坝发生在夜间，则不易发现，及时发布溃坝警报的可能性小，警报时间短。离大坝越近的地区，洪水到达所经历的时间越短，警报时间也越短。

若水库拥有较高的管理水平、较强的预警能力和通畅的警报发布设施，能够提前向下游居民发布准确的溃坝警报，警报时间达到1h以上，那么有助于风险人口及时安全地撤离，生命损失将会大大地减少，甚至不会造成任何生命损失。

警报传递可以是电子警报器、广播喇叭、蜂鸣器、电视、电话、短信等方式，在偏远的地区还可以采用人工传递警报、扩音喇叭喊话、吹哨子、敲打锣（鼓）、发射信号弹等方式。警报传递方式影响着风险人口接受警报的效果，警报时间受其影响。

3.5.4 溃坝经济损失计算

溃坝经济损失包括直接经济损失和间接经济损失^{［77，78］}。

直接经济损失包括水库工程损毁所造成的经济损失和洪水直接淹没所造成的可用货币计量的各类损失。直接淹没损失包括工业、农业、林业、牧业、副业、渔业、商业、交通、邮电、文教卫生、粮油贮存、工程设施、物资库存、农业机械、房屋、群众家产和专项损失这17类。

间接经济损失是指直接经济损失以外的可用货币计量的损失，主要包括由于采取各种措施（如防汛、抢险、避难、开辟临时交通线等）而增加的费用、骨干交通线路中断给有关工矿企业造成原材料中断而停工停产及产品积压的损失或运输绕道增加的费用、农产品减产给农产品加工企业和轻工业造成的损失等。

3.5.4.1 溃坝直接经济损失计算

溃坝直接经济损失可采用分类损失率法、单位面积综合损失法和人均综合损失法等方法计算。

（1）分类损失率法。采用分类损失率法时，溃坝直接经济损失D可按式（3.79）计算：

式中 R_i——第i个行政区的各类财产损失总值，万元；

R_ij——第i个行政区内、第j类财产的损失值，万元；

V_ijk——第i个行政区内、第k级淹没水深下第j类资产价值，万元；

η_ijk——第i个行政区内、第k级淹没水深下第j类资产损失率，根据溃坝洪水严重性、历时等因素确定，%；

n——行政区数；

m——资产种类数；

l——淹没水深等级数。

（2）单位面积综合损失法和人均综合损失法。采用单位面积综合损失法和人均综合损失法时，溃坝直接经济损失D可按式（3.80）或式（3.81）计算：

式中 A——溃坝洪水淹没范围，km²；

L_A——溃坝洪水淹没范围内单位面积损失值，万元/km²；

P_AR——溃坝洪水淹没范围内的风险人口，人；

L_P——风险人口人均损失值，万元/人。

3.5.4.2 溃坝间接经济损失计算

溃坝间接经济损失可采用系数折算法和调查分析法计算。

（1）系数折算法。采用系数折算法时，溃坝间接经济损失S可按式（3.82）计算：

式中 R_i——第i个行政区的直接经济损失总值，万元；

k_i——系数，可根据实际洪灾损失调查资料确定，缺少资料时，可取k_i=0.63；

n——行政区数。

（2）调查分析法。调查分析法应通过实地调查溃坝洪水淹没区社会经济受灾程度，在相关的社会经济统计资料基础上，运用数理统计及时间序列分析等方法估算受灾区的间接经济损失。

3.5.5 溃坝社会与环境影响计算

社会与环境影响涉及面很广，又很复杂，王仁钟等对社会与环境影响的各主要因素进行了量化，再综合成社会与环境影响指数^［34］。

社会影响要素：社会影响除生命损失外，主要包括政治影响（即对国家、社会安定的不利影响），因受伤或精神压力给人们造成身心健康的伤害，以及日常生活水平和生活质量的下降等，无法补救的文物古迹、艺术珍品和稀有动植物等的损失。

环境影响要素：环境影响主要包括河道形态的影响，生物及其生长栖息地（包括河流、湿地、表土和植被等）的丧失，人文景观（含公园与保护区）的破坏，以及易受影响或造成重大环境影响或污染的工业（包括核设施、化学储存设施、农药厂等）影响等。

本书中，溃坝社会与环境影响主要考虑溃坝洪水淹没范围内风险人口数量、城镇规模、基础设施重要性、文物古迹级别、河道形态破坏程度、动植物栖息地保护级别、自然景观级别、潜在污染企业规模等因素。溃坝社会与环境影响采用社会与环境影响指数计算。溃坝社会与环境影响指数定义为所考虑的社会与环境要素的乘积，按式（3.83）计算：

式中 I_SE——社会与环境影响指数；

C ₁——风险人口系数；

C ₂——城镇规模系数；

C ₃——基础设施重要性系数；

C ₄——文物古迹级别系数；

C ₅——河道形态破坏程度系数；

C ₆——动植物栖息地保护级别系数；

C ₇——自然景观级别系数；

C ₈——潜在污染企业规模系数。

社会与环境影响因素赋值在王仁钟^［34］的基础上做了改进，可查表3.22。