Marine gravity measurement error compensation method and system

Abstract

The invention discloses an ocean gravity measurement error compensation method and system, and belongs to the technical field of gravity measurement, and the method comprises the following steps: obtaining topographic parameters of a pre-measured ocean region and water flow fluctuation parameters of the pre-measured ocean region, and carrying out the comprehensive analysis to obtain a comprehensive complex index of the pre-measured ocean region; obtaining an ocean gravity measurement constraint parameter set; obtaining ocean gravity difference indexes, obtaining ocean gravity revision coefficients through matching, synchronously obtaining ocean tide parameters and a preset ocean tide influence prediction model, inputting the ocean tide parameters into the ocean tide influence prediction model to obtain ocean tide influence indexes, and obtaining ocean tide revision coefficients through matching; the ocean gravity correction data set is obtained based on the ocean gravity error correction coefficient and the ocean gravity data set analysis, and the problem that in the prior art, ocean gravity measurement is inaccurate due to complex submarine topography and tidal fluctuation interference is solved.

Description

Technical Field

[0001] The present invention relates to the field of gravity measurement technology, and in particular to a method and system for compensating an error in ocean gravity measurement. Background Art

[0002] With the deepening of marine scientific research, marine gravity measurement, as an important means of studying marine geology, seabed topography and ocean current changes, has gradually become one of the key technologies for marine monitoring. The existing marine gravity measurement error correction system is achieved through zero drift compensation, horizontal acceleration error compensation and other methods.

[0003] For example, the invention patent with announcement number: CN113341476B discloses a method for improving the spatial resolution of ocean gravity based on seabed topography-gravity combination, which includes: obtaining the theoretical gravity value at the seabed control point after performing altitude correction and intermediate layer correction on the normal gravity field of the geoid; determining the relationship among the free space gravity anomaly at the seabed control point, the gravity anomaly caused by water depth change and the free space gravity anomaly on the geoid; obtaining the basic formula model of the seabed topography-gravity combination method; solving to obtain the absolute gravity value at the discrete seabed control point; solving to obtain the gridded seabed absolute gravity value; obtaining the seabed topography-gravity combination method model through the existing gridded high-precision seabed topography data and gravity correction, and solving to obtain the gridded free space gravity anomaly on the geoid.

[0004] For example, the invention patent with announcement number: CN106405670B announces a gravity anomaly data processing method suitable for a strapdown marine gravimeter, including: using a Kalman filter to estimate the attitude error in the inertial navigation solution, correcting the attitude matrix and the vertical specific force component to obtain more accurate specific force information, and then calculating the gravity correction items based on the position, speed and altitude information provided by the PPP technology, and finally obtaining the gravity anomaly information along the route through a low-pass filter.

[0005] However, in the process of implementing the technical solution of the invention in the embodiments of the present application, the present application found that the above technology has at least the following technical problems: In the prior art, error compensation is only performed through zero drift compensation, horizontal acceleration error compensation and other methods, but these measures are still difficult to completely eliminate errors in the face of complex marine environments. In particular, factors such as the undulations of the seabed topography, the complexity of water flow fluctuations, and the periodic changes of tides will have a significant impact on gravity measurement. Therefore, there is currently a problem of inaccurate marine gravity measurement caused by complex seabed topography and tidal fluctuations. Summary of the invention

[0006] The embodiments of the present application solve the problem of inaccurate ocean gravity measurement caused by complex seabed topography and tidal fluctuation interference in the prior art by providing an ocean gravity measurement error compensation method and system, thereby achieving accurate analysis of ocean gravity.

[0007] The embodiment of the present application provides an ocean gravity measurement error compensation method, comprising the following steps: obtaining terrain parameters of a predicted ocean area and water flow fluctuation parameters of the predicted ocean area, and comprehensively analyzing to obtain comprehensive complex indicators of the predicted ocean area; obtaining an ocean gravity measurement constraint parameter set based on the analysis of the comprehensive complex indicators of the predicted ocean area; performing ocean gravity measurement based on the ocean gravity measurement constraint parameter set, obtaining an ocean gravity data set, analyzing to obtain an ocean gravity difference index, thereby matching to obtain an ocean gravity revision coefficient, synchronously obtaining ocean tide parameters and a preset ocean tide impact prediction model, and inputting the ocean tide parameters into the ocean tide impact prediction model to obtain an ocean tide impact index, thereby matching to obtain an ocean tide revision coefficient; performing analysis based on the ocean tide revision coefficient and the ocean gravity revision coefficient to obtain an ocean gravity error revision coefficient, and analyzing based on the ocean gravity error revision coefficient and the ocean gravity data set to obtain an ocean gravity correction data set.

[0008] Furthermore, the method of obtaining the topographic parameters of the ocean area to be measured and the water flow fluctuation parameters of the ocean area to be measured, and comprehensively analyzing to obtain the comprehensive complexity index of the ocean area to be measured, specifically comprises the following steps: obtaining the topographic parameters of the ocean area to be measured, the topographic parameters of the ocean area to be measured include the average seabed slope, the extreme difference of seabed terrain height, the average thickness of seabed rock formations and the average hardness of seabed rock formations; obtaining the water flow fluctuation parameters of the ocean area to be measured, the water flow fluctuation parameters of the ocean area to be measured include the maximum water flow velocity, the extreme difference of water flow velocity and the turbulence intensity at each screening point within a preset time period; obtaining the topographic reference set and the water flow fluctuation reference set preset in the database, and comparing them with the topographic parameters of the ocean area to be measured and the water flow fluctuation parameters of the ocean area to be measured, respectively, to obtain the comprehensive complexity index of the ocean area to be measured; the comprehensive complexity index of the ocean area to be measured is used to characterize the complexity of the ocean area to be measured; the topographic reference set includes a slope reference value, a terrain height extreme difference reference value, a seabed rock formation thickness reference value and a seabed rock formation hardness reference value; the water flow fluctuation reference set includes a water flow velocity reference value, a water flow velocity extreme difference reference value and a turbulence intensity reference value.

[0009] Furthermore, the comprehensive complexity index of the predicted ocean area includes the following specific analysis steps: analyzing based on the terrain parameters and terrain reference set of the predicted ocean area to obtain the terrain complexity index of the predicted ocean area; analyzing based on the water flow fluctuation parameters and water flow fluctuation reference set of the predicted ocean area to obtain the water flow fluctuation complexity index of the predicted ocean area; analyzing based on the terrain complexity index of the predicted ocean area and the water flow fluctuation complexity index of the predicted ocean area to obtain the comprehensive complexity index of the predicted ocean area; the terrain complexity index of the predicted ocean area is used to characterize the terrain complexity of the predicted ocean area; the water flow fluctuation complexity index of the predicted ocean area is used to characterize the water flow fluctuation complexity of the predicted ocean area.

