CN113656750A - Calculation method of magnetic induction intensity of ferromagnetic medium based on spatial wavenumber domain - Google Patents
Calculation method of magnetic induction intensity of ferromagnetic medium based on spatial wavenumber domain Download PDFInfo
- Publication number
- CN113656750A CN113656750A CN202111218085.7A CN202111218085A CN113656750A CN 113656750 A CN113656750 A CN 113656750A CN 202111218085 A CN202111218085 A CN 202111218085A CN 113656750 A CN113656750 A CN 113656750A
- Authority
- CN
- China
- Prior art keywords
- domain
- magnetic
- field
- spatial
- space
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000005291 magnetic effect Effects 0.000 title claims abstract description 300
- 238000004364 calculation method Methods 0.000 title claims abstract description 39
- 230000006698 induction Effects 0.000 title claims abstract description 34
- 230000005294 ferromagnetic effect Effects 0.000 title claims abstract description 15
- 230000002159 abnormal effect Effects 0.000 claims abstract description 131
- 230000005415 magnetization Effects 0.000 claims abstract description 24
- 230000002547 anomalous effect Effects 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims description 44
- 239000006185 dispersion Substances 0.000 claims description 8
- 238000004800 variational method Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 8
- 238000013508 migration Methods 0.000 abstract 1
- 230000005012 migration Effects 0.000 abstract 1
- 238000004590 computer program Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000005389 magnetism Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000006257 total synthesis reaction Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Operations Research (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
基于空间波数域的强磁介质的磁感应强度计算方法，剖分包含异常体的初始棱柱体模型，得到一系列的节点，对每个节点进行磁化率赋值，磁化率为标量；设定高斯参数，计算离散偏移波数，计算空间域背景场磁场强度，获得磁化强度模型；将空间域异常场磁位和磁化强度满足的方程转为空间波数域一维常微分方程；结合设定的空间波数域异常场磁位的边界条件，将空间波数域异常场磁位满足的边值问题模型等价转化为变分问题模型并求解，得到空间波数域异常场磁位、空间波数域异常场磁场强度后将其转换为空间域异常场磁位及空间域异常场磁场强度，根据输出的空间域异常场磁场强度求得空间域磁感应强度。本发明能更准确的对强磁性介质进行磁场数值模拟。
Based on the calculation method of magnetic induction intensity of ferromagnetic medium in the space wavenumber domain, the initial prism model containing the abnormal body is divided, a series of nodes are obtained, and the magnetic susceptibility is assigned to each node. Calculate the discrete migration wavenumber, calculate the magnetic field strength of the background field in the space domain, and obtain the magnetization model; convert the equation satisfying the magnetic potential and magnetization of the anomalous field in the space domain into a one-dimensional ordinary differential equation in the space wavenumber domain; combine the set space wavenumber domain The boundary condition of the magnetic potential of the anomalous field, the boundary value problem model satisfied by the magnetic potential of the anomalous field in the spatial wavenumber domain is equivalently transformed into a variational problem model and solved, and the magnetic potential of the anomalous field in the spatial wavenumber domain and the magnetic field strength of the anomalous field in the spatial wavenumber domain are obtained. It is converted into the magnetic potential of the anomalous field in the space domain and the magnetic field intensity of the anomalous field in the space domain, and the magnetic induction intensity in the space domain is obtained according to the output magnetic field intensity of the anomalous field in the space domain. The invention can carry out numerical simulation of the magnetic field for the ferromagnetic medium more accurately.
Description
Technical Field
The invention belongs to the technical field of strong magnet numerical simulation, and particularly relates to a magnetic induction intensity calculation method of a strong magnetic medium based on a space wavenumber domain.
Background
When the magnetic susceptibility value of the medium is greater than 0.1 SI, it is generally considered to be a ferromagnetic medium, and the self-demagnetizing field in the magnetic field formed by it is not negligible. At present, the magnetic induction intensity of the magnetic field of the ferromagnetic medium is usually calculated by neglecting a self-demagnetizing field, namely, the calculation is approximate as a weak magnetic condition, so that the magnetic induction intensity obtained by numerical simulation has larger deviation with an actual magnetic induction intensity value, and the magnetic measurement data is not favorably processed and explained.
For the calculation of the magnetic field of a ferromagnetic medium, research methods of related scholars are mainly divided into space-domain and frequency-domain methods. Among The spatial domain methods, The literature (volume a. The application of electronic computers to The calculation of effective magnetization [ J ]. geographic magnetization 1963.11(1): 51-58.) proposes to represent The effective magnetization of each volume element inside The magnet by a vector series and use successive approximation to solve, but when The magnetic susceptibility is high, The vector series converges slowly or even diverges, and The calculation efficiency is low. The literature (pursuis, m.b. j., and j.p. fill, a new iterative method for computing the magnetic field at high magnetic susceptibilities: geophisics, 2005.70: 53-62.) calculates the magnetic field of high susceptibility magnetic bodies on an iterative basis by means of a piecewise model defined by spherical voxels of uniform arbitrary diameter, but the high susceptibility calculation results are subject to large errors. In the frequency domain method, a new Gaussian FFT method is proposed based on a shift sampling method in documents (Wu, L.Y., and G.Tian, High-precision Fourier transforming of potential fields: Geophysics, 2014.79, No. 5, G59-G68, doi: 10.1190/geo 2014-0039.1), and is successfully applied to numerical simulation of gravity and magnetism, but is still only applicable to the case of low magnetic susceptibility. The literature (Likun, Daishikun, Chenganui, etc., three-dimensional numerical simulation of magnetic anomaly field integral solution in space and wave number mixed domain, geophysical report 2019.62 (11)) utilizes the characteristic that magnetic potential three-dimensional space domain integral is convolution to perform two-dimensional Fourier transform along the horizontal direction, and converts the three-dimensional integral problem satisfied by the magnetic potential in the space domain into the mutually independent vertical one-dimensional integral problem among different wave numbers, but is only suitable for the condition of low magnetic susceptibility.
Based on the current research situation, it is known that strong magnetism is rarely considered in the current magnetic induction intensity calculation method, and some calculation methods considering strong magnetism have low efficiency and low precision. Therefore, it is imperative to provide a high-efficiency and high-precision magnetic induction intensity calculation method suitable for strong magnets.
Disclosure of Invention
Aiming at the problems that the self-demagnetization effect is less considered in the calculation of the magnetic induction intensity of the strong magnet at present, and the research of the self-demagnetization effect is less considered, the precision is insufficient and the efficiency is low, the invention aims to provide a magnetic induction intensity calculation method of the strong magnetic medium based on the space wave number domain.
In order to achieve the technical purpose, the technical scheme provided by the invention is as follows:
in one aspect, the present invention provides a method for calculating magnetic induction intensity of a ferromagnetic medium based on a spatial wavenumber domain, including:
acquiring a three-dimensional target area containing an abnormal body, and establishing an initial prism model containing the three-dimensional target area;
for the initial prism model edgex、y、zThe directions are respectively divided at equal intervals inx、y、zObtaining a series of nodes in the direction; assigning a magnetic susceptibility to each node based on the magnetic susceptibility distribution data, the magnetic susceptibility being a scalar quantity, usingχRepresenting, at nodes of the anomaly part, each node being assigned according to a susceptibility value of the anomaly; the magnetic susceptibility of the nodes of the abnormal part is 0;
according to prism models and givenxDirection andygaussian parameter in direction, calculatingxAndydirection dispersion offset wavenumber;
calculating the strength of the main earth magnetic field at each node according to the main earth magnetic field model, and taking the strength of the main earth magnetic field as the strength of a spatial domain background field;
obtaining a calculation model of the magnetization intensity according to the intensity of the spatial domain background field magnetic field and the intensity of the spatial domain abnormal field magnetic field;
converting a three-dimensional Laplace equation which meets the magnetic potential and the magnetization intensity of the abnormal field in the spatial domain into a one-dimensional ordinary differential equation in the spatial wave number domain by utilizing two-dimensional Fourier transform;
based on a space wave number domain one-dimensional ordinary differential equation and in combination with a set boundary condition which needs to be met by the magnetic potential of the abnormal field of the space wave number domain, converting an edge value problem model which is met by the magnetic potential of the abnormal field of the space wave number domain into an equivalent variational problem model;
obtaining the magnetic potential of the space wave number domain abnormal field by solving a variational problem model;
obtaining the magnetic field intensity of the spatial wave number domain abnormal field based on the magnetic potential of the spatial wave number domain abnormal field, and converting the magnetic potential of the spatial wave number domain abnormal field and the magnetic field intensity of the spatial wave number domain abnormal field into the magnetic potential of the spatial wave number domain abnormal field and the magnetic field intensity of the spatial wave number domain abnormal field through inverse Fourier transform;
the total magnetic field intensity of the space domain is the sum of the magnetic field intensity of the abnormal field of the space domain and the magnetic field intensity of the background field of the space domain, and the total magnetic field intensity of the space domain is tightened by using a tightening operator to obtain the total magnetic field intensity of the space domain after tightening;
judging whether iterative convergence conditions are met or not, and if yes, outputting a space domain abnormal field magnetic potential and a space domain abnormal field magnetic field strength corresponding to the tightened space domain total magnetic field strength obtained through current calculation;
and solving according to the output spatial domain abnormal field magnetic field intensity to obtain spatial domain magnetic induction intensity.
In another aspect, the present invention provides a magnetic induction calculation apparatus for a ferromagnetic medium based on a spatial wave number domain, including:
the system comprises a first module, a second module and a third module, wherein the first module is used for acquiring a three-dimensional target area containing an abnormal body and establishing an initial prism model containing the three-dimensional target area;
a second module for aligning the initial prism model edgex、y、zThe directions are respectively divided at equal intervals inx、y、zObtaining a series of nodes in the direction; assigning a magnetic susceptibility to each node based on the magnetic susceptibility distribution data, the magnetic susceptibility being a scalar quantity, usingχIs shown therein in the abnormal body partNodes, each node being assigned according to a susceptibility value of the anomalous body; the magnetic susceptibility of the nodes of the abnormal part is 0;
a third module for determining the shape of the prism model and the givenxDirection andygaussian parameter in direction, calculatingxAndydirection dispersion offset wavenumber;
the fourth module is used for calculating the strength of the main earth magnetic field at each node according to the main earth magnetic field model and taking the strength of the main earth magnetic field as the strength of a spatial domain background field;
the fifth module is used for obtaining a calculation model of the magnetization intensity according to the ambient field magnetic field intensity of the spatial domain and the abnormal field magnetic field intensity of the spatial domain;
the sixth module is used for converting a three-dimensional Laplace equation which meets the magnetic potential and the magnetization intensity of the abnormal field in the spatial domain into a one-dimensional ordinary differential equation in the spatial wavenumber domain by utilizing two-dimensional Fourier transform;
the seventh module is used for converting an edge value problem model which is satisfied by the magnetic potential of the abnormal field in the space wave number domain into an equivalent variational problem model based on a one-dimensional ordinary differential equation in the space wave number domain and by combining the set boundary condition which is required to be satisfied by the magnetic potential of the abnormal field in the space wave number domain;
the eighth module is used for obtaining the magnetic potential of the space wave number domain abnormal field by solving the variational problem model;
the ninth module is used for solving the magnetic field intensity of the spatial wave number domain abnormal field based on the magnetic potential of the spatial wave number domain abnormal field, and converting the magnetic potential of the spatial wave number domain abnormal field and the magnetic field intensity of the spatial wave number domain abnormal field into the magnetic potential of the spatial wave number domain abnormal field and the magnetic field intensity of the spatial wave number domain abnormal field through inverse Fourier transform;
a tenth module, configured to apply a tightening operator to tighten the total magnetic field strength of the space domain to obtain a tightened total magnetic field strength of the space domain, where the total magnetic field strength of the space domain is a sum of an abnormal magnetic field strength of the space domain and a background magnetic field strength of the space domain;
an eleventh module, configured to determine whether an iterative convergence condition is currently satisfied, and if so, output a spatial domain abnormal field magnetic potential and a spatial domain abnormal field magnetic field strength corresponding to the tightened spatial domain total magnetic field strength obtained through current calculation;
and the twelfth module is used for solving according to the output space domain abnormal field magnetic field intensity to obtain the space domain magnetic induction intensity.
In another aspect, the present invention provides a computer device comprising a memory and a processor, the memory storing a computer program, the processor implementing the following steps when executing the computer program:
acquiring a three-dimensional target area containing an abnormal body, and establishing an initial prism model containing the three-dimensional target area;
for the initial prism model edgex、y、zThe directions are respectively divided at equal intervals inx、y、zObtaining a series of nodes in the direction; assigning a magnetic susceptibility to each node based on the magnetic susceptibility distribution data, the magnetic susceptibility being a scalar quantity, usingχRepresenting, at nodes of the anomaly part, each node being assigned according to a susceptibility value of the anomaly; the magnetic susceptibility of the nodes of the abnormal part is 0;
according to prism models and givenxDirection andygaussian parameter in direction, calculatingxAndydirection dispersion offset wavenumber;
calculating the strength of the main earth magnetic field at each node according to the main earth magnetic field model, and taking the strength of the main earth magnetic field as the strength of a spatial domain background field;
obtaining a calculation model of the magnetization intensity according to the intensity of the spatial domain background field magnetic field and the intensity of the spatial domain abnormal field magnetic field;
converting a three-dimensional Laplace equation which meets the magnetic potential and the magnetization intensity of the abnormal field in the spatial domain into a one-dimensional ordinary differential equation in the spatial wave number domain by utilizing two-dimensional Fourier transform;
based on a space wave number domain one-dimensional ordinary differential equation and in combination with a set boundary condition which needs to be met by the magnetic potential of the abnormal field of the space wave number domain, converting an edge value problem model which is met by the magnetic potential of the abnormal field of the space wave number domain into an equivalent variational problem model;
obtaining the magnetic potential of the space wave number domain abnormal field by solving a variational problem model;
obtaining the magnetic field intensity of the spatial wave number domain abnormal field based on the magnetic potential of the spatial wave number domain abnormal field, and converting the magnetic potential of the spatial wave number domain abnormal field and the magnetic field intensity of the spatial wave number domain abnormal field into the magnetic potential of the spatial wave number domain abnormal field and the magnetic field intensity of the spatial wave number domain abnormal field through inverse Fourier transform;
the total magnetic field intensity of the space domain is the sum of the magnetic field intensity of the abnormal field of the space domain and the magnetic field intensity of the background field of the space domain, and the total magnetic field intensity of the space domain is tightened by using a tightening operator to obtain the total magnetic field intensity of the space domain after tightening;
judging whether iterative convergence conditions are met or not, and if yes, outputting a space domain abnormal field magnetic potential and a space domain abnormal field magnetic field strength corresponding to the tightened space domain total magnetic field strength obtained through current calculation;
and solving according to the output spatial domain abnormal field magnetic field intensity to obtain spatial domain magnetic induction intensity.
In yet another aspect, the present invention also provides a computer readable storage medium having a computer program stored thereon, the computer program when executed by a processor implementing the steps of:
acquiring a three-dimensional target area containing an abnormal body, and establishing an initial prism model containing the three-dimensional target area;
for the initial prism model edgex、y、zThe directions are respectively divided at equal intervals inx、y、zObtaining a series of nodes in the direction; assigning a magnetic susceptibility to each node based on the magnetic susceptibility distribution data, the magnetic susceptibility being a scalar quantity, usingχRepresenting, at nodes of the anomaly part, each node being assigned according to a susceptibility value of the anomaly; the magnetic susceptibility of the nodes of the abnormal part is 0;
according to prism models and givenxDirection andygaussian parameter in direction, calculatingxAndydirection dispersion offset wavenumber;
calculating the strength of the main earth magnetic field at each node according to the main earth magnetic field model, and taking the strength of the main earth magnetic field as the strength of a spatial domain background field;
obtaining a calculation model of the magnetization intensity according to the intensity of the spatial domain background field magnetic field and the intensity of the spatial domain abnormal field magnetic field;
converting a three-dimensional Laplace equation which meets the magnetic potential and the magnetization intensity of the abnormal field in the spatial domain into a one-dimensional ordinary differential equation in the spatial wave number domain by utilizing two-dimensional Fourier transform;
based on a space wave number domain one-dimensional ordinary differential equation and in combination with a set boundary condition which needs to be met by the magnetic potential of the abnormal field of the space wave number domain, converting an edge value problem model which is met by the magnetic potential of the abnormal field of the space wave number domain into an equivalent variational problem model;
obtaining the magnetic potential of the space wave number domain abnormal field by solving a variational problem model;
obtaining the magnetic field intensity of the spatial wave number domain abnormal field based on the magnetic potential of the spatial wave number domain abnormal field, and converting the magnetic potential of the spatial wave number domain abnormal field and the magnetic field intensity of the spatial wave number domain abnormal field into the magnetic potential of the spatial wave number domain abnormal field and the magnetic field intensity of the spatial wave number domain abnormal field through inverse Fourier transform;
the total magnetic field intensity of the space domain is the sum of the magnetic field intensity of the abnormal field of the space domain and the magnetic field intensity of the background field of the space domain, and the total magnetic field intensity of the space domain is tightened by using a tightening operator to obtain the total magnetic field intensity of the space domain after tightening;
judging whether iterative convergence conditions are met or not, and if yes, outputting a space domain abnormal field magnetic potential and a space domain abnormal field magnetic field strength corresponding to the tightened space domain total magnetic field strength obtained through current calculation;
and solving according to the output spatial domain abnormal field magnetic field intensity to obtain spatial domain magnetic induction intensity.
Compared with the prior art, the invention has the advantages that:
1. the invention can more accurately carry out magnetic field numerical simulation on the ferromagnetic medium by considering the self-demagnetization effect.
2. For the condition that the right end item of the model has unknown parameters, the iterative method is used for solving, and a solving idea is provided for similar problems.
3. Further, in order to reduce the truncation effect of the fourier transform, an edge expansion process is adopted.
4. And a tightening operator is added, so that the algorithm can be stably converged.
5. The invention reduces the three-dimensional problem into one dimension through Fourier transformation, improves the calculation precision and the calculation efficiency, and has good algorithm parallelism.
6. The method has great significance for estimating the reserves of the iron ore deposit and inverting and explaining the actually measured data.
Drawings
FIG. 1 is a flow chart in one embodiment of the present invention;
FIG. 2 is a schematic representation of an initial prism model in an embodiment of the present invention;
FIG. 3 is a schematic representation of the initial prismatic model of FIG. 2, subdivided;
FIG. 4 is a schematic diagram of a target region and an anomaly in an embodiment of the invention;
FIG. 5 is a diagram illustrating a numerical solution, an analytical solution, and an absolute error of an isotropic sphere model according to an embodiment of the present invention; wherein (a), (b) and (c) each representsB ax The numerical solution of (a) is,B ax an analytical solution ofB ax The absolute error of the numerical solution and the analytic solution of (a); (d) respectively, e and f representB ay The numerical solution of (a) is,B ay an analytical solution ofB ay The absolute error of the numerical solution and the analytic solution of (a); (g) (h), (i) respectively representB az The numerical solution of (a) is,B az an analytical solution ofB az The absolute error of the numerical solution and the analytic solution of (a);
fig. 6 is an internal structural diagram of a computer device in an embodiment of the present invention.
Detailed Description
For the purpose of promoting a clear understanding of the objects, aspects and advantages of the embodiments of the invention, reference will now be made to the drawings and detailed description, wherein there are shown in the drawings and described in detail, various modifications of the embodiments described herein, and other embodiments of the invention will be apparent to those skilled in the art. The exemplary embodiments of the present invention and the description thereof are provided to explain the present invention and not to limit the present invention.
Referring to fig. 1, the method for calculating magnetic induction intensity of a ferromagnetic medium based on a spatial wavenumber domain provided by the invention comprises the following steps:
(S1) acquiring a three-dimensional target region including an anomaly, and creating an initial prism model including the three-dimensional target region, as shown in fig. 2, wherein the initial prism model is a cuboid or a cube, and the anomaly may be in any shape. The anomaly in fig. 2 is a cuboid, and the anomaly is located inside the initial prism model.
(S2) for the initial prism model edgex、y、zThe directions are respectively divided at equal intervals inx、y、zA series of nodes are obtained in the direction as shown in fig. 3.
Then, the magnetic susceptibility is assigned to each node according to the magnetic susceptibility distribution data, the magnetic susceptibility is scalar, and the magnetic susceptibility is used for assigning the magnetic susceptibility to each nodeχExpressed in SI.
It is understood that a person skilled in the art can select a suitable assignment manner according to the existing method. In one embodiment of the invention, for the nodes of the anomaly portion, each node is assigned according to the susceptibility value of the anomaly; the magnetic susceptibility of the nodes of the non-anomalous body part is 0;
(S3) according to the prism model and the givenxDirection andygaussian parameter in direction, calculatingxAndydirection dispersion offset wavenumber;
(S4) calculating the main earth magnetic field strength at each node as the spatial domain background field strength according to the main earth magnetic field model;
(S5) obtaining a calculation model of the magnetization intensity according to the spatial domain background field magnetic field intensity and the spatial domain abnormal field magnetic field intensity;
(S6) converting a three-dimensional Laplace equation which meets the magnetic potential and the magnetization intensity of the abnormal field in the spatial domain into a one-dimensional ordinary differential equation in the spatial wavenumber domain by utilizing two-dimensional Fourier transform;
(S7) converting the boundary value problem model satisfied by the abnormal magnetic potential of the space wave number domain into an equivalent variational problem model based on a one-dimensional ordinary differential equation of the space wave number domain and by combining the set boundary condition to be satisfied by the abnormal magnetic potential of the space wave number domain;
(S8) obtaining the magnetic potential of the space wave number domain abnormal field by solving a variational problem model;
(S9) obtaining the magnetic field intensity of the space wave number domain abnormal field based on the magnetic potential of the space wave number domain abnormal field, and converting the magnetic potential of the space wave number domain abnormal field and the magnetic field intensity of the space wave number domain abnormal field into the magnetic potential of the space wave number domain abnormal field and the magnetic field intensity of the space wave number domain abnormal field through inverse Fourier transform;
(S10) the total magnetic field strength of the space domain is the sum of the magnetic field strength of the abnormal field of the space domain and the magnetic field strength of the background field of the space domain, and the total magnetic field strength of the space domain is tightened by a tightening operator to obtain the total magnetic field strength of the space domain after tightening;
(S11) judging whether the iterative convergence condition is met or not, if not, returning (S8), and if so, outputting the space domain abnormal field magnetic potential and the space domain abnormal field magnetic field strength corresponding to the tightened space domain total magnetic field strength obtained by current calculation;
and (S12) solving according to the output spatial domain abnormal field magnetic field intensity to obtain spatial domain magnetic induction intensity.
It is understood that, in the step (S3) of the present invention, the gaussian parameter may be set by referring to the method for setting the gaussian parameter among the existing methods in the art, and the setting may be performedxAndyand calculating direction discrete offset wave number.
In step (S3) of an embodiment of the present invention, axNumber of Gaussian points of directionN x Interval [ -1,1 [ ]]Upper gauss pointt a Gaussian coefficient ofA a Whereina=1,2,...,N x (ii) a Given ayNumber of Gaussian points of directionN y Interval [ -1,1 [ ]]Upper gauss pointt b Gaussian coefficient ofA b Whereinb=1,2,...,N y 。
in step (S3) of an embodiment of the present invention,xandythe direction dispersion offset wavenumber is calculated by the following method:
in the formula,
wherein:k x to representxThe number of shifted wave-numbers of the directions,
to representxThe number of direction base waves is the number of direction base waves,NN x representing the initial prism model edgexThe number of nodes obtained by direction equal-interval subdivision,
representing the initial prism model edgexUnit interval length adopted in the direction equal interval division;k y to representyThe number of shifted wave-numbers of the directions,
to representyThe number of direction base waves is the number of direction base waves,NN y representing the initial prism model edgeyThe number of nodes obtained by direction equal-interval subdivision,
representing the initial prism model edgeyThe unit interval length adopted when the direction is equally spaced.
In the step (S4) of an embodiment of the present invention, the method includes:
calculating the intensity of the main earth magnetic field at each node according to the model IGRF of the main earth magnetic field, and taking the intensity as the intensity H of the background field magnetic field in the space domain at each node0It is the background field in numerical simulation, i.e. the magnetic field when there is no anomaly, and the unit is A/m. Magnetic field intensity H of spatial domain background field0The components of the three directions are respectively represented as H 0x 、H y0、H z0：
Wherein:
represents H0The L2 norm of (a),αis the magnetic tilt angle of the target area,βis the target area declination.
In step (S5) of an embodiment of the present invention, the calculation model of the magnetization is:
wherein HaThe magnetic field intensity generated by the abnormal body at each node, namely the magnetic field intensity of the abnormal field in the space domain is represented, the abnormal field in numerical simulation, namely the magnetic field generated by the abnormal magnetic susceptibility is represented in the unit of A/m, and the three components are respectively H ax 、H ay 、H az . The total magnetic field intensity H of the space domain at each node is the ambient field magnetic field intensity H of the space domain0Magnetic field intensity H corresponding to spatial domain abnormal fieldaAnd (4) summing.
In step (S6) of one embodiment of the present invention, the magnetic potential U of the spatial domain abnormal field a Three-dimensional Laplace equation satisfied by magnetization M
In the formula,
is a Laplace operator, is unfolded into
，
Showing the intensity of magnetization
And (5) solving divergence. The above equation (7) is developed as follows:
the two-dimensional Fourier transform of equation (7) is performed, and in order to reduce the truncation effect of the Fourier transform, an edge-expanding process, that is, an edge-expanding process is adoptedx，yThe directions are extended for a certain distance on the basis of the target area to avoid the influence of the truncation effect on the field value calculation.
Obtaining one-dimensional ordinary differential equation of space wave number domain through two-dimensional Fourier transform, and reservingzThe direction is the spatial domain:
wherein
Representing the magnetic potential of the anomalous field in the spatial wavenumber domain,
、
、
in order to be the wave-number domain magnetization,k x 、k y are respectively asx、yOffset wavenumber of directions. Equation (8) which is the magnetic potential of the anomalous field in the space wavenumber domain
A one-dimensional ordinary differential equation satisfied.
In step (S7), in order to obtain a definite solution of the control equation, appropriate boundary conditions are given. Taking the Z axis vertically downwards as the forward direction and taking the horizontal ground as the upper boundary Z in the target area under a Cartesian coordinate systemminTaking a sufficient distance from the underground to the abnormal body as a lower boundary ZmaxAnd the upper and lower boundary conditions meet:
an upper boundary:
（9）
lower bound:
（10）
wherein,
。
and (3) simultaneous equations (8), (9) and (10) to obtain an edge value problem model satisfied by the abnormal magnetic field magnetic potential in the space wavenumber domain:
converting an edge value problem model satisfied by the magnetic potential of the space wave number domain abnormal field into an equivalent variational problem model by using a variational method:
in a Cartesian coordinate system, unit subdivision is carried out along the z direction, and a quadratic interpolation function, namely a space wave number domain abnormal field magnetic potential, is adopted in each unit
Changing twice within a cell.
Magnetic potential of space wave number domain abnormal field
The right term of the satisfied one-dimensional ordinary differential equation comprises a background field and an abnormal field, and the abnormal field is unknown, so that an iterative solution is adopted, and a solution thought is provided for similar problems (namely, under the condition that the right term of the model has unknown parameters, an iterative method is proposed for solving).
Combining equation (6) and equation (7), it can be seen that the magnetization M is determined by the background field H0And an anomalous field H a Sum and product of magnetic susceptibility, and H = H0+H a And an abnormal field H a Unknown, so the first iteration assumes H a And (3) the sum of the abnormal field and the background field in the first iteration is replaced by the background field, so that the one-dimensional partial differential equation is changed into the one-dimensional ordinary differential equation to be solved to obtain a first abnormal field, and then the sum of the abnormal field and the background field is used as a right-end term total field to be solved for the next time.
And the problem solved each time is a variational problem, the variational problem is subjected to unit analysis and total synthesis item by item to obtain a five-diagonal equation consisting of all nodes, the quick and efficient solution can be realized by adopting a catch-up method, and the magnetic potential of the spatial wave number domain abnormal field at each node is obtained
。
In step (S9) of one embodiment of the present invention, the magnetic potential is based on the anomalous field in the space wave number domain
And obtaining the magnetic field intensity of the abnormal field in the space wave number domain by the following formula:
wherein,iis the number of the imaginary numbers,
presentation pairzThe differential is obtained by the differential analysis,
representing magnetic potential pairs of space wavenumber domain abnormal fieldszAnd (6) carrying out differentiation.
The magnetic potential U of the abnormal field in the spatial wave number domain and the magnetic field strength of the abnormal field in the spatial wave number domain can be obtained through Fourier inverse transformation a And spatial domain anomalous field strength H a 。
In the step (S10) of an embodiment of the present invention, a tightening operator is applied to the total magnetic field strength of the space domain to obtain the total magnetic field strength of the space domain after tightening, where the total magnetic field strength of the space domain is the sum of the magnetic field strength of the abnormal field of the space domain and the magnetic field strength of the background field of the space domain.
Wherein,jrepresenting the number of iterations, H j+1Is as followsjThe total magnetic field strength of the spatial domain calculated in +1 iteration,
are respectively the firstjAmbient field strength in the spatial domain and anomalous field strength in the spatial domain, H, in +1 iteration j Is shown asjCalculated in the sub-iterationThe total magnetic field strength of the space domain.
It is understood that the preset iteration termination condition refers to a preset model calculation constraint condition for constraining the whole model to converge in the performance calculation process, so that the model can output a result meeting the condition. In the present invention, the iteration termination condition in step (S11) may be set as:
when the above iteration convergence condition is satisfied, the iteration is stopped, wherein
Is shown asjThe total magnetic field strength of the tightened space domain obtained by calculation in the secondary iteration,
is shown asjAnd (5) calculating the total magnetic field intensity of the tightened space domain in +1 iteration.
Of course, in practical applications, a person skilled in the art may set other iteration termination conditions based on the prior art, the conventional technical means in the field, or the common general knowledge, and is not limited to the iteration termination conditions set in the above preferred embodiments of the present application.
Finally, in step (S12), the magnetic induction B is induced by the abnormal field a (unit is T) and the magnetic field intensity H of the spatial domain abnormal field a From the relationship of (1), the magnetic induction B can be obtained a And further to obtain B a Three components ofB ax ，B ay ，B az 。
Wherein,μis the absolute permeability of the ferromagnetic medium, in units of H/m.
Absolute magnetic permeabilityμAndχthe relationship between them is shown in equation (17),μ 0in order to obtain a magnetic permeability in a vacuum,
。
the accuracy and the efficiency of the method for calculating the magnetic induction intensity of the strong magnetic medium based on the space wave number domain are tested.
The test computer is configured to i5-4590, the main frequency is 3.30GHz, and the memory is 12 GB.
The target area is a three-dimensional prismatic structure, the size of which is 1000m × 1000m × 1000m,x，ythe direction of each side expanding is 3000 m. The anomaly in the target region is an isotropic sphere with a radius of 200m, the sphere center being at the model center (500 m, 500m, 500 m), which is schematically shown in fig. 4.
The anomalous body is an isotropic sphere, and the sphere magnetic susceptibility of the anomalous body is 2 SI. The background magnetic field intensity of the model region is 50000nT, the magnetic dip angle of the target region is 45 degrees, and the magnetic declination angle of the target region is 9 degrees.x，y，z The number of direction nodes is 101,101,101. Under the condition of single thread, edge expanding FFT is adopted, iteration is needed for 7 times when the convergence condition is achieved, the time is 56s, the occupied memory is 1.38GB, the occupied memory is low, and the efficiency is high. B thereof ax ，B ay ，B az The absolute error of the analytic solution is different from the field value by more than two orders of magnitude, and the numerical calculation requirement is met. FIG. 5 is a model calculationB ax ，B ay ，B az The numerical solution, analytical solution and absolute error of the components are shown in FIG. 5, wherein (a), (b) and (c) representB ax The numerical solution of (a) is,B ax an analytical solution ofB ax The absolute error of the numerical solution and the analytic solution of (a); in FIG. 5, (d), (e) and (f) representB ay The numerical solution of (a) is,B ay an analytical solution ofB ay The absolute error of the numerical solution and the analytic solution of (a); in FIG. 5, (g), (h) and (i) representB az The numerical solution of (a) is,B az an analytical solution ofB az And resolving the absolute error of the solution.
On the basis of considering the self-demagnetization effect, a three-dimensional problem is reduced to one dimension by adopting standard edge expanding FFT (fast Fourier transform), only the z direction is reserved, a one-dimensional finite element method is applied, and shape function quadratic interpolation is adopted in an element, so that iterative solution is carried out on a differential equation, and a tightening operator is applied to ensure the convergence of the algorithm, thereby greatly improving the calculation precision and the calculation efficiency, and the algorithm has good parallelism and occupies less memory. The invention can be used in the inversion algorithm of magnetic measurement data, and improves the accuracy and efficiency of inversion.
In one embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be as shown in fig. 6. The computer device includes a processor, a memory, a network interface, and a database connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The database of the computer device is used to store sample data. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to realize the above-mentioned method for calculating magnetic induction of a ferromagnetic medium based on a spatial wavenumber domain.
Those skilled in the art will appreciate that the architecture shown in fig. 6 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided, which includes a memory and a processor, the memory stores a computer program, and the processor implements the steps of the magnetic induction intensity calculation method based on the strong magnetic medium in the spatial wave number domain in the above embodiments when executing the computer program.
In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored, which when executed by a processor implements the steps of the method for calculating magnetic induction of a strong magnetic medium based on the spatial wavenumber domain in the above-described embodiments.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory, among others. Non-volatile memory can include read-only memory (ROM), Programmable ROM (PROM), Electrically Programmable ROM (EPROM), Electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Link DRAM (SLDRAM), Rambus Direct RAM (RDRAM), direct bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (10)
1. The method for calculating the magnetic induction intensity of the strong magnetic medium based on the space wave number domain is characterized by comprising the following steps of:
acquiring a three-dimensional target area containing an abnormal body, and establishing an initial prism model containing the three-dimensional target area;
for the initial prism model edgex、y、zThe directions are respectively divided at equal intervals inx、y、zObtaining a series of nodes in the direction; assigning a magnetic susceptibility to each node based on the magnetic susceptibility distribution data, the magnetic susceptibility being a scalar quantity, usingχRepresenting, at nodes of the anomaly part, each node being assigned according to a susceptibility value of the anomaly; the magnetic susceptibility of the nodes of the abnormal part is 0;
according to prism models and givenxDirection andygaussian parameter in direction, calculatingxAndydirection dispersion offset wavenumber;
calculating the strength of the main earth magnetic field at each node according to the main earth magnetic field model, and taking the strength of the main earth magnetic field as the strength of a spatial domain background field;
obtaining a calculation model of the magnetization intensity according to the intensity of the spatial domain background field magnetic field and the intensity of the spatial domain abnormal field magnetic field;
converting a three-dimensional Laplace equation which meets the magnetic potential and the magnetization intensity of the abnormal field in the spatial domain into a one-dimensional ordinary differential equation in the spatial wave number domain by utilizing two-dimensional Fourier transform;
based on a space wave number domain one-dimensional ordinary differential equation and in combination with a set boundary condition which needs to be met by the magnetic potential of the abnormal field of the space wave number domain, converting an edge value problem model which is met by the magnetic potential of the abnormal field of the space wave number domain into an equivalent variational problem model;
obtaining the magnetic potential of the space wave number domain abnormal field by solving a variational problem model;
obtaining the magnetic field intensity of the spatial wave number domain abnormal field based on the magnetic potential of the spatial wave number domain abnormal field, and converting the magnetic potential of the spatial wave number domain abnormal field and the magnetic field intensity of the spatial wave number domain abnormal field into the magnetic potential of the spatial wave number domain abnormal field and the magnetic field intensity of the spatial wave number domain abnormal field through inverse Fourier transform;
the total magnetic field intensity of the space domain is the sum of the magnetic field intensity of the abnormal field of the space domain and the magnetic field intensity of the background field of the space domain, and the total magnetic field intensity of the space domain is tightened by using a tightening operator to obtain the total magnetic field intensity of the space domain after tightening;
judging whether iterative convergence conditions are met or not, and if yes, outputting a space domain abnormal field magnetic potential and a space domain abnormal field magnetic field strength corresponding to the tightened space domain total magnetic field strength obtained through current calculation;
and solving according to the output spatial domain abnormal field magnetic field intensity to obtain spatial domain magnetic induction intensity.
2. The method according to claim 1, wherein the gaussian parameter comprises:
xnumber of Gaussian points of directionN x Interval [ -1,1 [ ]]Upper gauss pointt a Gaussian coefficient ofA a Whereina=1,2,...,N x ；
ynumber of Gaussian points of directionN y Interval [ -1,1 [ ]]Upper gauss pointt b Gaussian coefficient ofA b Whereinb=1,2,...,N y 。
3. the method of calculating the magnetic induction of a strong magnetic medium based on the spatial wave number domain according to claim 2,xandythe direction dispersion offset wavenumber is calculated by the following method:
in the formula,
wherein:k x to representxThe number of shifted wave-numbers of the directions,
to representxThe number of direction base waves is the number of direction base waves,NN x representing the initial prism model edgexThe number of nodes obtained by direction equal-interval subdivision,
representing the initial prism model edgexUnit interval length adopted in the direction equal interval division;k y to representyThe number of shifted wave-numbers of the directions,
to representyThe number of direction base waves is the number of direction base waves,NN y to representInitial prism model edgeyThe number of nodes obtained by direction equal-interval subdivision,
representing the initial prism model edgeyThe unit interval length adopted when the direction is equally spaced.
4. The method for calculating the magnetic induction of a strong magnetic medium based on the spatial wavenumber domain according to claim 1, 2 or 3, wherein the calculation model of the magnetization is:
wherein H a Representing the magnetic field strength produced by the anomaly at each node, i.e. its spatial domain anomalous field strength, H0And H represents the spatial domain background field magnetic field strength at each node, and H represents the spatial domain total magnetic field strength at each node.
5. The method according to claim 4, wherein the method is applied to the magnetic induction of the spatial-wavenumber-domain-based ferromagnetic medium for calculating the magnetic potential of the anomalous field in the spatial domainU a And the magnetization M satisfies the three-dimensional Laplace equation
Performing two-dimensional Fourier transform to obtain a space wave number domain one-dimensional ordinary differential equation:
wherein
In order to be the laplacian operator,
showing the intensity of magnetization
Solving divergence;
representing the magnetic potential of the anomalous field in the spatial wavenumber domain,
、
、
are respectively asx、y、zThe wave-number domain magnetization of the direction,k x 、k y are respectively asx、yOffset wavenumber of directions.
6. The method for calculating the magnetic induction intensity of the ferromagnetic medium based on the spatial wavenumber domain according to claim 5, wherein the boundary conditions that the magnetic potential of the abnormal field in the spatial wavenumber domain needs to satisfy are as follows:
an upper boundary:
lower bound:
wherein,
。
7. the method for calculating the magnetic induction intensity of the ferromagnetic medium based on the spatial wavenumber domain according to claim 6, wherein the boundary conditions that the magnetic potential of the abnormal field in the spatial wavenumber domain needs to satisfy and the one-dimensional ordinary differential equation in the spatial wavenumber domain jointly form an edge value problem model that the magnetic potential of the abnormal field in the spatial wavenumber domain satisfies; and converting the boundary value problem model satisfied by the magnetic potential of the abnormal field in the space wavenumber domain into an equivalent variational problem model by using a variational method.
8. The method according to claim 7, wherein the variational problem model is solved by a shape function method to obtain the magnetic potential of the anomalous field in the space wavenumber domain
。
9. The method according to claim 8, wherein the tightening is performed by using a tightening operator to obtain the total magnetic field strength of the tightened spatial domain, and the method comprises the following steps:
wherein,jrepresenting the number of iterations, H j+1Is as followsjThe total magnetic field strength of the spatial domain calculated in +1 iteration,
are respectively the firstjAmbient field strength in the spatial domain and anomalous field strength in the spatial domain, H, in +1 iteration j Is shown asjAnd calculating the total magnetic field intensity of the space domain in the secondary iteration.
10. The method for calculating the magnetic induction of a strong magnetic medium based on the spatial wavenumber domain according to claim 1, 2, 3, 5, 6, 7, 8 or 9, wherein the iteration termination condition is set as:
wherein
Is shown asjThe total magnetic field strength of the tightened space domain obtained by calculation in the secondary iteration,
is shown asjAnd (5) calculating the total magnetic field intensity of the tightened space domain in +1 iteration.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111218085.7A CN113656750B (en) | 2021-10-20 | 2021-10-20 | Calculation method of magnetic induction intensity of ferromagnetic medium based on spatial wavenumber domain |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111218085.7A CN113656750B (en) | 2021-10-20 | 2021-10-20 | Calculation method of magnetic induction intensity of ferromagnetic medium based on spatial wavenumber domain |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113656750A true CN113656750A (en) | 2021-11-16 |
| CN113656750B CN113656750B (en) | 2021-12-17 |
Family
ID=78484267
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111218085.7A Active CN113656750B (en) | 2021-10-20 | 2021-10-20 | Calculation method of magnetic induction intensity of ferromagnetic medium based on spatial wavenumber domain |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113656750B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114036805A (en) * | 2021-11-26 | 2022-02-11 | 中南大学 | Forward modeling method, device, equipment and medium for 3D steady heat conduction geothermal field |
| CN114036806A (en) * | 2021-11-26 | 2022-02-11 | 中南大学 | Three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium |
| CN114065586A (en) * | 2021-11-22 | 2022-02-18 | 中南大学 | Three-dimensional magnetotelluric space-wavenumber domain finite element numerical simulation method |
| CN115292973A (en) * | 2022-10-09 | 2022-11-04 | 中南大学 | A method and system for numerical simulation of three-dimensional magnetic field in spatial wavenumber domain with arbitrary sampling |
| CN115795231A (en) * | 2022-10-09 | 2023-03-14 | 中南大学 | Space wave number mixed domain three-dimensional high-intensity magnetic field iteration method and system |
| CN116029110A (en) * | 2022-12-26 | 2023-04-28 | 长江岩土工程有限公司 | Magnetic field numerical simulation method and system considering demagnetization and remanence |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5724270A (en) * | 1996-08-26 | 1998-03-03 | He Holdings, Inc. | Wave-number-frequency adaptive beamforming |
| US20030060981A1 (en) * | 1999-04-02 | 2003-03-27 | Conoco Inc. | Nonlinear constrained inversion method to determine base of salt interface from gravity and gravity tensor data |
| CN107748834A (en) * | 2017-11-22 | 2018-03-02 | 中南大学 | A fast and high-precision numerical simulation method for calculating the magnetic field of undulating observation surface |
| CN108197389A (en) * | 2018-01-04 | 2018-06-22 | 中南大学 | Quick, the high resolution numerical simulation method in two-dimentional ferromagnetic magnetic field |
| CN112859162A (en) * | 2019-11-28 | 2021-05-28 | 中国石油天然气集团有限公司 | Three-dimensional sound wave equation space domain forward modeling method and device based on pseudo-spectrum method |
-
2021
- 2021-10-20 CN CN202111218085.7A patent/CN113656750B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5724270A (en) * | 1996-08-26 | 1998-03-03 | He Holdings, Inc. | Wave-number-frequency adaptive beamforming |
| US20030060981A1 (en) * | 1999-04-02 | 2003-03-27 | Conoco Inc. | Nonlinear constrained inversion method to determine base of salt interface from gravity and gravity tensor data |
| CN107748834A (en) * | 2017-11-22 | 2018-03-02 | 中南大学 | A fast and high-precision numerical simulation method for calculating the magnetic field of undulating observation surface |
| CN108197389A (en) * | 2018-01-04 | 2018-06-22 | 中南大学 | Quick, the high resolution numerical simulation method in two-dimentional ferromagnetic magnetic field |
| CN112859162A (en) * | 2019-11-28 | 2021-05-28 | 中国石油天然气集团有限公司 | Three-dimensional sound wave equation space domain forward modeling method and device based on pseudo-spectrum method |
Non-Patent Citations (4)
| Title |
|---|
| JIANPING TAN等: ""The Scheme Design and Application of Large Gap Magnetic", 《2009 INTERNATIONAL CONFERENCE ON MEASURING TECHNOLOGY AND MECHATRONICS AUTOMATION》 * |
| 戴世坤等: ""重力异常场空间-波数混合域三维数值模拟"", 《地球物理学报》 * |
| 李昆等: ""空间波数混合域磁异常场积分解三维数值模拟"", 《地球物理学报》 * |
| 赵东东等: ""基于积分解的空间-波数混合域二度体磁异常数值模拟"", 《物探化探计算技术》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114065586A (en) * | 2021-11-22 | 2022-02-18 | 中南大学 | Three-dimensional magnetotelluric space-wavenumber domain finite element numerical simulation method |
| CN114036805A (en) * | 2021-11-26 | 2022-02-11 | 中南大学 | Forward modeling method, device, equipment and medium for 3D steady heat conduction geothermal field |
| CN114036806A (en) * | 2021-11-26 | 2022-02-11 | 中南大学 | Three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium |
| CN114036806B (en) * | 2021-11-26 | 2024-07-09 | 中南大学 | Three-dimensional ground temperature field numerical simulation method based on heat conductivity anisotropic medium |
| CN115292973A (en) * | 2022-10-09 | 2022-11-04 | 中南大学 | A method and system for numerical simulation of three-dimensional magnetic field in spatial wavenumber domain with arbitrary sampling |
| CN115795231A (en) * | 2022-10-09 | 2023-03-14 | 中南大学 | Space wave number mixed domain three-dimensional high-intensity magnetic field iteration method and system |
| CN115795231B (en) * | 2022-10-09 | 2023-08-04 | 中南大学 | A three-dimensional strong magnetic field iterative method and system in the space wavenumber mixed domain |
| CN116029110A (en) * | 2022-12-26 | 2023-04-28 | 长江岩土工程有限公司 | Magnetic field numerical simulation method and system considering demagnetization and remanence |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113656750B (en) | 2021-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113656750A (en) | Calculation method of magnetic induction intensity of ferromagnetic medium based on spatial wavenumber domain | |
| CN112800657B (en) | Gravity field numerical simulation method and device based on complex terrain and computer equipment | |
| CN112287534A (en) | NUFFT-based two-dimensional magnetic anomaly fast forward modeling method and device | |
| CN114021408B (en) | Two-dimensional strong magnetic field numerical simulation method, device, equipment and medium | |
| CN113962077B (en) | Three-dimensional anisotropic strong magnetic field numerical simulation method, device, equipment and medium | |
| Pollok et al. | Magnetic field prediction using generative adversarial networks | |
| CN114970289B (en) | Three-dimensional magnetotelluric anisotropy forward modeling numerical simulation method, equipment and medium | |
| CN114065585A (en) | Three-dimensional electrical source numerical simulation method based on coulomb specification | |
| Dai et al. | The forward modeling of 3D gravity and magnetic potential fields in space-wavenumber domains based on an integral method | |
| CN114004127B (en) | Two-dimensional main shaft anisotropic strong magnetic field numerical simulation method, device, equipment and medium | |
| CN109254327B (en) | Exploration method and exploration system for three-dimensional ferromagnetic body | |
| CN115292973B (en) | A method and system for numerical simulation of three-dimensional magnetic field in space wavenumber domain with random sampling | |
| Wang et al. | Efficient 2D modeling of magnetic anomalies using NUFFT in the Fourier domain | |
| CN107748834A (en) | A fast and high-precision numerical simulation method for calculating the magnetic field of undulating observation surface | |
| CN116029110A (en) | Magnetic field numerical simulation method and system considering demagnetization and remanence | |
| CN111950108B (en) | Method for calculating gravity gradient tensor of second-degree variable density body | |
| Chai | AE equation of potential field transformations in the wavenumber domain and its application | |
| CN104808249A (en) | High-precision geomagnetic field modeling method | |
| CN116244877B (en) | Three-dimensional magnetic field numerical simulation method and system based on 3D Fourier transform | |
| JP2984108B2 (en) | DC superposition inductance calculator | |
| CN120276058A (en) | Method, device, equipment and medium for rapidly calculating three-dimensional ground magnetic tensor | |
| JPH05100000A (en) | DC superimposition iron loss calculator | |
| Wang et al. | Magnetic forward modelling for high susceptibility at arbitrary heights based on compressed 2D and 3D fast discrete# xD; convolution techniques | |
| CN113139289B (en) | A forward and inversion method of magnetic measurement data under demagnetization effect based on integral equation | |
| CN107703539B (en) | Anti-alias seismic data interpolation method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |