40 results about "Time domain electromagnetics" patented technology

The invention relates to an active source magnetization-time domain electromagnetic induction polarization multi-parameter extraction method, and provides an active source magnetization-time domain electromagnetic induction polarization multi-parameter extraction method aiming at induction response caused by conductivity, polarization response caused by chargeable property and magnetization response caused by magnetic complex physical property in polymetallic ore detection. When the emission current is in a flat-topped stage, based on an active magnetic method, magnetic susceptibility extraction is carried out by using an Oxim method; when the emission current is in a turn-off stage, resistivity and polarization extraction is carried out by using a depth Q network method based on a time domain electromagnetic method by taking the extracted magnetic susceptibility as a constraint condition; and finally, multi-parameter extraction of resistivity, polarizability and magnetic susceptibility of the active source magnetization-time domain electromagnetic induction polarization method is realized.

The invention discloses a radiation source individual identification method, system and device based on space-time information separation, and relates to the field of radiation source individual identification, and the method comprises the steps: obtaining a to-be-detected time domain electromagnetic signal; inputting a to-be-detected time domain electromagnetic signal into the trained radiation source individual identification model to obtain a radiation source individual, and completing the identification of the radiation source individual; the radiation source individual identification model comprises a signal preprocessing and feature enhancement module, a multi-scale spatial feature coding module, a time sequence dynamic modeling module and a space-time cross attention fusion module which are connected in sequence; a loss function for training of the radiation source individual identification model adds triple loss on the basis of standard cross entropy loss. According to the invention, the accuracy and the recognition efficiency of radiation source individual recognition can be improved.

The invention relates to the technical field of geophysical electromagnetic data exploration, in particular to a three-dimensional time-frequency domain electromagnetic joint inversion method based on a finite element algorithm. The method comprises the following steps: respectively acquiring time domain and frequency domain electromagnetic observation data of a to-be-measured region; dividing a to-be-measured region to obtain a non-structural tetrahedral mesh to construct a time-frequency domain inversion initial model; performing forward modeling calculation on the time-frequency domain inversion initial model to obtain time domain and frequency domain electromagnetic response data; obtaining a time domain electromagnetic data residual error and a frequency domain electromagnetic data residual error, respectively obtaining corresponding adaptive weights, and establishing a joint inversion objective function; performing iterative solution on the joint inversion target function, and calculating an update vector of the inversion initial model; and iteratively updating the initial inversion model until a set condition is met, and completing the time-frequency domain electromagnetic joint inversion. According to the method, the accuracy of three-dimensional inversion is effectively improved through joint inversion, so that the practical effect of three-dimensional inversion is enhanced.

The invention discloses a transient electromagnetic data polarizability extraction method based on frequency domain resistivity, and the method comprises the following steps: designing measuring line and wire emission source positions in a test mining area, selecting observation parameters, and collecting field electromagnetic data, the field electromagnetic data comprising transient electromagnetic method electric field data and frequency domain electromagnetic method electric field data; performing inversion on actually measured frequency domain electromagnetic method electric field data by adopting an Occam inversion algorithm to obtain resistivity characteristics of different measuring points, and further obtaining resistivity distribution information; and taking the resistivity distribution information obtained by inversion as a constraint condition of transient electromagnetic method electric field data inversion, performing transient electromagnetic method polarizability inversion by analyzing a relationship between the resistivity and a time constant, and extracting the polarizability. By analyzing and processing the frequency domain electromagnetic data and the time domain electromagnetic data, the resistivity and the polarizability can be accurately extracted, parameter separation is achieved, and the problem of traditional time domain inversion coupling is solved.