[0010] Furthermore, the analysis of the comprehensive complexity indicators of the predicted ocean area is based on the analysis to obtain the ocean gravity measurement constraint parameter set, and the specific steps include: obtaining each comprehensive complexity indicator interval preset in the database and the ocean gravity measurement constraint reference parameter set corresponding to each comprehensive complexity indicator interval, and comparing with the comprehensive complexity indicator of the predicted ocean area; if the comprehensive complexity indicator of the predicted ocean area is within a preset comprehensive complexity indicator interval, then obtaining the ocean gravity measurement constraint reference parameter set corresponding to the comprehensive complexity indicator interval as the ocean gravity measurement constraint parameter set; the ocean gravity measurement constraint parameter set includes the number of ocean gravimeters deployed and the measurement constraint period.

[0011] Furthermore, the method of performing ocean gravity measurement based on the ocean gravity measurement constraint parameter set, obtaining an ocean gravity data set, and analyzing to obtain an ocean gravity difference index specifically comprises the following steps: deploying ocean gravity meters based on the ocean gravity measurement constraint parameter set, thereby performing ocean gravity measurement, and obtaining an ocean gravity data set, wherein the ocean gravity data set includes gravity acceleration measured by each ocean gravity meter at each time, and a center distribution distance and a vertical distribution height of each gravity meter; obtaining an ocean gravity fitting data set, and analyzing the data set with the ocean gravity data set to obtain an ocean gravity difference index; the ocean gravity fitting data set includes a gravity acceleration fitting value, a center distribution distance fitting value, and a vertical distribution height fitting value; and the ocean gravity difference index is used to characterize the degree of measurement difference of the ocean gravity meters.

[0012] Furthermore, the matching obtains the ocean gravity revision coefficient, and the specific steps include: obtaining each ocean gravity difference index interval preset in the database and the compensation reference factor corresponding to each ocean gravity difference index interval; based on the ocean gravity difference index and comparing with each ocean gravity difference index interval, if the ocean gravity difference index is within a preset ocean gravity difference index interval, then obtaining the compensation reference factor corresponding to the interval as the ocean gravity revision coefficient; the ocean gravity revision coefficient is used to revise the ocean gravity.

[0013] Furthermore, the ocean tidal parameters are input into the ocean tidal impact prediction model to obtain the ocean tidal impact index, and the ocean tidal revision coefficient is obtained by matching. The specific steps include: obtaining ocean tidal parameters, the ocean tidal parameters include the tidal amplitude, the maximum tidal flow velocity and the tidal water level change rate of the predicted ocean area in the measurement constraint period; inputting the ocean tidal parameters into the ocean tidal impact prediction model to obtain the ocean tidal impact index; the ocean tidal impact index is used to characterize the degree of influence of the ocean tide; obtaining each ocean tidal impact index interval preset in the database and the ocean tidal reference revision coefficient corresponding to each ocean tidal impact index interval, and comparing them with the ocean tidal impact index, if the ocean tidal impact index is within a certain ocean tidal impact index interval, then obtaining the ocean tidal reference revision coefficient corresponding to the ocean tidal impact index interval as the ocean tidal revision coefficient, the ocean tidal revision coefficient is used to reduce the influence of the ocean tide on the measurement of the ocean gravity acceleration.

[0014] Furthermore, the ocean tidal impact prediction model is specifically: ; In the formula, Indicates the ocean tide impact index, represents the tidal amplitude of the predicted ocean area during the measurement constraint period, represents the unit tidal amplitude influence factor, represents the maximum tidal velocity in the predicted ocean area during the measurement constraint period, Indicates the maximum flow velocity factor per unit tide, represents the rate of change of tidal level in the predicted ocean area during the measurement constraint period, Represents the factor affecting the rate of change of unit tidal water level.

[0015] Furthermore, the ocean gravity correction data set is obtained based on the analysis of the ocean gravity error revision coefficient and the ocean gravity data set, and the specific steps include: extracting the gravity acceleration measured by each ocean gravimeter in the ocean gravity data set, and performing discrete mean processing to obtain the specified gravity acceleration; based on the ocean gravity error revision coefficient, and analyzing the ocean gravity data set to obtain the gravity acceleration correction adjustment value, and combining the specified gravity acceleration to construct a specified gravity acceleration interval for the predicted ocean area; the specified gravity acceleration and the specified gravity acceleration interval for the predicted ocean area are jointly recorded as the ocean gravity correction data set.

[0016] The embodiment of the present application provides an ocean gravity measurement error compensation system, including: a comprehensive and complex analysis module, a constraint judgment module, a revision coefficient evaluation module and a correction execution module; wherein the comprehensive and complex analysis module is used to obtain the terrain parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area, and comprehensively analyze to obtain the comprehensive and complex indicators of the predicted ocean area; the constraint judgment module is used to obtain the ocean gravity measurement constraint parameter set based on the comprehensive and complex indicator analysis of the predicted ocean area; the revision coefficient evaluation module is used to perform ocean gravity measurement based on the ocean gravity measurement constraint parameter set, obtain the ocean gravity data set, analyze to obtain the ocean gravity difference index, thereby matching to obtain the ocean gravity revision coefficient, synchronously obtain the ocean tide parameters and the preset ocean tide influence prediction model, and input the ocean tide parameters into the ocean tide influence prediction model to obtain the ocean tide influence index, thereby matching to obtain the ocean tide revision coefficient; the correction execution module is used to analyze based on the ocean tide revision coefficient and the ocean gravity revision coefficient to obtain the ocean gravity error revision coefficient, and analyze based on the ocean gravity error revision coefficient and the ocean gravity data set to obtain the ocean gravity correction data set.

[0017] One or more technical solutions provided in the embodiments of the present application have at least the following technical effects or advantages: 1. The ocean gravity measurement error compensation method provided by the present invention obtains the terrain parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area, and comprehensively analyzes the comprehensive complex indicators of the predicted ocean area, thereby obtaining the ocean gravity measurement constraint parameter set, thereby realizing the accurate analysis of ocean gravity and ensuring that the resources invested in the measurement can fully fit the actual characteristics of the ocean, providing reliable guarantee for the effective utilization of resources, and effectively solving the problem of inaccurate ocean gravity measurement caused by complex seabed terrain and tidal fluctuation interference in the prior art.

[0018] 2. The present invention obtains an ocean gravity data set, analyzes it to obtain an ocean gravity revision coefficient, and then analyzes it in combination with an ocean tide revision coefficient to obtain an ocean gravity error revision coefficient. Then, the ocean gravity correction data set is obtained based on the ocean gravity error revision coefficient and the ocean gravity data set, thereby achieving accurate correction of the ocean gravity data.

[0019] 3. By obtaining the ocean tidal parameters and inputting them into the ocean tidal impact prediction model, the ocean tidal impact index is obtained, which is then matched with the preset ocean tidal correction coefficient, thereby realizing the compensation correction of the ocean tidal effect and reducing the interference of tidal changes on the ocean gravity measurement data. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] Figure 1A flow chart of a method for compensating ocean gravity measurement errors provided in an embodiment of the present application; Figure 2 A schematic diagram of the modules of the ocean gravity measurement error compensation system provided in an embodiment of the present application. DETAILED DESCRIPTION

[0021] The embodiments of the present application solve the problem of inaccurate ocean gravity measurement caused by complex seabed topography and tidal fluctuation interference in the prior art by providing an ocean gravity measurement error compensation method and system, thereby achieving accurate analysis of ocean gravity.

[0022] In order to better understand the above technical solution, the above technical solution will be described in detail below in conjunction with the accompanying drawings and specific implementation methods.

[0023] like Figure 1 As shown, it is a flow chart of the ocean gravity measurement error compensation method provided by an embodiment of the present application, and the method includes the following steps: obtaining the terrain parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area, and comprehensively analyzing to obtain the comprehensive complex index of the predicted ocean area; based on the comprehensive complex index analysis of the predicted ocean area, obtaining the ocean gravity measurement constraint parameter set; performing ocean gravity measurement based on the ocean gravity measurement constraint parameter set, obtaining the ocean gravity data set, analyzing to obtain the ocean gravity difference index, thereby matching to obtain the ocean gravity revision coefficient, synchronously obtaining the ocean tide parameters and the preset ocean tide influence prediction model, and inputting the ocean tide parameters into the ocean tide influence prediction model to obtain the ocean tide influence index, thereby matching to obtain the ocean tide revision coefficient; analyzing based on the ocean tide revision coefficient and the ocean gravity revision coefficient to obtain the ocean gravity error revision coefficient, and analyzing based on the ocean gravity error revision coefficient and the ocean gravity data set to obtain the ocean gravity correction data set.

[0024] In this embodiment, by analyzing and predicting the comprehensive and complex indicators of the ocean area, the error compensation measurement scheme can be adjusted according to the terrain and water flow fluctuation characteristics of different ocean areas, and the accurate analysis of ocean gravity is realized, and the resources invested in the measurement can fully fit the actual characteristics of the ocean, providing reliable guarantee for the effective use of resources. The accuracy of error compensation is thus improved, especially in applications in different seasons, different tides and ocean currents, ensuring the high consistency and accuracy of the measurement data.

[0025] Furthermore, the terrain parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area are obtained, and a comprehensive analysis is performed to obtain a comprehensive complexity index of the predicted ocean area. The specific steps include: obtaining the terrain parameters of the predicted ocean area, the terrain parameters of the predicted ocean area include the average seabed slope, the extreme difference of seabed terrain height, the average thickness of seabed rock layers and the average hardness of seabed rock layers; obtaining the water flow fluctuation parameters of the predicted ocean area, the water flow fluctuation parameters of the predicted ocean area include the maximum water flow velocity, the extreme difference of water flow velocity and the turbulence intensity of each screening point within a preset time period; obtaining the preset terrain reference set and water flow fluctuation reference set in the database, and comparing them with the terrain parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area, respectively, to obtain the comprehensive complexity index of the predicted ocean area; the comprehensive complexity index of the predicted ocean area is used to characterize the complexity of the predicted ocean area; the terrain reference set includes a slope reference value, a terrain height extreme difference reference value, a seabed rock layer thickness reference value and a seabed rock layer hardness reference value; the water flow fluctuation reference set includes a water flow velocity reference value, a water flow velocity extreme difference reference value and a turbulence intensity reference value.

[0026] In this embodiment, by analyzing the seabed topography and water flow fluctuation parameters, the geological and hydrological environment of the ocean area can be reflected in more detail, and the errors caused by ignoring these detail factors in the traditional measurement method can be reduced, which not only improves the accuracy of the data, but also can provide more targeted error correction schemes according to the characteristics of different ocean areas. By combining multiple key parameters and comparing them with the reference set in the database, the comprehensive complex index of the ocean area is obtained, the complexity of the ocean area is effectively quantified, and the analysis of ocean gravity data is made more efficient. Combined with the predicted seabed topography and water flow fluctuation parameters, this method can flexibly respond to changes in different regions, different seasons, and different tidal conditions. It has strong adaptability for compensating for ocean gravity measurement errors in complex environments and can be widely used in diverse ocean surveys and research projects.

[0027] It should be noted that the average seafloor slope, average hardness of seafloor rock formations, seafloor topography height and average thickness of seafloor rock formations can be obtained by querying relevant data from the Geological Survey Bureau. Water velocity and turbulence intensity can be obtained by querying relevant data from the marine management department. The extreme difference in water velocity can be obtained by calculating the difference between the maximum and minimum values ​​of the water velocity.

[0028] It should also be noted that through the comprehensive and complex indicators, the terrain and water flow characteristics of the ocean area can be accurately described, providing a scientific basis for the subsequent correction of gravity measurement errors, and reflecting the multidimensional complexity of the marine environment, especially in areas with more complex terrain and water flow fluctuations. The credibility of marine gravity measurement data is improved. The introduction of comprehensive and complex indicators enables marine gravity measurements to be implemented in a targeted manner according to specific marine environmental conditions, achieving accurate analysis of marine gravity and ensuring that the resources invested in measurement can fully fit the actual characteristics of the ocean, providing reliable guarantees for the effective use of resources.

[0029] Furthermore, the comprehensive complexity index of the ocean area to be measured is analyzed in the following specific steps: based on the terrain parameters and terrain reference set of the ocean area to be measured, the terrain complexity index of the ocean area to be measured is obtained; based on the water flow fluctuation parameters and water flow fluctuation reference set of the ocean area to be measured, the water flow fluctuation complexity index of the ocean area to be measured is obtained; based on the terrain complexity index of the ocean area to be measured and the water flow fluctuation complexity index of the ocean area to be measured, the comprehensive complexity index of the ocean area to be measured is obtained; the terrain complexity index of the ocean area to be measured is used to characterize the terrain complexity of the ocean area to be measured; the water flow fluctuation complexity index of the ocean area to be measured is used to characterize the water flow fluctuation complexity of the ocean area to be measured.

[0030] In this embodiment, it should be noted that by analyzing the terrain parameters of the predicted ocean area, the terrain complexity index of the predicted ocean area is obtained, taking into account the mutual influence relationship between these parameters, for example: areas with large seabed slopes are usually accompanied by large terrain height differences, because areas with large slopes mean that the seabed is undulating more, resulting in a larger change in terrain height, thereby indicating that the terrain of the ocean area is complex, resulting in a greater flow rate of seawater and a more intense flow of seawater. Therefore, the impact when measuring the ocean gravity acceleration will be more intense, resulting in an increased error. The thickness of the seabed rock layer affects the stability of the seabed topography. Thicker rock layers will suppress the slope changes of the seabed, making the slope gentler. Thin rock layers lead to steeper slope changes, resulting in larger terrain undulations. Harder rock layers can resist erosion and sedimentation, thereby forming steeper terrain, resulting in a larger terrain height difference. Soft rock layers are prone to sedimentation or erosion, forming a flatter terrain with a smaller terrain height difference.

