83 results about "Three dimensional inversion" patented technology

The invention discloses a non-linear conjugate gradient three-dimensional inversion method of a magnetotelluric field, which comprises the following steps of: a, establishing a three-dimensional model, and determining model parameters; b, acquiring impedance response through the model resistivity by using the staggered-grid finite difference of the magnetotelluric field as a forward modeling method; c, comparing the impedance which is acquired by forward modeling with an actually measured impedance to acquire an impedance data deviation; d, judging whether a target function value is enough small or not, finishing iteration when the target function value is small enough, otherwise entering the next step; e, calculating the gradient of a target function, and acquiring a search step-length and a search direction by using a non-linear conjugate gradient method as an inversion method so as to calculate a model modifier, wherein the resistivity calculated by using the impedance data error and current iteration is used as a preprocessing factor in the inversion method; and f, modifying the model, and returning to the step b. The non-linear conjugate gradient three-dimensional inversion method of the magnetotelluric field has the characteristics of high efficiency and high accuracy.

The invention provides a seismic inversion system and a seismic inversion method controlled by geologic facies. The seismic inversion system and the seismic inversion method controlled by the geologic facies comprises a geologic facies establishing module which establishes the geologic facies on a destination layer of a well point position, a seismic facies establishing module which establishes a seismic facies on a plane of the destination layer, a geologic facies constraint establishing module which establishes constraint to the geologic facies, a stochastic inversion module which chooses the number of constraint wells and conducts stochastic inversion controlled by the geologic facies, an analysis module which analyzes inversion data volume and further subdivides the geologic facies, and an output module which iterates repeatedly so as to obtain and output the three-dimensional inversion data volume. The seismic inversion system and the seismic inversion method controlled by the geologic facies solves the problem that due to the fact that geological sedimentary background and regularity in an inversion work area are not considered by the existing inversion method, and inversion results are often horizontally different from geological sedimentary regularity cognition. The seismic inversion system and the seismic inversion method controlled by the geologic facies has the advantages of being capable of accurately depicting oil deposit on the basis of control of the geologic facies.

The invention discloses a method for comprehensive recognition of igneous rocks by employing gravity, magnetism, electromagnetism, and earthquake data. The method includes steps: building templates of different igneous rocks and stratum physical properties; performing forward modeling of gravity, magnetism, electromagnetism, and earthquake, and making clear of response characteristics and rules; performing gravity and magnetism anomaly extraction of the planes of the igneous rocks and weak anomaly enhancement processing, and explaining and dividing distribution zones and range of the planes of the igneous rocks; performing electrical and seismic profile data igneous rock target processing explanation and making clear of the distribution rule of the igneous rocks in the vertical direction; performing combined inversion of gravity, magnetism, electromagnetism, and earthquake data profile, and building a vertical distribution frame of the igneous rocks; developing three-dimensional inversion of the physical properties of the igneous rocks with gravity and magnetism anomalies by employing the profile combined inversion result and drilling data as the constraint condition and with the combination of plane distribution of the igneous rocks; and recognizing the space distribution and lithological and lithofacies distribution rules of the igneous rocks in a comprehensive manner with the combination of physical property combined characteristics of different igneous rocks according to the three-dimensional distribution of the physical properties. According to the method, comprehensive recognition of the igneous rocks is realized by employing various data, and the reliability is higher compared with the recognition of the igneous rock by employing single data.

