Geomagnetic field modeling method, device, medium and computer program product

Abstract

The invention discloses a geomagnetic field modeling method and device, a medium and a computer program product, and relates to the technical field of geomagnetic field modeling, and the geomagnetic field modeling method comprises the steps: obtaining a scale factor matrix of a spherical basic current system; obtaining a first geomagnetic effect conversion matrix of the interpolation point; determining a magnetic field matrix of the interpolation point according to the scale factor matrix and the first geomagnetic effect conversion matrix; and constructing a geomagnetic field model based on the magnetic field matrix of the interpolation points. By constructing the geomagnetic field model, the magnetic field information of any point under the ionosphere in the global range can be determined, so that the blank of magnetic field data is filled, and the coverage range of global magnetic field data is expanded.

Description

Technical Field

[0001] The present application relates to the technical field of geomagnetic field modeling, and in particular to geomagnetic field modeling methods, devices, media and computer program products. Background Art

[0002] As one of the important physical fields of the Earth, the geomagnetic field has a significant impact on the Earth's environment and human activities. Initially, field observations were mainly carried out at geomagnetic observatories, and the distribution and change patterns of the geomagnetic field were analyzed through long-term accumulated data. With the rise of computer technology and numerical simulation methods, mathematical methods such as global magnetohydrodynamics, spherical harmonic analysis, and spherical cap harmonic analysis have gradually become important tools for geomagnetic field modeling.

[0003] At present, the spatial coverage of geomagnetic station observation data is insufficient, and most areas do not have geomagnetic stations, which makes it difficult to obtain high-precision magnetic field data on a global scale. Although the global magnetohydrodynamic model can make up for the lack of geomagnetic station observation data to a certain extent, it still has many limitations. For example: some complex physical processes are difficult to describe with a single global magnetohydrodynamic model, resulting in the model being unable to accurately predict changes in the geomagnetic field in some cases; spherical harmonic analysis requires a large number of data points to be used for global magnetic field analysis. These limitations will lead to reduced accuracy of global magnetic field data. Summary of the invention

[0004] The main purpose of this application is to provide a geomagnetic field modeling method, device, medium and computer program product, aiming to solve the technical problem of how to improve the accuracy of global magnetic field data.

[0005] To achieve the above purpose, the present application proposes a geomagnetic field modeling method, the method comprising: Get the scale factor matrix of the spherical basic current system SECS; Obtaining the first geomagnetic effect conversion matrix of the interpolation point; Determining the magnetic field matrix of the interpolation point according to the scale factor matrix and the first geomagnetic effect conversion matrix; A geomagnetic field model is constructed based on the magnetic field matrix of the interpolation points.

[0006] In one embodiment, the step of obtaining the first geomagnetic effect conversion matrix of the interpolation point includes: The radial component and the zenith component of the magnetic field at the interpolation point are acquired, and the first geomagnetic effect conversion matrix of the interpolation point is determined according to the radial component and the zenith component of the magnetic field at the interpolation point.

[0007] In one embodiment, the step of obtaining the radial component and the zenith component of the magnetic field of the interpolation point comprises: Acquire the angle information between the interpolation point and the pole of the SECS and the radial distance of the interpolation point; The radial component and the zenith component of the magnetic field of the interpolation point are determined according to the angle information and the radial distance.

[0008] In one embodiment, the step of obtaining the radial component and the zenith component of the magnetic field of the interpolation point comprises: Obtaining a non-divergent current density vector of the SECS and a radial distance between the interpolation points; The radial component and the zenith component of the magnetic field at the interpolation point are determined according to the unkempt current density vector and the radial distance.

[0009] In one embodiment, the step of obtaining a scale factor matrix of a spherical elementary current system SECS comprises: Obtain the second geomagnetic effect conversion matrix of the observation point; Obtaining a magnetic field disturbance vector of the observation point; A scale factor matrix of the SECS is determined according to the second geomagnetic effect conversion matrix and the magnetic field disturbance vector.

[0010] In one embodiment, the step of obtaining the second geomagnetic effect conversion matrix of the observation point includes: The radial component and the zenith component of the magnetic field at the observation point are obtained, and the second geomagnetic effect conversion matrix of the observation point is determined according to the radial component and the zenith component of the observation point.

[0011] In one embodiment, before the step of obtaining the scale factor matrix of the spherical basic current system SECS, the method further includes: Determine the grid area based on the distribution information of the observation points; On the grid area, the poles of the SECS are determined.

[0012] In addition, to achieve the above purpose, the present application also proposes a geomagnetic field modeling device, the geomagnetic field modeling device comprising: A scale factor acquisition module is used to obtain the scale factor matrix of the spherical basic current system SECS; A conversion matrix acquisition module, used to obtain a first geomagnetic effect conversion matrix of an interpolation point; A magnetic field matrix determination module, used to determine the magnetic field matrix of the interpolation point according to the scale factor matrix and the first geomagnetic effect conversion matrix; The model building module is used to build a geomagnetic field model based on the magnetic field matrix of the interpolation points.