[0031] The complex index of water flow fluctuation in the predicted ocean area is obtained by analyzing the water flow fluctuation parameters of the predicted ocean area, taking into account the mutual influence relationship between these parameters. For example, in areas with large water flow speed, the fluctuation of water flow is usually more violent, so the range of water flow speed will increase. Especially in areas with obvious tidal changes, the maximum water flow speed and the range of water flow speed are usually positively correlated. Turbulence intensity is usually more obvious in areas with large water flow speed. The greater the water flow speed, the stronger the turbulence phenomenon, and the turbulence intensity increases accordingly. The water flow in areas with high turbulence intensity changes more rapidly, resulting in a larger range of water flow speed. The increase in turbulence intensity indicates an increase in the instability of the water flow, which in turn leads to an increase in the speed range.

[0032] The terrain complexity index can accurately describe the terrain changes in the ocean area by quantitatively analyzing the average seafloor slope, the extreme difference in terrain height, the thickness of the rock layer and the hardness of the rock layer. Through this index, the topographic undulations and complexity of the seafloor can be more clearly identified, thereby providing more accurate basic data for subsequent marine gravity measurements and error compensation. In marine gravity measurements, the undulations of the terrain directly affect the accuracy of the measurement results. Through the terrain complexity index, the subtle changes in the seafloor topography can be deeply analyzed. During the seafloor survey, the terrain complexity index can be used to more efficiently select key measurement areas. This indicator helps researchers identify areas with large terrain changes, and provides data support for accurately locating gravity measurement points, selecting the best sampling location, and formulating reasonable exploration plans, thereby improving measurement efficiency and quality. The terrain complexity of different ocean areas varies greatly. After obtaining the terrain complexity index, customized analysis and compensation can be performed according to the specific characteristics of the area. For example, in areas with large seafloor slopes, the terrain complexity index can provide more guidance information, help adjust the measurement strategy, and improve the adaptability of the measurement strategy in complex environments.

[0033] By analyzing the parameters of water flow fluctuations, the complex characteristics of the ocean area can be fully and comprehensively reflected. It can not only provide multi-dimensional data of the seabed geology and hydrological environment, but also provide a unified representation of the characteristics of the entire ocean area, which helps to quickly identify the complexity of different sea areas, provide a quantitative basis for the correction of gravity measurement errors, and help reduce the error impact caused by ignoring a single factor. Especially when the terrain complexity and water flow fluctuations in the ocean area are significant, the comprehensive complexity index can make a more refined error correction through a unified analysis of the two.

[0034] The terrain complexity index of the predicted ocean area is obtained by: ; In the formula, Represents the terrain complexity index of the predicted ocean area, represents the average seafloor slope, Indicates the slope reference value, It means that the seabed topography is extremely different. Indicates the reference value of the extreme difference in terrain height. represents the average thickness of the seafloor rock layer, Indicates the reference value of the thickness of the seabed rock layer, Indicates the average hardness of the seabed rock. represents the reference value of the hardness of the seabed rock formation, e represents the natural constant, represents the seafloor slope influence coefficient, represents the influence coefficient of seabed topography, represents the influence coefficient of seabed rock thickness, Indicates the influence coefficient of seabed rock hardness.

[0035] The seabed slope influence coefficient, seabed topography influence coefficient, seabed rock thickness influence coefficient and seabed rock hardness influence coefficient can be obtained from the database, for example: obtain the average seabed slope stored in the database, and construct a seabed slope mapping set by combining the average seabed slope and the seabed slope influence coefficient corresponding to the average seabed slope, wherein there is a one-to-one or many-to-one correspondence in the seabed slope mapping set, and by inputting the real-time average seabed slope into the mapping set, the seabed slope influence coefficient corresponding to the average seabed slope can be obtained. Other influence coefficients such as the seabed topography influence coefficient, the seabed rock thickness influence coefficient and the seabed rock hardness influence coefficient can be obtained through their corresponding mapping sets in the same way, for example, the seabed topography influence coefficient needs to be obtained in the seabed topography mapping set, the seabed rock thickness influence coefficient needs to be obtained in the seabed rock thickness mapping set, and the seabed rock hardness influence coefficient needs to be obtained in the seabed rock hardness mapping set. It should be noted that the construction method of the seabed topography mapping set, the seabed rock thickness mapping set and the seabed rock hardness mapping set is the same as that of the seabed slope mapping set.

[0036] The complex index of water flow fluctuation in the predicted ocean area is obtained by: ; In the formula, Indicates the complexity index of water flow fluctuation in the predicted ocean area, represents the maximum water flow velocity at the i-th screening point, i represents the number of the screening point, , represents the total number of screening points, Indicates the reference value of water flow velocity, Indicates the extreme difference in water velocity at the i-th screening point, Indicates the extreme reference value of water velocity. represents the turbulence intensity at the i-th screening point, represents the reference value of turbulence intensity, represents the water velocity influence coefficient, It represents the influence coefficient of the extreme difference of water velocity, Represents the turbulence intensity influence coefficient.

[0037] The slope reference value, the extreme difference reference value of terrain height, the seabed rock thickness reference value and the seabed rock hardness reference value are obtained by extracting from the database, and the water flow velocity reference value, the extreme difference reference value of water flow velocity and the turbulence intensity reference value are obtained by extracting from the database.

[0038] It should be noted that the water flow velocity influence coefficient, the water flow velocity range reference value and the turbulence intensity reference value can be obtained by extracting from the database. The specific method is: construct a water flow velocity mapping set by obtaining the preset water flow velocity and the influence coefficient corresponding to the water flow velocity in the database, wherein there is a one-to-one correspondence or a many-to-one correspondence in the water flow velocity mapping set. The water flow velocity influence coefficient can be obtained by inputting the real-time water flow velocity into the mapping set. Other influence coefficients such as the water flow velocity range influence coefficient and the turbulence intensity influence coefficient can be obtained by constructing their corresponding mapping sets in the same way. For example, the water flow velocity range influence coefficient needs to be obtained from the water flow velocity range mapping set, and the turbulence intensity influence coefficient needs to be obtained from the turbulence intensity mapping set. It should be noted that the water flow velocity range mapping set and the turbulence intensity mapping set are constructed in the same way as the water flow velocity mapping set.

[0039] The comprehensive and complex indicators of the predicted ocean area are obtained in the following way: ; In the formula, Represents a comprehensive and complex index of the predicted ocean area, Represents the terrain complexity index of the predicted ocean area, It indicates the complexity index of water flow fluctuation in the predicted ocean area.

[0040] Furthermore, based on the analysis of the comprehensive and complex indicators of the predicted ocean area, an ocean gravity measurement constraint parameter set is obtained. The specific steps include: obtaining each comprehensive and complex indicator interval preset in the database and the ocean gravity measurement constraint reference parameter set corresponding to each comprehensive and complex indicator interval, and comparing them with the comprehensive and complex indicators of the predicted ocean area. If the comprehensive and complex indicators of the predicted ocean area are within a preset comprehensive and complex indicator interval, then the ocean gravity measurement constraint reference parameter set corresponding to the comprehensive and complex indicator interval is obtained as the ocean gravity measurement constraint parameter set; the ocean gravity measurement constraint parameter set includes the number of ocean gravimeters deployed and the measurement constraint period.

[0041] In this embodiment, the greater the comprehensive complexity index, the greater the corresponding number of deployments and the longer the constraint period. It should be noted that the marine gravity measurement constraint parameter set includes the number of marine gravimeter deployments and the measurement constraint period, wherein the measurement constraint period refers to the duration of the test.

[0042] By performing interval analysis on the comprehensive complex indicators of the ocean area and combining it with the corresponding reference parameter set of ocean gravity measurement constraints, the constraints in the process of ocean gravity measurement can be controlled more accurately. By matching the complex indicators of the predicted ocean area with the preset reference parameter set, the constraints of ocean gravity measurement can be flexibly adjusted for areas of different complexity, thereby improving the measurement accuracy. Moreover, the constraint parameters of ocean gravity measurement can be dynamically adjusted according to the actual conditions of different ocean areas, avoiding the limitations of the static measurement mode. In a complex ocean environment, the interval of comprehensive complex indicators provides personalized constraint parameters for each area, enhancing the adaptability of data processing, especially in areas with large environmental changes, which can effectively compensate for the error of ocean gravity measurement. By optimizing the number of ocean gravity meters and the measurement constraint period, the measurement resources can be reasonably allocated to ensure the comprehensiveness and efficiency of ocean gravity measurement. According to the different intervals of comprehensive complex indicators, the appropriate number of gravity meters and the measurement period are automatically selected, avoiding over-configuration or under-configuration, and maximizing the efficiency of resource use.

[0043] Through automated parameter matching and selection, the need for human intervention is reduced. In traditional marine gravity measurements, the selection of measurement constraint parameters usually depends on the experience of the operator, while the present invention can quickly and accurately select appropriate constraint parameters for each ocean area through a preset database, simplifying the operation process and improving the degree of automation and reliability. Through the automatic matching of comprehensive complex indicators and constraint parameter sets, efficient gravity measurement can be achieved in a wide range of ocean areas. Different areas obtain corresponding constraint parameter sets according to their complexity. The automated process reduces the complexity of manual selection and significantly improves the efficiency and coverage of marine gravity measurements in a wide range of areas. The present application can monitor and adjust the marine gravity measurement constraint parameter set in real time, and dynamically optimize the measurement parameters according to environmental changes and real-time data. For example, under the influence of ocean tidal changes, climate changes or other external factors, the marine gravity measurement constraint reference parameter set can be automatically adjusted through the obtained comprehensive complex indicator interval to ensure that the measurement data can adapt to the current environment.

[0044] Furthermore, ocean gravity measurement is performed based on the ocean gravity measurement constraint parameter set to obtain an ocean gravity data set, and an ocean gravity difference index is obtained by analysis. The specific steps include: deploying ocean gravity meters based on the ocean gravity measurement constraint parameter set, thereby performing ocean gravity measurement and obtaining an ocean gravity data set, the ocean gravity data set including the gravity acceleration measured by each ocean gravity meter at each time and the center distribution distance and vertical distribution height of each gravity meter; obtaining an ocean gravity fitting data set and analyzing it with the ocean gravity data set to obtain an ocean gravity difference index; the ocean gravity fitting data set includes gravity acceleration fitting values, center distribution distance fitting values ​​and vertical distribution height fitting values; the ocean gravity difference index is used to characterize the degree of measurement difference of the ocean gravity meter.

[0045] In this embodiment, it should be noted that the marine gravity fitting data set includes a gravity acceleration fitting value, a center distribution distance fitting value, and a vertical distribution height fitting value, wherein the gravity acceleration fitting value can be obtained by averaging the gravity acceleration measured by each marine gravimeter, obtaining the gravity acceleration of each marine gravimeter and averaging it to obtain the average gravity acceleration of each marine gravimeter, and then averaging the average gravity acceleration of each marine gravimeter measured each time to obtain the gravity acceleration fitting value. The center distribution distance fitting value refers to the average value of the center distribution distance measured by each gravity measuring instrument, and the vertical distribution height fitting value refers to the average value of the vertical distribution height measured by each gravity measuring instrument.

[0046] Gravitational acceleration can be measured by an absolute gravimeter. The center distribution distance refers to the distance between the locations of each ocean gravimeter and the center point ocean gravimeter. It can be measured by GPS positioning system. The vertical distribution height can be measured by an electronic level. The vertical distribution height refers to the vertical height of the gravimeter from the sea level.

[0047] The ocean gravity difference index is obtained by analyzing the ocean gravity data set, taking into account the mutual influence relationship between these parameters. For example, the gravity acceleration directly determines the gravity value at each measurement point, but it is affected by the position of the measuring instrument (i.e., the center distribution distance) and the change in height. The closer the gravimeter is located to the gravity source, the more accurate the measurement of gravity acceleration. Therefore, the shorter the center distribution distance, the more stable and accurate the measured gravity acceleration value is. In addition, the change in the vertical distribution height will also affect the measurement results. Due to the non-uniformity of the earth's gravity field, the measurement points at higher positions will have lower gravity acceleration as the height changes, so the height difference of the vertical distribution needs to be accurately measured and corrected to ensure the accuracy of the data.

[0048] By analyzing the relationship between gravity acceleration, center distribution distance and vertical distribution height, the error caused by the measurement position and height difference can be effectively reduced. Due to the complex marine environment, gravity measurement is often affected by multiple factors. By comprehensively considering the relationship between gravity acceleration, position and vertical height, the error caused by a single factor can be reduced to ensure the stability and consistency of the measurement data. The layout strategy of marine gravity measurement can be optimized. By rationally planning the location, layout and height difference of the measuring instruments, the interference factors in the measurement process can be reduced, ensuring that the measurement work is completed efficiently in a changing marine environment, thereby optimizing the utilization of measurement resources. For example, the measurement method can be adjusted according to the height difference and distribution distance of the measuring point, and the possible error source can be predicted in advance, so as to make targeted error corrections. This process of prediction and correction improves the efficiency and accuracy of error compensation.

[0049] The ocean gravity difference index is obtained by: ; In the formula, represents the ocean gravity difference index, represents the kth measured gravity acceleration of the jth ocean gravimeter, j represents the number of the ocean gravimeter, , represents the total number of ocean gravimeters, k represents the number of measurements, , represents the total number of measurements, represents the fitted value of gravity acceleration, represents the center distribution distance of the jth gravimeter, represents the central distribution distance fitting value, represents the vertical distribution height of the jth gravity meter, represents the vertical distribution height fitting value, represents the influence coefficient of gravity acceleration, represents the center distribution distance influence coefficient, Represents the vertical distribution height influence coefficient.

