Rock residual magnetization inclination angle calculation method and system based on magnetic susceptibility prior information

Abstract

This invention proposes a method and system for calculating the remanent magnetization dip angle of rocks based on prior information about magnetization. The method includes: dividing the target subsurface half-space into multiple two-dimensional semi-rectangular units by using magnetic profile data from pre-set ground measuring points as constraints; for each two-dimensional semi-rectangular unit, inverting the magnetic profile data to obtain the horizontal and vertical components of the total magnetization intensity; decomposing the horizontal and vertical components to obtain induced magnetization and remanent magnetization; and based on the vector relationship of the magnetization intensity components, using the horizontal and vertical components, induced magnetization, and remanent magnetization, and with magnetic susceptibility as prior information, obtaining the remanent magnetization dip angle of each two-dimensional semi-rectangular unit. This invention provides a method for calculating the remanent magnetization dip angle of rocks that effectively distinguishes between the induced magnetization direction and the remanent magnetization direction, achieving effective separation of their contributions.

Description

Technical Field

[0001] This invention relates to the field of geomagnetic exploration technology, and in particular to a method and system for calculating the residual magnetization dip angle of rocks based on magnetization prior information. Background Technology

[0002] The remanent magnetization vector of a rock records the complex tectonic evolution it has undergone throughout its long geological history, including paleomagnetic field information during initial formation, hydrothermal alteration, and even metamorphism. Reconstructing the remanent magnetization vector not only provides crucial evidence for identifying and dividing different geological units, but also offers important clues for revealing the evolutionary history and dynamic processes of rocks.

[0003] Traditional magnetization vector inversion methods can usually only obtain the total magnetization intensity vector that reflects the overall magnetism of the rock, and it is difficult to effectively separate the induced magnetization and residual magnetization components, thus limiting its application in geological interpretation and evolution analysis. Summary of the Invention

[0004] This invention provides a method and system for calculating the remanent magnetization dip angle of rocks based on magnetization prior information. Its main purpose is to provide a method for calculating the remanent magnetization dip angle of rocks that effectively distinguishes between the induced magnetization direction and the remanent magnetization direction, and achieves effective separation of their contributions.

[0005] In a first aspect, embodiments of the present invention provide a method for calculating the remanent magnetization dip angle of rocks based on magnetization prior information, including:

[0006] S1. Based on the magnetic profile data of the ground preset measuring points, the target underground half space is divided into two-dimensional half-grids, and the study area is divided into multiple two-dimensional half-rectangular units.

[0007] S2, for each two-dimensional semi-rectangular unit, the magnetic profile data is used as a constraint and inversion is performed to obtain the horizontal and vertical components of the total magnetization.

[0008] S3, decompose the horizontal component and the vertical component to obtain induced magnetization and residual magnetization;

[0009] S4. Based on the vector relationship of the magnetization intensity components, and according to the horizontal component and the vertical component, the induced magnetization and the remanent magnetization, the remanent magnetization tilt angle of each two-dimensional semi-rectangular unit is obtained using the magnetic susceptibility as prior information.

[0010] Further, step S2 is determined by the following formula:

[0011] ;

[0012] ;

[0013] ;

[0014] ;

[0015] ;

[0016] in, Indicates the first Magnetization intensity of the next iteration Indicates the first Magnetization intensity of the next iteration This represents the magnetic profile data. Indicates the first The forward modeling results of the magnetic profile data in the next iteration. Indicates the step size constant. This represents the data covariance matrix of the magnetic profile data. Indicates the first Jacobian matrix of the next iteration Indicates the first The weighting matrix of the next iteration. This represents the horizontal component. Represents the vertical component, This represents the total magnetization. Indicates the first Partial derivatives of the horizontal components in the next iteration Indicates the first Partial derivatives of the vertical component in the next iteration. Represents the identity matrix. Indicates the weighted parameters. Indicates the focus constraint parameters. This indicates a depth-weighted term.

