CN105569709B - The control method of High-geotemperature railway tunnel cooling ventilation - Google Patents
The control method of High-geotemperature railway tunnel cooling ventilation Download PDFInfo
- Publication number
- CN105569709B CN105569709B CN201610049329.6A CN201610049329A CN105569709B CN 105569709 B CN105569709 B CN 105569709B CN 201610049329 A CN201610049329 A CN 201610049329A CN 105569709 B CN105569709 B CN 105569709B
- Authority
- CN
- China
- Prior art keywords
- mrow
- tunnel
- msub
- mfrac
- temperature
- 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
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F1/00—Ventilation of mines or tunnels; Distribution of ventilating currents
- E21F1/003—Ventilation of traffic tunnels
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D27/00—Control, e.g. regulation, of pumps, pumping installations or pumping systems specially adapted for elastic fluids
Abstract
The invention discloses a kind of control method of High-geotemperature railway tunnel cooling ventilation; scheme provided by the invention; it is not only that the aeration-cooling scheme of tunnel machine proposes reliable foundation; the grain cooling machine low-temperature grain-storage technology mode proposed simultaneously using this programme, it can reach the purpose for reducing blower fan configuration and protection Fan Equipment：The continuous ventilating mode of existing tunnel blower fan use is substituted in the situation without train passage by fan operation, can effectively eliminate in the case that blower fan has train passage in tunnel and be risk of damage to；The high-power machinery used relative to existing tunnel is divulged information, and this programme can effectively avoid the electromechanical equipment in tunnel from being left unused for a long time, that is, reaches the purpose for improving ventilation equipment utilization rate in tunnel.Meanwhile this programme proposes one kind and realized in tunnel ventilation temperature-fall period, the minimum aeration control method of blower fan consumption power, using this method, the power consumption of tunnel ventilation cooling can be effectively reduced, reaches the purpose for saving railway traffic operating cost.
Description
Technical field
The present invention relates to length to be more than railway tunnel falling temperature technique means field, more particularly to a kind of High-geotemperature railway tunnel
The control method of road cooling ventilation.
Background technology
When tunnel primitive rock temperature is up to 35 DEG C or so, when relative humidity is up to more than 80%, the high-geotemperature problem in tunnel is non-
Chang Yanchong；In long-term operation, if effective temperature-reducing can not be carried out effectively, it will also deteriorate operating environment.Therefore, High-geotemperature tunnel
The influence that Climatic regionalization ensures the long-term operating environment in tunnel is the technical problem of urgent need to resolve.
It is a kind of effective temperature-reducing method to being aerated in tunnel,《Railway Tunnels Ventilation design specification》
(TB10068-2010) provided in：Passenger-cargo co electrified railway tunnelses of the length more than 15Km should set force ventilation.At present
Electric railway operation ventilation is often taken a breath in " Window time " using longitudinal ventilation mode, and it is by being installed on tunnel
Interior blower fan makes formation air-flow in tunnel be taken a breath.And in practical operation, for High-geotemperature super long tunnel (such as length of tunnel
More than 20Km), section " Window time " period, the interior mechanical blowing as caused by blower fan was tended not to temperature distinguished and admirable in tunnel
As defined in dropping in scope (less than 28 DEG C), operating environment in tunnel is caused not reach code requirement.
For the cooling problem in above-mentioned tunnel, solution of the prior art is：The ventilation improved in Window time is changed
Gas number；To tunnel by the way of continuous ventilating.The realization of the first solution is dependent on machine more in tunnel above
Electric equipment, so, electromechanical equipment input cost height and maintenance cost for tunnel cooling will necessarily be caused greatly while in majority
In the case of these electromechanical equipments be idle, be unfavorable for the utilization rate of equipment and installations of railway annex；The second way is in train process of passing through tunnel
When, particularly produced by the travel direction of train and blower fan air-flow it is in opposite direction when, the generation of running at high speed of train in tunnel
Piston Action Wind extreme influence will be produced to fan life, or even can directly damage blower fan.
The content of the invention
The problem of for the above-mentioned control method for lowering temp Shortcomings of tunnel ventilation in the prior art, the invention provides one kind
High-geotemperature railway tunnel cooling ventilation control method, by method proposes blower fan running scheme, can solve High-geotemperature
Tunnel because be unable to caused by draught capacity deficiency effective temperature-reducing, draught capacity it is too big caused by electromechanical equipment is idle, continuous ventilating
Caused by blower fan the problems such as damaging, realize by temperature in tunnel drop to as defined in scope, and reach and reduce blower fan configuration and protect
Protect the purpose of Fan Equipment.
In view of the above-mentioned problems, the control method of High-geotemperature railway tunnel cooling ventilation provided by the invention passes through following technology
Main points achieve the goal：The control method of High-geotemperature railway tunnel cooling ventilation, including the following steps that order is carried out：
S1：Daily in the case of without Piston Action Wind, operation is arranged at the blower fan in tunnel, makes to produce a certain speed in tunnel
The distinguished and admirable of angle value is aerated cooling to tunnel, and sets the distinguished and admirable duration in one day, if the tunnel after fan operation N
Temperature in road is reduced to a certain setting value, then record this it is aeration-cooling under the conditions of meet the ventilation cycle t that cooling requires, it is described
N is the integer more than or equal to 2；
S2：Calculate the ventilation energy consumption e of ventilation cycle t inner blowers；
S3：Obtain ventilation cycle matrix and ventilation energy consumption matrix：Repeat step S1 and S2, record meet the m that cooling requires
It is individual it is aeration-cooling under the conditions of ventilation cycle matrix T=(t1、t2..., tm) and each cycle ventilation energy consumption matrix E=(e1、
e2..., em), and in each S1 steps, distinguished and admirable velocity amplitude and/or duration numerical value distinguished and admirable in one day；
S4：Obtain the final ventilation scheme in tunnel：The ventilation cycle matrix T=(t that S3 steps are obtained1、t2..., tm) and
Divulge information energy consumption matrix E=(e1、e2..., em) be multiplied, obtain Ventilation Control discrimination matrix TE=(t1×e1、t2×e2..., tm×
em), take final ventilation scheme of the aeration-cooling condition corresponding to least member as tunnel in TE.
Above scheme is aimed in the case of without Piston Action Wind, is found and is met that tunnel cooling requires force ventilation scheme
Minimum power consumption.Above scheme can be put into practice on the spot in tunnel：The time and the wind within the time that setting blower ventilation is run daily
After speed, after verifying fan operation for a period of time, whether the temperature in tunnel is reduced to predetermined value, if meeting, cooling requires,
Record meets power when total days running, daily run time and the operation for the blower fan that cooling requires, above power and fortune
Value of the product of row time as Ventilation Control discrimination matrix, blower fan corresponding to the minimum value in above Ventilation Control discrimination matrix
Operational factor, as meet the final ventilation scheme of tunnel ventilation cooling conditions.
Because the air themperature that above force ventilation introduces into tunnel is influenceed by season, for the benefit of fan operation parameter
The accuracy of judgement, therefore the temperature in tunnel after fan operation N is reduced to a certain setting value and is set as fan operation parameter
Judge the cycle, so, preferably above ventilation cycle matrix T is remembered with year, then divulge information energy consumption matrix in parameter represent 1 year when
Between in the gross energy that is consumed of a certain fan parameter.
By above scheme, only the aeration-cooling scheme of tunnel machine does not propose reliable foundation, while using we
The grain cooling machine low-temperature grain-storage technology mode that case proposes, it can reach the purpose for reducing blower fan configuration and protection Fan Equipment：Pass through fan operation
The continuous ventilating mode of existing tunnel blower fan use is substituted in the situation without train passage, can effectively eliminate blower fan has in tunnel
It is risk of damage in the case of train passage；The high-power machinery used relative to existing tunnel is divulged information, and this programme can be effective
Avoid the electromechanical equipment in tunnel from being left unused for a long time, that is, reach the purpose for improving ventilation equipment utilization rate in tunnel.Meanwhile this
Scheme proposes one kind and realized in tunnel ventilation temperature-fall period, the minimum aeration control method of blower fan consumption power, using we
Method, the power consumption of tunnel ventilation cooling can be effectively reduced, reach the purpose for saving railway traffic operating cost.
As the further technical scheme of above control method：
The Window time run with reference to existing China railways and the cooling requirement of various regions, meet that tunnel cools for ease of obtaining
It is required that fan operation parameter, for distinguished and admirable velocity amplitude between 2 meter per seconds between 10 meter per seconds, the distinguished and admirable duration is small between 2
When/day is between 8 hours/day.
Railway tunnel has the characteristics that：In longitudinal depth in tunnel, tunnel wall each point temperature value does not wait, blower fan draws
The mechanical blowing original temperature entered is influenceed by seasonal climate, on tunnel longitudinal length, the heat transfer of tunnel wall each point and mechanical blowing
Speed, thus to shorten the cycle that optimal ventilation scheme obtains, can be by way of Computing, the control method is adopted
Carried out with computer simulation, before computer simulation, tunnel curve of ground temperature data are inputted into computer, year, temperature became outside hole
Change curve data, m distinguished and admirable velocity amplitude values, m distinguished and admirable duration values；
The computer simulation is carried out in the following way：By the convection heat transfer' heat-transfer by convection system of the tunnel wall obtained on the spot and air
The heat transfer coefficient input computer of number, country rock, using tunnel curve of ground temperature data as the original temperature value of tunnel wall, with hole outside
Year temperature Change curve data as tunnel portal air themperature boundary value, with the convective heat-transfer coefficient of tunnel wall and air,
The airflow velocity value at each time point is enclosed as cold and hot convection heat transfer' heat-transfer by convection relative velocity, the tunnel at each time point of air and tunnel wall
The difference of rock temperature value and air temperature value is that thermograde obtains heat output, and the heat output in N is summed, obtains N
Force ventilation heat output in year, finally obtain tunnel wall surface temperature drop.
In above method, if the airflow velocity of mechanical blowing is more than certain value, the heat output of tunnel wall and air is more than
Country rock so may be such that in tunnel that temperature declines, and is carried out using computer simulation to the heat output of tunnel wall, can be by tunnel
Wall is divided into multiple points, in ventilation process, to establish each point temperature for each point and change with time relation, pass through meter
The computing of calculation machine, the fan parameter for meeting that cooling requires is obtained in a short time.
Because three principal elements for influenceing tunnel temperature are respectively：Country rock temperature, blower fan force ventilation heat transfer, train fortune
Row Piston Action Wind is conducted heat, therefore Piston Action Wind heat transfer is considered in total heat transfer, can obtain the fan parameter simulation closer to actual value
Value, therefore：Also include the heat transfer conversion that driving produces Piston Action Wind, i.e. tunnel Piston Action Wind Calculation of Heat Transfer, the tunnel Piston Action Wind heat transfer
It is calculated as：Piston Action Wind curve data is inputted into computer, the tunnel wind speed angle value at each time point is as air and tunnel wall
The difference of cold and hot convection heat transfer' heat-transfer by convection relative velocity, the tunnel wall temperature value at each time point and air temperature value is obtained for thermograde
Piston Action Wind heat output, the Piston Action Wind heat output in N is summed, obtain the tunnel Piston Action Wind heat output in N, tunnel is lived
Hearsay heat and total heat transfer force ventilation heat output and that value is as tunnel palisades are filled in, finally obtains the final of tunnel wall
Temperature, the final temperature value are the reference value of a certain setting value.
Outside tunnel curve of ground temperature data, hole year temperature Change curve data and Piston Action Wind curve data acquisition pattern be as follows
Mode：
The tunnel curve of ground temperature data for survey the tunnel longitudinal direction ground temperature curve data that borehole data is drafted；
Temperature Change curve data was drawn and obtained by following formula year outside the hole：
In formula, taFor local year air-flow mean temperature；tmFor gas flow temperature amplitude, τ 1 is outside hole on year temperature Change curve
Certain point is in hourage residing then；
The Piston Action Wind curve data is drawn and obtained by following formula：
Wherein：V=vmWhen situation is that vehicle driving direction is consistent with blower ventilation direction；
V=-vmWhen situation is vehicle driving direction and blower ventilation in opposite direction；
When v=0 situations is in tunnels without vehicle.
Because length of tunnel is longer, when carrying out the conversion of tunnel surrounding temperature drop, along the radial direction in tunnel, wind-warm syndrome and enclose
Rock temperature can linear change, therefore set up complete mathematical description, be easy to more easily receive more accurate blower fan wind speed
And ventilation time, meet the blower fan wind speed of tunnel cooling conditions and the specific implementation of ventilation time as another acquisition,
Acquisition meets that the blower fan wind speed of tunnel cooling conditions and ventilation time are realized in the following way：By tunnel longitudinal direction unsteady state temperature
The governing equation (1) of field country rock part, governing equation (2) the composition tunnel longitudinal direction of the tunnel longitudinal direction distinguished and admirable part of Unsteady Temperature Field
The complete mathematical description of three-dimensional unstable state heat transfer：
The complete asynchronous long display difference equation of three-dimensional unstable state heat transfer is established using finite difference calculus, it is as follows：
The temperature of country rock nodeApproximate expression formula：
The temperature T of convective boundary node 1j,n+1Approximate expression formula：
The longitudinal distinguished and admirable approximate temperature in tunnelExpression formula：
To make Equations stable, equation (3), the calculating time step Δ τ of equation (4) and tunnel longitudinal space step-length
Δ x should meet equation (6) and equation (7) respectively：
More than it is various in, τ be setting time span, rjFor the radial distance of j-th of radial point in tunnel, R is tunnel
Radius is excavated, T is the temperature in country rock, and r is radial distance of any point along tunnel in country rock, and x is any point in country rock
Along tunnel longitudinal direction distance, ρ be air-flow density, A be tunnel area of section, TfFor tunnel gas flow temperature, TbFor tunnel wall
Temperature, v are airflow velocity, qsFor heat caused by equipment in tunnel, h is convection transfer rate, and U is tunnel cross section girth；
(i, j) be along i-th of tunnel axis direction, radially j-th of grid node,For country rock node (i, j) the n moment temperature
(℃)；Δ r is tunnel cross section radial distance step-length；Δ x is tunnel fore-and-aft distance step-length；Δ τ is time step；A is country rock
Thermal conductivity factor；FoFor Fourier number；cpThe specific heat at constant pressure of air-flow；For the wind-warm syndrome of i-node in n moment tunnels；
It it is 1 year by a ventilation cycle of temperature prediction in tunnel, using longitudinal ground temperature curve number value in tunnel as calculating
Start time tunnel surrounding longitudinal temperature initial value, year, temperature Change curve numerical value entered as a ventilation cycle inner tunnel outside hole
One's intention as revealed in what one says stream temperature initial value, Piston Action Wind change curve influences airflow velocity as driving in daily tunnel, during by daily ventilation
Between and corresponding airflow velocity, substitute into above difference equation, by equation (6) and equation (7) determination time step Δ τ and space
Step delta x, difference equation (3), equation (4), equation (5) are solved, obtain different year tunnel surrounding and airflow temperature
Along the distribution of longitudinal direction, temperature prediction in tunnel is carried out, obtains meeting the ventilation cycle t that cooling requires under the wind speed；
Outside tunnel curve of ground temperature data, hole year temperature Change curve data and Piston Action Wind curve data acquisition pattern be as follows
Mode：
The tunnel curve of ground temperature data for survey the tunnel longitudinal direction ground temperature curve data that borehole data is drafted；
Temperature Change curve data was drawn and obtained by following formula year outside the hole：
In formula, taFor local year air-flow mean temperature；tmFor gas flow temperature amplitude, τ 1 is outside hole on year temperature Change curve
Certain point is in hourage residing then；
The Piston Action Wind curve data is drawn and obtained by following formula：
Wherein v=vmWhen situation is that vehicle driving direction is consistent with blower ventilation direction；
V=-vmWhen situation is vehicle driving direction and blower ventilation in opposite direction；
When v=0 situations is in tunnels without vehicle.
As the specific implementation for being substantially shorter the solution time, area of computer aided can be used to check, you can by right
Difference equation (3), equation (4), equation (5) solve by the way of software MATLAB progress.
The invention has the advantages that：
Scheme provided by the invention, not the only aeration-cooling scheme of tunnel machine propose reliable foundation, use simultaneously
The grain cooling machine low-temperature grain-storage technology mode that this programme proposes, it can reach the purpose for reducing blower fan configuration and protection Fan Equipment：Pass through blower fan
The situation for running on no train passage substitutes the continuous ventilating mode that existing tunnel blower fan uses, and can effectively eliminate blower fan in tunnel
It is risk of damage in the case of inside having train passage；The high-power machinery used relative to existing tunnel is divulged information, and this programme can
Effectively avoid the electromechanical equipment in tunnel from being left unused for a long time, that is, reach the purpose for improving ventilation equipment utilization rate in tunnel.Together
When, this programme proposes one kind and realized in tunnel ventilation temperature-fall period, the minimum aeration control method of blower fan consumption power, uses
This method, the power consumption of tunnel ventilation cooling can be effectively reduced, reach the purpose for saving railway traffic operating cost.
Brief description of the drawings
Fig. 1 is the realization of one specific embodiment of control method of High-geotemperature railway tunnel cooling ventilation of the present invention
Procedural model；
Fig. 2 is the flow of one specific embodiment of control method of High-geotemperature railway tunnel cooling ventilation of the present invention
Figure.
Embodiment
With reference to embodiment, the present invention is described in further detail, but the present invention structure be not limited only to it is following
Embodiment.
Embodiment 1：
As shown in Fig. 2 the control method of High-geotemperature railway tunnel cooling ventilation, including the following steps that order is carried out：
S1：Daily in the case of without Piston Action Wind, operation is arranged at the blower fan in tunnel, makes to produce a certain speed in tunnel
The distinguished and admirable of angle value is aerated cooling to tunnel, and sets the distinguished and admirable duration in one day, if the tunnel after fan operation N
Temperature in road is reduced to a certain setting value, then record this it is aeration-cooling under the conditions of meet the ventilation cycle t that cooling requires, it is described
N is the integer more than or equal to 2；
S2：Calculate the ventilation energy consumption e of ventilation cycle t inner blowers；
S3：Obtain ventilation cycle matrix and ventilation energy consumption matrix：Repeat step S1 and S2, record meet the m that cooling requires
It is individual it is aeration-cooling under the conditions of ventilation cycle matrix T=(t1、t2..., tm) and each cycle ventilation energy consumption matrix E=(e1、
e2..., em), and in each S1 steps, distinguished and admirable velocity amplitude and/or duration numerical value distinguished and admirable in one day；
S4：Obtain the final ventilation scheme in tunnel：The ventilation cycle matrix T=(t that S3 steps are obtained1、t2..., tm) and
Divulge information energy consumption matrix E=(e1、e2..., em) be multiplied, obtain Ventilation Control discrimination matrix TE=(t1×e1、t2×e2..., tm×
em), take final ventilation scheme of the aeration-cooling condition corresponding to least member as tunnel in TE.
Above scheme is aimed in the case of without Piston Action Wind, is found and is met that tunnel cooling requires force ventilation scheme
Minimum power consumption.Above scheme can be put into practice on the spot in tunnel：The time and the wind within the time that setting blower ventilation is run daily
After speed, after verifying fan operation for a period of time, whether the temperature in tunnel is reduced to predetermined value, if meeting, cooling requires,
Record meets power when total days running, daily run time and the operation for the blower fan that cooling requires, above power and fortune
Value of the product of row time as Ventilation Control discrimination matrix, blower fan corresponding to the minimum value in above Ventilation Control discrimination matrix
Operational factor, as meet the final ventilation scheme of tunnel ventilation cooling conditions.
Because the air themperature that above force ventilation introduces into tunnel is influenceed by season, for the benefit of fan operation parameter
The accuracy of judgement, therefore the temperature in tunnel after fan operation N is reduced to a certain setting value and is set as fan operation parameter
Judge the cycle, so, preferably above ventilation cycle matrix T is remembered with year, then divulge information energy consumption matrix in parameter represent 1 year when
Between in the gross energy that is consumed of a certain fan parameter.
By above scheme, only the aeration-cooling scheme of tunnel machine does not propose reliable foundation, while using we
The grain cooling machine low-temperature grain-storage technology mode that case proposes, it can reach the purpose for reducing blower fan configuration and protection Fan Equipment：Pass through fan operation
The continuous ventilating mode of existing tunnel blower fan use is substituted in the situation without train passage, can effectively eliminate blower fan has in tunnel
It is risk of damage in the case of train passage；The high-power machinery used relative to existing tunnel is divulged information, and this programme can be effective
Avoid the electromechanical equipment in tunnel from being left unused for a long time, that is, reach the purpose for improving ventilation equipment utilization rate in tunnel.Meanwhile this
Scheme proposes one kind and realized in tunnel ventilation temperature-fall period, the minimum aeration control method of blower fan consumption power, using we
Method, the power consumption of tunnel ventilation cooling can be effectively reduced, reach the purpose for saving railway traffic operating cost.
In the present embodiment, the N is equal with t values.
Embodiment 2：
The present embodiment is further qualified on the basis of embodiment 1：As the further technical side of above control method
Case：
The Window time run with reference to existing China railways and the cooling requirement of various regions, meet that tunnel cools for ease of obtaining
It is required that fan operation parameter, for distinguished and admirable velocity amplitude between 2 meter per seconds between 10 meter per seconds, the distinguished and admirable duration is small between 2
When/day is between 8 hours/day.
Embodiment 3：
The present embodiment is further qualified on the basis of embodiment 1：Railway tunnel has the characteristics that：In the vertical of tunnel
To in depth, the mechanical blowing original temperature that tunnel wall each point temperature value does not wait, blower fan introduces is influenceed by seasonal climate, in tunnel
On longitudinal length, the heat transfer rate of tunnel wall each point and mechanical blowing, therefore to shorten the cycle that optimal ventilation scheme obtains,
Can be by way of Computing, the control method is carried out using computer simulation, before computer simulation, to calculating
It is individual distinguished and admirable that tunnel curve of ground temperature data, the outer year temperature Change curve data in hole, m distinguished and admirable velocity amplitude values, m are inputted in machine
Duration value；
The computer simulation is carried out in the following way：By the convection heat transfer' heat-transfer by convection system of the tunnel wall obtained on the spot and air
The heat transfer coefficient input computer of number, country rock, using tunnel curve of ground temperature data as the original temperature value of tunnel wall, with hole outside
Year temperature Change curve data as tunnel portal air themperature boundary value, with the convective heat-transfer coefficient of tunnel wall and air,
The airflow velocity value at each time point is enclosed as cold and hot convection heat transfer' heat-transfer by convection relative velocity, the tunnel at each time point of air and tunnel wall
The difference of rock temperature value and air temperature value is that thermograde obtains heat output, and the heat output in N is summed, obtains N
Force ventilation heat output in year, finally obtain tunnel wall surface temperature drop.
In above method, if the airflow velocity of mechanical blowing is more than certain value, the heat output of tunnel wall and air is more than
Country rock so may be such that in tunnel that temperature declines, and is carried out using computer simulation to the heat output of tunnel wall, can be by tunnel
Wall is divided into multiple points, in ventilation process, to establish each point temperature for each point and change with time relation, pass through meter
The computing of calculation machine, the fan parameter for meeting that cooling requires is obtained in a short time.
Further, because three principal elements for influenceing tunnel temperature are respectively：Country rock temperature, blower fan force ventilation pass
Heat, the heat transfer of train operation Piston Action Wind, therefore Piston Action Wind heat transfer is considered in total heat transfer, the blower fan closer to actual value can be obtained
Parameter simulation value, therefore：Also include the heat transfer conversion that driving produces Piston Action Wind, i.e. tunnel Piston Action Wind Calculation of Heat Transfer, the tunnel is lived
Plug hearsay heat is calculated as：Piston Action Wind curve data is inputted into computer, the tunnel wind speed angle value at each time point is as air and tunnel
The difference of the cold and hot convection heat transfer' heat-transfer by convection relative velocity of road wall, the tunnel wall temperature value at each time point and air temperature value is temperature
Gradient obtains Piston Action Wind heat output, and the Piston Action Wind heat output in N is summed, and obtains the tunnel Piston Action Wind heat transfer in N
Amount, tunnel Piston Action Wind heat output and total heat transfer force ventilation heat output and that value is as tunnel palisades, finally obtain tunnel
The final temperature of wall, the final temperature value are the reference value of a certain setting value.
Outside tunnel curve of ground temperature data, hole year temperature Change curve data and Piston Action Wind curve data acquisition pattern be as follows
Mode：
The tunnel curve of ground temperature data for survey the tunnel longitudinal direction ground temperature curve data that borehole data is drafted；
Temperature Change curve data was drawn and obtained by following formula year outside the hole：
In formula, taFor local year air-flow mean temperature；tmFor gas flow temperature amplitude, τ 1 is outside hole on year temperature Change curve
Certain point is in hourage residing then；
The Piston Action Wind curve data is drawn and obtained by following formula：
Wherein：V=vmWhen situation is that vehicle driving direction is consistent with blower ventilation direction；
V=-vmWhen situation is vehicle driving direction and blower ventilation in opposite direction；
When v=0 situations is in tunnels without vehicle.
Embodiment 4：
The present embodiment is further qualified on the basis of embodiment 1：Because length of tunnel is longer, tunnel surrounding is being carried out
When temperature drop converts, along the radial direction in tunnel, wind-warm syndrome and country rock temperature can linear change, therefore set up complete mathematics and retouch
State, be easy to more easily receive more accurate blower fan wind speed and ventilation time, meet tunnel cooling bar as another kind acquisition
The blower fan wind speed of part and the specific implementation of ventilation time, obtain the blower fan wind speed and ventilation time for meeting tunnel cooling conditions
Realize in the following way：By the governing equation (1) of tunnel longitudinal direction Unsteady Temperature Field country rock part, tunnel longitudinal direction unstable state temperature
Spend the complete mathematical description of governing equation (2) composition tunnel longitudinal direction three-dimensional unstable state heat transfer of the distinguished and admirable part in field：
The complete asynchronous long display difference equation of three-dimensional unstable state heat transfer is established using finite difference calculus, it is as follows：
The temperature of country rock nodeApproximate expression formula：
The temperature T of convective boundary node 1j,n+1Approximate expression formula：
The longitudinal distinguished and admirable approximate temperature in tunnelExpression formula：
To make Equations stable, equation (3), the calculating time step Δ τ of equation (4) and tunnel longitudinal space step-length
Δ x should meet equation (6) and equation (7) respectively：
More than it is various in, τ be setting time span, rjFor the radial distance of j-th of radial point in tunnel, R is tunnel
Radius is excavated, T is the temperature in country rock, and r is radial distance of any point along tunnel in country rock, and x is any point in country rock
Along tunnel longitudinal direction distance, ρ be air-flow density, A be tunnel area of section, TfFor tunnel gas flow temperature, TbFor tunnel wall
Temperature, v are airflow velocity, qsFor heat caused by equipment in tunnel, h is convection transfer rate, and U is tunnel cross section girth；
(i, j) be along i-th of tunnel axis direction, radially j-th of grid node,For country rock node (i, j) the n moment temperature
(℃)；Δ r is tunnel cross section radial distance step-length；Δ x is tunnel fore-and-aft distance step-length；Δ τ is time step；A is country rock
Thermal conductivity factor；FoFor Fourier number；cpThe specific heat at constant pressure of air-flow；For the wind-warm syndrome of i-node in n moment tunnels；
It it is 1 year by a ventilation cycle of temperature prediction in tunnel, using longitudinal ground temperature curve number value in tunnel as calculating
Start time tunnel surrounding longitudinal temperature initial value, year, temperature Change curve numerical value entered as a ventilation cycle inner tunnel outside hole
One's intention as revealed in what one says stream temperature initial value, Piston Action Wind change curve influences airflow velocity as driving in daily tunnel, during by daily ventilation
Between and corresponding airflow velocity, substitute into above difference equation, by equation (6) and equation (7) determination time step Δ τ and space
Step delta x, difference equation (3), equation (4), equation (5) are solved, obtain different year tunnel surrounding and airflow temperature
Along the distribution of longitudinal direction, temperature prediction in tunnel is carried out, obtains meeting the ventilation cycle t that cooling requires under the wind speed；
Outside tunnel curve of ground temperature data, hole year temperature Change curve data and Piston Action Wind curve data acquisition pattern be as follows
Mode：
The tunnel curve of ground temperature data for survey the tunnel longitudinal direction ground temperature curve data that borehole data is drafted；
Temperature Change curve data was drawn and obtained by following formula year outside the hole：
In formula, taFor local year air-flow mean temperature；tmFor gas flow temperature amplitude, τ 1 is outside hole on year temperature Change curve
Certain point is in hourage residing then；
The Piston Action Wind curve data is drawn and obtained by following formula：
Wherein v=vmWhen situation is that vehicle driving direction is consistent with blower ventilation direction；
V=-vmWhen situation is vehicle driving direction and blower ventilation in opposite direction；
When v=0 situations is in tunnels without vehicle.
As the specific implementation for being substantially shorter the solution time, area of computer aided can be used to check, you can by right
Difference equation (3), equation (4), equation (5) solve by the way of software MATLAB progress.
Such as Fig. 1, in the present embodiment, tunnel ventilation mode and fan operation parameter are ventilation model, tunnel curve of ground temperature
Outside data, hole year temperature Change curve data and Piston Action Wind curve data be it needs to be determined that parameter, heat conducting module, convection current-
Heat conducting module is country rock, the conduction model of tunnel wall, in ventilation process, obtains tunnel temperature field distribution, evaluation analysis mould
Block, that is, it is used to evaluate airflow temperature, energy consumption analysis etc. in the feasibility and tunnel of fan operation parameter, output module is used to export
The operational mode for the force ventilation that this control method finally determines.
Embodiment 5：
For ease of understanding, present embodiments provide and a kind of specifically obtain the implementation of blower fan optimized parameter：The manner
For a high temperature fracture with temperature up to 70 DEG C, length up to 34.5km, two-direction road 17 it is right/the ventilation wind in the High-geotemperature tunnel in day
The specific result of calculation of the opening time of speed and blower fan：
Calculated according to passenger-cargo carriage driving average speed by 90km/h, the running time in daily tunnel is about 13 hours.Examine
Consider two-direction road in tunnel early stage, two-way frequently Piston Action Wind is serious to axial flow blower performance impact, once blower fan work enters not
Stable region hump zone, shake is kicked into generation, blower fan will damage when serious；For jet blower, two-way frequent row in the case of unlatching
Car, also result in its damage.Thus the ventilation time not influenceed daily by driving only 11 hours, i.e., 660 minutes.
In this programme, record meets 5 aeration-cooling conditions that cooling requires, according to organization of driving in tunnel, by wind speed
Taken respectively 360 minutes, 540 minutes for the daily ventilation times of 2.4m/s；270 points are taken respectively by the daily ventilation times of wind speed 3.5m/s
Clock, 360 minutes；Taken 270 minutes by the daily ventilation times of wind speed 4.1m/s, the corresponding temperature for forming wind speed of being divulged information in 5 kinds of tunnels
Prediction, it is T=(10,5,4,3,2) to obtain the ventilation cycle matrix that temperature in tunnel is dropped to 28 DEG C by satisfaction under different wind speed, should
Each numeral in ventilation cycle matrix represents. how many years；Simultaneously by recording the ventilation energy consumption matrix in each cycle
(1502340,2253510,1662210,2216280,2016990), in the ventilation energy consumption matrix, each numerical value represents different logical
The electric energy that blower fan consumes every year under the conditions of wind, unit are kilowatt；Ventilation Control discrimination matrix for (15023400,11267550,
6648840,6648840,4033980), thus judge to open axial flow blower and jet blower 270 daily for 4.1m/s by wind speed
Minute ventilation is to meet the optimal selection of cooling ventilation.
Temperature in tunnel can be dropped to 28 DEG C by the selection in 2 years, and can guarantee that axle stream and jet blower in tunnel
The running interval time opens, and is effectively protected blower fan, and ventilation energy consumption is also minimum.
Above content is to combine the further description that specific preferred embodiment is made to the present invention, it is impossible to assert this
The embodiment of invention is confined to these explanations.For general technical staff of the technical field of the invention,
The other embodiment drawn in the case where not departing from technical scheme, should be included within the scope of the present invention.
Claims (7)
1. the control method of High-geotemperature railway tunnel cooling ventilation, it is characterised in that the following steps carried out including order：
S1：Daily in the case of without Piston Action Wind, operation is arranged at the blower fan in tunnel, makes to produce a certain velocity amplitude in tunnel
It is distinguished and admirable be aerated cooling to tunnel, and the distinguished and admirable duration in one day is set, if after fan operation N in tunnel
Temperature be reduced to a certain setting value, then record this it is aeration-cooling under the conditions of meet cooling require ventilation cycle t, the N be
Integer more than or equal to 2；
S2：Calculate the ventilation energy consumption e of ventilation cycle t inner blowers；
S3：Obtain ventilation cycle matrix and ventilation energy consumption matrix：Repeat step S1 and S2, record meet that the m that cooling requires is individual logical
Ventilation cycle matrix T=(t under wind cooling conditions1、t2..., tm) and each cycle ventilation energy consumption matrix E=(e1、e2...,
em), and in each S1 steps, distinguished and admirable velocity amplitude and/or duration numerical value distinguished and admirable in one day；
S4：Obtain the final ventilation scheme in tunnel：The ventilation cycle matrix T=(t that S3 steps are obtained1、t2..., tm) and ventilation
Energy consumption matrix E=(e1、e2..., em) be multiplied, obtain Ventilation Control discrimination matrix TE=(t1×e1、t2×e2..., tm×em),
Take final ventilation scheme of the aeration-cooling condition corresponding to least member as tunnel in TE.
2. the control method of High-geotemperature railway tunnel cooling ventilation according to claim 1, it is characterised in that distinguished and admirable speed
Angle value is between 2 meter per seconds between 10 meter per seconds, and the distinguished and admirable duration is between 2 hours/day between 8 hours/day.
3. the control method of High-geotemperature railway tunnel cooling ventilation according to claim 1, it is characterised in that the control
Method is carried out using computer simulation, and before computer simulation, tunnel curve of ground temperature data, year outside hole are inputted into computer
Temperature Change curve data, m distinguished and admirable velocity amplitude values, m distinguished and admirable duration values；
The computer simulation is carried out in the following way：By the convective heat-transfer coefficient of the tunnel wall obtained on the spot and air,
The heat transfer coefficient input computer of country rock, using tunnel curve of ground temperature data as the original temperature value of tunnel wall, with year outside hole
Temperature Change curve data is as tunnel portal air themperature boundary value, with the convective heat-transfer coefficient of tunnel wall and air, each
The airflow velocity value at time point is as air and cold and hot convection heat transfer' heat-transfer by convection relative velocity, the tunnel surrounding at each time point of tunnel wall
The difference of temperature value and air temperature value is that thermograde obtains heat output, and the heat output in N is summed, obtains N
Interior force ventilation heat output, finally obtain tunnel wall surface temperature drop.
4. the control method of High-geotemperature railway tunnel cooling ventilation according to claim 3, it is characterised in that also include row
Car produces the heat transfer conversion of Piston Action Wind, i.e., tunnel Piston Action Wind Calculation of Heat Transfer, the tunnel Piston Action Wind Calculation of Heat Transfer are：By Piston Action Wind
Curve data inputs computer, and the tunnel wind speed angle value at each time point is relative with the cold and hot convection heat transfer' heat-transfer by convection of tunnel wall as air
The difference of speed, the tunnel wall temperature value at each time point and air temperature value is that thermograde obtains Piston Action Wind heat output, to N
Piston Action Wind heat output in year is summed, and obtains the tunnel Piston Action Wind heat output in N, tunnel Piston Action Wind heat output and machinery
Total heat transfers ventilation heat output and that value is as tunnel palisades, finally obtain the final temperature of tunnel wall, the most final temperature
Angle value is the reference value of a certain setting value.
5. the control method of High-geotemperature railway tunnel cooling ventilation according to claim 4, it is characterised in that tunnel ground temperature
Outside curve data, hole year temperature Change curve data and Piston Action Wind curve data acquisition pattern be following manner：
The tunnel curve of ground temperature data for survey the tunnel longitudinal direction ground temperature curve data that borehole data is drafted；
Temperature Change curve data was drawn and obtained by following formula year outside the hole：
<mrow>
<msub>
<mi>t</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<msub>
<mi>t</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mn>8760</mn>
</mfrac>
<mo>&CenterDot;</mo>
<mi>&tau;</mi>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mi>&pi;</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
In formula, taFor local year air-flow mean temperature；tmFor gas flow temperature amplitude, τ 1 is certain point on year temperature Change curve outside hole
In hourage residing then；
The Piston Action Wind curve data is drawn and obtained by following formula：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>v</mi>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>v</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>v</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
Wherein：V=vmWhen situation is that vehicle driving direction is consistent with blower ventilation direction；
V=-vmWhen situation is vehicle driving direction and blower ventilation in opposite direction；
When v=0 situations is in tunnels without vehicle.
6. the control method of High-geotemperature railway tunnel cooling ventilation according to claim 1, it is characterised in that met
The blower fan wind speed and ventilation time of tunnel cooling conditions are realized in the following way：By tunnel longitudinal direction Unsteady Temperature Field country rock portion
Governing equation (1), governing equation (2) the composition tunnel of the tunnel longitudinal direction distinguished and admirable part of Unsteady Temperature Field divided is longitudinal three-dimensional non-steady
The complete mathematical description of state heat transfer：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mi>a</mi>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mi>r</mi>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>&GreaterEqual;</mo>
<mi>R</mi>
<mo>,</mo>
<mi>&tau;</mi>
<mo>></mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&rho;Ac</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>&tau;</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>v</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>h</mi>
<mi>U</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>b</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>q</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>></mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
The complete asynchronous long display difference equation of three-dimensional unstable state heat transfer is established using finite difference calculus, it is as follows：
The temperature of country rock nodeApproximate expression formula：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>F</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>a</mi>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>F</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>a</mi>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>F</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mi>a</mi>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mi>a</mi>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
The temperature T of convective boundary node 1j,n+1Approximate expression formula：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>F</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>r</mi>
</mrow>
<mrow>
<mn>4</mn>
<mi>R</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<mfrac>
<mrow>
<mi>a</mi>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>8</mn>
<mi>h</mi>
<mi>R</mi>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mi>&rho;</mi>
<mi>&Delta;</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>R</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mi>i</mi>
<mi>n</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>F</mi>
<mn>0</mn>
</msub>
<mi>R</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mi>a</mi>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<mrow>
<mn>2</mn>
<mi>&Delta;</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mn>8</mn>
<mrow>
<mn>4</mn>
<mi>R</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>a</mi>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mn>8</mn>
<mi>h</mi>
<mi>R</mi>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mi>&rho;</mi>
<mi>&Delta;</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>R</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
The longitudinal distinguished and admirable approximate temperature in tunnelExpression formula：
<mrow>
<msubsup>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mi>i</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>h</mi>
<mi>U</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&rho;Ac</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>q</mi>
<mi>s</mi>
</msub>
<mrow>
<msub>
<mi>&rho;Ac</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mi>v</mi>
<mfrac>
<mrow>
<msubsup>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
</mrow>
<mrow>
<mn>2</mn>
<mi>&Delta;</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mi>i</mi>
<mi>n</mi>
</msubsup>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msubsup>
<msub>
<mi>hUT</mi>
<mi>f</mi>
</msub>
<mi>i</mi>
<mi>n</mi>
</msubsup>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&rho;Ac</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>h</mi>
<mi>U</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&rho;Ac</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
To make Equations stable, equation (3), the calculating time step Δ τ of equation (4) and tunnel longitudinal space step delta x
Equation (6) and equation (7) should be met respectively：
<mrow>
<mi>&Delta;</mi>
<mi>&tau;</mi>
<mo>&le;</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<mi>a</mi>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>&tau;</mi>
<mo>&le;</mo>
<mfrac>
<mrow>
<mn>4</mn>
<mi>R</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>r</mi>
</mrow>
<mrow>
<mfrac>
<mrow>
<mn>2</mn>
<mi>a</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>R</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>a</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>R</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>a</mi>
<mi>&Delta;</mi>
<mi>r</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>8</mn>
<mi>h</mi>
<mi>R</mi>
</mrow>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mi>&rho;</mi>
<mi>&Delta;</mi>
<mi>r</mi>
</mrow>
</mfrac>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
More than it is various in, τ be setting time span, rjFor the radial distance of j-th of radial point in tunnel, R is tunnel excavation half
Footpath, T are the temperature in country rock, and r is radial distance of any point along tunnel in country rock, x be in country rock any point along tunnel
Longitudinal direction distance, ρ be air-flow density, A be tunnel area of section, TfFor tunnel gas flow temperature, TbFor tunnel wall surface temperature, v
For airflow velocity, qsFor heat caused by equipment in tunnel, h is convection transfer rate, and U is tunnel cross section girth；(i, j) is
Along i-th of tunnel axis direction, radially j-th of grid node,For country rock node (i, j) the n moment temperature (DEG C)；Δr
For tunnel cross section radial distance step-length；Δ x is tunnel fore-and-aft distance step-length；Δ τ is time step；A is the heat conduction system of country rock
Number；FoFor Fourier number；cpThe specific heat at constant pressure of air-flow；For the wind-warm syndrome of i-node in n moment tunnels；
It is 1 year by a ventilation cycle of temperature prediction in tunnel, longitudinal ground temperature curve number value in tunnel is started as calculating
Moment tunnel surrounding longitudinal temperature initial value, year, temperature Change curve numerical value as a ventilation cycle inner tunnel entered one's intention as revealed in what one says outside hole
Flow temperature initial value, Piston Action Wind change curve influences airflow velocity as driving in daily tunnel, by daily ventilation time and
Corresponding airflow velocity, above difference equation is substituted into, the time step Δ τ and spatial mesh size determined by equation (6) and equation (7)
Δ x, difference equation (3), equation (4), equation (5) are solved, obtain different year tunnel surrounding and airflow temperature edge is vertical
To distribution, carry out tunnel in temperature prediction, obtain under the wind speed meet cooling require ventilation cycle t；
Outside tunnel curve of ground temperature data, hole year temperature Change curve data and Piston Action Wind curve data acquisition pattern be such as lower section
Formula：
The tunnel curve of ground temperature data for survey the tunnel longitudinal direction ground temperature curve data that borehole data is drafted；
Temperature Change curve data was drawn and obtained by following formula year outside the hole：
<mrow>
<msub>
<mi>t</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<msub>
<mi>t</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mn>8760</mn>
</mfrac>
<mo>&CenterDot;</mo>
<mi>&tau;</mi>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mi>&pi;</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
In formula, taFor local year air-flow mean temperature；tmFor gas flow temperature amplitude, τ 1 is certain point on year temperature Change curve outside hole
In hourage residing then；
The Piston Action Wind curve data is drawn and obtained by following formula：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>v</mi>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>v</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>v</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
Wherein v=vmWhen situation is that vehicle driving direction is consistent with blower ventilation direction；
V=-vmWhen situation is vehicle driving direction and blower ventilation in opposite direction；
When v=0 situations is in tunnels without vehicle.
7. the control method of High-geotemperature railway tunnel cooling ventilation according to claim 6, it is characterised in that to difference side
Journey (3), equation (4), equation (5) solve and carried out using software MATLAB.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610049329.6A CN105569709B (en) | 2016-01-25 | 2016-01-25 | The control method of High-geotemperature railway tunnel cooling ventilation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610049329.6A CN105569709B (en) | 2016-01-25 | 2016-01-25 | The control method of High-geotemperature railway tunnel cooling ventilation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105569709A CN105569709A (en) | 2016-05-11 |
| CN105569709B true CN105569709B (en) | 2018-01-23 |
Family
ID=55880265
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610049329.6A Active CN105569709B (en) | 2016-01-25 | 2016-01-25 | The control method of High-geotemperature railway tunnel cooling ventilation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105569709B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108979696B (en) * | 2018-06-30 | 2019-12-10 | 华南理工大学 | optimized exhaust method for subway tunnel heat extraction system |
| CN112131639B (en) * | 2020-09-11 | 2022-07-08 | 中南大学 | Numerical simulation method for high-speed train over-high ground temperature tunnel |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0633700A (en) * | 1992-07-13 | 1994-02-08 | Toshiba Corp | Tunnel ventilation control system |
| JP3929734B2 (en) * | 2001-10-03 | 2007-06-13 | 三菱電機株式会社 | Tunnel ventilation simulation device and tunnel ventilation control device |
| CN101235723B (en) * | 2008-02-02 | 2010-08-18 | 西南交通大学 | Express highway section multi- tunnel gathering type intelligent aeration control method |
| CN101509390A (en) * | 2009-02-25 | 2009-08-19 | 湖南大学 | Tunnel ventilation control method based on fuzzy control |
| CN103838197A (en) * | 2012-11-25 | 2014-06-04 | 西安众智惠泽光电科技有限公司 | Ventilation equipment coordination control system for tunnel construction |
| TWI515363B (en) * | 2013-12-09 | 2016-01-01 | 財團法人工業技術研究院 | Ventilation and cooling device for wind turbine nacelle and control method thereof |
| CN104612739A (en) * | 2014-11-28 | 2015-05-13 | 中铁二十三局集团有限公司 | Small-section tunnel single-head long-distance ventilation construction method |
-
2016
- 2016-01-25 CN CN201610049329.6A patent/CN105569709B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105569709A (en) | 2016-05-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lin et al. | A study on underground tunnel ventilation for piston effects influenced by draught relief shaft in subway system | |
| CN105569709B (en) | The control method of High-geotemperature railway tunnel cooling ventilation | |
| CN105569707A (en) | Highway tunnel ventilation feedforward control method based on environmental forecasting | |
| CN101429870A (en) | Fan setting method for road tunnel multi-silo system | |
| He et al. | Optimization of outdoor design temperature for summer ventilation for undersea road tunnel using field measurement and statistics | |
| CN109611139B (en) | Natural draft tunnel ventilation mechanism test system | |
| CN105604599A (en) | System and method for independently carrying out natural ventilation of tunnel through solar air pump | |
| Khaleghi et al. | Analysis of unsteady airflow in a subway station influenced by train movement | |
| Beusen et al. | Energy savings in light rail through the optimization of heating and ventilation | |
| CN110321587A (en) | Rail vehicle additional air resistance due to tunnel calculation method based on numerical simulation | |
| CN206161574U (en) | Test device of high ground temperature tunnel country rock and distinguished and admirable heat transfer | |
| Wei et al. | Study on effects of the train-induced airflow on the temperature field of high-speed railway tunnels in cold regions | |
| CN106050292B (en) | A kind of vcehicular tunnel chimney draft system based on Far-infrared Heating | |
| CN107060863A (en) | A kind of subway tunnel relief hole blast control method and system | |
| He et al. | Energy saving in metro ventilation system based on multi-factor analysis and air characteristics of piston vent | |
| Sambolek | Model testing of road tunnel ventilation in normal traffic conditions | |
| Wang et al. | Measurement-based investigation of subway station tunnel thermal environment | |
| CN208546183U (en) | The more gradient running tunnel fire test systems of subway | |
| CN112241596A (en) | Urban power grid medium-long term production simulation method considering electric automobile | |
| CN202673321U (en) | Novel ventilation experimental model for road tunnel | |
| CN105715290A (en) | Method for arranging fan in highway tunnel with ramp at entrance and exit | |
| CN108979696A (en) | A kind of subway tunnel heat-extraction system optimization exhaust method | |
| CN203452830U (en) | Tunnel natural ventilation cooling structure | |
| CN104452508A (en) | Large-scale frozen roadbed cooling structure capable of strengthening convective cooling efficiency of rock blocks | |
| CN104420879B (en) | A kind of tunnel forced ventilation cooling device and cool-down method |
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 |