A method and system for verifying the accuracy of tidal force estimation

Abstract

This invention provides a method and system for verifying the accuracy of tidal force estimation. By calculating the tidal force at a given measuring station and the force exerted by the Earth on the station, the total change in gravitational acceleration at that station is calculated and matched with data collected by a gravimeter. Positioning is then performed using a gravitational positioning method. The accuracy of the tidal force calculation is evaluated based on the positioning accuracy, thus verifying the accuracy of the tidal force estimation. Due to error propagation, even small errors in the tidal force calculation will manifest as significant positioning errors in the positioning results. This invention, by verifying the tidal force through the positioning results, can greatly improve the verification accuracy.

Description

Technical Field

[0001] This invention relates to the field of measurement and positioning technology, and specifically to a technical solution for verifying the accuracy of tidal force estimation. Background Technology

[0002] Currently, there are many methods for verifying the accuracy of tidal forces, including tidal height observation, satellite altimetry, ocean dynamics models, and geophysical observation. However, these methods generally suffer from drawbacks such as numerous influencing factors, long time cycles, and high costs. Unlike conventional solutions, this invention proposes a method for verifying the accuracy of tidal force estimation using gravitational positioning (CN113310486A), considering that even small differences in tidal force can lead to significant differences in point coordinates. This method can effectively verify the accuracy of tidal force estimation and features short observation cycles, fewer influencing factors, and lower costs, making it a valuable method for verifying the accuracy of tidal force estimation. Summary of the Invention

[0003] To address the shortcomings of the existing technology, this invention proposes a technical solution for verifying the accuracy of tidal force.

[0004] To achieve the above objectives, the technical solution proposed in this invention is a method for verifying the accuracy of tidal force estimation. This method calculates the total gravitational acceleration change at a given measuring station by measuring the tidal force and the force exerted by the Earth on the station, and matches this calculation with data collected by a gravimeter. The station is then located using a gravitational positioning method. The accuracy of the tidal force calculation is evaluated based on the positioning accuracy, thus verifying the accuracy of the tidal force estimation. The implementation process includes the following steps:

[0005] (1) Calculate the tidal force exerted by the celestial body on the measuring station on Earth at time t;

[0006] (2) Calculate the resultant force of the gravitational force and centrifugal force exerted by the Earth on the measuring station on Earth at time t;

[0007] (3) Calculate the magnitude and direction of the total gravitational acceleration at the measuring stations on Earth;

[0008] (4) Conduct on-site gravity data collection at the measurement stations;

[0009] (5) Using the gravitational positioning method, the station is located by using the calculated gravitational acceleration and the gravitational acceleration collected in the field.

[0010] (6) Compare the two positioning results obtained in (5) to verify the positioning accuracy. Since a small tidal force error can cause a huge positioning error, the accuracy of the tidal force calculation can be verified based on the positioning error.

[0011] Furthermore, the method for calculating the tidal force of the celestial body on the measuring station on Earth in step (1) is as follows: the total tidal potential V of the celestial body on the measuring station on Earth is equal to V (1) and V (2) The sum of, where V (1) V represents the tidal range at a certain epoch of an celestial body at a certain measuring station on Earth. (2) This refers to the effect of Earth's oblateness on tidal levels.

[0012] Furthermore, the calculation method for the resultant force of the gravitational force and centrifugal force exerted by the Earth on the measuring station on Earth in step (2) is as follows: when considering the Earth's oblateness, the resultant force exerted by the Earth on the measuring station on Earth is calculated using the following formula:

[0013] F t =F g -F c

[0014] Among them, F t F represents universal gravitation. c This represents centrifugal force.

[0015] Furthermore, the formula for calculating the magnitude and direction of the total gravitational acceleration at each measuring station on Earth, taking into account tidal forces and the Earth's oblateness, is as follows:

[0016]

[0017] Where g is the total gravitational acceleration at the measuring station on Earth, and F t V represents universal gravitation, V is the total tidal potential exerted by the celestial body on the measuring station on Earth, and m is the mass of the measuring station.

