CN109492247A - 用于严寒地区隧道洞口深埋中心水沟长度的计算方法 - Google Patents
用于严寒地区隧道洞口深埋中心水沟长度的计算方法 Download PDFInfo
- Publication number
- CN109492247A CN109492247A CN201811108326.0A CN201811108326A CN109492247A CN 109492247 A CN109492247 A CN 109492247A CN 201811108326 A CN201811108326 A CN 201811108326A CN 109492247 A CN109492247 A CN 109492247A
- Authority
- CN
- China
- Prior art keywords
- tunnel
- length
- cold
- hole
- solar radiation
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Structural Engineering (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明涉及一种用于严寒地区隧道洞口深埋中心水沟长度的计算方法,分析位于寒冷及严寒地区隧道洞口的太阳辐射情况,计算隧道洞口的太阳辐射指标;将隧道洞口太阳辐射指标、洞口最冷月平均气温、已建成位于寒冷或严寒地区铁路隧道实测的排水通道结冰长度以及对应的保温水沟设置长度以及水沟结冰长度作为影响因素,拟定非线性曲面方程形式;利用多元非线性回归,计算出多元非线性曲面方程待定参数,得到寒冷及严寒地区隧道保温水沟设置长度计算公式。本发明克服了现有技术对保温水沟设置长度无法理论定性定量分析的缺点,提出了具有针对性的计算方法,提高了用于寒冷及严寒地区隧道洞口保温水沟设置长度的计算的准确性。
Description
技术领域
本发明涉及隧道防寒领域,具体涉及一种用于严寒地区隧道洞口深埋中心水沟长度的计算方法。
背景技术
根据我国已建成严寒地区隧道冻害调查概况,隧道冻害一般是从排水通道局部结冰开始的。对我国西北某新建高速铁路隧道进行统计,部分隧道建成后存在局部排水通道结冰的情况,严重威胁行车安全。特别是隧道洞口位于山体的阴坡且气温较低的区域,结冰问题更为突出,山体的阴、阳坡判断受到太阳辐射的直接影响。
完善有效的排水系统及可靠的防寒保温措施是防止冻害的关键。我国严寒地区隧道抗冻设防段的主要措施有防寒泄水洞、中心深埋水沟、中心深埋水沟和电伴热系统。隧道洞口中心深埋水沟长度设计时,通常依据设计人员的设计经验,缺乏定量的计算方法。
发明内容
本发明的目的是提供一种用于严寒地区隧道洞口深埋中心水沟长度的计算方法,综合考虑太阳辐射量和隧道洞口气温等影响因素,通过非线性回归确定严寒地区隧道中心深埋水沟的设置长度。
本发明所采用的技术方案为:
用于严寒地区隧道洞口深埋中心水沟长度的计算方法,其特征在于:
包括以下步骤:
步骤一:隧道洞口太阳辐射影响分析:
对位于寒冷及严寒地区隧道洞口的太阳辐射条件进行分析,通过隧道洞口所在坡面的方位角、坡面的自然坡度和地理纬度三个影响因素,计算得到隧道洞口的太阳辐射指标;
步骤二:拟定方程形式:
将隧道洞口太阳辐射指标、洞口最冷月平均气温作为自变量,已建成位于严寒地区铁路隧道实测的排水通道结冰长度以及对应的中心深埋水沟设置长度作为因变量,拟定多元非线性曲面方程形式;
步骤三:多元非线性回归:
计算得到多元非线性曲面方程待定的固定参数和待回归的收敛控制系数,对多元非线性曲面方程进行匹配,得到严寒地区隧道中心深埋水沟设置长度计算公式。
隧道洞口太阳辐射指标计算公式如下:
其中,Iθ为隧道洞口太阳辐射指标;IDθ为隧道洞口太阳直射辐射强度;Idθ为隧道洞口太阳散射辐射强度;IRθ为隧道洞口获得的地面反射辐射强度;IDN为法向太阳辐射强度,指太阳光线法线方向的太阳直射辐射强度,IDN=I0×Pm,I0为太阳常数,P为大气透明系数,m为大气质量,h为太阳高度角;i为太阳入射角, cosi=cosθsinh+sinθcoshcos(α-γ),γ为斜面方位角;α为太阳方位角;IDH为水平方向太阳直射辐射强度,IDH=IDNsinh=I0Pmsinh; IdH为水平方向太阳散射辐射强度θ为隧道洞口所在坡面的倾角;ρG为隧道所在坡面的地表反射率;Ih为水平方向的总辐射强度,
在计算太阳辐射指标时对部分变量进行赋值,其中:太阳方位角α为0,指时刻为12时的太阳方位角;P大气透明系数为0.722;隧道所在坡面的地表反射率ρG为0.2;太阳常数I0取1。
步骤二中,将隧道洞口太阳辐射指标、洞口最冷月平均气温作为自变量,已建成位于严寒地区铁路隧道实测的排水通道结冰长度以及对应的中心深埋水沟设置长度作为因变量,拟定多元非线性方程形式如下;
其中,T为洞口最冷月平均气温;Iθ为隧道洞口太阳辐射指标; L为高速铁路隧道的冻害长度;a、b、c为待回归的固定参数,d为待回归的收敛控制系数。
步骤三中,利用麦夸特法进行多元非线性回归,计算得到多元非线性曲面方程待定的固定参数和待回归的收敛控制系数,进而得到严寒地区隧道中心深埋水沟设置长度计算公式如下:
所述的严寒地区隧道中心深埋水沟设置长度计算公式的计算结果建立在隧道洞口最冷月平均气温、太阳辐射指标的基础上;严寒地区隧道中心深埋水沟设置长度还受到其它因素的影响,当引入洞口微地形、地下水发育程度、隧道坡度、自然风向、风速、列车行驶速度的其它影响因素时,需要在步骤二中增加上述自变量,重新拟定多元非线性方程形式,并在步骤三中增加待回归固定参数,重新进行多元非线性回归:
其中,m1m2...mn为影响严寒地区隧道中心深埋水沟设置长度计算公式的计算结果的其它因素,包括但不仅限于隧道洞口微地形、地下水发育程度、隧道坡度、自然风向、风速、列车行驶速度;为运算过程,表示各影响因素对计算结果的影响方式。
本发明具有以下优点:
本发明所述的一种用于严寒地区隧道洞口中心深埋水沟设置长度的计算方法,针对隧道洞口的特点,隧道洞口所在坡面的方位角、坡面的自然坡度和地理纬度三个影响因素计算隧道洞口的太阳辐射指标,对隧道洞口接收到的太阳辐射大小进行评价。分析隧道洞口太阳辐射指标、洞口最冷月平均气温对寒冷或严寒地区铁路隧道实测的排水通道结冰长度以及对应的保温水沟设置长度的影响,拟定了多元非线性曲面方程形式,通过多元非线性回归,提出了严寒地区隧道中心深埋水沟设置长度计算公式。本发明克服了现有技术对位于严寒地区隧道洞口中心深埋水沟设置长度无法理论定性定量分析的缺点,提出了具有针对性的计算方法,提高了严寒地区隧道洞口中心深埋水沟设置长度计算的准确性。
具体实施方式
下面结合具体实施方式对本发明进行详细的说明。
本发明涉及的一种用于严寒地区隧道洞口中心深埋水沟设置长度的计算方法,包括:
步骤一:隧道洞口太阳辐射影响分析。对位于寒冷及严寒地区隧道洞口的太阳辐射条件进行分析,通过隧道洞口所在坡面的方位角、坡面的自然坡度和地理纬度三个影响因素,计算得到隧道洞口的太阳辐射指标;
步骤二:拟定方程形式。将隧道洞口太阳辐射指标、洞口最冷月平均气温作为自变量,我国已建成位于严寒地区铁路隧道实测的排水通道结冰长度以及对应的中心深埋水沟设置长度作为因变量,拟定多元非线性曲面方程形式;
步骤三:多元非线性回归。计算得到多元非线性曲面方程待定的固定参数和待回归的收敛控制系数,对多元非线性曲面方程进行匹配,得到严寒地区隧道中心深埋水沟设置长度计算公式。
所述步骤一隧道洞口太阳辐射指标计算公式如下:
其中,Iθ为隧道洞口太阳辐射指标;IDθ为隧道洞口太阳直射辐射强度;Idθ为隧道洞口太阳散射辐射强度;IRθ为隧道洞口获得的地面反射辐射强度;IDN为法向太阳辐射强度,指太阳光线法线方向的太阳直射辐射强度,IDN=I0×Pm,I0为太阳常数,P为大气透明系数,m为大气质量,h为太阳高度角;i为太阳入射角,cosi=cosθsinh+sinθcoshcos(α-γ),γ为斜面方位角;α为太阳方位角;IDH为水平方向太阳直射辐射强度, IDH=IDNsinh=I0Pmsinh;IdH为水平方向太阳散射辐射强度θ为隧道洞口所在坡面的倾角;ρG为隧道所在坡面的地表反射率;Ih为水平方向的总辐射强度,
所述步骤一隧道洞口的太阳辐射受到太阳光线入射角、地形、云量、洞口微地形等因素的影响。我国直接测量太阳辐射的站点较少,很多地区并没有太阳辐射的观测值,需要间接计算获得,因而计算中需要选用形式简单,待定参数少的太阳辐射计算模型。为简化计算过程,本发明在计算太阳辐射指标时对部分变量进行赋值,其中:太阳方位角α为0,指时刻为12时的太阳方位角;P大气透明系数为0.722;隧道所在坡面的地表反射率ρG为0.2;太阳常数I0取1。
所述步骤二将隧道洞口太阳辐射指标、洞口最冷月平均气温作为自变量,我国已建成位于严寒地区铁路隧道实测的排水通道结冰长度以及对应的中心深埋水沟设置长度作为因变量,拟定多元非线性方程形式如下;
其中,T为洞口最冷月平均气温;Iθ为隧道洞口太阳辐射指标; L为高速铁路隧道的冻害长度;a、b、c为待回归的固定参数,d为待回归的收敛控制系数。
所述步骤三利用麦夸特法(Levenberg-Marquardt)进行多元非线性回归,计算得到多元非线性曲面方程待定的固定参数和待回归的收敛控制系数,进而得到严寒地区隧道中心深埋水沟设置长度计算公式如下:
所述的严寒地区隧道中心深埋水沟设置长度计算公式的计算结果建立在隧道洞口最冷月平均气温、太阳辐射指标的基础上。严寒地区隧道中心深埋水沟设置长度还受到其它因素的影响,当引入洞口微地形、地下水发育程度、隧道坡度、自然风向、风速、列车行驶速度等其它影响因素时,需要在步骤二中增加上述自变量,重新拟定多元非线性方程形式,并在步骤三中增加待回归固定参数,重新进行多元非线性回归:
其中,m1m2...mn为影响严寒地区隧道中心深埋水沟设置长度计算公式的计算结果的其它因素,包括但不仅限于隧道洞口微地形、地下水发育程度、隧道坡度、自然风向、风速、列车行驶速度。为运算过程,表示各影响因素对计算结果的影响方式。
所述严寒地区隧道中心深埋水沟设置长度计算公式的拟合过程建立在我国已建成位于严寒地区铁路隧道实测的排水通道结冰长度以及对应的中心深埋水沟设置长度的基础上,隧道位置处于我国西北地区高海拔温度较低的区域,我国东北等其它严寒地区隧道的适用性需要进一步论证;
受实测资料限制,严寒地区隧道中心深埋水沟设置长度计算公式仅适用于隧道洞口最冷月平均气温在-7.55℃至-18.0℃之间,且需要设置中心深埋水沟的隧道,最冷月平均气温低于-18.0℃隧道洞口的适应性需要进一步验证。
本发明所述的一种用于严寒地区隧道洞口中心深埋水沟设置长度的计算方法,针对隧道洞口的特点,隧道洞口所在坡面的方位角、坡面的自然坡度和地理纬度三个影响因素计算隧道洞口的太阳辐射指标,对隧道洞口接收到的太阳辐射大小进行评价。分析隧道洞口太阳辐射指标、洞口最冷月平均气温对寒冷或严寒地区铁路隧道实测的排水通道结冰长度以及对应的保温水沟设置长度的影响,拟定了多元非线性曲面方程形式,通过多元非线性回归,提出了严寒地区隧道中心深埋水沟设置长度计算公式。本发明克服了现有技术对位于严寒地区隧道洞口中心深埋水沟设置长度无法理论定性定量分析的缺点,提出了具有针对性的计算方法,提高了严寒地区隧道洞口中心深埋水沟设置长度计算的准确性。
本发明的内容不限于实施例所列举,本领域普通技术人员通过阅读本发明说明书而对本发明技术方案采取的任何等效的变换,均为本发明的权利要求所涵盖。
Claims (6)
1.用于严寒地区隧道洞口深埋中心水沟长度的计算方法,其特征在于:
包括以下步骤:
步骤一:隧道洞口太阳辐射影响分析:
对位于寒冷及严寒地区隧道洞口的太阳辐射条件进行分析,通过隧道洞口所在坡面的方位角、坡面的自然坡度和地理纬度三个影响因素,计算得到隧道洞口的太阳辐射指标;
步骤二:拟定方程形式:
将隧道洞口太阳辐射指标、洞口最冷月平均气温作为自变量,已建成位于严寒地区铁路隧道实测的排水通道结冰长度以及对应的中心深埋水沟设置长度作为因变量,拟定多元非线性曲面方程形式;
步骤三:多元非线性回归:
计算得到多元非线性曲面方程待定的固定参数和待回归的收敛控制系数,对多元非线性曲面方程进行匹配,得到严寒地区隧道中心深埋水沟设置长度计算公式。
2.根据权利要求1所述的用于严寒地区隧道洞口深埋中心水沟长度的计算方法,其特征在于:
隧道洞口太阳辐射指标计算公式如下:
其中,Iθ为隧道洞口太阳辐射指标;IDθ为隧道洞口太阳直射辐射强度;Idθ为隧道洞口太阳散射辐射强度;IRθ为隧道洞口获得的地面反射辐射强度;IDN为法向太阳辐射强度,指太阳光线法线方向的太阳直射辐射强度,IDN=I0×Pm,I0为太阳常数,P为大气透明系数,m为大气质量,h为太阳高度角;i为太阳入射角,cosi=cosθsinh+sinθcoshcos(α-γ),γ为斜面方位角;α为太阳方位角;IDH为水平方向太阳直射辐射强度,IDH=IDNsinh=I0Pmsinh;IdH为水平方向太阳散射辐射强度θ为隧道洞口所在坡面的倾角;ρG为隧道所在坡面的地表反射率;Ih为水平方向的总辐射强度,
3.根据权利要求2所述的用于严寒地区隧道洞口深埋中心水沟长度的计算方法,其特征在于:
在计算太阳辐射指标时对部分变量进行赋值,其中:太阳方位角α为0,指时刻为12时的太阳方位角;P大气透明系数为0.722;隧道所在坡面的地表反射率ρG为0.2;太阳常数I0取1。
4.根据权利要求3所述的用于严寒地区隧道洞口深埋中心水沟长度的计算方法,其特征在于:
步骤二中,将隧道洞口太阳辐射指标、洞口最冷月平均气温作为自变量,已建成位于严寒地区铁路隧道实测的排水通道结冰长度以及对应的中心深埋水沟设置长度作为因变量,拟定多元非线性方程形式如下;
其中,T为洞口最冷月平均气温;Iθ为隧道洞口太阳辐射指标;L为高速铁路隧道的冻害长度;a、b、c为待回归的固定参数,d为待回归的收敛控制系数。
5.根据权利要求4所述的用于严寒地区隧道洞口深埋中心水沟长度的计算方法,其特征在于:
步骤三中,利用麦夸特法进行多元非线性回归,计算得到多元非线性曲面方程待定的固定参数和待回归的收敛控制系数,进而得到严寒地区隧道中心深埋水沟设置长度计算公式如下:
6.根据权利要求5所述的用于严寒地区隧道洞口深埋中心水沟长度的计算方法,其特征在于:
所述的严寒地区隧道中心深埋水沟设置长度计算公式的计算结果建立在隧道洞口最冷月平均气温、太阳辐射指标的基础上;严寒地区隧道中心深埋水沟设置长度还受到其它因素的影响,当引入洞口微地形、地下水发育程度、隧道坡度、自然风向、风速、列车行驶速度的其它影响因素时,需要在步骤二中增加上述自变量,重新拟定多元非线性方程形式,并在步骤三中增加待回归固定参数,重新进行多元非线性回归:
其中,m1m2...mn为影响严寒地区隧道中心深埋水沟设置长度计算公式的计算结果的其它因素,包括但不仅限于隧道洞口微地形、地下水发育程度、隧道坡度、自然风向、风速、列车行驶速度;为运算过程,表示各影响因素对计算结果的影响方式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811108326.0A CN109492247B (zh) | 2018-09-21 | 2018-09-21 | 用于严寒地区隧道洞口深埋中心水沟长度的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811108326.0A CN109492247B (zh) | 2018-09-21 | 2018-09-21 | 用于严寒地区隧道洞口深埋中心水沟长度的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109492247A true CN109492247A (zh) | 2019-03-19 |
CN109492247B CN109492247B (zh) | 2022-12-27 |
Family
ID=65689369
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811108326.0A Active CN109492247B (zh) | 2018-09-21 | 2018-09-21 | 用于严寒地区隧道洞口深埋中心水沟长度的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109492247B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006234591A (ja) * | 2005-02-25 | 2006-09-07 | National Agriculture & Food Research Organization | 日射量演算方法及び装置 |
CN103761578A (zh) * | 2013-12-25 | 2014-04-30 | 国家电网公司 | 基于多元线性回归的太阳辐照预报方法 |
CN108316952A (zh) * | 2018-04-09 | 2018-07-24 | 铁科创恒新材料科技有限公司 | 高寒地区隧道保温防排水系统 |
-
2018
- 2018-09-21 CN CN201811108326.0A patent/CN109492247B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006234591A (ja) * | 2005-02-25 | 2006-09-07 | National Agriculture & Food Research Organization | 日射量演算方法及び装置 |
CN103761578A (zh) * | 2013-12-25 | 2014-04-30 | 国家电网公司 | 基于多元线性回归的太阳辐照预报方法 |
CN108316952A (zh) * | 2018-04-09 | 2018-07-24 | 铁科创恒新材料科技有限公司 | 高寒地区隧道保温防排水系统 |
Non-Patent Citations (1)
Title |
---|
朱国伟: "寒冷及严寒地区铁路隧道防排水设计探讨", 《铁道工程学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109492247B (zh) | 2022-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jun et al. | Field measurement of temperature inside tunnel in winter in Gangwon, Korea | |
CN106951612A (zh) | 土壤冻融过程中动态蓄水容量产流计算方法 | |
Yokobori et al. | Effect of land cover on air temperatures involved in the development of an intra-urban heat island | |
Zhao et al. | Field measurement of air temperature in a cold region tunnel in northeast China | |
Wu et al. | The thermal effect of differential solar exposure on embankments along the Qinghai–Tibet Railway | |
Lewis et al. | Surface energy balance and meltwater production for a Dry Valley glacier, Taylor Valley, Antarctica | |
Johanneck et al. | Comprehensive evaluation of effect of climate in mechanistic–empirical pavement design GuidePredictions | |
CN107761528A (zh) | 一种基于实时高温天气和交通量的沥青路面车辙预报方法 | |
CN109492248A (zh) | 用于寒冷及严寒地区隧道洞口保温水沟长度的计算方法 | |
CN115099162B (zh) | 一种复杂地形下风场的订正方法 | |
Tai et al. | Study on thermal performance of novel asymmetric crushed-rock-based embankment on the Qinghai-Tibet Railway in permafrost region | |
Ya-ling et al. | Calculation of temperature differences between the sunny slopes and the shady slopes along railways in permafrost regions on Qinghai–Tibet Plateau | |
Sokol et al. | First experience with the application of the METRo model in the Czech Republic | |
CN109492247A (zh) | 用于严寒地区隧道洞口深埋中心水沟长度的计算方法 | |
CN109214102A (zh) | 用于严寒地区隧道洞口防寒泄水洞长度的计算方法 | |
CN109268067A (zh) | 用于寒冷及严寒地区隧道防寒排水措施的判定方法 | |
Ižvolt et al. | Some experience and preliminary conclusions from the experimental monitoring of the temperature regime of subgrade structure | |
Lai et al. | Research situation and prospect for highway snowmelt deicing technology with electric heat tracing | |
Abild et al. | Application of the wind atlas method to extreme wind speed data | |
CN116401499A (zh) | 一种寒区运营长大隧道冬季洞内纵向气温场的计算方法 | |
CN110502847A (zh) | 基于热量收支平衡原理的多年冻土区热棒路基设计计算方法 | |
De Blasi | Scale dependence of hydrological effects from different climatic conditions on glacierized catchments | |
Saarelainen et al. | Climate change and freezing conditions in Finland | |
Naqvi et al. | Development, Construction, and Instrumentation of Pilot Freeze–Thaw Resistant Granular Roadways Test Cells | |
Ling et al. | Research on snow cover distance model for tunnel entrance/exit section in cold region under multi-parameter coupling conditions |
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 |