一种CO2地质封存区域的传感网络节点定位优化方法与流程

本发明涉及无线传感网络节点定位方法，尤其涉及一种co2地质封存区域的传感网络节点定位优化方法。

背景技术：

全球变暖问题已经引起国际的密切关注，二氧化碳作为最重要的温室气体亦是全球变暖的主要原因。目前，中国现已成为第二大co2排放国，并且将长期依赖矿物燃料特别是煤、石油和天然气等作为主要能源，是co2潜在的第一大排放国。虽然，大气中的co2可以通过陆地生态系统中的植被、微生物和土壤，以及海洋生态系统中的浮游生物吸收，但生物和其它固碳方式所产生的效果十分有限，因此还需要通过地质固碳等其它途径缓减co2气体排放的增长速度。一种较为直接、有效的方法是将固定点源产生的co2，捕集后长期储存于相对封闭的地质构造或海洋中，从而阻止或显著减少co2向大气中的人为排放，但是地质封存后的co2必须依靠地层的压力维持超临界流体态，一旦通过断层、裂隙或油气井等通道泄漏到地表，会在近地表形成比空气重的气云，对人类健康、生命安全及生态环境产生严重影响，使得对抗全球变暖的努力付之东流。为了确保co2能够长期安全封存在地下，就需要完善的监测技术提供保障，精准的监测设备能够提高对封存过程的认识并证实其有效性，验证co2地质封存持久性和安全性的关键是连续自动监测。传统的网络监测覆盖控制方案往往是在监测区域中或随机、或均匀地布设离散的传感节点，然而并没有从实际应用意义上解决网络覆盖的节点优化布设问题，造成了检测成本的增加和浪费。

与本发明最为相近的已有技术有：

专利号为200810201237.0公开了一种面向区域监测的无线传感器网络节点布设的优化方法，该方法利用粒子群优化算法和匈牙利算法从目标检测功能的角度来提高监测区域的覆盖率，高效地解决了无线传感器网络节点布设优化所面临的高维优化问题。然而上述无线传感器网络节点布设的优化方法面对检测区域复杂的地理、地质环境具有一定的局限性，且在该优化方法中，通过空投方式将传感器布撒在检测区域来确定初始位置，难以确定空投位置的科学性及准确性。

专利号为201310118083.x公开了一种无线传感器网络节点覆盖优化方法，该方法采用微分进化算法进行迭代改进，同时采用约束圆心范围、减少多个圆重叠的方法对群体进一步优化，使用尽量少的传感器完成尽可能高的区域覆盖率。但是针对二氧化碳地质封存监测传感器覆盖网而言，由于封存区域地质地理环境多元化，传感器覆盖网优化更多考虑的是检测区域内不同监测点的不同地理地质特征、监测点传感器网络节点密度以及可能出现的不同程度的干扰因素，该无线传感器网络节点覆盖优化方法在本方面具有不可实现性。

技术实现要素：

发明目的：为了克服现有技术中存在的不足，本发明提供一种co2地质封存区域的传感网络节点定位优化方法，具有较强的网络覆盖能力和网络连通性，同时可以减少节点冗余度和通信开销。

技术方案：为实现上述目的，本发明采用的技术方案为：

一种co2地质封存区域的传感网络节点定位优化方法，包括如下步骤：

步骤1)，通过分析co2地质封存区域的地质、地理和气象数据，获得co2地质封存区域的co2泄漏事件的影响因子集合并确定权重，然后运用gis空间分析技术得出co2地质封存区域敏感程度分布图；

步骤2)，根据co2地质封存区域的不同敏感等级，采用不同密度的网络覆盖控制算法布设传感监测节点；

步骤3)，对co2地质封存区域中部署的传感节点集合进行delaunay剖分，完成对覆盖网络的描述和优化表达。

