《1 引言》
1 引言
土石坝
土石坝工程
国外对此问题专门的研究不多见
综观以上研究成果, 存在的问题在于未能将土石方调配当成一个受时间、空间影响的动态调配系统, 仅考虑了其在空间“数量”上的平衡, 而未能考虑其在时间上的综合平衡;或者没有深入考虑调配中的一些约束条件及人为干预因素, 造成调配结果与实际操作有较大出入;在调配成果的表达上, 仅限于一般的文本输出, 形式过于简单, 不能给人以直观的形象表达;同时编写的计算机程序通用性不强。所有这些使其应用受到极大限制。
《2 研究方法》
2 研究方法
由于土石方调配系统是一个受时间-空间约束限制的复杂动态调配系统, 如果采用运筹学中只考虑空间因素, 即只考虑运距及供、需量的一般求解调配问题优化的方法则无法解决。因此, 本文采用了一种适合于求解动态调配系统的数学模型和求解方法。其思路为:将整个施工过程划分为若干时段进行递阶计算, 并将状态转移关系引入中转料场的平衡约束, 从而表达时-空转移条件, 通过时-空转移将动态问题静态化。再采用静态规划的方法求解并进行调配系统的综合平衡规划, 求出最优的调配方案。这种方法不但概念明确, 直观易行, 而且便于针对工程实际综合考虑各种影响因素及制约条件。此外, 它将枢纽施工全过程作为一个整体进行优化, 可协调各时期相互制约的关系, 便于监控全过程。
本文是在天津大学已有研究成果的基础上, 采用系统分析的方法, 根据工程实际施工进度, 以全过程系统调配成本最低为目标, 深入而细致地考虑了调配中的一些约束条件和人为干预因素, 建立调配系统优化模型;在成果的表达上, 本文借助于GIS强大的空间信息处理功能
《3 土石方动态调配数学模型》
3 土石方动态调配数学模型
《3.1 模型准备》
3.1 模型准备
在建立此问题的数学模型之前, 先对调配的料源作一说明。在这里, 将料源划分为供料源与受料源, 并依据其不同的性质, 供料源可细分为建筑物开挖供料源、直接开采料场供料源以及中转料场 (堆渣场) 供料源。受料源细分为坝体受料区受料源及中转料场受料源。可以看出, 中转料场既是供料源, 又是受料源。由于程序中的开挖方量均为有效利用量, 故在此不考虑弃渣的料, 也就不存在专门的弃渣场。
土石方调配的目的即为确定各施工时段由各供料源到各受料源的调配方量。任意时段土石方调配示意图如图1所示。从图中可以看出, 任意时段的土石方调配包括了运输路径的任意选择, 但不包括从一个中转料场到另一个中转料场的调配, 这是根据实际操作的要求而加入的一个约束条件。
《图1》
图1 任意时段土石方调配示意图 Fig.1 The sketch map of earth allocation in any stage
在建立土石方动态调配系统的数学模型时遵循了下面两个假定:一是所有的料场规划 (包括料场数量、位置及储量等) 已经完成, 即土石方调配是在既定的料场规划下进行的;二为直接开采料场任意时段开采的料没有数量与质量上的约束限制, 即认为直接开采料场可满足任何坝体受料区颗粒级配的要求, 且开采的量可满足任意需求量的要求。
《3.2 数学模型》
3.2 数学模型
为了描述方便, 先对各料源的数量作统一规定:k个建筑物开挖料源, l个直接开采料场, m个中转料场, n个坝体受料区。此问题的线性规划数学模型可简单描述为:土石方调配就是在满足各项施工进度和填筑材料质量与级配要求的前提下, 满足:1) 合理调配各建筑物、开采料场、中转料场的土石料;2) 全过程调配成本最低的目标, 以确定各时期、各区域相应的调配规划。
《3.2.1 目标函数 其目标函数为:》
3.2.1 目标函数 其目标函数为:
这里, Z为目标函数, 即总的调配费用 (万元) , C
《3.2.2 约束条件》
3.2.2 约束条件
1) 供量约束
其中, a
2) 需量约束
这里, b
3) 中转料场平衡约束 (时空转移条件)
其中, R
将式 (7) 代入式 (5) , 则可得式 (8) :
4) 优先料源约束。
优先料源约束指的是某些开挖出的土石料在优先满足坝体相应颗粒级配的部位填筑要求下, 多余的土石料在其他受料源进行统一调配。如有多个优先选择的料源, 则按照选定的优先级别进行调配。即
如地下工程开采的土石料优先用于主堆石区的填筑, 如有多余则满足次堆石区的填筑等, 依此类推。
5) 中转料场出料约束。
中转料场的出料约束指的是从中转料场调出的料要满足中转料场已有供料源库存料的要求, 即
这里, X
6) 不可能运输路径约束。
不可能运输路径有两层含义, 它既反映了各供料源与受料源的使用规划、开挖填筑项目的关系, 同时也隐含反映了道路系统及料物类型与开挖填筑的关系, 是进行土石方动态调配最重要也是最基础的数据。根据前面叙述过的供料源与受料源的详细划分, 可建立起一个二维矩阵, 其行表示供料源, 列表示受料源。行列相交的元素取值满足以下要求:不符合设计规定的运输路径取任意一负值如-1, 如开采出的石料不可能运输到粘土层进行填筑, 则其运输距离取为-1;对于符合设计规定的取实际运距, 如表1为某工程的土石方调配运距表。
对于某不可能的运输路径i→j, 其调配量为0, 即
《7) 中转料场调配结束的零约束。》
7) 中转料场调配结束的零约束。
所有中转料场在调配结束时其储量应为零, 即
其中, A
8) 非负约束。要求线性规划的所有变量为非负数, 即
表1 土石方调配运距表 Table 1 The table of haul distance of possible cut-and-fill pair m
《表1》
供料源 | 坝体受料区 | 中转料场 | ||||||||
ⅡA | ⅢA | ⅢB | ⅢC | ⅣA | ⅣB | 中转场1 | 中转场2 | 中转场3 | ||
(1) 溢洪道开挖 | -1 | -1 | -1 | -1 | 652 | -1 | 1 035 | -1 | -1 | |
(2) 泄洪洞明挖 | -1 | -1 | -1 | -1 | 726 | -1 | 816 | -1 | -1 | |
(3) 引水渠明挖 | -1 | -1 | -1 | -1 | 451 | -1 | -1 | -1 | 2 595 | |
(4) 出现平台明挖 | -1 | -1 | -1 | -1 | 466 | -1 | -1 | -1 | 2 202 | |
建 筑 物 开 挖 料 源 | (5) 尾水渠明挖 | -1 | -1 | -1 | -1 | 1 016 | -1 | -1 | -1 | 783 |
(6) 引水及竖井等洞挖 | -1 | 505 | 505 | -1 | 505 | -1 | -1 | 2 396 | -1 | |
(7) 主厂房及主变室 | -1 | 920 | 920 | -1 | 920 | -1 | -1 | 394 | -1 | |
(8) 厂房的其它洞室 | -1 | 1 020 | 1 020 | -1 | 1 020 | -1 | -1 | 425 | -1 | |
(9) 泄洪洞洞挖 | -1 | 883 | -1 | -1 | 883 | -1 | 1 180 | -1 | -1 | |
(10) 导流洞及灌浆洞洞挖 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | 1 438 | -1 | |
开采料场 | (11) 开采料场 | 4 716 | 4 716 | 4 716 | 4 716 | 4 716 | 4 716 | 3 658 | 2 546 | 2 341 |
(12) 中转场1 | 816 | 816 | 816 | 816 | 816 | 816 | -1 | -1 | -1 | |
中 转 料 场 | (13) 中转场2 | 516 | 516 | 516 | 516 | 516 | 516 | -1 | -1 | -1 |
(14) 中转场3 | 616 | 616 | 616 | 616 | 616 | 616 | -1 | -1 | -1 |
综合以上式 (1) ~式 (16) (式5, 7已合并为式8) , 便构成了土石方动态调配系统的数学模型。
《3.2.3 土石方动态调配程序设计》
3.2.3 土石方动态调配程序设计
上述土石方动态调配系统的数学模型是一个线性规划模型, 对线性规划的求解一般有单纯形法、图解法以及改进的单纯形法, 本文采用改进的单纯形法对问题进行求解
《4 可视化的土石方动态调配》
4 可视化的土石方动态调配
鉴于以往调配成果表达上的不足, 本文采用目前强大的空间地形处理工具地理信息系统技术对其进行了处理。地理信息系统 (GIS) 是一门近年来迅速发展起来的、介于地球科学与信息科学之间的交叉学科, 也是一门地学空间数据与计算机技术相结合的新型空间信息技术
成果可视化表达的第一步是建立整个施工场地的三维数字地形模型
土石方调配涉及到与地理空间位置或特征相关联的空间地物, 包括:地表和地下的各种人工和自然现象, 依据地理实体的分布特征、空间特征及其相互关系, 可将空间对象抽象为三维空间中的点、线、面、体等四种不同的基本图形元素
最后运用动画技术显示整个调配过程。由优化调配计算程序得出的调配成果按照一定格式写入GIS数据库中, 并进行处理, 形成各条调配路径, 由GIS以图形的形式表现出来。通过众多这种画面的连续播放形成了动画, 将调配的各时段依次联接起来, 于是就演示了土石方调配的整个时段, 达到直观、形象的目的。
《图2》
图2 土石方动态调配规划程序流程图 Fig.2 The program flow chart of earth allocation
《5 工程实例》
5 工程实例
某混凝土面板堆石坝工程最大坝高185.5 m, 坝顶长423.75 m, 宽高比为2.28, 属狭窄河床高面板堆石坝。坝体总填筑方量为961.1×104 m3。
本工程有土料场四个, 石料场最终确定为一个, 其运距较远, 但料性及容量均能满足要求。同时可利用料为343.8×104 m3, 尚需从石料场开采堆石料593.6×104 m3 (压实方) , 考虑各种损耗并折算成自然方为577.13×104 m3。
根据工程施工进度的按排, 坝体填筑分五期进行, 考虑前期准备及中间停歇, 坝体共分八期进行施工, 坝体分区如图3所示。从施工进度表上确定各时期的开挖量与填筑量, 并依实际工程的要求建立调配运距表 (如表1所示) 及确定各种约束条件, 然后利用软件进行优化计算 (软件的部分应用界面如图4和图5所示) , 即可确定最优的调配方案。调配结果表明直接上坝率偏低, 仅为37%, 其原因为开挖利用料的主要填筑区开始填筑时间较迟, 开挖料需转运堆存, 造成直接上坝率较低。因数据量过大, 详细调配结果在此从略。
将调配成果在GIS环境下进行处理, 可获得可视化的输出, 图6为时段四的调配过程示意图。
《6 结语》
6 结语
土石方调配是土石坝设计与施工中必须考虑的一项重要内容, 是进行其他各项工作的基础。本文在现有研究成果的基础上, 采用系统分析的方法,
《图3》
图3 某堆石坝坝体分区示意图 Fig.3 The sketch map of zones of the rock-fill dam
《图4》
图4 软件运距约束对话框 Fig.4 The dialogue of hauling distance of the software
《图5》
图5 软件料源优先约束对话框 Fig.5 The dialogue of priority constraints of accepting sources of the software
将一个受时间-空间约束的动态问题静态化, 利用线性规划对问题进行求解, 简化了问题的复杂性。同时, 建立的数学模型灵活性强, 可进行一些人为的干预;调配成果的可视化表达能使调配过程直观形象地展现出来;开发的软件不仅可以满足土石坝施工中料物调配的需要, 还适用于其他相类似工程的土石方调配, 如港口工程、公路工程等, 因此, 具有较强的通用性。
《图6》
图6 施工时段四的调配过程示意图 Fig.6 The demo graph of earth allocation process in 4th stage
由于影响土石方调配的因素众多, 且随机事件时有发生, 如何更加真实地体现其受到的众多约束以及人为干预因素, 即如何更进一步增强软件的通用性, 将是今后进行此项研究工作的重要内容。