CN113312771B - Calculation method and application of limited wind speed of side key smoke exhaust of tunnel - Google Patents

Calculation method and application of limited wind speed of side key smoke exhaust of tunnel Download PDF

Info

Publication number
CN113312771B
CN113312771B CN202110598123.XA CN202110598123A CN113312771B CN 113312771 B CN113312771 B CN 113312771B CN 202110598123 A CN202110598123 A CN 202110598123A CN 113312771 B CN113312771 B CN 113312771B
Authority
CN
China
Prior art keywords
wind speed
tunnel
fire
dimensionless
smoke exhaust
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
Application number
CN202110598123.XA
Other languages
Chinese (zh)
Other versions
CN113312771A (en
Inventor
姜学鹏
刘婵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University of Science and Technology WHUST
Original Assignee
Wuhan University of Science and Technology WHUST
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University of Science and Technology WHUST filed Critical Wuhan University of Science and Technology WHUST
Priority to CN202110598123.XA priority Critical patent/CN113312771B/en
Publication of CN113312771A publication Critical patent/CN113312771A/en
Application granted granted Critical
Publication of CN113312771B publication Critical patent/CN113312771B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

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)
  • Ventilation (AREA)

Abstract

The invention relates to a method for calculating limited wind speed of side key smoke exhaust of a tunnel and application thereof, comprising the following specific steps of: determining a limited wind speed ucDefining, determining influencing factors, establishing a limited wind speed ucA relationship to an influencing factor; determining basic dimension, and establishing limited wind speed ucA dimensional relation with the influencing factor; determining basic physical quantity, obtaining dimensionless items of influence factors according to the pi theorem, and deducing to obtain the limited wind speed ucThe dimensionless calculation formula of (1); obtaining the numerical value of the smoke backflow length L of the fire under different working conditions through numerical simulation, obtaining the limited wind speed through the smoke backflow length, drawing the simulation result into a scatter diagram, and determining the influence factors on the limited wind speed ucThe influence of (c); performing data fitting on the result to obtain the value of each unknown coefficient in the dimensionless relational expression, and further establishing the dimensionless limited wind speed ucThe calculation formula of (2). The method has the advantages of being simple, suitable for different tunnels, greatly influencing the efficient utilization of the smoke exhaust capacity of the tunnels and providing guidance for fire smoke control and fire rescue of the side key smoke exhaust tunnels.

Description

一种隧道侧部重点排烟受限风速的计算方法及应用A calculation method and application of limited wind speed for key smoke exhaust on the side of tunnel

技术领域technical field

本发明涉及隧道火灾烟气控制技术领域,具体涉及一种隧道侧部重点排烟受限风速的计算方法及应用。The invention relates to the technical field of tunnel fire smoke control, in particular to a calculation method and application of a limited wind speed for key smoke exhaust at the side of a tunnel.

背景技术Background technique

火灾是隧道运营的最大风险之一。由于隧道具有狭长且封闭的特性,导致隧道相比普通地面建筑而言,具有独特空间结构和火灾特性。隧道一旦发生火灾,烟气无法及时排出,极易造成人员窒息或中毒,隧道内车辆来不及驶离就会形成二次伤害,造成群死群伤的重大灾害性事故。控制受限风速有利于高效利用隧道的排烟能力,将烟气控制在一定范围内从而保证人员疏散与救援。Fire is one of the greatest risks to tunnel operations. Compared with ordinary ground buildings, tunnels have unique spatial structure and fire characteristics due to the narrow and closed characteristics of tunnels. In the event of a fire in the tunnel, the smoke cannot be discharged in time, which will easily cause suffocation or poisoning of people. If the vehicles in the tunnel are too late to leave, secondary injuries will occur, resulting in major catastrophic accidents involving mass deaths and injuries. Controlling the limited wind speed is conducive to the efficient use of the smoke exhaust capacity of the tunnel, and the smoke is controlled within a certain range to ensure personnel evacuation and rescue.

隧道火灾场景下受限风速uc的定义:烟气逆流长度对变动的排烟口流速不敏感时对应的最小纵向诱导风速称为受限风速uc,烟气逆流长度指火源在隧道行进方向上与排烟口的距离,如图2,受限风速uc是保障烟气逆流长度不因排烟口流速降低而增加的纵向诱导风的合适风速,是控制烟气不扩散蔓延的控制指标。目前的研究主要探索了临界风速的计算方法,临界风速是刚好使烟气逆流长度为零的纵向风速,可确保火源上游没有烟气,但采用临界风速控制烟气蔓延,一方面由于临界风速理论上大于受限风速,故不能使烟气向最近的上游排烟口排烟;另一方面,侧部点式排烟模式火灾临界风速的预测模型与侧部无量纲排烟量、无量纲火源热释放速率、无量纲排烟口距火源的距离相关,因此用临界风速模型预测的纵向风速控制烟气的逆流,仍受排烟量的变化的影响。The definition of the restricted wind speed uc in the tunnel fire scenario: the minimum longitudinal induced wind speed corresponding to the flue gas reverse flow length is not sensitive to the fluctuating exhaust flow velocity is called the restricted wind speed uc , and the flue gas reverse flow length refers to the fire source traveling in the tunnel. The distance from the exhaust port in the direction, as shown in Figure 2, the limited wind speed u c is the appropriate wind speed to ensure that the length of the reverse flow of the flue gas does not increase due to the decrease of the flow rate of the exhaust port. index. The current research mainly explores the calculation method of the critical wind speed. The critical wind speed is the longitudinal wind speed that just makes the counterflow length of the flue gas zero, which can ensure that there is no flue gas upstream of the fire source, but the critical wind speed is used to control the spread of the flue gas. On the one hand, due to the critical wind speed Theoretically, it is greater than the limited wind speed, so the smoke cannot be exhausted to the nearest upstream smoke outlet; The heat release rate of the fire source and the distance between the dimensionless smoke outlet and the fire source are related. Therefore, the longitudinal wind speed predicted by the critical wind speed model to control the reverse flow of the flue gas is still affected by the change of the smoke exhaust volume.

目前,已有的受限风速计算方法多是针对排烟口位于隧道顶部的情况,根据烟气逆流长度与隧道高度的比值来确定受限风速,应用侧部重点排烟模式的隧道中排烟口位于隧道侧壁,相比于顶部排烟,此时烟气的受力状态与运动状态已发生改变,国内外尚无针对侧部重点排烟模式隧道火灾受限风速的计算方法。At present, most of the existing limited wind speed calculation methods are aimed at the situation where the smoke exhaust port is located at the top of the tunnel. The mouth is located on the side wall of the tunnel. Compared with the smoke exhaust at the top, the stress state and motion state of the smoke have changed at this time. There is no calculation method for the limited wind speed of the tunnel fire in the key smoke exhaust mode at home and abroad.

发明内容SUMMARY OF THE INVENTION

本发明所要解决的技术问题是提供一种隧道侧部重点排烟火灾受限风速的计算方法及应用,旨在解决上述技术问题。The technical problem to be solved by the present invention is to provide a calculation method and application of the limited wind speed of the key smoke exhaust fire at the side of the tunnel, aiming at solving the above technical problems.

本发明解决上述技术问题的技术方案如下:The technical scheme that the present invention solves the above-mentioned technical problems is as follows:

一种隧道侧部重点排烟火灾受限风速的计算方法,包括以下具体步骤:A method for calculating the limited wind speed of a key smoke exhaust fire at the side of a tunnel, comprising the following specific steps:

S1:确定隧道火灾场景下受限风速uc的影响因素,建立所述受限风速uc与所述影响因素的关系式:f(uc,Q,ρ0,CP,T0,g,HD)=0;S1: Determine the influencing factors of the limited wind speed uc in the tunnel fire scenario, and establish the relationship between the limited wind speed uc and the influencing factors: f ( uc ,Q,ρ 0 , C P ,T 0 ,g , HD )=0;