The invention discloses a method for inversing a wavenumber domain of underground three-dimensional magnetic distribution by magnetic total tensor gradient and an application thereof. The method comprises the steps: known prior information of a study area and observation magnetic total tensor gradient anomality are read; inversion parameters, an inversion depth range, step length, and the like areset; constrained inversion is performed on the observation magnetic total tensor gradient anomality to obtain a magnetic three-dimensional mesh model; forward is performed on the observation magnetictotal tensor gradient anomality to obtain theoretical magnetic total tensor gradient anomality; subtraction is performed on the theoretical anomaly and the observation anomaly to obtain a theoreticalanomaly residual error; and if a minimum norm solution of the anomaly residual error does not satisfy limits of error, the inversion is performed on the anomaly residual error to obtain an anomaly residual error model, the magnetic three-dimensional mesh model is updated, the steps are repeated until the anomaly residual error satisfies requirements, and a final magnetic three-dimensional mesh model is an inversion solution. The method can reflect three-dimensional distribution of magnetic properties (magnetic susceptibility or magnetization intensity) of underground space geologic bodies, and has the application value for improving the three-dimensional inversion efficiency and fine geological simulation of large data bodies.

The invention provides a rapid three-dimensional inversion method for transient electromagnetism based on different constraint conditions. According to the technical scheme of the invention, the data processing and three-dimensional forward calculation method based on the transient electromagnetic moment inversion is adopted. In this way, the problem that the current three-dimensional inversion is large in data volume and complex in forward calculation can be solved. According to the method, firstly, a constraint condition for a time constant vector during the inversion process is constructed. After that, the iteration on the inversion process is conducted based on the optimization algorithm, so that the inversion process is optimized to be close to an actual underground structure. According to the technical scheme of the invention, based on the above method, the rapid three-dimensional inversion process of an underground abnormal target can be realized by an ordinary computer. Therefore, the method has an important application prospect during the interpretation process of real-time transient electromagnetic data.

The invention, which belongs to the evaluation field of geological disasters of water conservancy projects, especially relates to a measurement arrangement system for exploring a riverway karst cave and implementation method thereof, particularly to a multi-dimensional detection method for diagnosing a karst body hidden in a riverway by using a parallel electrical method. The multi-dimensional arrangement apparatus for exploring riverway karst body distribution comprises a test system and a test instrument. The test system consists of two multi-core cable lines arranged at two river embankments in parallel, infinity B electrodes corresponding to the single multi-core cable line, and potential reference N electrodes. 64 electrodes are arranged at each multi-core cable line; and each electrode channel is connected with the ground by a copper rod. The test instrument is a parallel electric method instrument. According to the multi-dimensional detection apparatus and method of exploring a riverway karst cave, because of utilization of the parallel electric method detection technology, the work efficiency of the conventional resistivity method is improved; with an interpreting technique based on combination of a two-dimensional apparent resistivity image and three-dimensional inversion imaging, the karst cave exploration precision is improved and defects in the prior art are overcome.

The application provides a three-dimensional gravity field fast and high-precision inversion algorithm in a spherical coordinate system based on a three-dimensional Gauss-legendre numerical integration formula (3D-GLQ). The core content comprises the following two points that the first point lies in the aspect of three-dimensional gravity field forward modeling in the spherical coordinate system,a nuclear matrix equivalent storage strategy is proposed aiming at the problem of low computational efficiency of the traditional 3D-GLQ three-dimensional gravity field forward modeling algorithm, thenuclear matrix equivalence is combined with the existing adaptive splitting strategy, the computational efficiency is increased by two orders of magnitude under the premise of ensuring that the maximum relative error of the forward modeling calculation is less than 0.1%, and the memory occupancy is greatly reduced at the same time; and the second point lies in that the above fast and high-precision forward modeling algorithm is applied to the gravity field three-dimensional inversion, and the gravity field inversion objective function and depth weighting function in the spherical coordinate system are constructed, so that the computational efficiency and credibility of the large-scale gravity field three-dimensional inversion in the spherical coordinate system are greatly improved.

