CN104102815B - 一种滨海核电厂厂区积水深度计算方法 - Google Patents
一种滨海核电厂厂区积水深度计算方法 Download PDFInfo
- Publication number
- CN104102815B CN104102815B CN201410275826.9A CN201410275826A CN104102815B CN 104102815 B CN104102815 B CN 104102815B CN 201410275826 A CN201410275826 A CN 201410275826A CN 104102815 B CN104102815 B CN 104102815B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mfrac
- msub
- msup
- rainwater
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 42
- 238000000205 computational method Methods 0.000 title abstract 2
- 239000013535 sea water Substances 0.000 claims abstract description 16
- 239000002352 surface water Substances 0.000 claims abstract description 16
- 230000008878 coupling Effects 0.000 claims abstract description 4
- 238000010168 coupling process Methods 0.000 claims abstract description 4
- 238000005859 coupling reaction Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000003746 surface roughness Effects 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 238000005086 pumping Methods 0.000 claims description 3
- 238000007711 solidification Methods 0.000 claims description 3
- 230000008023 solidification Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000012279 drainage procedure Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 abstract 1
- 230000002093 peripheral effect Effects 0.000 abstract 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种滨海核电厂厂区积水深度计算方法，包括：(Ⅰ)根据核电厂厂区降雨情况、外海越浪、地面坡度和厂区地表固化情况,建立可表征地表水深与雨水管网下泄量的关系的地表径流模型；(Ⅱ)针对雨水管网排水过程，考虑雨水受外海海水的顶托作用，建立雨水管网水动力模型；(Ⅲ)对厂区进行计算网格划分，通过使雨水管网下泄量和雨水在管网系统中的排出量相等，实现地表径流模型和管网水动力学模型的耦合计算，确定厂区不同区域的积水深度。本发明通过对厂区降雨和厂区外围海水越浪造成的影响进行组合计算，所得到的积水深度的数值更加科学合理，据此可有针对性地对核电厂重要构筑物提出防洪预案和防洪措施，以提高厂区的防洪能力。
Description
技术领域
本发明属于积水深度计算技术领域，具体涉及一种滨海核电厂厂区积水深度计算方法。
背景技术
滨海核电厂根据厂区内构筑物的重要性不同，针对不同区域的雨水管网系统采用不同设计标准。核岛所在的主厂区采用千年一遇暴雨强度设计，可能最大暴雨强度进行校核。厂前区雨水管网系统按照一般大型工矿企业的设计标准设计。厂区构筑物屋顶设有屋面排水系统，来收集、输送和排泄雨水，厂区地面通过设置不同坡度引导雨水地面径流，最终厂区雨水通过地面设置的雨水口进入地下雨水管网系统，并最终排入大海。
由于恶劣的天气系统，如风暴或台风，除带来大量的降雨外，同时还可能会引起外海潮位的增长，外海在台风的影响下，波浪会巨大，台风会导致部分海水越过厂区的防护堤进入厂区，这部分越浪海水会在厂区地面产生漫流。由于厂区的雨水系统设计仅考虑降雨产生的地面雨水的排出，而未考虑到越浪海水的排出，同时也未考虑此时外海潮位很大，雨水管网出口位于海面以下，雨水的排水受到外海潮位的顶托，造成厂区雨水的排水不畅，从而可能导致厂区雨水以及越浪海水的滞留，引起厂区局部水淹。
发明内容
本发明的目的是提供一种同时考虑降雨(含雨水管网排放)以及越浪海水漫流情况的滨海核电厂厂区积水深度计算方法，为厂区防水淹提供基础数据。
为达到以上目的，本发明采用的技术方案是：一种滨海核电厂厂区积水深度计算方法，包括以下步骤：
(Ⅰ)根据核电厂厂区降雨情况、外海越浪、地面坡度和厂区地表固化情况,建立地表径流模型，所述地表径流模型用于表征地表水深与雨水管网下泄量的关系；
(Ⅱ)针对雨水管网排水过程，考虑雨水受外海海水的顶托作用，建立雨水管网水动力模型，所述雨水管网水动力模型用于表征地表水深与雨水在管网系统中的排出量的关系；
(Ⅲ)对厂区进行计算网格划分，通过使雨水管网下泄量和雨水在管网系统中的排出量相等，实现步骤(Ⅰ)的地表径流模型和步骤(Ⅱ)的管网水动力学模型在每个时间步长内的耦合计算，确定出厂区不同区域的积水深度，即地表水深。
进一步，步骤(Ⅰ)中，所述地表径流模型包括地表水流连续性方程(1)和地表水流动力学方程，地表水流动力学方程包括X方向的动量方程(2)和Y方向的动量方程(3)：
其中：
h表示地表水深；p表示降雨量；qg表示雨水井下泄量；u、v分别表示空间变量x、y方向的垂向平均速度；t表示时间变量；z表示越浪海水在厂区地面引起水面高程；γt表示紊动粘性系数；n表示地表糙率；
步骤(Ⅱ)中，所述雨水管网水动力模型包括雨水管网中水流连续性方程(4)和雨水管网中水流动量方程(5)以及管道中水流连续性方程方程(6)和管道中水流动量方程(7)：
其中：
Qi表示第i个连接管道出入该节点的流量；
Qj表示直接进入节点的变化的入流量或节点的抽水量或溢流量或渗漏量；
Anj表示节点的水平断面面积；
t表示时间变量；
Yj表示节点水深；
Zi、Pi、vi分别表示第i个管道出口的底高、大气压和速度；
Hj表示每单位容积入流的净能量输入，以水头的形式表示；
Ki表示第i个管道入口或出口的损失系数；
γ表示液体的容重；g表示重力加速度；
Q表示雨水在管网系统中的排出量；
x表示水流沿着管道长度方向移动的距离；
q表示单位长度管道的旁侧入流量；
A表示过水断面的面积；
h1表示管道断面平均水深；
θ为管道与水平面的夹角；so为管道坡度，且so＝sinθ；
β表示动量变动改正系数；K和K'为非水力静态分布的改正系数；
T表示作用于过水断面的内部压力产生的作用力；
Ux表示x方向的旁侧入流流速；
sf表示摩阻比降。
再进一步，式(7)中，sf采用下述式(8)计算：
其中，g表示重力加速度，n表示地表糙率，Q表示雨水在管网系统中的排出量，h1表示管道断面平均水深。
本发明通过对厂区降雨和厂区外围海水越浪造成的影响进行组合计算，得到同时考虑了外海洪水影响的不同降雨强度下厂区的不同区域的积水深度，所得到的积水深度的数值更加科学合理，能够有针对性地对核电厂重要构筑物提出防洪预案和防洪措施，可大大提高厂区的防洪能力。
具体实施方式
下面结合具体实施方式对本发明作进一步描述。
本发明所提供的一种滨海核电厂厂区积水深度计算方法，包括以下步骤：
第一步，根据核电厂厂区降雨、外海越浪、地面坡度和厂区地表固化情况(是指与地表径流系数有关的因素，如硬化比例，绿化情况，反映了降雨在地表汇流和下渗的特征),建立用于表征地表水深与雨水管网下泄量的关系的地表径流模型。
本发明中，所建立的地表径流模型包括地表水流连续性方程(1)和地表水流动力学方程，地表水流动力学方程包括X方向的动量方程(2)和Y方向的动量方程(3)：
其中：p表示降雨量；qg表示雨水井下泄量(即雨水管网下泄量)；u、v分别表示空间变量x、y方向的垂向平均速度；t表示时间变量；z表示越浪海水在厂区地面引起水面高程；h表示地表水深；γt表示紊动粘性系数；n表示地表糙率(可直接根据室外给排水规范选取)。
可根据计算区域情况确定厂区地表边界类型、边界越浪海水和降雨强度p，通过建筑物的产汇流计算及海水在地面漫流情况确定雨水在地面引起地表水深h、越浪海水在地面引起水面高程z h和雨水井下泄量qg。
上述偏微分方程组是一个抛物型偏微分方程，可以用六点差分格式求解。
第二步，针对雨水管网排水过程，考虑雨水受外海海水的顶托作用，建立用于表征地表水深与雨水在管网系统中的排出量的关系的雨水管网水动力模型。
本发明中，所建立的雨水管网水动力模型包括雨水管网中水流连续性方程(4)和雨水管网中水流动量方程(5)以及管道中水流连续性方程方程(6)和管道中水流动量方程(7)：
其中：
Qi表示第i个连接管道出入该节点的流量；Qj表示直接进入节点的变化的入流量或节点的抽水量或溢流量或渗漏量；Anj表示节点的水平断面面积；t表示时间变量；Yj表示节点水深；Zi、Pi、vi分别表示第i个连接管道出口的底高、大气压和速度；Hj表示每单位容积入流的净能量输入，以水头的形式表示；Ki表示第i个连接管道入口或出口的损失系数；γ表示液体的容重；g表示重力加速度；Q表示雨水在管网系统中的排出量；x表示水流沿着管道长度方向移动的距离；q表示单位长度管道的旁侧入流量；A表示过水断面的面积；h1表示管道断面平均水深；θ为管道与水平面的夹角；so为管道坡度，且so＝sinθ；β表示动量变动改正系数；K和K'为非水力静态分布的改正系数；T表示作用于过水断面的内部压力产生的作用力；Ux表示x方向的旁侧入流流速；sf表示摩阻比降。
其中，式(7)中的sf采用下述式(8)计算：
其中，g表示重力加速度，n表示地表糙率，Q表示雨水在管网系统中的排出量，h1表示管道断面平均水深。
在这一步中，可根据雨水管网设计图纸资料，确定雨水管网计算节点的深度Yj、节点的水平断面面积Anj、连接管道出口的底高Zi、大气压Pi和速度vi，雨水在管网系统中的排出量Q等，计算雨水在管网中的流动情况。
第三步，对厂区进行计算网格划分，通过使雨水管网下泄量qg和雨水在管网系统中的排出量Q相等，实现第一步的地表径流模型和第二步的管网水动力学模型在每个时间步长内的耦合计算，确定出厂区不同区域的积水深度数值。
根据滨海核电厂海域地理环境，厂区可分为多个区段，由厂区各个区段所处的位置分别确定防洪护堤设计标准。厂区护堤护面结构在不同的堤段采用不同的结构形式，以节省工程投资。
厂区内不同区段的积水深度跟各个区段所处的位置有关，因此计算积水深度时要对厂区进行计算网格划分，对不同的区段，使用不同的初始条件和边界条件求解。
上述实施例只是对本发明的举例说明，本发明也可以以其它的特定方式或其它的特定形式实施，而不偏离本发明的要旨或本质特征。因此，描述的实施方式从任何方面来看均应视为说明性而非限定性的。本发明的范围应由附加的权利要求说明，任何与权利要求的意图和范围等效的变化也应包含在本发明的范围内。
Claims (2)
1.一种滨海核电厂厂区积水深度计算方法，包括以下步骤：
(Ⅰ)根据核电厂厂区降雨情况、外海越浪、地面坡度和厂区地表固化情况,建立地表径流模型，所述地表径流模型用于表征地表水深与雨水管网下泄量的关系；
所述地表径流模型包括地表水流连续性方程(1)和地表水流动力学方程，地表水流动力学方程包括X方向的动量方程(2)和Y方向的动量方程(3)：
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mi>v</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>p</mi>
<mo>-</mo>
<msub>
<mi>q</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mi>u</mi>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mi>v</mi>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>g</mi>
<mi>h</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mi>gn</mi>
<mn>2</mn>
</msup>
<mi>h</mi>
<msqrt>
<mrow>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mi>u</mi>
</mrow>
<msup>
<mi>h</mi>
<mrow>
<mn>4</mn>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msup>
</mfrac>
<mo>=</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>h&gamma;</mi>
<mi>t</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>h&gamma;</mi>
<mi>t</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mi>v</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mi>u</mi>
<mi>v</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mi>v</mi>
<mi>v</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>g</mi>
<mi>h</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mi>gn</mi>
<mn>2</mn>
</msup>
<mi>h</mi>
<msqrt>
<mrow>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mi>v</mi>
</mrow>
<msup>
<mi>h</mi>
<mrow>
<mn>4</mn>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msup>
</mfrac>
<mo>=</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>h&gamma;</mi>
<mi>t</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>v</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>h&gamma;</mi>
<mi>t</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>v</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，
h表示地表水深；p表示降雨量；qg表示雨水井下泄量；u、v分别表示空间变量x、y方向的垂向平均速度；t表示时间变量；z表示越浪海水在厂区地面引起水面高程；γt表示紊动粘性系数；n表示地表糙率；
(Ⅱ)针对雨水在管网的排出过程，同时考虑雨水受外海海水的顶托作用，建立雨水管网水动力模型，所述雨水管网水动力模型用于表征地表水深与雨水在管网系统中的排出量的关系；
所述雨水管网水动力模型包括雨水管网中水流连续性方程(4)和雨水管网中水流动量方程(5)以及管道中水流连续性方程方程(6)和管道中水流动量方程(7)：
<mrow>
<msub>
<mi>&Sigma;Q</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mi>j</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msub>
<mi>dY</mi>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mi>j</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msub>
<mi>dY</mi>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>&Sigma;</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<mi>g</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mi>&gamma;</mi>
</mfrac>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mi>j</mi>
</msub>
<msub>
<mi>H</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mo>&Sigma;</mo>
<msub>
<mi>Q</mi>
<mi>j</mi>
</msub>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mfrac>
<mrow>
<msup>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<mi>g</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>A</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>Q</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>q</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mi>g</mi>
<mi>A</mi>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>Q</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>g</mi>
<mi>A</mi>
</mrow>
</mfrac>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&beta;</mi>
<mi>A</mi>
</mfrac>
<msup>
<mi>Q</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>cos</mi>
<mi>&theta;</mi>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>Kh</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>-</mo>
<msup>
<mi>K</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mfrac>
<mn>1</mn>
<mi>A</mi>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>A</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>s</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>&gamma;</mi>
<mi>A</mi>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>qU</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<mi>g</mi>
<mi>A</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中：
Qi表示第i个连接管道出入节点的流量；
Qj表示直接进入节点的变化的入流量或节点的抽水量或溢流量或渗漏量；
Anj表示节点的水平断面面积；
Yj表示节点水深；
Zi、Pi、vi分别表示第i个管道出口的底高、大气压和速度；
Hj表示每单位容积入流的净能量输入，以水头的形式表示；
Ki表示第i个管道入口或出口的损失系数；
γ表示液体的容重；g表示重力加速度；
Q表示雨水在管网系统中的排出量；
x表示水流沿着管道长度方向移动的距离；
q表示单位长度管道的旁侧入流量；
A表示过水断面的面积；
h1表示管道断面平均水深；
θ为管道与水平面的夹角；so为管道坡度，且so＝sinθ；
β表示动量变动改正系数；K和K'为非水力静态分布的改正系数；
T表示作用于过水断面的内部压力产生的作用力；
Ux表示x方向的旁侧入流流速；
sf表示摩阻比降；
(Ⅲ)对厂区进行计算网格划分，通过使雨水管网下泄量和雨水在管网系统中的排出量相等，实现步骤(Ⅰ)的地表径流模型和步骤(Ⅱ)的管网水动力学模型在每个时间步长内的耦合计算，确定出厂区不同区域的积水深度，即地表水深。
2.根据权利要求1所述的一种滨海核电厂厂区积水深度计算方法，其特征在于，式(7)中，sf采用下述式(8)计算：
<mrow>
<msub>
<mi>s</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>gn</mi>
<mn>2</mn>
</msup>
<msup>
<mi>Q</mi>
<mn>2</mn>
</msup>
</mrow>
<msubsup>
<mi>h</mi>
<mn>1</mn>
<mrow>
<mn>4</mn>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msubsup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，g表示重力加速度，n表示地表糙率，Q表示雨水在管网系统中的排出量，h1表示管道断面平均水深。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410275826.9A CN104102815B (zh) | 2014-06-19 | 2014-06-19 | 一种滨海核电厂厂区积水深度计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410275826.9A CN104102815B (zh) | 2014-06-19 | 2014-06-19 | 一种滨海核电厂厂区积水深度计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104102815A CN104102815A (zh) | 2014-10-15 |
| CN104102815B true CN104102815B (zh) | 2017-10-24 |
Family
ID=51670960
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410275826.9A Active CN104102815B (zh) | 2014-06-19 | 2014-06-19 | 一种滨海核电厂厂区积水深度计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104102815B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108151967B (zh) * | 2017-11-16 | 2021-05-18 | 中国核电工程有限公司 | 一种核电厂核岛厂房液体泄漏率测量装置 |
| CN111369102B (zh) * | 2020-02-04 | 2021-01-19 | 中国水利水电科学研究院 | 内涝风险点的提取方法及装置 |
| CN111985082B (zh) * | 2020-07-23 | 2024-04-26 | 西安理工大学 | 一种高守恒性水文水动力耦合模拟方法 |
| CN118296715A (zh) * | 2024-06-06 | 2024-07-05 | 中交天津港湾工程研究院有限公司 | 核电厂区积水仿真计算溯源方法、装置、终端及存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101692273A (zh) * | 2009-09-28 | 2010-04-07 | 北京工业大学 | 城市排水管网在线水力模型建模方法 |
| CN101964020A (zh) * | 2010-09-21 | 2011-02-02 | 西安交通大学 | 一种江河流域水动力模型取水量的预测方法 |
| CN103473435A (zh) * | 2012-06-06 | 2013-12-25 | 邵丹娜 | 城市雨水排水调蓄计算模型及方法 |
-
2014
- 2014-06-19 CN CN201410275826.9A patent/CN104102815B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101692273A (zh) * | 2009-09-28 | 2010-04-07 | 北京工业大学 | 城市排水管网在线水力模型建模方法 |
| CN101964020A (zh) * | 2010-09-21 | 2011-02-02 | 西安交通大学 | 一种江河流域水动力模型取水量的预测方法 |
| CN103473435A (zh) * | 2012-06-06 | 2013-12-25 | 邵丹娜 | 城市雨水排水调蓄计算模型及方法 |
Non-Patent Citations (2)
| Title |
|---|
| 在役核电厂防水淹改造及探讨;陈坚刚 等;《中国核科学技术进展报告（第三卷）——中国核学会2013年学术年会论文集第2册（铀矿冶分卷、核能动力分卷（上））》;20130911;第198-202页 * |
| 极端降雨下核电厂雨水排水系统的设计和校核;潘俊杰 等;《中国给水排水》;20130531;第29卷(第10期);第31-34页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104102815A (zh) | 2014-10-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104732073B (zh) | 地表水‑地下水耦合模拟的计算方法 | |
| CN102930357B (zh) | 岩溶隧道暗河暴雨涌水洪峰峰值及峰值时间的预测方法 | |
| CN113723024B (zh) | 一种适用于滨海地区的“溪流”-“河道”-“河口”分布式洪水过程模拟方法 | |
| CN109919372A (zh) | 一种基于全时空的城市暴雨积水评估建模方法 | |
| CN104102815B (zh) | 一种滨海核电厂厂区积水深度计算方法 | |
| CN103147464A (zh) | 一种地下室抗浮板的泄压结构及其施工方法 | |
| CN113792448B (zh) | 一种河道与泛区凌汛壅水-溃堤-淹没耦合模拟方法 | |
| CN107916646A (zh) | 一种雨洪资源挖潜量估算方法 | |
| CN115828786A (zh) | 一种次洪过程中淤地坝系减沙量的计算方法 | |
| CN108625337B (zh) | 一种确定潮流界以下沙质河床河段整治水位的方法 | |
| CN109736259B (zh) | 山区河流复式滩地利用范围及防洪安全位置的水力学确定方法 | |
| CN107816031A (zh) | 利用高含沙洪水的沉积物固定流动沙丘的方法 | |
| CN113962475A (zh) | 一种下穿岩溶天坑隧道暴雨期最大涌水量峰值估计方法 | |
| CN104195979A (zh) | 河道交汇水流停滞区楔锥体及其构建方法与应用 | |
| CN102720266A (zh) | 虹吸式蓄排水隧道清淤装置 | |
| CN203174745U (zh) | 雨水竖向分区排水系统 | |
| CN102778260B (zh) | 一种潜水完整井或井群动态降水过程中排水量及排水时间的计算方法 | |
| CN102768050B (zh) | 一种小流域全自动三级测流装置及其使用方法 | |
| CN103981838B (zh) | 一种护坡营造方法 | |
| Huang et al. | Numerical study on the impact of Gongji Road Rain Pump on the waterlogging in Huinan, Pudong District | |
| CN113779814A (zh) | 一种大尺度台风洪涝模拟计算方法 | |
| CN202141477U (zh) | 一种小流域全自动三级测流装置 | |
| Dong et al. | Sustainable development of water resources and hydraulic engineering in China | |
| CN112884279B (zh) | 一种沟道梯级土地整治工程暴雨洪水安全风险评估方法 | |
| Yin et al. | Effectiveness of restoring a saltmarsh in the eroded steep coast |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |