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

Abstract

The present application discloses a geoelectric field modeling method, device, medium and computer program product, which relates to the field of geoelectric field modeling technology, including: obtaining the regional electric field signal of the point to be determined; obtaining the local electric field signal of the point to be determined; determining the total electric field signal of the point to be determined based on the regional electric field signal and the local electric field signal; and constructing a geoelectric field model based on the total electric field signal. The present application calculates the total electric field signal by constructing a geoelectric field model, and divides the total electric field signal into a regional electric field signal not affected by a magnetic storm and a local electric field signal affected by a magnetic storm, thereby improving the accuracy of the total electric field signal calculated by the geoelectric field model.

Description

Technical Field

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

[0002] With the deepening of understanding of solar activity and its impact on the Earth's magnetic field, geoelectric field modeling technology has gradually become a research hotspot. With the advancement of technology, geoelectric field modeling technology has begun to involve more complex interactions between the earth's electrical structure and the geomagnetic field, as well as how to more accurately simulate the impact of geomagnetic storms on the geoelectric field.

[0003] At present, a variety of methods have been developed for geoelectric field modeling technology, including the plane wave method, the complex image method, and the fast Hankel function change method. These methods can achieve relatively accurate modeling in areas with flat terrain and relatively simple geoelectric structures. However, for areas with complex geoelectric structures, the application of these methods is limited because they need to ensure that there are no sudden changes in the geoelectric structure in the horizontal direction in the area. If this cannot be guaranteed, the accuracy of the electric field signal calculated using the geoelectric field model will be reduced. In addition, these methods also rely on the accuracy of the geoelectric conductivity model. Once there is an error in the model, it will directly affect the accuracy of the electric field signal obtained according to the model. Summary of the invention

[0004] The main purpose of the present application is to provide a geoelectric field modeling method, device, medium and computer program product, aiming to solve the technical problem of low accuracy of electric field signals determined based on geoelectric field models.

[0005] To achieve the above objectives, the present application proposes a geoelectric field modeling method, the method comprising:

[0006] Obtain the regional electric field signal of the point to be found;

[0007] Acquire the local electric field signal of the point to be determined;

[0008] Determining the total electric field signal of the point to be determined according to the regional electric field signal and the local electric field signal;

[0009] A geoelectric field model is constructed based on the total electric field signal.

[0010] In one embodiment, the step of obtaining the regional electric field signal of the point to be determined includes:

[0011] Determining the regional magnetic field signal of the point to be determined according to the magnetic field signal of the preset observation point and the first horizontal magnetic transfer function;

[0012] The regional electric field signal of the point to be determined is determined according to the regional magnetic field signal of the point to be determined and the local electromagnetic transfer function.

[0013] In one embodiment, the step of obtaining the regional electric field signal of the point to be determined includes:

[0014] According to the magnetic field signal of the preset observation point and the electromagnetic transfer function between stations, the regional electric field signal of the point to be determined is determined.

[0015] In one embodiment, the step of obtaining the local electric field signal of the point to be determined includes:

[0016] Obtain the local magnetic field signal of the preset geomagnetic station;

[0017] Determining the local magnetic field signal of the point to be determined according to the local magnetic field signal of the geomagnetic station;

[0018] The local electric field signal of the point to be determined is determined according to the local magnetic field signal of the point to be determined and the local electromagnetic transfer function.

[0019] In one embodiment, the step of obtaining a local magnetic field signal of a preset geomagnetic station includes:

[0020] Determine the regional magnetic field signal of the geomagnetic station according to the magnetic field signal of the preset observation point and the second horizontal magnetic transfer function;

[0021] The total magnetic field signal of the geomagnetic station is acquired, and the local magnetic field signal of the geomagnetic station is determined according to the regional magnetic field signal of the geomagnetic station and the total magnetic field signal.

[0022] In one embodiment, after the step of constructing a geoelectric field model based on the total electric field signal, the method further includes:

[0023] Based on the geoelectric field model, obtaining a calculated value of the total electric field signal of the point to be determined;

[0024] Determining an error between the calculated value and the obtained measured value of the total electric field signal of the point to be determined;

[0025] According to the error, the selection of the geomagnetic station is adjusted.

[0026] In one embodiment, after the step of constructing a geoelectric field model based on the total electric field signal, the method further includes:

[0027] Performing an inverse Fourier transform on the total electric field signal;

[0028] A visualization module based on the geoelectric field model is used to represent the transformed total electric field signal.

[0029] In addition, to achieve the above-mentioned purpose, the present application also proposes a geoelectric field modeling device, the geoelectric field modeling device comprising:

[0030] A regional signal acquisition module is used to obtain the regional electric field signal of the point to be found;

[0031] A local signal acquisition module, used to acquire the local electric field signal of the point to be determined;

[0032] An electric field signal determination module, used to determine the total electric field signal of the point to be determined based on the regional electric field signal and the local electric field signal;

[0033] The model building module is used to build a geoelectric field model based on the total electric field signal.

[0034] In addition, to achieve the above-mentioned purpose, the present application also proposes a geoelectric field modeling device, which includes: 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 geoelectric field modeling method as described above.

[0035] 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 geoelectric field modeling method described above are implemented.

[0036] 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, the steps of the geoelectric field modeling method described above are implemented.

[0037] One or more technical solutions proposed in the present application have at least the following technical effects: separately obtaining the regional electric field signal and the local electric field signal of the point to be determined, which helps to reduce the error or deviation that may be caused by a single signal source and improve the accuracy of the electric field signal acquisition; then, based on the regional electric field signal and the local electric field signal, the total electric field signal of the point to be determined is determined, and the accuracy of the determined electric field signal can be further improved by comprehensively considering the local electric field signal affected by the magnetic storm and the regional electric field signal not affected by the magnetic storm; based on the total electric field signal, a geoelectric field model is constructed, and the electric field signals of various locations around the world can be calculated, thereby providing a more accurate prediction of the electric field distribution. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] 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.

[0039] 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.

[0040] Figure 1 A schematic diagram of a flow chart provided for the first embodiment of the electric field modeling method of the present application;

[0041] Figure 2 A schematic diagram of the electric field signal changing with time provided in the second embodiment of the present application;

[0042] Figure 3 A schematic diagram of the module structure of the geoelectric field modeling program provided in Example 2 of the present application;

[0043] Figure 4 A schematic diagram of a brief flow chart of the geoelectric field modeling method provided in the second embodiment of the present application;

[0044] Figure 5 This is a schematic diagram of the module structure of the electric field modeling device according to an embodiment of the present application;

[0045] Figure 6 Schematic diagram of the device structure of the hardware operating environment involved in the method for modeling the ground electric field in the embodiment of the present application. DETAILED DESCRIPTION

[0046] 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.

[0047] 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.

[0048] In this embodiment, for the convenience of description, the following description is made with the terminal as the execution subject.

[0049] In conventional technology, the plane wave method may not be able to accurately reflect the changing behavior of the electromagnetic field at the lateral mutation of conductivity when dealing with complex geoelectric structures, especially in the modeling of complex geoelectric structures where analytical expressions cannot be derived; the complex mirror method becomes complicated and difficult to achieve fast calculation when dealing with geoelectric structures with multiple interfaces or irregular shapes. This method relies on the concept of mirror points, which may be difficult to find suitable mirror points in actual geological structures, limiting its application under complex geological conditions; although the fast Hankel change method has relatively high computational efficiency, it often has problems with unstable computational accuracy when dealing with high-frequency signals. These methods are difficult to ensure the accuracy and reliability of the modeling results.

[0050] The present application provides a solution, by obtaining the magnetic field information of a preset observation point that is not affected by a magnetic storm, and then using the electromagnetic transfer function and the horizontal magnetic transfer function to determine the regional electric field signal of the point to be determined. This step can effectively utilize limited observation data, and through frequency domain analysis and transfer function calculation, quickly derive the response of the geoelectric field and the geomagnetic field at different frequencies, thereby improving the calculation efficiency; at the same time, by obtaining the local magnetic field signal of the observation site, the local magnetic field signal of the point to be determined is determined. This step mainly uses the magnetic field signal of the known point to determine the magnetic field signal of the unknown point; then using the electromagnetic transfer function, the local magnetic field signal is converted into a local electric field signal, and then the total electric field signal of the point to be determined can be determined according to the regional electric field signal and the local electric field signal. Comprehensively considering the local electric field signal affected by the magnetic storm and the regional electric field signal not affected by the magnetic storm can further reduce the error and improve the calculation accuracy; based on the total electric field signal, a geoelectric field model is constructed, and the electric field signals of various locations around the world can be calculated, thereby providing a more accurate prediction of the electric field distribution; for the geoelectric field model, the relevant parameters can also be adjusted according to the error between the calculated value and the measured value, thereby further improving the calculation accuracy of the geoelectric field model.

[0051] It should be noted that the execution subject of this embodiment may 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 geoelectric field modeling device as an example of the execution subject to illustrate this embodiment and the following embodiments.

[0052] Based on this, the present application embodiment provides a geoelectric field modeling method, referring to Figure 1 , Figure 1 FIG. 1 is a flow chart of a first embodiment of an electric field modeling method of the present application.

[0053] In this embodiment, the geoelectric field modeling method includes steps S10 to S40:

[0054] Step S10, obtaining the regional electric field signal of the point to be determined;

[0055] It should be noted that the regional electric field signal refers to the electric field data calculated at a specified geographical location (the point to be determined), which is not affected by magnetic storms or is less affected by magnetic storms, including the intensity, direction and other relevant data of the electric field.

[0056] Exemplarily, the regional electric field signal at the location point is calculated by acquiring magnetic field data and inter-station electromagnetic transfer functions of magnetic field observation sites that are not affected by or are less affected by magnetic storms. The inter-station electromagnetic transfer function can be determined by historical electric field data and historical magnetic field data between the two points.

[0057] Step S20, obtaining the local electric field signal of the point to be determined;

[0058] It should be noted that the local electric field signal refers to the electric field disturbance data after the specified geographical location (the point to be determined) is affected by the magnetic storm, which also includes information such as the intensity and direction of the electric field.

[0059] Exemplarily, a geomagnetic station that is close to the point to be determined and has accurate observation data is selected, and the electric field observation data during the magnetic storm is collected from the selected geomagnetic station, including the geomagnetic horizontal component, vertical component, etc. The geographical location difference between the point to be determined and the geomagnetic station is taken into account to simulate and predict the impact of the magnetic storm on the electric field of the point to be determined, thereby determining the local electric field signal of the point to be determined.

[0060] It can be understood that by dividing the magnetic source into regional and local, the calculation accuracy of the electric field signal can be effectively improved by respectively obtaining the regional electric field signals that are not affected by the magnetic storm or are less affected by the magnetic storm and the local electric field signals that are affected by the magnetic storm at a specified geographical location.

[0061] Step S30, determining the total electric field signal of the point to be determined according to the regional electric field signal and the local electric field signal;

[0062] It should be noted that the total electric field signal refers to the electric field portion of the electromagnetic field generated by charge distribution or current flow, usually expressed in the form of electric field strength, which can be captured by sensors or measuring devices and converted into processable digital signals.

[0063] Exemplarily, the total electric field signal can be divided into a regional electric field signal that is not affected by the magnetic storm or is less affected by the magnetic storm and a local electric field signal that is affected by the magnetic storm, wherein the regional electric field signal can be converted from the regional magnetic field signal using the electromagnetic transfer function; similarly, the local electric field signal can be converted from the local magnetic field signal using the electromagnetic transfer function. The specific conversion process can refer to the following electric field signal determination formula:

[0064]

[0065] In this formula, represents the total electric field signal in the frequency domain at point A (the point to be determined), represents the regional electric field signal in the frequency domain at point A, represents the local electric field signal in the frequency domain at point A, represents the local electromagnetic transfer function at point A, represents the regional magnetic field signal in the frequency domain at point A, represents the local magnetic field signal in the frequency domain at point A. Among them, the local electromagnetic transfer function The acquisition can be calculated by considering the physical laws between the magnetic field and the electric field (such as Maxwell's equations) and the medium characteristics (such as electrical conductivity, magnetic permeability, etc.), or it can be inferred from the historical electric field data and magnetic field data of point A. This embodiment does not make specific limitations on this.

[0066] It can be understood that by dividing the electric field signal into regional electric field signals that are not affected by magnetic storms or are less affected by magnetic storms and local electric field signals that are affected by magnetic storms, the regional electric field signals and the local electric field signals are calculated separately, thereby improving the accuracy and reliability of the calculated electric field signals.

[0067] Step S40: constructing a geoelectric field model based on the total electric field signal.

[0068] It should be noted that the geoelectric field model is a mathematical or physical model used to describe and simulate the electric field distribution at a certain location on the earth's surface or underground. The model is usually composed of parameters such as electric field strength, direction, and distribution law, and can be presented in the form of mathematical formulas, charts, or three-dimensional visualization.

[0069] Exemplarily, according to the method of acquiring the total electric field signal, a geoelectric field model is constructed. The total electric field signal can be obtained by adding the regional electric field signal and the local electric field signal. Therefore, the geoelectric field model at least includes the addition equation, and the parameters that need to be acquired for the geoelectric field model include the regional electric field signal and the local electric field signal. In addition, the regional electric field signal can be converted according to the regional magnetic field signal. Therefore, the geoelectric field model at least includes the addition equation and the magnetic field conversion equation. The parameters that need to be acquired for the geoelectric field model include the regional magnetic field signal and the local electric field signal. Similarly, the method of acquiring these two parameters can also be further refined.

[0070] This embodiment provides a geoelectric field modeling method. By constructing a geoelectric field model, the total electric field signal of a specified geographical location can be calculated in real time and accurately, thereby improving calculation efficiency. At the same time, by dividing the total electric field signal into a regional electric field signal and a local electric field signal, the calculation result of the total electric field signal is also more accurate.

[0071] In a feasible implementation, step S10 may include:

[0072] Step S11, determining a regional magnetic field signal of a point to be determined according to a magnetic field signal of a preset observation point and a first horizontal magnetic transfer function;

[0073] It should be noted that the preset observation point refers to an observation site near the point to be determined that is not affected by the local magnetic storm. The magnetic field signal can be captured and recorded by sensors and other measuring equipment. After obtaining the magnetic field data of the site, operations such as screening, unit conversion, removal of linear trends and frequency screening can be performed for subsequent calculations.

[0074] In addition, it should be noted that the horizontal magnetic transfer function is a mathematical tool or model that converts magnetic field signals from one position to another, so that the magnetic field distribution of an unknown position (the point to be determined) can be inferred based on the magnetic field data of a known observation point. The first horizontal magnetic transfer function refers to the magnetic field data conversion function between the point to be determined and the observation point.

[0075] For example, a nearby observation site n that is not affected by the local magnetic storm is selected, and the specific regional magnetic field calculation method can refer to the following formula:

[0076]

[0077] In this formula, represents the change of magnetic field in the time domain at point A (the point to be determined), represents the horizontal magnetic transfer function from point n to point A, represents the magnetic field variation in the time domain at point n, Represents Fourier transform, which is mainly responsible for converting magnetic field data in the time domain into magnetic field data in the frequency domain. Represents the inverse Fourier transform, which is used to convert the magnetic field data in the frequency domain back to the time domain. Since the data provided by the observation station is in the time domain, and the tensor value of the horizontal magnetic transfer function corresponds one-to-one to the frequency, the magnetic field data conversion between different locations can only be performed in the frequency domain. Therefore, the magnetic field data provided by the observation station needs to be Fourier transformed before it can be used for the magnetic field data calculation of point A. At this time, the magnetic field signal in the frequency domain at point A has been obtained, but in order to facilitate drawing, the inverse Fourier transform can be used to convert the magnetic field signal in the frequency domain back to the time domain.

[0078] Step S12, determining the regional electric field signal of the point to be determined based on the regional magnetic field signal of the point to be determined and the local electromagnetic transfer function.

[0079] It should be noted that the electromagnetic transfer function refers to a function that describes the conversion relationship between the magnetic field and the electric field, and is used to calculate the electric field signal based on the known magnetic field signal. The local electromagnetic transfer function is responsible for converting the local magnetic field signal into the local electric field signal.

[0080] In this embodiment, the magnetic field signal of the point to be determined is determined by the horizontal magnetic transfer function and the magnetic field data of the observation point. Compared with other methods based on historical data and inference based on geographical characteristics, the regional magnetic field signal calculated by this method is more accurate and reliable. Then, the local electromagnetic transfer function is used to convert the regional magnetic field signal of the point to be determined into the regional electric field signal of the point to be determined, thereby improving the calculation accuracy of the regional electric field signal.

[0081] In a feasible implementation manner, step S10 includes:

[0082] Step S13, determining the regional electric field signal of the point to be determined based on the magnetic field signal of the preset observation point and the electromagnetic transfer function between stations.

[0083] It should be noted that the observation point refers to a nearby observation station near the point to be determined that is not affected by the local magnetic storm; the inter-station electromagnetic transfer function is responsible for converting the magnetic field signal at a certain point into the electric field signal at another point.

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

[0085]

[0086] In this formula, represents the electric field in the frequency domain at point A, represents the electromagnetic transfer function from point B to point A, represents the magnetic field in the frequency domain at point B.

[0087] In this implementation, by directly converting the magnetic field signal of the observation point into the regional electric field signal of the point to be calculated, the error of signal conversion can be reduced to a certain extent, thereby improving the calculation accuracy of the regional electric field signal.

[0088] In a feasible implementation manner, step S20 includes:

[0089] Step S21, obtaining a local magnetic field signal of a preset geomagnetic station;

[0090] It should be noted that geomagnetic stations refer to observation stations equipped with high-precision geomagnetic sensors, which are specially used to measure the geomagnetic field and can provide continuous and long-term geomagnetic field data. The initial geomagnetic stations can be preset according to the distance between the point to be determined and the geomagnetic station, the density of the geomagnetic station, etc., and then the selection of geomagnetic stations can be adjusted according to the relevant errors during the verification of the model accuracy.

[0091] Exemplarily, the local magnetic field signal can be determined by subtracting the regional magnetic field signal from the total magnetic field signal measured by the geomagnetic station.

[0092] Step S22, determining the local magnetic field signal of the point to be determined based on the local magnetic field signal of the geomagnetic station.

[0093] For example, since the spherical current basic system interpolation algorithm can be used to calculate the magnetic field signal of the position point according to the magnetic field signal of the known point, the local magnetic field signal of the point to be determined can be determined according to the local magnetic field signal of the geomagnetic station with the help of the spherical current basic system interpolation algorithm, and the following calculation formula can be referred to:

[0094]

[0095] In this formula, represents the local magnetic field signal in the time domain at point A, SECS represents the spherical current basic system interpolation algorithm, It represents the local magnetic field signal of the nearby geomagnetic stations used for interpolation calculation. The number of geomagnetic stations can be one or more. Within a certain spatial range, multiple geomagnetic stations adjacent to the point to be interpolated are selected as reference stations to improve the accuracy of the calculated local magnetic field signal of point A.

[0096] Exemplarily, a global magnetohydrodynamic model can be used to determine the local magnetic field signal of the point to be determined. When using the global magnetohydrodynamic model to calculate the local magnetic field signal of the point to be determined, in addition to the parameters required to process the local magnetic field signal of the geomagnetic station, some space weather parameters (such as solar wind parameters, etc.) need to be provided. By simulating the interaction process between the solar wind and the magnetosphere, the local magnetic field signal of the point to be determined can be calculated.

[0097] Step S23, determining the local electric field signal of the point to be determined according to the local magnetic field signal of the point to be determined and the local electromagnetic transfer function.

[0098] In this embodiment, the local magnetic field signal of the point to be determined is determined based on the local magnetic field signal of the geomagnetic station that is disturbed by the magnetic storm, which improves the accuracy of the calculation of the local magnetic field signal, thereby improving the accuracy of the local electric field signal and further improving the accuracy of the overall electric field signal.

[0099] In a feasible implementation, step S21 includes:

[0100] Step S211, determining a regional magnetic field signal of a geomagnetic station according to a magnetic field signal of a preset observation point and a second horizontal magnetic transfer function;

[0101] It should be noted that the preset observation point refers to a site near the geomagnetic observatory that is not affected by the magnetic storm. Its magnetic field signal can be captured by magnetic field measurement equipment and is usually expressed in parameters such as magnetic field intensity and direction. The second horizontal magnetic transfer function is responsible for realizing the conversion of magnetic field signals between the observation point and the geomagnetic observatory.

[0102] When the horizontal magnetic function between the point to be determined and a certain observation point is known, the magnetic field signal of the observation point is obtained, the horizontal magnetic transfer function is applied for calculation, and the regional magnetic field signal of the point to be determined is output.

[0103] Step S212, acquiring the total magnetic field signal of the geomagnetic station, and determining the local magnetic field signal of the geomagnetic station based on the regional magnetic field signal and the total magnetic field signal of the geomagnetic station.

[0104] It should be noted that the high-precision geomagnetic sensors deployed at the geomagnetic stations can measure the total magnetic field data of the location of the geomagnetic stations, including the total magnetic field strength and direction.

[0105] For example, the following formula may be used for calculation:

[0106]

[0107] In this formula, represents the local magnetic field signal in the time domain at site k, represents the total magnetic field signal in the time domain at site k, represents the regional magnetic field signal in the time domain at site k.

[0108] In this embodiment, by accurately calculating the local magnetic field signal of the geomagnetic station affected by the magnetic storm, the accuracy of the local magnetic field signal of the point to be determined based on the local magnetic field signal of the geomagnetic station can be improved, thereby improving the accuracy of the total electric field signal.

[0109] 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 later. On this basis, after step S40, it also includes:

[0110] Step A10, obtaining a calculated value of a total electric field signal of a point to be determined based on a geoelectric field model;

[0111] Step A20, determining the error between the calculated value and the measured value of the total electric field signal of the point to be determined;

[0112] It should be noted that the measured value refers to the physical quantity directly observed by the measuring equipment; the error refers to the difference between the calculated value and the actual measured value, which reflects the accuracy of the mathematical model or algorithm in simulating or predicting actual phenomena, and is an important basis for optimizing model parameters and adjusting model structure.

[0113] Step A30, adjusting the selection of geomagnetic stations according to the error.

[0114] It should be noted that a certain degree of error is normal, but if the error is relatively large, the error can be reduced and the calculation accuracy can be improved by adjusting the selection of geomagnetic stations.

[0115] For example, if it is found that the difference between the calculated value and the measured value is large, the observation data of some geomagnetic stations that are too far away can be abandoned, and only the geomagnetic stations near the point to be determined can be selected for determining the relevant electric field signals and magnetic field signals.

[0116] In this embodiment, by measuring the difference between the measured value and the calculated value of the electric field signal, the accuracy of the geoelectric field model can be verified; at the same time, according to the difference, the relevant parameters of the geoelectric field model are adjusted to reduce the error between the calculated value and the measured value, thereby improving the accuracy and reliability of the geoelectric field model.

[0117] In a feasible implementation manner, after step S40, the method further includes:

[0118] Step E10, performing inverse Fourier transform on the total electric field signal;

[0119] It should be noted that the inverse Fourier transform is a method of converting frequency domain signals back to time domain signals. By converting frequency domain signals into time domain signals, the change of electric field signals over time can be intuitively observed and plotted.

[0120] Step E20, characterizing the transformed total electric field signal based on the visualization module of the geoelectric field model.

[0121] It should be noted that the visualization module is mainly responsible for displaying information such as electric field signals and electric field distribution in the form of images, so as to intuitively display the transformation of electric field signals, as well as the difference information between the calculated value and the measurement of the electric field signal or other information that is easy to analyze.

[0122] For example, please refer to Figure 2 , Figure 2 A schematic diagram of the calculated and measured values ​​of the electric field signal of an ESK site changing with time in two components in different directions is provided, wherein the horizontal axis represents time in minutes, and the vertical axis represents the intensity information of the electric field signal in millivolts per kilometer. Figure 2 A in the figure shows the trend of the electric field intensity of the total electric field signal in the Ex direction over time. Figure 2 Figure B shows the trend of the electric field strength of the total electric field signal in the Ey direction over time. The black line represents the measured value of the total electric field signal strength in the Ex and Ey directions, and the red line represents the calculated value of the total electric field signal strength in the Ex and Ey directions based on the geoelectric field model. Figure 2It can be found that the measured and calculated values ​​of the total electric field signal intensity in the Ex and Ey directions have a high degree of coincidence, which shows that the calculation results of the geoelectric field model are accurate.

[0123] In this implementation, the regional electric field signal and the local electric field signal in the frequency domain are inversely transformed by Fourier transform and converted from the frequency domain back to the time domain to facilitate drawing. Then, through the visualization module of the geoelectric field model, the change of the electric field signal over time is presented in the form of an image, which is convenient for researchers to analyze and study.

[0124] 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.

[0125] The geoelectric field modeling program provided by the embodiment of the present application may include:

[0126] The calculation module is used to schedule each submodule to realize electric field simulation, magnetic field calculation and error calculation;

[0127] Grid module, used for interpolation calculation;

[0128] A visualization module, used to provide visualization information for the calculation module;

[0129] The data management module is used to save data information and provide data information to the computing module and the network module.

[0130] It should be noted that visualization information refers to information related to electric field signals in the format of visually presented pictures; data information refers to the basic information required to realize electric field simulation, magnetic field calculation and error calculation, and the calculated result data for visualization display, including the data required to complete the interpolation calculation and the calculated magnetic field data, electric field data, etc.

[0131] For example, please refer to Figure 3 , Figure 3A schematic diagram of the module structure of a geoelectric 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. The efficient simulation and analysis of the geoelectric field is realized through modular design. The calculation module 101 is the core of the program, which is mainly responsible for scheduling other submodules, performing magnetic field calculation and error calculation, and integrating various data information, such as visualization image parameter information and electric field signal calculation results, so as to realize electric field simulation. The network module 102 is responsible for interpolation calculation. When performing interpolation calculation, the geographic information and public parameters in the data management module 104 are read. The geographic information includes the coordinate information of the point to be calculated and the observation point, etc. The public parameters include the list of geomagnetic stations, the ionosphere radius and the earth radius, etc., and the calculated related magnetic field data is sent to the data management module for storage, or sent to the calculation module 101, so that it determines the electric field signal according to the magnetic field data. The visualization module 103 mainly provides visualization information for the calculation module 101, including data visualization and geographic location visualization. Data visualization refers to displaying the trajectory of the acquired electric field signal changing with a certain factor in a graphical manner; geographic location visualization refers to displaying the electric field signals of all calculated 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, common parameters, and calculation results, and retrieves and provides corresponding data according to the needs of the calculation module and the grid module.

[0132] In this embodiment, by constructing a geoelectric field modeling program, information such as the distribution and intensity of the electric field signal can be calculated more accurately, and the electric field data can be displayed in a visual form, thereby completing the simulation of the changes in the geoelectric field.

[0133] For example, to help understand the implementation process of the geoelectric field modeling method obtained by combining this embodiment with the above-mentioned embodiment 1, please refer to Figure 4 , Figure 4 A brief flowchart of a geoelectric field modeling method is provided, specifically:

[0134] Step S201, read and process magnetic field data, filter the magnetic field according to geomagnetic stations and time, perform unit conversion, remove linear trends and frequency filtering, etc.

[0135] Step S202, reading and processing the electromagnetic transfer function and the horizontal magnetic transfer function, and reading the relevant information of the electromagnetic transfer function and the horizontal magnetic transfer function that may be used according to the location information of the point to be determined and the observation stations in its vicinity.

[0136] Step S203, determining the electric field signal based on the geoelectric field model. After acquiring relevant data, the geoelectric field model can calculate the electric field signal of the point to be determined.

[0137] Step S204, outputting visualization results. For the convenience of analysis and display, after calculating the relevant electric field signals, the geoelectric field model will output visualization results of the electric field signals, such as a line graph showing the changes of the electric field signals over time.

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

[0139] The present application also provides a geoelectric field modeling device, please refer to Figure 5 , the geoelectric field modeling device comprises:

[0140] The regional signal acquisition module 10 is used to acquire the regional electric field signal of the point to be determined;

[0141] A local signal acquisition module 20 is used to acquire a local electric field signal of a desired point;

[0142] The electric field signal determination module 30 is used to determine the total electric field signal of the point to be determined based on the regional electric field signal and the local electric field signal;

[0143] The model building module 40 is used to build a geoelectric field model based on the total electric field signal.

[0144] The geoelectric field modeling device provided in the embodiment of the present application adopts the geoelectric field modeling method in the above embodiment, which can solve the technical problem of low accuracy of the electric field signal determined based on the geoelectric field model. Compared with the prior art, the beneficial effects of the geoelectric field modeling device provided in the present application are the same as the beneficial effects of the geoelectric field modeling method provided in the above embodiment, and other technical features in the geoelectric field modeling device are the same as the features disclosed in the above embodiment method, which will not be repeated here.

[0145] An embodiment of the present application provides a geoelectric 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 geoelectric field modeling method in the above-mentioned embodiment one.

[0146] Reference below Figure 6, which shows a schematic diagram of the structure of a geoelectric field modeling device suitable for implementing the embodiments of the present application. The geoelectric field modeling device in the embodiments 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 geoelectric 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.

[0147] like Figure 6 As shown, the geoelectric 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 geoelectric 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 geoelectric field modeling device to communicate with other devices wirelessly or by wire to exchange data. Although the figure shows a geoelectric field modeling device with various systems, 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 have alternatively.

[0148] 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.

[0149] The geoelectric field modeling device provided in the embodiment of the present application adopts the geoelectric field modeling method in the above embodiment, which can solve the technical problem of low accuracy of the electric field signal determined based on the geoelectric field model. Compared with the prior art, the beneficial effects of the geoelectric field modeling device provided in the present application are the same as the beneficial effects of the geoelectric field modeling method provided in the above embodiment, and the other technical features in the geoelectric field modeling device are the same as the features disclosed in the method of the previous embodiment, which will not be repeated here.

[0150] 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.

[0151] 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.

[0152] 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 geoelectric field modeling method in the above-mentioned embodiment.

[0153] 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 that can 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.

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

[0155] The above-mentioned computer-readable storage medium carries one or more programs. When the above-mentioned one or more programs are executed by the geoelectric field modeling device, the geoelectric field modeling device: obtains the regional electric field signal of the point to be determined; obtains the local electric field signal of the point to be determined; determines the total electric field signal of the point to be determined based on the regional electric field signal and the local electric field signal; and constructs a geoelectric field model based on the total electric field signal.

[0156] 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).

[0157] 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.

[0158] 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.

[0159] 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 geoelectric field modeling method, and can solve the technical problem of low accuracy of electric field signals determined based on geoelectric field models. 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 geoelectric field modeling method provided in the above-mentioned embodiment, and will not be described in detail here.

[0160] An embodiment of the present application also provides a computer program product, including a computer program, which implements the steps of the geoelectric field modeling method as described above when executed by a processor.

[0161] The computer program product provided in the embodiment of the present application can solve the technical problem of low accuracy of electric field signals determined based on the geoelectric field model. 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 geoelectric field modeling method provided in the above embodiment, which will not be described in detail here.

[0162] 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 under the technical concept of the present application and the contents of the specification and drawings 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 geoelectric field modeling method, characterized in that: The geoelectric field modeling method comprises: Acquire a regional electric field signal of the point to be determined, wherein the regional electric field signal is electric field data of the point to be determined when it is not affected by the magnetic storm; Obtain the local magnetic field signal of the preset geomagnetic station; Determining the local magnetic field signal of the point to be determined according to the local magnetic field signal of the geomagnetic station; Determine the local electric field signal of the point to be determined according to the local magnetic field signal of the point to be determined and the local electromagnetic transfer function, wherein the local electric field signal is the electric field disturbance data of the point to be determined after being affected by the magnetic storm; Determining the total electric field signal of the point to be determined according to the regional electric field signal and the local electric field signal; constructing a geoelectric field model based on the total electric field signal; Wherein, the step of obtaining the regional electric field signal of the point to be determined comprises: Determine the regional magnetic field signal of the point to be determined based on the magnetic field signal of the preset observation point and the first horizontal magnetic transfer function; determine the regional electric field signal of the point to be determined based on the regional magnetic field signal of the point to be determined and the local electromagnetic transfer function; Alternatively, the regional electric field signal of the point to be determined is determined based on the magnetic field signal of a preset observation point and the inter-station electromagnetic transfer function.

2. The geoelectric field modeling method as claimed in claim 1, characterized in that: The step of obtaining the local magnetic field signal of the preset geomagnetic station includes: Determine the regional magnetic field signal of the geomagnetic station according to the magnetic field signal of the preset observation point and the second horizontal magnetic transfer function; The total magnetic field signal of the geomagnetic station is acquired, and the local magnetic field signal of the geomagnetic station is determined according to the regional magnetic field signal of the geomagnetic station and the total magnetic field signal.

3. The geoelectric field modeling method according to claim 1, characterized in that: After the step of constructing a geoelectric field model based on the total electric field signal, the method further includes: Based on the geoelectric field model, obtaining a calculated value of the total electric field signal of the point to be determined; Determining an error between the calculated value and the obtained measured value of the total electric field signal of the point to be determined; According to the error, the selection of the geomagnetic station is adjusted.

4. The geoelectric field modeling method according to claim 1, characterized in that: After the step of constructing a geoelectric field model based on the total electric field signal, the method further includes: Performing an inverse Fourier transform on the total electric field signal; A visualization module based on the geoelectric field model is used to represent the transformed total electric field signal.

5. A geoelectric 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 geoelectric field modeling method according to any one of claims 1 to 4.

6. 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 geoelectric field modeling method according to any one of claims 1 to 4 are implemented.

7. 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 geoelectric field modeling method according to any one of claims 1 to 4 are implemented.

Images