社団法人 物理探査学会
第141回(2019年度秋季)学術講演会


MT法の1-Dスパースインバージョン

講演要旨(和文)
MTのインバージョンでは,未知パラメータの平滑化拘束を用いて非線形最小二乗法を解く方法が一般的である.平滑化拘束を用いれば解を得られる反面,その解は滑らかに変化する地下構造に限定される.この問題を解決するため、スパースインバージョンが有効であると考えられる.本研究では,非線形最小二乗法を解くMT法にスパースインバージョンを適用し,1次元スパースインバージョンを開発した.スパースインバージョンでは,L1ノルムを最小にするという条件を逆問題に与えることで,スパースな解が選択される.開発したプログラムでは,L1ノルム最適化には近接勾配アルゴリズムとしてFISTAを用いた.地下構造モデルとしては,2層,3層モデルを想定し,模擬データを用いてスパースインバージョンを行い,本プログラムの有効性を検討した.

講演要旨(英文)
Magnetotelluric (MT) method is a powerful geophysical tool for mapping and sounding subsurface resistivity distributions. This method is widely used in the exploration for geothermal resources. MT 1-D inversion is usually performed by using nonlinear least squares method with smoothness constrained requirement. However, we often have to take into consideration subsurface structure which has sharp boundaries such as a fault. In order to solve this problem, we developed 1-D sparse inversion program based on Fast Iterative Shrinkage Thresholding Algorithm (FISTA). In this program, we applied the L1 norm minimization to nonlinear least squares method. Some examples of MT 1-D sparse inversion using synthetic data are shown in this paper.