[0013] In addition, to achieve the above-mentioned purpose, the present application also proposes a geomagnetic field modeling device, which includes: a memory, a processor, and a computer program stored in the memory and executable on the processor, and the computer program is configured to implement the steps of the geomagnetic field modeling method as described above.

[0014] In addition, to achieve the above-mentioned purpose, the present application also proposes a storage medium, which is a computer-readable storage medium, and a computer program is stored on the storage medium. When the computer program is executed by a processor, the steps of the geomagnetic field modeling method described above are implemented.

[0015] In addition, to achieve the above-mentioned purpose, the present application also provides a computer program product, which includes a computer program, and when the computer program is executed by a processor, it implements the steps of the geomagnetic field modeling method as described above.

[0016] One or more technical solutions proposed in the present application have at least the following technical effects: by obtaining the scale factor matrix of the spherical basic current system SECS, the degree of influence of the current source in SECS on the magnetic field of the interpolation point can be determined, which provides a basis for subsequent magnetic field calculations and helps to accurately simulate and predict changes in the earth's magnetic field; by obtaining the first geomagnetic effect conversion matrix of the interpolation point, the geomagnetic effect of the interpolation point can be calculated according to the relative position of SECS and the interpolation point under a fixed current coefficient of SECS, including key information such as magnetic field strength and direction, which provides accurate magnetic field information for subsequent determination of the magnetic field matrix; and then, the magnetic field matrix of the interpolation point can be determined according to the scale factor matrix and the first geomagnetic effect conversion matrix, and the current distribution of SECS and the geomagnetic effect of the interpolation point are linked, thereby realizing accurate conversion from current distribution to magnetic field distribution of the interpolation point, making the magnetic field calculation result of the interpolation point more accurate; based on the magnetic field matrix of the interpolation point, a geomagnetic field model is constructed, and the magnetic field of the unknown point can be predicted through the geomagnetic field model, thereby accurately reflecting the distribution of the earth's magnetic field. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0018] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, for ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0019] Figure 1 A schematic diagram of a flow chart provided for the first embodiment of the geomagnetic field modeling method of the present application; Figure 2 A schematic diagram of a spherical basic current system without dispersion provided in Example 1 of the present application; Figure 3 A schematic diagram of the module structure of the geomagnetic field modeling program provided in Example 2 of the present application; Figure 4 A schematic diagram of the change of magnetic field data over time provided in Example 2 of the present application; Figure 5 This is a schematic diagram of the module structure of the geomagnetic field modeling device according to an embodiment of the present application; Figure 6 Schematic diagram of the device structure of the hardware operating environment involved in the geomagnetic field modeling method in the embodiment of the present application. DETAILED DESCRIPTION

[0020] It should be understood that the specific embodiments described herein are only used to explain the technical solutions of the present application and are not used to limit the present application.

[0021] In order to better understand the technical solution of the present application, a detailed description will be given below in conjunction with the accompanying drawings and specific implementation methods.

[0022] In this embodiment, for ease of description, the following description is made with the geomagnetic field modeling device as the execution subject.

[0023] Because in conventional technology, global magnetohydrodynamic models also rely on external input data such as space weather parameters, the accuracy and availability of these data also limit the scope of application of the model. Although spherical harmonic analysis can be used for global magnetic field analysis, it requires a large number of data points and cannot meet the resolution of local analysis. Spherical cap harmonic analysis can be better used for local magnetic field analysis, but the calculation amount is larger and the calculation cost is higher. These problems limit the large-area coverage of global magnetic field data and the improvement of magnetic field accuracy.

[0024] The present application provides a solution, which obtains the scale factor matrix of SECS through the relevant magnetic field data of the observation point, so that the magnetic field strength of different regions can be calibrated, which effectively reduces the error introduced by the uneven distribution of the earth's magnetic field data and improves the accuracy of the overall model; further, according to the angular information of the interpolation point and the pole of SECS and the radial distance of the interpolation point, the radial component and zenith component of the magnetic field of the interpolation point are determined, so as to construct the first geomagnetic effect conversion matrix of the interpolation point, that is, according to the current distribution of SECS and the relative position of SECS and the interpolation point, the magnetic field component intensity of the interpolation point is determined, which provides accurate magnetic field distribution information for the subsequent determination of the magnetic field matrix; the determination of the radial component and zenith component of the magnetic field of the interpolation point can also be based on SEC The radial distance between the interpolation point and the divergence-free current density vector of S is determined, which is essentially to determine the size of the magnetic field component of the interpolation point based on the relative position and the current distribution of SECS; then, the magnetic field matrix of the interpolation point can be determined based on the scale factor matrix and the first geomagnetic effect conversion matrix. Due to the different relative positions of SECS and the interpolation points, the influence of the geomagnetic effect of the interpolation point inferred based on the current distribution of SECS will be different. Therefore, based on the spatial relationship between SECS and the interpolation points, the change of the geomagnetic effect caused by the current distribution of SECS is reasonably inferred, ensuring the accurate evaluation of the geomagnetic effect of the interpolation point; based on the magnetic field matrix of the interpolation point, a geomagnetic field model is constructed. Through the geomagnetic field model, the magnetic field of the unknown point can be predicted, thereby accurately reflecting the distribution of the earth's magnetic field.

