Research Of Complex Resistivity 2.5D Electromagnetic Forward And Inversion With Topography

The resistivity and induced polarization of geologic body generally show the characteristics of complex resistivity. Some scholars have expressed the complex resistivity by some mathematical models, for example the Cole-Cole model. It has been proved that the electrical anomaly of geologic body can be quantitative estimated by Cole-Cole model parameters (direct resistivityρ0, chargeability m , time constantτand frequency dependencec ). It is very important to study the methods to obtain the complex parameters of geologic body.There are a lot of mountains in china, mountainous area take part in two over three. In order to enhance the veracity and reasonably of the complex resistivity data, we must consider the effect of topography. At present among data interpretation of electromagnetic method, the computation quantity of 3D forward problem and inversion is enormous. In the mountainous area, the efficiency of 3D forward problem and inversion cannot obtain the enhancement but the cost increase. The 2D forward problem and inversion have small computation quantity and achieve many achievements of research. But because of the source is 3D or finite, the 2D forward problem and inversion are only approximate simulation. Therefor, 2.5D mode(l3D or finite source over an earth whose conductivity varies in 2D)simulation have particular predominance. Compared with 3D model, 2.5D model reduces the computation quantity, which can be fulfilled easily on common computers. Compared with 2D model, it is easier for 2.5D model with topography to approach field geologic condition. so the research on electromagnetic forward and inversion of 2.5D model with topography has important theoretical and practical significance. In practical, electromagnetic effect and induced polarization exist at the same time and influence each other. The electromagnetic coupling will become the interference factor on induced polarization method, while the induced polarization will become the interference factor on frequency domain and time domain electric methods. so it is hard to separate them. Considering all these problems, we develop the research on electromagnetic forward and inversion of 2.5D complex resistivity with topography.In this paper, we use the finite element method to compute the 2.5D electromagnetic response with topography, In finite element method solution, the solve region must be divided into small units, in this paper, we use the anomalous quadrangle grid to divide the region, and the realistic topography is realized by different altitude of node, furthermore the grid are small and denser in the middle area where the source and anomalous body distribute, while the grid are big and sparser at the edge of area. As a result, not only there is enough space for dividing area, but also complex parts of the geoelectric cross section gather in the center of dividing area. We deduce the frequency domain differential equation of electric field and magnetic field coupling according to Maxwell equation for this block homogeneous medium, then isoparametric elements and Galerkin method are used to disperse the differential equation, we obtain the linear equations which various nodes function value satisfies. Introducing Cole-Cole model into the 2.5D electromagnetic model,and introducing the pseudo-delta function to solve the source place irregularity. Moreover storage of system matrix employ the univariate nonzero element contraction which can greatly reduce the memory occupancy, the Cholesky preconditioned biconjugate gradient (ICBCG) method is used to solve the global system matrix equations. The integration is carried by continued fraction method in Fourier inverse transform.By model calculating, this article analyses the effect of source frequency, topography, and Cole-Cole model parameters (ρ0, m ,τ,c ) on electric and magnetic field at survey site. It concludes that: 1. the value of horizontal electric field is changing from positive to negative gradually with source frequency; the real component of electric field is a positive anomaly above the anomalous body when the direct resistivity of anomalous body is high resistivity, the performance is a negative anomaly whenρ0is low resistivity, the imaginary part of electric field is opposite to real part; The real component of horizontal electric field decreases with the chargeability above anomalous body, the imaginary part of electric field increases with chargeability at left side of anomalous body, but decreases at the right; The horizontal electric field decrease with time-constant above anomalous body; The real component of horizontal electric field increases with the frequency dependence above anomalous body, but the imaginary part decreases with it; The abnormity of terrain increases with the frequency; When identical frequency the real part of electric field is bigger than the imaginary part. The amplitude of horizontal electric field at the middle of topography almost don't change with direct resistivity, chargeability, time constant and frequency dependence. 2. the vertical magnetic field is a negative anomaly above the anomalous body when the direct resistivity of anomalous body is high resistivity, the performance is a positive anomaly whenρ0is low resistivity; The real component of vertical magnetic field decreases with the chargeability at left side of anomalous body, but increases at the right; the imaginary part of vertical magnetic field increases with chargeability at left side of anomalous body, but decreases at the right; The real part of vertical magnetic field increases with time-constant above anomalous body, while imaginary part decreases with time-constant; The real component of vertical magnetic field increases with the frequency dependence above anomalous body, but the imaginary part is opposite to real part. The amplitude of vertical magnetic field at the middle of topography almost don't change with direct resistivity, chargeability, time constant and frequency dependence of the anomalous body.As for the inversion of induced polarization parameter, this paper develops the partial derivative matrix that reflects the rate of real measurement field to the Cole-Cole model parameters by reciprocity theory. It is achieved that the complex resistivity 2.5D electromagnetic inversion of dipole source with topography by using classical damped least square method. In the inversion, the convergence rate of dc resistivity is fastest , the convergence rate of dependence is the second, and the convergence rate of time constant is slowest. On the other side, due to strong correlativity between time constant and chargeability, the time constant can't be accurate inversed when all the complex resistivity parameters take part in the inversion. For this problem, considering the chargeability can be measured or estimated before the inversion calculation. So we fixed the values of the chargeability in the inversion, and then inversed other resistivity parameters. Through the method we can obtain accurate result of inversion. The results of theoretical calculation example indicate that the method about inversing resistivity parameters of 2.5D geologic body from electromagnetic information is feasible.The production of research further developed the theory of forward and inversion for complex resistivity method. It has important theoretical and practical significance.