[0018] Moreover, the gravitational positioning method in step (5) is implemented by combining the gravitational force generated by celestial bodies, the gravitational force generated by the Earth, and other inertial forces, which causes the gravitational acceleration of the location to be located to change. By periodically monitoring the change in gravitational acceleration of the location to be located, the positioning result of the location to be located is obtained by reverse calculation.

[0019] Moreover, the positioning accuracy in step (6) is expressed as planar accuracy, elevation accuracy and three-dimensional accuracy. Planar accuracy and elevation accuracy are the coordinate errors of a specified point in the horizontal plane and vertical direction, while three-dimensional accuracy is a combination of planar accuracy and elevation accuracy.

[0020] Furthermore, horizontal accuracy is assessed using circular probability error or root mean square error.

[0021] Furthermore, vertical accuracy assessment is expressed using elevation accuracy assessment, employing either absolute or relative accuracy.

[0022] On the other hand, the present invention also provides a tidal force estimation accuracy verification system, including a processor and a memory, wherein the memory is used to store program instructions, and the processor is used to call the stored instructions in the memory to execute a tidal force estimation accuracy verification method as described above.

[0023] On the other hand, the present invention also provides a tidal force estimation accuracy verification system, including a readable storage medium on which a computer program is stored, and when the computer program is executed, it implements a tidal force estimation accuracy verification method as described above.

[0024] This invention provides a method for verifying the accuracy of tidal force calculations. The calculated tidal force is matched with data collected by a gravimeter, and the location is determined using a gravitational positioning method. The accuracy of the tidal force calculation is evaluated based on the positioning accuracy. The gravitational positioning method is based on the combined effects of Earth's gravity, celestial gravity, and other inertial forces, which cause changes in the magnitude and direction of gravitational acceleration at the location to be located. The location is determined by periodically monitoring these changes in magnitude and direction. The positioning result is then compared with the actual coordinates to verify the accuracy of the tidal force calculation. Due to error propagation, even small errors in the tidal force calculation will manifest as significant positioning errors in the positioning results. Therefore, verifying the tidal force using the positioning results greatly improves the verification accuracy.

[0025] Compared with existing technologies:

[0026] This invention utilizes the principle that even small errors in tidal force can lead to significant deviations in the magnitude and direction of gravitational acceleration, resulting in substantial inaccuracies in positioning results. It employs a gravitational method for positioning, using two data sources for calculation, and verifies the accuracy of the tidal force estimation based on the differences in the positioning results. This method is characterized by its short timeframe, low cost, wide applicability, and fewer influencing factors, making it significantly different from traditional methods for verifying the accuracy of tidal force estimation.

[0027] The solution of this invention is simple to implement and highly practical. It solves the problems of high cost and inconvenience in practical application of existing technologies, improves user experience and enriches scientific research methods, and has important value. Attached Figure Description

[0028] Figure 1 This is a flowchart of a method according to an embodiment of the present invention. Detailed Implementation

[0029] The technical solution of the present invention will be described in detail below with reference to the accompanying drawings and embodiments.

[0030] This invention proposes that tidal force is a force with a relatively clear law of action, which can be estimated by physical laws such as the law of universal gravitation. This makes the study of tidal force have certain scientific significance and application value.

[0031] This invention utilizes the minute calculation deviation of the tidal force exerted by celestial bodies on the Earth. Although this deviation is small, it can still affect the Earth's gravitational field, causing changes in the magnitude and direction of gravitational acceleration. Since changes in gravitational acceleration can lead to significant deviations in positioning results, this invention leverages this tidal force error, employing a gravitational method for positioning. This invention uses two data sources for estimation and comparison to verify the accuracy of the tidal force estimation.

[0032] Compared to traditional methods for verifying the accuracy of tidal forces, this invention offers significant advantages. Traditional methods require a large amount of observational data, resulting in high verification costs, while this invention requires a shorter timeframe, is less expensive, and is applicable to a wider range of scenarios. Furthermore, the verification results of traditional methods are susceptible to various influencing factors, while the method employed in this invention reduces these factors, thereby improving the accuracy and reliability of the verification.

