CN117647847A - Crust magnetic structure inversion method based on depth constraint of inside and magnetic substrate - Google Patents
Crust magnetic structure inversion method based on depth constraint of inside and magnetic substrate Download PDFInfo
- Publication number
- CN117647847A CN117647847A CN202311675344.8A CN202311675344A CN117647847A CN 117647847 A CN117647847 A CN 117647847A CN 202311675344 A CN202311675344 A CN 202311675344A CN 117647847 A CN117647847 A CN 117647847A
- Authority
- CN
- China
- Prior art keywords
- depth
- magnetic
- crust
- inversion
- basin
- 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.)
- Pending
Links
- 230000005291 magnetic effect Effects 0.000 title claims abstract description 177
- 238000000034 method Methods 0.000 title claims abstract description 74
- 239000000758 substrate Substances 0.000 title claims abstract description 63
- 239000011435 rock Substances 0.000 claims abstract description 44
- 238000011160 research Methods 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- LOUPRKONTZGTKE-WZBLMQSHSA-N Quinine Chemical compound C([C@H]([C@H](C1)C=C)C2)C[N@@]1[C@@H]2[C@H](O)C1=CC=NC2=CC=C(OC)C=C21 LOUPRKONTZGTKE-WZBLMQSHSA-N 0.000 claims abstract description 30
- 235000001258 Cinchona calisaya Nutrition 0.000 claims abstract description 15
- LOUPRKONTZGTKE-UHFFFAOYSA-N cinchonine Natural products C1C(C(C2)C=C)CCN2C1C(O)C1=CC=NC2=CC=C(OC)C=C21 LOUPRKONTZGTKE-UHFFFAOYSA-N 0.000 claims abstract description 15
- 229960000948 quinine Drugs 0.000 claims abstract description 15
- 238000011835 investigation Methods 0.000 claims abstract description 9
- 238000004590 computer program Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 6
- 230000010287 polarization Effects 0.000 claims description 5
- 238000009933 burial Methods 0.000 claims description 3
- 206010054949 Metaplasia Diseases 0.000 claims description 2
- 230000015689 metaplastic ossification Effects 0.000 claims description 2
- 101150077190 sinI gene Proteins 0.000 claims 1
- 230000008859 change Effects 0.000 abstract description 5
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 4
- 239000000463 material Substances 0.000 abstract description 4
- 239000011707 mineral Substances 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 7
- 238000005553 drilling Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000002844 melting Methods 0.000 description 2
- 230000008018 melting Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- SZVJSHCCFOBDDC-UHFFFAOYSA-N iron(II,III) oxide Inorganic materials O=[Fe]O[Fe]O[Fe]=O SZVJSHCCFOBDDC-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- 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
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a crust magnetic structure inversion method based on the depth constraint of an inside and a magnetic substrate, which comprises the following steps: performing tangent method and external quinine method calculation on the actually measured aeromagnetic section of the research area to obtain the depth of the magnetic substrate; obtaining the Curie depth of the research area through a one-dimensional steady-state heat conduction equation based on the ground temperature gradient of the basin in the research area, the heat yield of the crust rock and the heat conductivity of the crust rock; and under the constraint of the depth of the magnetic substrate and the depth of the Curie face, performing depth weighted magnetic anomaly three-dimensional inversion to obtain the magnetic susceptibility of the basin crust. The method effectively relieves the problem of longitudinal resolution of susceptibility inversion, accurately constrains transverse change of the magnetic layer, provides more reliable results for detecting the crust material structure of the basin, and further provides mineral resource investigation for service science research.
Description
Technical Field
The invention relates to the technical field of geophysics, in particular to a crust magnetic structure inversion method based on the depth constraint of an inside and a magnetic substrate.
Background
The longitudinal resolution of the aeromagnetic anomaly is limited, and the corresponding magnetic susceptibility three-dimensional inversion has strong multi-resolution, so that the understanding of the crust magnetic structure is restricted. The magnetic structure of the crust is of great importance for understanding the composition, origin and mineral resource effects of the crust. The conventional crust susceptibility three-dimensional inversion has no constraint of upper and lower limits of the magnetic layer, resulting in a stronger magnetic distribution in the calculated result in nonmagnetic or weakly magnetic sedimentary layers and rock rings below curie, which is clearly different from the fact.
Disclosure of Invention
The invention provides a crust magnetic structure inversion method based on the depth constraint of an inside and a magnetic substrate, which mainly aims to effectively relieve the problem of longitudinal resolution of magnetic susceptibility inversion, accurately constrain the transverse change of a magnetic layer and provide a more reliable result for the crust material structure detection of a basin.
In a first aspect, an embodiment of the present invention provides a method for inversion of a crust magnetic structure based on depth constraints of an in-house and magnetic substrate, including:
performing tangent method and external quinine method calculation on the actually measured aeromagnetic section of the research area to obtain the depth of the magnetic substrate;
obtaining the Curie depth of the research area through a one-dimensional steady-state heat conduction equation based on the ground temperature gradient of the basin in the research area, the heat yield of the crust rock and the heat conductivity of the crust rock;
and under the constraint of the depth of the magnetic substrate and the depth of the Curie face, performing depth weighted magnetic anomaly three-dimensional inversion to obtain the magnetic susceptibility of the basin crust.
Further, the step of obtaining the depth of the magnetic substrate includes the steps of:
performing tangent method and external quinine method calculation on the aeromagnetic section to obtain a minimum buried depth value of the magnetic body;
removing shallow magnetic rock mass which does not accord with preset rock type information in the minimum burial depth value of the magnetic body;
and carrying out interpolation treatment on the minimum embedded depth value of the removed magnetic body to obtain the magnetic substrate depth.
Further, the step of obtaining the curie depth of the investigation region by a one-dimensional steady state heat conduction equation based on the geothermal gradient of the basin in the investigation region, the thermal yield of the crust rock and the thermal conductivity of the crust rock comprises:
constructing the one-dimensional steady-state heat conduction equation based on the ground temperature gradient of the basin in the research area, the heat yield of the crust rock and the heat conductivity of the crust rock, and obtaining the temperature structures of the crust and the top of the rock circle;
performing vertical linear interpolation on the temperature structures at the tops of the crust and the rock ring to obtain Curie point depths of preset discrete points;
and carrying out minimum curvature interpolation on the Curie point depth of each preset discrete point on a plane to obtain the Curie depth of the research area.
Further, the one-dimensional steady-state heat conduction equation is calculated as follows:
where z represents a depth variable, T (z) represents a temperature at the depth z, H 0 Represents the heat yield rate at the shallow part of the crust, h r Representing the decay index of the thermal yield with depth, k (T (z)) representing the thermal conductivity,representing the ground temperature gradient.
Further, the step of obtaining the magnetic susceptibility of the basin crust by performing depth weighted three-dimensional inversion of the magnetic anomaly under the constraint of the magnetic base depth and the curie face depth comprises:
an objective function is established for depth weighted inversion of susceptibility structure detection in the basin,
and taking the depth of the magnetic substrate as an upper interface constraint, the depth of the Curie surface as a lower interface constraint, taking the statistical data of the magnetic susceptibility measured in the field as a physical constraint, and carrying out aeromagnetic anomaly depth weighted inversion according to the objective function to obtain the magnetic susceptibility of the basin crust.
Further, the calculation formula of the objective function is as follows:
wherein μ represents a regularization factor, d represents observed data, m represents model parameters, G represents a forward operator of uniformly magnetized cubic magnetic anomalies, W d Represents a diagonal weighting matrix, W m Representing a model weighting matrix consisting of the sum of the zero-order and first-order finite differences of the model, Z represents a depth weighting matrix that decays exponentially with depth.
Further, before the step of establishing the objective function for depth weighted inversion of susceptibility structure detection in the basin, the method further comprises:
and performing pole-changing treatment on the grid magnetic anomaly actually measured in the research area to obtain a variable dip angle pole-changing result of the basin.
Further, the calculation formula of the polarization treatment is as follows:
S(u,v)=S T (u,v)·H(u,v),
α 0 =cos I 0 ·cos D 0 ,
β 0 =cos I 0 ·sin D 0 ,
γ 0 =sin I 0 ,
wherein S (u, v) represents post-metaplasiaFourier transform, S T (u, v) represents the Fourier transform of the grid magnetic anomaly, H (u, v) represents the Fourier transform of the polarization factor, I represents the imaginary unit, (u, v) represents the frequency domain argument, I 0 Represents geomagnetic inclination angle, D 0 Representing the geomagnetic declination.
In a second aspect, embodiments of the present invention provide an in-house and magnetic-basement depth-constrained basin crust magnetic computing system comprising:
the magnetic substrate depth calculation module is used for carrying out tangential method and external quinine method calculation on the actually measured aeromagnetic section of the research area to obtain the magnetic substrate depth;
the Curie depth calculation module is used for obtaining the Curie depth of the research area through a one-dimensional steady-state heat conduction equation based on the ground temperature gradient of the basin in the research area, the heat yield of the crust rock and the heat conductivity of the crust rock;
and the magnetic susceptibility calculation module is used for carrying out depth weighted magnetic anomaly three-dimensional inversion under the constraint of the magnetic substrate depth and the Curie depth to obtain the magnetic susceptibility of the basin crust.
In a third aspect, embodiments of the present invention provide a computer storage medium storing a computer program which, when executed by a processor, implements the steps of a method of inversion of a crustal magnetic structure based on in-house and magnetic substrate depth constraints described above.
According to the inversion method of the crust magnetic structure based on the depth constraint of the inside and the magnetic substrate, the depth of the magnetic substrate is obtained through calculation by an efficient manual tangent method and an external quinine method; then, based on heat flow data and heat conductivity and heat yield data existing in the basin, obtaining the depth of 580 ℃ inside by solving a one-dimensional steady-state heat conduction equation; and under the constraint of the depth of the magnetic substrate and the depth of the Curie face, performing depth weighted magnetic anomaly three-dimensional inversion to obtain the magnetic susceptibility three-dimensional structure of the basin crust. The method effectively relieves the problem of longitudinal resolution of susceptibility inversion, accurately constrains transverse change of the magnetic layer, provides more reliable results for detecting the crust material structure of the basin, and further provides mineral resource investigation for service science research.
Drawings
FIG. 1 is a flow chart of a method for inverting a crustal magnetic structure based on depth constraints of an inside and a magnetic substrate provided by an embodiment of the invention;
FIG. 2 is a schematic diagram of a tangent method according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a crustal magnetic structure inversion system based on depth constraints of an in-house and magnetic substrate according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of a computer device according to an embodiment of the present invention.
The achievement of the objects, functional features and advantages of the present invention will be further described with reference to the accompanying drawings, in conjunction with the embodiments.
Detailed Description
Embodiments of the present application are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements or elements having like or similar functions throughout. The embodiments described below by referring to the drawings are exemplary only for the purpose of explaining the present application and are not to be construed as limiting the present application.
In order to better understand the solution of the present application, the following description will make clear and complete descriptions of the technical solution of the embodiment of the present application with reference to the accompanying drawings in the embodiment of the present application. It will be apparent that the described embodiments are only some, but not all, of the embodiments of the present application. All other embodiments, which can be made by those skilled in the art based on the embodiments herein without making any inventive effort, are intended to be within the scope of the present application.
In the embodiment of the application, at least one refers to one or more; plural means two or more. In the description of the present application, the words "first," "second," "third," and the like are used solely for the purpose of distinguishing between descriptions and not necessarily for the purpose of indicating or implying a relative importance or order. Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the present application, the meaning of "a plurality" is two or more, unless explicitly defined otherwise.
Reference in the specification to "one embodiment" or "some embodiments" or the like means that a particular feature, structure, or characteristic described in connection with the embodiment is included in one or more embodiments of the application. Thus, the terms "comprising," "including," "having," and variations thereof herein mean "including but not limited to," unless expressly specified otherwise.
Fig. 1 is a flowchart of a method for inverting a crust magnetic structure based on depth constraint of an inside and a magnetic substrate according to an embodiment of the present invention, as shown in fig. 1, including:
s110, performing tangent method and external quinine method calculation on the actually measured aeromagnetic section of the research area to obtain the depth of the magnetic substrate;
in some embodiments, the step S110 includes:
step 1, carrying out tangent method and external quinine method calculation on the aeromagnetic section to obtain a minimum buried depth value of a magnetic body;
based on forward abnormal characteristics of the uniformly magnetized two-degree thick plate-shaped body, the aeromagnetic section before the chemical pole is calculated by a manual tangent method and an external quinine method, and the minimum buried depth point of the underground magnetic body is obtained by an Euler deconvolution automatic calculation method. FIG. 2 is a schematic diagram showing the basis of a tangent method according to an embodiment of the present invention, wherein for the symmetrical magnetic anomaly curves shown in FIG. 2, tangents of two wings of the curve are respectively made, and 0 value lines are respectively crossed at p 1 And p 4 Two points; and then making a horizontal tangent line at the maximum value Mag_Max of the curve, wherein the vertical projection of the intersection point of the horizontal tangent line and the tangent line of two wings of the curve is positioned at p of a 0-value line 2 And p 3 Two points. Meter for minimum buried depth h of magnetic bodyThe calculation formula is as follows:
h=(p 4 -p 3 )/K h , (1)
K h =f[(p 4 -p 1 )/(p 3 -p 2 )], (2)
wherein K is h The statistical table can be calculated by tangent method.
Step 2, removing shallow magnetic rock mass which does not accord with preset rock type information in the minimum buried depth value of the magnetic body;
in the embodiment of the invention, the preset rock type information refers to removing depth points which obviously do not accord with geological facts in the step 1 according to rock types and susceptibility characteristics of different depths confirmed by exploration seismic and drilling data. For example, the minimum burial depth points of the volcanic rock obtained by a tangent method and an external quinine method can be removed through drilling the depth of the sedimentary volcanic rock encountered by drilling. By this method, discrete depth points of the magnetic base of the entire basin can be obtained.
And step 3, performing interpolation processing on the minimum embedded depth value of the removed magnetic body to obtain the magnetic substrate depth.
And (3) obtaining a magnetic substrate depth contour line through a minimum curvature interpolation method according to the discrete depth point of the minimum buried depth of the magnetic body obtained in the step (2), and obtaining the magnetic substrate depth.
S120, obtaining the Curie depth of the research area through a one-dimensional steady-state heat conduction equation based on the ground temperature gradient of the basin in the research area, the heat yield of the crust rock and the heat conductivity of the crust rock;
in some embodiments, step S120 includes:
step 4, constructing the one-dimensional steady-state heat conduction equation based on the ground temperature gradient of the basin in the research area, the heat yield of the crust rock and the heat conductivity of the crust rock, and obtaining the temperature structures of the crust and the top of the rock circle;
the one-dimensional steady-state heat conduction equation describes a state in which heat flowing in and out of the inside of a substance is balanced in one dimension, and the expression is as follows:
where z represents a depth variable, T (z) represents a temperature at the depth z, H 0 Represents the heat yield rate at the shallow part of the crust, h r An index of decay of the thermal yield with depth, k (T (z)) representing the thermal conductivity, is an amount that varies with temperature,the ground temperature gradient is represented, and the change amount of the temperature per unit distance is represented.
The earth surface heat flow value can be obtained through the actually measured heat conductivity and the earth temperature gradient of the well drilling in the basin and the product of the heat conductivity and the earth temperature gradient, the heat flow value is taken as constraint, meanwhile, the heat conductivity and the heat yield obtained through statistics of the well drilling sample, the crust and the rock ring mantle are solved, and the temperature structures of the crust and the rock ring top can be obtained.
Step 5, performing linear interpolation in the vertical direction on the temperature structures at the tops of the crust and the rock ring to obtain Curie point depths of all preset discrete points; and carrying out minimum curvature interpolation on the Curie point depth of each preset discrete point on a plane to obtain the Curie depth of the research area.
Generally, the demagnetizing temperature of pure magnetite is about 580 degrees, so the curie temperature of a typical continental region is set to 580 degrees. And (3) performing linear interpolation in the vertical direction on the temperature structures of the crust and the top of the rock ring generated in the step (4), and obtaining the depth of the Curie temperature of each discrete point, namely the Curie point depth.
The discrete curie point depths are then interpolated with minimal curvature on a plane to obtain the curie point depth, i.e., the curie plane depth, of the entire basin.
S130, under the constraint of the depth of the magnetic substrate and the depth of the Curie face, performing depth weighted magnetic anomaly three-dimensional inversion to obtain the magnetic susceptibility of the basin crust.
In some embodiments, step S130 includes:
step 6, performing pole-changing treatment on the grid magnetic anomaly actually measured in the research area to obtain a variable dip angle pole-changing result of the basin;
grid magnetic anomalies are the basis of inversion of the three-dimensional magnetic structure of the crust, and the pole-changing magnetic anomalies can improve inversion speed based on the geometric grid equivalent compression technology, so pole-changing treatment on the aeromagnetic anomalies is necessary. In order to obtain a more reliable electrode melting result, a point-to-point inclination angle electrode melting method is adopted. The formula of the aeromagnetization pole is as follows:
S(u,v)=S T (u,v)·H(u,v), (4)
α 0 =cos I 0 ·cos D 0 , (6)
β 0 =cos I 0 ·sin D 0 , (7)
γ 0 =sin I 0 , (8)
wherein S (u, v) represents the Fourier transform after the metaplasia, S T (u, v) represents the Fourier transform of the grid magnetic anomaly, H (u, v) represents the Fourier transform of the polarization factor, I represents the imaginary unit, (u, v) represents the frequency domain argument, I 0 Represents geomagnetic inclination angle, D 0 Representing the geomagnetic declination.
The method for changing the dip angle pole is characterized in that Fourier transformation H (u, v) of pole changing factors of corresponding points is calculated according to the geomagnetic dip angle and geomagnetic deflection angle of each point, then the Fourier transformation H (u, v) is brought into a formula (4), inverse Fourier transformation is carried out, the process is circulated, and finally the dip angle pole changing result of the whole basin is obtained.
Step 7, establishing a depth weighted inversion objective function aiming at the detection of the susceptibility structure in the basin;
as an embodiment, the calculation formula of the objective function is as follows:
wherein, the right side of the medium formula in the formula (9) comprises two terms, wherein the first term is a data fitting term; the second term is a regularization term. μ represents a regularization factor for balancing the data fitting term and the regularization term; d represents the observed data, m represents the model parameters, G represents the forward operator of the uniformly magnetized cubic magnetic anomaly, W d Representing a diagonal weighting matrix, the elements of which are the inverse of the estimated observed data noise standard deviation; w (W) m And Z represents a depth weighting matrix which decays exponentially with depth to prevent the inversion model from concentrating at a shallow position.
And 8, taking the depth of the magnetic substrate as an upper interface constraint, taking the depth of the Curie surface as a lower interface constraint, taking the statistical data of the magnetic susceptibility measured in the field as a physical constraint, and carrying out aeromagnetic abnormal depth weighted inversion according to the objective function to obtain the magnetic susceptibility of the basin crust.
And (3) taking the depth of the magnetic substrate obtained in the step (3) as an upper interface constraint, taking the depth of the Curie surface obtained in the step (5) as a lower interface constraint, and carrying out aeromagnetic abnormal depth weighted inversion on the statistical data of the magnetic susceptibility measured in the field as physical property constraint to obtain the magnetic susceptibility distribution of the magnetic layer in the basin.
Specifically, in the inversion process, the magnetic susceptibility of the model units shallower than the depth of the magnetic substrate and the model units deeper than the Curie depth are forcedly set to 0. In the inversion process, if the value of the model unit between the depth of the magnetic substrate and the depth of the Curie surface is smaller than zero or larger than the maximum value of the magnetic susceptibility of the field actually measured ferromagnetic specimen, the magnetic susceptibility is forced to be in a reasonable variation range.
According to the inversion method of the crust magnetic structure based on the depth constraint of the inside and the magnetic substrate, the depth of the magnetic substrate is obtained through calculation by an efficient manual tangent method and an external quinine method; then, based on heat flow data and heat conductivity and heat yield data existing in the basin, obtaining the depth of 580 ℃ inside by solving a one-dimensional steady-state heat conduction equation; and under the constraint of the depth of the magnetic substrate and the depth of the Curie face, performing depth weighted magnetic anomaly three-dimensional inversion to obtain the magnetic susceptibility three-dimensional structure of the basin crust. The method effectively relieves the problem of longitudinal resolution of susceptibility inversion, accurately constrains transverse change of the magnetic layer, provides more reliable results for detecting the crust material structure of the basin, and further provides mineral resource investigation for service science research.
Fig. 3 is a schematic structural diagram of an inversion system of a crust magnetic structure based on depth constraint of an inside and a magnetic substrate according to an embodiment of the present invention, as shown in fig. 3, the system includes a magnetic substrate depth calculation module 310, a curie depth calculation module 320, and a magnetic power calculation module 330, wherein:
the magnetic substrate depth calculation module 310 is configured to perform tangent method and external quinine method calculation on the actually measured aeromagnetic section in the research area to obtain a magnetic substrate depth;
the curie depth calculation module 320 is configured to obtain the curie depth of the investigation region by using a one-dimensional steady-state heat conduction equation based on the ground temperature gradient of the basin in the investigation region, the heat yield of the crust rock, and the heat conductivity of the crust rock;
the magnetic susceptibility calculation module 330 is configured to perform depth weighted three-dimensional inversion of the magnetic anomaly under the constraint of the magnetic base depth and the curie face depth to obtain the magnetic susceptibility of the basin crust.
The implementation process of the system embodiment corresponding to the method embodiment is the same as that of the method embodiment, and reference is made to the method embodiment for details, so that the system embodiment is not repeated.
The various modules in the above-described earth crust magnetic structure inversion system based on the in-house and magnetic substrate depth constraints may be implemented in whole or in part in software, hardware, and combinations thereof. The above modules may be embedded in hardware or may be independent of a processor in the computer device, or may be stored in software in a memory in the computer device, so that the processor may call and execute operations corresponding to the above modules.
Fig. 4 is a schematic structural diagram of a computer device according to an embodiment of the present invention, where the computer device may be a server, and an internal structure diagram of the computer device may be as shown in fig. 4. 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 includes a computer storage medium, an 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 media. The database of the computer device is used for storing data, such as magnetic substrate depth and Curie surface depth, generated or acquired during the inversion method of the crust magnetic structure based on the depth constraint of the inside and the magnetic substrate. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program when executed by the processor implements a method for inversion of a crust magnetic structure based on an in-house and magnetic substrate depth constraint.
In one embodiment, a computer device is provided that includes a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the steps of one of the above embodiments of a method for inversion of a crustal magnetic structure based on an in-house and magnetic substrate depth constraint. Alternatively, the processor, when executing the computer program, performs the functions of the modules/units in this embodiment of the earth's crust magnetic structure inversion system based on the depth constraints of the in-house and magnetic substrates.
In one embodiment, a computer storage medium is provided, on which a computer program is stored, which when executed by a processor, implements the steps of a method for inversion of a crustal magnetic structure based on in-house and magnetic substrate depth constraints in the above embodiment. Alternatively, the computer program, when executed by the processor, performs the functions of the modules/units in this embodiment of an earth's crust magnetic structure inversion system based on the in-house and magnetic substrate depth constraints.
Those skilled in the art will appreciate that implementing all or part of the above described methods may be accomplished by way of a computer program stored on a non-transitory computer readable storage medium, which when executed, may comprise the steps of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the various embodiments provided herein may include non-volatile and/or volatile memory. The nonvolatile 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), memory bus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), among others.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-described division of the functional units and modules is illustrated, and in practical application, the above-described functional distribution may be performed by different functional units and modules according to needs, i.e. the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-described functions.
The above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention, and are intended to be included in the scope of the present invention.
Claims (10)
1. The earth crust magnetic structure inversion method based on the depth constraint of the inside and the magnetic substrate is characterized by comprising the following steps:
performing tangent method and external quinine method calculation on the actually measured aeromagnetic section of the research area to obtain the depth of the magnetic substrate;
obtaining the Curie depth of the research area through a one-dimensional steady-state heat conduction equation based on the ground temperature gradient of the basin in the research area, the heat yield of the crust rock and the heat conductivity of the crust rock;
and under the constraint of the depth of the magnetic substrate and the depth of the Curie face, performing depth weighted magnetic anomaly three-dimensional inversion to obtain the magnetic susceptibility of the basin crust.
2. The inversion method of the crustal magnetic structure based on the constraints of the depth of the inside and the magnetic substrate according to claim 1, wherein the step of obtaining the depth of the magnetic substrate by performing tangential method and external quinine method calculation on the actually measured aeromagnetic section of the research area comprises the following steps:
performing tangent method and external quinine method calculation on the aeromagnetic section to obtain a minimum buried depth value of the magnetic body;
removing shallow magnetic rock mass which does not accord with preset rock type information in the minimum burial depth value of the magnetic body;
and carrying out interpolation treatment on the minimum embedded depth value of the removed magnetic body to obtain the magnetic substrate depth.
3. The method of inversion of a crustal magnetic structure based on the constraints of depth of the inside and magnetic base of claim 1, wherein said step of obtaining the curie depth of said investigation region by one-dimensional steady state thermal conduction equation based on the earth temperature gradient of the basin in said investigation region, the thermal yield of the crustal rock and the thermal conductivity of said crustal rock comprises:
constructing the one-dimensional steady-state heat conduction equation based on the ground temperature gradient of the basin in the research area, the heat yield of the crust rock and the heat conductivity of the crust rock, and obtaining the temperature structures of the crust and the top of the rock circle;
performing vertical linear interpolation on the temperature structures at the tops of the crust and the rock ring to obtain Curie point depths of preset discrete points;
and carrying out minimum curvature interpolation on the Curie point depth of each preset discrete point on a plane to obtain the Curie depth of the research area.
4. The inversion method of the crustal magnetic structure based on the depth constraints of the inside and magnetic substrates according to claim 3, wherein the one-dimensional steady-state heat conduction equation is calculated as follows:
where z represents a depth variable, T (z) represents a temperature at the depth z, H 0 Represents the heat yield rate at the shallow part of the crust, h r Representing the decay index of the thermal yield with depth, k (T (z)) representing the thermal conductivity,representing the ground temperature gradient.
5. The method of inversion of a magnetic structure of the earth crust based on depth constraints of the magnetic substrate and the inside of the earth crust according to claim 1, wherein the step of obtaining the magnetic susceptibility of the earth crust by depth weighted three-dimensional inversion of the magnetic anomaly under the constraints of the depth of the magnetic substrate and the depth of the curie face comprises:
establishing a depth weighted inversion objective function aiming at the detection of the susceptibility structure in the basin;
and taking the depth of the magnetic substrate as an upper interface constraint, the depth of the Curie surface as a lower interface constraint, taking the statistical data of the magnetic susceptibility measured in the field as a physical constraint, and carrying out aeromagnetic anomaly depth weighted inversion according to the objective function to obtain the magnetic susceptibility of the basin crust.
6. The method of inversion of a crustal magnetic structure based on depth constraints of an in-house and magnetic base of claim 5, wherein said objective function is calculated as:
wherein μ represents a regularization factor, d represents observed data, m represents model parameters, G represents a forward operator of uniformly magnetized cubic magnetic anomalies, W d Represents a diagonal weighting matrix, W m Representing a model weighting matrix consisting of the sum of the zero-order and first-order finite differences of the model, Z represents a depth weighting matrix that decays exponentially with depth.
7. The method of inversion of a crustal magnetic structure based on the depth constraints of the inside and magnetic base of claim 5, further comprising, prior to said step of establishing an objective function for depth weighted inversion of susceptibility structure detection in the basin:
and performing pole-changing treatment on the grid magnetic anomaly actually measured in the research area to obtain a variable dip angle pole-changing result of the basin.
8. The inversion method of the crust magnetic structure based on the depth constraint of the inside and the magnetic substrate according to claim 7, wherein the calculation formula of the polarization treatment is as follows:
S(u,v)=S T (u,v)·H(u,v),
α 0 =cosI 0 ·cosD 0 ,
β 0 =cosI 0 ·sinD 0 ,
γ 0 =sinI 0 ,
wherein S (u, v) represents the Fourier transform after the metaplasia, S T (u, v) represents the Fourier transform of the grid magnetic anomaly, H (u, v) represents the Fourier transform of the polarization factor, I represents the imaginary unit, (u, v) represents the frequency domain argument, I 0 Represents geomagnetic inclination angle, D 0 Representing the geomagnetic declination.
9. An inversion system for a magnetic structure of a crust based on depth constraints of an in-house and magnetic substrate, comprising:
the magnetic substrate depth calculation module is used for carrying out tangential method and external quinine method calculation on the actually measured aeromagnetic section of the research area to obtain the magnetic substrate depth;
the Curie depth calculation module is used for obtaining the Curie depth of the research area through a one-dimensional steady-state heat conduction equation based on the ground temperature gradient of the basin in the research area, the heat yield of the crust rock and the heat conductivity of the crust rock;
and the magnetic susceptibility calculation module is used for carrying out depth weighted magnetic anomaly three-dimensional inversion under the constraint of the magnetic substrate depth and the Curie depth to obtain the magnetic susceptibility of the basin crust.
10. A computer storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the earth crust magnetic structure inversion method based on the depth constraints of the in-house and magnetic substrate as claimed in any one of claims 1 to 8.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311675344.8A CN117647847A (en) | 2023-12-07 | 2023-12-07 | Crust magnetic structure inversion method based on depth constraint of inside and magnetic substrate |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311675344.8A CN117647847A (en) | 2023-12-07 | 2023-12-07 | Crust magnetic structure inversion method based on depth constraint of inside and magnetic substrate |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117647847A true CN117647847A (en) | 2024-03-05 |
Family
ID=90047586
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311675344.8A Pending CN117647847A (en) | 2023-12-07 | 2023-12-07 | Crust magnetic structure inversion method based on depth constraint of inside and magnetic substrate |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117647847A (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140365193A1 (en) * | 2011-12-20 | 2014-12-11 | Edward Karl Biegert | Method to constrain a basin model with curie depth |
RU2680262C1 (en) * | 2017-12-07 | 2019-02-19 | Федеральное государственное бюджетное учреждение науки Институт земного магнетизма, ионосферы и распространения радиоволн им. Н.В. Пушкова Российской академии наук (ИЗМИРАН) | Curie isotherm depth of location determining method |
CN110515136A (en) * | 2019-07-03 | 2019-11-29 | 吉林大学 | A kind of terrestrial heat flow estimation method based on boundary surface of gravitation and magnetism inverting |
CN116360000A (en) * | 2023-02-17 | 2023-06-30 | 中国自然资源航空物探遥感中心 | Aeromagnetic detection technology for crust material structure |
CN117148457A (en) * | 2023-08-29 | 2023-12-01 | 长安大学 | Magnetic layer magnetization modulus calculation method and system |
-
2023
- 2023-12-07 CN CN202311675344.8A patent/CN117647847A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140365193A1 (en) * | 2011-12-20 | 2014-12-11 | Edward Karl Biegert | Method to constrain a basin model with curie depth |
RU2680262C1 (en) * | 2017-12-07 | 2019-02-19 | Федеральное государственное бюджетное учреждение науки Институт земного магнетизма, ионосферы и распространения радиоволн им. Н.В. Пушкова Российской академии наук (ИЗМИРАН) | Curie isotherm depth of location determining method |
CN110515136A (en) * | 2019-07-03 | 2019-11-29 | 吉林大学 | A kind of terrestrial heat flow estimation method based on boundary surface of gravitation and magnetism inverting |
CN116360000A (en) * | 2023-02-17 | 2023-06-30 | 中国自然资源航空物探遥感中心 | Aeromagnetic detection technology for crust material structure |
CN117148457A (en) * | 2023-08-29 | 2023-12-01 | 长安大学 | Magnetic layer magnetization modulus calculation method and system |
Non-Patent Citations (1)
Title |
---|
马杰;王万银;罗新刚;姚攀;赵志刚;邱之云;纪晓琳;鲁宝亮;冯旭亮;: "基于磁力资料的南海及邻区磁性基底和居里面深度特征研究", 地球物理学报, no. 10, 15 October 2018 (2018-10-15) * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pilkington et al. | Mitigating remanent magnetization effects in magnetic data using the normalized source strength | |
Miensopust et al. | Lithospheric structures and Precambrian terrane boundaries in northeastern Botswana revealed through magnetotelluric profiling as part of the Southern African Magnetotelluric Experiment | |
Stavrev et al. | Magnetic field transforms with low sensitivity to the direction of source magnetization and high centricity [Link] | |
CN108508490B (en) | Magnetic tensor gradient data equilibrium boundary identification method based on analytic signals | |
Gerovska et al. | Estimating the magnetization direction of sources from southeast Bulgaria through correlation between reduced‐to‐the‐pole and total magnitude anomalies | |
CN106405664B (en) | A kind of magnetic anomaly normalizing pole method | |
CN113255230B (en) | Gravity model forward modeling method and system based on MQ radial basis function | |
CN111856599B (en) | Magnetic measurement data equivalent source pole and type conversion method based on PDE | |
CN110687610B (en) | Gravity and magnetic data correlation analysis-based field source positioning and attribute identification method | |
CN111337992A (en) | Method for obtaining depth of field source based on downward continuation of bit field data | |
Fedi et al. | Decorrugation and removal of directional trends of magnetic fields by the wavelet transform: application to archaeological areas | |
CN110346835A (en) | Magnetotelluric forward modeling method, forward modeling system, storage medium and electronic equipment | |
Yang et al. | 3D gravity inversion with optimized mesh based on edge and center anomaly detection | |
Li et al. | A paradigm shift in magnetic data interpretation: Increased value through magnetization inversions | |
CN117647847A (en) | Crust magnetic structure inversion method based on depth constraint of inside and magnetic substrate | |
Rybin et al. | Deep Geoelectric Structure of the Earth’s Crust and the Upper Mantle of the Pamir–Alai Zone | |
CN112684498A (en) | Reservoir fracture prediction method and system based on wide-azimuth seismic data | |
Yang et al. | Inversion of high-amplitude magnetic total field anomaly: an application to the Mengku iron-ore deposit, northwest China | |
CN113075749B (en) | Method and system for locating favorable ore-forming space of sandstone uranium deposit in anticline development area | |
CN112328955B (en) | Method for processing heavy magnetic data, storage medium and device | |
Smith et al. | An enhanced method for source parameter imaging of magnetic data collected for mineral exploration | |
de Souza et al. | Using parity decomposition for interpreting magnetic anomalies from dikes having arbitrary dip angles, induced and remanent magnetization | |
CN112651102A (en) | Bit field full-automatic extreme point depth estimation method under undulating terrain | |
Chen et al. | Potential field data interpolation by Taylor series expansion | |
Fang et al. | A new axial constraint with variable angle for the inversion of gravity data |
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 |