[0025] It should be noted that the execution subject of this embodiment can be a computing service device with data processing, network communication and program running functions, such as a tablet computer, a personal computer, a mobile phone, etc., or an electronic device capable of realizing the above functions, etc. The following takes the geomagnetic field modeling device as an example of the execution subject to illustrate this embodiment and the following embodiments.

[0026] Based on this, the present application embodiment provides a method for modeling the geomagnetic field, referring to Figure 1 , Figure 1 FIG. 1 is a flow chart of a first embodiment of a local magnetic field modeling method.

[0027] In this embodiment, the geomagnetic field modeling method includes steps S10 to S40: Step S10, obtaining a scale factor matrix of a spherical basic current system SECS; It should be noted that SECS (Spherical Elementary Current Systems) is a mathematical method used to represent and analyze the Earth's magnetic field and ionospheric current system, which can represent the divergence and curl of the current. The poles of SECS can be any point in the specified area. Since the specified area can be divided into several regular or irregular grids, it is ensured that the poles of SECS fall on the grid points. The poles of SECS can also be called grid points; in other words, after the specified area is gridded, the grid points are determined as the poles of SECS to establish a spherical elementary current system.

[0028] In addition, it should be noted that the scale factor is a coefficient factor that reflects the size of the spherical basic current system and is used to determine the degree of influence of the spherical basic current system on the magnetic field of the interpolation point.

[0029] For example, the specific form of the scale factor matrix can refer to the following formula:

[0030] The I in this formula represents the scale factor matrix, represents the scale factor of the lth SECS, Indicates the total number of SECS poles, that is, the total number of grid points. In a custom grid area, you can use The current distribution of each SECS is used to infer the magnetic field distribution at the interpolation point.

[0031] Step S20, obtaining a first geomagnetic effect conversion matrix of the interpolation point; It should be noted that the interpolation point refers to a virtual point introduced to estimate or predict the unknown data value located near the known data point, which is mainly used to fill the gaps in the magnetic field data; the geomagnetic effect conversion matrix is ​​mainly used to display the intensity of the magnetic field components of the interpolation point in the radial direction and the zenith direction, that is, when the SECS scale factor is a fixed value (such as 1A), the magnetic field intensity information of the interpolation point in all directions.

[0032] For example, please refer to the following formula:

[0033] In this formula, T represents the first geomagnetic effect conversion matrix, T k,l,θ It means that when I is a fixed value (such as 1A), the disturbance of the zenith angle direction of the kth point due to the lth SECS is k,l,φ It means that when I is a fixed value, the disturbance of the azimuth direction of the kth point due to the lth SECS is kRepresents the total number of interpolation points. The points to be sought here are interpolation points. By performing mathematical calculations on the magnetic field disturbances in all directions of the interpolation points, the magnetic field data of the interpolation points in all directions can be obtained.

[0034] Step S30, determining the magnetic field matrix of the interpolation point according to the scale factor matrix and the first geomagnetic effect conversion matrix; It should be noted that the magnetic field matrix represents the state of the geomagnetic field at a specific location (interpolation point). The matrix contains information such as the magnitude of the magnetic field component at that location and can be obtained by performing matrix multiplication operations on the scale factor matrix and the first geomagnetic effect conversion matrix.

[0035] For example, the formula To determine the magnetic field matrix Z of the interpolation point, where T contains the magnetic field disturbance value of the interpolation point when the scale factor is a fixed value, which is generally set to 1A. However, due to the different sizes of each SECS and the relative positions of its poles (grid points) to the interpolation point, the geomagnetic effects of the interpolation points inferred from the current distribution of the SECS will also be different. It is necessary to determine the size of the SECS of each grid point based on the scale factor matrix, thereby determining the direction and size of the magnetic field at the interpolation point.

[0036] It can be understood that by combining the scale factor matrix and the first geomagnetic effect conversion matrix, the magnetic field of the interpolation point can be predicted more accurately and the gaps in the spatial magnetic field data can be filled.

[0037] Step S40, constructing a geomagnetic field model based on the magnetic field matrix of the interpolation points.

[0038] It should be noted that the geomagnetic field model refers to a mathematical model used to simulate the magnetic field data of any point under the ionosphere on the earth. It can describe and predict the spatial distribution of the magnetic field and the changes of the magnetic field over time. The model usually includes a series of equations and parameters to simulate the various characteristics of the geomagnetic field.

[0039] Exemplarily, based on the magnetic field matrix of the interpolation point, the method for obtaining the magnetic field matrix (parametric equation) is clarified, and the geomagnetic field model is determined according to the data of the magnetic field matrix. Specifically, in order to obtain the magnetic field matrix of the interpolation point, the geomagnetic field model needs to calculate the scale factor matrix and the geomagnetic effect conversion matrix of the interpolation point according to the parametric equation; the scale factor matrix can be obtained by obtaining the magnetic field data of the observation point and the geomagnetic effect conversion matrix of the observation point, and then performing singular value decomposition on the parametric equation; and the geomagnetic effect conversion matrix of the interpolation point and the observation point can be determined after calculation based on the relative position of the required point (interpolation point or observation point) and the grid point and related formulas.

