一种基于种群进化算法电阻率溶洞识别方法技术

本发明专利技术公开了一种基于种群进化算法的电阻率溶洞识别方法，该方法主要步骤如下：S1.构建溶洞有限元模型，利用电阻抗法正问题相关公式求得结构测量点的电势分布数据；S2.构建溶洞结构的目标函数，即待优化的目标函数；S3.利用种群进化算法不断的迭代优化目标函数，满足停止准则后，保存最优解，便可识别出溶洞位置。该方法对比传统的灵敏度或梯度矩阵方法可以准确地检测出各种溶洞的位置，无需初值和梯度信息，并且对噪声不敏感。本发明专利技术通过单个溶洞和多个溶洞两个实施例说明了该方法的有效性和鲁棒性，与粒子群算法相比，本发明专利技术即使在一定的测量噪声下也具有较好的探测精度，具有良好的工程应用能力。

A method of resistivity cave recognition based on population evolution algorithm

【技术实现步骤摘要】
一种基于种群进化算法电阻率溶洞识别方法
本专利技术涉及一种基于种群进化算法电阻率溶洞识别方法，属于地下溶洞探测

技术介绍
目前，在石灰岩地下溶洞发育的地区，建筑工程场地的地下溶洞，往往位于地下水位以下，建设各方为了保证建筑物基础的安全，不惜投入大量的人力、财力及时间进行石灰岩地区建设工程场地地下溶洞的勘查工作，采用传统的工程地质钻探仅为一孔之见，或采用传统的电法勘探，由于地下溶洞埋深较大或地下溶洞规模较小而导致其电阻率与围岩石灰岩电阻率差异较小而无法准确地探查出来，给建筑物留下严重的安全隐患。地下溶洞会对房屋建筑、地铁等建造产生一定危害，及时探测出地下溶洞并加以处理，对结构施工、使用等中的安全有重要意义。电阻率法探测溶洞问题是反问题的识别研究中的优化问题，其基本思想是：确定探测区域，电极方案，溶洞的出现会引起测量的电势分布数据的改变，进而反应到土层视电阻率的改变，可以利用这些变化对溶洞位置进行定位。即通过定义一个关于溶洞结构的目标函数，然后利用各种优化方法来实现探测溶洞的位置。传统的电阻率溶洞探测方法引入了正则化方法，采用灵敏度或梯度矩阵来识别溶洞位置，具有以下缺点：对初值敏感，对噪声敏感，容易陷入局部最优解，从而识别效果较差。
技术实现思路
本专利技术的目的在于提供一种具有实用性、有效性和准确性的跨孔电阻率溶洞识别方法，该方法可以准确地检测出各种溶洞的位置，并且对噪声不敏感，具有良好的工程应用能力。本专利技术的技术方案：一种跨孔电阻率溶洞识别方法，包括以下步骤：步骤一：建立溶洞土层有限元模型，确定电极方案，利用电阻抗法得到测量点的电势分布数据；步骤二：构建溶洞结构的目标函数，即优化的目标函数，该目标函数如下：c＝[c1；c2；...；cm](1)其中，g(c)为目标函数，||(.)||表示向量的2范数，即c是各单元视电阻率，Su表示测量电势点集合上的数据，是第i组电极方案测量的电势分布数据，是第i组电极方案识别出来的电势分布数据，是s组测量数据组成的矩阵，R(c)是s组识别数据组成的矩阵；步骤三：利用种群进化算法不断优化目标函数，满足停止准则后，最终可以得到溶洞的识别位置。上述方法的步骤三中利用种群进化算法对目标函数优化的具体过程如下：S1，种群参数和初始化：控制参数包括迭代最大次数Itermax、种群大小数Np、决策变量数n、决策变量的上界Xmax和下界Xmin：Xp，q＝Xmin+rand·(Xmax-Xmin)，p＝1，2，...，Np，q＝1，2，...，n(1)其中，rand是在[0,1]范围内的均匀分布的随机数，候选解的适应值f＝(f1，f2，...，fNp)根据以下目标函数计算：fp＝obj(Xp，1，Xp，2，…，Xp，n)，p＝1，2，…，Np(2)S2，新位置生成阶段：在这个阶段，种群候选解被认为会围绕自已的初始位置移动到一个更好的位置，对于候选解p，从种群(p≠1)中随机选择另一个候选解1，用于生成潜在的搜索方向，根据候选解p新旧位置之间的距离，通过以下两个模型更新位置：(i)模型1：候选解p的新位置直接向候选解1学习，探索半径随迭代动态变化：(ii)模型2：在原位置附近产生候选解p的新位置：S3，空间搜索增强阶段：为了进一步加强对每个维度的深入搜索，每次迭代都采取以下步骤：(i)在总体中找出最佳和最差的候选解；(ii)通过改变一个维度的值，同时保持其他维度的值；(iii)将新生成的解与原解的适应度值进行比较，保留较好的解；(iv)对其他方面分别重复步骤(ii)和(iii)，新生成的解决方案由以下公式产生：S4，停止准则阶段：保存当前最优解，如果满足最大迭代次数，算法将终止；否则，将重复步骤S2和S3。与现有技术相比，本专利技术技术方案的有益效果是：本专利技术通过溶洞结构测量点的电势分布数据构建目标函数，利用元启发式算法来识别溶洞位置，对初值不敏感，对噪声不敏感，且不容易陷入局部最优解，具有更好的效率和精度。附图说明图1为溶洞识别问题归化为优化问题的流程图；图2为种群进化算法的实现流程示意图；图3为本专利技术实施例1中单个溶洞的有限元模型；图4为本专利技术实施例2中多个溶洞的有限元模型；图5为本专利技术实施例1中采用本专利技术所述方法在无噪音情况下的检测结果；图6为本专利技术实施例1中采用本专利技术所述方法在有0.1％噪音情况下的检测结果；图7为本专利技术实施例1中采用粒子群算法在无噪音情况下的检测结果；图8为本专利技术实施例1中采用粒子群算法在有0.1％噪音情况下的检测结果；图9为本专利技术实施例1中的有0.1％噪音情况下单个溶洞的算法收敛曲线对比图。图10为本专利技术实施例2中采用本专利技术所述方法在无噪音情况下的检测结果；图11为本专利技术实施例2中采用本专利技术所述方法在有噪音情况下的检测结果；图12为本专利技术实施例2中采用粒子群算法在无噪音情况下的检测结果；图13为本专利技术实施例2中采用粒子群算法在有0.1％噪音情况下的检测结果；图14为本专利技术实施例2中有0.1％噪音情况下多个溶洞的算法收敛曲线对比图。具体实施方式为了使本专利技术的目的、技术方案和优点更加清楚，下面结合附图和实施例对本专利技术作进一步的详细说明。参见图1，本专利技术的一种基于种群进化算法的电阻率溶洞识别方法，具体过程分为两步：(1)目标函数溶洞模型简化为矩形探测区域。电流点与电势点不一致，但数目相等。使用有限元进行正问题求解，有如下方程：Ku＝f其中，u＝[u1；u2；…；un]包含了所有测量节点电势，f＝[f1；f2；…；fn]则表征电流向量，不同的电极方案，f不同。K为刚度阵，与各单元的视电阻率c＝[c1；c2；…；cm]线性相关其中，Kj为单位电阻率时的单元刚度矩阵。当土层中没有溶洞时，所有单元的视电阻率取c0＝10，当出现溶洞时，所在单元的视电阻率会随之降低，本专利技术中出现溶洞的单元视电阻率取3，即cj＝3(j＝1,2,3,...,m)即代表第j单元出现溶洞。通过上述有限元的公式，c与u存在函数关系，我们可以通过识别出c推出计算的电势分布数据基于电势分布数据构建的目标函数如下：c＝[c1；c2；...；cm](2)其中，c是各单元视电阻率，Su表示测量电势点集合上的数据，是第i组电极方案测量的电势分布数据，是第i组电极方案识别出来的电势分布数据，是s组测量数据组成的矩阵，R(c)是s组计算数据组成的矩阵。当测量的电势分布数据与识别的电势分布数据完全吻合时，数值最小，溶洞识别问题等价为优化问题，种群的某个个体位置即为一种视电阻率分布情况，当目标函数达到极小值时，通过识别出来的视电阻率c，即最优解，便能反映出溶洞的位本文档来自技高网...