[0017] Further, step S3 is determined by the following formula:

[0018] ;

[0019] in, This represents the horizontal component of the induced magnetization. This represents the vertical component of the induced magnetization. This represents the horizontal component of the residual magnetization. This represents the vertical component of the residual magnetization. Indicates the inclination angle of the Earth's magnetic field. This indicates the deflection angle of the geomagnetic field. The tilt angle representing the residual magnetization, The deflection angle represents the residual magnetization.

[0020] Further, step S4 is determined by the following formula:

[0021] ;

[0022] in, Indicates the residual magnetization tilt angle, This represents the magnetic susceptibility. Indicates the strength of the Earth's magnetic field. The deflection angle represents the total magnetization.

[0023] Furthermore, the magnetic profile data is obtained through field survey lines or by extracting corresponding profiles from gridded data.

[0024] Furthermore, the lateral length of the two-dimensional semi-rectangular unit is equal to the preset measurement point spacing of the magnetic measurement profile, and the vertical length of the two-dimensional semi-rectangular unit is determined by the reference measurement point spacing.

[0025] Furthermore, the lateral direction refers to the horizontal direction parallel to the magnetic measurement profile;

[0026] The vertical direction is the depth direction perpendicular to the magnetic profile.

[0027] Secondly, embodiments of the present invention provide a rock remanent magnetization dip calculation system based on magnetization prior information, comprising:

[0028] The subdivision module is used to perform two-dimensional semi-grid subdivision of the target underground half-space based on the magnetic profile data of the ground preset measuring points, and to divide the study area into multiple two-dimensional semi-rectangular units.

[0029] The inversion module is used to perform inversion on each two-dimensional semi-rectangular unit, using the magnetic profile data as a constraint, to obtain the horizontal and vertical components of the total magnetization.

[0030] The decomposition module is used to decompose the horizontal component and the vertical component to obtain induced magnetization and residual magnetization;

[0031] The calculation module is used to obtain the remanent magnetization tilt angle of each two-dimensional semi-rectangular unit based on the vector relationship of the magnetization intensity components, according to the horizontal component and the vertical component, the induced magnetization and the remanent magnetization, and using the magnetic susceptibility as prior information.

[0032] Thirdly, embodiments of the present invention provide a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the steps of the above-described method for calculating the remanent magnetization dip angle of rocks based on magnetization prior information.

[0033] Fourthly, embodiments of the present invention provide a computer storage medium storing a computer program, which, when executed by a processor, implements the steps of the above-described method for calculating the remanent magnetization dip angle of rocks based on magnetization prior information.

[0034] This invention proposes a method and system for calculating the remanent magnetization dip angle of rocks based on magnetization a priori information. This method calculates the remanent magnetization dip angle of rocks based on magnetization a priori information. Building upon two-and-a-half-dimensional magnetization vector inversion, it constructs a vector relationship between the horizontal and vertical components of the total magnetization intensity, combined with known magnetic susceptibility information, thereby effectively extracting the remanent magnetization dip angle. In geological contexts with significant remanent magnetization, this method can accurately distinguish between the induced magnetization direction and the remanent magnetization direction, achieving effective separation of their contributions.

[0035] This invention has good applicability and stability under strong remanent magnetization conditions, and can provide technical means for the interpretation of complex magnetic anomalies, thereby improving the application effect of magnetic data in resource exploration. It has important theoretical significance and application value. Attached Figure Description

[0036] Figure 1 A flowchart is provided for an embodiment of the present invention to illustrate a method for calculating the remanent magnetization dip angle of rocks based on magnetization prior information;

[0037] Figure 2 This is a schematic diagram of a rock remanent magnetization dip angle calculation system based on magnetization prior information, provided in an embodiment of the present invention.

[0038] The realization of the objective, functional features and advantages of the present invention will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0039] The embodiments of this application are described in detail below. Examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain this application, and should not be construed as limiting this application.

