Please wait a minute...
E-mail Alert Rss
 
物探与化探  2022, Vol. 46 Issue (3): 743-749    DOI: 10.11720/wtyht.2022.1238
  方法研究·信息处理·仪器研制 本期目录 | 过刊浏览 | 高级检索 |
基于Python语言的瑞利面波频散反演实行方案
吴为治(), 娄利, 王鹏, 王宾
安徽省勘查技术院,安徽 合肥 230041
A Python-based scheme of Rayleigh-wave dispersion inversion
WU Wei-Zhi(), LOU Li, WANG Peng, WANG Bin
Geological Exploration Technology Institute of Anhui Province,Heifei 230041,China
全文: PDF(2912 KB)   HTML
输出: BibTeX | EndNote (RIS)      
摘要 

本文基于 Python 语言构建了瑞利面波频散反演工作流程:①通过pysurf96软件包实现水平层状模型的频散曲线正演;②建立描述频散曲线拟合度的目标函数;③利用scikit-opt软件包中的启发式算法实现频散曲线反演。提出并解决了调用函数过程中所遇到的问题,实验测试表明反演结果可靠,同时具备一定的计算效率。由此实现了搭建基于Python语言的面波频散反演地下层状结构研究平台,为其他研究者依靠开源软件进行反演运算提供了方法支持。最后,通过利用彭一波研究海拉尔盆地背景噪声时提取的频散曲线对地壳以及上地幔结构进行了反演,得到不错的效果。

服务
把本文推荐给朋友
加入引用管理器
E-mail Alert
RSS
作者相关文章
吴为治
娄利
王鹏
王宾
关键词 Python语言瑞利面波频散反演启发式算法    
Abstract

This study developed a workflow of Rayleigh-wave dispersion inversion using the Python programming language, and the detail is as follows.First,carry out the forward modeling of the dispersion curves on a horizontal layered model using the pysurf96 software package.Second,create an objective function used to describe the fitting degree of the dispersion curves.Third,complete the dispersion curve inversion using the heuristic algorithm in the scikit-opt software package.The problems encountered in the function call in the workflow were proposed and solved.The results show that the Python-based dispersion curve inversion of Rayleigh wave in multilayered media is reliable and offers a certain computational efficiency.In this way,this study built a Python-based inversion platform of underground layered structures using the wave dispersion,thus providing a method for other researchers to do inversion using open-source software.Finally,this study carried out the inversion of the crust and upper mantle structures using the dispersion curves extracted from the study of Yi-bo Peng on the noise in the Hailar Basin,achieving ideal results.