【技术保护点】
1.一种基于种群进化算法的电阻率溶洞识别方法，其特征在于包括以下步骤：/n步骤一：建立溶洞土层有限元模型，确定电极方案，利用电阻抗法得到测量点的电势分布数据；/n步骤二：构建溶洞结构的目标函数，即优化的目标函数，该目标函数如下：/n

【技术特征摘要】
1.一种基于种群进化算法的电阻率溶洞识别方法，其特征在于包括以下步骤：
步骤一：建立溶洞土层有限元模型，确定电极方案，利用电阻抗法得到测量点的电势分布数据；
步骤二：构建溶洞结构的目标函数，即优化的目标函数，该目标函数如下：



其中，g(c)为目标函数，||(.)||表示向量的2范数，即c是各单元视电阻率，Su表示测量电势点集合上的数据，是第i组电极方案测量的电势分布数据，是第i组电极方案识别出来的电势分布数据，是s组测量数据组成的矩阵，R(c)是s组识别数据组成的矩阵；
步骤三：利用种群进化算法不断优化目标函数，满足停止准则后，最终可以得到溶洞的识别位置。


2.根据权利要求1所述的基于种群进化算法的电阻率溶洞识别方法，其特征在于：步骤三中利用种群进化算法的对目标函数优化的具体过程如下：
S1，种群参数和初始化：控制参数包括迭代最大次数Itermax、种群大小数Np、决策变量数n、决策变量的上界Xmax和下界Xmin：
Xp，q＝Xmin+rand·(Xmax-Xmin)，p＝1，2，...，Np，q＝1，2，...，n(1)
其中，rand是在[0,1]范围内的均匀分布...

【专利技术属性】
技术研发人员：李晓阳，孙瑞甲，卢春燕，罗盈洲，
申请(专利权)人：中国建筑第四工程局有限公司，中建四局第一建筑工程有限公司，
类型：发明
国别省市：广东;44