A quality assessment method for aeromagnetic diurnal variation data

Abstract

The present invention belongs to the field of aeromagnetic survey technology, and specifically relates to a method for evaluating the quality of aeromagnetic daily variation data, comprising the following steps: obtaining aeromagnetic daily variation data and time attribute information of aeromagnetic survey lines; processing aeromagnetic daily variation data: taking time as a reference, calculating the three-minute change rate of aeromagnetic daily variation data, and constructing a three-minute change rate curve; using nonlinear filtering, obtaining the noise jitter curve of aeromagnetic daily variation data measurement by calculation; and visually expressing the three-minute change rate curve and the jitter curve. The present invention provides a method for evaluating the quality of aeromagnetic daily variation data, by calculating the three-minute change rate and noise jitter curve of aeromagnetic daily variation data, performing quality evaluation and visual expression on the aeromagnetic daily variation data, by which the method can intuitively judge whether a re-flight is required when a magnetic storm occurs, and can more intuitively and richly express the quality information of aeromagnetic daily variation data, thereby improving the efficiency of aeromagnetic survey operations.

Description

Technical Field

[0001] The invention belongs to the technical field of aeromagnetic survey, and in particular relates to a quality assessment method for aeromagnetic daily variation data. Background Art

[0002] In the process of aeromagnetic survey, the measurement of diurnal magnetic variation is one of the important auxiliary measures for aeromagnetic survey. The data obtained is the basis for the correction of diurnal magnetic variation. The quality of the data is directly related to the accuracy of aeromagnetic data processing and indirectly affects the effect of aeromagnetic interpretation results. Therefore, there is a clear requirement in the Technical Specifications for Aeromagnetic Survey (DZ / T 0142-2010) that a diurnal variation observation station must be established when conducting aeromagnetic survey. The purpose of aeromagnetic diurnal variation observation is to obtain the data of the change of the geomagnetic field over time during the flight, so as to eliminate the influence of diurnal variation and short-period disturbance of the geomagnetic field, thereby improving the processing accuracy of aeromagnetic data.

[0003] In 2000-2002, Xue Dianjun and others developed the AirProbe software based on the "Aerial Geophysical Data Processing Software System under WINDOWS System" project of the National Land and Resources Survey, which included visual magnetic diurnal variation data editing software, which can perform numerical statistics on the noise of aeromagnetic diurnal variation data after the flight measurement. In 2012, with the support of the "Research on Aerogeophysical Methods and Data Processing and Interpretation System" project of the National Land and Resources Survey, the project team improved and upgraded the visual magnetic diurnal variation data editing software based on C# language and plug-in technology. In recent years, with the development of high-precision aeromagnetic survey work, the refinement of measurement tasks and the application of unattended working mode, the requirements for data quality have become more stringent.

[0004] At present, there are few studies on the quality evaluation methods of magnetic diurnal variation data at home and abroad. The quality of aeromagnetic diurnal variation data is mainly controlled by ground surveyors who interpret and mark the simulated paper rolls output in real time during magnetic diurnal variation observations. This method relies too much on manual judgment and lacks an intuitive way to express data quality. Summary of the invention

[0005] In view of the above problems, the technical solution of the present invention is proposed. The present invention provides a method for evaluating the quality of aeromagnetic diurnal variation data, which is used to solve the problem that there are few studies on the quality evaluation methods of aeromagnetic diurnal variation data at home and abroad, and the quality of aeromagnetic diurnal variation data is mainly controlled by ground surveyors to interpret and mark the simulated paper rolls output in real time during the magnetic diurnal variation observation, which is overly dependent on manual judgment and lacks an intuitive expression method for data quality.

[0006] The present invention is implemented by the following technical solutions:

[0007] A method for evaluating the quality of aeromagnetic diurnal variation data comprises the following steps:

[0008] S100, obtaining aeromagnetic diurnal variation data and aeromagnetic survey line time attribute information;