所述影响因素包括火源热释放速率Q、空气密度ρ0、空气定压比热Cp、空气温度T0、重力加速度g、隧道高度HD,物理参数的数量n为7,火源热释放速率Q的单位为kg.m2/s3,空气密度ρ0的单位为kg/m3,空气定压比热Cp的单位为m2/s2.K,空气温度T0的单位为K,重力加速度g单位为m/s2,隧道高度HD的单位为m;The influencing factors include fire source heat release rate Q, air density ρ 0 , air constant pressure specific heat C p , air temperature T 0 , gravitational acceleration g, tunnel height H D , the number of physical parameters n is 7, and the fire source heat The unit of release rate Q is kg.m 2 /s 3 , the unit of air density ρ 0 is kg/m 3 , the unit of air constant pressure specific heat C p is m 2 /s 2 .K, the unit of air temperature T 0 is K, the unit of gravitational acceleration g is m/ s 2 , and the unit of tunnel height HD is m;

S2:根据所述影响因素的单位确定基本量纲,将所述影响因素由所述基本量纲表示,建立所述受限风速uc与所述影响因素的量纲关系式;S2: Determine a basic dimension according to the unit of the influencing factor, represent the influencing factor by the basic dimension, and establish a dimensional relationship between the limited wind speed uc and the influencing factor;

S3:根据所述影响因素的单位确定基本量纲,将所述影响因素由所述基本量纲表示,建立所述受限风速uc与所述影响因素的量纲关系式;S3: Determine a basic dimension according to the unit of the influencing factor, represent the influencing factor by the basic dimension, and establish a dimensional relationship between the limited wind speed uc and the influencing factor;

所述步骤S3中,所述基本量纲包括质量M、时间t、长度L、温度T,基本量纲数η为4;In the step S3, the basic dimension includes mass M, time t, length L, and temperature T, and the basic dimension number η is 4;

所述步骤S3中,所述受限风速uc与所述影响因素的关系式为:In the step S3, the relationship between the limited wind speed uc and the influencing factors is:

f(Lt-1,ML2t-3,ML-3,L2t-2T-1,T,Lt-2,L)=0;f(Lt -1 , ML 2 t -3 , ML -3 , L 2 t -2 T -1 , T, Lt -2 , L)=0;

S4:确定所述影响因素的基本物理量,根据π定理得到所述影响因素的无量纲参数的数量为n-η=3,再根据π定理中循环量的选取原则,选取与长度L直接相关的变量HD,选取与温度T直接相关的变量T0,选取与时间t直接相关的变量g,选取与质量M直接相关的变量ρ0,作为循环量,用这4个循环量与其它n-η个物理参数中的其它所有参数依次组合成无量纲参数Π1、Π2、Π3,将所述步骤S3中的所述关系式转变为无量纲关系式,再获得所述受限风速uc的无量纲计算公式;S4: Determine the basic physical quantity of the influencing factor, obtain the number of dimensionless parameters of the influencing factor according to the π theorem as n-η=3, and then select the one directly related to the length L according to the selection principle of the circulation quantity in the π theorem Variable H D , select the variable T 0 directly related to the temperature T, select the variable g directly related to the time t, select the variable ρ 0 directly related to the mass M, as the circulation amount, use these 4 circulation amounts and other n- All other parameters in the n physical parameters are sequentially combined into dimensionless parameters Π 1 , Π 2 , Π 3 , the relational expression in the step S3 is converted into a dimensionless relational expression, and then the restricted wind speed u is obtained. The dimensionless calculation formula of c ;

所述步骤S4中,所述影响因素的无量纲参数Π1、Π2、Π3为:In the step S4, the dimensionless parameters Π 1 , Π 2 , and Π 3 of the influencing factors are:

Figure GDA0003574147540000031
Figure GDA0003574147540000031

所述受限风速uc的无量纲计算公式为:The dimensionless calculation formula of the restricted wind speed u c is:

Figure GDA0003574147540000032
Figure GDA0003574147540000032

Figure GDA0003574147540000033
uc *为无量纲风速,Q**为无量纲功率;which is
Figure GDA0003574147540000033
u c * is the dimensionless wind speed, Q ** is the dimensionless power;

S5:FDS数值模拟部分:建立隧道的缩尺几何模型,设置虚拟火源为丙烷燃烧,模拟不同的火源热释放速率,设置排烟口流速v从小到大逐步递增,得到FDS软件拟合出的不同烟气逆流长度L,当烟气逆流长度L随排烟口流速v的递增而减少的趋势变得缓慢时,此时的排烟口流速v在隧道内产生的诱导风速为该火灾条件下的实际受限风速uc的数值,并测量记录受限风速uc的影响因素的数值,隧道内的一侧壁上沿纵长方向设有两个排烟口,两个排烟口之间布置火源,火源位于单向三车道中间车道的中心线上,采用稳态火或t平方火,所述缩尺几何模型内沿顶壁的纵长方向上均匀间隔设有多个温度测点,提取FDS软件拟合出的受限风速uc的影响因素,所述缩尺几何模型的一处断面设置有多个流速测点,提取FDS软件拟合出的受限风速ucS5: FDS numerical simulation part: establish a scaled geometric model of the tunnel, set the virtual fire source as propane combustion, simulate different heat release rates of the fire source, set the flow velocity v of the exhaust port to gradually increase from small to large, and get the FDS software fitting When the flue gas counterflow length L decreases slowly with the increase of the flue gas flow velocity v, the induced wind speed generated in the tunnel by the flue gas flow velocity v at this time is the fire condition The value of the actual restricted wind speed uc at the lower end of the tunnel, and the values of the influencing factors of the restricted wind speed uc are measured and recorded. One side wall in the tunnel is provided with two smoke exhaust ports along the longitudinal direction, the The fire source is located on the center line of the one-way three-lane middle lane, and the steady-state fire or t-square fire is used. In the scaled geometric model, there are multiple temperature uniform intervals along the longitudinal direction of the top wall. Measuring points, extract the influencing factors of the restricted wind speed uc fitted by the FDS software, a section of the scaled geometric model is provided with a plurality of flow velocity measuring points, and extract the restricted wind speed uc fitted by the FDS software;

S6:实体模拟:建立隧道的缩尺几何模型,在缩尺几何模型中设置不同的火灾条件,包括用燃烧器作为火源,模拟不同的火源热释放速率,设置不同的排烟口流速v,并用来控制不同烟气逆流长度L,并在该火灾条件下改变纵向风速的大小,并观察烟气逆流长度不受排烟口流速影响的时刻,测量该时刻的实际受限风速uc的数值,并测量记录受限风速uc的影响因素的数值,隧道内的一侧壁上沿纵长方向设有两个排烟口,两个排烟口之间布置火源,火源位于单向三车道中间车道的中心线上,所述缩尺几何模型内沿顶壁的纵长方向上均匀间隔安装多个温度传感器,测量受限风速uc的影响因素,所述缩尺几何模型的一处断面安装风速仪,测量受限风速ucS6: Entity simulation: establish a scaled geometric model of the tunnel, and set different fire conditions in the scaled geometric model, including using the burner as the fire source, simulating different heat release rates of the fire source, and setting different flow velocity v of the exhaust port , and used to control different flue gas counterflow lengths L, and change the vertical wind speed under the fire conditions, and observe the moment when the flue gas counterflow length is not affected by the flow velocity of the exhaust port, and measure the actual restricted wind speed u c at this moment. The numerical value, and the numerical value of the influencing factors of the restricted wind speed uc are measured and recorded. There are two smoke exhaust ports on one side wall in the tunnel along the longitudinal direction, and the fire source is arranged between the two smoke exhaust ports. To the center line of the middle lane of the three lanes, a plurality of temperature sensors are installed at uniform intervals along the longitudinal direction of the top wall in the scaled geometric model to measure the influencing factors of the restricted wind speed uc . An anemometer is installed on one section to measure the restricted wind speed u c ;

S7:分别根据步骤S5和S6得到的不同火灾条件下受限风速uc和影响因素的数值,将模拟结果分别绘制成散点图,若拟合线图形相近则选用步骤5得到的散点图进行数据拟合;S7: According to the values of the restricted wind speed uc and the influencing factors under different fire conditions obtained in steps S5 and S6, respectively, draw the simulation results into scatter plots. If the fitted line graphs are similar, select the scatter plot obtained in step 5. perform data fitting;