The invention discloses a polluted soil detection method based on a three-dimensional high-density electrical resistivity method. The polluted soil detection method comprises the following steps that a facet-shaped detection unit zone is arranged in a polluted region, measuring lines are arranged in the facet-shaped detection unit zone and include a plurality of detection electrodes distributed in the facet-shaped detection unit zone in a uniformly-spaced dot matrix mode and single or multiple serially connected detection cables sequentially connected with the detection electrodes; electrical resistivity data of the detection electrodes of the measuring lines are acquired; three-dimensional inversion processing is conducted on the electrical resistivity data to obtain corresponding electrical resistivity profile data; the polluted region resistivity data measured by the measuring lines and electrical resistivity data of a non-polluted region are compared, and meanwhile the pollutant distribution in the polluted region is judged by combining with the electrical resistivity profile data. The polluted soil detection method has the advantages the time condition, power supply condition, electrode grounding condition and the like for data acquisition have consistency, and accordingly obtained electrical resistivity values have comparability and contrast.

The invention discloses a focusing sounding type three-dimensional induced polarization advanced detection system and method loaded on a TBM. The advanced detection system comprises an electrode system, an electrode telescoping device, a multi-path rotation device, a cable system and an induced polarization detection instrument; the electrode system is arranged on a TBM cutterhead, a shield, a trolley and a tunnel side wall, wherein electrodes on the cutterhead, the shield and the trolley are loaded on the TBM through the electrode telescoping device; and the induced polarization detection instrument controls the electrode telescoping device to realize automatic telescoping control over the electrodes, telescopic state detection and contact state detection, and the multi-path rotation device is adopted to realize uninterrupted rotation of measurement signals, hydraulic oil paths and cutterhead water paths and realize load integration of the induced polarization detection system and theTBM. The system and the method have the advantages that the problem of induced polarization advanced detection of a source of water and mud inrush disaster of a TBM construction tunnel can be solved,the automation degree of measurement is high, the safety, the rapidness and the reliability are achieved, three-dimensional inversion imaging can be carried out on a water-contained body in the frontof a tunnel face, and therefore the safe and efficient construction requirements of the TBM construction tunnel are met.

A three-dimensional inversion method for airborne transient electromagnetic conductivity based on the Bornn iterative method, involving electromagnetic exploration. Provide a set of three-dimensional conductivity inversion algorithm based on the frequency domain to improve the calculation speed of the inversion; the algorithm uses the fitting of the secondary field value to perform the inversion calculation, thereby removing the influence of the ground response of the direct wave on the signal and improving the inversion calculation speed. Acting precision. Firstly, the time-domain signal of the total field is obtained from the transient electromagnetic system, and then the time-domain response of the corresponding theoretical uniform half-space is calculated to extract the secondary field signal required for the inversion. First use the electric field integral equation to solve the electric field distribution in the calculation area, and calculate the secondary field value through the definition of the contrast function. Then use the Born iterative method to establish the inversion process, fit the secondary field values ​​in the frequency domain, and iteratively solve the contrast function, so as to obtain the distribution of the medium in the underground calculation area, and realize the three-dimensional inversion of the conductivity of the transient electromagnetic system. play.

The invention relates to a method for a wire power source multichannel receiving ground nuclear magnetic resonance and a system thereof. The method comprises the following the steps that: 1) a transmitter is used to supply alternating current with Larmor frenquency to the underground through two embedded emission electrodes, and the distance of the electrodes and the current are changed to generate different excited magnetic fields; 2) the excited magnetic fields are broken, and signals which are in exponential decaying laws are received through a high-sensitive probe in the process that excited protons are precessed into a normal balanced state; 3) the signals are transmitted to a receiver through signal transmission lines, and received signals of the receiver are sent to a signal processing and imaging system to carry out treatments to the signals, such as pretreatment, denoising, inversion, imaging, etc. The method for the wire power source multichannel receiving ground nuclear magnetic resonance and the system thereof have the advantages that: small volume of meticulous exploration underground can be provided; the method and the system thereof have the characteristics of high efficiency of multichannel receiving work, large explorating space, and is favorable for three-dimensional inversion imaging and so on; the wire power source are excited so that construction is flexible and simple; the efficiency of the multichannel receiving work is high; and an aim of exploration is realized through changing emission current.

