1. 引言
龙舟坪镇位居清江中下游、东临宜都,北接宜昌,地处“一坝两库”(隔河岩电站大坝,高坝洲与隔河岩库区),是长阳土家族自治县县城所在地,是全县政治、经济、文化中心。境内318国道、清江黄金水道横贯东西,王渔省道纵穿南北。地势西高东低,最高海拔1740米,最低海拔80米,年平均气温18℃。龙舟坪镇管辖13个社区、14个行政村,镇人民政府驻津洋口村(如图1所示)。全镇国土面积340平方公里,山林面积9344公顷,耕地5257公顷,其中水田443公顷。近年来全镇紧紧围绕建设“清江旅游重镇、宜昌卫星名镇、湖北经济强镇”的目标,以旅游服务、商贸流通、建筑安装、酒精食品、蔬菜瓜果、水产养殖、花木药材为发展重点,逐步形成了“沿江兴旅游、城区强工贸、国道创特色、城郊建支柱”的经济发展格局[1]。在地下燃气管网建设方面,面临着资金投入大、技术要求高、协调难度大、安全隐患多、后期维护难、经济效益低、村民接受度低、政策支持不足等问题和挑战,需要从规划、建设与运营等多个环节进行系统性的设计和管理、确保为农村居民提供安全可靠的燃气供应[2]。
注:该图基于自然资源部标准底图服务网站下载的审图号为GS(2016)008号的标准地图制作,底图无修改。
Figure 1. Geographic location of 13 communities and 14 administrative villages in Longzhouping in satellite image
图1. 卫星影像图中龙舟坪13个社区14个行政村地理位置
本文从降低地下燃气管道建设成本出发,应用图论中的最小支撑树理论模型,结合GIS卫星影像测量技术,研究地下燃气管网的规划布局,力求使得管网总长最短的规划路径。
2. 最小支撑树模型建模
2.1. 问题描述
长阳县龙舟坪镇要沿13个社区14个行政村居民点规划布局规划地下燃气管网,依照沿公路网形成的各乡镇居民点分布图,完成地下燃气管网的规划布局,要求提供解决地下燃气管网路线优化方案,使得燃气管网地下铺设总路径最短。
2.2. 最小支撑树建模
设无向图G为龙舟坪各村(社区)形成的无向赋权图,权值
表示两村(社区) i与j之间的距离。根据无向赋权图构造邻接矩阵
(对称矩阵),27是镇辖村(社区)个数,每个村(社区)就是一个节点。对图
,
表示顶点集合,E表示边的集合,
。
求图G从顶点1开始的最小支撑树,其数学规划模型如下:
1) 目标函数为
2) 约束条件,对于起始点1,至少有一条出去的路,即
3) 对于其他节点,恰有一条路进入,即
4) 对于图G不形成圈,约束为
最小支撑树问题的0-1整数规划模型为[3] [4]
2.3. 龙舟坪各乡村(社区)空间地理数据的搜集
为方便图论网络规划建模,先对龙舟坪镇各村(社区)标号[5],编号如下表1。
Table 1. No. of villages (communities) in Longzhouping Town
表1. 龙舟坪镇各村(社区)编号
序号 |
地名 |
序号 |
地名 |
序号 |
地名 |
1 |
湖口社区 |
10 |
后山社区 |
19 |
三渔冲村 |
2 |
四冲湾社区 |
11 |
西寺坪社区 |
20 |
晒鼓坪村 |
3 |
龙门社区 |
12 |
板桥铺社区 |
21 |
厚丰溪村 |
4 |
花坪社区 |
13 |
白氏坪社区 |
22 |
两河口村 |
5 |
枫竹园社区 |
14 |
刘家坳村 |
23 |
王家棚村 |
6 |
津洋口社区 |
15 |
黄家坪村 |
24 |
土地坡村 |
7 |
丹水社区 |
16 |
王子石村 |
25 |
全伏山村 |
8 |
何家坪社区 |
17 |
合子坳村 |
26 |
郑家榜村 |
9 |
鸣凤社区 |
18 |
胡家棚村 |
27 |
观坪林场 |
在天地影像图中统一以各村(社区)党员活动中心WGS-84坐标系下的坐标位置为基点,沿道路测量得到龙舟坪镇各村(社区)距离,构成邻接矩阵[6] [7],如下表2 (限于篇幅只列出了有限的第1到27行、第1到9列)。
Table 2. Adjacency matrix formed by the distance between villages (communities) in Longzhouping Town (Colums 1~15) (Unit: km)
表2. 龙舟坪镇各村(社区)距离构成的邻接矩阵表(第1到15列) (单位:km)
起点终点 |
湖口社区 |
四冲湾区 |
龙门社区 |
花坪社区 |
枫竹园区 |
津洋口区 |
丹水社区 |
何家坪区 |
鸣凤社区 |
湖口社区 |
|
|
|
|
|
|
|
|
|
四冲湾区 |
1.2 |
|
|
|
|
|
|
|
|
龙门社区 |
1.3 |
0.64 |
|
|
|
|
|
|
|
花坪社区 |
2.3 |
1.1 |
1.3 |
|
|
|
|
|
|
枫竹园区 |
4.6 |
3.7 |
4.2 |
3.1 |
|
|
|
|
|
津洋口区 |
4.5 |
3.9 |
4.5 |
3.7 |
1.3 |
|
|
|
|
丹水社区 |
4.7 |
4.3 |
5 |
4.5 |
2.5 |
1.2 |
|
|
|
何家坪区 |
3.1 |
2.3 |
2.9 |
2.1 |
1.5 |
1.6 |
2.4 |
|
|
鸣凤社区 |
2.6 |
3.8 |
3.8 |
4.9 |
6.9 |
6.5 |
6.2 |
5.4 |
|
后山社区 |
1.3 |
1 |
1.6 |
1.8 |
3.3 |
3.2 |
3.5 |
1.8 |
3.7 |
西寺坪区 |
4.3 |
3.2 |
3.6 |
2.3 |
1.5 |
2.6 |
3.8 |
1.9 |
6.9 |
板桥铺区 |
4.4 |
5.6 |
5.5 |
6.7 |
8.6 |
8.1 |
7.7 |
7.1 |
1.8 |
白氏坪区 |
5.5 |
6.4 |
6.7 |
7.4 |
7.7 |
6.9 |
5.9 |
6.7 |
4.1 |
刘家坳村 |
8 |
8.3 |
8.9 |
9 |
7.8 |
6.5 |
5.3 |
7.3 |
7.8 |
黄家坪村 |
7.6 |
7.4 |
8 |
7.5 |
5.1 |
4 |
3.1 |
5.4 |
8.6 |
王子石村 |
9.5 |
9.2 |
9.9 |
9.2 |
6.5 |
5.6 |
4.9 |
7.2 |
10.7 |
合子坳村 |
6.6 |
6 |
6.6 |
5.7 |
2.7 |
2.1 |
2.3 |
3.7 |
8.5 |
胡家棚村 |
8.1 |
7.4 |
7.9 |
6.9 |
3.8 |
3.6 |
3.9 |
5.1 |
10.1 |
三渔冲村 |
7.5 |
6.5 |
7 |
5.8 |
2.9 |
3.5 |
4.5 |
4.4 |
9.8 |
晒鼓坪村 |
6.4 |
5.3 |
5.6 |
4.4 |
2.2 |
3.3 |
4.5 |
3.5 |
8.8 |
厚丰溪村 |
10.6 |
9.5 |
9.9 |
8.6 |
6.3 |
7.1 |
8.1 |
7.8 |
13.1 |
两河口村 |
10.3 |
9.2 |
9.6 |
8.4 |
5.8 |
6.5 |
7.3 |
7.3 |
12.7 |
王家棚村 |
14.2 |
13.2 |
13.7 |
12.5 |
9.7 |
10 |
10.5 |
11.1 |
16.5 |
土地坡村 |
16.5 |
15.5 |
16 |
14.7 |
11.9 |
12.3 |
12.9 |
13.4 |
18.8 |
全伏山村 |
19.7 |
18.7 |
19.1 |
17.9 |
15.1 |
15.5 |
16 |
16.6 |
22 |
郑家榜村 |
20.7 |
19.7 |
20.1 |
18.8 |
16.1 |
16.5 |
17.1 |
17.6 |
23 |
观坪林场 |
23.7 |
22.6 |
23 |
21.7 |
19.2 |
19.7 |
20.4 |
20.7 |
26.1 |
3. 模型分析与计算
求解最小支撑树的算法有很多种,如Prim算法与Kruskal算法[6],考虑到上述模型为0-1整数规划模型[6] [7],应用破圈法求图 G 的最小支撑树[7],在LINGO环境下编程得到实现[4] [7],计算结果如下图2,总长度54.84 km。
Figure 2. Calculation results of the minimum support tree theory model
图2. 最小支撑树理论模型计算结果
4. 计算结果的可行性分析
由于清江流域的存在,该方案的多处连接跨过了水域:8号何家坪社区和5号枫竹园社区;11号西寺坪社区和5号枫竹园社区;8号何家坪社区和10号后山社区;1号湖口社区和9号板桥铺社区;11号西寺坪社区和20号晒鼓坪村;15号黄家坪村和14号刘家坳村;6号津洋口社区和7号丹水社区;24号土地坡村和25号全伏山村;19号三渔冲村和22号两河口村;26号郑家榜村和27号观坪林场。为了规划布局的可行性,可以将龙舟坪镇各村分为7组,7组分别为:
第一组:1,2,3,4,7,7,10
第二组:5,6,15,16,17,18,19,20
第三组:9,12,13,14
第四组:8,11
第五组:21,22,23,24
第六组:25,26
第七组:27
地下燃气管网的铺设不仅需要避开清江流域,而且尽可能地考虑到每个天然气供应站的负载平衡,以及覆盖村子数相对集中与均衡,可修正规划方案如下图3,总长40.64 km,此方案需要多建设5个天然气站点:第一组与第三组合计1站;第二组1站;第四组1站;第五组1站;第六组与第七组合计1站。虽然增加了建设成本,但是方案的可行性增强了,且能更好保障管道铺设过程中的施工人员安全。图4是应用GIS可视化技术在天地图中展示的规划布局方案[8]。
Figure 3. Revised gas pipeline planning scheme
图3. 修正后的燃气管道规划方案
注:该图基于自然资源部标准底图服务网站下载的审图号为GS(2016)008号的标准地图制作,底图无修改。
Figure 4. The revised planning scheme is visualized in 2D GIS on the sky map
图4. 修正后的规划方案在天地图中GIS二维可视化
5. 结论
应用图论最小支撑树理论模型与0-1整数规划的破圈算法,结合天地图GIS测量技术得到的龙舟坪镇各村(社区)地下燃气管网GIS规划布局初步方案。可以考虑以下几个因素对模型继续改进:
① 考虑地形因素:龙舟坪镇的地形复杂,包括山地、丘陵、河流等。不同的地形对管网的施工难度和成本有很大影响。可以在模型中加入地形系数,以反映不同地形的施工难度和成本差异。例如,在丘陵施工难度更大,成本更高,地形系数可以设置为较高值;而在相对平缓的地区施工难度较小,成本较低,地形系数可以设置为较低值。
② 考虑人口密度因素[9]:龙舟坪镇不同区域的人口密度差异较大。人口密度对管网的需求和成本有影响。可以在模型中加入人口密度系数,以反映不同区域的管网需求和成本差异。例如,人口密度高的区域,对管网的需求更大,成本也更高,湖口社区,四冲湾社区,龙门社区,花坪社区,后山社区,这5个社区离得很近,导致该地区人口密度大,且建筑物多,公路复杂,人口密度系数可以设置为较高值;而人口密度低的区域对管网的需求较小,成本也较低,白氏坪社区,每户之间距离较远,人口密度系数可以设置为较低值。
③ 考虑环境保护因素[10]:在乡村地区环境保护至关重要,管网建设可能对环境产生影响,包括对植被、水源和土壤的破坏。可以在模型中加入环境保护系数,以反映不同区域的环保要求和成本差异。例如在生态敏感区域,环保要求更高,成本也更高,环境保护系数可以设置为较高值,观坪林场、郑家榜村就需要保护其森林资源,板桥铺村、枫竹园社区距离清江较近,在铺设时需要保护水资源;在普通区域环保要求较低,成本也较低,环境保护系数可以设置为较低值。
总之,按照地下燃气管网行业施工标准,乡村社区间管道距离测量数据肯定还与实际距离存在误差,但只需要求解最短布线路径,具体的管道规划布局方案还需考虑上述多因素考,本研究具有重要参考价值,不影响总体规划布局结果。
致 谢
长阳土家族自治县城乡建设规划局。
NOTES
*通讯作者。