CN-121997663-A - Magnetic anomaly high-efficiency simulation method based on hybrid algorithm

CN121997663ACN 121997663 ACN121997663 ACN 121997663ACN-121997663-A

Abstract

The invention discloses a magnetic anomaly high-efficiency simulation method based on a hybrid algorithm, which comprises the steps of establishing a Cartesian coordinate system for a target area, acquiring information, observation point information and magnetic field intensity of any grid unit of the target area, dividing the grid unit by adopting a distance self-adaptive algorithm based on the information and the observation point information of the grid unit to acquire a plurality of near-field area grid units and a plurality of far-field area grid units, calculating nuclear matrix elements corresponding to any far-field area grid unit by adopting a point magnetic source algorithm, calculating the nuclear matrix elements corresponding to any near-field area grid unit by adopting a numerical integration method, synthesizing the nuclear matrix elements corresponding to the far-field area grid unit and the nuclear matrix elements corresponding to the near-field area grid unit, and obtaining a magnetic anomaly vector. Through the scheme, the method has the advantages of simple logic, accuracy, reliability and the like, and has high practical value and popularization value in the technical field of physical simulation.

Inventors

WANG CHANGNAN
Mei Ensheng
Min Kefan
Tang Gaolin
YANG GUANGYUAN
ZOU WEN
YU XIANGLONG
TAN WENJIE
LIANG JUN
PANG YIJIE
GUAN HUA

Assignees

四川里伍铜业股份有限公司

Dates

Publication Date: 20260508
Application Date: 20260128

Claims (7)

1. The magnetic anomaly high-efficiency simulation method based on the hybrid algorithm is characterized by comprising the following steps of: Establishing a Cartesian coordinate system for a target area, and acquiring information, observation point information and magnetic field intensity of any grid unit of the target area; Dividing the grid cells by adopting a distance self-adaptive algorithm based on the information of the grid cells and the information of the observation points to obtain a plurality of near-field region grid cells and a plurality of far-field region grid cells; Calculating the nuclear matrix element corresponding to any far-field area grid unit by adopting a point magnetic source algorithm, and calculating the nuclear matrix element corresponding to any near-field area grid unit by adopting a numerical integration method; and synthesizing the nuclear matrix elements corresponding to the far-field region grid cells and the nuclear matrix elements corresponding to the near-field region grid cells, and obtaining the magnetic anomaly vector.
2. The hybrid algorithm-based efficient simulation method of magnetic anomalies according to claim 1, wherein the information of the grid cells is Wherein, the method comprises the steps of, Representing the coordinates of the jth grid cell in the x-direction of the Cartesian coordinate system; Representing the coordinates of the jth grid cell in the y-direction of the Cartesian coordinate system; representing the coordinates of the jth grid cell in the z-direction of the Cartesian coordinate system; representing the width of the jth grid cell in the x-direction of the Cartesian coordinate system; representing the width of the jth grid cell in the y-direction of the Cartesian coordinate system; representing the width of the jth grid cell in the z-direction of the Cartesian coordinate system; indicating the magnetic susceptibility value of the jth grid cell.
3. The hybrid algorithm-based efficient simulation method of magnetic anomalies, as set forth in claim 2, wherein the observation point information is Wherein, the method comprises the steps of, Representing the coordinates of the ith observation point in the x direction of the Cartesian coordinate system; Representing the coordinates of the ith observation point in the y direction of the Cartesian coordinate system; Representing the coordinates of the i-th observation point in the z-direction of the cartesian coordinate system.
4. The method for efficiently simulating magnetic anomalies based on a hybrid algorithm according to claim 3, wherein the method for efficiently simulating magnetic anomalies based on information of grid cells and information of observation points and dividing the grid cells by a distance adaptive algorithm to obtain a plurality of near-field region grid cells and a plurality of far-field region grid cells comprises: if the grid cell satisfies: And (2) and If not, dividing the grid cell into a near field region grid cell; Wherein, the Representing the distance between the i-th observation point and the center point of the j-th grid cell.
5. The method of claim 4, further comprising dividing a j-th grid unit of the far field region grid units into k×m×n square units according to the length a, the width b and the height c, and calculating a nuclear matrix element corresponding to any far field region grid unit by using a point magnetic source algorithm, wherein the expression is: ; Wherein, the Representing a kernel function matrix element; T represents the magnetic field vector of the observation point to the center of the unit; representing the minimum side length of the far field area grid unit in a Cartesian coordinate system; Minimum side length of the far field area grid unit in a Cartesian coordinate system The expression of (2) is: wherein, the method comprises the steps of, A minimum side length of the grid cell representing the jth far field region in a Cartesian coordinate system; the expression of the magnetic field vector t of the observation point to the center of the unit is as follows: ; wherein, the method comprises the steps of, Representing a unit magnetic field component in an x-direction of a Cartesian coordinate system; representing a unit magnetic field component in the y-direction of the Cartesian coordinate system; representing a unit magnetic field component in a z-direction of a Cartesian coordinate system; representing the distance between the ith observation point and the jth cube unit in the x direction of the Cartesian coordinate system, Representing the distance between the ith observation point and the jth cube unit in the y direction of the Cartesian coordinate system, The distance between the ith observation point and the jth cube unit in the z direction of the Cartesian coordinate system is represented by k, wherein k represents the quotient of the length of the cube unit in the x direction of the Cartesian coordinate system divided by the minimum size, m represents the quotient of the length of the cube unit in the y direction of the Cartesian coordinate system divided by the minimum size, and n represents the quotient of the length of the cube unit in the z direction of the Cartesian coordinate system divided by the minimum size.
6. The efficient simulation method of magnetic anomalies based on the hybrid algorithm of claim 5, wherein the numerical integration method is adopted to calculate the nuclear matrix element corresponding to any near-field region grid unit, and the expression is as follows: wherein, the method comprises the steps of, Representing a unit magnetic field component in an x-direction of a Cartesian coordinate system; representing a unit magnetic field component in the y-direction of the Cartesian coordinate system; Represents the unit magnetic field component in the z-direction of the Cartesian coordinate system; Representing a sign of a constant-integral operation.
7. The efficient simulation method of magnetic anomaly based on a hybrid algorithm of claim 6, wherein the synthesis of the core matrix elements corresponding to the far field region grid cells and the core matrix elements corresponding to the near field region grid cells is performed, and the magnetic anomaly vector is obtained, and the expression is: wherein, the method comprises the steps of, Representing a magnetic anomaly vector; Represents the number of grid cells, j E ; Representing the number of observation points, i E 。