[0009] S200, processing of aeromagnetic diurnal variation data:

[0010] Based on time, the three-minute change rate of the aeromagnetic daily variation data is calculated and a three-minute change rate curve is constructed;

[0011] Calculate the characteristic value of the curvature of the aeromagnetic diurnal variation data curve, use the characteristic value to nonlinearly iterate and approximate the aeromagnetic diurnal variation data curve, fit the filtered abnormal curve, and make a difference between the aeromagnetic diurnal variation data and the filtered magnetic diurnal variation data to obtain the measurement noise jitter curve;

[0012] S300, visually expressing the three-minute rate of change curve and the measurement noise jitter curve according to the time attribute information of the aeromagnetic survey line.

[0013] As a further solution of the present invention, the aeromagnetic diurnal variation data is the aeromagnetic diurnal variation observation data outputted in real time by a diurnal magnetic variation observation station when flying on a survey line.

[0014] As a further solution of the present invention, the calculation method of the three-minute change rate of the magnetic daily variation data is:

[0015] Taking time as the reference and time M as the calculation point, the calculation is performed using a change rate calculation model, and the change rate calculation model is:

[0016] d M =M Max -M Min

[0017] M Max =[M i-1.5 , M i+1.5 ]

[0018] M Min =[M i-1.5 , M i+1.5 ]

[0019] Where: M i is the magnetic measurement value at time i, M Max 、M Min Indicates the maximum and minimum magnetic measurement values ​​of the current calculation point M within a 3-minute time range, d M It is the change of magnetic measurement amplitude within 3 minutes.

[0020] As a further solution of the present invention, the method for constructing the three-minute change rate curve is:

[0021] The horizontal coordinate of the coordinate system is time, and the amplitude change of the magnetic measurement value is d. MThe three-minute change rate curve of the magnetic diurnal variation data is constructed as the ordinate of the coordinate system.

[0022] As a further solution of the present invention, the method for calculating the measurement noise jitter curve comprises the following steps:

[0023] S210, given a filter window width WD, the magnetic diurnal variation data curve is taken with the filter window width WD as the interval, and the characteristic value S is calculated using the following formula:

[0024] S i =△T i -0.5·(△T i-1 +△T i+1 )

[0025] Where: S i represents the eigenvalue of the i-th point, △T i represents the magnetic measurement value of the i-th point;

[0026] S220, judge S i ·S i+1 <0, S i ·S i-1 <0 whether they are satisfied at the same time, if so, proceed to step S230, if not, proceed to step S240;

[0027] S230, calculating the trend abnormal value to obtain the filtered magnetic daily variation data curve. The calculation formula of the trend abnormal value is as follows:

[0028] TF i =△T i -0.5 SK i

[0029] Among them: TF i represents the trend outlier of the i-th point;

[0030] S240, moving the filter window forward with the filter window width WD as the step length, and proceeding to step S100;

[0031] S250, calculate the measurement noise jitter curve of the aeromagnetic daily variation data, the calculation formula is as follows:

[0032] D j =△T i -TF i

[0033] Where: D j Represents the measurement noise jitter of the i-th point.

[0034] As a further solution of the present invention, the method of visual expression comprises the following steps:

[0035] S310, obtaining a three-minute rate of change curve and a measurement noise jitter curve;

[0036] S320, processing the three-minute change rate curve, the processing method is:

[0037] Match the time information of the noise jitter curve with the time attribute information of the aeromagnetic survey line, fill the time overlap segment on the three-minute change rate curve with color, and mark the survey line number; mark the part with a three-minute change rate exceeding 5nT in red;

[0038] S330, outputting a measured noise jitter curve and a processed three-minute rate of change curve to form an information quality report.

[0039] A system for executing the method for evaluating the quality of aeromagnetic diurnal variation data, the system comprising a processing module, a storage module, and a computer program stored in the storage module and executable on the processing module, wherein:

[0040] A storage module, used for storing a three-minute rate of change curve and a measurement noise jitter curve;

[0041] The processing module is used to implement the steps of the above-mentioned aeromagnetic diurnal variation data quality assessment method when executing the program.

[0042] The technical solution provided by the present invention has the following beneficial effects: the present invention provides a method for quality assessment of aeromagnetic diurnal variation data, by calculating the three-minute change rate of the magnetic diurnal variation data and the measurement noise jitter curve, the quality of the magnetic diurnal variation data is assessed, and a visualization method can be used to intuitively determine whether a re-flight is required when a magnetic storm occurs, and the magnetic diurnal variation data information can be expressed more intuitively and richly, thereby improving the efficiency of aeromagnetic survey operations, saving flight costs in a certain sense, and solving practical problems in production. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] The accompanying drawings are used to provide a further understanding of the present invention and constitute a part of the specification. Together with the embodiments of the present invention, they are used to explain the present invention and do not constitute a limitation of the present invention. In the accompanying drawings:

[0044] Figure 1 The flowchart of the method for evaluating the quality of aeromagnetic diurnal variation data according to Embodiment 1 of the present invention is shown in FIG.

[0045] Figure 2 It is a schematic diagram before and after fitting of the magnetic diurnal variation curve in the method for evaluating the quality of aeromagnetic diurnal variation data in Example 1 of the present invention.

[0046] Figure 3 The three-minute rate of change curve and the measurement noise jitter diagram are visualized expressions of the aeromagnetic diurnal variation data quality assessment method of Example 1 of the present invention.

[0047] Figure 4 It is a flow chart of a method for calculating a measurement noise jitter curve in a method for evaluating aeromagnetic diurnal variation data quality according to Embodiment 1 of the present invention. DETAILED DESCRIPTION

[0048] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0049] It should be noted that various aspects of the embodiments within the scope of the appended claims are described below. It should be apparent that the aspects described herein may be embodied in a wide variety of forms, and any specific structure and / or function described herein is merely illustrative. Based on the present disclosure, it should be understood by those skilled in the art that an aspect described herein may be implemented independently of any other aspect, and two or more of these aspects may be combined in various ways. For example, any number of aspects described herein may be used to implement the device and / or practice the method. In addition, other structures and / or functionalities other than one or more of the aspects described herein may be used to implement this device and / or practice this method.

[0050] It should be understood that although the terms first, second, etc. may be used to describe various information in the embodiments of the present invention, the information should not be limited to these terms. These terms are only used to distinguish the same type of information from each other.

[0051] In the process of aeromagnetic survey, the measurement of diurnal magnetic variation is one of the important auxiliary measures for aeromagnetic survey. The data obtained is the basis for the correction of diurnal magnetic variation. The quality of the data is directly related to the accuracy of aeromagnetic data processing and indirectly affects the effect of aeromagnetic interpretation results. Therefore, there is a clear requirement in the Technical Specifications for Aeromagnetic Survey (DZ / T 0142-2010) that a diurnal variation observation station must be established when conducting aeromagnetic survey. The purpose of aeromagnetic diurnal variation observation is to obtain the data of the change of the geomagnetic field over time during the flight, so as to eliminate the influence of diurnal variation and short-period disturbance of the geomagnetic field, thereby improving the processing accuracy of aeromagnetic data.

[0052] At present, there are few studies on the quality evaluation methods of magnetic diurnal variation data at home and abroad. The quality of aeromagnetic diurnal variation data is mainly controlled by ground surveyors who interpret and mark the simulated paper rolls output in real time during magnetic diurnal variation observations. This method relies too much on manual judgment and lacks an intuitive way to express data quality.

[0053] Based on this, see Figure 1The embodiment of the present application provides a method for evaluating the quality of aeromagnetic diurnal variation data, which is applied to the quality evaluation and visual expression of aeromagnetic diurnal variation data, and includes the following steps:

