Multi-source data-based optimal depth reference surface geodetic height model construction method
Technical Field
The invention relates to the technical field of ocean measurement, in particular to a method for constructing an optimal depth datum plane geodetic height model based on multi-source data.
Background
Marine mapping is a very important fundamental task, and all development and construction work in the ocean is independent of the support of the various fundamental data provided by marine mapping. Thus, marine mapping is receiving attention from governments around the world as a basis and guideline for marine development. Sea duct measurement is a precursor to marine surveying, which is an internationally common special technique whose main task is to measure the depth of the relevant sea area to ensure the safety of vessel voyage. At present, the development of ocean development and national defense construction cannot be satisfied by acquiring a small amount of information through the traditional measurement means, and along with the development of various measurement platforms such as special measurement ships, airplanes, satellites and the like, the information such as complete ocean geometric fields, physical fields and the like can be obtained by acquiring data through sensors. Traditional paper sea chart is evolving into dynamic real-time sea service with advanced technologies such as digital sea chart, database and information system as cores. The references of marine mapping can be divided into planar and vertical references of spatial position information and gravitational and magnetic references of geophysical survey information, where the vertical references of marine mapping are referred to as depth reference planes, which are depth data start planes. Under the influence of tide, ocean current, wind and wave and other factors, the sea surface can vibrate, the height of the sea surface is fluctuant along with the fluctuation of the tide and the tide, and sometimes the highest drop can reach tens of meters. These differences vary with changes in tidal water duration and tidal interval, particularly in certain sea areas. In order to correct the deviation caused by tides in the measured water depth data, a stable initial surface needs to be determined, and the instantaneous water depth at a certain point observed at different moments is calculated to the initial surface, namely the depth reference surface. Whereas the depth reference plane defines criteria that are close to but not lower than the lowest tide level that may actually occur, i.e. they are all defined as being close to but not lower than the lowest tide level that may actually occur.
China has wide ocean area, and different depth references are adopted in different periods. The data obtained by the existing observation technology usually takes a reference ellipsoid as a starting surface, and in order to ensure the safety of offshore operation, depth data needs to be converted into a reference standard. When measuring ocean depths, an average sea level and a depth datum plane are usually used as reference datum, and in practical situations, the reference datum used in measurement tasks of different times and different units is difficult to be consistent. Because of tidal fluctuation, coastal areas at sea-land junctions are difficult to accurately measure, and the local average sea level is mainly adopted as a reference standard; in some ports and channels, a depth datum plane is generally adopted as a reference datum for ensuring navigation safety. For decades, although China has extensively explored and measured the surrounding sea area, due to the long-term lack of uniform high-precision ocean measurement references, the fusion of data exchange acquired by different departments and different times, the splicing of related maps and land elevation data and ocean depth data is difficult to realize. This greatly reduces the efficiency of the measurement and the value of the data used, wasting a lot of labor.
If a unified and continuous vertical datum is determined, conversion and output of sounding data on different reference datum can be easily achieved, and ocean and land geographic information can be more easily spliced and combined. Therefore, the method for constructing the ground high model of the optimal depth datum plane based on the multi-source data is provided.
Disclosure of Invention
The invention aims at overcoming the defects of the prior art, and provides an optimal depth datum plane geodetic high model construction method based on multi-source data, so as to solve the problems of the background art.
In order to achieve the above purpose, the present invention provides the following technical solutions: the method for constructing the ground high model of the optimal depth datum plane based on the multi-source data comprises the following specific steps:
s1: constructing an average sea surface height model;
extracting data points of average sea surface high lattice points of the south sea area based on a CNES_CLS2015 model, and establishing a depth datum plane geodetic height model of the south sea area based on a CGCS2000 reference ellipsoid;
s2: constructing a south sea optimal tide model;
carrying out harmonic analysis on tide level data of a long-term tide station, evaluating the precision of four tide models, screening out tide branches with highest precision, and combining the optimal tide branches to realize optimization of the tide models;
s3: constructing a depth reference surface L value model of the south China sea area based on the refined optimal model, and correcting the depth reference L value model by utilizing the L value of the long-term tide station;
s4: and constructing a depth reference surface geodetic model of the research area by adopting a model difference method based on the average sea surface height model and the depth reference surface model.
As a preferred embodiment of the present invention, in S1 and S2, the resolution of the average sea surface high model and the resolution of the depth reference plane L value model need to be the same.
As a preferable technical scheme of the invention, the optimization method for the tide model in the step S2 comprises the following steps: and (3) carrying out precision verification on the plurality of tide models according to the tide checking station, comparing and analyzing the tide with the highest precision of different tide models, and then combining the tide with the highest precision into a new tide model.
As a preferable technical scheme of the invention, the optimization method for the tide model in the step S2 is based on the difference ratio relation optimization method as follows:
the variation of the difference ratio relation between different adjustment and constant of two adjacent tide stations has stronger correlation and consistency, namely the difference ratio relation does not change along with the position variation, and the adjustment constant of the tide stations in the middle and short period is corrected by the difference ratio relation of different tide dividing adjustment constants of the tide stations in the long period; according to the method, according to the precision of the tide model, a tide with better precision is selected as a basic tide; the method comprises the following specific steps:
the method comprises the steps of selecting a tide with highest precision as a basic tide, and setting the amplitude and the retarded angle relation of the basic tide in a tide station to be expressed by the following formula:
wherein h' represents the amplitude ratio of the remaining moisture to the basic moisture; g' represents the delay angle difference between the rest moisture and the basic moisture; the superscript i indicates the ith moisture division to be refined; the superscript S indicates the selected basic moisture; o represents a tide station;
the difference ratio relation of the tide stations is applied to lattice points of the tide model nearby, and the calculation formula is as follows:
in the middle of
Newly calculating the amplitude of lattice points for the tidal model; />
Newly calculating the delay angle of the lattice point for the tide model;
the amplitude and the delay angle of the basic tide of the original grid point of the tide model are respectively. />
As a preferable technical scheme of the invention, the model difference method used for constructing the depth datum plane geodetic high model in the S4 is as follows:
firstly, analyzing the interrelationship of a reference ellipsoid, an average sea level and a sea chart depth standard used in the model construction process, wherein H is tidal height and H s Distance h from depth datum plane to reference ellipsoid g The distance from the average sea level to the reference ellipsoid is given by L, which is the rising reading reference value;
if another h mss Is the average sea level height of a point, thenThe depth reference surface geodetic height calculation formula is as follows: h is a s =h mss -L。
As a preferred technical solution of the present invention, the constructing the average sea surface elevation model in S1 includes the following operations:
s11: a depth reference surface geodetic height model construction flow;
s12: an average picture high contour map;
s13: the principle of tidal model refinement;
s14: a depth datum plane L value contour map;
s15: relationships between various references;
s16: model difference principle;
s17: depth reference surface geodetic model contour map.
The beneficial effects of the invention are as follows: the reference ellipsoid is used as a final reference standard, a depth reference surface model established based on an average sea surface high model and a high-precision tide model established by multi-source data, namely satellite height measurement data, is used for establishing an optimal depth reference surface geodetic high model of a research area by adopting a model difference method, and a conversion relation of various standards is determined, so that a foundation is laid for unifying sea-land vertical standards.
Drawings
FIG. 1 is a flow chart of the depth reference surface geodetic high model construction of the present invention;
FIG. 2 is an average picture high contour plot of the present invention;
FIG. 3 is a schematic diagram of a tidal model refinement of the present invention;
FIG. 4 is a depth reference plane L value contour map of the present invention;
FIG. 5 is a graph of the relationship between various references of the present invention;
FIG. 6 is a schematic diagram of the model difference method of the present invention;
FIG. 7 is a contour map of a depth referencing geodetic model of the present invention;
Detailed Description
The preferred embodiments of the present invention will be described in detail below with reference to the attached drawings so that the advantages and features of the present invention can be more easily understood by those skilled in the art, thereby making clear and defining the scope of the present invention.
Examples: referring to fig. 1-7, the present invention provides a technical solution: the embodiment takes a south sea area as a research object, combines satellite altitude data, tide model data and average sea surface altitude model data, reduces an average sea surface and a depth reference surface to a reference ellipsoid by establishing a depth reference surface geodetic altitude model of a research area, and provides a basis for sea Liu Chuizhi reference conversion, and the construction method specifically comprises the following steps:
s1: constructing an average sea surface height model; extracting data points of average sea surface high lattice points of the south sea area based on a CNES_CLS2015 model, and establishing a depth datum plane geodetic height model of the south sea area based on a CGCS2000 reference ellipsoid; the specific operation is as follows:
s11: a depth reference surface geodetic height model construction flow;
s12: an average picture high contour map;
s13: the principle of tidal model refinement;
s14: a depth datum plane L value contour map;
s15: relationships between various references;
s16: model difference principle;
s17: a depth datum surface geodetic model contour map;
s2: constructing a south sea optimal tide model; carrying out harmonic analysis on tide level data of a long-term tide station, evaluating the precision of four tide models, screening out tide branches with highest precision, and combining the optimal tide branches to realize optimization of the tide models;
s3: constructing a depth reference surface L value model of the south China sea area based on the refined optimal model, and correcting the depth reference L value model by utilizing the L value of the long-term tide station;
s4: and constructing a depth reference surface geodetic model of the research area by adopting a model difference method based on the average sea surface height model and the depth reference surface model.
In S1 and S3, the resolution of the average sea surface elevation model and the depth reference plane L value model need to be the same.
In S2, the tide model is optimized, and the optimal tide combination optimization method is as follows:
the method mainly comprises the steps of verifying the precision of a plurality of tide models according to a tide station, comparing and analyzing tide branches with highest precision of different tide models, and then combining the tide branches with highest precision into a new tide model; the method carries out accuracy verification on three different tide models TPXO_7.2, NAO99b, TPXO_yellow and FES2014b, and the result shows that the accuracy of O1 moisture division of the TPXO_7.2 model is highest; the M2 moisture separation precision of the NAO99b model is highest; the moisture separation precision of Q1, P1, K1 and K2 of the TPXO_yellow model is highest; the FES2014b tidal model has the highest accuracy of the moisture separation of N2 and S2.
In S2, the tidal model is optimized, based on the difference-to-ratio relationship optimization method as follows:
the variation of the difference ratio relation between different modulations and constants of two adjacent tide stations has stronger correlation and consistency, namely the difference ratio relation does not change along with the position variation. The harmonic constants of the middle-short-term tide stations can be corrected by the difference ratio relation of different tide-dividing harmonic constants of the long-term tide stations. Based on the idea, the grid points of the tidal model in a certain range are corrected by utilizing long-term tide checking data around the south China sea and satellite altitude data so as to achieve the purpose of model refinement; according to the method, according to the precision of the tide model, a tide with better precision is selected as a basic tide; the method comprises the following specific steps:
the method comprises the steps of selecting a tide with highest precision as a basic tide, and setting the amplitude and the retarded angle relation of the basic tide in a tide station to be expressed by the following formula:
wherein h' represents the amplitude ratio of the remaining moisture to the basic moisture; g' represents the delay angle difference between the rest moisture and the basic moisture; the superscript i indicates the ith moisture division to be refined; the superscript S indicates the selected basic moisture; o represents the tide station.
The difference ratio relation of the tide stations is applied to lattice points of the tide model nearby, and the calculation formula is as follows:
in the middle of
Newly calculating the amplitude of lattice points for the tidal model; />
Newly calculating the delay angle of the lattice point for the tide model;
the amplitude and the delay angle of the basic tide of the original grid point of the tide model are respectively.
The model difference method principle used for constructing the depth datum surface geodetic height model in S4 is as follows:
firstly, analyzing the interrelationship of reference ellipsoids, average sea level and sea chart depth references used in the model construction process, wherein H is tide height and H as shown in figure 5 s Distance h from depth datum plane to reference ellipsoid g The distance from the mean sea level to the reference ellipsoid, L, is the read-up reference value there (relative to the mean sea level);
if another h mss For an average sea surface height (relative to a reference ellipsoid) at a point, the depth reference surface geodetic height at that point is calculated as follows: h is a s =h mss -L。
The reference ellipsoid is used as a final reference standard, a depth reference surface model established based on an average sea surface high model and a high-precision tide model established by multi-source data, namely satellite height measurement data, is used for establishing an optimal depth reference surface geodetic high model of a research area by adopting a model difference method, and a conversion relation of various standards is determined, so that a foundation is laid for unifying sea-land vertical standards.
The foregoing examples illustrate only a few embodiments of the invention, which are described in detail and are not to be construed as limiting the scope of the invention. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the invention, which are all within the scope of the invention.