Processing math: 27%
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

山区铁路“路基-环境”的耦合优化调控

鲍学英 沈杜华 李亚娟 贺振霞 张程皓 陈慧鑫

鲍学英, 沈杜华, 李亚娟, 贺振霞, 张程皓, 陈慧鑫. 山区铁路“路基-环境”的耦合优化调控[J]. 西南交通大学学报. doi: 10.3969/j.issn.0258-2724.20230340
引用本文: 鲍学英, 沈杜华, 李亚娟, 贺振霞, 张程皓, 陈慧鑫. 山区铁路“路基-环境”的耦合优化调控[J]. 西南交通大学学报. doi: 10.3969/j.issn.0258-2724.20230340
BAO Xueying, SHEN Duhua, LI Yajuan, HE Zhenxia, ZHANG Chenghao, CHEN Huixin. Coupling Optimization and Regulation of Roadbed and Environment in Mountainous Railways[J]. Journal of Southwest Jiaotong University. doi: 10.3969/j.issn.0258-2724.20230340
Citation: BAO Xueying, SHEN Duhua, LI Yajuan, HE Zhenxia, ZHANG Chenghao, CHEN Huixin. Coupling Optimization and Regulation of Roadbed and Environment in Mountainous Railways[J]. Journal of Southwest Jiaotong University. doi: 10.3969/j.issn.0258-2724.20230340

山区铁路“路基-环境”的耦合优化调控

doi: 10.3969/j.issn.0258-2724.20230340
基金项目: 国家自然科学基金资助项目(51768034);中央引导地方科技发展资金项目(22ZY1QA005)
详细信息
    通讯作者:

    鲍学英(1974—),女,宁夏中卫人,教授,博士,从事绿色铁路及工程管理方面的研究;E-mail:813257032@qq.com

  • 中图分类号: U215