[0054] S100, obtaining aeromagnetic diurnal variation data and aeromagnetic survey line time attribute information;

[0055] S200, processing of aeromagnetic diurnal variation data:

[0056] Based on time, the three-minute change rate of the aeromagnetic daily variation data is calculated and a three-minute change rate curve is constructed;

[0057] Calculate the characteristic value of the curvature of the aeromagnetic diurnal variation data curve, use the characteristic value to nonlinearly iterate and approximate the aeromagnetic diurnal variation data curve to obtain the filtered magnetic diurnal variation data, and make a difference between the aeromagnetic diurnal variation data and the filtered magnetic diurnal variation data to obtain the measurement noise jitter curve;

[0058] S300, visually expressing the three-minute rate of change curve and the measurement noise jitter curve according to the time attribute information of the aeromagnetic survey line.

[0059] In the process of aeromagnetic survey, the measurement of diurnal magnetic variation is one of the important auxiliary methods for aeromagnetic survey. The data obtained is the basis for the correction of diurnal magnetic variation. The quality of the data is directly related to the accuracy of aeromagnetic data processing and indirectly affects the effect of aeromagnetic interpretation results. When conducting aeromagnetic survey, a diurnal magnetic variation observation station must be established. The purpose of aeromagnetic diurnal variation observation is to obtain the data of the change of the geomagnetic field over time during the flight. In step S100 of the embodiment of the present invention, the aeromagnetic diurnal variation data curve is obtained by using the aeromagnetic diurnal variation observation data outputted instantly by the magnetometer of the magnetic diurnal variation observation station.

[0060] In step S200 of the present invention, the calculation method of the three-minute change rate of the daily magnetic variation data is:

[0061] Taking time as the reference and time M as the calculation point, the calculation is performed using a change rate calculation model, and the change rate calculation model is:

[0062] d M =M Max -M Min

[0063] M Max =[M i-1.5 , M i+1.5 ]

[0064] M Min =[M i-1.5 , M i+1.5 ]

[0065] Where: M iis the magnetic measurement value at time i, M Max 、M Min Indicates the maximum and minimum magnetic measurement values ​​of the current calculation point M within a 3-minute time range, d M It is the change of magnetic measurement amplitude within 3 minutes;

[0066] It should be noted that, in the embodiment of the present invention, time 0 is taken as the starting point, that is, i≥0, when 0≤i≤1.5, M i-1.5 Represents the magnetic measurement value at time 0, that is, M i-1.5 =M0, when i>1.5, M i-1.5 Indicates the magnetic measurement value at time (i-1.5);

[0067] Also, if Figure 4 As shown, in step S200, the method for constructing the three-minute change rate curve is:

[0068] With time T as the horizontal coordinate of the coordinate system, the amplitude change d of the magnetic measurement value M The three-minute rate of change curve of the daily magnetic variation data is constructed using the ordinate of the coordinate system, such as Figure 3 shown.

[0069] like Figure 4 As shown, in step S200 of the embodiment of the present invention, the method for calculating the measured noise jitter curve includes the following steps:

[0070] S210, given a filter window width WD, the magnetic diurnal variation data curve is taken with the filter window width WD as the interval, and the characteristic value S is calculated using the following formula:

[0071] S i =△T i -0.5·(△T i-1 +△T i+1 )

[0072] Where: S i represents the eigenvalue of the i-th point, △T i represents the magnetic measurement value of the i-th point;

[0073] S220, judge S i ·S i+1 <0, S i ·S i-1 <0 whether they are satisfied at the same time, if so, proceed to step S230, if not, proceed to step S240;

[0074] S230, calculating the trend abnormal value to obtain the filtered magnetic daily variation data curve. The calculation formula of the trend abnormal value is as follows:

[0075] TF i =△T i-0.5 SK i

[0076] Among them: TF i represents the trend outlier of the i-th point;