[0040] To enable those skilled in the art to better understand the solutions of this application, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this application, and not all of them. All other embodiments obtained by those skilled in the art based on the embodiments of this application without creative effort are within the scope of protection of this application.

[0041] In the embodiments of this application, "at least one" refers to one or more; "multiple" refers to two or more. In the description of this application, terms such as "first," "second," and "third" are used only for descriptive purposes and should not be construed as indicating or implying relative importance or order. Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of indicated technical features. Therefore, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of this application, "multiple" means two or more, unless otherwise explicitly specified.

[0042] References such as “one embodiment” or “some embodiments” as described in this specification mean that one or more embodiments of this application include a specific feature, structure, or characteristic described in connection with that embodiment. Therefore, the terms “comprising,” “including,” “having,” and variations thereof, as used in this specification, mean “including, but not limited to,” unless otherwise specifically emphasized.

[0043] In practical geophysical exploration, magnetic survey data typically contains superimposed information on both induced magnetization and remanent magnetization components. This provides the possibility of extracting the remanent magnetization vector based on specific prior information. Therefore, this invention provides a method for calculating the remanent magnetization dip angle of rocks based on prior magnetization information. Figure 1 A flowchart of a method for calculating the remanent magnetization dip of rocks based on magnetization prior information is provided in this embodiment of the invention, as shown below. Figure 1 As shown, the method includes:

[0044] S1. Based on the magnetic profile data of the ground preset measuring points, the target underground half space is divided into two-dimensional half-grids, and the study area is divided into multiple two-dimensional half-rectangular units.

[0045] First, geomagnetic observations are conducted along a certain survey line, and a complete set of measurement data is recorded in sequence according to the preset measurement point locations. This set of magnetic profile data includes supporting information such as the preset measurement point coordinates, observed magnetic field values, and elevation.

[0046] In one implementation, the magnetic profile data is obtained by field survey lines or by extracting corresponding profiles from gridded data.

[0047] In practice, the magnetic profile data can be obtained from the actual field survey line. Staff will lay out a straight survey line on the ground in the field, set up a preset measuring point at a certain distance, and use a magnetometer to measure the geomagnetic field data at each preset measuring point. The magnetic profile data of this survey line will be directly recorded.

[0048] The magnetic profile data can also be extracted from the gridded data. In the field, a large number of survey lines are usually laid out in an area to obtain a large number of discrete survey points. Then, through interpolation and gridding, magnetic anomaly data in the form of a two-dimensional regular grid is generated. This magnetic anomaly data is similar to a magnetic field distribution map. Each grid point has a corresponding magnetic field value. In this two-dimensional grid data, a straight path (virtual survey line) is manually selected, and the magnetic field, coordinates and other data corresponding to the grid points passed through this path are extracted and combined to form an equivalent magnetic profile data.

[0049] Therefore, in this embodiment of the invention, geomagnetic measurement data (i.e., magnetic profile data) corresponding to a survey line is required first. This data can be obtained directly from a field measurement along a line, or it can be used by cutting out data from a virtual survey line from a pre-made gridded magnetic field data of the entire area.

[0050] Next, the target underground half-space is divided into two-dimensional semi-grids. Specifically, the underground half-space can be divided into two-dimensional semi-grids based on the length and point spacing of the magnetic survey profile data, resulting in multiple two-dimensional semi-rectangular units. The lateral spacing of these two-dimensional semi-rectangular units is consistent with the point spacing of the profile, while the vertical spacing is usually set based on the reference point spacing, thereby dividing the study area into regular two-dimensional semi-rectangular units.

[0051] It should be noted that, in magnetic exploration, the term "subsurface semi-space" refers to the space above the ground that is filled with air, while the subsurface extends infinitely below. However, in actual practice, only the finite depth range to be studied is considered.

[0052] Furthermore, the lateral length of the two-dimensional semi-rectangular unit is equal to the preset measurement point spacing of the magnetic measurement profile, and the vertical length of the two-dimensional semi-rectangular unit is determined by the reference measurement point spacing.

