Joint inversion of electromagnetic and seismic data based on structural constraints using variational born iteration method

An efficient 2-D joint full-waveform inversion method for electromagnetic and seismic data in a layered medium background is developed. The joint inversion method based on the integral equation (IE) method is first proposed in this paper. In forward computation, the IE method is employed, which usually has smaller discretized computation domain and less cumulative error compared with the finite-difference method. In addition, fast Fourier transform is used to accelerate the convolution between Green's functions and induced sources due to the shift invariance property of the layered Green's functions in the horizontal direction. In the inversion model, the crossgradient function is incorporated into the cost function of the separate inversion to enforce the structure similarity between electric conductivity and seismic-wave velocity. We use the improved variational Born iteration method and two different iteration strategies to minimize the cost function and reconstruct the contrasts. Several typical models in geophysical applications are used to validate our joint inversion method, and the numerical simulation results show that joint inversion can improve the inversion results when compared with those from the separate inversion.