[0033] See Figure 1 This invention proposes a method for verifying the accuracy of tidal force estimation. By calculating the tidal force at a given measuring station and the force exerted by the Earth on the station, the total change in gravitational acceleration at that station is calculated and matched with data collected by a gravimeter. The station is then located using a gravitational positioning method. The accuracy of the tidal force calculation is evaluated based on the positioning accuracy, thus verifying the accuracy of the tidal force estimation. The implementation process includes the following steps:

[0034] (1) Calculate the tidal force exerted by the celestial body on the measuring station on Earth at time t.

[0035] (2) Calculate the resultant force of the gravitational force and centrifugal force exerted by the Earth on the measuring station on Earth at time t.

[0036] (3) Calculate the magnitude and direction of the total gravitational acceleration at the measuring station on Earth.

[0037] (4) Conduct on-site gravity data collection at the measurement stations.

[0038] (5) Using the universal gravitation positioning method, the station is located by using the calculated gravitational acceleration and the gravitational acceleration collected in the field.

[0039] (6) Compare the positioning results of the two to verify the positioning accuracy. Since a small error in tidal force can cause a huge positioning error, the accuracy of the tidal force calculation can be verified based on the positioning error.

[0040] For ease of implementation and reference, the following detailed description of the preferred implementation of the embodiments is provided:

[0041] First, calculate the change in gravitational acceleration for the first data source.

[0042] Select the point to be determined, set it as the station point, and calculate its initial coordinate values. Therefore, the geocentric distance R between the celestial body and the measuring station can be calculated. J Calculate the zenith distance Z between the celestial body and the station. J When calculating the tidal force exerted on a measuring station on Earth by a celestial body of a certain epoch, the tidal potential V at that measuring station on Earth is calculated. (1) It can be represented as:

[0043]

[0044] Among them, GM J P represents the product of the gravitational constant and the mass of celestial body J. n () represents the nth-order Legendre function, R J and r represent the geocentric distances of the celestial body and the measuring station, respectively, and Z J This is the geocentric zenith distance between the celestial body and the observation station.

[0045] The influence of Earth's oblateness on tidal level V (2) for:

[0046]

[0047] Where J2 represents the second-order zonal harmonic coefficient of the Earth's gravitational field, (α J δ J ), These represent the geocentric longitude and geocentric latitude of the celestial body and the station in the ITRS (International Terrestrial Reference System), respectively, where 'a' is the major radius of the Earth's reference ellipsoid. For example, the fully normalized Legendre function of order n and degree m, This represents a first-order, zero-degree fully normalized associated Legendre function.

[0048] Then the total tidal potential V of celestial bodies at terrestrial stations is equal to V. (1) and V (2) sum:

[0049] V = V (1) +V (2)

[0050] The relationship between tidal force and tidal potential V is denoted as F. V :

[0051]

[0052] Where l represents the normal direction of the tidal location where the measuring station is located.

[0053] The tidal depth V can also be calculated and verified using the following formula:

[0054]

[0055] Among them, GM J R is the product of the gravitational constant and the mass of celestial body J. J and r represent the geocentric distances of the celestial body and the measuring station, respectively, and Z J P is the geocentric zenith distance between the celestial body and the station. n () is the nth-order Legendre function, J2 represents the second-order spherical harmonic coefficient of the Earth's gravitational field, (α) J δ J ), denoted as geocentric longitude and geocentric latitude of the celestial body and the station in the ITRS (International Terrestrial Reference System), respectively, and a is the major radius of the Earth's reference ellipsoid.

[0056] When considering the Earth's oblateness, the resultant force F exerted by the Earth on the measuring stations on Earth can be calculated using the following formula. t :

[0057] F t =F g -F c

[0058] Among them, F g F represents universal gravitation. c This represents centrifugal force. Specifically, the formula for universal gravitation is:

[0059]