进一步的，所述步骤1)中，确定权重的具体步骤为：根据提取的封存区域co2泄露事件的影响因子集合，构建co2地质封存区域环境敏感影响因子评价指标体系；所述评价指标体系由目标层、准则层和判别层组成，其中一级评估指标为监测环境敏感度a；一级评估指标包括的二级评估指标有：地质储层b1、地形地貌b2、气象风场b3；二级评估指标包括的三级评估指标有：埋藏深度c1、断裂活动c2、储层渗透度c3、储层孔隙度c4、地热条件c5、坡度c6、坡向c7、矿井位置c8、土地利用c9、地表覆盖c10、土壤类型c11、盛行风力c12、盛行风向c13；依据评价指标体系的影响因子层次结构，采用层析分析法构建判断矩阵进行计算比较，经层层迭代计算后得到最后一级指标层对一级指标层的总排序权重。

进一步的，所述步骤1)中，运用gis空间分析技术得出co2地质封存区域敏感程度分布图，具体包括如下步骤：

a，分析地质调查和矿井专题地图数据，得到co2地质封存区域煤储层及其围岩的孔隙度、渗透度、地热条件、埋藏深度和断裂活动影响因子的值，利用gis技术将各单项影响因子通过加权叠加进行空间叠加，获得co2地质封存区域的co2潜在泄漏通道，并对其进行缓冲区分析得出co2泄漏的影响程度和空间分布状况图层；

b，基于co2地质封存区域已获取的数字高程模型数据，确认co2地质封存区域的地理位置及范围，并采集co2地质封存区域的自然地理特征；运用获取的数字高程模型数据进行坡度和坡向分析，得出地形坡度图以及地形坡向图；并对地形坡度图和地形坡向图分别进行栅格重分类，得到坡度分级图和坡向分类图；其中分类标准为：在地形坡度图中将小于15°定义为缓坡，大于等于15°定义为陡坡；在地形坡向图中将坡向分为向南、向北、向东、向西、向东南、向西南、向东北和向西北8个方向；

基于遥感数据提取co2地质封存区域的土地利用现状、土壤资源类型和植被覆盖现状，并现场确认co2地质封存区域、城镇居民点以及其他人为co2排放源的相对位置，对其进行距离分析，并将分析结果按照四个距离层次进行栅格重分类，得出人为co2排放源影响范围图层；

c，提取co2地质封存区域内的最小、平均和最大风速，对co2地质封存区域内主导风向进行编码，结合地形坡度图以及地形坡向图生成迎风坡、背风坡和顺风坡的风场空间分布图层；然后以矿井位置为原点，以主导风向为轴线，在矿井下风向画扇形区域作为泄露扩散蓄积区分布，所述扇形区域应包括整个co2地质封存区域；

d，运用arcgis地理信息处理软件，将所述co2泄漏的影响程度和空间分布状况图层、所述人为co2排放源影响范围图层、所述迎风坡、背风坡和顺风坡的风场空间分布图层运用gis技术按照所述评价指标体系所得权重进行加权叠加，分别计算出每个评价单元的泄露监测敏感度综合指数，得出co2地质封存区域敏感程度分布图，并对其进行五个等级重分类，得出co2地质封存区域敏感等级空间分布图，分别为高敏感、较高敏感、一般敏感、较低敏感和低敏感。

进一步的，所述步骤2)具体为：假设每个传感节点均实行全方位监测，将其覆盖范围作为感知半径为r的圆形区域，且设各传感节点都具有相同的发射功率，即所有的传感节点的检测半径r均相等；以传感节点间距离a间接表示覆盖密度，根据co2地质封存区域敏感等级不同，以传感节点为中心，以传感节点间距离a作为格网边长，六个方向新增传感节点对co2地质封存区域进行规则三角网剖分，并满足条件：若co2地质封存区域敏感等级越高，则传感节点间距离a越小。

进一步的，所述步骤3)包括如下具体步骤：

a，以co2注入井作为簇头节点构建初始delaunay三角网剖分，求解簇头节点的初始voronoi域；

