[发明专利]瞬变电磁三维有限差分正演方法及系统

权利要求书

1.一种瞬变电磁三维有限差分正演方法，其特征在于，包括：

根据地电模型的模型参数，设置时间步长为初始时间步长；

基于所述初始时间步长，进行预设次数的时间步进迭代，并执行三维有限差分正演模拟；

将所述时间步长设置为当前时刻的目标时间步长；

基于所述目标时间步长，继续进行时间步进迭代，并执行三维有限差分正演模拟，直至时间步进迭代的总时长达到预设时长，三维有限差分正演模拟过程结束；

其中，所述目标时间步长基于当前时刻所述地电模型中地下介质的视电导率确定。

2.根据权利要求1所述的瞬变电磁三维有限差分正演方法，其特征在于，所述目标时间步长具体基于如下方式确定：

基于当前时刻的所述视电导率，确定当前时刻的备选目标时间步长；

基于所述备选目标时间步长，确定所述目标时间步长。

3.根据权利要求2所述的瞬变电磁三维有限差分正演方法，其特征在于，所述基于当前时刻的所述视电导率，确定当前时刻的备选目标时间步长，具体包括：

基于当前时刻的所述视电导率、所述地下介质的最小磁导率、所述地电模型中离散的最小网格间距以及当前时刻进行的时间步进迭代中电流完全关断后的延迟时间，确定所述备选目标时间步长。

4.根据权利要求2所述的瞬变电磁三维有限差分正演方法，其特征在于，所述基于所述备选目标时间步长，确定所述目标时间步长，具体包括：

若判断获知当前时刻的视电导率增大，则基于第一预设倍数的前一时刻的时间步长，以及所述备选目标时间步长之间的最小值，确定所述目标时间步长；

若判断获知当前时刻的视电导率减小，则基于第二预设倍数的前一时刻的时间步长，以及所述备选目标时间步长之间的最大值，确定所述目标时间步长；

其中，所述第一预设倍数大于所述第二预设倍数。

5.根据权利要求1-4中任一项所述的瞬变电磁三维有限差分正演方法，其特征在于，所述视电导率基于每个时间量级范围内进行预设次数的二分法搜索计算得到；所述二分法搜索计算的停止条件包括：

基于当前搜索计算得到的所述视电导率，执行三维有限差分正演模拟得到的模拟结果数据与所述地电模型的观测数据之间的相对偏差小于预设阈值。

6.根据权利要求5所述的瞬变电磁三维有限差分正演方法，其特征在于，所述二分法搜索计算的停止条件具体包括：

其中，d_app为所述模拟结果数据，d_obs为所述观测数据，D为所述预设阈值。

7.根据权利要求1-4中任一项所述的瞬变电磁三维有限差分正演方法，其特征在于，所述根据地电模型的模型参数，设置时间步长为初始时间步长，具体包括：

基于空气中的磁导率、所述地电模型中离散的最小网格间距以及所述地下介质的最小电导率，确定所述初始时间步长。

8.一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至7任一项所述瞬变电磁三维有限差分正演方法的步骤。

9.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至7任一项所述瞬变电磁三维有限差分正演方法的步骤。

10.一种计算机程序产品，包括计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至7任一项所述瞬变电磁三维有限差分正演方法的步骤。