[0053] Wherein, the lateral direction is the horizontal direction parallel to the magnetic measurement profile;

[0054] The vertical direction is the depth direction perpendicular to the magnetic profile.

[0055] It should also be noted that, in this embodiment of the invention, the preset measuring point spacing refers to the horizontal distance between two adjacent preset measuring points on the magnetic measurement profile.

[0056] In this embodiment of the invention, "horizontal" refers to the horizontal direction parallel to the magnetic measurement profile. The lateral spacing of the two-dimensional semi-rectangular units is equal to the point spacing of the preset measurement points on the magnetic measurement profile, so as to ensure that the sampling density of the two-dimensional semi-rectangular units matches that of the measured data.

[0057] In this embodiment of the invention, "vertical" refers to the depth direction that is perpendicularly downward. The vertical spacing of the two-dimensional semi-rectangular units is not necessarily equal to the point spacing, but is generally selected based on the reference point spacing. It can be equal to the reference point spacing, or slightly smaller or slightly larger than the reference point spacing, and can be adaptively adjusted according to the actual situation.

[0058] Therefore, in this embodiment of the invention, the underground research area is divided into regularly shaped two-dimensional semi-rectangular units. Under the assumption of two-dimensionality and semi-rectangularity, these rectangular units extend infinitely or with a fixed thickness along the direction of the vertical profile, forming two-dimensional semi-rectangular units for subsequent inversion calculations such as magnetic susceptibility and remanent magnetization inclination.

[0059] In this embodiment of the invention, by matching the point spacing laterally, the sampling density of the underground two-dimensional semi-rectangular unit can correspond to the observation density on the surface, avoiding calculation errors caused by scale mismatch. Furthermore, when the geological body (such as dikes, layered rock masses, etc.) in the study area has a long strike and minimal variation along the strike, using two-dimensional semi-rectangular units can significantly reduce the amount of computation, replacing complex three-dimensional modeling.

[0060] S2, for each two-dimensional semi-rectangular unit, the magnetic profile data is used as a constraint and inversion is performed to obtain the horizontal and vertical components of the total magnetization.

[0061] In this embodiment of the invention, for each two-dimensional semi-rectangular unit, magnetic profile data is used as a constraint and inversion is performed to obtain the horizontal and vertical components of the total magnetization of each two-dimensional semi-rectangular unit.

[0062] It should be noted that the magnetization vector refers to the magnetization intensity of a rock, which is a vector quantity comprising induced magnetization and remanent magnetization. Induced magnetization is caused by the Earth's magnetic field and is related to magnetic susceptibility, while remanent magnetization is the inherent magnetism retained during the rock's formation and has parameters such as magnitude, dip angle, and deflection angle. In two-and-a-half-dimensional profile problems, the focus is on the magnetization intensity components within the profile.

[0063] Inversion refers to using known surface magnetic profile data as observations to calculate the magnetization vector of each two-dimensional semi-rectangular subsurface unit, including the magnitude, direction, and remanent magnetization. Simply put, it involves using surface magnetic profile data to infer the strength and orientation of the magnetism of each small subsurface rock mass.

[0064] In this embodiment of the invention, the inversion of the two-dimensional semi-rectangular magnetization vector means that the physical properties and shape of each two-dimensional semi-rectangular unit remain unchanged along the direction of the vertical cross-section, and only change within the cross-sectional (lateral and vertical) plane. Therefore, this embodiment of the invention does not require a complete three-dimensional calculation, but only inverts the magnetization vector within the cross-sectional plane; and the magnetization vector includes directional information (i.e., remanent magnetization tilt angle), and is not just an inversion of scalar magnetic susceptibility.