The invention provides a method for identifying a mountain front deep part construction layer through utilizing gravity interface inversion and is a gravity density interface inversion method for difficulty in depicting deep part stratum structure attributes during mountain front oil gas target exploration. The method comprises steps that gravity potential field data pre-processing is carried outto form basic data before data inversion; different target area lithology and stratum physical property rules are analyzed, and physical property base of density interface inversion is made to be clear; a variable density physical property model is established, and an inversion physical property model is provided; gravity-earthquake combined step by step fine delamination is carried out, and superposition field source separation precision is enhanced; drilling well-earthquake constraints are added to the space domain, and interface constraints are provided by skeleton profile combined inversion; gravity density interface three-dimensional inversion is developed, and the geological structure attributes of the mountain front deep part construction layer is made to be clear. The method is advantaged in that distribution characteristics of the mountain front deep part construction layer is effectively depicted, geological basis is provided for mountain front deep oil gas exploration, andoil gas exploration risks are reduced.

The invention relates to a method for building a horizontal layered soil electrical model, which comprises the steps of building a composite layered soil model according to magnetotelluric three-dimensional inversion soil resistivity data, selecting the earth surface potential to act as a fitting index, fitting the earth surface potential of a composite model acquired by using a finite element method and the earth surface potential of a horizontal model deduced by using a Green function method by adopting a least square method so as to acquire an objective function of the horizontal layered soil electrical model, carrying out optimization solving on the objective function by adopting a genetic algorithm, acquiring a soil resistivity parameter and the corresponding layer depth of each layer in the horizontal layered soil electrical model, and thus realizing conversion from composite layering to horizontal layering of the model. The method provided by the invention reduces the complexity of the model, improves the reliability of the model, and can provide more accurate soil resistivity model data for analysis of an earth electrode.

The invention discloses a sea magneto-telluric field non-linear conjugate gradient three-dimensional parallel inversion method, comprising the steps of: adding a sea layer with a fixed resistivity value into an initial model, and deriving the interpolation expression and computational process of the magneto-telluric field value with sea layer model each node to realize directive calculation; and deriving the computational expression and solution process of the model modifiers (a computational process mainly comprises variable parameters such as an objective function, an objective function gradient, a searching direction, a searching step size, etc.) with sea layer model each node to realize directive calculation. The method is suitable for magneto-telluric field non-linear conjugate gradient three-dimensional inversion of sea or other water areas.

The invention discloses a three-dimensional while-drilling induced polarization rock precise imaging device and method. The method is characterized in that a drilling bit in drilling serves as a power supply electrode A which is placed into a hole; a working face is provided with a three-dimensional array measuring electrode; an electrode N and an electrode B are arranged behind the working face; and while-drilling real-time data measurement and real-time three-dimensional inversion imaging are realized, and three-dimensional refined real-time imaging detection of unfavorable geology before the working face is realized. The method can realize while-drilling real-time measurement without interfering drilling construction; the power supply electrode is arranged in the drill hole, so that interference from metal bodies and the like behind a tunnel can be prevented; the drilling bit is served as the power supply electrode creatively; and the advance drilling machine is small in change, and in-hole induced polarization detecting precision and work efficiency are improved.

The present invention relates to a geological integrated identification method and device based on gravity-magnetic-electric-seismic three-dimensional joint inversion. The method comprises: obtaining a gravity-magnetic-electric-seismic joint inversion result and the physical property statistics of a research area; performing cluster analysis of the gravity-magnetic-electric-seismic joint inversion result, and obtaining the corresponding clustering result; and dividing the inhomogeneity in the clustering result into different stratums according to the physical property statistics of the research area, and performing clustering mean statistics to obtain the geological integrated interpretation result of the gravity-magnetic-electric-seismic joint inversion. Through the process mentioned above, the geological integrated identification method and device based on the gravity-magnetic-electric-seismic three-dimensional joint inversion overcome the influence of the human factors and obtain a reliable interpretation model through adoption of the limited geology and earthquake constraint conditions so as to realize the objective and integrated evaluation of the three-dimensional inversion result and especially have the high applicability for the complex geologic structure condition. Especially, in the condition of large interpretation workload, the geological integrated identification method based on the gravity-magnetic-electric-seismic three-dimensional joint inversion can greatly reduce the workload. The geological integrated identification method and device based on the gravity-magnetic-electric-seismic three-dimensional joint inversion is rapid, efficient and economic, and has wide market prospects and good popularization and usage values.