Description

Magnetic anomaly high-efficiency simulation method based on hybrid algorithm Technical Field The invention relates to the technical field of physical simulation, in particular to a magnetic anomaly high-efficiency simulation method based on a hybrid algorithm. Background In the field of physical simulation technology, calculating the spatial magnetic anomaly distribution caused by a magnetic medium is a fundamental and key technology. The technology provides explanation basis for a plurality of fields such as engineering and environmental geophysics (such as unexplosive object detection and archaeological investigation), material science and industrial detection (such as nondestructive inspection and magnetic characteristic imaging), planetary science and space physics (such as planetary shell magnetic structure inversion) and the like by constructing a linear forward model of 'magnetic source-magnetic field' response. To facilitate modeling of complex morphologies of magnetic media, the target region is typically split into a series of regular (e.g., cuboid) grid cells, and the susceptibility within each grid cell is assumed to be uniform. The core calculation process comprises the steps of precisely calculating the magnetic field contribution generated by each rectangular grid unit at each observation point through a numerical integration method based on a potential field theory to form a kernel function matrix, and finally obtaining magnetic anomalies through multiplication of the kernel function matrix and vectors representing susceptibility distribution. For example, the disclosure technology of the disclosure number CN119575490A, named as a self-constrained magnetic anomaly property inversion method and system based on a correlation coefficient, and the disclosure technology of the disclosure number CN120912802A, named as a geophysical heavy magnetic anomaly source model construction and rapid forward modeling method. Although the forward modeling method based on the global numerical integration can ensure the precision, the forward modeling method has obvious efficiency bottleneck when being applied to large-scale practical application, and the root cause of the forward modeling method is that the calculation complexity is high, namely, the calculation of each element of the kernel function matrix needs to perform triple numerical integration (needing eight values) on the grid unit body. The core disadvantage of the method is that the method cannot be optimized according to the physical rule, namely, the method adopts the same integral flow for all grid cells, the characteristic that far field can be approximated to a point magnetic source is not considered, and huge calculation resource waste is caused by the indiscriminate processing mode. Therefore, it is needed to propose a magnetic anomaly efficient simulation method based on a hybrid algorithm, which is simple in logic, accurate and reliable. Disclosure of Invention Aiming at the problems, the invention aims to provide a magnetic anomaly high-efficiency simulation method based on a hybrid algorithm, which adopts the following technical scheme: a magnetic anomaly high-efficiency simulation method based on a hybrid algorithm comprises the following steps: Establishing a Cartesian coordinate system for a target area, and acquiring information, observation point information and magnetic field intensity of any grid unit of the target area; Dividing the grid cells by adopting a distance self-adaptive algorithm based on the information of the grid cells and the information of the observation points to obtain a plurality of near-field region grid cells and a plurality of far-field region grid cells; Calculating the nuclear matrix element corresponding to any far-field area grid unit by adopting a point magnetic source algorithm, and calculating the nuclear matrix element corresponding to any near-field area grid unit by adopting a numerical integration method; and synthesizing the nuclear matrix elements corresponding to the far-field region grid cells and the nuclear matrix elements corresponding to the near-field region grid cells, and obtaining the magnetic anomaly vector. Compared with the prior art, the invention has the following beneficial effects: (1) The invention is based on the information of the grid cells and the information of the observation points, and adopts a distance self-adaptive algorithm to divide the grid cells to obtain a plurality of near-field region grid cells and a plurality of far-field region grid cells, which are used for distinguishing the point magnetic source calculation region and the integral calculation region. (2) According to the invention, the nuclear matrix element corresponding to any far-field area grid unit is calculated by adopting a point magnetic source algorithm, and the nuclear matrix element corresponding to any near-field area grid unit is calculated by adopting a numerical integration method, so that