S8:对步骤5得到的散点图的结果散点进行数据拟合,得到所述步骤S4中的所述无量纲计算公式中的未知系数k2的取值,将所得的uc *=k1Q**的计算值绘制散点图并进行数据拟合,得到所述步骤S4中的所述无量纲计算公式中的未知系数k1的取值,进而建立无量纲式的所述受限风速uc的计算公式:

Figure GDA0003574147540000041
S8: perform data fitting on the result scatter of the scatter diagram obtained in step 5, obtain the value of the unknown coefficient k 2 in the dimensionless calculation formula in step S4, and set the obtained u c * = k The calculated value of 1 Q ** is drawn into a scatter plot and data fitting is performed to obtain the value of the unknown coefficient k 1 in the dimensionless calculation formula in the step S4, and then establish the limit of the dimensionless formula. The calculation formula of wind speed u c :
Figure GDA0003574147540000041

在上述技术方案的基础上,本发明还可以做如下改进。On the basis of the above technical solutions, the present invention can also be improved as follows.

进一步,一种隧道侧部重点排烟受限风速的计算方法的应用,包括以下步骤:Further, an application of the method for calculating the limited wind speed of key smoke exhausting at the side of the tunnel includes the following steps:

步骤1.根据隧道侧向重点排烟火灾烟气逆流长度的计算公式L*=76(v*)-1.67(Q*)0.56(l*)-0.23及获取的火灾信息:Q、ρ0、Cp、T0、g、HD、L、l、WD计算发生火灾的隧道的烟气逆流长度L,L*=L/HD,l*=l/HD,Q*=Q/(HD 3/2g1/2WDT0Cpρ0),v*为无量纲排烟速率,其中WD为隧道宽度,v为排烟口流速,l为排烟口到火源的间距;Step 1. According to the calculation formula L * =76(v * ) -1.67 (Q * ) 0.56 (l * )- 0.23 and the obtained fire information: Q, ρ0 , C p , T 0 , g, HD , L, l, WD calculate the length L of the backflow of smoke in the fire tunnel, L * =L/ HD , l * =l/ HD , Q * =Q/ (H D 3/2 g 1/2 W D T 0 C p ρ 0 ), v * is the dimensionless smoke exhaust rate, where W D is the tunnel width, v is the flow velocity of the exhaust port, and l is the exhaust port to the fire source spacing;

步骤2.若L大于l,计算隧道侧部重点排烟受限风速uc,且使该隧道事故段两端站厅的风机同时向事故段排风形成纵向风,并使纵向风的流速大于uc,直至隧道内温度传感器或CO报警器检测出,烟气逆流长度L=l,再减低纵向风的流速至ucStep 2. If L is greater than 1, calculate the limited wind speed u c for key smoke exhaust on the side of the tunnel, and make the fans in the station halls at both ends of the accident section of the tunnel exhaust air to the accident section at the same time to form longitudinal wind, and make the velocity of the longitudinal wind greater than u c , until the temperature sensor in the tunnel or the CO alarm detects that the flue gas counterflow length L=l, and then reduce the flow velocity of the longitudinal wind to u c ;

步骤3.若隧道事故段两端站厅的风机排出的风使L小于l,则降低纵向风的流速至小于uc,直至烟气逆流长度L=l,再提高纵向风的流速至ucStep 3. If the wind discharged by the fans in the station halls at both ends of the tunnel accident section makes L less than l, reduce the flow rate of the longitudinal wind to less than uc until the length of the countercurrent flow of the flue gas L= l , and then increase the flow rate of the longitudinal wind to u c .

与现有技术相比,本发明的有益效果是:Compared with the prior art, the beneficial effects of the present invention are:

方法简单,可根据隧道实际情况设置参数,适用于采用侧部重点排烟方式的不同隧道。方法科学有效,应用了π定理和量纲分析,更具理论基础。尚无侧部重点排烟模式下受限风速的计算方法。本方法在量纲关系推导过程中考虑了主要的火源热释放速率的影响因素,得到的结果具创新性和实际工程意义。隧道火灾受限风速对高效利用隧道排烟能力产生极大的影响,通过预测计算隧道火灾受限风速能够为侧部重点排烟下隧道火灾的消防救援提供指导。The method is simple, and the parameters can be set according to the actual conditions of the tunnel, and it is suitable for different tunnels with side-focused smoke exhaust methods. The method is scientific and effective, applying the π theorem and dimensional analysis, and has a more theoretical basis. There is no calculation method for the restricted wind speed in the side-focused smoke exhaust mode. In this method, the main influencing factors of the heat release rate of the fire source are considered in the process of derivation of the dimensional relationship, and the obtained results have innovative and practical engineering significance. The limited wind speed of tunnel fire has a great impact on the efficient utilization of tunnel smoke exhaust capacity. Predicting and calculating the limited wind speed of tunnel fire can provide guidance for fire rescue in tunnel fire under the key smoke exhaust on the side.

附图说明Description of drawings

图1为本发明的隧道侧部重点排烟受限风速的计算方法流程图;Fig. 1 is the flow chart of the calculation method of the limited wind speed of key smoke exhausting at the side of the tunnel of the present invention;

图2为本发明中侧部重点排烟模式下受限风速示意图;FIG. 2 is a schematic diagram of the limited wind speed in the side key smoke exhaust mode of the present invention;

图3为本发明中建立的缩尺寸模型隧道示意图;3 is a schematic diagram of a reduced-scale model tunnel established in the present invention;

图4为本发明中受限风速与火源热释放速率的关系图;Fig. 4 is the relation diagram of restricted wind speed and fire source heat release rate in the present invention;

图5为本发明中受限风速数值模拟结果图;Fig. 5 is the numerical simulation result diagram of restricted wind speed in the present invention;

图6为FDS软件的Velocity测点直接测出诱导风速的截面图。Figure 6 is a cross-sectional view of the induced wind speed directly measured by the Velocity measuring point of the FDS software.

附图中,各标号所代表的部件列表如下:In the accompanying drawings, the parts list represented by each number is as follows:

1、单向三车道隧道,2、排烟口,3、温度测点,4、流速测点。1. One-way three-lane tunnel, 2. Smoke outlet, 3. Temperature measuring point, 4. Flow velocity measuring point.

具体实施方式Detailed ways

以下结合附图及具体实施例对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。The principles and features of the present invention will be described below with reference to the accompanying drawings and specific embodiments. The examples are only used to explain the present invention, but not to limit the scope of the present invention.

如图1至图5所示,本发明提供一种隧道侧部重点排烟火灾受限风速的计算方法,包括以下具体步骤:As shown in Figures 1 to 5, the present invention provides a method for calculating the limited wind speed of a key smoke exhaust fire on the side of a tunnel, comprising the following specific steps:

S1:确定隧道火灾场景下受限风速uc的定义:烟气逆流长度对变动的排烟口流速不敏感时对应的最小纵向诱导风速称为受限风速uc,烟气逆流长度指火源在隧道行进方向上与排烟口的距离。S1: Determine the definition of the restricted wind speed u c in the tunnel fire scenario: the minimum longitudinal induced wind speed corresponding to the flue gas backflow length is not sensitive to the fluctuating flue gas flow velocity is called the restricted wind speed uc , and the flue gas backflow length refers to the fire source The distance from the exhaust port in the direction of travel of the tunnel.

S2:确定隧道火灾场景下受限风速uc的影响因素,建立所述受限风速uc与所述影响因素的关系式:f(uc,Q,ρ0,CP,T0,g,HD)=0;S2: Determine the influencing factors of the limited wind speed uc in the tunnel fire scenario, and establish the relationship between the limited wind speed uc and the influencing factors: f ( uc ,Q,ρ 0 , C P ,T 0 ,g , HD )=0;

