matlab瑞利波频散曲线代码瑞利波蒙特卡罗反演
一种联合反演R波频散曲线的代码
MC_R_CyU_VariableThickness.m
此代码执行瑞利波频散曲线(相位和群联合)的蒙特卡罗反演(马尔可夫链
MC)。
它搜索适合观察数据的一维高概率模型。
该代码使用
Herrmann
地震学代码中的
surf96
()。
这组代码应该在尝试运行反演之前安装。
与任何
MC
代码一样,运行它需要时间!
速度取决于您的计算机和
Matlab
版本。
根据迭代次数,反演可能需要长达一天的时间,但经验表明
100000
次迭代就可以解决问题。
要有耐心!
请阅读:Bosch
(1999;
2001
and
2005)
和
Mosegaard
&
Sambridge
(2002)
了解方法和背景。
代码从一个模型开始(要修改的层和要考虑的层数应按照说明进行设置),并且很难找到解决方案。
这是通过选择一层或一组层(随机)并改变层的速度或厚度来完成的。
图层不允许出现或消失。
如果“扰动”减少了错误,则它被接受,否则它有被接受或拒绝的概率。
如果被拒绝,则重复先前的模型。
参数:
IT:迭代次数
M
2021-08-20 16:58:13
126KB
系统开源
1