[0077] S240, moving the filter window forward with the filter window width WD as the step length, and proceeding to step S100;

[0078] S250, calculate the measurement noise jitter curve of the aeromagnetic daily variation data, the calculation formula is as follows:

[0079] D j =△T i -TF i

[0080] Where: D j Represents the measurement noise jitter of the i-th point.

[0081] It is understandable that in step S220, S i+1 and S i-1 represents a point one filter window width away from point i, then S i+1 and S i-1 When the S values ​​of point i, point i+1, and point i-1 are alternately positive and negative, which are two filter window widths apart, it means that the curves of these three points have the same bending direction and have the same regional anomaly. This serves as the basis for the existence of local anomalies, that is, there is noise interference in the area and filtering is required.

[0082] In addition, the filtering window width WD is generally the distance between the half-extreme values ​​of the abnormal curve, and the abnormal curve in the embodiment of the present invention is the aeromagnetic diurnal variation data curve.

[0083] Please refer again Figure 2 , as a specific application of the embodiment of the present invention, Figure 2 In the calculation, point K is used as the calculation center, and K with a window width WD away from point K is selected. M1 and K P1， Calculate points K, K M1 and K P1 The characteristic values ​​of are:

[0084] SK1=△T(K M1 )-0.5·(△T(K M2 )+△T(K))

[0085] SK2=△T(K)-0.5·(△T(K M1 )+△T(K P1 ))

[0086] SK3=△T(K P1)-0.5·(△T(K)+△T(K P2 ))

[0087] Assume SK l , SK2, SK3 simultaneously satisfy the conditions SK3·SK2<0, SK2·SK1<0, calculate K Ml , K, K Pl Three-point trend outliers:

[0088] TF(K M1 )=△T(K M1 )-0.5·SK1

[0089] TF(K)=△T(K)-0.5·SK2

[0090] TF(K P1 )=△T(K P1 )-0.5·SK2

[0091] After the calculation is completed, the filter window is moved forward with the filter window width WD as the step length, and the above steps are repeated. If SK l , SK2, and SK3 do not simultaneously satisfy the conditions SK3·SK2<0, SK2·SK1<0, no calculation is required, and the above steps are repeated until the calculation of the entire curve is completed. K Ml , K, K Pl The abnormal characteristic values ​​of the three points satisfy SK3·SK2<0, SK2·SK1<0 at the same time, indicating that the curve is in SK l The bending directions in SK2 and SK3 are the same, and they have the same regional anomaly, that is, there is noise interference in the area and filtering is required.

[0092] In the embodiment of the present invention, the method of visual expression in step S300 includes the following steps:

[0093] S310, obtaining a three-minute rate of change curve and a measurement noise jitter curve;

[0094] S320, processing the three-minute change rate curve, the processing method is:

[0095] Match the time information of the noise jitter curve with the time attribute information of the aeromagnetic survey line, fill the time overlap segment on the three-minute change rate curve with color, and mark the survey line number; mark the part with a three-minute change rate exceeding 5nT in red;

[0096] S330, outputting a measured noise jitter curve and a processed three-minute rate of change curve to form an information quality report.

[0097] Please combine Figure 3In the embodiment of the present invention, it should be noted that in order to make the attached Figure 3 The drawings required by the application documents must be met. Figure 3 When filling the time overlap segment with color, gray is used. In actual applications, other colors can be used to fill the time overlap segment to make it stand out. Figure 3 The time overlap period is 13:00-16:00, and the marked measurement line numbers include 11695, 11715, 11705, 11725, 11735 and 11745. Among them, the measurement line numbers with a three-minute change rate exceeding 5nT are 11695 and 11735. Since the curve graph of the application document can only be in black, the attached Figure 3 It is not possible to display it in red, but in actual application, the part of the curve where the rate of change exceeds 5nT in three minutes is displayed in red;