The invention relates to a gravity gradient data three-dimensional inversion method based on partial smooth constraints. The method comprises the following steps of firstly, calculating second-order partial derivatives of the gravitational potential of a cuboid unit in x, y and z directions of a three-dimensional space to obtain a forward modeling formula of a cuboid element, giving a model, and performing forward modeling calculation according to the forward modeling formula to obtain observation data; manually assigning values to the roughness constraint matrix according to certain prior information to obtain a roughness matrix, and overcoming a skin effect by adopting a weighting mode of integral sensitivity; and finally, writing a target equation containing parameters such as an integral sensitivity matrix, a roughness matrix, regularization parameters, observation data and the like, and solving the target equation by adopting a conjugate gradient method to obtain an inversion result. The technology is suitable for three-dimensional gravity and gravity gradient data inversion under certain geological or geophysical prior conditions.

The present invention discloses an apparent magnetization intensity three-dimensional inversion method based on full-space regularization downward continuation data. A three-dimensional regularization downward continuation operator is used to realize different depths of continuation to obtain a stable magnetic abnormality full space filed, the minimal value of a target equation is solved, the regularization downward continuation operator is substituted and the positive extension field values of different depths are calculated. A full space field inversion target function is constructed. A multi-dimensional search golden section algorithm is used to solve the minimum value of the inversion target function, and the inversion speed is raised. The method has the advantages that the full depth stable downward continuation is realized, and the full space filed with reliable accuracy can be obtained.

The invention provides a PDE-based magnetic measurement data equivalent source reduction-to-the-pole and type conversion method, which comprises the steps: S1, obtaining magnetic field data d0, and establishing a topographic relief curved surface according to topographic height information of an area where the magnetic field data are located; S2, determining a spatial range of mesh generation, andperforming non-uniform multilayer mesh generation; S3, according to the magnetic inclination angle, the magnetic declination angle and the magnetic induction intensity of the geomagnetic field, carrying out PDE three-dimensional inversion calculation with a depth regularization factor, a positive value constraint term and a regularization term on the magnetic field data d0 to obtain a multi-layerequivalent source model of a magnetic anomalous body; and S4, giving different geomagnetic pole parameters, and performing non-linear PDE magnetic field three-dimensional forward modeling calculationby using a multi-layer equivalent source model to obtain one or more of reduction-to-the-pole data, magnetic anomaly total field data, magnetic anomaly component data and magnetic gradient tensor data. The method has the beneficial effects that the underground magnetic anomaly data in a complex environment can be subjected to self-adaptive and high-precision reduction-to-the-pole and data type conversion calculation.

The invention provides a three-dimensional magnetotelluric tipper inversion method using nonlinear conjugate gradients and belongs to the technical field of magnetotelluric sounding. The method comprises steps as follows: S1, inputting data and performing ranking; S2, distributing the parameter of each frequency point to all cpus for tipper inversion and computing; S3, stopping parallel computing and combining the data, and comparing the difference between an invention result and the input date; S4, distributing the data difference and the parameter of each frequency point to all cpus for parallel computing of objective function gradients, computing the search direction p and the step length a and accordingly obtaining a model modification amount; S5, updating a module according to the model modification amount and stopping iteration when an objective function is one hour enough, otherwise, using a new model to repeat the steps from S2 to S5. According to the inversion method, tipper three-dimensional inversion is realized on the basis of tipper data by introducing a nonlinear conjugate gradient method, and the invention result approaches to real underground three-dimensional electrical structure information greatly; the method adopts a parallel computing structure, and the computing efficiency is improved.

