采用自适应泰森多边形参数化的瞬变电磁二维贝叶斯反演方法

    公开(公告)号:CN114896564B

    公开(公告)日:2024-08-02

    申请号:CN202210566416.4

    申请日:2022-05-23

    IPC分类号: G06F17/18 G06F17/16 G06N7/01

    摘要: 本发明公开了一种采用自适应泰森多边形参数化的瞬变电磁二维贝叶斯反演方法,包括:设定相关反演参数;随机生成初始模型m0;设定模型的修改方式;随机挑选一种修改方式,生成新的建议模型m1;判断建议模型m1是否满足先验准则;如果模型m1满足先验准则,则进行一维正演模拟计算所有测点处模型m1的瞬变电磁响应F(m);根据瞬变电磁响应F(m)计算数据的似然度p(d|m1);若似然度p(d|m1)满足MH接受准则,则接受该新建议模型,并设置m0=m1。本发明尤其适用于瞬变电磁剖面二维反演,能在反演过程中自动调整模型参数的个数及尺寸去更好地适应观测数据的空间分辨率,从而避免了对模型参数与正则化参数的先验选择,同时贝叶斯方法能提供反演模型的不确定信息并更准确地确定实测数据的有效探测深度。

    一种基于舒尔补方法的瞬变电磁三维快速正演方法

    公开(公告)号:CN113569447A

    公开(公告)日:2021-10-29

    申请号:CN202110761321.3

    申请日:2021-07-06

    IPC分类号: G06F30/23 G06F17/11 G06F17/16

    摘要: 本发明公开了一种基于舒尔补方法的瞬变电磁三维快速正演方法,该方法将舒尔补方法成功引入到时域有限体积法瞬变电磁三维正演中,将整个模拟区域划分为背景区域与异常区域,其中背景区域在多正演过程中电导率不发生变化只有异常区域的电导率更新。依据区域划分位置,采用舒尔补分解可将正演方程分解为不同的部分,异常区域电导率更新后,更新模型的正演仅需求解异常区域所对应的舒尔补方程,这将大幅度提高瞬变电磁三维正演的效率。本发明尤其适用于对地下介质电导率分布已有较好了解,而未知电导率区域只占模拟区域的较小部分的情况,能降低正演过程中仅有小区域电导率变化时的系数矩阵方程规模,从而提高多模型瞬变电磁正演模拟效率。

    采用自适应泰森多边形参数化的瞬变电磁二维贝叶斯反演方法

    公开(公告)号:CN114896564A

    公开(公告)日:2022-08-12

    申请号:CN202210566416.4

    申请日:2022-05-23

    IPC分类号: G06F17/18 G06F17/16 G06N7/00

    摘要: 本发明公开了一种采用自适应泰森多边形参数化的瞬变电磁二维贝叶斯反演方法,包括:设定相关反演参数;随机生成初始模型m0;设定模型的修改方式;随机挑选一种修改方式,生成新的建议模型m1;判断建议模型m1是否满足先验准则;如果模型m1满足先验准则,则进行一维正演模拟计算所有测点处模型m1的瞬变电磁响应F(m);根据瞬变电磁响应F(m)计算数据的似然度p(d|m1);若似然度p(d|m1)满足MH接受准则,则接受该新建议模型,并设置m0=m1。本发明尤其适用于瞬变电磁剖面二维反演,能在反演过程中自动调整模型参数的个数及尺寸去更好地适应观测数据的空间分辨率,从而避免了对模型参数与正则化参数的先验选择,同时贝叶斯方法能提供反演模型的不确定信息并更准确地确定实测数据的有效探测深度。

    一种基于舒尔补方法的瞬变电磁三维快速正演方法

    公开(公告)号:CN113569447B

    公开(公告)日:2023-11-03

    申请号:CN202110761321.3

    申请日:2021-07-06

    IPC分类号: G06F30/23 G06F17/11 G06F17/16

    摘要: 本发明公开了一种基于舒尔补方法的瞬变电磁三维快速正演方法,该方法将舒尔补方法成功引入到时域有限体积法瞬变电磁三维正演中,将整个模拟区域划分为背景区域与异常区域,其中背景区域在多正演过程中电导率不发生变化只有异常区域的电导率更新。依据区域划分位置,采用舒尔补分解可将正演方程分解为不同的部分,异常区域电导率更新后,更新模型的正演仅需求解异常区域所对应的舒尔补方程,这将大幅度提高瞬变电磁三维正演的效率。本发明尤其适用于对地下介质电导率分布已有较好了解,而未知电导率区域只占模拟区域的较小部分的情况,能降低正演过程中仅有小区域电导率变化时的系数矩阵方程规模,从而提高多模型瞬变电磁正演模拟效率。