[0065] In the actual inversion process, the measured magnetic profile data are used as constraints, and two-dimensional semi-rectangular units are used as the model space to be solved. The unknown is the magnetization vector of each two-dimensional semi-rectangular unit. The specific inversion process is as follows: first, set an initial value for the magnetization vector of each two-dimensional semi-rectangular unit; use a two-dimensional semi-forward modeling algorithm to calculate the theoretical magnetic anomaly of the surface corresponding to this initial model; compare the difference between the theoretical value and the measured magnetic profile data in step S1; continuously adjust the magnetization vector of each two-dimensional semi-rectangular unit through optimization algorithms to reduce the difference; after iterative convergence, an optimal distribution of magnetization vectors of each underground two-dimensional semi-rectangular unit is obtained.

[0066] As one implementation method, based on the magnetic profile data from step 1 and the two-dimensional semi-subdivided units established in step 2, a two-dimensional semi-magnetization vector inversion is performed. The inversion is calculated using the following iterative formula:

[0067] ;

[0068] in, Indicates the first Magnetization intensity of the next iteration Indicates the first Magnetization intensity of the next iteration This represents the magnetic profile data. This represents the forward modeling result of the magnetic profile data. Indicates the step size constant. This represents the data covariance matrix of the magnetic profile data. Represents the Jacobian matrix. This represents a weighted matrix.

[0069] The specific variables are as follows:

[0070] ;

[0071] ;

[0072] ;

[0073] ;

[0074] in, This represents the partial derivative of the horizontal component. This represents the partial derivative of the vertical component. Represents the identity matrix. Indicates the weighted parameters. Indicates the focus constraint parameters. This indicates a depth-weighted term.

[0075] The ultimate goal of this invention is to determine the remanent magnetization dip angle of a rock. The magnetization intensity vector includes two parts: the induced magnetization vector and the remanent magnetization vector. After obtaining the total magnetization intensity through this vector inversion, the induced and remanent parts are separated by combining the magnetization prior information, and the remanent magnetization dip angle can be further solved.

[0076] S3, decompose the horizontal component and the vertical component to obtain induced magnetization and residual magnetization;

[0077] In this embodiment of the invention, various magnetization intensities are decomposed into horizontal and vertical components. The total magnetization intensity obtained by inversion is combined with the pre-given magnetization a priori information to calculate the induced magnetization component. The two are subtracted to obtain the components of the residual magnetization. Then, using trigonometric function relationships, a mathematical formula for directly calculating the residual magnetization tilt angle is derived.

[0078] As one implementation method, step S3 is determined by the following formula:

[0079] ;

[0080] in, This represents the horizontal component of the induced magnetization. This represents the vertical component of the induced magnetization. This represents the horizontal component of the residual magnetization. This represents the vertical component of the residual magnetization. Indicates the inclination angle of the Earth's magnetic field. This indicates the deflection angle of the geomagnetic field. The tilt angle representing the residual magnetization, The deflection angle represents the residual magnetization.

[0081] In this embodiment of the invention, the magnetic susceptibility is fixed as prior information, and the induced magnetization is directly calculated instead of being used as an unknown quantity in the inversion, thereby reducing the ambiguity of the inversion.

[0082] S4. Based on the vector relationship of the magnetization intensity components, and according to the horizontal component and the vertical component, the induced magnetization and the remanent magnetization, the remanent magnetization tilt angle of each two-dimensional semi-rectangular unit is obtained using the magnetic susceptibility as prior information.

[0083] In one implementation method, step S4 is determined by the following formula:

[0084] ;

[0085] in, Indicates the residual magnetization tilt angle, This represents the magnetic susceptibility. Indicates the strength of the Earth's magnetic field. The deflection angle represents the total magnetization intensity.

[0086] in, The magnetic susceptibility is a fundamental magnetic parameter of rocks, representing their ability to be magnetized under the influence of an external magnetic field. It can be measured from rock outcrops or borehole cores in the field and is a priori information that needs to be input into the calculation.

[0087] Indicates the strength of the Earth's magnetic field. The total magnetization deflection angle is represented by the total magnetization deflection angle. In the two-dimensional semi-magnetization vector inversion, it is usually assumed that the direction of the survey line is consistent with the direction of the total magnetization deflection angle.