The invention relates to a method and a system for three-dimensional detection on an earth electrode field and surrounding soil resistivity. The method comprises the following steps of (1) combining an AMT (audio magnetotelluric) method and an MT (magnetotelluric) method to arrange measuring points in and around the initial range of the earth electrode field; (2) using a high-precision measuring instrument to measure the earth coordinate and elevation of each arranged measuring point, so as to obtain the coordinate data of each measuring point; (3) using an MT instrument system to collect field data of each measuring point according to the requirements of MT and AMT technical specifications, so as to obtain time sequence data of each measuring point; (4) according to the collected coordinate data and time sequence data of each measuring point at field, performing data processing and three-dimensional inversion, so as to obtain the three-dimensional resistivity structure meeting the actual underground geology conditions. The method can be widely applied to the selection of the earth electrode field and the detection of the soil resistivity.

The invention provides a high-precision regional troposphere refractive index three-dimensional inversion method. On the basis of sounding meteorological data accumulated by the national meteorological bureau for a long time, a troposphere refractive index profile model adapting to a test task area in China is established, so that the three-dimensional inversion precision of the troposphere refractive index is improved, and the method comprises the following steps: calculating the troposphere refractive index profile on the basis of the sounding meteorological data; through statistical analysis, establishing a statistical relational expression between each height refractive index of the troposphere of each sounding site and the ground parameters; establishing an inversion model of different height refractive indexes of the troposphere of the sounding site; establishing a high-precision refractive index inversion model by utilizing geographical location information, time information anda ground meteorological information driving model of a target site and adopting a double-harmonic spline interpolation method. Therefore, the inversion precision can be further improved, the inversion effect on any point in the area is better, the number of driving elements is small, and the method has higher generalization performance.

The invention provides an upper continuation and lower continuation method for a multilayer equivalent source of magnetic survey data. The method comprises the steps of S1, inputting magnetic field data, and establishing a topographic relief curved surface according to topographic height information of an area where a relief observation curved surface is located; S2, determining a spatial range ofmesh generation, performing continuous structured non-uniform mesh generation, and further determining an inversion model solving space; S3, according to the magnetic inclination angle, the magneticdeclination angle and the magnetic induction intensity of the geomagnetic field, performing integral equation three-dimensional inversion calculation with a depth regularization factor, a positive value constraint term and a regularization term on the magnetic field data based on an inversion model solving space to obtain a multi-layer equivalent source model of a magnetic anomalous body; and S4,performing magnetic field forward modeling calculation based on an integral equation by utilizing the multilayer equivalent source model to obtain magnetic field data after upper continuation and / or lower continuation generated by the magnetic anomalous body. The method has the beneficial effects that the required upper continuation or lower continuation data can be adaptively, quickly, efficiently and accurately generated for the magnetic anomalous body.

The invention relates to a method for obtaining three-dimensional distribution positions of uranium anomalies in a sand body through a well-seismic combination method. The method comprises the following steps: reading three-dimensional earthquake, drilling and logging data; calibrating a uranium abnormal target layer section and a fine well earthquake; carrying out well data rock physical sensitivity analysis of a uranium anomaly target layer section, fine structure interpretation and inversion initial model establishment of the uranium anomaly target layer section, high-resolution three-dimensional well seismic sand body inversion and precision evaluation of the uranium anomaly target layer section, rapidly obtaining the corresponding relation between uranium anomaly and sand body space by superposing a gamma logging curve and a three-dimensional inversion data body and the like. The specific distribution of uranium anomalies in a sand body is solved, and a reliable drill distributionbasis is provided for detailed exploration. Compared with an analysis method only based on drilling data, the method provided by the invention can obtain a well-to-seismic combination method including rock physical sensitivity analysis and precision evaluation, and has better pertinence, reliability and precision guarantee. A large amount of drilling workload is greatly reduced, the working efficiency of detailed survey is improved, and the field working cost is greatly reduced.