b，将簇头节点作为初始生长点，根据初始生长点所处区域的监测敏感度级别选择不同的传感节点距离a作为格网布设，以正三角形作为格网划分单元，以初始生长点为中心从六个方向计算拟增传感节点点集的空间位置，依次判断拟增传感节点是否在voronoi中，将落在voronoi域中的传感节点添加进生长点集中，并将新增传感节点作为新的生长点再次进行敏感度判断，直至voronoi域所有生长节点的拟增传感节点点集都在voronoi域外。

c，结合簇头节点和感知节点，求解delaunay监测网络优化覆盖控制方案。有益效果：本发明方法根据co2地质封存区域的地质、地理和气象等数据，分析出泄漏事件的影响因子集合，采用可变密度节点布设优化不同封存区域内的无线网络传感节点部署方案，具有较强的网络覆盖能力和网络连通性，同时可以减少节点冗余度和通信开销。能够降低布设成本、提高监测质量和延长生命时间，提高泄漏监测预警的实时性、预见性和有效性。

附图说明

图1是本发明一种co2地质封存区域的传感网络节点定位优化方法的流程图。

图2是本发明确定co2地质封存区域泄露敏感度分区技术路线图。

图3是本发明确定co2地质封存区域泄露敏感度分区图。

图4是本发明可变密度传感节点部署模型高密度节点布设效果图。

图5是本发明可变密度传感节点部署模型较高密度节点布设效果图。

图6是本发明可变密度传感节点部署模型一般密度节点布设效果图。

图7是本发明可变密度传感节点部署模型较低密度节点布设效果图。

图8是本发明可变密度传感节点部署模型低密度节点布设效果图。

图9是本发明可变密度的网络优化覆盖算法流程图。

图10是本发明仿真区域煤层co2封存区域监测网络优化部署场景效果图。

注：图中●表示传感节点布设位置。

具体实施方式

下面结合附图对本发明作更进一步的说明。

本发明以沁水盆地煤层co2注入区5000m*4000m作为监测仿真区域，以100m*100m分辨率栅格进行监测敏感度分析，以14个ch4开采井作为簇头节点进行常规监测，传感节点的感知半径r＝100m，按照监测敏感程度优化部署算法布设监测场景。

如图1所示，本发明一种co2地质封存区域的传感网络节点定位优化方法主要流程图，首先，通过分析co2地质封存区域的地质、地理和气象数据，获得co2地质封存区域的co2泄漏事件的影响因子集合并确定权重，然后运用gis空间分析技术得出co2地质封存区域敏感程度分布图。其次，设计不同密度的监测传感网络节点覆盖控制方案，根据封存区域的不同敏感等级采用不同密度的网络覆盖控制算法或密集或稀疏地布设传感监测节点。最后，对co2地质封存区域中部署的传感节点集合进行delaunay剖分，完成对覆盖网络的描述和优化表达。

如表1所示，确定权重的具体步骤为：根据提取的封存区域co2泄露事件的影响因子集合，构建co2地质封存区域环境敏感影响因子评价指标体系。评价指标体系由目标层、准则层和判别层组成，一级评估指标为监测环境敏感度a；一级评估指标包括的二级评估指标有：地质储层b1、地形地貌b2、气象风场b3；二级评估指标包括的三级评估指标有：埋藏深度c1、断裂活动c2、储层渗透度c3、储层孔隙度c4、地热条件c5、坡度c6、坡向c7、矿井位置c8、土地利用c9、地表覆盖c10、土壤类型c11、盛行风力c12、盛行风向c13。

表1co2地质封存区域泄露监测敏感度影响因子体系表

各影响因子对封存区域敏感程度的影响并不相同，依据评价指标体系影响因子层次结构，采用层次分析法构建判断矩阵进行计算比较，经层层迭代计算后得到最后一级指标层对一级指标层的总排序权重，如表2所示。

表2三级指标层(c)的总排序权重

如图2所示，运用gis空间分析技术得出co2地质封存区域敏感程度分布图，具体包括如下步骤：

a，分析地质调查和矿井专题地图数据，得到co2地质封存区域煤储层及其围岩的孔隙度、渗透度、地热条件、埋藏深度和断裂活动影响因子的值，利用gis技术将各单项影响因子通过加权叠加进行空间叠加，获得co2地质封存区域的co2潜在泄漏通道，并对其进行缓冲区分析得出co2泄漏的影响程度和空间分布状况图层。