所述影响因素包括火源热释放速率Q、空气密度ρ0、空气定压比热Cp、空气温度T0、重力加速度g、隧道高度HD,物理参数的数量n为7,火源热释放速率Q的单位为kg.m2/s3,空气密度ρ0的单位为kg/m3,空气定压比热Cp的单位为m2/s2.K,空气温度T0的单位为K,重力加速度g单位为m/s2,隧道高度HD的单位为m。The influencing factors include fire source heat release rate Q, air density ρ 0 , air constant pressure specific heat C p , air temperature T 0 , gravitational acceleration g, tunnel height H D , the number of physical parameters n is 7, and the fire source heat The unit of release rate Q is kg.m 2 /s 3 , the unit of air density ρ 0 is kg/m 3 , the unit of air constant pressure specific heat C p is m 2 /s 2 .K, the unit of air temperature T 0 is K, the unit of gravitational acceleration g is m/ s 2 , and the unit of tunnel height HD is m.

S3:根据所述影响因素的单位确定基本量纲,将所述影响因素由所述基本量纲表示,建立所述受限风速uc与所述影响因素的量纲关系式;S3: Determine a basic dimension according to the unit of the influencing factor, represent the influencing factor by the basic dimension, and establish a dimensional relationship between the limited wind speed uc and the influencing factor;

所述步骤S3中,所述基本量纲包括质量M、时间t、长度L、温度T,基本量纲数η为4;In the step S3, the basic dimension includes mass M, time t, length L, and temperature T, and the basic dimension number η is 4;

所述步骤S3中,所述受限风速uc与所述影响因素的关系式为:In the step S3, the relationship between the limited wind speed uc and the influencing factors is:

f(Lt-1,ML2t-3,ML-3,L2t-2T-1,T,Lt-2,L)=0;f(Lt -1 , ML 2 t -3 , ML -3 , L 2 t -2 T -1 , T, Lt -2 , L)=0;

S4:确定所述影响因素的基本物理量,根据π定理得到所述影响因素的无量纲参数的数量为n-η=3,再根据π定理中循环量的选取原则,选取与长度L直接相关的变量HD,选取与温度T直接相关的变量T0,选取与时间t直接相关的变量g,选取与质量M直接相关的变量ρ0,作为循环量,用这4个循环量与其它n-η个物理参数中的其它所有参数依次组合成无量纲参数Π1、Π2、Π3,将所述步骤S3中的所述关系式转变为无量纲关系式,再获得所述受限风速uc的无量纲计算公式;S4: Determine the basic physical quantity of the influencing factor, obtain the number of dimensionless parameters of the influencing factor according to the π theorem as n-η=3, and then select the one directly related to the length L according to the selection principle of the circulation quantity in the π theorem Variable H D , select the variable T 0 directly related to the temperature T, select the variable g directly related to the time t, select the variable ρ 0 directly related to the mass M, as the circulation amount, use these 4 circulation amounts and other n- All other parameters in the n physical parameters are sequentially combined into dimensionless parameters Π 1 , Π 2 , Π 3 , the relational expression in the step S3 is converted into a dimensionless relational expression, and then the restricted wind speed u is obtained. The dimensionless calculation formula of c ;

所述步骤S4中,所述影响因素的无量纲参数Π1、Π2、Π3为:In the step S4, the dimensionless parameters Π 1 , Π 2 , and Π 3 of the influencing factors are:

Figure GDA0003574147540000071
Figure GDA0003574147540000071

所述受限风速uc的无量纲计算公式为:The dimensionless calculation formula of the restricted wind speed u c is:

Figure GDA0003574147540000072
Figure GDA0003574147540000072

Figure GDA0003574147540000081
uc *为无量纲风速,Q**为无量纲功率;which is
Figure GDA0003574147540000081
u c * is the dimensionless wind speed, Q ** is the dimensionless power;

S5:FDS数值模拟部分:建立隧道的缩尺几何模型,设置虚拟火源为丙烷燃烧,模拟不同的火源热释放速率。排烟口尺寸为0.3m×0.1m,设置排烟口流速v从小到大逐步递增,得到FDS软件拟合出的不同烟气逆流长度L,当烟气逆流长度L随排烟口流速v的递增而减少的趋势变得缓慢时,此时的排烟口流速v在隧道内产生的诱导风速为该火灾条件下的实际受限风速uc的数值,并测量记录受限风速uc的影响因素的数值,隧道的缩尺几何模型按1:10缩尺,尺寸为长×宽×高为50m×1.1m×0.45m;单向三车道隧道1内的一侧壁上沿纵长方向设有两个排烟口2,两个排烟口之间布置火源,火源位于单向三车道中间车道的中心线上,采用稳态火或t平方火,尺寸长×宽0.6m×0.2m;所述火源的功率为9.5~94.9kW,对应的全尺寸火源热释放速率为3~30MW,所述缩尺几何模型内沿顶壁的纵长方向上均匀间隔设有多个温度测点,提取FDS软件拟合出的受限风速uc的影响因素,所述缩尺几何模型的一处断面设置有多个流速测点,提取FDS软件拟合出的受限风速ucS5: FDS numerical simulation part: establish a scaled geometric model of the tunnel, set the virtual fire source as propane combustion, and simulate the heat release rate of different fire sources. The size of the exhaust port is 0.3m×0.1m, and the flow velocity v of the exhaust port is set to gradually increase from small to large, and the different flue gas counterflow lengths L fitted by the FDS software are obtained. When the trend of increasing and decreasing becomes slow, the induced wind speed in the tunnel generated by the flow velocity v of the exhaust vent at this time is the value of the actual restricted wind speed uc under the fire condition, and the influence of the restricted wind speed uc is measured and recorded. The numerical value of the factor, the scaled geometric model of the tunnel is scaled at 1:10, and the dimensions are length × width × height 50m × 1.1m × 0.45m; There are two smoke exhaust ports 2, and the fire source is arranged between the two smoke exhaust ports. The fire source is located on the center line of the middle lane of the one-way three-lane. Steady-state fire or t-square fire is adopted, and the size is 0.6m × 0.2 m; the power of the fire source is 9.5-94.9kW, the corresponding heat release rate of the full-scale fire source is 3-30MW, and the scaled geometric model is provided with a plurality of temperatures evenly spaced along the longitudinal direction of the top wall Measuring points, extract the influencing factors of the restricted wind speed uc fitted by the FDS software, a section of the scaled geometric model is provided with a plurality of flow velocity measuring points, and extract the restricted wind speed uc fitted by the FDS software;

S5:实体实验部分:建立隧道的缩尺几何模型,在缩尺几何模型中设置不同的火灾条件,包括用燃烧器作为火源,模拟不同的火源热释放速率。排烟口尺寸为0.3m×0.1m,设置不同的排烟口流速v,并用来控制不同烟气逆流长度L,并在该火灾条件下改变纵向风速的大小,并观察烟气逆流长度不受排烟口流速影响的时刻,测量该时刻的实际受限风速uc的数值,并测量记录受限风速uc的影响因素的数值,隧道的缩尺几何模型按1:10缩尺,尺寸为长×宽×高为10m×1.1m×0.45m;单向三车道隧道1内的一侧壁上沿纵长方向设有两个排烟口2,两个排烟口之间布置火源,火源位于单向三车道中间车道的中心线上,尺寸长×宽15cm×15cm;所述火源的功率为9.5~94.9kW,对应的全尺寸火源热释放速率为3~30MW。所述缩尺几何模型内沿顶壁的纵长方向上均匀间隔安装多个温度传感器,测量受限风速uc的影响因素,所述缩尺几何模型的一处断面安装风速仪,测量受限风速ucS5: Entity experiment part: establish a scaled geometric model of the tunnel, and set different fire conditions in the scaled geometric model, including using the burner as the fire source to simulate different heat release rates of the fire source. The size of the exhaust port is 0.3m×0.1m, and different flow velocity v of the exhaust port is set to control the length L of the reverse flow of the flue gas, and change the size of the longitudinal wind speed under the fire conditions, and observe that the length of the reverse flow of the flue gas is not affected. At the moment when the flow velocity of the exhaust port is affected, measure the value of the actual restricted wind speed uc at that moment, and measure and record the value of the influencing factors of the restricted wind speed uc . The scaled geometric model of the tunnel is scaled at 1:10, and the size is The length×width×height is 10m×1.1m×0.45m; one side wall of the one-way three-lane tunnel 1 is provided with two smoke exhaust ports 2 along the longitudinal direction, and the fire source is arranged between the two smoke exhaust ports. The fire source is located on the center line of the one-way three-lane middle lane, with a size of 15cm×15cm in length and width; the power of the fire source is 9.5-94.9kW, and the corresponding full-size fire source heat release rate is 3-30MW. In the scaled geometric model, a plurality of temperature sensors are installed at uniform intervals along the longitudinal direction of the top wall to measure the influencing factors of the limited wind speed u c . An anemometer is installed on a section of the scaled geometric model, and the measurement is limited wind speed u c ;