The invention provides an average value model establishment method suitable for small disturbance stability research of a string type photovoltaic grid-connected system, and the method comprises the steps: firstly, carrying out the mathematical modeling of a string type photovoltaic converter SPWM modulation process; then, respectively establishing an alternating current side equivalent model and a direct current side equivalent model of the average value model of the converter; and finally, establishing a string photovoltaic independent grid-connected system through a two-stage boosting transformer in PSCAD / EMTDC electromagnetic transient simulation software, and comparing time domain electromagnetic transient simulation results of a detailed model and an average value model which are adopted by converter modulation in the same small disturbance instability scene, so as to verify the accuracy of the average value model. According to the method, the average value model of the string type photovoltaic converter is established, the problems of low simulation efficiency and overlarge calculation energy consumption caused by detailed modeling in the converter modulation process can be effectively avoided, and the method is suitable for small disturbance stability research under string type photovoltaic sudden increase of a complex alternating current system.

The invention relates to the field of geophysical exploration electromagnetic method detection, in particular to a time domain electromagnetic electrical source emission system. The objective of the invention is to realize impedance matching of an electrical source equivalent load and effective suppression of power supply ripples generated by a direct-current power supply bus. Comprising a high-power direct-current power supply, a transmitting host and a fractional order ripple suppression circuit, and the high-power direct-current power supply provides three-phase alternating current; the transmitting host comprises an H-bridge inverter circuit, the H-bridge inverter circuit is connected with a transmitting cable and is alternately conducted under the control of a group of complementary pulse width modulation signals, and transmitting current is generated in a loop formed by the transmitting cable and the ground; the fractional-order ripple suppression circuit comprises a matching capacitor and a first fractional-order device which are connected in series, and the matching capacitor and the first fractional-order device are connected to the two ends of the H-bridge inverter circuit in parallel. And the performance of the time-domain electromagnetic electrical source emission system can be improved.

The application discloses a semi-airborne full-time electromagnetic forward method based on an Octree grid, which performs spatial discretization on time-domain Maxwell equations, thereby effectively reducing unknowns of equations to be solved. For full-time electromagnetic field numerical simulation, the application combines the high calculation precision of the second-order backward Euler method with the high OFF-Time calculation efficiency of the model reduction algorithm, adopts the second-order BDF and the model reduction algorithm to respectively solve ON-Time and OFF-Time time-domain control equations, forms a joint algorithm to realize simulation of semi-airborne time-domain electromagnetic full-time response, and finally completes the forward calculation process. Through numerical simulation verification, in the full-time transient electromagnetic forward process, the application can keep small occupied memory and also has higher calculation efficiency, and the calculation precision is not lower than that of other methods.

The invention relates to the technical field of aviation electromagnetic detection, in particular to a rapid time domain aviation electromagnetic data variable-dimension Bayesian inversion method, system and device and a medium. According to the method, a mapping relation between geoelectric model parameters and time domain electromagnetic response data is constructed by adopting a one-dimensional convolutional neural network according to a deep learning network, and rapid prediction of the time domain aviation electromagnetic response data in a variable-dimension Bayesian inversion process is realized; secondly, during parallel calculation, only one Markov chain is operated in the preheating period, and N Markov chains are executed in parallel after the preheating period is ended, so that the parallel strategy can avoid the problems that the N Markov chains are executed in parallel from the beginning of the preheating period, and the sampling of the Markov chains is not convergent and the calculation time is wasted due to a relatively poor initial state; the condition that Markov chain sampling is not convergent can be found in time in the sampling process of the preheating period, so that the next operation is stopped, and wrong calculation is avoided.

The present invention relates to the field of geophysical exploration technology, and in particular to a physical simulation method, device and system for a time domain electromagnetic method. The system comprises an electromagnetic signal transmitter, a DC stable power supply, a receiver and a transient electromagnetic signal acquisition instrument. The receiver is used to carry the surrounding rock part and the abnormal body part of the manufacturing material of the physical simulation model; the DC stable power supply is used to power the electromagnetic signal transmitter, and the electromagnetic signal transmitter is used to transmit a bipolar rectangular pulse current. The pulse width and off time of the bipolar rectangular pulse current are determined according to the time constant of the abnormal body part; and according to the parameter characteristics of the surrounding rock part, the diffusion distance of the transmitted signal in the surrounding rock part over time is determined, and a data collection time period is selected according to the diffusion distance. Then, the receiver synchronously collects electric field data and magnetic field data within the data collection time period through a receiving electrode and a flux gate arranged on its surface, and forwards the data to the transient electromagnetic signal acquisition instrument for analysis.