[0098] In addition, the three-minute change rate of the daily aeromagnetic data exceeding 5nT is often caused by a magnetic storm during the aeromagnetic survey, which results in the data quality not meeting the requirements of the aeromagnetic survey specifications. A magnetic storm usually lasts for more than ten hours around the world. As long as it has no effect on the final aeromagnetic survey results, the occurrence of a magnetic storm can be ignored. Figure 3 It can be seen that the 3-minute rate of change curve exceeds 5nT at 4 points, which is consistent with the data characteristics of a magnetic storm. Considering that a magnetic storm occurred during the data measurement of this flight, some areas of the survey line were filled.

[0099] In addition, during the aeromagnetic survey, the entire flight includes route flight and line flight. Route flight is a non-effective flight outside the general survey area, and line flight is used to refer to an effective flight carried out strictly in accordance with the planned survey line. The final measurement result is based on the line flight. Therefore, in the embodiment of the present invention, the aeromagnetic daily variation data is the aeromagnetic daily variation observation data output by the magnetometer in real time during the line flight. Please refer to Figure 3 , Figure 3 Only two of the four magnetic storms occurred on the survey lines, so only survey lines 11695 and 11735 needed to be re-surveyed;

[0100] The above application shows that when a magnetic storm occurs, it is not necessary to re-fly. You should first determine whether the magnetic storm causes the three-minute change rate of the aeromagnetic daily variation data to exceed 5nT, and whether the time period when the three-minute change rate exceeds the standard coincides with the survey line. Then, determine the location and frequency of the three-minute change rate exceeding the standard, and consider multiple factors such as the comprehensive flight cost to decide whether to re-fly.

[0101] The present invention also discloses a system for executing the method for evaluating the quality of aeromagnetic diurnal variation data, the system comprising a processing module, a storage module, and a computer program stored in the storage module and executable on the processing module, wherein:

[0102] A storage module, used for storing a three-minute rate of change curve and a measurement noise jitter curve;

[0103] The processing module is used to implement the steps of the above-mentioned aeromagnetic diurnal variation data quality assessment method when executing the program.

[0104] In addition, some embodiments may include a storage medium having a program for executing the methods described in this specification on a computer, on which at least one instruction, at least one program, code set or instruction set is stored, and when the at least one instruction, at least one program, code set or instruction set is loaded and executed by a processor, the steps in the above-mentioned method embodiments are implemented. Examples of the computer-readable recording medium include hardware devices specially constructed to store and execute program commands: magnetic media such as hard disks, floppy disks and tapes, optical recording media such as CD-ROMs and DVDs, magneto-optical media such as floppy disks, and ROMs, RAMs, flash memories, etc. Examples of program commands may include: machine language codes written by a compiler and high-level language codes executed by a computer using an interpreter, etc.

[0105] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through at least one instruction, at least one program, code set or instruction set, and the at least one instruction, at least one program, code set or instruction set can be stored in a non-volatile computer-readable storage medium. When the at least one instruction, at least one program, code set or instruction set is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory.

[0106] In summary, the present invention provides a method for quality assessment of aeromagnetic daily variation data. The quality of the aeromagnetic daily variation data is assessed by calculating the three-minute change rate of the aeromagnetic daily variation data and the measurement noise jitter curve. The visualization method can be used to intuitively determine whether a re-flight is required when a magnetic storm occurs. The aeromagnetic daily variation data information can be expressed more intuitively and richly, thereby improving the efficiency of aeromagnetic survey operations, saving flight costs in a certain sense, and solving practical problems in production.