[0040] This embodiment provides a geomagnetic field modeling method. By constructing a geomagnetic field model, the geomagnetic field of all points under the ionosphere in a region can be calculated, thereby effectively filling the gaps in magnetic field data and improving the coverage of global magnetic field data; at the same time, by combining the scale factor matrix and the first geomagnetic effect conversion matrix, the current distribution of SECS is linked to the geomagnetic effect of the interpolation point, which is conducive to constructing an accurate geomagnetic field model, so that the magnetic field data of the interpolation point predicted based on the geomagnetic field model is more accurate.

[0041] In a feasible implementation manner, step S20 includes: Step S21, obtaining the radial component and the zenith component of the interpolation point, and determining the first geomagnetic effect conversion matrix of the interpolation point according to the radial component and the zenith component of the interpolation point.

[0042] It should be noted that the radial component refers to the projection of the magnetic field vector in the radial direction, and the radial direction refers to the direction of the line connecting the interpolation point to the center of the earth; the zenith component refers to the projection of the magnetic field vector in the zenith direction, and the zenith direction refers to the direction perpendicular to the line connecting the interpolation point to the center of the earth.

[0043] From the radial component and zenith component to the first geomagnetic effect conversion matrix, the coordinate system is mainly converted. The radial component and zenith component are converted from the spherical basic current system to the geographic spherical coordinate system using the preset or calculated conversion matrix to display the magnetic field data more intuitively. The converted radial component and zenith component can be used to construct the first geomagnetic effect conversion matrix.

[0044] In this implementation, the radial component and the zenith component are converted into magnetic field data of a geographic spherical coordinate system through coordinate system conversion, which makes it easier to compare and uniformly display global magnetic field data.

[0045] In a feasible implementation manner, the step of obtaining the radial component and the zenith component of the magnetic field of the interpolation point in step S21 includes: Step A211, obtaining the angle information between the interpolation point and the pole of SECS and the radial distance of the interpolation point; It should be noted that the angle information refers to the angle between the line connecting the interpolation point and the center of the earth and the line connecting the SECS pole and the center of the earth; the radial distance refers to the straight-line distance from the interpolation point to the center of the earth. The relative position relationship between the interpolation point and the SECS pole can be determined by combining the radial distance and angle information.

[0046] Step A212, determining the radial component and the zenith component of the magnetic field at the interpolation point according to the angle information and the radial distance.

[0047] Exemplarily, the radial component of the magnetic field may be determined by referring to the following formula:

[0048] In this formula, represents the radial component of the magnetic field at the interpolation point, r represents the radial distance, Indicates angle information. represents the reference radius of the ionosphere, represents the vacuum permeability, Indicates the scale factor of SECS. The size of the radial component is mainly related to the relative position (angle information and radial distance) between the interpolation point and the SECS pole.

[0049] Exemplarily, the zenith component of the magnetic field may be determined with reference to the following formula:

[0050] In this formula, It represents the component of the magnetic field at the interpolation point in the zenith direction (zenith component). Its magnitude is mainly related to the relative position between the interpolation point and the pole of SECS. The meaning of other characters is the same as the formula for determining the radial component, so they will not be repeated here.

[0051] In this embodiment, by determining the radial component and the zenith component in the spherical basic coordinate system, the magnetic field disturbance in the azimuth and zenith angle directions of the interpolation point inferred according to the current distribution of SECS is accurately obtained, thereby improving the accuracy of the magnetic field data.

[0052] In a feasible implementation manner, the step of obtaining the radial component and the zenith component of the magnetic field of the interpolation point in step S21 includes: Step B211, obtaining the divergence-free current density vector of the SECS and the radial distance of the interpolation point; It should be noted that SECS can represent the divergence and curl of the current. Since the curl-free part of the current does not produce any magnetic effect below the ionosphere, this property can be used to determine the strength of the magnetic field based on the divergence-free current below the ionosphere.

[0053] The divergence-free current density vector of SECS can be determined by analyzing the current density distribution of all points on the SECS sphere. Specifically, the divergence-free current density vector of any point on the SECS sphere can be calculated by considering the angular relationship between the point and the interpolation point and applying the corresponding physical formula.

[0054] For example, please refer to Figure 2 , Figure 2A schematic diagram of a spherical basic current system without divergence is provided. The circle represents the spherical basic current system of a certain SECS pole, the arrow represents the divergence-free current density vector, and the thickness of the arrow represents the magnitude of the divergence-free current density vector. The thicker the arrow, the larger the divergence-free current. The method of determining the divergence-free current density vector of a spherical point based on the angle relationship between the SECS spherical point and the SECS pole can refer to the following formula:

[0055] In this formula, represents the undissipated current density vector, It represents the scale factor of the basic system. It is usually set to a fixed value (such as 1A) when calculating the radial component and the zenith component. represents the radial component of the spherical basic current coordinate system, represents the azimuthal unit vector in the spherical fundamental current coordinate system.

[0056] Step B212, determining the radial component and the zenith component of the magnetic field at the interpolation point according to the undeviated current density vector and the radial distance.