The application discloses a kind of super electrical big platform time domain electromagnetic environment and electromagnetic interference simulation prediction method, comprising the following steps: Fourier transform is carried out to time domain excitation signal, and the frequency spectrum curve corresponding to periodic time domain excitation signal is obtained;The obtained frequency spectrum curve is dispersed based on transmission path tracking, and the specific discrete frequency point to be calculated is obtained;Super electrical big platform is decomposed, and full-wave frequency domain electromagnetic environment and electromagnetic interference are solved in parallel;Frequency domain electromagnetic environment and electromagnetic interference are simulated, and simulation data are corrected, and the sampling frequency point electromagnetic environment and electromagnetic interference simulation data are obtained;The sampling frequency point electromagnetic environment and electromagnetic interference simulation data of the obtained super electrical big platform are multiplied with the corresponding frequency point of the frequency spectrum data of excitation source time domain signal, then inverse Fourier transform is carried out, and the super electrical big platform time domain electromagnetic environment and electromagnetic interference coupling signal are obtained.The method provided in the application solves the super electrical big platform time domain electromagnetic environment simulation prediction problem.

The invention discloses a vehicle permanent magnet synchronous motor vibration transmission capability analysis method, which establishes a simulation calculation model of a permanent magnet synchronous motor, and specifically comprises the following steps: step 1) calculating a cogging torque peak-to-peak value of the permanent magnet synchronous motor; step 2) calculating a counter electromotive force harmonic amplitude ratio of the permanent magnet synchronous motor; 3) calculating the time domain electromagnetic force of the permanent magnet synchronous motor; 4) calculating the tangential force density of the permanent magnet synchronous motor; 5) calculating the radial force density of the permanent magnet synchronous motor; 6) extracting stator structure parameters of the permanent magnet synchronous motor; and 7) analyzing the vibration transmission capability of the permanent magnet synchronous motor. According to the method, under the condition that high-frequency electromagnetic noise simulation is not carried out, rapid evaluation of the electromagnetic noise suppression capability of a motor design scheme can be realized, so that the development period is shortened, and the cost is saved.

The invention relates to the technical field of time-domain electromagnetic geological exploration, in particular to a time-domain superconducting electromagnetic multi-effect multi-parameter extraction method and system, and the method comprises the steps: calculating time-domain superconducting electromagnetic induction-polarization-magnetization multi-effect response, and constructing a data error fitting term based on the multi-effect response and the residual error of measured data; constructing a constraint matrix of multi-parameter correlation coefficients including resistivity, polarizability and magnetic susceptibility, and constructing a multi-parameter statistical constraint term for describing statistical correlation among multiple parameters based on the constraint matrix; constructing a regularization objective function; an improved Gaussian-Newton optimization algorithm is adopted, the regularization objective function is minimized, the optimal value of the regularization objective function is iteratively solved, and the resistivity, the polarizability and the magnetic susceptibility are obtained through synchronous inversion. According to the invention, multi-parameter efficient and accurate extraction imaging in the deep mineral exploration and geological structure identification field is realized.

A time domain electromagnetic tensor based logging method, relating to the field of induction logging. A transmitting coil system Tx and a receiving coil system Rx are each composed of three orthogonal direction coils, and the transmitting coil system and the receiving coil system form a borehole probe; a time domain electromagnetic field is excited by means of two orthogonal transmitting coil systems Tx, and a secondary field signal is received by two orthogonal receiving coil systems Rx located above transmitting coils and having a distance D; during each measurement, three directions are used for transmitting, nine component data is received at the same time, and each component is a series of induction potentials in one-to-one correspondence with measurement time sequences, and a potential or resistance change, from near to far, of the rock stratum where a measurement point is located is obtained after time-depth conversion.