Coupling Optimization and Regulation of Roadbed and Environment in Mountainous Railways

  • 摘要:

    为实现山区铁路“路基-环境”高质量协调发展,提出一种“路基-环境”耦合优化调控方法. 首先,对于铁路“路基-环境”相容共生进行定义,构建绿色要素指标体系,并结合耦合魔方游戏模型明晰耦合调控框架;其次,利用耦合协调度模型、施压承载模型及各关键要素间的函数关系共同构建“路基-环境”调控优化目标函数及约束条件,将路基工程绿色关键要素作为主控变量,采取智能优化算法(CSA)进行求解,得到相容共生状态下的各主控变量最优解;最后,以某山区铁路路基为例进行实证分析. 结果表明:当各主控变量路堑开挖尺寸、路堤填筑尺寸、支挡结构设计、支挡结构布设、工程防护结构设计、边坡工程防护布设、植物防护结构设计、边坡植物防护布设、生态声屏障结构设计、生态声屏障布设依次分别优化36.83%,43.14%,49.93%,68.91%,69.98%,68.91%,23.42%,68.91%,19.64%,19.60%比例时,可以实现铁路“路基-环境”从初级协调状态向中级协调状态的演化,研究结果验证了构建的“路基-环境”调控优化模型的合理性以及CSA算法求最优解的有效性,为实现山区铁路路基工程绿色化建设提供了科学参考.

     

  • 山区地势地形、气候极其复杂,环境敏感区众多,铁路路基工程的建设活动会给周边区域环境造成严重损害,同时,区域环境的脆弱性也会反馈于路基工程的建设活动并对其产生约束,彼此间相互作用、相互影响、相互制约[1],通俗表现为一方利益受损,另一方利益受益的状态,但在双方博弈过程中,会存在协调平衡点,此时,“路基-环境”会产生共同创造、相互促进的关系,促使整体系统利益最大化. 因此,如何探析路基工程与环境之间的耦合状况,并通过优化调控方式促使“路基-环境”向协调平衡状态进行演变,对铁路路基绿色高效建设具有重要意义.

    “路基-环境”绿色高效内涵不仅包含环境效益最大化原则,还容纳路基工程设计资源投入最小化的目的. 目前,大多数学者对于“路基-环境”绿色高效的研究较为局限于单一方面,甚至其中的一小部分,缺乏系统性思维. 在环境效益最大化原则,Su等[2]通过对不同根系与土壤之间的关系进行力学试验,为生态护坡技术的应用提供了参考;Sukmak等[3]研究石榴石残渣和轮胎衍生骨料改善不良地质的可行性,减少废弃材料对环境的损害;郑雨茜等[4]以铁路路基爆破为研究对象,从绿色施工理念的角度提出动态“三位一体”综合爆破效益评估模型. 在设计资源投入最小化目的上,李雨浓等[5]基于G1法和改进数据包络分析(DEA)模型对铁路工程材料投入进行评价,以此为依据减少材料的投入成本;鲍学英等[6]构建技术接口的关键要素均衡优化模型,分析质量、进度、成本、安全要素间的关系,并对关键要素进行均衡优化以减少管理成本. 因此,从系统整体性思维分析“路基-环境”协调发展,寻找铁路路基与环境之间和谐共存的利益平衡点具有一定的前景.

    “路基-环境”绿色高效的实施必然离不开工程智能优化决策,学者进行了广泛研究. 李伟等[7]建立一个综合考虑建设和运营费用、环境影响代价及线路约束的线形优化模型,并基于改进的遗传算法生成最优选线方案;段晓晨等[8]将三维非线性智能控制技术运用实际工程,通过建立张家庄大桥运营维护成本预测模型获取最优成本值,为桥梁运营阶段的管理提供借鉴和参考;薛新功等[9]对20世纪60年代至今国内外铁路线路专业智能优化策略进行总结,详细介绍各类方法的适用性及优缺点. 上述研究集中在桥梁、选线等专业领域,缺乏对铁路路基与环境之间的相互影响的分析决策.

    综上,本研究以山区铁路路基为研究对象,对“路基-环境”相容共生进行定义,构建绿色要素体系,利用耦合魔方游戏模型明确“路基-环境”调控优化的途径与方法,识别绿色关键要素作为主控变量,构建以“路基-环境”耦合协调度最优为目标、各绿色要素函数关系约束、“满载”条件约束下的耦合优化模型,基于MATLAB,采取变色龙群算法进行优化求解,得到最优状态下主控变量的最优值,通过模型最优值与工程实际值进行对比管理,促使路基工程绿色关键要素不断向最优值靠近,最终实现铁路“路基-环境”的相容共生.

    根据肖建庄等[10]对工程结构可持续性设计理论的架构,将“路基-环境”相容共生表征为路基工程建设活动对环境产生的影响维持在环境承载力所界定的合理阈值范围内,即达到一种“满载而不过载”的微妙平衡状态下,彼此间协调状态达到最优. 因此,为达到“路基-环境”相容共生的目的,本文依托铁路“路基-环境”交互影响特点,挖掘影响相容共生的绿色要素,并在此基础上提出耦合魔方游戏模型建立目标调控框架.

    山区铁路路基工程建设活动与环境之间存在错综复杂的联系,彼此间相互作用、相互影响、相互制约,会形成大量由一对一、一对多、多对多影响路径所构成的交互影响网络,网络中任意节点的变化都会牵一发而动全身,对铁路“路基-环境”相容共生产生影响,因此,参考文献[11-13]及《铁路路基设计规范》(TB 10001—2016)[14]、《铁路工程环境保护设计规范》(TB 10501—2016)[15],从微观层面“路基-环境”交互影响关系进行分析,如图1所示. 图中,黄色、紫色和黑色文字分别表示路基工程绿色要素、环境绿色要素和“路基-环境”影响路径.

    图  1  铁路“路基-环境”交互影响关系
    Figure  1.  Interaction between railway roadbed and environment

    在此基础上,对“路基-环境”交互影响网络内影响因素进行集成,从基体结构、地基处理、支挡结构、边坡防护、声屏障、防排水等设计环节提取路基工程绿色要素,分析其对地质环境、大气环境、水环境、生态环境、土壤环境、声环境的影响,识别环境绿色要素,构建铁路“路基-环境”绿色要素体系,如图2所示. 图中:z1为路基表土清理时单位体积产生的粉尘平均浓度;z11为从隧道工程弃渣场运送弃渣填筑可利用率;z12为路基挖方用作填方的弃渣利用率;z14z16z18分别为设计单位所设计的路基支挡结构、边坡工程防护、边坡植物防护的结构形式是否与周围环境相协调且功能有效,设计合理性能否达到要求,例如,在有滑坡体的路基地段,宜采用抗滑桩等支挡结构进行防护;z15z17z19为在实际施工中施工单位是否严格遵循技术要求与组织管理,从而确保施工质量达到设计标准,例如,边坡植物防护中,严格遵循植物种类选配标准,极大提高边坡稳定性、增强水土保持能力;z20表示安装声屏障前后在某特定位置上的声压级之差;z22z23z24z25表示路基截排水结构在单位时间内所能通过的最大截排水量;z27z28为路基施工区大气环境PM10(CO)浓度值;z30为在爆破影响下介质质点往复运动的速度;z31为单位面积土壤在单位时间内侵蚀量的大小;z32为通过措施减少水土流失面积与未实施措施时水土流失面积之比.

    图  2  铁路“路基-环境”绿色要素指标体系
    Figure  2.  Index system for green elements of railway roadbed and environment

    基于上述相容共生定义及绿色要素体系,提出铁路“路基-环境”耦合魔方游戏模型[16],建立铁路“路基-环境”调控优化目标框架,如图3所示. 将“路基-环境”的整个耦合过程映射为魔方的旋转操作过程,在一定时间内,将混乱状态的魔方还原成有序状态(混乱状态表现为魔方的各个小立方体颜色错乱,有序状态表现为魔方每一面颜色相同). 其中,每个小立方体代表一个相对完整的“路基-环境”耦合子系统,关节轴起到连接各个小立方体纽带作用,确保它们能够在旋转操作中相互作用、相互影响;路基工程的建设活动为转动魔方的动力,支撑魔方不断转动;属于水、土壤、声、大气、生态、地质等环境的绿色要素代表小立方体的每一面,表现为如何达到魔方有序状态的标识;路基与环境耦合形成的一面有序状态则表示单一环境与工程达到高度协调水平;由于本文提出的是一种“有损魔方”的理想模型,受制于魔方几何约束、物理约束的结构特点,会随着魔方的旋转产生使用损坏,当旋转次数达到魔方可以承受的限值区间,魔方就会失去使用价值,同理,环境的承载力表现为魔方可以承受转动次数区间,转动次数超出了魔方承受范围,认为魔方损坏.

    图  3  魔方游戏模型
    Figure  3.  Rubik’s cube game model

    魔方游戏的最终得分由三方面组成:还原度、错误率以及转动效率. 还原度指魔方的有序状态,即路基与环境之间的高度耦合协调水平;错误率是转动魔方过程中失误转动次数与魔方可承受次数的比值,表现为在旋转作用下魔方受损程度,即路基工程建设活动对环境的影响效果;转动效率是魔方从无序到有序所花费的步骤,即路基工程投入转化环境效益的能力. 游戏结果一般会出现3种状况:1)玩家在规定的时间内,还原度最低,错误率高,转动效率低;2)玩家超出规定时间的限制,还原度最高,错误率最低,转动效率高;3)玩家在规定的时间,通过利用关键小方块改进方法,还原度提升,错误率较低,转动效率较高. 因此,第3种结果得分最优,结合路基工程实践,应在满足铁路“路基-环境”满载的条件下,通过调控绿色关键要素达到协调水平最优,共同实现路基工程与环境之间的相容共生.

    基于耦合魔方游戏模型的构建过程,有必要识别绿色关键要素,将其作为耦合优化模型的主控变量,并以此采取施压承载模型对“路基-环境”的耦合效率进行测算,保证铁路“路基-环境”处于满载状态,结合耦合协调度模型对协调水平进行测度,促使“路基-环境”耦合协调状态有向高度协调状态演变的趋势.

    为消除量纲,使数据具有可比性,对要素进行极差归一化处理.

    xi,g={ximinximaxximinxi,正向要素,maxxiximaxximinxi,负向要素, (1)

    式中:xi为第i项绿色要素的实际值,xi,g为第i项绿色要素的归一化值.

    在“路基-环境”交互影响关系网络中,路基绿色要素对环境绿色要素的影响路径总体呈现一对一、一对多甚至多对多的情况,具有复杂化、虚拟性的特点,难以识别关键要素. 为保障后续优化调控顺利实施,借助鲍学英等[17]构建的复杂网络模型理论,将影响路径分为路基、环境系统内部影响和“路基-环境”系统间影响,通过二维云模型、Spearman相关系数对绿色关键要素进行识别.

    2.2.1   自身属性重要度

    Spearman相关系数是一种利用单调方程来评价2个变量之间相关性的方法,可以在同一系统内有效得到绿色要素的自身重要度,如式(2)所示:

    Wi=Mj=1ρij,ρij=16Kk=1d2KK(K21), (2)

    式中:Wi为路基、环境由于自身内部影响下第i项绿色要素的属性重要度,M为第i项绿色要素与其他绿色要素有关联关系的要素个数,ρij为第i项绿色要素与第j项绿色要素之间的Spearman相关系数,dk为在第k位专家打分中2个变量中的等级差,K指专家个数.

    2.2.2   影响间重要度

    二维云模型是一种可以将定性概念与定量数据进行转换的方法,通过表征路基建设活动对环境的影响后果Eh和概率Eg可以直观看出“路基-环境”之间的影响关系,确定绿色要素影响间重要度. 具体计算如式

    {Ex=1KKk=1uk,En=π2×1KKk=1|ukEx|,He=|S2E2n|,S2=1K1Kk=1(ukEx)2, (3)
    Li=11+(EhjEh)2+(EgjEg)2,Vi={CmaxiLi若靠近上界,Cmixi+Li若靠近下界, (4)

    式中:ExEnHeS2分别为样本期望、熵、超熵和样本方差;uk为第k位专家的评估分值;Li为第i项绿色要素二维云贴进度;Cimax C_i^{{\mathrm{mix}}} 分别为 i 项绿色要素贴进度所处区间范围的上限和下限; {V_i} 为第 i 项绿色要素二维云评价值; {E_{{\text{hj}}}} {E_{{\text{gj}}}} 分别为标准风险云后果等级和概率等级的期望值,具体云评价等级标准依照参考文献[17-18].

    2.2.3   综合重要度

    综合重要度由自身属性重要度及影响间重要度构成,综合重要度越大,代表该绿色要素越重要. 第 i 项绿色要素的综合重要度为

    {I_i} = \left( {{W_i}/\sum\limits_{i = 1}^{f} {{W_i}} + {V_i}/\sum\limits_{i = 1}^{f} {{V_i}} } \right)\Biggr/2,f=m,n, (5)

    式中:{V_i} 为第i 项绿色要素二维云评价值; \displaystyle\sum\limits_{i = 1}^{f} {{W_i}} \displaystyle\sum\limits_{i = 1}^{f} {{V_i}} 分别为自身属性重要度总值和影响路径重要度的总值; m 为路基工程系统的绿色要素个数,如图2中绿色要素 {x_1} {x_{19}},m=19 n 为环境系统的绿色要素个数,如图2中绿色要素 {x_{20}} {x_{26}},n=7 .

    实现铁路“路基-环境”相容共生的首要目的是他们彼此间的协调状态最优,而耦合协调度模型可在多个系统或子系统之间进行协调和交互的过程中分析事物的协调发展水平. 因此,通过耦合协调度模型对协调状态的测度,可以明晰“路基-环境”耦合协调状态是否有进一步提升空间[19-20],关系如式(6)所示

    \left\{ \begin{array}{l} \min D = 1 - \left(\dfrac{{{T_{\text{a}}} + {T_{\text{b}}}}}{2} \left(\dfrac{{{T_{\text{a}}}}}{{{T_{\text{b}}}}}\times [1 - \left( {{T_{\text{b}}} - {T_{\text{a}}}} \right)]\right)^{1/2}\right)^{1/2} ,\\ {T_{\text{a}}} = \displaystyle\sum\limits_{i = 1}^m { {{\omega _i}x_{i,{\text{g}}}^{}} } , \\ {T_{\text{b}}} = \displaystyle\sum\limits_{i = 1}^n { {{\omega _i}x_{i,{\text{g}}}^{}} } , \end{array} \right. (6)

    式中: {T_{\text{a}}} {T_{\text{b}}} 分别为路基工程、环境的综合发展水平,表示在绿色要素相互作用的影响下系统发展水平的高低; D 为耦合协调度,具体标准划分如表1所示[21-22],表示为路基工程与环境之间在绿色要素作用下耦合协调状态的强弱; \omega _{_i}^{} 为第 i 项绿色要素的权值.

    表  1  耦合协调度划分标准
    Table  1.  Classification criteria for ccoupling coordination degree
    第1分类划分标准第2分类及划分标准
    协调发展类0.9<D≤1.0优质协调
    0.8<D≤0.9良好协调
    0.7<D≤0.8中级协调
    0.6<D≤0.7初级协调
    过渡发展类0.5<D≤0.6勉强发展
    0.4<D≤0.5濒临失调
    失调衰退类0.3<D≤0.4轻度失调
    0.2<D≤0.3中度失调
    0.1<D≤0.2严重失调
    0<D≤0.1极度失调
    下载: 导出CSV 
    | 显示表格
    2.4.1   “路基-环境”影响关系约束

    参考《铁路路基设计规范》(TB 10001—2016)[14]、《铁路工程环境保护设计规范》(TB 10501—2016)[15]、《铁路声屏障工程设计规范》(TB 10505—2019)[23]、《土壤侵蚀分类分级标准》(sl190—2007)[24]、《环境空气质量标准》(GB 3095—2012)[25]、《声环境质量标准》(GB 3096—2008)[26]、《地表水环境质量标准》(GB 3838—2002)[27]、《爆破安全规程》(GB 6722—2014)[28]及文献[29],建立铁路“路基-环境”绿色要素函数关系,并明确环境绿色要素的限制范围,具体关系如式(7)所示.

    \left\{ \begin{gathered} {{\textit{z}}_{27}} = {{\textit{z}}_1} + y\left(1 + 2.5\left(\frac{{{{\textit{z}}_{10}}}}{b} - 1\right)\right), \\ {{\textit{z}}_{28}} = \frac{{({t_{\text{s}}}({c_1} - {{\textit{z}}_3}){e_{\text{s}}} + {t_{\text{y}}}{{\textit{z}}_{13}}{e_{\text{y}}}({c_2} - {{\textit{z}}_2}))\varphi }}{{{f_{\text{c}}}{\text{π}} {r^2}}}, \\ {{\textit{z}}_{29}} = ({s_{\text{c}}}{W_{\text{c}}}/100 + {s_{\text{w}}}{W_{\text{w}}}/100) (1 - {{\textit{z}}_{26}}), \\ {{\textit{z}}_{30}} = q{\left(\frac{{{\textit{z}}_{10}^{1/3}}}{R}\right)^\alpha }, \\ {{\textit{z}}_{31}} = 0.067{W^{1.627}_{\text{w}}}{L_{\text{S}}}C(0.6508 - 0.3436\times \\ \quad \lg({{\textit{z}}_{18}}{{\textit{z}}_{19}}/A)), \\ {{\textit{z}}_{32}} = ({{\textit{z}}_{16}}{{\textit{z}}_{17}} + {{\textit{z}}_{14}}{{\textit{z}}_{15}} + {{\textit{z}}_{17}}{{\textit{z}}_{18}})/A, \\ {{\textit{z}}_{33}} = ({{\textit{z}}_2} + {{\textit{z}}_3}){n_{\text{d}}} + \sum {{n_{\text{q}}}} + 20\lg \left(\frac{{0.006{\textit{z}}_{10}^{0.52}}}{{2 \times {{10}^{ - 5}}}}\right), \\ {{\textit{z}}_{34}} = {n_{\text{y}}} - {{\textit{z}}_{20}}, \\ {{\textit{z}}_{35}}_{} = \frac{{{\textit{z}}_{18}^{}{{\textit{z}}_{19}} + {{\textit{z}}_{21}}}}{{{{\textit{z}}_4}{{\textit{z}}_5} + {p_{\text{q}}}{{\textit{z}}_4} + {{\textit{z}}_7}{{\textit{z}}_8} + {p_{\text{t}}}{{\textit{z}}_7}}}, \\ {{\textit{z}}_{25}} \leqslant 50,{{\textit{z}}_{28}} \leqslant 4,{{\textit{z}}_{29}} \leqslant 20,8 \leqslant {{\textit{z}}_{30}} \leqslant 12,{{\textit{z}}_{31}} \leqslant 500, \\ {{\textit{z}}_{32}} \geqslant 85\% ,{{\textit{z}}_{33}} \leqslant 70,{{\textit{z}}_{34}} \leqslant 70,{{\textit{z}}_{35}} \geqslant 25\%, \end{gathered} \right. (7)

    式中: y 为原来的爆破产尘浓度(mg/m3), b 为原来爆破的炸药量(t), {t_{\text{s}}} 为施工机械工作台班数(台班), {t_{\text{y}}} 为运输机械运渣次数(次数), c_{1} c_{2} 分别为施工机械、运输机械总数目(台), {e_{\text{s}}} 为整体施工机械燃油消耗量均值(L), {e_{\text{y}}} 为整体运输机械燃油消耗量均值(L), \varphi 为燃油CO排放系数(mg/L), {f_{\text{c}}} 为路基长度与监测范围的比值, r 为CO监测点半径(m), {W_{\text{c}}} 为收集处理后废水(L), {W_{\text{w}}} 为未收集处理径流(L),s_{\mathrm{c}} {s_{\text{w}}} 分别为收集、未收集废水里固体悬浮物(SS)含量(mg/L), q \alpha 分别为爆破点至保护对象间的地形、地质条件有关的系数和衰减系数, R 为爆破后安全允许距离(m), {L_{\text{S}}} 为地形因子, C 为土壤可蚀性因子, A 为路基边坡面积(m2); {n_{\text{d}}} 、{n_{\text{q}}}、 {n_{\text{y}}}分别为电能机械设备、燃油机械设备、机车运行噪声平均值(dB), {p_{\text{q}}} 为路堑边坡坡面长(m), {p_{\text{t}}} 为路堤边坡坡面长(m).

    2.4.2   “满载”约束

    在魔方游戏模型中,魔方作为一个三维物体,具备空间特征,因此,将状态空间模型与魔方相结合[30],构成三维空间约束模型,为展现旋转过程中错误率对魔方的施压承载效果,引入超曲面概念,将三维空间映射成XOY二维超曲面状态,以X轴表示魔方横向约束,Y轴表示为魔方纵向约束,Z轴表示为魔方游戏的转动错误率,如图4所示. XYZ表现为最佳的错误率曲面. 点Ac表现为玩家在规定时间前提下错误率最佳状态,点Bc表现为玩家为降低错误率损失了时间,点Cc表现为玩家过于强调时间,错误率高于魔方承受的次数限值,魔方损坏,得分扣除.

    图  4  魔方游戏施压承载模型
    Figure  4.  Rubik’s cube game pressure bearing model

    因此,将魔方游戏施压承载模型应用于铁路“路基-环境”,确保路基工程建设活动对环境的影响始终处于XYZ曲面,使两系统在相互耦合过程中达到一种临界状态,即在不破坏环境条件下,路基工程设计资源投入最大化,构建“路基-环境”满载约束. 关系如式(8)所示.

    \left\{ \begin{array}{l} 1 - c \leqslant {X_{{\text{rgs}}}} = \dfrac{{{F_{l}} \cos \;\theta }}{{{F_{s}}}} \leqslant 1 + c , \\ {F_u} = { \sqrt{\displaystyle\sum\limits_{i = 1}^f {{w_i}x^{2}_{i,{\mathrm{g}}}} }},f = m,n,\;u \in \{{l,{\mathrm{s}}}\}, \\ \cos \;\theta = \dfrac{{\displaystyle\sum\nolimits_{i = 1}^{m} \displaystyle\sum\nolimits_{i = 1}^{n} {x^2_{i,{\text{g}}}} }}{\sqrt{{{ {\displaystyle\sum\limits_{i = 1}^m {x{{_{i,{\text{g}}}^2}^{}}\displaystyle\sum\limits_{i = 1}^n {x_{i,{\text{g}}}^2} } } }}}}, \end{array} \right. (8)

    式中: {F_{l}} 为环境理想承载力,表示为路基工程建设活动对环境的影响限值; {F_{s}} 为环境实际影响值,表示为路基工程建设活动对环境的实际影响效果; \theta 为随机旋转角度; {X_{{\text{rgs}}}} 为绿色耦合效率指数,表征路基工程与环境之间的施压承载状态,用于判断是否存在设计资源的浪费; c 为路基工程建设方的决策偏向,值越大,表示越偏向于环境保护,越小越,表示越偏向于工程设计投入高效利用,本文取值0.1; {X_{{\text{rgs}}}} \leqslant 1 - c ,表示路基工程建设活动强度小于环境承载力(可载); 1 - c \leqslant {X_{{\text{rgs}}}} \leqslant 1 + c ,表示路基工程建设活动强度刚好处于环境承载力之内(满载); 1 + c \leqslant {X_{{\text{rgs}}}} ,表示路基工程建设活动强度大于环境承载力(超载).

    变色龙群优化算法(CSA)是Braik于2021年提出的一种改进变色龙算法(CA). 该算法能够通过模拟变色龙猎食动态行为进行数学建模和计算,相较于传统的模拟退火算法(SA)、粒子群算法等,具有搜索能力强,收敛速度快等特征[31-32],且能够有效处理单、多目标复杂问题,符合本文所构建的优化调控模型. 主要步骤如下:

    步骤1 搜索猎物

    变色龙的搜寻食物行为如式(9)所示.

    {y}_{t + 1,\lambda ,\eta }^{}=\left\{\begin{array}{ll}{y}_{t,\lambda ,\eta }^{} + {p}_{1}\left({P}_{t,\lambda ,\eta }^{}-{G}_{t,\eta }^{}\right){r}_{2} +\\ \quad{p}_{2}\left({G}_{t,\eta }^{}-{y}_{t,\lambda ,\eta }^{}\right){r}_{1},\quad{\text{当}}{r}_{q}\geqslant{P}_{\text{p}}, \\ {y}_{t,\lambda ,\eta }^{} + \mu \left(\left({u}_{\text{b}}-{l}_{\text{b}}\right){r}_{3} + {l}_{\text{b}}\right)\\ \quad \mathrm{sgn}(r_{\text{a}}-0.5),\quad{\text{当}}{r}_{q} \lt {P}_{\text{p}}, \end{array} \right. (9)

    式中: y_{t,\lambda ,\eta }^{} t 次迭代第 \lambda 只变色龙在 \eta 维迭代后的位置;{p}_{1} {p}_{2} 为控制算法搜索能力的正系数; {P_{\text{p}}} 为感知猎物的概率;{r}_{q} {r}_{{\mathrm{a}}} 为(0,1)间的随机数; {{\mathrm{sgn}}} (r_{\text{a}} - 0.5) 代表变色龙旋转的方向,一般取值为±1; P_{t,\lambda ,\eta }^{} 为第 t 次迭代第 \lambda 只变色龙在 \eta 维迭代后的最佳位置; G_{t,\eta }^{} 为第 t 次迭代在 \eta 维中的全局最优位置; {r_1} , {r_2} , {r_3} 为0到1之间产生的随机数值; {u_{\text{b}}} {l_{\text{b}}} 分别为搜索空间的上界和下界; \mu 为通过 t 表达的变色龙搜索能力, \mu = {{\mathrm{e}}^{{{( - \delta t/{T_{\text{c}}})}^3}}} ,其中, \delta 为敏感系数, {T_{\text{c}}} 为最大迭代次数.

    步骤2 旋转眼睛搜索

    变色龙能利用眼睛360°旋转获取食物的具体位置,并且也会根据猎物的位置变化而更新自己的位置,更新位置关系如式(10)所示.

    \left\{\begin{array}{l} {y}_{t + 1,\lambda }^{}={m}_{\text{b}}({y}_{t,\lambda }-{y}_{t,{\lambda },{\text{q}}}) + {y}_{t,{\lambda },{\text{q}}},\\ {m}_{\text{b}}={R}_{\text{b}}(\theta ,{V}_{\text{z}1},{V}_{\text{z2}})\\ \theta ={r}_{\text{q}}\mathrm{sgn}(r_{\text{a}}-0.5)\text{π},\end{array}\right. (10)

    式中: {y}_{t + 1,\lambda }^{} 为第 \lambda 只变色龙在第 t + 1 次迭代位置; y_{t,{\lambda ,{\text{q}}}} 为第 \lambda 只变色龙第 t 次迭代的中心位置; {{\boldsymbol{m}}_{\text{b}}} 为变色龙眼睛的旋转矩阵; {{\boldsymbol{R}}_{\text{b}}} (·)为旋转矩阵的生成函数; {{\boldsymbol{V}}_{{\text{z}}1}} {{\boldsymbol{V}}_{{\text{z2}}}} 为坐标空间形成的2个正交向量.

    步骤3 吐舌捕食猎物

    在更新位置后,一旦猎物离自己较近,可以用吐舌的方式进行捕捉猎物,并且这时候的位置称为最优解,具体关系如式(11)所示.

    \left\{ \begin{array}{l} v_{t + 1,\lambda ,\eta }^{} = v_{t,\lambda ,\eta }^{} + {\varpi _1}{r_1}\left( {{G_{t,\eta }} - {y_{t,\lambda ,\eta }}} \right) + {\varpi _2}{r_2}\left( {{P_{t,\lambda ,\eta }} - {y_{t,\lambda ,\eta }}} \right), \\ y_{t + 1,\lambda ,\eta }^{} = {y_{t,\lambda ,\eta }} + \dfrac{{ {{{\left( {v_{t,\lambda ,\eta }^{}} \right)}^2} - {{\left( {v_{t - 1,\lambda ,\eta }^{}} \right)}^2}} }}{{2a}},\\ a = 2590(1 - {{\mathrm{e}}^{ - \lg \;t}}), \end{array} \right. (11)

    式中: v_{t,\lambda ,\eta }^{} 为第 t 次迭代的速度; {\varpi _1} {\varpi _2} 分别为 {G_{t,\eta }} {P_{t,\lambda ,\eta }} 的认知因子; a 为变色龙的加速度,随着 t 的变化而非线性变化.

    整体的变色龙群优化算法流程如图5所示.

    图  5  变色龙群算法流程
    Figure  5.  Flow chart of CSA

    选取山区某铁路路基为研究对象,全线路基包含隧道出口段和进口段两部分,共1178 m,地处区域地势地形极其复杂,高差大,工程施工多为高路堤深路堑形式,施工难度巨大. 此外,由于地处峡谷边缘区,峡谷风剧烈并伴随强降水,易引发水土流失,路基边坡防护设计任务艰巨,路基区域附近多湿地、暗河等水体,植物、动物及附近居民需求高,环境保护要求高,脆弱的环境因素会对路基工程的建设活动产生重大影响,有必要挖掘关键要素,并对其进行优化调控.

    结合现场实测数据及相关工程文件,依据鲍学英等[15]对于绿色关键要素的识别方法,构建多层复杂关联网络,通过式(1) ~ (5)分别得到路基工程及环境系统绿色要素的综合重要度,路基工程内部绿色要素 {x_1} {x_{19}} 的综合重要度依次为0.06、0.03、0.11、0.09、0.03、0.02、0.03、0.02、0.04、0.04、0.08、0.06、0.10、0.11、0.04、0.04、0.04、0.04、0.03,环境绿色要素 {x_{20}} {x_{26}} 的综合重要度依次为0.07、0.21、0.08、0.10、0.28、0.05、0.21,因此,识别影响“路基-环境”相容共生的绿色关键要素为路堑开挖尺寸、路堤填筑尺寸、支挡结构设计、支挡结构布设、工程防护结构设计、边坡工程防护布设、植物防护结构设计、边坡植物防护布设、生态声屏障结构设计、生态声屏障布设.

    通过式(1)、(6)求得路基工程与环境的耦合协调度值为0.672,通过式(8)求得绿色耦合效率值为1.65,表示该“路基-环境”处于超载状态,耦合协调度也相对较低,处于初级协调,需要降低工程上投入,并将绿色关键要素作为主控变量进行优化调控,促使“路基-环境”保持满载的状态下耦合协调度向上提升.

    基于上述分析,结合式(6) ~ (8),设置求解模型的目标函数及约束条件,基于MATLAB平台并使用式(9) ~ (11)编写变色龙群算法程序,设置相关参数,种群数量为400个,最大迭代次数为5000次. 同时,为了对比CSA算法的准确性和科学性,将SA和新型麻雀搜索算法(SSA)进行对比,在同一平台运行20次,得到最终的耦合协调度最优值及相关的变量参数. 取其中具有代表性的一次结果,如图6所示.

    图  6  CSA、SA、SSA算法对比结果
    Figure  6.  Comparison results of CSA, SA, and SSA

    相比SA及SSA算法结果,CSA算法的求解效果最接近于最优值,这是由于在建模过程中,CSA算法模拟了变色龙眼睛360°旋转特征,使猎物被发现的概率增加,并能使用高速发射的粘性舌头高效率捕捉猎物,具备较强的局部最优脱困能力,精准获取全局最优值,从而求解效果较好.

    CSA算法求得的山区铁路“路基-环境”绿色耦合效率值为1.04、耦合协调度最优值为0.77,表明路基工程建设活动对环境的影响处于环境承载力临界区间,“路基-环境”协调水平由初级协调提升到中级协调状态.

    图7所示,整体主控变量优化情况大致分为3个等级. {x_{10}} {x_{11}} {x_{12}} {x_{14}} 优化比例较高,达到70%左右,是影响该“路基-环境”相容共生最主要绿色关键要素,应优先对其进行调控; {x_3} {x_4} {x_9} 优化比例达到45%左右,是次要需要进行调控的要素; {x_{13}} {x_{15}} {x_{16}} 优化比例最小,处于20%区间,表明在这些绿色关键要素领域路基建设活动较为成熟,调控幅度较小. 基于此,在满足工程可靠性且不浪费设计资源的前提下,通过以上主控变量的分级优化调控,促使“路基-环境”耦合协调状态向相容共生转变.

    图  7  主控变量最优解优化比例
    Figure  7.  Optimization proportion of optimal solution for main control variables

    1) 将铁路“路基-环境”复杂的耦合状况通过耦合魔方游戏模型进行展现,解释了铁路“路基-环境”的最优耦合状态为满载条件下协调水平最优,提出了铁路“路基-环境”调控框架;运用耦合协调度模型和施压承载模型对某山区铁路路基工程与环境之间的耦合状况进行分析,表明该“路基-环境”协调水平处于初级协调和“超载”状态,存在一定的优化空间.

    2) 构建以山区铁路“路基-环境”协调水平最优为目标,各绿色要素函数关系及满载条件约束的耦合调控优化模型,采取变色龙群算法对主控变量进行优化求解,得到满载条件下最优耦合协调值和路堑开挖尺寸、路堤填筑尺寸、支挡结构设计、支挡结构布设、工程防护结构设计、边坡工程防护布设、植物防护结构设计、边坡植物防护布设、生态声屏障结构设计、生态声屏障布设主控变量的优化比例,为山区铁路“路基-环境”调控优化提供了思路.

    3) 将变色龙群优化算法和传统算法SA及新型算法SSA进行结果对比,结果表明CSA算法全局搜索能力最优,论证了该模型的有效性与科学性.

  • 图 1  铁路“路基-环境”交互影响关系

    Figure 1.  Interaction between railway roadbed and environment

    图 2  铁路“路基-环境”绿色要素指标体系

    Figure 2.  Index system for green elements of railway roadbed and environment

    图 3  魔方游戏模型

    Figure 3.  Rubik’s cube game model

    图 4  魔方游戏施压承载模型

    Figure 4.  Rubik’s cube game pressure bearing model

    图 5  变色龙群算法流程

    Figure 5.  Flow chart of CSA

    图 6  CSA、SA、SSA算法对比结果

    Figure 6.  Comparison results of CSA, SA, and SSA

    图 7  主控变量最优解优化比例

    Figure 7.  Optimization proportion of optimal solution for main control variables

    表  1  耦合协调度划分标准

    Table  1.   Classification criteria for ccoupling coordination degree

    第1分类划分标准第2分类及划分标准
    协调发展类0.9<D≤1.0优质协调
    0.8<D≤0.9良好协调
    0.7<D≤0.8中级协调
    0.6<D≤0.7初级协调
    过渡发展类0.5<D≤0.6勉强发展
    0.4<D≤0.5濒临失调
    失调衰退类0.3<D≤0.4轻度失调
    0.2<D≤0.3中度失调
    0.1<D≤0.2严重失调
    0<D≤0.1极度失调
    下载: 导出CSV
  • [1] 李艳鸽,姜楠,韩征,等. “绿色铁路设计” 核心理念的界定及应用探讨——以某西部铁路工程为例[J]. 铁道科学与工程学报,2022,19(11): 3439-3446.

    LI Yange, JIANG Nan, HAN Zheng, et al. Discussion on definition and application of the core concept of Green Railway Design―a case study of a Tibet railway project[J]. Journal of Railway Science and Engineering, 2022, 19(11): 3439-3446.
    [2] SU L J, HU B L, XIE Q J, et al. Experimental and theoretical study of mechanical properties of root-soil interface for slope protection[J]. Journal of Mountain Science, 2020, 17(11): 2784-2795. doi: 10.1007/s11629-020-6077-4
    [3] SUKMAK P, SUKMAK G, HORPIBULSUK S, et al. Improved mechanical properties of cement-stabilized soft clay using garnet residues and tire-derived aggregates for subgrade applications[J]. Sustainability, 2021, 13(21): 11692.1-11692.22.
    [4] 郑雨茜,鲍学英. 绿色施工下的铁路路基爆破效益分析[J]. 爆破,2019,36(1): 147-154. doi: 10.3963/j.issn.1001-487X.2019.01.022

    ZHENG Yuxi, BAO Xueying. Blasting benefits analysis of railway subgrade under green construction[J]. Blasting, 2019, 36(1): 147-154. doi: 10.3963/j.issn.1001-487X.2019.01.022
    [5] 李雨浓,鲍学英. 川藏铁路绿色施工节材措施综合效果评价研究[J]. 铁道科学与工程学报,2021,18(6): 1613-1621.

    LI Yunong, BAO Xueying. Study on comprehensive effect evaluation of material-saving measures for green construction of Sichuan-Tibet railway[J]. Journal of Railway Science and Engineering, 2021, 18(6): 1613-1621.
    [6] 鲍学英,李亚娟,胡所亭,等. 铁路桥隧工程技术接口的关键要素均衡优化[J]. 浙江大学学报(工学版),2022,56(3): 558-568.

    BAO Xueying, LI Yajuan, HU Suoting, et al. Tradeoff optimization of key elements of technical interface of railway bridge-tunnel engineering[J]. Journal of Zhejiang University (Engineering Science), 2022, 56(3): 558-568.
    [7] 李伟,蒲浩,赵海峰,等. 基于分步编码改进遗传算法的铁路智能选线[J]. 西南交通大学学报,2013,48(5): 831-838. doi: 10.3969/j.issn.0258-2724.2013.05.008

    LI Wei, PU Hao, ZHAO Haifeng, et al. Intelligent railway aignment optimization based on stepwise encoding genetic algorithm[J]. Journal of Southwest Jiaotong University, 2013, 48(5): 831-838. doi: 10.3969/j.issn.0258-2724.2013.05.008
    [8] 段晓晨,喇海霞,胡天明,等. 桥梁工程运维成本三维非线性智能控制研究[J]. 铁道工程学报,2020,37(9): 102-107. doi: 10.3969/j.issn.1006-2106.2020.09.018

    DUAN Xiaochen, LA Haixia, HU Tianming, et al. Research on the 3D nonlinear intelligent control for bridge engineering operation and maintenance cost[J]. Journal of Railway Engineering Society, 2020, 37(9): 102-107. doi: 10.3969/j.issn.1006-2106.2020.09.018
    [9] 薛新功,李伟,蒲浩. 铁路线路智能优化方法研究综述[J]. 铁道学报,2018,40(3): 6-15. doi: 10.3969/j.issn.1001-8360.2018.03.002

    XUE Xingong, LI Wei, PU Hao. Review on intelligent optimization methods for railway alignment[J]. Journal of the China Railway Society, 2018, 40(3): 6-15. doi: 10.3969/j.issn.1001-8360.2018.03.002
    [10] 肖建庄,夏冰,肖绪文. 工程结构可持续性设计理论架构[J]. 土木工程学报,2020,53(6): 1-12.

    XIAO Jianzhuang, XIA Bing, XIAO Xuwen. Theoretical framework for sustainability design of engineering structures[J]. China Civil Engineering Journal, 2020, 53(6): 1-12.
    [11] 王思源,赵敏敏,闫晶,等. 川藏铁路西藏昌都段生态保护重要性评价[J]. 现代地质,2021,35(1): 234-243.

    WANG Siyuan, ZHAO Minmin, YAN Jing, et al. Evaluation on the importance of ecological protection in Changdu section of the Sichuan-Tibet railway[J]. Geoscience, 2021, 35(1): 234-243.
    [12] 彭建兵,崔鹏,庄建琦. 川藏铁路对工程地质提出的挑战[J]. 岩石力学与工程学报,2020,39(12): 2377-2389.

    PENG Jianbing, CUI Peng, ZHUANG Jianqi. Challenges to engineering geology of Sichuan: Tibet railway[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(12): 2377-2389.
    [13] 蔡迪文,张克存,安志山,等. 青藏铁路格拉段风动力环境及其对铁路沙害的影响[J]. 中国沙漠,2017,37(1): 40-47. doi: 10.7522/j.issn.1000-694X.2015.00272

    CAI Diwen, ZHANG Kecun, AN Zhishan, et al. Wind energy environments and its impacts on railway sand hazards along gerlha section of the Qinghai-Tibet railway[J]. Journal of Desert Research, 2017, 37(1): 40-47. doi: 10.7522/j.issn.1000-694X.2015.00272
    [14] 国家铁路局. 铁路路基设计规范:TB 10001—2016[S]. 北京:中国铁道出版社,2017.
    [15] 国家铁路局. 铁路工程环境保护设计规范:TB 10501—2016[S]. 北京:中国铁道出版社,2016.
    [16] 刘海猛,方创琳,李咏红. 城镇化与生态环境“耦合魔方” 的基本概念及框架[J]. 地理学报,2019,74(8): 1489-1507. doi: 10.11821/dlxb201908001

    LIU Haimeng, FANG Chuanglin, LI Yonghong. The coupled human and natural cube: a conceptual framework for analyzing urbanization and eco-environment interactions[J]. Acta Geographica Sinica, 2019, 74(8): 1489-1507. doi: 10.11821/dlxb201908001
    [17] 鲍学英,沈杜华,郭海东,等. 山区铁路工程 “路基-环境” 绿色关键要素识别研究[J]. 安全与环境学报,2024,24(2): 787-795.

    BAO Xueying, SHEN Duhua, GUO Haidong, et al. Study on the identification of “subgrade-environment” green key elements in mountain railway engineering[J]. Journal of Safety and Environment, 2024, 24(2): 787-795.
    [18] 莫俊文,滕仓国,李甲,等. 基于熵权—二维云模型的高铁建设工程系统韧性评价[J]. 铁道科学与工程学报,2022,19(1): 26-33.

    MO Junwen, TENG Cangguo, LI Jia, et al. Resilience evaluation of high-speed railway construction engineering system based on entropy weight-two dimensional cloud model[J]. Journal of Railway Science and Engineering, 2022, 19(1): 26-33.
    [19] 鲍学英,肖岚月,万炳彤,等. 节材视角下山区铁路与材料资源协调发展研究[J]. 铁道科学与工程学报,2023,20(11): 4367-4376.

    BAO Xueying, XIAO Lanyue, WAN Bingtong, et al. Coordinated development of mountain railways and material resources from perspective of material saving[J]. Journal of Railway Science and Engineering, 2023, 20(11): 4367-4376.
    [20] 任慧君. 黄河流域9省区能源产业链与资源环境耦合关系研究[D]. 北京:中国矿业大学(北京),2020.
    [21] 王恒,方兰. 中国水-能源-粮食纽带系统安全水平与全要素生产率时空耦合协调关系分析[J]. 水资源保护,2023,39(1): 150-157.

    WANG Heng, FANG Lan. Spatial-temporal coupling coordination relationship between the security level of water-energy-food nexus system and total factor productivity in China[J]. Water Resources Protection, 2023, 39(1): 150-157.
    [22] 贺振霞,鲍学英,胡所亭,等. 铁路工程多系统技术接口复合网络信息质量耦合协调分析[J]. 土木工程学报,2022,55(11): 118-130.

    HE Zhenxia, BAO Xueying, HU Suoting, et al. Coupling and coordination analysis of information quality of railway engineering multi-system technical interface composite network[J]. China Civil Engineering Journal, 2022, 55(11): 118-130.
    [23] 国家铁路局. 铁路声屏障工程设计规范:TB 10505—2019[S]. 北京:中国铁道出版社,2019.
    [24] 中华人民共和国水利部. 土壤侵蚀分类分级标准:SL 190—2007[S]. 北京:中国水利水电出版社,2008.
    [25] 国家质量监督检验检疫总局,中国国家标准化管理委员会. 环境空气质量标准:GB 3095—2012[S]. 北京:中国环境科学出版社,2016.
    [26] 环境保护部. 声环境质量标准:GB 3096—2008[S]. 北京:中国环境科学出版社,2008.
    [27] 国家环境保护总局,国家质量监督检验检疫总局. 地表水环境质量标准:GB 3838—2002[S]. 北京:中国环境科学出版社,2002.
    [28] 国家质量监督检验检疫总局,中国国家标准化管理委员会. 爆破安全规程:GB 6722—2014[S]. 北京:中国标准出版社,2015.
    [29] 刘坚,李泽华,白和强,等. 露天矿爆破粉尘排放量影响因素分析及试验研究[J]. 爆破,2017,34(4): 169-174. doi: 10.3963/j.issn.1001-487X.2017.04.030

    LIU Jian, LI Zehua, BAI Heqiang, et al. Research on influence factors of blasting dust emission of open-pit mine[J]. Blasting, 2017, 34(4): 169-174. doi: 10.3963/j.issn.1001-487X.2017.04.030
    [30] 鲍学英,郑雨茜,王起才. 绿色施工下的深路堑施工群综合承载度分析[J]. 浙江大学学报(工学版),2019,53(3): 482-491.

    BAO Xueying, ZHENG Yuxi, WANG Qicai. Construction group comprehensive bearing capacity analysis of deep cutting under green construction[J]. Journal of Zhejiang University (Engineering Science), 2019, 53(3): 482-491.
    [31] BRAIK M S. Chameleon Swarm Algorithm: a bio-inspired optimizer for solving engineering design problems[J]. Expert Systems with Applications, 2021, 174: 114685.1-114685.25.
    [32] KARYPIS G, HAN E H, KUMAR V. Chameleon: hierarchical clustering using dynamic modeling[J]. Computer, 1999, 32(8): 68-75. doi: 10.1109/2.781637
  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  52
  • HTML全文浏览量:  26
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-07-13
  • 修回日期:  2023-10-11
  • 网络出版日期:  2025-02-25

目录

/

返回文章
返回