[0060] Where G is the gravitational constant, M is the mass of the Earth, m is the mass of the measuring station, r is the distance from the Earth's center to the measuring station, and f is the Earth's oblateness. To measure the latitude of the station.

[0061] The formula for centrifugal force is:

[0062]

[0063] Where ω is the angular velocity of Earth's rotation, r is the distance from the measuring station to the Earth's axis, and f is the Earth's oblateness. Here is the latitude of the measuring station. Substituting these two formulas into the resultant force formula yields the resultant force exerted by the Earth on the measuring station on Earth. These formulas are of great significance in many studies in geophysics and earth sciences, and are therefore frequently cited and used. For ease of implementation and reference, they are provided here.

[0064] The formula for calculating the magnitude and direction of the total gravitational acceleration at each measuring station on Earth (considering tidal force V and Earth's oblateness) is as follows:

[0065]

[0066] Where g is the total gravitational acceleration at the measuring station on Earth.

[0067] After calculating the gravitational acceleration, data from different times can be selected to calculate the gravitational acceleration at different times, thereby calculating the change in gravitational acceleration, which includes changes in magnitude and direction.

[0068] Next, we calculate the change in gravitational acceleration from the second data source.

[0069] The invention employs a gravitational positioning method to measure changes in gravitational acceleration at a test site. The positioning method is as follows: based on the combined effects of gravitational forces generated by celestial bodies, gravitational forces generated by the Earth, and other inertial forces, the gravitational acceleration at the location to be located changes. By periodically monitoring the changes in gravitational acceleration at the location to be located, the positioning result of the location is obtained by reverse calculation.

[0070] This embodiment specifically describes a method for measuring the change in the direction of gravitational acceleration at a location to be located, which includes the following steps:

[0071] Based on the inclinometer's accuracy level, different time intervals are selected to observe the inclinometer's angles θ1, θ2, θ3, θ4 at four times t1, t2, t3, t4, and the differences Δθ1, Δθ2, Δθ3 between adjacent observations are calculated, where:

[0072] Δθ1=θ2-θ1

[0073] Δθ2=θ3-θ2

[0074] Δθ3=θ4-θ3

[0075] Assuming the positions of the celestial bodies and the coordinates of the location to be determined can be accurately obtained at each moment, the directions of the tidal forces V1, V2, V3, and V4 at each moment can be calculated using the law of universal gravitation, and the differences between adjacent directions ΔV1, ΔV2, and ΔV3 can be calculated, where:

[0076] ΔV1=V2-V1

[0077] ΔV2=V3-V2

[0078] ΔV3=V4-V3

[0079] Assuming the location to be determined is the actual location of the inclinometer, then Δθ1, Δθ2, and Δθ3 should be equal to ΔV1, ΔV2, and ΔV3, respectively. Since the directional differences are non-linear, they can be solved using the following system of equations:

[0080] f(Δθ1,Δθ2,Δθ3)=ΔV1

[0081] f(Δθ2,Δθ3,ΔV1)=ΔV2

[0082] f(Δθ3,ΔV1,ΔV2)=ΔV3

[0083] Where t1, t2, t3, and t4 are known precise observations, Δθ1, Δθ2, and Δθ3 are known observations, and f() is a function of the direction difference with respect to θ. This system of equations has three unknowns, and therefore can be solved using Taylor expansion, iteration, or other methods.

[0084] Furthermore, errors in the magnitude and direction of the tidal force exerted by celestial bodies on the measuring station on Earth can lead to significant deviations in the positioning results. Therefore, this invention utilizes the relationship between tidal force and tidal position for positioning. Specifically, the tidal position can be calculated using the following formula: V = V(1) + V(2), where V(1) represents the tidal position of the measuring station on Earth, and V(2) represents the influence of the Earth's oblateness on the tidal position. V(1) can be calculated based on the universal gravitation formula and known information such as the location of the celestial body and the measuring station.

[0085] The magnitude change of gravitational acceleration can be easily obtained using a gravimeter, which will not be elaborated upon in this invention. The preferred embodiment is:

[0086] The change in the direction of gravitational acceleration at the monitoring station is measured using an inclinometer.

[0087] The magnitude of the change in gravitational acceleration at the monitoring station is measured using a gravimeter.

[0088] When periodically monitoring changes in gravitational acceleration at the location of the monitoring station, an atomic clock is used to obtain the current time.

[0089] After obtaining the magnitude and direction changes of gravitational acceleration from the two data sources, the gravitational positioning method can be used to locate the measuring station separately. By monitoring the changes in gravitational acceleration at the location to be located, the positioning result can be inferred. The magnitude and direction changes of gravitational acceleration are estimated based on the initial solution of the station's own coordinates. Combined with the observed data, the positioning result is obtained through iterative linearization.

[0090] When monitoring the change in the magnitude of gravitational acceleration, let time t be... i The magnitude of gravity is g i The change in the magnitude of gravitational acceleration can be expressed as Δg i =g i+1 -g i Let i = 1…k; assuming that according to the coordinates… Given the time and the position of the relevant celestial body, we can calculate the value of time t at each time. i Time to t i+1 The change in the magnitude of gravitational acceleration ΔA at time t. i If i = 1…k, where k is an integer not less than 3, then the coordinates of the point can be deduced by reverse calculation.

[0091] The specific formula for the change in the direction of gravitational acceleration is as follows:

[0092]

[0093] in, Let t1, t2, ..., t be the coordinates of the measurement stations. k+1 ΔV at each time point i =V i+1 -V i For the estimated data of the change in direction of gravitational acceleration, f1() is the difference in the direction of acceleration ΔV. i about The function, i = 1…k, where k is an integer not less than 3.

[0094] The specific formula for the change in the magnitude of gravitational acceleration is as follows:

[0095]

[0096] in, Let t1, t2, ..., t be the coordinates of the measurement stations. k+1 ΔA represents the values ​​at each time point. i =A i+1 -A i For the estimated data of the change in magnitude of gravitational acceleration, f2() is the difference in magnitude of acceleration ΔA. i about The function, i = 1…k, where k is an integer not less than 3.

[0097] By periodically monitoring the changes in gravitational acceleration at the location to be determined, the positioning result can be deduced. This is achieved by calculating estimated data on the magnitude and direction of gravitational acceleration changes based on the initial solution of the coordinates, and then combining this with observational data. The equations governing the magnitude and direction of gravitational acceleration are iteratively linearized and solved until the error converges to a preset limit.

[0098] Linearize the two sets of equations together, take the previously obtained initial values ​​and perform iterative calculations. Stop the iteration when the preset limit is met, and you can get the coordinate results of the measurement station.

[0099] Preferably:

[0100] The initial solution of its own coordinates is based on the current approximate position obtained from an inertial positioning device or a geomagnetic field positioning device.

[0101] The data sources for the changes in gravitational acceleration include field observations and calculations that take tidal forces into account.

[0102] The method for calculating the magnitude and direction of changes in gravitational acceleration includes the following steps:

[0103] (1) To calculate the acceleration and direction caused by the Earth, we need to use the initial solution of its own coordinates, the mass of the Earth and the gravitational constant.

[0104] (2) To calculate the acceleration and direction caused by the Earth's rotation, we need to use the initial solution of its own coordinates and the Earth's rotation speed.

[0105] (3) To calculate the acceleration and direction in the Earth-Moon system, we need to use the initial solutions of time, lunar mass, lunar coordinates, gravitational constant, geocentric coordinates and self-coordinates.

[0106] (4) To calculate the acceleration and direction in the Sun-Earth system, we need to use the initial solutions of time, solar mass, solar coordinates, gravitational constant, geocentric coordinates and self-coordinates.

[0107] (5) To calculate the magnitude and direction of acceleration caused by tidal forces and other influences, it is necessary to use the initial solutions of time, solar mass, solar coordinates, lunar mass, lunar coordinates, other celestial bodies (Venus, Mars, etc.) mass, other celestial bodies coordinates, gravitational constant, geocentric coordinates and self-coordinates.

[0108] (6) Combine the above acceleration vectors to form the magnitude and direction of acceleration, and calculate the corresponding difference.

[0109] (7) Iteratively linearize the solution by using the equations for the change in direction and / or the change in magnitude of gravitational acceleration until the error is less than the set error limit, and obtain the positioning result.

[0110] Based on the key technology proposed in this invention, the coordinate results from the two data sources are compared to verify the accuracy of the tidal force calculation. The accuracy of the positioning result can be represented by three parameters: planar accuracy, elevation accuracy, and three-dimensional accuracy. Among them, planar accuracy and elevation accuracy are the coordinate errors of a specified point in the horizontal and vertical directions, while three-dimensional accuracy is a combination of planar accuracy and elevation accuracy.

[0111] The obtained positions of the measuring stations and celestial bodies are in a geocentric coordinate system, which requires coordinate system transformation for accuracy evaluation. The coordinates of the measuring stations in the geocentric coordinate system are transformed to obtain their plane coordinates and elevations. In practice, existing technologies can be used for this transformation, which will not be elaborated upon in this invention. After obtaining the plane coordinates and elevations from the two data sources, the differences between them are evaluated, which can be divided into plane evaluation and elevation evaluation:

[0112] Horizontal accuracy can be assessed using CEP (Circular Error Probability) or RMSE (Root Mean Square Error). CEP refers to the average distance of all measurement points within the error radius with a 50% probability. Assume there are N measurement points, and the coordinates of the i-th measurement point are (x... i ,yi), its actual coordinates are (x' i y' i If ), then CEP can be calculated as:

[0113]

[0114] RMSE is the square root of the average of the sum of squares of the deviations at all measurement points. Again, assuming there are N measurement points, and the coordinates of the i-th measurement point are (x... i y i Its actual coordinates are (x') i y' i If ), then RMSE can be calculated as:

[0115]

[0116] Vertical accuracy is generally assessed using elevation accuracy, typically employing both absolute and relative accuracy methods. In practical applications, the Height Accuracy Index (HPE) is commonly used to evaluate vertical accuracy. The calculation method for HPE is similar to that for CEP, but it requires converting horizontal errors into vertical errors, i.e.:

[0117]

[0118] Among them, h i h' is the elevation of the measurement point. i N represents the actual elevation of the corresponding point, and N is the total number of measurement points.

[0119] The minute errors in tidal force estimation, after iterative calculation using the aforementioned equations, can lead to significant differences in point coordinates. Therefore, the accuracy of the tidal force estimation can be assessed simply by considering the values ​​of CEP, RMSE, and HPE. This method amplifies minute errors in tidal force estimation into more significant point errors, facilitating the detection of these errors and verifying the accuracy of the tidal force estimation.

[0120] In specific implementation, the method proposed in the technical solution of this invention can be automatically executed by those skilled in the art using computer software technology. System devices for implementing the method, such as computer-readable storage media storing the corresponding computer program of the technical solution of this invention and computer equipment including the computer program running the corresponding computer program, should also be within the protection scope of this invention.

[0121] In some possible embodiments, a tidal force estimation accuracy verification system is provided, including a processor and a memory, wherein the memory is used to store program instructions, and the processor is used to call the stored instructions in the memory to execute a tidal force estimation accuracy verification method as described above.

[0122] In some possible embodiments, a tidal force estimation accuracy verification system is provided, including a readable storage medium on which a computer program is stored, wherein when the computer program is executed, it implements a method for verifying the accuracy of tidal force estimation as described above.

[0123] In some possible embodiments, multiple inclinometers can be used to perform differential calculations to eliminate certain errors and obtain more accurate and reliable results on the change in the direction of gravitational acceleration.

[0124] The specific embodiments described herein are merely illustrative of the spirit of the invention. Those skilled in the art to which this invention pertains may make various modifications or additions to the described specific embodiments or use similar methods to substitute them, without departing from the spirit of the invention or exceeding the scope defined by the appended claims.

Claims

1. A method for verifying the accuracy of tidal force estimation, characterized in that: By calculating the tidal force and the Earth's force on the station, the total change in gravitational acceleration at that station is calculated and matched with data collected by a gravimeter. The station is then located using a gravitational positioning method. The accuracy of the tidal force calculation is evaluated based on the positioning accuracy, thus verifying the accuracy of the tidal force estimation. The process includes the following steps. (1) Calculate the tidal force exerted by the celestial body on the measuring station on Earth at time t; (2) Calculate the resultant force of the gravitational force and centrifugal force exerted by the Earth on the measuring station on Earth at time t; (3) Calculate the magnitude and direction of the total gravitational acceleration at the measuring stations on Earth; (4) Conduct on-site gravity data collection at the measurement stations; (5) Using the gravitational positioning method, the station is located by using the calculated gravitational acceleration and the gravitational acceleration collected in the field. (6) Compare the two positioning results obtained in (5) to verify the positioning accuracy. Since a small tidal force error can cause a huge positioning error, the accuracy of the tidal force calculation can be verified based on the positioning error.

2. The method for verifying the accuracy of tidal force estimation according to claim 1, characterized in that: The method for calculating the tidal force of the celestial body on the station on the earth in step (1) is that the total tidal potential V of the celestial body on the station on the earth is equal to the sum of V (1) and V (2) , wherein V (1) is the tidal potential of the celestial body on the station on the earth at a certain epoch, and V (2) is the influence of the earth's flattening on the tidal potential.

3. The method for verifying the accuracy of tidal force estimation according to claim 1, characterized in that: The calculation method for the resultant force of the gravitational force and centrifugal force exerted by the Earth on the measuring station on Earth in step (2) is as follows: When considering the Earth's oblateness, the resultant force exerted by the Earth on the measuring station on Earth is calculated using the following formula: F t =F g -F c Among them, F t F represents universal gravitation. c This represents centrifugal force.

4. The method for verifying the accuracy of tidal force estimation according to claim 1, characterized in that: The formula for calculating the magnitude and direction of the total gravitational acceleration at each measuring station on Earth, taking into account tidal forces and the Earth's oblateness, is as follows: Where g is the total gravitational acceleration at the measuring station on Earth, and F t V represents universal gravitation, V is the total tidal potential exerted by the celestial body on the measuring station on Earth, and m is the mass of the measuring station.

5. The method for verifying the accuracy of tidal force estimation according to claim 1, characterized in that: The gravitational positioning method in step (5) is implemented by combining the gravitational force generated by celestial bodies, the gravitational force generated by the Earth, and other inertial forces to cause a change in the gravitational acceleration of the location to be located. The positioning result of the location is obtained by periodically monitoring the change in the gravitational acceleration of the location to be located.

6. The method for verifying the accuracy of tidal force estimation according to claim 1, 2, 3, 4, or 5, characterized in that: The accuracy of the positioning result in step (6) is represented by planar accuracy, elevation accuracy and three-dimensional accuracy. Planar accuracy and elevation accuracy are the coordinate errors of the specified point in the horizontal plane and vertical direction, while three-dimensional accuracy is a combination of planar accuracy and elevation accuracy.

7. The method for verifying the accuracy of tidal force estimation according to claim 6, characterized in that: Horizontal accuracy is assessed using circular probability error or root mean square error.

8. The method for verifying the accuracy of tidal force estimation according to claim 6, characterized in that: Vertical accuracy is assessed using elevation accuracy, which can be either absolute or relative.

9. A system for verifying the accuracy of tidal force estimation, characterized in that: It includes a processor and a memory, the memory being used to store program instructions, and the processor being used to call the stored instructions in the memory to execute the tidal force estimation accuracy verification method as described in any one of claims 1-8.

10. A system for verifying the accuracy of tidal force estimation, characterized in that: The method includes a readable storage medium on which a computer program is stored, and when the computer program is executed, it implements a method for verifying the accuracy of tidal force estimation as described in any one of claims 1-8.

Images