The invention provides a time domain finite element fast algorithm for electromagnetic force coupling simulation of a loudspeaker, which comprises the following steps of: performing three-dimensional modeling on a micro loudspeaker, setting parameter information of a solving region, and extracting a calculation model; performing discretization processing on the micro loudspeaker model by using a tetrahedral mesh; an electromagnetic force coupling discrete process is deduced, a hybrid finite element fast time domain solver is constructed, a discrete model and solving parameters of the miniature loudspeaker are imported into the solver, and the solver is initialized; and carrying out coil geometric analysis, calculating the current density in the coil, taking the current density as an excitation source of a magnetic field, calling a hybrid finite element fast time domain solver, calculating to obtain the displacement of the vibrating diaphragm, and finally obtaining the time domain vibration curve of the vibrating diaphragm. According to the time domain hybrid finite element fast algorithm, the tree-cotree technology and the Lanczos modal order reduction technology are introduced, the method can be used for analyzing the time domain electromagnetic force coupling simulation problem of the complex micro loudspeaker, and the method has the advantage of being high in solving efficiency under the condition that the solving precision is met.

The application discloses a borehole time domain electromagnetic far detection device and method, which comprises a ground system, an armored cable, a remote transmission system and a time domain electromagnetic logging instrument; the remote transmission system is connected with the time domain electromagnetic logging instrument and is arranged in a borehole; the time domain electromagnetic logging instrument comprises an electronic circuit, an acquisition module and a transmission module which are sequentially arranged from top to bottom; the electronic circuit is connected with the acquisition module and the transmission module for detection; the acquisition module comprises a first array type three-component magnetic receiving sensor, a second array type three-component magnetic receiving sensor, an electric receiving sensor, a third array type three-component magnetic receiving sensor and a fourth array type three-component magnetic receiving sensor which are sequentially arranged from top to bottom; and the transmission module comprises a magnetic source and an electric source which are sequentially arranged from top to bottom. The time domain magnetic and electric hybrid transmission source is adopted for transmission, the array type magnetic sensor is adopted for reception and the electric sensor is adopted for reception, so that the source distance of the detector is shortened, and the borehole far detection is realized.

The invention relates to a superparamagnetic effect identification and separation method in time domain electromagnetic response. The method comprises the following steps: establishing a double-feature superparamagnetic effect identification model based on a fuzzy support vector machine algorithm; fitting the induction-magnetization effect in power law attenuation in the late period by adopting a Levenberg-Marquardt algorithm to obtain a quasi-superparamagnetic effect, and performing de-equalization and whitening on the quasi-superparamagnetic effect and the total effect as input signals of an independent component analysis algorithm; and selecting a mutual information minimization criterion and a natural gradient algorithm to update and optimize the separation matrix according to the statistical characteristics of the transient electromagnetic signal until the convergence condition is met, and finally separating to obtain the induction response and the superparamagnetic effect. According to the method, the superparamagnetic effect can be effectively identified and accurately separated, the inversion interpretation precision of the time-domain electromagnetic method in superficial superparamagnetic layer exploration is effectively improved, and technical support is provided for the development of the time-domain electromagnetic method in refined detection.

The present application relates to a kind of system response correction methods of time domain electromagnetic detection system.The present application method is used to correct the influence of system response on effective signal, removes the influence of inherent characteristics of receiving sensor on effective signal by training adaptive output compensation of the theoretical and actual response of correction device;And by the real measured current waveform and adaptive current sensor compensation, the real transmitting current waveform is obtained, and the response of earth under unit step excitation is obtained according to its inverse function.The present application can remove the influence of system response on time domain electromagnetic detection system signal, improve signal quality, reduce the blind area time of system detection, improve the detection range of system.

