严寒地区铁路隧道中心深埋水沟设置长度计算方法

严寒地区铁路隧道中心深埋水沟设置长度计算方法陈霄汉【期刊名称】《《铁道标准设计》》【年(卷)，期】2019(063)004【总页数】5页(P118-122)【关键词】铁路隧道;太阳辐射;统计分析;中心深埋水沟;非线性拟合【作者】陈霄汉【作者单位】陕西省铁道及地下交通工程重点实验室(中铁一院)西安710043【正文语种】中文【中图分类】U453.6引言根据我国已建成严寒地区隧道冻害调查概况，隧道冻害一般是从排水通道局部结冰开始的，完善有效的排水系统及可靠的防寒保温措施是防止冻害的关键［1-4］。目前，国内隧道防寒措施主要有保温水沟、中心深埋水沟、防寒泄水洞和电伴热丝系统等。隧道冻害的发生，往往与洞口气温、洞口微地形、洞口朝向、隧道纵坡、地下水发育程度、自然风向、风速、列车行驶速度等诸多因素有关。严寒地区隧道中心深埋水沟设置长度暂无完善的定量计算方法，设计时通常根据工程类比确定。随着已建成的严寒地区隧道数量不断增长，可能出现冻害现象的隧道数量将逐渐增多，设置合适长度的中心深埋水沟可以经济、合理、有效地解决隧道冻害问题。张国柱、夏才初［5-6］等人通过建立隧道传热模型，计算寒区隧道围岩轴向、径向温度的理论解，分析隧道洞内及围岩温度分布规律，确定寒区隧道防寒保暖段的设防长度。日本学者黑川羲范［7］回归分析264座隧道的实测数据，提出利用隧道洞口气温计算隧道内-5°C位置距洞口长度的经验公式。张生［8］分析我国建设部房建行业和日本国铁的研究成果后，提出适用于季节性冻土隧道洞口保温段长度、衬砌保温及水沟保温的设计计算方法。隧道设计中通常不考虑太阳辐射的大小，其对隧道的影响在运营后才会体现。对已建成位于严寒地区的铁路隧道排水通道局部结冰情况进行统计，位于山体阴坡的隧道洞口的排水通道更易结冰，威胁行车安全。根据运营后的现场实测，对隧道洞口中心深埋水沟的设置长度进行优化，结果见表1。表1优化后中心深埋水沟设置长度统计隧道坡向洞口最冷月平均气温/C中心深埋水沟设置长度/ml号洞口阳坡-13.351050.002号洞口阴坡-13.641400.003号洞口阳坡-17.521300.004号洞口阴坡-17.481650.005号洞口阳坡-17.441300.006号洞口阴坡-13.511450.007号洞口阳坡-13.481100.008号洞口阴坡-13.011450.009号洞口阳坡-12.961100.0010号洞口阴坡-12.361600.0011号洞口阳坡-12.361100.0012号洞口阴坡-12.321700.00如表1所示，最冷月平均气温大致相当的隧道洞口，位于山体阴坡的隧道洞口发生排水通道局部结冰的情况更为严重，需要设置更长的中心深埋水沟。山体的阴、阳坡是由于坡面受到的太阳辐射不同而形成，为了准确分析隧道洞口的太阳辐射情况，建立了隧道洞口太阳辐射评价方法，在黑川羲公式的基础上，引入太阳辐射指标分析隧道洞口的太阳辐射情况，结合最冷月平均气温和已建成隧道优化后的中心深埋水沟设置长度进行分析并提出经验公式。1太阳辐射强度影响因素太阳辐射是指太阳以电磁波的形式向夕卜传递能量，表示太阳辐射强弱的物理量，称为太阳辐射强度，太阳辐射强度的大小受到天气状况、地理纬度、海拔高低、照射角度等因素的影响。太阳辐射强度影响因素［9-12］如下。1.1太阳高度角太阳高度角是指观测点与地心连线的地平面和该点太阳光入射方向的夹角。太阳高度角为90°时，太阳辐射强度最大；太阳斜射地面角度越大，太阳辐射强度相对越小。太阳高度角计算公式如下sinh=sin^sin8+cos^cos8(1)式中h——太阳高度角；中——隧道洞口的地理纬度；6——赤纬角(太阳赤纬)，地球赤道平面与太阳和地球中心的连线之间的夹角。1.2太阳入射角太阳入射角是指太阳直射光线与壁面法线之间的夹角。太阳入射角由地表自然坡度、坡面方位角、太阳高度角和太阳方位角确定。cosi=cos0sinh+sin0coshcos(a-y)(2)式中i——太阳入射角；0——自然坡度；Y——坡面方位角，指坡面法线在水平面上的投影与正南方向(当地子午线)的夹角，—般以偏东为负，偏西为正；a——太阳方位角，指太阳至地面上某给定点的连线在水平面上的投影与正南方向的夹角。一般以偏东为负，偏西为正。1.3大气质量大气质量反映了大气对地球表面接收太阳光的影响程度，指太阳光线穿过地球大气层的路程与太阳在天顶位置时光线穿过地球大气层的路程之比。⑶式中m 大气质量。1.4水平面太阳辐射强度法向太阳辐射强度法向太阳辐射强度是指与太阳光线相垂直的表面上(即太阳光线法线方向)的太阳直射辐射强度IDN=I0xPm⑷式中P 大气透明系数；I0——太阳常数，太阳常数是指在日地平均距离(D=1.496x108km)上，大气顶界垂直于太阳光线的单位面积每秒钟接受的太阳辐射［13-15］。水平面直射太阳辐射强度水平面直射太阳辐射强度是指，未被地球大气层吸收、反射及折射仍保持原来的方向直达地球表面在水平面上的太阳辐射IDH=IDNsinh=I0Pmsinh(5)式中idh——水平面直射太阳辐射强度。⑶水平面散射辐射强度水平面散射辐射强度是指，经过大气和云层的反射、折射、散射作用，改变了原来的传播方向达到地球表面的、并无特定方向在水平面上的这部分太阳辐射(6)式中idH——水平面散射辐射强度。⑷水平面总辐射强度水平面总辐射强度指水平面散射辐射与直射辐射的总和⑺式中Ih——水平面总辐射强度。2隧道洞口太阳辐射评价2.1太阳辐射强度计算模型隧道洞口的山体有一定的自然坡度，因此，隧道洞口的太阳辐射强度按倾斜面进行计算。隧道洞口太阳直射辐射强度(8)式中ID0——隧道洞口太阳直射辐射强度。隧道洞口太阳散射辐射强度⑼式中Id0——隧道洞口太阳散射辐射强度。⑶隧道洞口地面反射辐射强度(10)式中ir0__隧道洞口地面反射辐射强度；pG——隧道洞口的地表反射率。常见典型地表面的地表反射率［16］详见表2。表2典型地表面的地表反射率地表状态反射率地表状态反射率沙漠0.24~028干草地0.15~0.25干燥裸地0.1~0.2湿草地0.14~0.26湿裸地0.08~0.09新雪0.81干湿土0.14冰面0.69⑷隧道洞口总辐射强度I0=ID0+Id0+IR0(11)式中I0——隧道洞口太阳总辐射强度。2.2隧道洞口太阳辐射指标为直观反映隧道洞口的太阳辐射强弱，现引入太阳辐射指标作为严寒地区隧道洞口太阳辐射评价指标。太阳辐射指标由隧道洞口总辐射强度计算公式简化而来，针对隧道洞口的太阳辐射特点，对部分固定的计算参数进行赋值。以北半球冬至日中午12时作为计算时间，确定赤纬角和太阳方位角，大气透明系数取我国西北高海拔地区12月的均值［17］,拟定隧道洞口的地表面为干燥裸地。太阳常数为固定值，为简化计算取值为1。太阳辐射指标计算结果的大小仅表示相同时间内、单位面积上的辐射量的大小，不表示实际接收到的太阳辐射值。太阳辐射指标固定参数见表3。表3太阳辐射指标固定参数固定参数取值赤纬角S23.45太阳方位角a0大气透明系数P(12月)0.722地表反射率pG0.2太阳常数I01.0太阳辐射指标计算(12)式中IQ——隧道洞口太阳辐射指标；固定部分参数的隧道洞口总辐射强度。经简化，严寒地区隧道洞口太阳辐射指标由隧道洞口的地理纬度、坡面方位角和山体自然坡度计算得到。选择表1中我国已建成位于严寒地区铁路隧道的洞口作为研究对象，太阳辐射指标计算参数见表4。表4太阳辐射指标计算参数隧道地理纬度/(°)坡面方位角/(°)自然坡度/(°)1号洞口-37.19-8.6039.712号洞口-37.33157.8813.163号洞口-37.73-23.5821.434号洞□-37.79173.2121.215号洞口-37.80-7.0419.856号洞口-37.88160.6419.267号洞口-37.88-20.7730.978号洞口-37.91155.3022.639号洞口-37.91-24.6928.2110号洞口-37.96172.6933.8911号洞口-37.96-5.7922.6612号洞口-38.03176.3440.78由公式(1)~公式(12)计算得到的太阳辐射指标结果如表5所示。表5太阳辐射指标计算结果统计隧道洞口IQ隧道洞口IQ1号0.567号0.502号0.238号0.153号0.469号0.494号0.1510号0.045号0.4611号0.466号0.1712号0.01经计算，12个隧道洞口的太阳辐射指标均值在山体阴坡为0.13，阳坡为0.49，与山体坡向一致，可以准确反映隧道洞口的太阳辐射情况。3中心深埋水沟设置长度计算3.1方程形式日本学者黑川羲范提出的黑川羲公式［7］如下所示L=65.6x(5-T)0.818(13)式中T——洞口最冷月平均气温，°C;L——防寒措施设置长度，m。在黑川羲公式形式的基础上引入太阳辐射指标，拟定公式形式如下(14)式中T——洞口最冷月平均气温，°C;L——中心深埋水沟设置长度，m；a、b、c、d——待回归的计算参数，其中d为收敛控制系数。3.2防寒措施设置长度计算公式计算参数非线性回归方法为麦夸特法(Levenberg-Marquardt)[18-21],基本原理是通过迭代程序来计算残差的平方和，当残差平方和达到最小值时，迭代过程结束，所得到的拟合方程就是曲线拟合的最终结果。研究对象为已建成的位于严寒地区铁路隧道，T值、L值分别取表1中洞口最冷月平均气温和优化后的中心深埋水沟设置长度，IQ值取表5的太阳辐射指标计算结果。将隧道洞口太阳辐射指标、洞口最冷月平均气温作为自变量，优化后的严寒地区中心深埋水沟设置长度作为因变量，进行非线性回归后得到结果如表6所示。表6防寒措施长度计算回归统计计算参数abcd回归值1074-7.6250.2240.926优化算法麦夸特法^▼£门6£丐*30"3汕)收敛判断标准10-10均方差(RMSE)12.147残差平方和(SSE)1770.544相关系数R0.999相关系数之平方(R2)0.997决定系数(DC)0.997卡方系数(Chi-Square)0.680F统计(F-Statistic)3710.332根据回归结果，严寒地区中心深埋水沟设置长度计算公式如下(15)式中L——严寒地区中心深埋水沟设置长度，m。严寒地区中心深埋水沟设置长度计算公式曲面如图1所示。图1中心深埋水沟设置长度曲面4结论与建议在黑川羲公式的基础上，考虑了隧道洞口的太阳辐射对中心深埋水沟设置长度的影响，基于已建成的位于严寒地区铁路隧道优化后的中心深埋水沟设置长度，结合太阳辐射指标、最冷月平均气温，采用麦夸特法非线性回归得到中心深埋水沟设置长度计算公式。由于公式回归过程中资料有限，今后需要在以下方面进一步研究。实际工程中，计算结果需要结合隧址区地下水发育程度、洞口风向、洞口风速、列车行驶速度等因素做出适当优化。非线性拟合过程建立在我国西北高海拔严寒地区铁路隧道优化后的防寒措施的基础上，公路隧道以及东北等其他严寒地区隧道的适用性需要进一步论证。⑶由于实测资料有限，公式适用于洞口最冷月平均气温在-7.55~-18.0°C的铁路隧道，最冷月平均气温低于-18.0C隧道洞口的适应性需要进一步验证。参考文献：【相关文献】夏勇，马志富.严寒地区隧道防排水改进措施探讨[J].铁道标准设计，2012(S1):27-29,73.刘利玮，何文社，路仕洋，等.当金山隧道防寒泄水洞适宜埋设位置数值分析[J].铁道标准设计，2018,62(4):130-135.宁少武.严寒地区客运专线隧道冻害预防技术和管理措施探讨[J].铁道标准设计，2012(S1):147-149.马志富.铁路隧道防排水设计探讨[J].铁道标准设计，2011(5):76-79.张国柱，夏才初，殷卓.寒区隧道轴向及径向温度分布理论解[J].同济大学学报(自然科学版)，2010,38(8):1117-1122,1160.夏才初，范东方，李强，等.寒区隧道保温层铺设长度的计算方法[J].同济大学学报(自然科学版)，2016,44(9):1363-1370.黑川羲范(Japan).卜^^儿^以以防止刀现状上冏题点[J].卜^初比地下，1980,11(12):13-21.张生.高寒隧道保温设计研究及热工计算方法初探[J].现代隧道技术，2011,48(3):59-63.张春桂，文明章.利用卫星资料估算福建晴空太阳辐射[J].自然资源学报，2014,29(9):1496-1507.和清华，谢云.我国太阳总辐射气候学计算方法研究[J].自然资源学报，2010,25(2):308-319.林媛