地震波阻抗反演可以直接提供地下介质的弹性信息，为后续的煤层识别，水通道预测以及陷落柱识别等提供可靠的资料。在面对复杂地质条件时，现有的波阻抗反演方法精度较低，不能满足实际勘探的需求。反射率法通过求解弹性波动方程获得弹性参数信息，能够模拟全波场响应以及地震波的各种传播效应，能更精确地考虑透射损失、多次波等引起的地震波振幅和相位变化，精度较高，同时计算效率高于其他波动方程求解方法。因此，基于贝叶斯框架，建立了反射率法的反演目标函数并通过引入柯西分布作为先验分布提高反演结果的分辨率，形成了一种基于贝叶斯的反射率法反演方法。将该方法应用于布尔台煤矿地震数据中，获得的波阻抗反演结果与测井资料匹配度较好，相比于常规方法分辨率更高，精度高，能够有效识别薄煤层和深部煤层分布，有效验证了新方法的可行性和有效性。该方法能为利用地震反演技术识别煤层及陷落柱等小构造，预测顶底板水分布提供有力的资料。
2022-02-22 13:28:11
4.71MB
1
描述： bulirsch-Stoer 单步 ODE 传播器已在 MATLAB 中Juergen Dietel 的数值 ODE 求解器的 MEX 改编形式。 这很可能是 Mathworks FileExchange 上发布的最快的单步数值传播器。
所有函数都是用 C++ 编写的，必须编译才能在 MATLAB 上运行。 此代码已经过验证，可以使用 MS Visual Studio 2010 正确编译。
请使用以下命令序列进行编译： >> mex BulirschMex.cpp bulirsch.cpp base_r.cpp t_dgls.cpp vmblock.cpp
有关使用 Burlisch-Stoer 与 ODE113 的示例轨道传播比较，请参阅 testDriver.m。
2022-02-13 21:35:31
31KB
1