Key wordsPython library    Rayleigh wave    dispersion curve inversion    heuristic algorithm
收稿日期: 2021-08-12      修回日期: 2022-03-24      出版日期: 2022-06-20
ZTFLH:  P631  
基金资助:安徽蒙城地球物理国家野外科学观测站联合开放基金(MENGO-202111)
作者简介: 吴为治(1989- ),男,2015年毕业于中国科学技术大学,主要从事地质工程方面的工作。Email: 1781367830@qq.com
引用本文:   
吴为治, 娄利, 王鹏, 王宾. 基于Python语言的瑞利面波频散反演实行方案[J]. 物探与化探, 2022, 46(3): 743-749.
WU Wei-Zhi, LOU Li, WANG Peng, WANG Bin. A Python-based scheme of Rayleigh-wave dispersion inversion. Geophysical and Geochemical Exploration, 2022, 46(3): 743-749.
链接本文:  
https://www.wutanyuhuatan.com/CN/10.11720/wtyht.2022.1238      或      https://www.wutanyuhuatan.com/CN/Y2022/V46/I3/743
层厚/m vs/(m·s-1) vp/(m·s-1) 密度/(kg·m-3)
10 200 400 2000
半空间 400 800 2000
Table 1  两层层状介质模型参数
Fig.1  基于pysurf96计算的两层层状介质模型的瑞利面波频散曲线
Fig.2  不同启发式算法反演的结果
a—差分进化算法;b—遗传算法;c—粒子群算法;d—模拟退火算法
Fig.3  频率为2 Hz瑞利面波相速度敏感核
Fig.4  不同启发式算法的效率和稳定性
a—差分进化算法;b—遗传算法;c—粒子群算法;d—模拟退火算法
Fig.5  遗传算法反演海拉尔盆地面瑞利波频散结果
[1] 彭一波, 姜明明, 艾印双. 基于Python语言的ObsPy软件包从地震背景噪声中提取瑞利面波经验格林函数的实行方案[J]. 地球物理学进展, 2019, 34(3):919-927.
[1] Peng Y B, Jiang M M, Ai Y S. Efficient scheme to extract Green functions of Rayleigh wave from seismic noise via a Python library for seismology ObsPy[J]. Progress in Geophysics, 2019, 34(3):919-927.
[2] 夏江海. 高频面波方法[M]. 武汉: 中国地质大学出版社, 2015.
[2] Xia J H. High-frequency surface waves method[M]. Wuhan: China University of Geosciences Press, 2015.
[3] Knopoff L. A matrix method for elastic wave problems[J]. Bulletin of the Seismological Society of America, 1964, 54(1):431-438.
doi: 10.1785/BSSA0540010431
[4] Chen X. A systematic and efficient method of computing normal modes for multilayered half-space[J]. Geophysical Journal International, 1993, 115(2):391-409.
doi: 10.1111/j.1365-246X.1993.tb01194.x
[5] Herrmann R B. Computer programs in seismology:An evolving tool for instruction and research[J]. Seismological Research Letters, 2013, 84(6):1081-1088.
doi: 10.1785/0220110096
[6] Yang T, Yi W. A Study of leaky mode wave and Zig-Zag dispersion curve in Rayleigh-wave exploration[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2005, 27(4):295-300.
[7] Pan Y, Xia J, Zeng C. Verification of correctness of using real part of complex root as Rayleigh-wave phase velocity with synthetic data[J]. Journal of Applied Geophysics, 2013, 88(1):94-100.
doi: 10.1016/j.jappgeo.2012.09.012
[8] 凡友华, 刘雪峰, 陈晓非. 面波频散反演地下层状结构的拟牛顿法[J]. 物探与化探, 2006, 30(5):456-459.
[8] Fan Y H, Liu X F, Chen X F. The Quasi Newton method in the inversion of the dispersion curve of Rayleigh wave in multilayered media[J]. Geophysical and Geochemical Exploration, 2006, 30(5):456-459.
[9] 石耀霖, 金文. 面波频散反演地球内部构造的遗传算法[J]. 地球物理学报, 1995, 38(2):189-198.
[9] Shi Y L, Jin W. Genetic algorithms inversion of lithospheric structure from surface wave dispersion[J]. Chinese Journal of Geophysics, 1995, 38(2):189-198.
[10] 翟佳羽, 赵园园, 安丁酉. 面波频散反演地下层状结构的蚁群算法[J]. 物探与化探, 2010, 34(4):476-481.
[10] Zhai J Y, Zhao Y Y, An D Y. The ant colony algorithm for the inversion of the disperson curve of Rayleigh wave in multilayered media[J]. Geophysical and Geochemical Exploration, 2010, 34(4):476-481.
[11] Zhang J, Zhang H, Chen E, et al. Real-time earthquake monitoring using a search engine method[J]. Nature Communications, 2014, 5(1):1-9.
[12] 李英康, 高锐, 姚聿涛, 等. 大兴安岭造山带及两侧盆地的地壳速度结构[J]. 地球物理学进展, 2014, 29(1):73-83.
[12] Li Y K, Gao R, Yao Y T, et al. The srust velocity structrue of DA Hinggan Ling orgenic belt and the basins on both sides[J]. Progress in Geophysics, 2014, 29(1):73-83.
[13] 李皎皎, 黄金莉, 刘志坤. 用背景噪声和地震面波反演东北地区岩石圈速度结构[J]. 地震, 2012, 32(4):22-32.
[13] Li J J, Huang J L, Liu Z K. Inversion of lithospheric velocity structure in Northeast China using ambient noise and seismic surface wave[J]. Earthquake, 2012, 32(4):22-32.
[14] Zha Y, Webb S C. Crustal shear velocity structure in the Southern Lau Basin constrained by seafloor compliance[J]. Journal of Geophysical Research:Solid Earth, 2016, 121(5):3220-3237.
doi: 10.1002/2015JB012688
[1] 许艺煌, 黄真萍, 程志伟, 陈少博, 陈振明. 高密度电阻率法在弃渣堆积体分布调查中的应用[J]. 物探与化探, 2020, 44(2): 435-440.
[2] 熊章强, 唐圣松, 张大洲. 瑞利面波数值模拟中的PML吸收边界条件[J]. 物探与化探, 2009, 33(4): 453-457.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址:北京市学院路29号航遥中心 邮编:100083
电话:010-62060192;62060193 E-mail:whtbjb@sina.com