[0057] For example, mathematical operations such as applying Maxwell's equations can be used to perform surface integration on information such as the divergence-free current density vector of the SECS and the first radial distance of the interpolation point to determine the radial component and zenith component of the magnetic field at the interpolation point, thereby further determining the magnetic field distribution at the interpolation point.

[0058] In this implementation, by converting the current data into magnetic field data, the magnitude of the magnetic field component of the interpolation point is essentially determined according to the relative position and the current distribution of the SECS, so that the accuracy of the magnetic field data can be improved.

[0059] Based on the first embodiment of the present application, in the second embodiment of the present application, the same or similar contents as those in the above-mentioned embodiment 1 can be referred to the above introduction, and will not be repeated in the following. On this basis, step S10 includes: Step S11, obtaining a second geomagnetic effect conversion matrix of the observation point; It should be noted that the second geomagnetic effect conversion matrix is ​​the magnetic field component data of the observation point determined according to the relative position between the pole of SECS and the observation point.

[0060] Exemplarily, the second geomagnetic effect conversion matrix is ​​composed of the geomagnetic disturbance in the zenith angle direction and the geomagnetic disturbance in the azimuth angle direction, and can refer to the following formula:

[0061] In this formula, T' represents the second geomagnetic effect conversion matrix, the point to be calculated here is the observation point, n obs Indicates the number of observation points.

[0062] It can be understood that by obtaining the second geomagnetic effect conversion matrix of the SECS poles, the change of the geomagnetic effect can be predicted more accurately, so that the magnetic field strength of the interpolation point in different directions can be predicted according to the positional relationship between the SECS poles and the interpolation points.

[0063] Step S12, obtaining the magnetic field disturbance vector of the observation point; It should be noted that the observation point refers to a location with a geomagnetic observatory that can directly measure the local magnetic field data. Figure 1 , Figure 1 A distribution map of observation points in Europe is provided; the magnetic field disturbance vector is a physical quantity that describes the change of the magnetic field. It indicates the degree of deviation of parameters such as magnetic field intensity and direction from a certain reference state or average state at a certain observation point, and usually includes two components: intensity and direction.

[0064] For example, please refer to the following formula:

[0065] In this formula, Z' represents the magnetic field disturbance vector, Z k,θ represents the magnetic field disturbance in the zenith angle direction measured at the kth observation point, Z k,φ It represents the magnetic field disturbance in the azimuth direction measured at the kth observation point. These magnetic field disturbance data can be obtained by reading the relevant files of the observation station, obtaining the latitude and longitude information and magnetic field data information of the observation point, and after filtering, unit conversion, removing linear trends and frequency filtering of the magnetic field data information, the magnetic field disturbance data of each observation point in the specified direction can be obtained, thereby obtaining the magnetic field disturbance vector.

[0066] Step S13, determining the scale factor matrix of SECS according to the second geomagnetic effect conversion matrix and the magnetic field disturbance vector.

[0067] It should be noted that due to , the scale factor matrix can be obtained by using methods such as singular value decomposition.

[0068] In this embodiment, by using the second geomagnetic effect conversion matrix and the magnetic field disturbance vector, the scale factor matrix of the appropriate SECS can be determined, that is, a suitable scale factor is assigned to each SECS, so as to better infer the magnetic field data of the interpolation point according to the current distribution of the SECS and improve the accuracy of the magnetic field data.

[0069] In a feasible implementation, step S11 includes: Step S111, obtaining the radial component and the zenith component of the magnetic field at the observation point, and determining the second geomagnetic effect conversion matrix of the observation point according to the radial component and the zenith component of the observation point.

[0070] The process of determining the second geomagnetic effect conversion matrix is ​​the same as that of determining the first geomagnetic effect conversion matrix. The main difference between the two lies in the different points to be determined. The first geomagnetic effect conversion matrix is ​​mainly performed for interpolation points, while the second geomagnetic effect conversion matrix is ​​mainly performed for observation points.

[0071] Similarly, there are at least two ways to obtain the radial component and zenith component of the magnetic field at the observation point. The first way to obtain is to determine it based on the angle information between the observation point and the pole of SECS (the angle between the line connecting the observation point and the center of the earth and the line connecting the pole of SECS and the center of the earth) and the radial distance of the observation point (the straight-line distance from the observation point to the center of the earth). The specific determination method is the same as the determination method in step A212, which will not be repeated here; the second way to obtain is to determine it based on the non-divergent current density vector of SECS and the radial distance of the observation point. The specific determination method is the same as the determination method in step B212, which will not be repeated here. In addition, the process of determining the second geomagnetic effect conversion matrix of the observation point based on the radial component and the zenith component of the observation point is the same as step S21, so the determination process of the second geomagnetic effect conversion matrix will not be repeated.

[0072] In a feasible implementation manner, before step S10, the method further includes: Step S01, determining a grid area according to distribution information of observation points; It should be noted that distribution information refers to statistical characteristics such as the positional relationship, density, and coverage of observation points on the ground; the grid area refers to the target analysis area. In order to facilitate analysis and research, the target analysis area can be divided into several regular or irregular small areas based on actual needs and the distribution of observation points.

[0073] For example, in areas where observation stations (observation points) are densely populated, the grid area can be appropriately expanded. Even if the grid area is large, the magnetic field data of the interpolation points can be accurately predicted.

[0074] Step S02, determining the poles of SECS on the grid area.

[0075] For example, the grid points can be directly determined as the poles of SECS, which means that each grid point will become a reference point for magnetic field data calculation or interpolation; or suitable positions can be selected as SECS poles at any location other than the grid points according to the specific magnetic field data distribution characteristics and analysis requirements, so as to facilitate subsequent magnetic field data interpolation calculations.

[0076] In the global and local modeling of the geomagnetic field, by determining the poles of SECS, the distribution and changes of the geomagnetic field can be simulated more accurately; and by determining the poles of SECS in a divided grid area, the magnetic field information of the interpolation point can be better predicted based on the magnetic field information of the observation point.

[0077] In this implementation, by reasonably dividing the grid area, the magnetic field data of the observation point can be better understood and utilized, thereby improving the accuracy of the calculation of the magnetic field data of the unknown point (interpolation point).

[0078] An embodiment of the present application also provides a computer program product, including a computer program, which implements the steps of the above-mentioned geomagnetic field modeling method when executed by a processor.

[0079] The computer program provided by the embodiment of the present application may include: The calculation module is used to schedule each submodule to realize magnetic field simulation; Grid module, used to determine grid area and perform interpolation calculations; A visualization module, used to provide visualization information for the calculation module; The data management module is used to save data information and provide data information to the computing module and the network module.

[0080] It should be noted that visualization information refers to the picture format information of visual presentation, such as displaying magnetic field data in the form of line graphs, vector fields, etc.; data information refers to the basic information required to realize magnetic field simulation and the calculated result data for visualization display, including the data required to complete the interpolation calculation and the calculated magnetic field data.

[0081] For example, please refer to Figure 3 , Figure 3A schematic diagram of the module structure of a geomagnetic field modeling program is provided. The geomagnetic field modeling program 100 includes a calculation module 101, a network module 102, a visualization module 103 and a data management module 104, and realizes efficient simulation and analysis of the geomagnetic field through modular design. The calculation module 101 is the core of the program, which is mainly responsible for scheduling other submodules and integrating various data information, such as visualization image parameter information and magnetic field calculation results, so as to complete the complete display of the simulated magnetic field. The network module 102 is responsible for determining the grid area and performing interpolation calculations. When performing interpolation calculations, the geographic information and common parameters in the data management module 104 will be read, such as the coordinate information of the interpolation point and the grid point, as well as the common parameters such as vacuum permeability and the reference radius of the ionosphere, and the calculated related magnetic field data will be sent to the data management module for storage. The visualization module 103 mainly provides visualization information for the calculation module 101, including data visualization and geographic location visualization. Data visualization refers to the display of the trajectory of the magnetic field data of the interpolation point obtained with a certain factor in a graphical manner, such as a line graph, etc.; geographic location visualization refers to the display of the magnetic field data of all interpolation points at a specific time in the form of longitude and latitude or other geographic locations. The data management module 104 is responsible for storing data information such as geographic location information, public parameters, and calculation results, and retrieves and provides corresponding data according to the requirements of the calculation module and the grid module.

[0082] For example, please refer to Figure 4 , Figure 4 Provides a line diagram of the magnetic field data of the interpolation point changing with time. Figure 4 The horizontal axis of the two line graphs in the figure represents time in minutes; the vertical axis represents the magnetic field strength in nT (nanometer); Mx and My represent the components of the magnetic field data of the interpolation point in different directions. During a magnetic storm, the magnetic field disturbance of the interpolation point will change all the time, and the strength of the magnetic field in different directions will also be affected, and will change differently with time and the strength of the magnetic storm.

[0083] In this embodiment, by constructing a geomagnetic field modeling program, information such as the distribution and intensity of the geomagnetic field can be obtained more accurately. After obtaining the position information of the interpolation point, it is possible to respond quickly to it, calculate and display the magnetic field data in a visual form, thereby completing the simulation and prediction of changes in the Earth's magnetic field.

[0084] It should be noted that the above examples are only used to understand the present application and do not constitute a limitation on the geomagnetic field modeling method of the present application. More forms of simple transformations based on this technical concept are all within the scope of protection of the present application.

[0085] The present application also provides a geomagnetic field modeling device, please refer to Figure 5 , the geomagnetic field modeling device comprises: A scale factor acquisition module 10 is used to acquire a scale factor matrix of a spherical basic current system SECS; A conversion matrix acquisition module 20, used to acquire a first geomagnetic effect conversion matrix of an interpolation point; A magnetic field matrix determination module 30, used to determine the magnetic field matrix of the interpolation point according to the scale factor matrix and the first geomagnetic effect conversion matrix; The model building module 40 is used to build a geomagnetic field model based on the magnetic field matrix of the interpolation points.

[0086] The geomagnetic field modeling device provided in the embodiment of the present application adopts the geomagnetic field modeling method in the above embodiment, which can solve the technical problem of how to improve the accuracy of global magnetic field data. Compared with the prior art, the beneficial effects of the geomagnetic field modeling device provided in the present application are the same as the beneficial effects of the geomagnetic field modeling method provided in the above embodiment, and the other technical features in the geomagnetic field modeling device are the same as the features disclosed in the above embodiment method, which will not be repeated here.

[0087] An embodiment of the present application provides a geomagnetic field modeling device, which includes: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the geomagnetic field modeling method in the above-mentioned embodiment one.

[0088] Reference below Figure 6 , which shows a schematic diagram of the structure of a geomagnetic field modeling device suitable for implementing the embodiment of the present application. The geomagnetic field modeling device in the embodiment of the present application may include but is not limited to mobile terminals such as mobile phones, laptop computers, digital broadcast receivers, PDAs (Personal Digital Assistants), PADs (Portable Application Descriptions), PMPs (Portable Media Players), vehicle-mounted terminals (such as vehicle-mounted navigation terminals), etc., and fixed terminals such as digital TVs, desktop computers, etc. Figure 6 The geomagnetic field modeling device shown is merely an example and should not bring any limitation to the functions and scope of use of the embodiments of the present application.

[0089] like Figure 6As shown, the geomagnetic field modeling device may include a processing device 1001 (such as a central processing unit, a graphics processor, etc.), which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM: Read Only Memory) 1002 or a program loaded from a storage device 1003 to a random access memory (RAM: Random Access Memory) 1004. In RAM1004, various programs and data required for the operation of the geomagnetic field modeling device are also stored. The processing device 1001, ROM1002 and RAM1004 are connected to each other through a bus 1005. An input / output (I / O) interface 1006 is also connected to the bus. Generally, the following systems can be connected to the I / O interface 1006: an input device 1007 including, for example, a touch screen, a touch pad, a keyboard, a mouse, an image sensor, a microphone, an accelerometer, a gyroscope, etc.; an output device 1008 including, for example, a liquid crystal display (LCD: Liquid Crystal Display), a speaker, a vibrator, etc.; a storage device 1003 including, for example, a magnetic tape, a hard disk, etc.; and a communication device 1009. The communication device 1009 can allow the geomagnetic field modeling device to communicate wirelessly or wired with other devices to exchange data. Although the geomagnetic field modeling device with various systems is shown in the figure, it should be understood that it is not required to implement or have all the systems shown. More or fewer systems can be implemented or provided alternatively.

[0090] In particular, according to the embodiments disclosed in the present application, the process described above with reference to the flowchart can be implemented as a computer software program. For example, the embodiments disclosed in the present application include a computer program product, which includes a computer program carried on a computer-readable medium, and the computer program includes a program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from a network through a communication device, or installed from a storage device 1003, or installed from a ROM 1002. When the computer program is executed by the processing device 1001, the above-mentioned functions defined in the method of the embodiment disclosed in the present application are executed.

[0091] The geomagnetic field modeling device provided in the embodiment of the present application adopts the geomagnetic field modeling method in the above embodiment, which can solve the technical problem of how to improve the accuracy of global magnetic field data. Compared with the prior art, the beneficial effects of the geomagnetic field modeling device provided in the present application are the same as the beneficial effects of the geomagnetic field modeling method provided in the above embodiment, and the other technical features in the geomagnetic field modeling device are the same as the features disclosed in the method of the previous embodiment, which will not be repeated here.

[0092] It should be understood that the various parts disclosed in this application can be implemented by hardware, software, firmware or a combination thereof. In the description of the above embodiments, specific features, structures, materials or characteristics can be combined in any one or more embodiments or examples in a suitable manner.

[0093] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art who is familiar with the present technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

[0094] An embodiment of the present application provides a computer-readable storage medium having computer-readable program instructions (ie, a computer program) stored thereon, and the computer-readable program instructions are used to execute the geomagnetic field modeling method in the above-mentioned embodiment.

[0095] The computer-readable storage medium provided in the embodiment of the present application may be, for example, a USB flash drive, but is not limited to electrical, magnetic, optical, electromagnetic, infrared, or semiconductor systems or devices, or any combination of the above. More specific examples of computer-readable storage media may include, but are not limited to: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM: Random Access Memory), a read-only memory (ROM: Read Only Memory), an erasable programmable read-only memory (EPROM: Erasable Programmable Read Only Memory or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM: CD-Read Only Memory), an optical storage device, a magnetic storage device, or any suitable combination of the above. In this embodiment, the computer-readable storage medium may be any tangible medium containing or storing a program, which may be used by or in combination with an instruction execution system or device. The program code contained on the computer-readable storage medium may be transmitted using any appropriate medium, including but not limited to: wires, optical cables, RF (Radio Frequency: Radio Frequency), etc., or any suitable combination of the above.

[0096] The computer-readable storage medium may be included in the geomagnetic field modeling device; or may exist independently without being assembled into the geomagnetic field modeling device.

[0097] The above-mentioned computer-readable storage medium carries one or more programs. When the above-mentioned one or more programs are executed by the geomagnetic field modeling device, the geomagnetic field modeling device: obtains the scale factor matrix of the spherical basic current system SECS; obtains the first geomagnetic effect conversion matrix of the interpolation point; determines the magnetic field matrix of the interpolation point according to the scale factor matrix and the first geomagnetic effect conversion matrix; and constructs a geomagnetic field model based on the magnetic field matrix of the interpolation point.

[0098] Computer program code for performing the operations of the present application may be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages ​​such as Java, Smalltalk, C++, and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code may be executed entirely on the user's computer, partially on the user's computer, as a separate software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0099] The flow chart and block diagram in the accompanying drawings illustrate the possible architecture, function and operation of the system, method and computer program product according to various embodiments of the present application. In this regard, each square box in the flow chart or block diagram can represent a module, a program segment or a part of a code, and the module, the program segment or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two square boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or flow chart, and the combination of the square boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0100] The modules involved in the embodiments described in this application may be implemented by software or hardware, wherein the name of the module does not constitute a limitation on the unit itself in some cases.

[0101] The readable storage medium provided in the embodiment of the present application is a computer-readable storage medium, which stores computer-readable program instructions (i.e., computer programs) for executing the above-mentioned geomagnetic field modeling method, and can solve the technical problem of how to improve the accuracy of global magnetic field data. Compared with the prior art, the beneficial effects of the computer-readable storage medium provided in the present application are the same as the beneficial effects of the geomagnetic field modeling method provided in the above-mentioned embodiment, and will not be repeated here.

[0102] An embodiment of the present application also provides a computer program product, including a computer program, which implements the steps of the above-mentioned geomagnetic field modeling method when executed by a processor.

[0103] The computer program product provided in the embodiment of the present application can solve the technical problem of how to improve the accuracy of global magnetic field data. Compared with the prior art, the beneficial effects of the computer program product provided in the present application are the same as the beneficial effects of the geomagnetic field modeling method provided in the above embodiment, which will not be repeated here.

[0104] The above descriptions are only some embodiments of the present application, and are not intended to limit the patent scope of the present application. All equivalent structural changes made using the contents of the present application specification and drawings under the technical concept of the present application, or direct / indirect applications in other related technical fields are included in the patent protection scope of the present application.

Claims

1. A geomagnetic field modeling method, characterized in that: The geomagnetic field modeling method comprises: Get the scale factor matrix of the spherical basic current system SECS; Obtaining the first geomagnetic effect conversion matrix of the interpolation point; Determining the magnetic field matrix of the interpolation point according to the scale factor matrix and the first geomagnetic effect conversion matrix; A geomagnetic field model is constructed based on the magnetic field matrix of the interpolation points.

2. The geomagnetic field modeling method according to claim 1, characterized in that: The step of obtaining the first geomagnetic effect conversion matrix of the interpolation point comprises: The radial component and the zenith component of the magnetic field at the interpolation point are acquired, and the first geomagnetic effect conversion matrix of the interpolation point is determined according to the radial component and the zenith component of the magnetic field at the interpolation point.

3. The geomagnetic field modeling method according to claim 2, characterized in that: The step of obtaining the radial component and the zenith component of the magnetic field of the interpolation point comprises: Acquire the angle information between the interpolation point and the pole of the SECS and the radial distance of the interpolation point; The radial component and the zenith component of the magnetic field of the interpolation point are determined according to the angle information and the radial distance.

4. The geomagnetic field modeling method according to claim 2, characterized in that: The step of obtaining the radial component and the zenith component of the magnetic field of the interpolation point comprises: Obtaining a non-divergent current density vector of the SECS and a radial distance between the interpolation points; The radial component and the zenith component of the magnetic field at the interpolation point are determined according to the unkempt current density vector and the radial distance.

5. The geomagnetic field modeling method according to claim 1, characterized in that: The step of obtaining the scale factor matrix of the spherical basic current system SECS comprises: Obtain the second geomagnetic effect conversion matrix of the observation point; Obtaining a magnetic field disturbance vector of the observation point; A scale factor matrix of the SECS is determined according to the second geomagnetic effect conversion matrix and the magnetic field disturbance vector.

6. The geomagnetic field modeling method according to claim 5, characterized in that: The step of obtaining the second geomagnetic effect conversion matrix of the observation point comprises: The radial component and the zenith component of the magnetic field at the observation point are obtained, and the second geomagnetic effect conversion matrix of the observation point is determined according to the radial component and the zenith component of the observation point.

7. The geomagnetic field modeling method according to claim 1, characterized in that: Before the step of obtaining the scale factor matrix of the spherical basic current system SECS, the method further includes: Determine the grid area based on the distribution information of the observation points; On the grid area, the poles of the SECS are determined.

8. A geomagnetic field modeling device, characterized in that: The device comprises: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the computer program is configured to implement the steps of the geomagnetic field modeling method according to any one of claims 1 to 7.

9. A storage medium, characterized in that: The storage medium is a computer-readable storage medium, and a computer program is stored on the storage medium. When the computer program is executed by a processor, the steps of the geomagnetic field modeling method according to any one of claims 1 to 7 are implemented.

10. A computer program product, characterized in that The computer program product comprises a computer program, and when the computer program is executed by a processor, the steps of the geomagnetic field modeling method according to any one of claims 1 to 7 are implemented.

Images