S6:根据步骤S5得到的不同火灾条件下受限风速uc和影响因素的数值,将模拟结果绘制成散点图;S6: According to the numerical values of the restricted wind speed uc and the influencing factors under different fire conditions obtained in step S5, the simulation results are drawn into a scatter diagram;

S7:对所有散点图的结果散点进行非线性拟合,得到所述步骤S4中的所述无量纲计算公式中的各未知系数k1、k2的取值,进而建立无量纲式的所述受限风速uc的计算公式:

Figure GDA0003574147540000091
S7: Perform nonlinear fitting on the resulting scatter points of all scatter plots, obtain the values of the unknown coefficients k 1 and k 2 in the dimensionless calculation formula in the step S4, and then establish a dimensionless formula The calculation formula of the restricted wind speed uc :
Figure GDA0003574147540000091

上述火灾动力模拟软件(FDS)是现有技术。The above-mentioned fire dynamic simulation software (FDS) is the prior art.

在步骤S5中,FDS数值模拟部分,单向三车道隧道1为1:10缩尺寸长×宽×高为50m×1.1m×0.45m;单向三车道隧道1内的一侧壁上设有两个排烟口2,且其内顶壁上均匀间隔设有多个温度测点3,其隧道某处断面设置有多个流速测点4。In step S5, in the FDS numerical simulation part, the one-way three-lane tunnel 1 is a 1:10 reduced size length × width × height of 50m × 1.1m × 0.45m; one side wall of the one-way three-lane tunnel 1 is provided with The two smoke exhaust ports 2 are provided with a plurality of temperature measuring points 3 evenly spaced on the inner top wall thereof, and a plurality of flow velocity measuring points 4 are arranged on a certain section of the tunnel.

在步骤S5中,实体实验部分,单向三车道隧道1为1:10缩尺寸长×宽×高为10m×1.1m×0.45m;单向三车道隧道1内的一侧壁上设有两个排烟口2,且其内顶壁上均匀间隔设有多个温度测点3,其隧道某处断面设置有多个流速测点4,每个温度测点3处通过本领域技术人员所能想到的方式安装有直径1mm,探头0.1m的k型热电偶,每个热电偶传感器均通过线路与32通道无纸记录仪连接,温度传感器检测对应温度测点3处的温度,并将对应的温度信号发送给控制器,控制器接收对应的温度信号储存。每个流速测点4处通过本领域技术人员所能想到的方式安装有流速测点,采用型号(Anemomaster KA23)的热线风速仪检测对应流速测点4处的风速,并将对应的风速信号发送给控制器,控制器接收对应的流速信号储存。上述控制器、温度传感器和风速传感器均采用现有技术,各个温度传感器、风速传感器与控制器之间的控制线路同样为现有技术。在上述条件的基础上,依次改变火源热释放速率、排烟口流速的取值形成不同的烟气逆流长度,在每个火源热释放速率、排烟口流速、烟气逆流长度的火灾条件下,改变纵向风速观察并测量得到不同火灾条件下的受限风速数值,具体模拟结果见表1:In step S5, in the physical experiment part, the one-way three-lane tunnel 1 is a 1:10 reduced size length × width × height of 10m × 1.1m × 0.45m; one side wall of the one-way three-lane tunnel 1 is provided with two There are a number of smoke exhaust ports 2, and a plurality of temperature measuring points 3 are evenly spaced on the inner top wall, and a plurality of flow velocity measuring points 4 are arranged on a certain section of the tunnel. A k-type thermocouple with a diameter of 1mm and a probe of 0.1m is installed in a conceivable way. Each thermocouple sensor is connected to a 32-channel paperless recorder through a line. The temperature signal is sent to the controller, and the controller receives the corresponding temperature signal and stores it. Each flow velocity measuring point 4 is installed with a flow velocity measuring point in a way that those skilled in the art can think of. A hot-wire anemometer of the model (Anemomaster KA23) is used to detect the wind speed at the corresponding flow velocity measuring point 4, and the corresponding wind speed signal is sent. To the controller, the controller receives the corresponding flow rate signal and stores it. The above-mentioned controller, temperature sensor and wind speed sensor all adopt the prior art, and the control circuit between each temperature sensor, wind speed sensor and the controller is also the prior art. On the basis of the above conditions, the values of the heat release rate of the fire source and the flow rate of the smoke exhaust port are changed in turn to form different flue gas backflow lengths. Under different fire conditions, the restricted wind speed values under different fire conditions were obtained by changing the longitudinal wind speed observation and measurement. The specific simulation results are shown in Table 1:

表1 FDS数值模拟结果:Table 1 FDS numerical simulation results:

Figure GDA0003574147540000101
Figure GDA0003574147540000101

模拟过程以第一组为例:第一组Q=9.5kW时,改变的排烟口流速是3m/s到11m/s,最后导致烟气逆流长度是从5.95m降到0.75m,这个下降的过程发现当烟气逆流长度达到1.05m/s时,再增大排烟口流速,烟气逆流长度L下降的很缓慢即变得不敏感。那这个1.05m对应的排烟口流速就是8m/s,因排烟口流速8m/s导致事故隧道两端均会产生一个诱导风速,这个诱导风速的值就是Q=9.5kW对应的受限风速,该诱导风速值0.51m/s是通过FDS的Velocity测点直接测出,该测点如图6所示。其余组的模拟过程类似。The simulation process takes the first group as an example: when the first group Q=9.5kW, the flue gas flow velocity changed is 3m/s to 11m/s, and finally the length of the flue gas counterflow is reduced from 5.95m to 0.75m. This decrease It is found that when the flue gas counterflow length reaches 1.05m/s, and the flow velocity of the exhaust port is increased again, the flue gas counterflow length L decreases very slowly and becomes insensitive. Then the flow rate of the smoke outlet corresponding to this 1.05m is 8m/s. Because the flow rate of the smoke outlet is 8m/s, an induced wind speed will be generated at both ends of the accident tunnel. The value of this induced wind speed is the restricted wind speed corresponding to Q=9.5kW , the induced wind speed value of 0.51m/s is directly measured by the Velocity measuring point of FDS, which is shown in Figure 6. The simulation process for the remaining groups was similar.

表2实体实验结果:Table 2 Entity experimental results:

Figure GDA0003574147540000111
Figure GDA0003574147540000111

因为实验条件有限,导致实际隧道不够长,只能做几组较小火源热释放速率的实验,而且实验的风机功率最大只能让排烟口流速达到11m/s。所以整个实验数据较少。还因为实验的时候是冬天,导致数值模拟与实验室的环境温度差别较大,最终实验结果与FDS模拟结果有误差。第一组实验模拟多次误差较小,其余组实验趋势相同。该实验验证了FDS模拟的结果正确性。根据表1和表2的数值模拟结果,绘制散点图,得到受限风速uc与各影响因素的变化关系,具体见图4,FDS模拟结果和实体实验结果趋势相同,因此得到的函数相似:Due to the limited experimental conditions, the actual tunnel is not long enough, so only a few sets of experiments with small fire source heat release rate can be done, and the maximum fan power in the experiment can only make the exhaust flow velocity reach 11m/s. So the whole experimental data is less. Also, because the experiment was in winter, there was a big difference between the numerical simulation and the ambient temperature in the laboratory, and the final experimental results were inaccurate with the FDS simulation results. The first group of experiments has a small error of multiple simulations, and the other groups of experiments have the same trend. This experiment verifies the correctness of the results of the FDS simulation. According to the numerical simulation results in Table 1 and Table 2, draw a scatter diagram to obtain the relationship between the restricted wind speed uc and various influencing factors, as shown in Figure 4. The FDS simulation results and the physical experimental results have the same trend, so the obtained functions are similar :

