发明内容
为了解决上述技术问题,本发明提出一种山区铁路施工便道线路网智能布设方法,所述布设方法用于实现山区铁路施工便道线路网智能布设。本发明通过不断伸展路线来形成施工便道线路网,并实时协调多种约束,从而获取最优线路网方案。
本发明解决上述技术问题的技术方案是:
一种山区铁路施工便道线路网智能布设方法,包括以下步骤:
S1:获取目标研究区域,采集工程信息,所述工程信息包括工程结构信息、基础地理信息以及施工节点信息;
S2:将步骤S1所述研究区域划分为若干个网格,并添加网格属性;
S3:基于步骤S1采集的信息,构建山区铁路施工便道线路网优化模型;具体步骤如下:
所述山区铁路施工便道线路网优化模型包括设计变量、约束条件以及目标函数;所述设计变量包括道路平面线形信息以及道路纵断面线形信息;所述平面线形信息包括平面交点坐标、平面圆曲线半径长度、平面缓和曲线长度;所述纵断面线形信息包括变坡点里程、高程以及竖曲线半径长度;所述约束条件至少包括禁区约束、高程约束、结构物约束以及道路几何约束;
所述约束条件的具体检测步骤为:
(1)禁区约束检测
检测路径中的网格是否属于禁区,即判断网格的Is_ForbidZone属性值是否为1;
(2)高程约束检测
检测路径中每相邻路径网格的纵坡坡度是否小于最大纵坡坡度:
其中ρ为所述局部路径的纵坡坡度;ρmax为最大纵坡坡度;xw1,yw1,zw1为网格Ww1的中心空间坐标;xw2,yw2,zw2为网格Ww2的中心空间坐标;所述ΔH为网格Ww1和网格Ww2之间的高度差;所述ΔS为网格Ww1和网格Ww2之间的距离;所述网格Ww1和网格Ww2为路径中相邻局部路径网格;
(3)结构物约束检测
检测路径的路基、桥梁、隧道结构是否满足约束:
桥梁设置分界高度HBridge和隧道设置分界高度HTunnel以及路径与地面线间的高度差ΔHD,判断出路径需构建的路基、桥梁或隧道结构,接着对所述结构物进行约束检测:
其中lT为所述现有最优路径隧道长;lTmax为最大隧道长;hS为所述现有最优路径路基边坡高;hSmax为最大路基边坡高;lS为所述现有最优路径便桥单孔跨径;lsmax为最大便桥单孔跨径;hP为所述现有最优路径便桥墩台高;hPmax为最大便桥墩台高;
(4)道路几何约束检测
检测路径的曲线半径、圆曲线长度、缓和曲线长度、坡段长度、纵坡、平均纵坡、竖曲线是否满足设计规范要求;
所述目标函数为建设施工便道工程费;所述目标函数对应的优化目标为:所需工程费用最小;其中,所述工程费用包括:土石方工程费用、修建便桥工程费用、修建隧道工程费用、土地占用费用以及道路铺面费用;
S4:求解施工干线,所述施工干线用于连接主要施工节点与既有道路网;
S5:基于步骤S4得到的施工干线路径,更新网格属性;
S6:基于步骤S5更新后的既有道路网,求解施工支线,所述施工支线用于连接次要施工节点与既有道路网;
S7:基于步骤S4得到的施工干线以及步骤S6得到的施工支线,求解施工便道线路网,将所述施工通道线路网路径曲线化,最终得到最优施工通道线路网。
本发明的有益效果:
1、本发明的布设方法采用二次建网法,首先生成施工干线,接着将施工干线并入进既有道路网,生成施工支线,最后形成施工便道线路网。与人工经验相比,提高了设计效率及准确性。
2、本发明将弗洛伊德算法应用到施工便道路径探索中去,通过动态规划不断更新路径矩阵,全面考虑多种路径的可能性,并实时充分检测各种约束,使生成的施工便道线路网更加准确可靠。
3、本发明将新生成的施工干线接入进既有路网,能够更加充分的利用资源,同时使设计经济合理。
4、本发明具有智能化、自动化程度高的优点,具有较好的推广价值。
具体实施方式
下面将结合附图,对本发明的技术方案做进一步的详细说明
如图1所示,图1为本发明的流程图。本发明的山区铁路施工便道线路网智能布设方法包括以下步骤:
S1:获取目标研究区域,采集工程信息;
工程信息包括工程结构信息、基础地理信息以及施工节点信息;
工程结构信息包括:(a)施工便桥桥面宽;(b)桥梁设置分界高度HBridge;(c)隧道跨径长度;(d)隧道设置分界高度HTunnel;(e)路基宽度;(f)路基边坡坡度;
基础地理信息包括:(a)用地类型(农田、建设用地);(b)地价信息;(c)地物信息(既有道路);(d)禁区类型(地质灾害频发区、自然保护区);
施工节点信息包括主要施工节点以及次要施工节点;其中,主要施工节点包括隧道进口端洞门位置处、隧道出口端洞门位置处以及竖井、斜井施工节点,次要施工节点包括项目部、弃渣场以及拌合站;
S2:如图2所示,将研究区域划分为若干个网格,并添加网格属性;
正方形网格预定边长为d,步骤S1获取的整个研究区域可划分为I*J个正方形规则网格,共计N个网格(即N=I*J),自所述研究区域左上角起,各个网格依次记为Wi(i=1,2,3,...,N),将获取的基础地理信息以及施工节点信息添加进网格属性中,相应各网格中心空间坐标为(xi,yi,zi),研究区域内N个网格Wi(i=1,2,3,...,N)的属性可表示为矩阵,其对应的行向量为:Di(xi,yi,zi,ZoneTypei,Is_ForbidZonei,Is_RoadZonei,Is_FPointi,Is_SPointi);
S3:构建山区铁路施工便道线路网优化模型。
该山区铁路施工便道线路网优化模型包括设计变量、约束条件以及目标函数;
设计变量:
设计变量包括平面交点坐标、平面圆曲线半径长度、平面缓和曲线长度、变坡点里程、高程以及竖曲线半径长度,采用以下变量表示设计变量:
平面交点X坐标列向量:X=[X1,X2,X3,...,Xnp]T
平面交点Y坐标列向量:Y=[Y1,Y2,Y3,...,Ynp]T
平面圆曲线半径列向量:R=[R1,R2,R3,...,Rnp]T
平面缓和曲线列向量:l=[l1,l2,l3,...,lnp]T
变坡点里程列向量:K=[K1,K2,K3,...,Knz]T
变坡点高程列向量:H=[H1,H2,H3,...,Hnz]T
竖曲线半径列向量:RV=[RV1,RV2,RV3,...,RVnz]T
式中:np为平面交点数量;nz为纵断面变坡点数量;
约束条件:
(1)禁区约束检测
检测路径中的网格是否属于禁区,即判断网格的Is_ForbidZone属性值是否为1;
(2)高程约束检测
检测路径中每相邻路径网格的纵坡坡度是否小于最大纵坡坡度:
其中ρ为所述局部路径的纵坡坡度;ρmax为最大纵坡坡度;(xw1,yw1,zw1)为网格Ww1的中心空间坐标;(xw2,yw2,zw2)为网格Ww2的中心空间坐标;所述ΔH为网格Ww1和网格Ww2之间的高度差;所述ΔS为网格Ww1和网格Ww2之间的距离;所述网格Ww1和网格Ww2为路径中相邻局部路径网格;
(3)结构物约束检测
检测路径的路基、桥梁、隧道结构是否满足约束:
首先,基于步骤S1中所述桥梁设置分界高度HBridge和隧道设置分界高度HTunnel以及路径与地面线间的高度差ΔHD,判断出路径需构建的路基、桥梁或隧道结构,接着对所述结构物进行约束检测:
其中lT为所述现有最优路径隧道长;lTmax为最大隧道长;hS为所述现有最优路径路基边坡高;hSmax为最大路基边坡高;lS为所述现有最优路径便桥单孔跨径;lsmax为最大便桥单孔跨径;hP为所述现有最优路径便桥墩台高;hPmax为最大便桥墩台高;
(4)道路几何约束检测
检测路径的曲线半径、圆曲线长度、缓和曲线长度、坡段长度、纵坡、平均纵坡、竖曲线是否满足川藏铁路施工道路和施工供电工程勘察设计暂行规定的设计规范要求;
(5)运输时间约束检测
检测路径的运输时间是否满足约束:
其中t为货车单次运输时间;将路径的道路划分为K段,每段道路长度为Lk,每段道路货车行驶速度为vk;Tmax为货车单次运输最大等待时长;
目标函数:
目标函数为建设施工便道工程费CCR;目标函数对应的优化目标为:所需工程费用最小;其中,工程费用CCR包括:土石方工程费用CES、修建便桥工程费用CBR、修建隧道工程费用CTU、土地占用费用CR以及道路铺面费用CPR;
(1)土石方工程费用CES
式中:
为挖方单价;
为挖方体积;
为填方单价;
为填方体积;
为将土方移至垃圾填埋场的单位成本;
为移至垃圾填埋场的土方体积;
为从借用坑中移土的单位成本;
为从借用坑中移土的土方体积;
(2)修建便桥工程费用CBR
式中:UBRi为第i座桥梁每延米的费用;LBRi为第i座桥梁的长度;UAi为第i座桥梁的桥台每延米费用;NAi为第i座桥梁的桥台数量;
(3)修建隧道工程费用CTU
式中:UTUi为第i座隧道每延米的费用;LTUi为第i座隧道的长度;UD为单个隧道洞门的费用;nTU为隧道的数量;
(4)土地占用费用CR
CR=UR×nR×d2 (3.2.4)
式中:UR为单位面积土地占用费;nR为道路占据的网格数;d为网格边长;
(5)道路铺面费用CPR
CPR=UPR×WR×LR (3.2.5)
式中:UPR为单位面积的道路铺面费;WR为道路路面宽度;LR为道路路面长度;
综上,抽象出数学表达,得到如下施工便道线路网优化模型:
式中:
(1)X,Y,R,l,K,H,RV为设计变量;
(2)CCR=CES+CBR+CTU+CR+CPR为模型目标函数,详见(3.1)所述;
(3) si(x,y,z)≤0为禁区约束;hi(x,y,z)≤0为高程约束;ci(h,l)≤0为结构物约束;
ci(x,y,z,R,l,K)≤0为道路几何约束;ti(l,v)≤0为运输时间约束,详见(3.2)
所述;
S4:使用弗洛伊德(Floyd)算法求解施工便道干线布设,生成施工便道干线方案,所述施工便道干线用于连接主要施工节点与既有道路网,主要实现步骤如图3所示,具体步骤如下:
S4-1:初始化最小路径代价值矩阵Val和最短路径矩阵Path;
初始化N阶方阵Val的主对角线元素为0,非主对角元素为∞;
初始化N阶方阵Path的主对角线元素为-1,非主对角元素为相应元素第二个下角标值;
所述最小路径代价值矩阵Val为N阶方阵,其每个元素Val[i][j]表示网格Wi和Wj之间的最小路径代价值,其中,所述代价值为对应路径施工便道的工程费;所述最短路径矩阵Path为N阶方阵,其每个元素Path[i][j]表示起点为Wi、终点为Wj的两个网格之间的最短路径要经过网格的下角标值;
在计算所述Val矩阵代价值时,当两个网格的Is_RoadZone属性值均为1时,即两个网格处均为既有道路,取其代价值为0;
S4-2:将网格W1加入进路径中;
对新加网格W1进行禁区约束检测:
基于W1网格Is_ForbidZone属性值,若步骤S4-2所述网格W1为非禁区,则进入步骤S4-3;若步骤S4-2所述网格W1为禁区,则跳转至步骤S4-4;
S4-3:更新所述Val和Path两个矩阵;
比较路径Wi→Wj和路径Wi→W1→Wj的代价值,取路径代价值较小者为网格Wi到网格Wj的最短路径;所述路径Wi→Wj为起点为Wi、终点为Wj的路径;所述路径Wi→W1→Wj为起点为Wi,经网格W1然后到达终点Wj的路径;
接着,按上述最短路径更新Path矩阵;
进一步地,步骤S4-3包括高程约束检测、结构物约束检测、道路几何约束检测以及运输时间约束检测,当四个约束检测均通过时,才可更新所述Val和Path两个矩阵,进入步骤S4-4;否则跳过步骤S4-3,进入步骤S4-4;
S4-4:首次取m=2,将网格Wm加入进路径中;
得到路径Wi→…→Wm和路径Wm→…→Wj,所述路径Wi→…→Wm为起点为Wi、终点为Wm的两个网格之间的最短路径(路径所经网格的下标需不小于1且不大于m-1);所述路径Wm→…→Wj为起点为Wm、终点为Wj的两个网格之间的最短路径(路径所经网格的下标需不小于1且不大于m-1);
进一步地,步骤S4-4包括对新加网格Wm进行禁区约束检测:
基于网格Is_ForbidZone属性值,若步骤S4-4所述网格Wm为非禁区,则进入步骤S4-5;若步骤S4-4所述网格Wm为禁区,则跳转至步骤S4-6;
S4-5:更新所述Val和Path两个矩阵;
比较路径Wi→…→Wm→…→Wj和步骤S4-4前已得到的起点为Wi、终点为Wj的两个网格间的最短路径的代价值(路径所经网格的下标需不小于1且不大于m-1),取路径代价值较小者为网格Wi到网格Wj的最短路径(路径所经网格的下标需不小于1且不大于m);
接着,按上述最短路径更新Path矩阵;
进一步地,步骤S4-5包括高程约束检测、结构物约束检测、道路几何约束检测以及运输时间约束检测,当四个约束检测均通过时,才可更新所述Val和Path两个矩阵,进入步骤S4-6;否则跳过步骤S4-5,进入步骤S4-6;其中,基于步骤S4-5中按最短路径代价值更新的Path矩阵,对现有最优路径进行约束检测,所述现有最优路径指起点为网格Wi、终点为网格Wj的最小代价值路径,如图4所示,得到所述现有最优路径的步骤为:
(1)若Path[i][j]=-1,则没有从网格Wi到网格Wj的路线;
(2)若Path[i][j]=j,则从网格Wi直接到网格Wj即为最优路径;
(3)若Path[i][j]=t≠j,则从网格Wi到网格Wj的最优路径要先从网格Wi经网格Wt,然后基于Path矩阵查找Path[t][j]是否等于j,若Path[t][j]=j,则从网格Wi到网格Wj的最优路径接着从网格Wt经网格Wj;若Path[t][j]≠j,则循环查找直至最优路径走到网格Wj结束;
S4-6:将步骤S4-4中的m加1,重复步骤S4-4和步骤S4-5,直至研究区域内所有网格Wi(i=1,2,3,...,N)均被加入路径中;
S4-7:通过比较路径代价值,追溯出每个主要施工节点对应的最优施工便道干线路径,如图5所示;
进一步地,步骤S4-7中寻找每个主要施工节点对应的最优施工便道干线路径的具体步骤为:
基于所述Val和Path两个矩阵,最优干线路径依次以主要施工节点网格Wp为起点,比较Val[p][j](j=1,2,3,…,N,andj≠p)的值,取最小者为对应最优干线路径,其中,所述主要施工节点网格指Is_FPoint属性值为1的网格;所述最优干线路径指连接主要施工节点与既有道路,即所述网格Wj的Is_RoadZone属性值应为1;同理,根据本发明前文提到的寻找最优路径的方法,根据Path矩阵找到对应最优干线路径依次经过的网格,此处不再赘述。
S5:更新网格属性;
为充分利用道路资源,将步骤S4得到的所有施工便道干线加入进既有道路网,如图6所示,使布设更加经济合理,具体步骤为:基于步骤S4得到的施工便道干线路径,更新所述干线路径经过的网格,即Is_RoadZone属性值更新为1;
S6:如图7所示,使用弗洛伊德(Floyd)算法求解施工便道支线布设,所述施工便道支线用于连接次要施工节点与既有道路网,生成施工便道支线方案的具体步骤与步骤S4所述弗洛伊德算法求解施工干线方法相同,此处不再赘述;
S7:如图8所示,基于步骤S4得到的施工通道干线以及步骤S6得到的施工通道支线,求解施工便道线路网,采用弦切支距法将所述施工通道线路网路径曲线化,最终得到最优施工通道线路网。