[0088] In addition, the estimated value of the total magnetization deflection can be obtained in advance by introducing other magnetization direction estimation methods.

[0089] In this embodiment of the invention, based on the formula for calculating the residual magnetization inclination derived in the previous step, the residual magnetization inclination of each two-dimensional semi-rectangular unit divided underground is calculated by substituting its corresponding magnetic data, and finally the distribution result of the residual magnetization inclination of all two-dimensional semi-rectangular units in the entire study area is obtained.

[0090] This invention, in its embodiments, combines a complete two-dimensional semi-magnetic inversion process, outlining its advantages in six dimensions: computational accuracy, solution stability, computational efficiency, data applicability, engineering practicality, and target relevance.

[0091] (1) Introducing a priori constraints on magnetization greatly alleviates the ambiguity of magnetic inversion.

[0092] Conventional vector magnetic inversion simultaneously solves for magnetic susceptibility, remanent magnetization, and tilt angle, resulting in numerous unknowns, multiple solutions, and distorted results. In contrast, this invention uses magnetic susceptibility as prior known information to directly calculate the induced magnetization component, treating only remanent magnetization-related parameters as unknowns, thus reducing the number of unknowns in the model. This effectively narrows the solution space, reduces inversion iteration divergence and local optima problems, and improves the reliability of remanent magnetization tilt angle calculation.

[0093] (2) Two-dimensional semi-mesh modeling is adopted to balance computational accuracy and computational efficiency.

[0094] The embodiments of this invention are adapted to the scenario of single magnetic profile data in the field. The two-and-a-half assumption is applicable to most exploration targets such as rock masses, veins, and strata with stable strike extension, and has a high degree of geological fit. Compared with full three-dimensional grid inversion, the number of grids is greatly reduced, the calculation time for forward modeling and iterative inversion is significantly reduced, and the requirements for computer hardware are lower.

[0095] Furthermore, in this embodiment of the invention, the horizontal spacing of the underground grid matches the distance between the surface measuring points, and the spatial sampling scale is unified, avoiding calculation deviations caused by scale mismatch, and the model discretization is more reasonable.

[0096] (3) The data source is flexible and highly adaptable to actual production in the field.

[0097] There are two reliable sources for input magnetic profiles:

[0098] It can directly use the original measured profile data collected in the field along the survey line, and is suitable for conventional field magnetic exploration operations;

[0099] Virtual profiles can be extracted from regional gridded magnetic survey data, making full use of existing large-area exploration databases, eliminating the need for repeated fieldwork and saving exploration costs.

[0100] (4) The component vector decomposition method can directly output the residual magnetization tilt angle quantitatively.

[0101] This invention establishes an explicit dip angle calculation formula based on the vector decomposition relationship of total magnetization intensity, eliminating the need for complex quadratic fitting and empirical diagrams. It achieves full-grid, unit-by-unit quantitative calculation, outputting a continuous and refined residual magnetization dip angle distribution across the entire subsurface region of the profile, rather than a single-point qualitative judgment. It distinguishes between induced magnetization and residual magnetization contributions, accurately extracting residual magnetic signals related to paleomagnetism in rocks, serving the interpretation of tectonic, paleomagnetic, and mineral magnetic properties.

[0102] (5) The process is standardized and modularized, which facilitates its procedural implementation.

[0103] The entire process is logically closed-loop and clearly structured: data input → mesh generation → vector inversion → magnetization decomposition → dip angle calculation by mesh, with each step being independent and controllable; each step has clear input and output, making it easy to write computer programs for automated calculation, and can be integrated into gravity and magnetic processing and interpretation software to process multiple magnetic profiles in batches, making it suitable for large-scale production applications.

[0104] (6) It has stronger quantitative interpretation capabilities and improves the geological interpretation level of magnetic anomalies.