[0050] It should be noted that the gravity acceleration influence coefficient, the center distribution distance influence coefficient and the vertical distribution height influence coefficient can be obtained through the database, for example: a gravity acceleration mapping set is constructed by obtaining the gravity acceleration preset in the database and the influence coefficient corresponding to the gravity acceleration, wherein the gravity acceleration mapping set has a one-to-one correspondence or a many-to-one correspondence, and the gravity acceleration influence coefficient can be obtained by inputting the real-time gravity acceleration into the mapping set. Other influence coefficients such as the center distribution distance influence coefficient and the vertical distribution height influence coefficient can be obtained through their corresponding mapping sets in the same way. For example, the center distribution distance influence coefficient needs to be obtained in the center distribution distance mapping set, and the vertical distribution height influence coefficient needs to be obtained in the vertical distribution height mapping set. It should be noted that the construction method of the center distribution distance mapping set and the vertical distribution height mapping set is the same as that of the gravity acceleration mapping set.

[0051] Furthermore, the ocean gravity revision coefficient is obtained by matching, and the specific steps include: obtaining each ocean gravity difference index interval preset in the database and the compensation reference factor corresponding to each ocean gravity difference index interval; based on the ocean gravity difference index and comparing with each ocean gravity difference index interval, if the ocean gravity difference index is within a preset ocean gravity difference index interval, then obtaining the compensation reference factor corresponding to the interval as the ocean gravity revision coefficient; the ocean gravity revision coefficient is used to revise the ocean gravity.

[0052] In this embodiment, it should be noted that the larger the ocean gravity difference index is, the larger the ocean gravity revision coefficient is. By obtaining the revision coefficient based on the accurate matching of the ocean gravity difference index and the preset interval, the accuracy and credibility of the ocean gravity measurement data are improved. In this process, the acquisition of the revision coefficient not only performs customized compensation for specific difference indicators, but also effectively ensures that in a complex and changeable ocean environment, the data correction can accurately adapt to the actual measurement conditions, so that each ocean gravity measurement result can dynamically adjust the correction factor according to the real-time difference index, and improve the flexibility and accuracy of error compensation. The process of automated matching optimizes the data processing flow, avoids tedious manual calculations and judgments, shortens the data correction cycle, and provides efficient data support for large-scale ocean gravity measurements. On this basis, the corrected data can be more accurately fused and comprehensively analyzed, which improves the accuracy and reliability of subsequent processing steps, especially in data fusion and multi-source information integration, the corrected data promotes the accuracy of the analysis results.

[0053] Furthermore, the ocean tidal parameters are input into the ocean tidal impact prediction model to obtain the ocean tidal impact index, and the ocean tidal revision coefficient is obtained by matching. The specific steps include: obtaining the ocean tidal parameters, the ocean tidal parameters include the tidal amplitude, the maximum tidal flow velocity and the tidal water level change rate of the predicted ocean area in the measurement constraint period; inputting the ocean tidal parameters into the ocean tidal impact prediction model to obtain the ocean tidal impact index; the ocean tidal impact index is used to characterize the degree of influence of the ocean tide; obtaining each ocean tidal impact index interval preset in the database and the ocean tidal reference revision coefficient corresponding to each ocean tidal impact index interval, and comparing them with the ocean tidal impact index, if the ocean tidal impact index is within a certain ocean tidal impact index interval, then obtaining the ocean tidal reference revision coefficient corresponding to the ocean tidal impact index interval as the ocean tidal revision coefficient, the ocean tidal revision coefficient is used to reduce the influence of the ocean tide on the measurement of the ocean gravity acceleration.

[0054] In this embodiment, it should be noted that the tidal amplitude refers to the vertical distance between the lowest tide level and the highest tide level of the seawater in the tidal phenomenon, that is, the maximum fluctuation value of the tidal wave, which can be measured by a tidal observation instrument (such as a tidal meter), and the tidal meter determines the amplitude of the tide by recording the change of the water level. The maximum tidal flow velocity refers to the maximum value of the water flow velocity in a tidal cycle. The maximum tidal flow velocity can be obtained by measuring the flow velocity through a water flow velocity instrument (such as ADCP, acoustic Doppler current profiler). The tidal water level change rate represents the speed of water level change, which can be obtained through continuous water level measurement. The commonly used equipment is a water level recorder (such as a tidal meter).

[0055] Through a comprehensive analysis of tidal amplitude, maximum tidal velocity, and tidal level change rate, the impact of ocean tides on measurement results can be accurately assessed. The dynamic characteristics of the tide determine whether the error correction in ocean gravity measurement is accurate enough. Revision combined with the prediction model not only ensures the accuracy of ocean gravity measurement, but also improves the adaptability to tidal changes, making the error more flexible and accurate. This process reduces the interference of error sources on measurement results, thereby improving the reliability of ocean measurement data. By fully understanding the impact of tides, the compensation method can dynamically respond to tidal changes, effectively improve the quality of ocean gravity measurement data, and promote progress in the fields of marine scientific research and resource exploration.

[0056] Furthermore, the ocean tide impact prediction model is as follows: ; In the formula, Indicates the ocean tide impact index, represents the tidal amplitude of the predicted ocean area during the measurement constraint period, represents the unit tidal amplitude influence factor, represents the maximum tidal velocity in the predicted ocean area during the measurement constraint period, Indicates the maximum flow velocity factor per unit tide, represents the rate of change of tidal level in the predicted ocean area during the measurement constraint period, Represents the factor affecting the rate of change of unit tidal water level.

[0057] In this embodiment, the ocean tidal impact index is obtained by analyzing the ocean tidal parameters, taking into account the mutual influence relationship between these parameters, for example: the tidal amplitude determines the maximum height change of the water level within the tidal cycle, which is directly related to the flow rate of the water flow. When the tidal amplitude increases, the fluctuation range of the water level expands, and the flow rate of the water flow also increases accordingly, so there is a positive correlation between the maximum tidal flow rate and the tidal amplitude. The tidal water level change rate is closely related to the tidal amplitude and the maximum tidal flow rate, because the tidal water level change rate is affected by the tidal amplitude, and the water level change rate has different rates between high tide and low tide. In summary, the tidal amplitude, the maximum tidal flow rate and the tidal water level change rate are interdependent and interact with each other, and they are mutually causal through physical relationships. Changes in the tidal amplitude and the maximum flow rate usually directly lead to changes in the water level change rate.