The invention discloses a current element three-dimensional inversion method based on an algebraic elastic network regularization method, which comprises two functional form formulas of expressing a current element column vector by I, substituting an algebraic elastic network regularization operator into the formula 1 to obtain a formula. Beneficial effects of the method are that: according to thealgebraic elastic network regularization method, the direction and position calculation result of current element distribution is basically kept consistent with actual distribution, the ill-conditioned problem embodied in the calculation process is solved, and the algebraic elastic network regularization method in the algebraic regularization method has high accuracy and stability in the aspect of current element three-dimensional imaging calculation.

The invention discloses a ground nuclear magnetic resonance three-dimensional inversion method. First, dimensionality reduction is performed on a three-dimensional forward model through straightening transformation, and the three-dimensional inversion problem is abstracted to a matrix equation solving model; a constrained GMERR algorithm is provided for solving the SNMR signal three-dimensional inversion problem; an intermediate solution of the water content is obtained through the GMERR algorithm first, the value range of the intermediate solution is then forcibly converted to range from 0 to 100% through constraint conditions, and the optimal solution is obtained through multiple iterative operations. The ground nuclear magnetic resonance three-dimensional inversion method is high in inversion precision at different signal-to-noise ratios, good in stability and high in anti-noise capacity.

The invention discloses a multi-channel transient electromagnetic inversion method based on a footprint technology and a transient electromagnetic survey device, and the method can greatly reduce thethree-dimensional inversion space corresponding to a single data point, improves the multi-channel transient electromagnetic three-dimensional inversion efficiency, and comprises the following steps:S1, transmitting a pseudo-random coded current; S2, receiving a pseudo-random response containing geological information; S3, obtaining a ground pulse response; S4, loading an inversion area and meshgeneration model parameters; S5, constructing a minimum parameter functional; S6, constructing a footprint weighting matrix, and constructing a footprint weighting matrix; and S7, obtaining the distribution condition of the geoelectric structure in the inversion area. The device comprises a transmitting probe, a receiver and a working platform.

The invention provides a PDE-based magnetic measurement data equivalent source upper continuation and lower continuation method, which comprises the steps of S1, obtaining magnetic field data d0, andestablishing a topographic relief curved surface according to topographic height information of an area where the magnetic field data is located; s2, determining a spatial range of mesh generation, and performing non-uniform multilayer mesh generation; s3, according to the magnetic inclination angle, the magnetic declination angle and the magnetic induction intensity of the geomagnetic field, carrying out PDE three-dimensional inversion calculation with a depth regularization factor, a positive value constraint term and a regularization term on the magnetic field data d0 to obtain a multi-layer equivalent source model of the magnetic anomalous body; and S4, performing PDE-based magnetic field forward modeling calculation by using the multilayer equivalent source model to obtain upper continuation and / or lower continuation magnetic field data generated by the magnetic anomalous body. The invention has the beneficial effect that self-adaptive and high-precision upper continuation or lower continuation can be carried out on underground magnetic anomaly data in a complex environment.

The invention discloses a transient electromagnetic and induced polarization effect decoupling method. The method is characterized by dividing a receiving signal into a direct current component and analternating current component and taking the alternating current component as a mixing signal; using an independent component analysis algorithm to realize transient electromagnetic and induced polarization effect decoupling; extracting a simple transient electromagnetic signal and an induced polarization signal; and then, carrying out three-dimensional inversion and interpretation. By using themethod of the invention, low-resistance geological abnormal information in geological mining detection can be effectively extracted, and an induced polarization effect is suppressed, and the comprehensive interpretation effect of transient electromagnetic and induced polarization in ore body exploration can be effectively increased.