The present application relates to the technical field of geophysical electromagnetic data exploration, in particular to a three-dimensional time-frequency domain electromagnetic joint inversion method based on finite element algorithm. The method obtains time domain and frequency domain electromagnetic observation data of the region to be measured respectively; the region to be measured is profiled to obtain a non-structural tetrahedral mesh to construct a time-frequency domain inversion initial model; the time-frequency domain inversion initial model is forward calculated to obtain time domain and frequency domain electromagnetic response data; time domain electromagnetic data residuals and frequency domain electromagnetic data residuals are obtained, and their corresponding adaptive weights are obtained respectively to establish a joint inversion objective function; the joint inversion objective function is iteratively solved to calculate the update vector of the inversion initial model; the inversion initial model is iteratively updated until the set condition is reached, and the time-frequency domain electromagnetic joint inversion is completed. The present application effectively improves the accuracy of three-dimensional inversion through joint inversion, thereby enhancing the practical effect of three-dimensional inversion.

The invention discloses a multi-device cooperative time domain electromagnetic hiding system and method for an electromagnetic positioning system, and aims to improve the concealment and survivability of a positioned device cluster in a complex electromagnetic environment. According to the invention, through signal synthesis among multiple devices and dynamic adjustment of phase and time parameters, an electromagnetic waveform with confusion can be generated so as to mislead a positioning result of an electromagnetic positioning system. The positioned equipment cluster comprises a plurality of mobile carriers, signal sources and a control center, the control center is responsible for commanding the mobile carriers, the signal sources carried by the mobile carriers transmit electromagnetic signals, and the positioning precision of electromagnetic positioning equipment can be effectively reduced by changing parameters such as phase, amplitude and forwarding delay of the signals. According to the technical scheme, the electromagnetic hiding effect of the positioned device cluster on the electromagnetic positioning system can be remarkably improved, and the method has important application value in the aspects of device safety and attack defense.

The invention relates to a vibration calculation method of an embedded spoke type permanent magnet motor for flywheel energy storage, which comprises the following steps of: respectively establishing equivalent magnetic network models for stator, air gap and rotor areas of the embedded spoke type permanent magnet motor according to the distribution characteristics of magnetic lines in the embedded spoke type permanent magnet motor; based on a Maxwell stress tensor method, calculating the time domain electromagnetic force density and the centralized electromagnetic force of the stator tooth surface of the spoke-embedded permanent magnet motor through the air gap magnetic flux density calculated by the equivalent magnetic network model; solving a vibration transfer function set of the embedded spoke type permanent magnet motor containing the electromagnetic force modulation effect through an analytical method; and finally, calculating the vibration response of the embedded spoke type permanent magnet motor through the combination of the vibration transfer function set and the centralized electromagnetic force in the frequency domain. According to the method, rapid calculation of the vibration response of the embedded spoke type permanent magnet motor is realized.

The invention provides a transient electromagnetic forward modeling method and device containing an induced polarization effect and a storage medium, and relates to the technical field of electrical digital data processing. The method comprises the following steps: constructing a geologic model; applying an excitation field to the geologic model; determining a cut-off time T1 and a cut-off frequency F1 within a time period from T0 to T2 after the excitation field disappears; obtaining time domain electromagnetic data in a time period from T0 to T1 and frequency domain electromagnetic data in a frequency section from Fmin to F1 based on the first mathematical model and a Maxwell equation set; based on the frequency domain electromagnetic data, the time domain electromagnetic data and a second mathematical model, obtaining response electromagnetic data containing an induced polarization effect at any position in the geologic model within a time period from T0 to T2; according to the transient electromagnetic forward modeling method containing the induced polarization effect, the problem that in the prior art, due to the fact that calculation is difficult in the frequency domain, transient electromagnetic three-dimensional rapid forward modeling containing the induced polarization effect is difficult to achieve can be solved.