[0058] It should be noted that the unit tidal amplitude influence factor, the unit tidal maximum flow velocity influence factor and the unit tidal water level change rate influence factor can be obtained through the database. The unit tidal amplitude influence factor indicates the degree of influence of each unit tidal amplitude on the ocean tide, the unit tidal maximum flow velocity influence factor indicates the degree of influence of each unit tidal maximum flow velocity on the ocean tide, and the unit tidal water level change rate influence factor indicates the degree of influence of each unit tidal water level change rate on the ocean tide.

[0059] It should also be noted that the unit tidal amplitude influence factor can be constructed into a tidal amplitude mapping set by using the tidal amplitude stored in the database and the influence factor corresponding to the tidal amplitude, wherein there is a one-to-one or many-to-one correspondence in the tidal amplitude mapping set, and the unit tidal amplitude influence factor can be obtained by inputting the real-time tidal amplitude into the tidal amplitude mapping set. Other influence factors such as the unit tidal maximum flow velocity influence factor and the unit tidal water level change rate influence factor can be obtained by constructing their corresponding mapping sets in the same way, for example: the unit tidal maximum flow velocity influence factor needs to be obtained in the tidal flow velocity mapping set, and the unit tidal water level change rate influence factor needs to be obtained in the tidal water level change rate mapping set. It should be noted that the tidal flow velocity mapping set and the tidal water level change rate mapping set are constructed in the same way as the tidal amplitude mapping set.

[0060] Furthermore, an ocean gravity correction data set is obtained based on the ocean gravity error revision coefficient and the ocean gravity data set analysis. The specific steps include: extracting the gravity acceleration measured by each ocean gravimeter in the ocean gravity data set, and performing discrete mean processing to obtain a specified gravity acceleration; based on the ocean gravity error revision coefficient, and analyzing the ocean gravity data set to obtain a gravity acceleration correction adjustment value, and combining the specified gravity acceleration to construct a specified gravity acceleration interval for the predicted ocean area; the specified gravity acceleration and the specified gravity acceleration interval for the predicted ocean area are jointly recorded as the ocean gravity correction data set.

[0061] By analyzing the ocean gravity error correction coefficient and the ocean gravity data set to obtain the ocean gravity correction data set, the gravity measurement errors caused by environmental factors, terrain factors, etc. can be effectively eliminated, thereby improving the accuracy and reliability of ocean gravity data.

[0062] It should be noted that the gravity acceleration measured by each ocean gravimeter in the ocean gravity data set is extracted, and discrete mean processing is performed to obtain the specified gravity acceleration, wherein discrete mean processing refers to obtaining the gravity acceleration measured by each ocean gravimeter, removing the maximum and minimum values, and taking the average value to obtain the reference mean of the gravity acceleration measured by each ocean gravimeter, and taking the average value after removing the maximum and minimum values ​​of all the data obtained from each measurement, that is, obtaining the discrete mean.

[0063] like Figure 2 As shown, it is a module schematic diagram of the marine gravity measurement error compensation system provided by the embodiment of the present application. The marine gravity measurement error compensation system provided by the embodiment of the present application includes: a comprehensive and complex analysis module, a constraint judgment module, a revision coefficient evaluation module and a correction execution module; wherein the comprehensive and complex analysis module is used to obtain the terrain parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area, and comprehensively analyze to obtain the comprehensive and complex indicators of the predicted ocean area; the constraint judgment module is used to analyze the comprehensive and complex indicators of the predicted ocean area to obtain the marine gravity measurement constraint parameter set; the revision coefficient evaluation ... The invention relates to a method for performing ocean gravity measurement based on an ocean gravity measurement constraint parameter set, obtaining an ocean gravity data set, analyzing and obtaining an ocean gravity difference index, thereby matching and obtaining an ocean gravity revision coefficient, synchronously obtaining ocean tidal parameters and a preset ocean tidal impact prediction model, and inputting the ocean tidal parameters into the ocean tidal impact prediction model to obtain an ocean tidal impact index, thereby matching and obtaining an ocean tidal revision coefficient; a correction execution module is used to perform analysis based on the ocean tidal revision coefficient and the ocean gravity revision coefficient to obtain an ocean gravity error revision coefficient, and to obtain an ocean gravity correction data set based on the analysis of the ocean gravity error revision coefficient and the ocean gravity data set.

[0064] To summarize, this embodiment obtains the terrain parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area, and conducts comprehensive analysis to obtain the comprehensive complex indicators of the predicted ocean area, thereby obtaining a set of ocean gravity measurement constraint parameters, thereby realizing accurate analysis of ocean gravity, and effectively solving the problem of inaccurate ocean gravity measurement caused by complex seabed topography and tidal fluctuation interference in the prior art.

[0065] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0066] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0067] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0068] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0069] Although the preferred embodiments of the present invention have been described, those skilled in the art may make other changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the present invention.

[0070] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include these modifications and variations.

Claims

1. A method for compensating marine gravity measurement errors, characterized in that: The following steps are involved: Obtain the topographic parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area, and conduct comprehensive analysis to obtain the comprehensive and complex indicators of the predicted ocean area; Based on the comprehensive and complex index analysis of the predicted ocean area, the constraint parameter set of ocean gravity measurement is obtained; Perform ocean gravity measurement based on the ocean gravity measurement constraint parameter set, obtain an ocean gravity data set, analyze and obtain an ocean gravity difference index, thereby matching and obtaining an ocean gravity revision coefficient, synchronously obtain ocean tidal parameters and a preset ocean tidal impact prediction model, and input the ocean tidal parameters into the ocean tidal impact prediction model to obtain an ocean tidal impact index, thereby matching and obtaining an ocean tidal revision coefficient; Based on the analysis of the ocean tide correction coefficient and the ocean gravity correction coefficient, the ocean gravity error correction coefficient is obtained, and based on the analysis of the ocean gravity error correction coefficient and the ocean gravity data set, the ocean gravity correction data set is obtained.

2. The method for compensating ocean gravity measurement errors according to claim 1, characterized in that: The steps of obtaining the topographic parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area and comprehensively analyzing to obtain the comprehensive complex index of the predicted ocean area include: Acquiring topographic parameters of the predicted ocean area, wherein the topographic parameters of the predicted ocean area include an average seafloor slope, an extreme difference in seafloor topography height, an average thickness of seafloor rock formations, and an average hardness of seafloor rock formations; Obtaining water flow fluctuation parameters of the predicted ocean area, wherein the water flow fluctuation parameters of the predicted ocean area include the maximum water flow velocity, the water flow velocity range and the turbulence intensity at each screening point within a preset time period; Obtaining a preset terrain reference set and a water flow fluctuation reference set in the database, and comparing them with the terrain parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area, respectively, to obtain a comprehensive complex index of the predicted ocean area; The comprehensive complexity index of the predicted ocean area is used to characterize the complexity of the predicted ocean area; The terrain reference set includes a slope reference value, a terrain height extreme difference reference value, a seabed rock layer thickness reference value, and a seabed rock layer hardness reference value; The water flow fluctuation reference set includes a water flow velocity reference value, a water flow velocity extreme difference reference value and a turbulence intensity reference value.