b，基于co2地质封存区域已获取的数字高程模型数据，确认co2地质封存区域的地理位置及范围，并采集co2地质封存区域的自然地理特征。运用获取的数字高程模型数据进行坡度和坡向分析，得出地形坡度图以及地形坡向图；并对地形坡度图和地形坡向图分别进行栅格重分类，得到坡度分级图和坡向分类图；其中分类标准为：在地形坡度图中将小于15°定义为缓坡，大于等于15°定义为陡坡；在地形坡向图中将坡向分为向南、向北、向东、向西、向东南、向西南、向东北和向西北8个方向。

基于遥感数据提取co2地质封存区域的土地利用现状、土壤资源类型和植被覆盖现状，并现场确认co2地质封存区域、城镇居民点以及其他人为co2排放源的相对位置，对其进行距离分析，并将分析结果按照四个距离层次进行栅格重分类，得出人为co2排放源影响范围图层。

c，提取co2地质封存区域内的最小、平均和最大风速，对co2地质封存区域内主导风向进行编码，结合地形坡度图以及地形坡向图生成迎风坡、背风坡和顺风坡的风场空间分布图层。然后以矿井位置为原点，以主导风向为轴线，在矿井下风向画扇形区域作为泄露扩散蓄积区分布，所述扇形区域应包括整个co2地质封存区域。

d，运用arcgis10.2地理信息处理软件，将co2泄漏的影响程度和空间分布状况图层、人为co2排放源影响范围图层、迎风坡、背风坡和顺风坡的风场空间分布图层运用gis技术按照所述评价指标体系所得权重进行加权叠加，分别计算出每个评价单元的泄露监测敏感度综合指数，监测区域中一个栅格为一个评价单元，得出co2封存区域敏感程度分布图，并对其进行五个等级重分类，得出co2封存区域敏感等级空间分布图，分别为高敏感、较高敏感、一般敏感、较低敏感和低敏感，定性地给出煤层co2注入区泄露监测敏感度定量分区方案，如图3所示。

根据封存区域敏感程度不同，设计不同的网络覆盖控制算法或密集、或稀疏地布设传感节点，co2地质封存区域敏感等级越高，传感节点间距离a越小。假设每个传感节点均实行全方位监测，将其覆盖范围作为感知半径为r的圆形区域，且设各传感节点都具有相同的发射功率，即所有的传感节点的检测半径r均相等；以传感节点间距离a间接表示覆盖密度，根据co2封存区域敏感程度不同，以传感节点为中心，以传感节点间距离a作为格网边长，六个方向新增传感节点对co2封存区域进行规则三角网剖分，并满足条件：co2封存区域敏感等级越高，传感节点间距离a越小。即对于高敏感区域采用高密度节点布设(a＝r)，如图4所示；对于较高敏感区域采用较高密度节点布设如图5所示；对于中敏感区域其采用一般密度节点布设如图6所示；对于较低敏感区域和低敏感区域，分别采用较低密度节点布设(a＝2r)和低密度节点布设(a＝4r)，如图7、8所示，如需了解未覆盖区域的co2浓度，可通过空间插值法得出。

如图9所示，以co2注入井作为簇头节点构建初始delaunay三角网剖分，求解簇头节点的初始voronoi域；将簇头节点作为初始生长点，根据初始生长点所处区域的监测敏感度级别选择不同的传感节点距离a作为格网布设，以正三角形作为格网划分单元，以初始生长点为中心从六个方向计算拟增传感节点点集的空间位置，依次判断拟增传感节点是否在voronoi中，将落在voronoi域中的传感节点添加进生长点集中，并将新增传感节点作为新的生长点再次进行敏感度判断，直至voronoi域所有生长节点的拟增传感节点点集都在voronoi域外。结合簇头节点和感知节点，求解delaunay监测网络优化覆盖控制方案，结果如图10所示。

以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。