图4为受限风速uc与火源热释放速率Q的关系图,由图4可以获知受限风速uc随火源热释放速率Q的增大先增大后趋于稳定;Figure 4 is a graph showing the relationship between the restricted wind speed uc and the fire source heat release rate Q. From Figure 4, it can be known that the restricted wind speed uc first increases and then tends to be stable with the increase of the fire source heat release rate Q;

利用0rigin软件对图4进行数据拟合,可得uc *与Q**的函数关系为

Figure GDA0003574147540000112
即0<Q**≤0.52时,k2=1;当Q**>0.52时,k2=0;Using Origin software to fit the data in Figure 4, the functional relationship between u c * and Q ** can be obtained as
Figure GDA0003574147540000112
That is, when 0<Q ** ≤0.52, k 2 =1; when Q ** >0.52, k 2 =0;

为了确定上式中系数k1和k1'的值,将所得的uc *=k1Q**的计算值绘制于图5,由图5可见数值模拟结果均在一条直线附近波动,k1为0.34,相关系数为0.99,k1'为0.4。表明了无量纲受限风速uc *的数值计算公式和数值模拟结果的一致性。将k1=0.34代入上式可得无量纲受限风速uc *的计算公式为:In order to determine the values of the coefficients k 1 and k 1 ' in the above formula, the calculated values of u c * = k 1 Q ** are plotted in Fig. 5. It can be seen from Fig. 5 that the numerical simulation results fluctuate around a straight line, k 1 is 0.34, the correlation coefficient is 0.99, and k 1 ' is 0.4. It is shown that the numerical formula of the dimensionless confined wind speed u c * is consistent with the numerical simulation results. Substituting k 1 =0.34 into the above formula, the calculation formula of dimensionless restricted wind speed u c * is:

Figure GDA0003574147540000121
Figure GDA0003574147540000121

式中:uc *:无量纲受限风速,Q**:无量纲火源热释放速率。In the formula: u c * : Dimensionless confined wind speed, Q ** : Dimensionless fire source heat release rate.

通过此方法,可以快速得到不同火源热释放速率下受限风速的数值。Through this method, the value of the limited wind speed under different heat release rates of the fire source can be quickly obtained.

本发明的有益效果是:可根据隧道实际情况设置参数,适用于采用侧部重点排烟方式的不同隧道。方法科学有效,应用了定理和量纲分析,更具理论基础。暂无隧道火灾无量纲受限风速uc *的计算方法,本方法在量纲关系推导过程中考虑了最重要的影响因素火源热释放速率,得到的结果具创新性和实际工程意义。隧道火灾受限风速对高效利用隧道排烟能力产生极大的影响,通过预测计算隧道火灾受限风速能够为侧部重点排烟下隧道火灾的消防救援提供指导。The beneficial effect of the invention is that the parameters can be set according to the actual conditions of the tunnel, and it is suitable for different tunnels adopting the key smoke exhaust method on the side. The method is scientific and effective, applying theorems and dimensional analysis, and has a more theoretical basis. There is currently no calculation method for the dimensionless restricted wind speed u c * in tunnel fires. This method considers the most important factor in the derivation of the dimension relationship, the heat release rate of the fire source, and the obtained results are innovative and of practical engineering significance. The limited wind speed of tunnel fire has a great impact on the efficient utilization of tunnel smoke exhaust capacity. Predicting and calculating the limited wind speed of tunnel fire can provide guidance for fire rescue in tunnel fire under the key smoke exhaust on the side.

一种隧道侧部重点排烟受限风速的计算方法的应用,An application of the calculation method for the limited wind speed of key smoke exhaust on the side of the tunnel,

步骤1.根据专利CN111027176A的隧道侧向重点排烟火灾烟气逆流长度的计算公式L*=76(v*)-1.67(Q*)0.56(l*)-0.23及获取的火灾信息:Q、ρ0、Cp、T0、g、HD、L、l、WD计算发生火灾的隧道的烟气逆流长度L,L*=L/HD,l*=l/HD,Q*=Q/(HD 3/2g1/2WDT0Cpρ0),v*为无量纲排烟速率,v*=v/(HD 1/2g1/2),其中WD为隧道宽度,v为排烟口流速,l为排烟口到火源的间距;Step 1. According to the calculation formula L * =76(v * ) -1.67 (Q * ) 0.56 (l * )- 0.23 and the obtained fire information: Q, ρ 0 , C p , T 0 , g, HD , L, l, WD calculate the length L of the backflow of smoke in the tunnel where the fire occurred, L * =L/ HD , l * =l/ HD , Q * =Q/( HD 3/2 g 1/2 W D T 0 C p ρ 0 ), v * is the dimensionless smoke exhaust rate, v * =v/( HD 1/2 g 1/2 ), where WD is the width of the tunnel, v is the flow velocity of the exhaust port, and l is the distance from the exhaust port to the fire source;

步骤2.若L大于l,计算隧道侧部重点排烟受限风速uc,且使该隧道事故段两端站厅的风机同时向事故段排风形成纵向风,并使纵向风的流速大于uc,直至隧道内温度传感器或CO报警器检测出,烟气逆流长度L=l,再减低纵向风的流速至ucStep 2. If L is greater than 1, calculate the limited wind speed u c for key smoke exhaust on the side of the tunnel, and make the fans in the station halls at both ends of the accident section of the tunnel exhaust air to the accident section at the same time to form longitudinal wind, and make the velocity of the longitudinal wind greater than u c , until the temperature sensor in the tunnel or the CO alarm detects that the flue gas counterflow length L=l, and then reduce the flow velocity of the longitudinal wind to u c ;

步骤3.若隧道事故段两端站厅的风机排出的风使L小于l,则降低纵向风的流速至小于uc,直至烟气逆流长度L=l,再提高纵向风的流速至ucStep 3. If the wind discharged by the fans in the station halls at both ends of the tunnel accident section makes L less than l, reduce the flow rate of the longitudinal wind to less than uc until the length of the countercurrent flow of the flue gas L= l , and then increase the flow rate of the longitudinal wind to u c .

以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection of the present invention. within the range.

Claims (2)