[0105] Traditional magnetic methods often only invert scalar values ​​of magnetic susceptibility, failing to distinguish between induced and remanent magnetization. This invention, however, fully preserves the vector information of magnetization intensity and quantitatively calculates the remanent magnetization dip, enabling it to: determine the direction of the paleomagnetic field during the formation of rock masses, differentiate between mineralized / stratigraphic magnetic anomalies, analyze the magnetic alteration of rock masses by tectonic movements, and expand the geological interpretation depth of magnetic data.

[0106] Figure 2 This is a schematic diagram of a rock remanent magnetization dip calculation system based on magnetization prior information, provided in an embodiment of the present invention. Figure 2 As shown, the system includes:

[0107] The subdivision module 210 is used to perform two-dimensional semi-grid subdivision of the target underground half-space based on the magnetic profile data of the ground preset measuring points, and to divide the study area into multiple two-dimensional semi-rectangular units.

[0108] Inversion module 220 is used to invert each two-dimensional semi-rectangular unit using the magnetic profile data as a constraint to obtain the horizontal and vertical components of the total magnetization.

[0109] The decomposition module 230 is used to decompose the horizontal component and the vertical component to obtain induced magnetization and residual magnetization;

[0110] The calculation module 240 is used to obtain the remanent magnetization tilt angle of each two-dimensional semi-rectangular unit based on the vector relationship of the magnetization intensity components, according to the horizontal component and the vertical component, the induced magnetization and the remanent magnetization, and using the magnetic susceptibility as prior information.

[0111] This embodiment is a system embodiment corresponding to the above method embodiment. Its specific implementation process is the same as that of the above method embodiment. For details, please refer to the above method embodiment. This system embodiment does not make any specific limitations.

[0112] The modules in the aforementioned rock remanent magnetization dip calculation system based on magnetization prior information can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device, or stored in the computer device's memory as software, so that the processor can call and execute the corresponding operations of each module.

[0113] In one embodiment, a computer device is provided, which may be a server. The computer device includes a processor, memory, a network interface, and a database connected via a system bus. The processor provides computational and control capabilities. The memory includes a computer storage medium and internal memory. The computer storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the computer storage medium. The database stores data generated or acquired during the execution of a method for calculating the remanent magnetization dip of rocks based on magnetization a priori information, such as magnetic profile data, and the horizontal and vertical components of total magnetization. The network interface communicates with external terminals via a network connection. When the computer program is executed by the processor, it implements a method for calculating the remanent magnetization dip of rocks based on magnetization a priori information.

[0114] In one embodiment, a computer device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the steps of a method for calculating the remanent magnetization dip angle of rocks based on magnetization a priori information as described in the above embodiment. Alternatively, when the processor executes the computer program, it implements the functions of each module / unit in this embodiment of a system for calculating the remanent magnetization dip angle of rocks based on magnetization a priori information.

[0115] In one embodiment, a computer storage medium is provided, on which a computer program is stored. When executed by a processor, the computer program implements the steps of the rock remanent magnetization dip angle calculation method based on magnetization a priori information described in the above embodiment. Alternatively, when executed by a processor, the computer program implements the functions of each module / unit in the rock remanent magnetization dip angle calculation system based on magnetization a priori information described in this embodiment.

[0116] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments described above. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may 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), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), RAMbus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0117] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the above-described division of functional units and modules is used as an example. In practical applications, the above functions can be assigned to different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above.

[0118] The above-described embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be included within the protection scope of the present invention.

Claims

1. A method for calculating the remanent magnetization dip angle of rocks based on prior magnetization information, characterized in that, include: S1. Based on the magnetic profile data of the ground preset measuring points, the target underground half space is divided into two-dimensional half-grids, and the study area is divided into multiple two-dimensional half-rectangular units. S2, for each two-dimensional semi-rectangular unit, the magnetic profile data is used as a constraint and inversion is performed to obtain the horizontal and vertical components of the total magnetization. S3, decompose the horizontal component and the vertical component to obtain induced magnetization and residual magnetization; S4. Based on the vector relationship of the magnetization intensity components, and according to the horizontal component and the vertical component, the induced magnetization and the remanent magnetization, the remanent magnetization tilt angle of each two-dimensional semi-rectangular unit is obtained using the magnetic susceptibility as prior information.

