岩土力学 ›› 2022, Vol. 43 ›› Issue (3): 843-856.doi: 10.16285/j.rsm.2021.0872

• 数值分析 • 上一篇    

基于变步长加速搜索的微震源定位方法

贾宝新1, 2,李峰1,潘一山3,周琳力1   

  1. 1. 辽宁工程技术大学 土木工程学院,辽宁 阜新 123000; 2. 辽宁工程技术大学 辽宁省矿山沉陷灾害防治重点实验室,辽宁 阜新 123000;3. 辽宁大学 环境学院,辽宁 沈阳 110036
  • 收稿日期:2021-06-10 修回日期:2021-12-30 出版日期:2022-03-22 发布日期:2022-03-23
  • 作者简介:贾宝新,男,1978年生,博士,教授,主要从事矿山灾害力学等方面的研究。
  • 基金资助:
    国家自然科学基金(No.51774173);辽宁工程技术大学学科创新团队资助项目(No.LNTU20TD08);辽宁省“兴辽英才计划”项目(No.XLYC2007163);辽宁“百千万人才工程”培养经费。

Microseismic source locating method based on variable step size accelerated search

JIA Bao-xin1, 2, LI Feng1, PAN Yi-shan3, ZHOU Lin-li1   

  1. 1. School of Civil Engineering, Liaoning Technical University, Fuxin, Liaoning 123000, China; 2. Liaoning Key Laboratory of Mine Subsidence Disaster Prevention and Control, Liaoning Technical University, Fuxin, Liaoning 123000, China; 3. School of Environment, Liaoning University, Shenyang, Liaoning 110036, China
  • Received:2021-06-10 Revised:2021-12-30 Online:2022-03-22 Published:2022-03-23
  • Supported by:
    This work was supported by the National Natural Science Foundation of China(51774173), the Discipline Innovation Team of Liaoning Technical University (LNTU20TD08), the “Rejuvenating Liaoning Talents Plan” Project of Liaoning Province(XLYC2007163) and the Liaoning BaiQianWan Talents Program.

摘要: 微震定位方法是微震监测技术的重要组成部分,其关键是定位震源位置。利用空间网格划分并计算网格交点目标函数值,对微震定位目标函数二维及三维空间分布进行了分析,并据此获取了目标函数连续且极小值唯一、单轴收敛范围逐步减小、各轴收敛范围不一的规律。利用以上规律及模式搜索法、网格搜索法的优缺点,探索出了基于连续比较模块、变步长模块、加速模块的变步长加速搜索法。通过模拟算例与工程数据下收敛稳定性、结果精确度、计算速度以及参数初始值影响程度4个指标的效果对比,结果表明:模拟算例下,对比模拟退火算法、遗传算法,变步长加速搜索法的目标函数值标准差、定位误差标准差、波速误差标准差均为0;该算法的定位误差平均值分别为其余二者的0.7%、1.9%;该算法的计算时间平均值分别为其余二者的6.9%、33.2%。该算法单独更改各参数对定位误差的影响在0.005~0.025 m之间;减小搜索步长下限可有效提高结果精确度,并增加相应的计算时间。在规定初至到时与目标函数模型及检波器位置坐标下,搜索算法对定位精度无实质影响。

关键词: 微震定位, 空间图像, 变步长加速搜索法, 模拟退火算法, 遗传算法

Abstract: Microseismic locating method is an important part of microseismic monitoring technology, the key of which is to locate the hypocenter. We analyze the two-dimensional and three-dimensional spatial distributions of microseismic locating objective functions by using spatial gridding and calculating the objective function values of grid intersections. Accordingly, we find that the objective function is continuous and the minimum value is unique, the convergence range of single axis decreases gradually, and the convergence range of each axis varies. Using the above findings and the advantages and disadvantages of pattern search method and grid search method, we propose the variable step size accelerated search method based on continuous comparison module, the variable step size module and the acceleration module. Through the comparison of four indexes between the simulation example and the engineering data: the convergence stability, the accuracy of the results, the speed of calculation and the degree of influence of initial values of parameters, we get the results that: in the simulation example, compared with that of simulated annealing algorithm and genetic algorithm, standard deviations of objective function value, locating error and wave velocity error of the variable step size accelerated search method are all 0; the average locating error of the variable step size accelerated search method is 0.7% and 1.9% of that of other two algorithm, respectively; the average calculation time of this algorithm is 6.9% and 33.2% of that of the other two algorithm, respectively. The influence of the algorithm changing each parameter individually on the locating error is between 0.005 m and 0.025 m; reducing the lower limit of the step size for search can effectively improve the accuracy of the results but increase calculation time. When the initial arrival time, objective function model and coordinates of the geophone position are specified, the search algorithm has no substantial impact on the locating accuracy.

