Parallel Calculation Of3D FDTD Helicopter Transient Electromagnetic Response Based On The GPU

With the rapid development of China’s economy, the demand for mineral resourcesalso increases dramatically. The contradiction between supply and demand has becomeincreasingly intensified, where the demand of the bulk minerals such as iron, manganese,chromium, copper, aluminum, potassium salt is particularly prominent. Helicoptertransient electromagnetic method is an exploration technology which is developed becauseof the special environment, for example, forest-swamp landscape area and desert-Gobilandscape area. In China, the shape of the mineral resource’s distribution of is diverse.Most of them present as lenticular, lentiform and vein type. However, the actual situationis just the opposite. The theoretical investigation of helicopter transient electromagneticmethod is mainly based on the uniform half-space model or layered earth model. On thebasis of that, one-dimensional or two-dimensional approximate inversion surely has awide gap comparing with the actual structure of the ore body. In order to reduce the error,we use a three-dimensional anomalous body model to approximate lenticular, lentiform etal mineral distribution. The following problems exist in3D finite difference numericalcalculation method. Firstly, there are no absorbing boundary conditions in the truncatedborders of computational domain. Artificially, we let the value of truncated borders beequal to zero. Do like this, if computational domain is not large enough, the field will bedivergence in the late stage. Secondly, the calculation time is long. So it is unable to meetthe rapid demand of the3D electromagnetic inversion.Based on the problems above, the thesis is in the support of863major tackle keyproblems in science and technology project "the bird time domain helicopter-borneelectromagnetic exploration system", and national natural science foundation project "research about the3D geological identify key technology based on the full waveformaviation time-domain electromagnetic detection". In the basis of serial computing of3Dfinite-difference transient electromagnetic response which doesn’t have the absorbingboundary conditions, we use heterogeneous platforms based on CPU and GPU to researchthe parallel computing of3D finite difference helicopter transient electromagneticresponse with Mur absorbing boundary conditions. Major researches and innovations areas follows:1．According to the expressing which the loop is on the surface and receiving isunderground, we derive time domain initial field formula of step incentives which the loopis in the air and receiving is underground. By means of the small amount of approximationmethod, making the formula only have one Bessel function. Applying the MAXWELLfirst and second equation, through differential discrete method, I derive pulse initialelectromagnetic field which is needed in the iteration field computation. At last, realizeinitial field parallel computation. When the grid is101×101×50, acceleration ratio is21.04.When the grid is201×201×150, acceleration ratio is13.2. The denser the grid is, the moretime can be saved by parallel calculation.2． Research the calculation method of bicubic spline interpolation under thenot-a-knot boundary condition: decompose the bicubic spline interpolation into two cubicspline interpolations, interpolating along the x direction first, and then y direction. At last,the row or column of two-dimensional array to be interpolated is assigned to Thread in theGPU to realize parallel calculation of interpolation. Ultimately, the parallel computingresult and interpolation function in the matlab are compared; the error is controlledbetween10-4%~10-8%.3．Using3D finite differential time domain to discrete curl equations in theMAXWELL equation，I derived x, y, z components of electric and magnetic field in thecomputational domain. Adopting3D Mur absorbing boundary conditions, making use offirst order and second order approximation derive electric field formula on the truncatedboundary. Solving the problem, when the rock resistivity is larger, using no absorbingboundary conditions to calculate3D electromagnetic response will produce the divergence in the late time. Effectively extend the electromagnetic response calculation time andimprove accuracy. Applying upward-continuation theory solves the problem aboutmagnetic field component calculation in the air. At last, achieve parallel computing ofiteration field, acceleration ratio is12.03and7.32respectively, and supply the foundationfor3D inversion of anomaly.4. Transmitter-receiver loop applies the central loop device. I studies helicoptertransient electromagnetic response characteristics of layered earth model, plate-like bodywith or without covering layer model and3D abnormal body model. Supply forwardsample data for inverse explain of3D target. According to the abnormal morphology ofprofile curve of plate-like body in the homogeneous half-space model, the tendency of theplate-like can be roughly determined. Single peak is a horizontal plate-like body,equivalent double peak is an upright plate-like body, and the non-equivalent double peakis an inclined plate-like body. The meters changing which the aircraft deviates from thenormal flight altitude will cause false anomaly, and therefore it need to be corrected. Bythe thinking of distinguished time window, I research the detection capability of helicoptertransient electromagnetic detection system. Only changing the conductivity of theanomalous body, the other parameters are constant, then the greater the conductance ratiobetween the anomaly and surrounding rock, the more easily detected the anomaly. Whenonly change the depth of anomaly, the shallower the anomaly is, the more easily detectedthe anomaly.5. Combining the G-S inverse Laplace transform with F.N.Kong241’s Hankeltransform calculates helicopter transient electromagnetic response of the layered earthwhen transmitted current is a ramp function and is completely switched off. By theconvolution between the first derivative of trapezoidal wave and the3D anomaly responsewhen transmitted current is a step function, I calculate3D anomaly helicopter transientelectromagnetic response when the transmitted current is a trapezoidal form.The innovation of the paper are as follows：1. Put forward the bicubic spline interpolation computing method under not-a-knotboundary condition. Solve the problem that we can’t calculate the first and second order derivative of magnetic field value at two end points of the computing area and solveinterpolation calculation problem in the upward-continuation, to lay the groundwork forthe GPU-based iterative field of parallel computing.2. Put forward applying the truncated boundary’s electric field formula to computeiteration field of3D finite difference helicopter transient electromagnetic response.Solving the problem, when the rock resistivity is larger, using no absorbing boundaryconditions to calculate3D electromagnetic response will produce the divergence problemin the late time. Effectively improve the accuracy of numerical computation ofelectromagnetic lately.3. By analyzing numerical calculation features of initial field, electromagnetic fieldwithin the region, truncated boundary’s electric field and magnetic field in the air, putforward applying parallel computing method based on GPU to initial field and iterationfield of electromagnetic response. Complete the parallel calculation of helicopter transientelectromagnetic response of3D finite differential method based on the GPU, savingcomputing time.