1.一种隧道侧部重点排烟火灾受限风速的计算方法,包括以下具体步骤:1. A method for calculating the limited wind speed of a key smoke exhaust fire on the side of a tunnel, comprising the following specific steps: S1:确定隧道火灾场景下受限风速uc的影响因素,建立所述受限风速uc与所述影响因素的关系式:f(uc,Q,ρ0,CP,T0,g,HD)=0;S1: Determine the influencing factors of the limited wind speed uc in the tunnel fire scenario, and establish the relationship between the limited wind speed uc and the influencing factors: f ( uc ,Q,ρ 0 , C P ,T 0 ,g , HD )=0; 所述影响因素包括火源热释放速率Q、空气密度ρ0、空气定压比热Cp、空气温度T0、重力加速度g、隧道高度HD,物理参数的数量n为7,火源热释放速率Q的单位为kg.m2/s3,空气密度ρ0的单位为kg/m3,空气定压比热Cp的单位为m2/s2.K,空气温度T0的单位为K,重力加速度g单位为m/s2,隧道高度HD的单位为m;The influencing factors include fire source heat release rate Q, air density ρ 0 , air constant pressure specific heat C p , air temperature T 0 , gravitational acceleration g, tunnel height H D , the number of physical parameters n is 7, and the fire source heat The unit of release rate Q is kg.m 2 /s 3 , the unit of air density ρ 0 is kg/m 3 , the unit of air constant pressure specific heat C p is m 2 /s 2 .K, the unit of air temperature T 0 is K, the unit of gravitational acceleration g is m/ s 2 , and the unit of tunnel height HD is m; S2:根据所述影响因素的单位确定基本量纲,将所述影响因素由所述基本量纲表示,建立所述受限风速uc与所述影响因素的量纲关系式;S2: Determine a basic dimension according to the unit of the influencing factor, represent the influencing factor by the basic dimension, and establish a dimensional relationship between the limited wind speed uc and the influencing factor; S3:所述基本量纲包括质量M、时间t、长度L、温度T,基本量纲数η为4;S3: the basic dimension includes mass M, time t, length L, temperature T, and the basic dimension number η is 4; 所述步骤S2中,所述受限风速uc与所述影响因素的关系式为:In the step S2, the relationship between the limited wind speed uc and the influencing factor is: f(Lt-1,ML2t-3,ML-3,L2t-2T-1,T,Lt-2,L)=0;f(Lt -1 , ML 2 t -3 , ML -3 , L 2 t -2 T -1 , T, Lt -2 , L)=0; S4:确定所述影响因素的基本物理量,根据π定理得到所述影响因素的无量纲参数的数量为n-η=3,再根据π定理中循环量的选取原则,选取与长度L直接相关的变量HD,选取与温度T直接相关的变量T0,选取与时间t直接相关的变量g,选取与质量M直接相关的变量ρ0,作为循环量,用这4个循环量与其它n-η个物理参数中的其它所有参数依次组合成无量纲参数Π1、Π2、Π3,将所述步骤S3中的所述关系式转变为无量纲关系式,再获得所述受限风速uc的无量纲计算公式;S4: Determine the basic physical quantity of the influencing factor, obtain the number of dimensionless parameters of the influencing factor according to the π theorem as n-η=3, and then select the one directly related to the length L according to the selection principle of the circulation quantity in the π theorem Variable H D , select the variable T 0 directly related to the temperature T, select the variable g directly related to the time t, select the variable ρ 0 directly related to the mass M, as the circulation amount, use these 4 circulation amounts and other n- All other parameters in the n physical parameters are sequentially combined into dimensionless parameters Π 1 , Π 2 , Π 3 , the relational expression in the step S3 is converted into a dimensionless relational expression, and then the restricted wind speed u is obtained. The dimensionless calculation formula of c ; 所述步骤S4中,所述影响因素的无量纲参数Π1、Π2、Π3为:In the step S4, the dimensionless parameters Π 1 , Π 2 , and Π 3 of the influencing factors are:
Figure FDA0003574147530000021
Figure FDA0003574147530000021
所述受限风速uc的无量纲计算公式为:The dimensionless calculation formula of the restricted wind speed u c is:
Figure FDA0003574147530000022
Figure FDA0003574147530000022
Figure FDA0003574147530000023
uc *为无量纲风速,Q**为无量纲功率;
which is
Figure FDA0003574147530000023
u c * is the dimensionless wind speed, Q ** is the dimensionless power;
S5:FDS数值模拟部分:建立隧道的缩尺几何模型,设置虚拟火源为丙烷燃烧,模拟不同的火源热释放速率,设置排烟口流速v从小到大逐步递增,得到FDS软件拟合出的不同烟气逆流长度L,当烟气逆流长度L随排烟口流速v的递增而减少的趋势变得缓慢时,此时的排烟口流速v在隧道内产生的诱导风速为该火灾条件下的实际受限风速uc的数值,并测量记录受限风速uc的影响因素的数值,隧道内的一侧壁上沿纵长方向设有两个排烟口,两个排烟口之间布置火源,火源位于单向三车道中间车道的中心线上,采用稳态火或t平方火,所述缩尺几何模型内沿顶壁的纵长方向上均匀间隔设有多个温度测点,提取FDS软件拟合出的受限风速uc的影响因素,所述缩尺几何模型的一处断面设置有多个流速测点,提取FDS软件拟合出的受限风速ucS5: FDS numerical simulation part: establish the scaled geometric model of the tunnel, set the virtual fire source as propane combustion, simulate the heat release rate of different fire sources, set the flow velocity v of the exhaust port to increase gradually from small to large, and get the FDS software fitting When the flue gas counterflow length L decreases slowly with the increase of the flue gas flow velocity v, the induced wind speed generated in the tunnel by the flue gas flow velocity v at this time is the fire condition The value of the actual restricted wind speed uc at the lower end of the tunnel, and the values of the influencing factors of the restricted wind speed uc are measured and recorded. One side wall in the tunnel is provided with two smoke exhaust ports along the longitudinal direction, the The fire source is located on the center line of the one-way three-lane middle lane, and the steady-state fire or t-square fire is used. In the scaled geometric model, there are multiple temperature uniform intervals along the longitudinal direction of the top wall. Measuring points, extract the influencing factors of the restricted wind speed uc fitted by the FDS software, a section of the scaled geometric model is provided with a plurality of flow velocity measuring points, and extract the restricted wind speed uc fitted by the FDS software; S6:实体模拟:建立隧道的缩尺几何模型,在缩尺几何模型中设置不同的火灾条件,包括用燃烧器作为火源,模拟不同的火源热释放速率,设置不同的排烟口流速v,并用来控制不同烟气逆流长度L,并在该火灾条件下改变纵向风速的大小,并观察烟气逆流长度不受排烟口流速影响的时刻,测量该时刻的实际受限风速uc的数值,并测量记录受限风速uc的影响因素的数值,隧道内的一侧壁上沿纵长方向设有两个排烟口,两个排烟口之间布置火源,火源位于单向三车道中间车道的中心线上,所述缩尺几何模型内沿顶壁的纵长方向上均匀间隔安装多个温度传感器,测量受限风速uc的影响因素,所述缩尺几何模型的一处断面安装风速仪,测量受限风速ucS6: Entity simulation: establish a scaled geometric model of the tunnel, and set different fire conditions in the scaled geometric model, including using the burner as the fire source, simulating different heat release rates of the fire source, and setting different smoke exhaust flow velocity v , and used to control different flue gas counterflow lengths L, and change the vertical wind speed under the fire conditions, and observe the moment when the flue gas counterflow length is not affected by the flow velocity of the exhaust port, and measure the actual restricted wind speed u c at this moment. The numerical value, and the numerical value of the influencing factors of the restricted wind speed uc are measured and recorded. There are two smoke exhaust ports on one side wall in the tunnel along the longitudinal direction, and the fire source is arranged between the two smoke exhaust ports. To the center line of the middle lane of the three lanes, a plurality of temperature sensors are installed at uniform intervals along the longitudinal direction of the top wall in the scaled geometric model to measure the influencing factors of the restricted wind speed uc . An anemometer is installed on one section to measure the restricted wind speed u c ; S7:分别根据步骤S5和S6得到的不同火灾条件下受限风速uc和影响因素的数值,将模拟结果分别绘制成散点图,若拟合线图形相近则选用步骤5得到的散点图进行数据拟合;S7: According to the values of the restricted wind speed uc and the influencing factors under different fire conditions obtained in steps S5 and S6, respectively, draw the simulation results into scatter plots. If the fitted line graphs are similar, select the scatter plot obtained in step 5. perform data fitting; S8:对步骤5得到的所有散点图的结果散点进行数据拟合,得到所述步骤S4中的所述无量纲计算公式中的未知系数k2的取值,将所得的
Figure FDA0003574147530000031
的计算值绘制散点图并进行数据拟合,得到所述步骤S4中的所述无量纲计算公式中的未知系数k1的取值,进而建立无量纲式的所述受限风速uc的计算公式:
Figure FDA0003574147530000032
S8: Perform data fitting on the resulting scatter points of all the scatter plots obtained in step 5, and obtain the value of the unknown coefficient k 2 in the dimensionless calculation formula in the step S4.
Figure FDA0003574147530000031
Draw a scatter plot of the calculated value and perform data fitting to obtain the value of the unknown coefficient k 1 in the dimensionless calculation formula in the step S4, and then establish the dimensionless formula of the restricted wind speed u c Calculation formula:
Figure FDA0003574147530000032
2.如权利要求1所述的隧道侧部重点排烟受限风速的计算方法的应用,包括以下步骤:2. the application of the calculation method of the key smoke exhaust restricted wind speed of the tunnel side as claimed in claim 1, comprises the following steps: 步骤1.根据隧道侧向重点排烟火灾烟气逆流长度的计算公式L*=76(v*)-1.67(Q*)0.56(l*)-0.23及获取的火灾信息:Q、ρ0、Cp、T0、g、HD、L、l、WD计算发生火灾的隧道的烟气逆流长度L,L*=L/HD,l*=l/HD,Q*=Q/(HD 3/2g1/2WDT0Cpρ0),v*为无量纲排烟速率,其中WD为隧道宽度,v为排烟口流速,l为排烟口到火源的间距;Step 1. According to the calculation formula L * =76(v * ) -1.67 (Q * ) 0.56 (l * )- 0.23 and the obtained fire information: Q, ρ0 , C p , T 0 , g, HD , L, l, WD calculate the length L of the backflow of smoke in the fire tunnel, L * =L/ HD , l * =l/ HD , Q * =Q/ (H D 3/2 g 1/2 W D T 0 C p ρ 0 ), v * is the dimensionless smoke exhaust rate, where W D is the tunnel width, v is the flow velocity of the exhaust port, and l is the exhaust port to the fire source spacing; 步骤2.若L大于l,计算隧道侧部重点排烟受限风速uc,且使该隧道事故段两端站厅的风机同时向事故段排风形成纵向风,并使纵向风的流速大于uc,直至隧道内温度传感器或CO报警器检测出,烟气逆流长度L=l,再减低纵向风的流速至ucStep 2. If L is greater than 1, calculate the limited wind speed u c for key smoke exhaust on the side of the tunnel, and make the fans in the station halls at both ends of the accident section of the tunnel exhaust air to the accident section at the same time to form longitudinal wind, and make the velocity of the longitudinal wind greater than u c , until the temperature sensor in the tunnel or the CO alarm detects that the flue gas counterflow length L=l, and then reduce the flow velocity of the longitudinal wind to u c ; 步骤3.若隧道事故段两端站厅的风机排出的风使L小于l,则降低纵向风的流速至小于uc,直至烟气逆流长度L=l,再提高纵向风的流速至ucStep 3. If the wind discharged by the fans in the station halls at both ends of the tunnel accident section makes L less than l, reduce the flow rate of the longitudinal wind to less than uc until the length of the countercurrent flow of the flue gas L= l , and then increase the flow rate of the longitudinal wind to u c .
CN202110598123.XA 2021-05-31 2021-05-31 Calculation method and application of limited wind speed of side key smoke exhaust of tunnel Active CN113312771B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110598123.XA CN113312771B (en) 2021-05-31 2021-05-31 Calculation method and application of limited wind speed of side key smoke exhaust of tunnel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110598123.XA CN113312771B (en) 2021-05-31 2021-05-31 Calculation method and application of limited wind speed of side key smoke exhaust of tunnel