The invention discloses a large-loop source transient electromagnetic fast forward modeling method based on GPU acceleration, and the method comprises the steps: firstly carrying out the parameter setting and preprocessing, then segmenting a large-loop source into a north edge, a south edge, an east edge and a west edge, and carrying out the discretization of each edge into N electric dipoles through employing a Gaussian-Legendre quadrature method; constructing a loop edge-frequency sampling point-wave number-time channel four-level parallel computing architecture to compute each edge; wherein each side in the GPU adopts a frequency-wave number-time parallel computing architecture, and parallel computing is performed on each electric dipole, so that the frequency-time domain response of each electric dipole is obtained; and finally, after calculation of all the electric dipoles is completed, vectorization calculation and summation are carried out on transient responses of all the electric dipoles on the GPU, a large loop source total field response at the current observation point is obtained, and a final time domain electromagnetic field value is output. According to the method, the advantages of large-scale parallel computing of the GPU can be efficiently and fully played, so that the computing efficiency is effectively improved on the premise of ensuring the computing precision.

The application relates to a time-domain electromagnetic space-time double-gradient resistivity imaging method, device, computer equipment, computer readable storage medium and computer program product. The method comprises the following steps: acquiring transient electromagnetic method data to be processed, wherein the data to be processed comprises electric field component attenuation data of m*n measuring points, and a far zone condition is met between the m*n measuring points and a long wire source; acquiring a spatial gradient expression of apparent resistivity and a time gradient expression of apparent resistivity; performing finite difference fusion based on the spatial gradient expression of resistivity and the time gradient expression of apparent resistivity to obtain a gradient resistivity expression; and generating a resistivity distribution map of a subsurface medium according to the data to be processed and the gradient resistivity expression. The abnormal body boundary description accuracy is improved.

The application discloses a data storage method of time domain electromagnetic signals, comprising the following steps: initiating a transfer transaction from a sending account to a receiving account of a block chain transfer transaction, and placing time domain electromagnetic signal data in a parameter input area; performing a transfer operation to complete the on-chain storage of the time domain electromagnetic signal data. The data storage method of the time domain electromagnetic signals on the block chain chain based on the transfer transaction writes the time domain electromagnetic signal data in the form of the transfer transaction into the block chain through the decentralized idea, thereby avoiding the phenomenon of "super-privilege user tampering" of the centralized storage, and improving the security and reliability of the time domain electromagnetic signal data storage.

The present disclosure belongs to the technical field of geophysical exploration signal processing, and particularly relates to a time-domain electromagnetic polarization parameter extraction method based on an improved physical guided neural network, comprising generating a sample set; constructing a loss function constrained by physical information; using the sample set, minimizing the loss function, and using an improved whale optimization algorithm to optimize the hyperparameters of the neural network to obtain a time-domain electromagnetic polarization parameter neural network inversion model; wherein the improved whale optimization algorithm introduces chaotic mapping to generate an initial population in the initialization stage, introduces a dynamic opposite learning strategy in the bubble net attack stage, and introduces an adaptive t-distribution disturbance variation in the search prey stage; obtaining measured data and preprocessing the same, using the time-domain electromagnetic polarization parameter neural network inversion model to process the preprocessed measured data, and extracting polarization medium model parameters. The method improves the accuracy and robustness of polarization parameter interpretation.

This invention discloses a rapid detection method for shallow ore bodies based on the fusion of time-domain electromagnetic method and AI, comprising: acquiring time-domain electromagnetic induction data of each unit within the detection area and extracting voltage decay curves to obtain preliminary ore body electromagnetic response characteristics; analyzing electromagnetic signal abrupt change points between adjacent units based on the preliminary ore body electromagnetic response characteristics to determine ore body boundary location information; processing the ore body boundary location information using a clustering algorithm to generate a three-dimensional geological structure model; obtaining clear electromagnetic response signals based on the three-dimensional geological structure model; identifying aliased electromagnetic response components in overlapping areas based on the clear electromagnetic response signals to determine existing ore body structures; when multiple ore body structures are determined to exist, acquiring response components and decomposing them into independent electromagnetic signal components using a signal separation algorithm to obtain time constant difference characteristics of each ore body; constructing location modeling relationships based on the time constant difference characteristics and generating a spatial distribution map of overlapping ore bodies.