3. The method for compensating ocean gravity measurement errors according to claim 2, characterized in that: The comprehensive and complex indicators of the predicted ocean area include: Based on the terrain parameters and terrain reference set of the predicted ocean area, the terrain complexity index of the predicted ocean area is obtained; Based on the analysis of the water flow fluctuation parameters and the water flow fluctuation reference set in the predicted ocean area, the water flow fluctuation complexity index in the predicted ocean area is obtained; Based on the terrain complexity index and the water flow fluctuation complexity index of the predicted ocean area, the comprehensive complexity index of the predicted ocean area is obtained; The terrain complexity index of the predicted ocean area is used to characterize the terrain complexity of the predicted ocean area; The water flow fluctuation complexity index of the predicted ocean area is used to characterize the complexity of the water flow fluctuation in the predicted ocean area.

4. The method for compensating ocean gravity measurement errors according to claim 1, characterized in that: The comprehensive and complex index analysis based on the predicted measurement of the ocean area is used to obtain the ocean gravity measurement constraint parameter set, and the specific steps include: Obtain each preset comprehensive and complex index interval in the database and the ocean gravity measurement constraint reference parameter set corresponding to each comprehensive and complex index interval, and compare them with the comprehensive and complex index of the predicted ocean area. If the comprehensive and complex index of the predicted ocean area is within a preset comprehensive and complex index interval, obtain the ocean gravity measurement constraint reference parameter set corresponding to the comprehensive and complex index interval as the ocean gravity measurement constraint parameter set; The ocean gravity measurement constraint parameter set includes the number of ocean gravity meters deployed and the measurement constraint period.

5. The method for compensating ocean gravity measurement errors according to claim 1, characterized in that: The method of performing ocean gravity measurement based on the ocean gravity measurement constraint parameter set, obtaining an ocean gravity data set, and analyzing to obtain an ocean gravity difference index specifically includes the following steps: Deploy ocean gravimeters based on the ocean gravity measurement constraint parameter set, thereby performing ocean gravity measurement and obtaining an ocean gravity data set, wherein the ocean gravity data set includes gravity acceleration measured by each ocean gravimeter at each time and a center distribution distance and a vertical distribution height of each gravity measuring instrument; Obtain the ocean gravity fitting data set and analyze it with the ocean gravity data set to obtain the ocean gravity difference index; The ocean gravity fitting data set includes a gravity acceleration fitting value, a center distribution distance fitting value, and a vertical distribution height fitting value; The ocean gravity difference index is used to characterize the degree of measurement difference of the ocean gravimeter.

6. The method for compensating ocean gravity measurement errors according to claim 1, characterized in that: The matching obtains the ocean gravity revision coefficient, and the specific steps include: Obtaining each ocean gravity difference index interval preset in the database and the compensation reference factor corresponding to each ocean gravity difference index interval; Based on the ocean gravity difference index and compared with each ocean gravity difference index interval, if the ocean gravity difference index is within a preset ocean gravity difference index interval, the compensation reference factor corresponding to the interval is obtained as the ocean gravity revision coefficient; The ocean gravity revision coefficient is used to revise the ocean gravity.

7. The method for compensating ocean gravity measurement errors according to claim 1, characterized in that: The ocean tidal parameters are input into the ocean tidal impact prediction model to obtain the ocean tidal impact index, and the ocean tidal revision coefficient is obtained by matching. The specific steps include: Acquiring ocean tidal parameters, the ocean tidal parameters including tidal amplitude, maximum tidal flow velocity and tidal water level change rate of the predicted ocean area in the measurement constraint period; Input the ocean tide parameters into the ocean tide impact prediction model to obtain the ocean tide impact index; The ocean tide impact index is used to characterize the impact degree of ocean tide; Obtain each ocean tidal impact index interval preset in the database and the ocean tidal reference revision coefficient corresponding to each ocean tidal impact index interval, and compare with the ocean tidal impact index; if the ocean tidal impact index is within a certain ocean tidal impact index interval, obtain the ocean tidal reference revision coefficient corresponding to the ocean tidal impact index interval as the ocean tidal revision coefficient, and the ocean tidal revision coefficient is used to reduce the influence of ocean tides on the measurement of ocean gravity acceleration.

8. The method for compensating ocean gravity measurement errors according to claim 1, characterized in that: The ocean tide impact prediction model is specifically: ； In the formula, Indicates the ocean tide impact index, represents the tidal amplitude of the predicted ocean area during the measurement constraint period, represents the unit tidal amplitude influence factor, represents the maximum tidal velocity in the predicted ocean area during the measurement constraint period, Indicates the maximum flow velocity factor per unit tide, represents the rate of change of tidal level in the predicted ocean area during the measurement constraint period, Represents the factor affecting the rate of change of unit tidal water level.

9. The method for compensating ocean gravity measurement errors according to claim 1, characterized in that: The ocean gravity correction dataset is obtained based on the ocean gravity error correction coefficient and the ocean gravity dataset analysis, and the specific steps include: Extract the gravity acceleration measured by each ocean gravity meter in the ocean gravity data set, and perform discrete mean processing to obtain the specified gravity acceleration; Based on the ocean gravity error correction coefficient, the gravity acceleration correction adjustment value is obtained by analyzing the ocean gravity data set, and combined with the specified gravity acceleration to construct the specified gravity acceleration interval of the predicted ocean area; The specified gravity acceleration and the specified gravity acceleration interval of the predicted ocean area are jointly recorded as the ocean gravity correction data set.

10. A system using the method for compensating ocean gravity measurement errors as claimed in any one of claims 1 to 9, characterized in that: include: Comprehensive complex analysis module, constraint determination module, revision coefficient evaluation module and revision execution module; The comprehensive and complex analysis module is used to obtain the topographic parameters of the predicted ocean area and the water flow fluctuation parameters of the predicted ocean area, and to obtain the comprehensive and complex indicators of the predicted ocean area through comprehensive analysis; The constraint determination module is used to obtain a set of marine gravity measurement constraint parameters based on comprehensive and complex index analysis of the predicted ocean area; The revision coefficient evaluation module is used to perform ocean gravity measurement based on the ocean gravity measurement constraint parameter set, obtain an ocean gravity data set, analyze and obtain an ocean gravity difference index, thereby matching and obtaining an ocean gravity revision coefficient, synchronously obtain ocean tide parameters and a preset ocean tide impact prediction model, and input the ocean tide parameters into the ocean tide impact prediction model to obtain an ocean tide impact index, thereby matching and obtaining an ocean tide revision coefficient; The correction execution module is used to analyze based on the ocean tide correction coefficient and the ocean gravity correction coefficient to obtain the ocean gravity error correction coefficient, and to analyze based on the ocean gravity error correction coefficient and the ocean gravity data set to obtain the ocean gravity correction data set.

Images