基于Python语言的瑞利面波频散反演实行方案

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

A Python-based scheme of Rayleigh-wave dispersion inversion

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.