Publications (2)

Publication Number Publication Date
CN113312771A CN113312771A (en) 2021-08-27
CN113312771B true CN113312771B (en) 2022-05-20

Family

ID=77376276

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110598123.XA Active CN113312771B (en) 2021-05-31 2021-05-31 Calculation method and application of limited wind speed of side key smoke exhaust of tunnel

Country Status (1)

Country Link
CN (1) CN113312771B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118940502A (en) * 2024-07-23 2024-11-12 武汉科技大学 A prediction method for critical wind speed in an upward single slope tunnel with a water spray system
CN119226671B (en) * 2024-11-28 2025-06-13 中交第一公路勘察设计研究院有限公司 Tunnel fire source detection offset distance calculation, adaptive correction method and equipment

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108022496A (en) * 2018-01-05 2018-05-11 清华大学 Tunnel fire test device and test method of a kind of top containing crossbeam
CN110909474A (en) * 2019-11-27 2020-03-24 武汉科技大学 A method for determining the effective air volume of the tunnel shaft in the subway section
CN111027176A (en) * 2019-11-07 2020-04-17 武汉科技大学 A calculation method for the counterflow length of the smoke in the lateral key smoke exhaust fire in the tunnel
CN111075493A (en) * 2019-12-05 2020-04-28 西南交通大学 A design method for controlling wind speed of flue gas in plateau tunnel fire

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0819608D0 (en) * 2008-10-24 2008-12-03 Mosen Ltd Improved tunnel ventilation device
BR112013001871B1 (en) * 2010-07-27 2019-09-10 Josip Pavetic tunnel ventilation system under normal and fire conditions and tunnel ventilation and fire fighting under normal and fire conditions
CN111462609B (en) * 2020-04-21 2022-03-29 合肥工业大学 Fire burning and ventilation control system of spiral tunnel group

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108022496A (en) * 2018-01-05 2018-05-11 清华大学 Tunnel fire test device and test method of a kind of top containing crossbeam
CN111027176A (en) * 2019-11-07 2020-04-17 武汉科技大学 A calculation method for the counterflow length of the smoke in the lateral key smoke exhaust fire in the tunnel
CN110909474A (en) * 2019-11-27 2020-03-24 武汉科技大学 A method for determining the effective air volume of the tunnel shaft in the subway section
CN111075493A (en) * 2019-12-05 2020-04-28 西南交通大学 A design method for controlling wind speed of flue gas in plateau tunnel fire

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Study on induced airflow velocity of point smoke extraction in road tunnel fires;Xuepeng Jiang等;《Tunnelling and Underground Space Technology》;20171021;第637-643页 *
侧部点式排烟隧道火灾临界风速研究;姜学鹏等;《中国安全科学学报》;20210331;第105-111页 *
隧道火灾临界风速计算公式的比较分析;程婧园;《消防科学与技术》;20130315(第03期);第252-255页 *

Also Published As

Publication number Publication date
CN113312771A (en) 2021-08-27

Similar Documents

Publication Publication Date Title
CN113312771B (en) Calculation method and application of limited wind speed of side key smoke exhaust of tunnel
CN111027176B (en) A Calculation Method for the Backflow Length of Smoke in a Tunnel&#39;s Lateral Focused Smoke Exhaust Fire
CN106710421B (en) Tunnel fire simulation experiment device with adjustable slope under longitudinal wind effect
Gao et al. The airborne transmission of infection between flats in high-rise residential buildings: Particle simulation
CN110909474B (en) A method for determining the effective air volume of air supply in subway tunnel shafts
CN106055880A (en) Goaf spontaneous combustion prediction method
CN104657578A (en) Method for obtaining intensity of natural wind inside extra-long inclined-shaft free tunnel through meteorological data
CN107145737B (en) A Reverse Identification Algorithm of Pollution Sources in Unsteady Flow Field Based on Markov Chain
CN104834778A (en) Method for optimizing control parameters of ventilation and air conditioning system of subway station
CN106384165B (en) Smoke exhaust efficiency evaluation method based on orthogonality analysis method
CN103631995A (en) A method and system for simulating the temperature field in a vehicle compartment
CN115659805B (en) A method to establish a mathematical model of air pollution based on big data
CN113887054A (en) A linear temperature measurement inversion method for goaf based on GA-BP neural network
Abouelhamd et al. Identification of probabilistic size of breathing zone during single inhalation phase in semi-outdoor environmental scenarios
CN104933230A (en) Atmospheric pressure influence considered mine goaf temperature field simulation method
CN119005041A (en) Prediction method for flue gas countercurrent length of uplink single-slope tunnel of water adding spray system
CN109753679B (en) Visual monitoring method and system for dust accumulation blockage of air preheater
CN113482693B (en) Tunnel smoke-proof partition division simulation device and method, calculation model and application
Liu et al. Effects of baffles on heat and contaminant removal for tunnel smoke extraction: Evaluation criteria regarding removal efficiency and flow resistance
CN117291064A (en) Method and device for solving failure of ventilation network air quantity and air quality under condition of external fire
Kashef et al. Some effects on natural ventilation system for subway tunnel fires
CN206726624U (en) A kind of emulation platform of Highway Tunnel Ventilation System
CN206038344U (en) Exhaust port heat extraction efficient measurement system
Lanchava et al. Fire development study on physical models of transport tunnels
Langmans et al. Experimental analysis of cavity ventilation behind residential rainscreen cladding systems

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