2. The method for calculating the remanent magnetization dip angle of rocks based on prior magnetization information according to claim 1, characterized in that, Step S2 is determined by the following formula: ； ； ； ； ； in, Indicates the first Magnetization intensity of the next iteration Indicates the first Magnetization intensity of the next iteration This represents the magnetic profile data. Indicates the first The forward modeling results of the magnetic profile data in the next iteration. Indicates the step size constant. This represents the data covariance matrix of the magnetic profile data. Indicates the first Jacobian matrix of the next iteration Indicates the first The weighting matrix of the next iteration. This represents the horizontal component. Represents the vertical component, This represents the total magnetization. Indicates the first Partial derivatives of the horizontal components in the next iteration Indicates the first Partial derivatives of the vertical component in the next iteration. Represents the identity matrix. Indicates the weighted parameters. Indicates the focus constraint parameters. This indicates a depth-weighted term.

3. The method for calculating the remanent magnetization dip angle of rocks based on prior magnetization information according to claim 2, characterized in that, Step S3 is determined by the following formula: ； in, This represents the horizontal component of the induced magnetization. This represents the vertical component of the induced magnetization. This represents the horizontal component of the residual magnetization. This represents the vertical component of the residual magnetization. Indicates the inclination angle of the Earth's magnetic field. This indicates the deflection angle of the geomagnetic field. The tilt angle representing the residual magnetization, The deflection angle represents the residual magnetization.

4. The method for calculating the remanent magnetization dip angle of rocks based on prior magnetization information according to claim 3, characterized in that, Step S4 is determined by the following formula: ； in, Indicates the residual magnetization tilt angle, This represents the magnetic susceptibility. Indicates the strength of the Earth's magnetic field. The deflection angle represents the total magnetization.

5. The method for calculating the remanent magnetization dip angle of rocks based on prior magnetization information according to claim 1, characterized in that, The magnetic profile data is obtained through field survey lines or by extracting corresponding profiles from gridded data.

6. The method for calculating the remanent magnetization dip angle of rocks based on magnetization prior information according to claim 1, characterized in that, The lateral length of the two-dimensional semi-rectangular unit is equal to the preset measurement point spacing of the magnetic measurement profile, and the vertical length of the two-dimensional semi-rectangular unit is determined by the reference measurement point spacing.

7. The method for calculating the remanent magnetization dip angle of rocks based on prior magnetization information according to claim 6, characterized in that, The lateral direction is the horizontal direction parallel to the magnetic measurement profile; the vertical direction is the depth direction perpendicular to the magnetic measurement profile.

8. A system for calculating the remanent magnetization dip angle of rocks based on prior magnetization information, characterized in that, include: The subdivision module is used to perform two-dimensional semi-grid subdivision of the target underground half-space based on the magnetic profile data of the ground preset measuring points, and to divide the study area into multiple two-dimensional semi-rectangular units. The inversion module is used to perform inversion on each two-dimensional semi-rectangular unit, using the magnetic profile data as a constraint, to obtain the horizontal and vertical components of the total magnetization. The decomposition module is used to decompose the horizontal component and the vertical component to obtain induced magnetization and residual magnetization; The calculation module is used to obtain the remanent magnetization tilt angle of each two-dimensional semi-rectangular unit based on the vector relationship of the magnetization intensity components, according to the horizontal component and the vertical component, the induced magnetization and the remanent magnetization, and using the magnetic susceptibility as prior information.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the rock remanent magnetization dip angle calculation method based on magnetization prior information as described in any one of claims 1 to 7.

10. A computer storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the steps of the rock remanent magnetization dip angle calculation method based on magnetization prior information as described in any one of claims 1 to 7.

Images