Key words: microseismic locating, spatial image, variable step size accelerated search method, simulated annealing algorithm, genetic algorithm

中图分类号: TD315
[1] 王新龙, 聂利青, 蔡国军, 张宁, 赵泽宁, 刘薛宁, 宋登辉, . 基于孔压静力触探技术的SVR优化算法评估土体液性指数[J]. 岩土力学, 2024, 45(S1): 645-653.
[2] 肖晓春, 丁振, 丁鑫, 徐军, 樊玉峰, 李子阳, . 基于P波到时拾取的分析法与遗传算法联合定位方法与应用[J]. 岩土力学, 2024, 45(7): 2195-2207.
[3] 李林, 孙砖芹. 考虑基质吸力影响的非饱和路堤三维稳定性上限分析[J]. 岩土力学, 2024, 45(4): 1014-1025.
[4] 董福瑞, 王述红, 侯钦宽. 基于主成分分析的多参数岩体结构面 优势分组方法研究[J]. 岩土力学, 2022, 43(9): 2457-2470.
[5] 王伟, 陈国庆, 郑水全, 张广泽, 王栋, . 考虑张拉-剪切渐进破坏的边坡矢量和法研究[J]. 岩土力学, 2019, 40(S1): 468-476.
[6] 崔学杰, 晏鄂川, 陈 武. 基于改进遗传算法的岩体结构面产状聚类分析[J]. 岩土力学, 2019, 40(S1): 374-380.
[7] 蒋若辰, 徐奴文, 戴峰, 周家文. 基于快速行进迎风线性插值的微震定位算法研究[J]. 岩土力学, 2019, 40(9): 3697-3708.
[8] 马春辉, 杨杰, 程琳, 李婷, 李雅琦, . 基于量子遗传算法与多输出混合核相关向量机的堆石坝材料参数自适应反演研究[J]. 岩土力学, 2019, 40(6): 2397-2406.
[9] 郭 超, 高永涛, 吴顺川, 成子桥, 张诗淮, 韩龙强, . 基于三维快速扫描算法与到时差数据库技术的层状介质震源定位方法研究[J]. 岩土力学, 2019, 40(3): 1229-1238.
[10] 温树杰,梁 超,宋亮亮,刘 刚,. 基于最小势能法的三维临界滑裂面搜索方法[J]. , 2018, 39(7): 2708-2714.
[11] 王开禾,罗先启,沈 辉,张海涛. 围岩力学参数反演的GSA-BP神经网络模型及应用[J]. , 2016, 37(S1): 631-638.
[12] 王 伟 ,杨 敏 ,上官士青,. 控制差异沉降的桩筏基础桩径优化分析方法[J]. , 2015, 36(S2): 178-184.
[13] 李 忠 ,杨 俊 , . 基于MPGA的复杂应力状态边坡稳定性分析[J]. , 2015, 36(5): 1488-1495.
[14] 李 宁 ,王李管 ,贾明涛 ,陈建宏 ,谭正华,. 改进遗传算法和支持向量机的岩体结构面聚类分析[J]. , 2014, 35(S2): 405-411.
[15] 邱道宏 ,李术才 ,薛翊国 ,田 昊 ,闫茂旺,. 基于数字钻进技术和量子遗传-径向基函数神经网络的围岩类别超前识别技术研究[J]. , 2014, 35(7): 2013-2018.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!