岩土力学 ›› 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] 崔学杰, 晏鄂川, 陈 武. 基于改进遗传算法的岩体结构面产状聚类分析[J]. 岩土力学, 2019, 40(S1): 374-380.
[2] 王伟, 陈国庆, 郑水全, 张广泽, 王栋, . 考虑张拉-剪切渐进破坏的边坡矢量和法研究[J]. 岩土力学, 2019, 40(S1): 468-476.
[3] 蒋若辰, 徐奴文, 戴峰, 周家文. 基于快速行进迎风线性插值的微震定位算法研究[J]. 岩土力学, 2019, 40(9): 3697-3708.
[4] 马春辉, 杨杰, 程琳, 李婷, 李雅琦, . 基于量子遗传算法与多输出混合核相关向量机的堆石坝材料参数自适应反演研究[J]. 岩土力学, 2019, 40(6): 2397-2406.
[5] 郭 超, 高永涛, 吴顺川, 成子桥, 张诗淮, 韩龙强, . 基于三维快速扫描算法与到时差数据库技术的层状介质震源定位方法研究[J]. 岩土力学, 2019, 40(3): 1229-1238.
[6] 温树杰,梁 超,宋亮亮,刘 刚,. 基于最小势能法的三维临界滑裂面搜索方法[J]. , 2018, 39(7): 2708-2714.
[7] 王开禾,罗先启,沈 辉,张海涛. 围岩力学参数反演的GSA-BP神经网络模型及应用[J]. , 2016, 37(S1): 631-638.
[8] 王 伟 ,杨 敏 ,上官士青,. 控制差异沉降的桩筏基础桩径优化分析方法[J]. , 2015, 36(S2): 178-184.
[9] 李 忠 ,杨 俊 , . 基于MPGA的复杂应力状态边坡稳定性分析[J]. , 2015, 36(5): 1488-1495.
[10] 李 宁 ,王李管 ,贾明涛 ,陈建宏 ,谭正华,. 改进遗传算法和支持向量机的岩体结构面聚类分析[J]. , 2014, 35(S2): 405-411.
[11] 邱道宏 ,李术才 ,薛翊国 ,田 昊 ,闫茂旺,. 基于数字钻进技术和量子遗传-径向基函数神经网络的围岩类别超前识别技术研究[J]. , 2014, 35(7): 2013-2018.
[12] 文 斌 ,吴青柏 ,蒋观利 ,张 鹏 . 模拟退火优化算法的冻土热传导参数反分析[J]. , 2013, 34(8): 2401-2408.
[13] 林永生,陈胜宏. 基于有限元计算的边坡三维滑裂面搜索[J]. , 2013, 34(4): 1191-1196.
[14] 袁文华. 人工冻土黏弹塑性本构参数反分析研究[J]. , 2013, 34(11): 3091-3095.
[15] 李南生,唐 博,谈风婕,谢利辉. 基于统一强度理论的土石坝边坡稳定分析遗传算法[J]. , 2013, 34(1): 243-249.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 姚仰平,侯 伟. 土的基本力学特性及其弹塑性描述[J]. , 2009, 30(10): 2881 -2902 .
[2] 徐金明,羌培,张鹏飞. 粉质黏土图像的纹理特征分析[J]. , 2009, 30(10): 2903 -2907 .
[3] 向天兵,冯夏庭,陈炳瑞,江 权,张传庆. 三向应力状态下单结构面岩石试样破坏机制与真三轴试验研究[J]. , 2009, 30(10): 2908 -2916 .
[4] 石玉玲,门玉明,彭建兵,黄强兵,刘洪佳. 地裂缝对不同结构形式桥梁桥面的破坏试验研究[J]. , 2009, 30(10): 2917 -2922 .
[5] 夏栋舟,何益斌,刘建华. 土-结构动力相互作用体系阻尼及地震反应分析[J]. , 2009, 30(10): 2923 -2928 .
[6] 徐速超,冯夏庭,陈炳瑞. 矽卡岩单轴循环加卸载试验及声发射特性研究[J]. , 2009, 30(10): 2929 -2934 .
[7] 张力霆,齐清兰,魏静,霍倩,周国斌. 淤填黏土固结过程中孔隙比的变化规律[J]. , 2009, 30(10): 2935 -2939 .
[8] 张其一. 复合加载模式下地基失效机制研究[J]. , 2009, 30(10): 2940 -2944 .
[9] 易 俊,姜永东,鲜学福,罗 云,张 瑜. 声场促进煤层气渗流的应力-温度-渗流压力场的流固动态耦合模型[J]. , 2009, 30(10): 2945 -2949 .
[10] 陶干强,杨仕教,任凤玉. 崩落矿岩散粒体流动性能试验研究[J]. , 2009, 30(10): 2950 -2954 .