[0107] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for evaluating the quality of aeromagnetic diurnal variation data, characterized in that: The following steps are involved: Obtain aeromagnetic diurnal variation data and time attribute information of aeromagnetic survey lines; Processing of aeromagnetic diurnal data: Based on time, the three-minute change rate of the aeromagnetic daily variation data is calculated and a three-minute change rate curve is constructed; Calculate the characteristic value of the curvature of the aeromagnetic diurnal variation data curve, use the characteristic value to nonlinearly iterate and approximate the aeromagnetic diurnal variation data curve to obtain the filtered magnetic diurnal variation data, and make a difference between the aeromagnetic diurnal variation data and the filtered magnetic diurnal variation data to obtain the measurement noise jitter curve; The three-minute rate of change curve and the measurement noise jitter curve are visualized according to the time attribute information of the aeromagnetic survey line; The calculation method of the three-minute change rate of the aeromagnetic daily variation data is as follows: Taking time as the reference and time M as the calculation point, the calculation is performed using a change rate calculation model, and the change rate calculation model is: d M =M Max -M Min M Max =[M i-1.5 ,M i+1.5 ] M Min =[M i-1.5 ,M i+1.5 ] Where: M i is the magnetic measurement value at time i, M Max 、M Min Indicates the maximum and minimum magnetic measurement values ​​of the current calculation point M within a 3-minute time range, d M It is the change of magnetic measurement amplitude within 3 minutes.

2. The method for evaluating the quality of aeromagnetic diurnal variation data according to claim 1, characterized in that: The aeromagnetic diurnal variation data are aeromagnetic diurnal variation observation data outputted in real time by a magnetic diurnal variation observation station when flying on a survey line.

3. The method for evaluating the quality of aeromagnetic diurnal variation data according to claim 1, characterized in that: The construction method of the three-minute change rate curve is: The horizontal coordinate of the coordinate system is time, and the amplitude change of the magnetic measurement value is d. M The three-minute change rate curve of the magnetic diurnal variation data is constructed as the ordinate of the coordinate system.

4. The method for evaluating the quality of aeromagnetic diurnal variation data according to claim 1, characterized in that: The method for calculating the measurement noise jitter curve comprises the following steps: S210, given a filter window width WD, the magnetic diurnal variation data curve is taken with the filter window width WD as the interval, and the characteristic value S is calculated using the following formula: S i =△T i -0.5·(△T i-1 +△T i+1 ) Where: S i represents the eigenvalue of the i-th point, △T i represents the magnetic measurement value of the i-th point; S220, judge S i ·S i+1 <0, S i ·S i-1 <0 whether they are satisfied at the same time, if so, proceed to step S230, if not, proceed to step S240; S230, calculating the trend abnormal value to obtain the filtered magnetic diurnal variation data curve. The calculation formula of the trend abnormal value is as follows: TF i =△T i -0.5 SK i Among them: TF i represents the trend outlier of the i-th point; S240, moving the filter window forward with the filter window width WD as the step length, and proceeding to step S210; S250, calculate the measurement noise jitter curve of the aeromagnetic daily variation data, the calculation formula is as follows: D j =△T i -TF i Where: D j Represents the measurement noise jitter of the i-th point.

5. The method for evaluating the quality of aeromagnetic diurnal variation data according to claim 1, characterized in that: The method of visual expression comprises the following steps: Obtain three-minute rate of change curve and measurement noise jitter curve; The three-minute change rate curve is processed as follows: Match the time information of the noise jitter curve with the time attribute information of the aeromagnetic survey line, fill the time overlap segment on the three-minute change rate curve with color, and mark the survey line number; mark the part with a three-minute change rate exceeding 5nT in red; Output the measured noise jitter curve and the processed three-minute rate of change curve to form an information quality report.

6. A system for executing the method for evaluating the quality of aeromagnetic diurnal variation data as claimed in claim 5, characterized in that: The system comprises a processing module, a storage module, and a computer program stored in the storage module and executable on the processing module, wherein: A storage module, used for storing a three-minute rate of change curve and a measurement noise jitter curve; A processing module is used to implement the steps of the method for evaluating the quality of aeromagnetic diurnal variation data as claimed in any one of claims 1 to 5 when executing the program.

Images