The present invention discloses a time-varying electromagnetic field implicit finite volume calculation method, equipment, medium and system, which belongs to the field of numerical solution of time-domain electromagnetic finite volume method, including: iterative solution of Maxwell equations time-varying electromagnetic field by implicit double time step method based on time iterative advancement of conserved electromagnetic field and spatial flux residual; the outer layer of the simulation model is a physical time step cycle until the calculation converges; the inner layer of the simulation model is a virtual time step sub-iteration cycle until the sub-iteration converges; in each virtual time sub-iteration process, each grid block and each grid unit is sequentially subjected to spatial flux calculation and implicit iterative solution calculation, and the next level virtual time sub-iteration step number of conserved electromagnetic field values ​​is updated. The present invention can solve large-scale electromagnetic scattering problems of arbitrarily complex shapes, high-frequency electrical large-size targets.

The invention discloses a machine learning-based time domain electromagnetic transmitter performance evaluation method, which comprises the following steps of: extracting key parameter data from an operation log of an electromagnetic transmitter, and preprocessing the key parameter data to obtain a standardized data set; constructing a multi-dimensional feature set according to the standardized data set, and dividing the multi-dimensional feature set into an input feature vector and a target output variable; training a regression prediction model according to the multi-dimensional feature set to obtain a trained performance prediction model; and predicting the to-be-evaluated data according to the performance prediction model to obtain a prediction result, and performing feature importance analysis on the prediction result to obtain an analysis result. According to the method, intelligent evaluation of transmitter performance under the conditions of multiple devices and multiple working conditions can be achieved, the limitation of a traditional experience or single-index method is overcome, the precision and interpretability of performance analysis are improved, and a new data driving scheme is provided for instrument operation optimization and abnormality diagnosis.

The application discloses a GPU acceleration-based large-loop-source transient electromagnetic fast forward method, which comprises the following steps: firstly, parameter setting and pretreatment are performed; then, the large-loop-source is divided into four edges of north, south, east and west, and each edge is discretized into N electric dipoles by using a Gauss-Legendre quadrature method; a four-level parallel computing architecture of loop edge-frequency sampling point-wave number-time channel is constructed to calculate each edge; wherein, the parallel computing architecture of frequency-wave number-time is used in each edge in the GPU to perform parallel calculation on each electric dipole, so that the frequency-time domain response of each electric dipole is obtained; finally, after the calculation of all electric dipoles is completed, the vectorization calculation and summation of the transient responses of all electric dipoles are performed on the GPU to obtain the total field response of the large-loop-source at the current observation point, and the final time domain electromagnetic field value is output. The above method can efficiently and fully exert the large-scale parallel computing advantage of the GPU, so that the calculation efficiency is effectively improved under the premise of ensuring the calculation precision.

The invention discloses a shallow ore body rapid detection method based on fusion of a time domain electromagnetic method and AI, and the method comprises the steps: obtaining time domain electromagnetic induction data of each unit in a detection region, and extracting a voltage attenuation curve to obtain initial ore body electromagnetic response features; analyzing electromagnetic signal mutation points between adjacent units based on the preliminary ore body electromagnetic response characteristics, and determining ore body boundary position information; processing the ore body boundary position information by adopting a clustering algorithm to generate a three-dimensional geologic structure model; obtaining a clear electromagnetic response signal based on the three-dimensional geological structure model; on the basis of the clear electromagnetic response signals, aliasing electromagnetic response components in the overlapping area are identified, and an existing ore body structure is determined; when it is judged that multiple ore body structures exist, response components are obtained and decomposed into independent electromagnetic signal components through a signal separation algorithm, and time constant difference characteristics of all ore bodies are obtained; and constructing a position modeling relation based on the time constant difference characteristics, and generating an overlapped ore body spatial distribution diagram.