본문 바로가기 주메뉴 바로가기
검색 검색영역닫기 검색 검색영역닫기 ENGLISH 메뉴 전체보기 메뉴 전체보기

논문

Efficient electric resistivity inversion using adjoint state of mixed finite-element method for Poisson`s equation

https://doi.org/10.1016/j.jcp.2005.09.007

  • 저자Taeyoung Ha,Sukjoon Pyun and Changsoo Shin
  • 학술지Journal of Computational Physics 214/1
  • 등재유형
  • 게재일자(2006)


We propose an electric resistivity inversion method that is similar to the reverse time migration technique applied to seismic data. For calculating model responses and inversion, we use the mixed finite-element method with the standard P1  P0 pair for triangular decompositions, which makes it possible to compute both the electric potential and the electric field vector economically. In order to apply the adjoint state of the Poisson equation in the resistivity inverse problem, we introduce an apparent electric field defined as the dot product between the computed electric field vector and a weighting factor and then defining a virtual source to compute the partial derivative of the electric field vector. We exploit the adjoint state (the symmetry of Green’s function) of matrix equations derived from solving the Poisson equation by the mixed finite-element method, for the calculation of the steepest descent direction of our objective function. By computing the steepest descent direction by a dot product of backpropagated residual and virtual source, we can avoid the cumbersome and expensive process of computing the Jacobian matrix directly. We calibrate our algorithm on a synthetic of a buried conductive block and obtain an image that is compatible with the limits of the resistivity method.


We propose an electric resistivity inversion method that is similar to the reverse time migration technique applied to seismic data. For calculating model responses and inversion, we use the mixed finite-element method with the standard P1  P0 pair for triangular decompositions, which makes it possible to compute both the electric potential and the electric field vector economically. In order to apply the adjoint state of the Poisson equation in the resistivity inverse problem, we introduce an apparent electric field defined as the dot product between the computed electric field vector and a weighting factor and then defining a virtual source to compute the partial derivative of the electric field vector. We exploit the adjoint state (the symmetry of Green’s function) of matrix equations derived from solving the Poisson equation by the mixed finite-element method, for the calculation of the steepest descent direction of our objective function. By computing the steepest descent direction by a dot product of backpropagated residual and virtual source, we can avoid the cumbersome and expensive process of computing the Jacobian matrix directly. We calibrate our algorithm on a synthetic of a buried conductive block and obtain an image that is compatible with the limits of the resistivity method.

이 페이지에서 제공하는 정보에 대해 만족하십니까?