CN107883931A - Tide rectification method and device, storage medium based on structure instantaneous water level model - Google Patents
Tide rectification method and device, storage medium based on structure instantaneous water level model Download PDFInfo
- Publication number
- CN107883931A CN107883931A CN201711009326.0A CN201711009326A CN107883931A CN 107883931 A CN107883931 A CN 107883931A CN 201711009326 A CN201711009326 A CN 201711009326A CN 107883931 A CN107883931 A CN 107883931A
- Authority
- CN
- China
- Prior art keywords
- tidal
- water level
- instantaneous
- tide
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
- G01C13/008—Surveying specially adapted to open water, e.g. sea, lake, river or canal measuring depth of open water
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Hydrology & Water Resources (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of tide rectification method and device, storage medium based on structure instantaneous water level model, wherein, methods described includes step：Numerical simulation is carried out to the astronomical tidal level surveyed in area first with ocean numerical models, obtain comprising the instantaneous astronomical tide numerical model including the measurement activity duration, then according to survey tidal observation data structure that the tidal station that area is distributed observes it is instantaneous more than water level correction model and instantaneous residual water level correction model, finally according to the instantaneous astronomical tide numerical model, instantaneous remaining water level correction model and instantaneous residual water level correction model structure instantaneous water level model, the accurate correction of area's tidal level is surveyed by the instantaneous water level model realization of structure, obtain accurate water-depth measurement data, so as to solve prior art when carrying out water-depth measurement, tidal station distribution because surveying area is less or surveys the problem of tidal current field complexity in area causes tidal level not corrected accurately.
Description
Technical field
The present invention relates to water-depth measurement field, more particularly to tide rectification method and dress based on structure instantaneous water level model
Put, storage medium.
Background technology
Now, water-depth measurement has become the central task of hydrographic survey, while water-depth measurement is also to make sea chart, obtain
The most effective approach of sea-floor relief change, the means of water-depth measurement have a lot, including airborne laser water-depth measurement and satellite remote sensing
Depth extraction etc., however, for coastal waters littoral sea water-depth measurement, still walked at present with ship based on aerial survey amount, no matter pass through
Which kind of means, it is largely high for the Instantaneous Sea Level at measurement moment that achievement is surveyed in water-depth measurement, different big containing at different moments
The small water level as caused by tide influences, and is measured different from ocean, in coastal waters, particularly Coastal sea area, because the depth of water is shallower, and tide
Nighttide turns into the principal element for restricting bathymetric surveying precision, and due to shallow water effect, the tidal range of littoral sea also much larger than
Ocean marine site, therefore, in high-precision coastal waters coastwise sounding operation, tide rectification turns into an essential operation
Content.
Instantaneous Sea Level influenceed by tidal effect and has behavioral characteristics, in order to be stablized during hydrographic survey
Underwater topography depth of water numerical value, it is necessary to eliminate the influence of tide, i.e., by the instantaneous depth measurement value naturalization in sea area to one it is stable known to
The process that depth datum gets on is referred to as tide rectification, because the precision of tide rectification is directly affected finally into the precision of figure, because
Each coastal state in this world attaches great importance to tide rectification, because the sea area situation and technical requirements of various countries are different, because
The method of this used tide rectification is also different, but its general principle is all based on measuring the water level sight of sea area tidal station
Measurement information or tide prediction value carry out temporal interpolation and most of regional space interpolation, these interpolation algorithms is several using space
What interpolation, however, when surveying tidal station of area and being distributed insufficient or when surveying area's change of tidal current field complexity, survey the spatial variations of area's tidal level
The result that may be come out with geometry interpolation is not consistent, or even larger error occurs.
Therefore, prior art has yet to be improved and developed.
The content of the invention
In view of above-mentioned the deficiencies in the prior art, it is an object of the invention to provide the tidal level based on structure instantaneous water level model
Correcting method and device, storage medium, it is intended to less due to surveying the distribution of tidal station of area when solving existing method progress water-depth measurement
Or survey area's change of tidal current field complexity so that the technical problem that tidal level can not be corrected accurately.
Technical scheme is as follows：
A kind of tide rectification method based on structure instantaneous water level model, wherein, including step：
Numerical simulation is carried out to the astronomical tidal level surveyed in area using ocean numerical models, obtains existing comprising the measurement activity duration
Interior instantaneous astronomical tide numerical model；
According to surveying the instantaneous remaining water level correction model of tidal observation data structure that the tidal station that area is distributed observes and instantaneously remaining water
Position correction model；
According to the instantaneous astronomical tide numerical model, it is instantaneous more than water level correction model and instantaneous residual water level correction model structure
Instantaneous water level model is built, the accurate correction of area's tidal level is surveyed by the instantaneous water level model realization.
Described tide rectification method, wherein, the step is entered using ocean numerical models to surveying the astronomical tidal level in area
Row numerical simulation, obtain, comprising the instantaneous astronomical tide numerical model including the measurement activity duration, specifically including：
Determine simulated domain and extract simulated domain water front data, make the grid that becomes more meticulous of simulated domain；
Measure and calculate and survey area's bathymetric data, the bathymetric data is loaded into the grid that becomes more meticulous, obtains ocean
Domain condition needed for numerical model；
Open on boundary point one year tide prediction number as simulation duration is provided using the harmonic constant forecast of partial tide
According to foundation selects the demand of pattern, and data are handled, obtain the open boundaryconditions of corresponding format；
Select and the primary condition before the numerical model startup of ocean is set, high accuracy is carried out to the astronomical tide process surveyed in area
Numerical simulation, obtain instantaneous astronomical tide numerical model.
Described tide rectification method, wherein, the step is according to the tidal observation data structure for surveying the tidal station observation that area is distributed
Instantaneous remaining water level correction model is built to specifically include：
Based on the shot and long term tidal station tidal observation data for surveying area's distribution, using harmonic analysis principle, mediation meter is carried out to tidal level
Calculate, so as to obtain partial tide harmonic constant；
Harmonic constant is modified, obtains the higher harmonic constant of stability, and carry out based on stable harmonic constant
Tide prediction；
Forecast tidal level value is subtracted using tidal station actual measurement tidal level, obtains the remaining water level value of the tidal station, then test based on all
The remaining water level value at tide station, using interpolation algorithm, obtain surveying remaining water level model in area.
Described tide rectification method, wherein, the step is according to the tidal observation data structure for surveying the tidal station observation that area is distributed
Building instantaneous residual water level correction model also includes：
The remainder error in the simulation error and remaining water level extraction process of ocean numerical models in itself is extracted, is calculated using interpolation
Method, build instantaneous residual water level correction model.
Described tide rectification method, wherein, the instantaneous remaining water level correction model expression formula is R (t)=h (t)-MSL-
T (t)-Δ, wherein, R (t) represents Yu Shuiwei sequences, and h (t) represents measured water level sequence, and MSL represents mean sea level, and T (t) is represented
Astronomical tidal level, Δ represent observation error.
Described tide rectification method, wherein, the instantaneous residual water level correction model expression formula is Ei(t)=hi(t)-
Ti(t)-Ri(t), wherein, Ei(t) the remaining tide station modified value corresponding to tidal station i is represented；hi(t) tidal station i actual measurement tides are represented
Place value, Ti(t) tidal station i simulation astronomical tide place value, R are representedi(t) tidal station i remaining water level value is represented.
Described tide rectification method, wherein, the instantaneous water level model expression is Tj(t)=hj(t)+Rj(t)+Cj
(t), wherein, Tj(t) it is tide rectification value at depth measurement point j required by t, hj(t) be t at depth measurement point j astronomical tidal level
The analogue value, Rj(t) be t at depth measurement point j remaining water level correction value, Cj(t) be t at depth measurement point j residual water level amendment
Value.
A kind of storage medium, wherein, a plurality of instruction is stored with, the instruction is suitable to be loaded by processor and performed above-mentioned
The step of described in meaning one based on the tide rectification method for building instantaneous water level model.
A kind of tide rectification device based on structure instantaneous water level model, wherein, including processor, it is adapted for carrying out each finger
Order；And storage device, suitable for storing a plurality of instruction, the instruction is suitable to be loaded by processor and perform above-mentioned any one institute
The step of stating the tide rectification method based on structure instantaneous water level model.
Beneficial effect：The invention discloses it is a kind of based on structure instantaneous water level model tide rectification method, first with
Ocean numerical models carry out numerical simulation to the astronomical tidal level surveyed in area, obtain comprising the instantaneous astronomy including the measurement activity duration
Damp numerical model, instantaneous remaining water level correction model and instantaneously then is built according to surveying the tidal observation data that the tidal station that area is distributed observes
Residual water level correction model, finally according to the instantaneous astronomical tide numerical model, instantaneous remaining water level correction model and instantaneous residue
Water level correction model builds instantaneous water level model.The present invention can be realized effectively by the instantaneous water level model of structure surveys area's tidal level
Accurate correction, so as to solve prior art when carrying out water-depth measurement, because survey area tidal station distribution it is less or survey area
Tidal current field complexity the problem of causing tidal level not corrected accurately.
Brief description of the drawings
Fig. 1 is a kind of flow chart of the tide rectification method preferred embodiment based on structure instantaneous water level model of the present invention.
Fig. 2 is a kind of structural frames of the tide rectification device preferred embodiment based on structure instantaneous water level model of the present invention
Figure.
Embodiment
The present invention provides a kind of tide rectification method and device, storage medium based on structure instantaneous water level model, to make
The purpose of the present invention, technical scheme and effect are clearer, clear and definite, and the present invention is described in more detail below.It should be appreciated that
The specific embodiments described herein are merely illustrative of the present invention, is not intended to limit the present invention.
Referring to Fig. 1, Fig. 1 is a kind of tide rectification method preferred embodiment based on structure instantaneous water level model of the present invention
Flow chart, as illustrated, it includes step：
S10, using ocean numerical models to survey area in astronomical tidal level carry out numerical simulation, obtain comprising measurement operation when
Between including instantaneous astronomical tide numerical model；
S20, build instantaneous remaining water level correction model and instantaneously remain according to surveying the tidal observation data that the tidal station that area is distributed observes
Remaining water level correction model；
S30, according to the instantaneous astronomical tide numerical model, it is instantaneous more than water level correction model and instantaneous residual water level amendment mould
Type builds instantaneous water level model, and the accurate correction of area's tidal level is surveyed by the instantaneous water level model realization.
Specifically, prior art carry out water-depth measurement when, when survey tidal station of area distribution it is insufficient when or survey area's tidal current field
When changing complicated, the result that surveying the spatial variations of area's tidal level may come out with geometry interpolation is not consistent, or even occurs larger
Error, tidal level can not be corrected accurately, cause the inaccuracy of water-depth measurement result.
To solve the above problems, the invention discloses it is a kind of based on structure instantaneous water level model tide rectification method, its
In, methods described includes step：Instantaneous astronomical tide model is built using ocean numerical models, is seen using the tidal station for surveying area's distribution
The instantaneous remaining water level correction model of tidal observation data structure of survey and instantaneous residual water level correction model, comprehensive instantaneous astronomical tide model,
Instantaneous remaining water level correction model and instantaneous residual water level correction model structure instantaneous water level model, and it is used for tide rectification.This hair
Bright when solving existing method and carrying out water-depth measurement, the tidal station distribution because surveying area is less or to survey the tidal current field in area complicated, makes
Tidal level can not be corrected accurately, the problem of so as to cause water-depth measurement result inaccuracy.
In the step S10, numerical simulation is carried out to the astronomical tidal level surveyed in area using ocean numerical models, wrapped
Instantaneous astronomical tide numerical model including the activity duration containing measurement.
Specifically, astronomical tide numerical simulation is the height carried out by ocean numerical models to the astronomical tide process surveyed in area
Accurate values are simulated, and simulation partial tide number is M2、S2、N2、K2、K1、O1、P1、Q1And M4Nine most important partial tides, specific method is such as
Under：
From ocean numerical models：Conventional numerical model has FVCOM, POM, ROM and MIKE21 etc..
Determine simulated domain and mesh generation：Simulated domain, which will cover, surveys area, and mould is visually extracted using satellite remote-sensing image
Intend region water front data, utilize the grid that becomes more meticulous of SMS software development simulated domains.The process of mesh generation in zoning
In, the angle quality control of grid has carried out local cypher, water front average mark resolution between 30 ° -120 °, to the water front of complexity
About 0.002 °, about 0.01 ° of -0.04 ° of gradual change of border average mark resolution is opened, centre has island part to use local cypher.
Sea-floor relief inverting：Prepare to survey area's bathymetric data, the quality of bathymetric data is very big on the influence of astronomical tide simulation precision.
All Charted depth datas in the survey area published from navigation support portion of naval, the bathymetric data is to be based on theoretical depth benchmark
Water-filling is entered in face, it is necessary to the local theoretical depth reference plane and the interpolation L of mean sea level that are provided according to survey head of district phase tidal station
The Reference Transforming of deep data.Most bathymetric data Reference Transforming is mean sea level at last.Using SMS softwares, by bathymetric data plus
It is downloaded in the grid made, finally gives astronomical tide numerical simulation required domain (domain) condition.
Open boundaryconditions：Open boundaryconditions are the most important parts of numerical simulation.Utilize M in OTPS programs2、S2、N2、K2、
K1、O1、P1、Q1And M4The harmonic constant forecast of nine partial tides provides the same tide prediction in open on boundary point one year of simulation duration
Data, the time, which includes, surveys area's activity duration, and time interval is 1 hour, and tidal plane of reference uses mean sea level, then, according to
From numerical model demand, border will be opened and make corresponding format requirement.
Ocean numerical models start：Pattern also needs to set primary condition before activation, and pattern uses cold start-up, and temperature is set
10 DEG C (annual surface temperatures) are set to, salinity is arranged to 35 ‰, using Smagorinsky turbulent closure models as horizontal mixed
Conjunction scheme, coefficient of friction are arranged to 0.001, and horizontal mixed coefficint is 0.1, and the modeling time sets the time with open boundaryconditions
Consistent, pattern exports once result per hour, and instantaneous astronomical tide numerical model is built according to the result data of output.
Further, in the present embodiment, the step S20, the tidal observation number observed according to the tidal station for surveying area's distribution
According to the instantaneous remaining water level correction model of structure and instantaneous residual water level correction model.
Instantaneous remaining water level correction model is based on the shot and long term tidal station tidal observation data for surveying area's distribution, utilizes harmonic analysis original
Reason, carries out mediation calculating, so as to obtain partial tide harmonic constant to tidal level；Related algorithm is recycled to be modified harmonic constant,
Obtain the higher harmonic constant of stability；Tide prediction is carried out based on stable harmonic constant again, so as to obtain forecasting tidal level；Most
Afterwards, forecast tidal level value is subtracted using tidal station actual measurement tidal level, so as to obtain the remaining water level value of the tidal station；All tidal observations are based on again
The remaining water level value stood, using the interpolation algorithm of adaptation, it can obtain surveying remaining water level model in area.
Specifically, remaining water level refers to due to being caused sea by random meteorologic factor of the short cycle such as wind, air pressure, precipitation
Disturbance, causes short-term water level abnormality phenomenon occur, its mathematic(al) representation can be classified as：
R (t)=h (t)-MSL-T (t)-Δ formula (1),
Wherein, R (t) represents Yu Shuiwei sequences, and h (t) represents measured water level sequence, and MSL represents mean sea level, and T (t) represents astronomical tide
Position, Δ represent observation error (negligible).MSL and T (t) is obtained by the harmonic analysis to measured water level.As special with signal
The remaining water level part of sign, it is stronger due to remaining water level in coastal waters and bank although not can determine that its clear and definite function model
Spatial coherence, it has the rule in space scale than tide with bigger gradual change, can be achieved when in a suitable approach remaining
The reconstruct of water level spatial field.
1), partial tide harmonic constant computational methods：Tidal level is mainly acted on lower produced by the power to lead tide of the celestial bodies such as the moon, the sun
Tide can by multiple cosine terms stacking pattern express, wherein each cosine term represents a partial tide, its calculation formula is：
In formula：a0Mean sea level is represented, m represents partial tide number, aj=fjHjcos[gj-(v0+u)j],bj=fjHjsin[gj-
(v0+u)j], fjRepresent node factor, σjRepresent partial tide angular speed, (v0+u)jPartial tide first phase is represented, t represents instantaneous time, Hj、gj
As partial tide harmonic constant.Harmonic analysis is exactly to go to approach actual tidal level ζ (t) by algorithm.During using choosing corresponding observation
The intermediate time of section is as the timeorigin calculated.By the principle of least square, it is necessary to make：
When D for minimum, with this determination system measure aj, bj.If take the data of 369 days to be used as to analyze, T=8856 hours in formula.
(2) formula is substituted into (3) formula to obtain：
Continuous variable is quadratured, substituted using the observation of discretization, wherein t=k Δs t=k (generally takes 1 hour),
Then K=-N ,-N+1 ..., 0 ... N-1, N, if taking 369 days, T=369*24=8856=2N, so N=4428,2N+1
=8857.For convenience's sake, when summing below end to end two respectively take half, this is from the meaning of least square method, not
Lose general.
For convenience, a0As the special partial tide processing that angular speed is zero.When end to end two respectively take half when,
Formula (4) can be write as：
According to the principle of least square method, D is sought to aj, bjPartial derivative, and make its be equal to zero, obtain formula (6)：
M be partial tide sum, Σ represent summation when end to end two respectively take half.Meanwhile because formula (7)：
∑sinσjkcosσjK=∑ cos σjksinσjK=0；
So formula (6) can be write as：
Obtained coefficient is brought into formula (6), two Symmetric Equations can be obtained.There is m+1 to previous equation group
Individual equation, a can be solved0, a1..., am, the latter equation group has m equation, can solve b0, b1..., bm。
After obtaining a and b, our cans calculate R and θ according to following formula：
After obtaining R and θ, just by intermediate date at the time of, calculate f, V0+ u, and then can obtains harmonic constant H, g.
2) tide prediction, is carried out according to harmonic constant：The tune that each partial tide is tried to achieve using the harmonic analysis of tide is described above
With the algorithm of constant.The tide prediction of any time can also be carried out using the partial tide harmonic constant tried to achieve.
Wherein, m is the number of partial tide, and f is the node factor of forecast day partial tide；(V0+ u) be forecast day partial tide Greenwich time
Between initial phase angle；A0It is the height of yearly mean sea level or multi-year mean sea level.
3), instantaneous remaining water level correction model：From formula (1), the extraction accuracy of remaining water level depends on astronomical tidal level
Precision.Wherein astronomical tide tidal level is to forecast to obtain according to the partial tide harmonic constant of tidal station, therefore the stabilization of partial tide harmonic constant
Property determine remaining water level extraction accuracy.In actual water-depth measurement operation, some long-term short-term tidal observations of tidal station are generally only had
Data, or construction party can lay some short-term stations as progress tidal level observation, but short-term tidal observation data are directly reconciled
Often ratio of precision is relatively low for the obtained harmonic constant of analysis.Therefore, using water level extracting method more than following two, can effectively solve
This problem.
(1) there is water level extraction more than the short-term station of stable harmonic constant
According to the research of harmonic constant stability, the tune obtained using tide gauge (progress per hour) harmonic analysis in 1 year
There is fabulous stability with constant.Therefore there is the Survey of tide and other times of certain a year and a day for unified website
One month observational data, the harmonic constant that can be stablized by the tidal level harmonic analysis of 1 year, then utilize stable mediation
Constant carries out Interpretation Method of Area Rainfall to short-term station observation period, and the water level of forecast is finally subtracted using the water level value of elapsed time section
Value, formula are as follows：R (t)=h (t)-MSL-TL(t) formula (14), wherein, h (t) represents short-term actual measurement tide level data, MSL tables
Show the mean sea level of the point, TL(t) represent with the forecast water level value stablized obtained by harmonic constant progress Interpretation Method of Area Rainfall.
(2) without stablizing water level extraction more than the short-term station of harmonic constant
Tidal level observation data are referred to as mid-term observational data within one month, can only ensure the basic cycle for meeting tidal fluctuations, with
Result obtained by observational data harmonic analysis in 1 year is there is obvious deviation, while this deviation is tested in the neighbouring of simultaneous observation
Tide has certain correlation between station.According to this correlation, the harmonic constant obtained by the short-term tidal station that is related to is carried out
Improve.
Wherein, H, g are respectively the amplitude and delay angle for studying partial tide；Second subscript m and L represent the moon and data point in the whole year respectively
Analyse result.And treat and correct tidal station B there is same change expression formula：
It is assumed that at two stations, the result of the whole month data harmonic analysis with for many years or analysis result in the whole year difference it is equal, i.e.,：Δ
HAm=Δ HBm, Δ gAm=Δ gBm；Had according to formula (15) and (16)：
Wherein, A tidal level long-established station, B are tidal level short-term station.
The partial tide harmonic constant obtained by short cycle tide gauge is refined using the algorithm, so as to accurately be tied
Fruit.Then Interpretation Method of Area Rainfall is done using revised result to extract remaining water level, and formula is as follows：
R (t)=h (t)-MSL-Tm(t) formula (18), wherein, h (t) represents short-term tidal station tide gauge, and MSL is represented
Mean sea level, Tm(t (represents forecast water level obtained by revised harmonic constant.
According to the amendment value sequence of water level extraction algorithm, just available each tidal station more than above two, using instead away from
It is as follows from weighted formula：
Wherein, Rj(t) it is the remaining water level correction value of t at point j, Ri(t) the remaining water that tidal station i t is extracted is represented
Position correction value；J is the locus of depth measurement point, and i is the sequence number of tidal station；λiIt is water level H t tidal station i at depth measurement point j
Weight, p represent exponential quantity.
Further, in the present embodiment, the step S20 also includes：
Remainder error in the simulation error and remaining water level extraction process of extraction pattern in itself, utilize interpolation algorithm, structure
Instantaneous residual water level correction model.
Specifically, the residue that residual water level error is included in the simulation error and remaining water level extraction process of pattern in itself is missed
Difference composition.Residual water level error includes the error with certain spatial coherence, while also there is nonlinear error.
The calculation formula of residual water level is as follows：Ei(t)=hi(t)-Ti(t)-Ri(t) formula (21), wherein, Ei(t) tidal observation is represented
The remaining tide station modified value stood corresponding to i；hi(t) tidal station i actual measurement tidal level values, T are representedi(t) tidal station i simulation day is represented
Literary tidal level value, Ri(t) water level value more than expression tidal station i.
Because, there is linearity error, therefore the whole area that surveys has a kind of error sequence of stabilization in residual water level.To surveying
Each tidal station residual water level value in area carries out relationship analysis, can obtain the relative coefficient between each tidal station.
Wherein, ρ represents the linearly related characteristic value spent closely between variable X and Y.
The time difference method water level interpolation method after more conventional improvement is used herein for the interpolation of residual water level：
Formula (23), wherein Cj(t) at depth measurement point j t residual water level corrected value, μiFor tidal station i coefficient correlation
Weighted value, Sj(t) it is by water level interpolation algorithm between tidal station j and most related tidal station.DABFor A, B two in time difference method
Tidal station, DABFor the distance between A, B, DAj、DBjFor depth measurement point j to A, the distance at the stations of B two, SA(t)、SB(t) it is the station of A, B two
Water level sequence, Δ tAB(t) it is the when difference sequence at the station of A, B two, Δ tAj(t)、ΔtBj(t) it is depth measurement point j to A, the time difference at the stations of B two
Sequence.The time difference at the station of A, B two is mainly had instantaneity by time difference method improved model instead, preferably to reflect remaining water
Position error change.
Comprehensive instantaneous astronomical tide model, instantaneous remaining water level correction model and instantaneous residual water level correction model build instantaneous water
Bit model, it can obtain surveying area's instantaneous water level model：Tj(t)=hj(t)+Rj(t)+Cj(t) formula (25) wherein, Tj(t) it is depth measurement
Tide rectification value at point j required by t, hj(t) be t at depth measurement point j the astronomical tidal level analogue value, Rj(t) it is depth measurement point
The remaining water level correction value of t, C at jj(t) be t at depth measurement point j residual water level correction value.
The instantaneous water level model built by the present invention can effectively realize the accurate correction for surveying area's tidal level, obtain accurately
Water-depth measurement data, so as to solve prior art when carrying out water-depth measurement, the tidal station distribution because surveying area is less or surveys
The problem of tidal current field complexity in area causes tidal level not corrected accurately.
Based on the above method, present invention also offers a kind of storage medium, wherein, it is stored with a plurality of instruction, the instruction
Suitable for loaded by processor and performed it is above-mentioned any one based on structure instantaneous water level model tide rectification method the step of.
Further, the present invention also provides a kind of tide rectification device based on structure instantaneous water level model, wherein, such as scheme
Shown in 2, including processor 10, it is adapted for carrying out each instruction；And storage device 20, suitable for storing a plurality of instruction, the instruction is suitable
In loaded by processor 10 and perform it is above-mentioned any one based on structure instantaneous water level model tide rectification method the step of.
In summary, the invention discloses it is a kind of based on structure instantaneous water level model tide rectification method, first with
Ocean numerical models carry out numerical simulation to the astronomical tidal level surveyed in area, obtain comprising the instantaneous astronomy including the measurement activity duration
Damp numerical model, instantaneous remaining water level correction model and instantaneously then is built according to surveying the tidal observation data that the tidal station that area is distributed observes
Residual water level correction model, finally according to the instantaneous astronomical tide numerical model, instantaneous remaining water level correction model and instantaneous residue
Water level correction model builds instantaneous water level model.The present invention can be realized effectively by the instantaneous water level model of structure surveys area's tidal level
Accurate correction, so as to solve prior art when carrying out water-depth measurement, because survey area tidal station distribution it is less or survey area
Tidal current field complexity the problem of causing tidal level not corrected accurately.
It should be appreciated that the application of the present invention is not limited to above-mentioned citing, for those of ordinary skills, can
To be improved or converted according to the above description, all these modifications and variations should all belong to the guarantor of appended claims of the present invention
Protect scope.
Claims (9)
- A kind of 1. tide rectification method based on structure instantaneous water level model, it is characterised in that including step：Numerical simulation is carried out to the astronomical tidal level surveyed in area using ocean numerical models, obtained comprising including the measurement activity duration Instantaneous astronomical tide numerical model；According to survey tidal observation data structure that the tidal station that area is distributed observes it is instantaneous more than water level correction model and instantaneous residual water level repair Positive model；Wink is built according to instantaneous astronomical tide numerical model, instantaneous remaining water level correction model and the instantaneous residual water level correction model Shi Shuiwei models, the accurate correction of area's tidal level is surveyed by the instantaneous water level model realization.
- 2. tide rectification method according to claim 1, it is characterised in that the step is using ocean numerical models to surveying Astronomical tidal level in area carries out numerical simulation, obtains including the instantaneous astronomical tide numerical model including the measurement activity duration, specifically Including：Determine simulated domain and extract simulated domain water front data, make the grid that becomes more meticulous of simulated domain；Measure and calculate and survey area's bathymetric data, the bathymetric data is loaded into the grid that becomes more meticulous, obtains ocean numerical value Domain condition needed for pattern；Open on boundary point one year tide prediction data as simulation duration are provided using the harmonic constant forecast of partial tide, according to According to the demand from pattern, data are handled, obtain the open boundaryconditions of corresponding format；Select and the primary condition before the numerical model startup of ocean is set, high-precision number is carried out to the astronomical tide process surveyed in area Value simulation, obtains instantaneous astronomical tide numerical model.
- 3. tide rectification method according to claim 1, it is characterised in that the tidal station that the step is distributed according to area is surveyed The instantaneous remaining water level correction model of tidal observation data structure of observation specifically includes：Based on the shot and long term tidal station tidal observation data for surveying area's distribution, using harmonic analysis principle, mediation calculating is carried out to tidal level, from And obtain partial tide harmonic constant；Harmonic constant is modified, obtains the higher harmonic constant of stability, and tidal level is carried out based on stable harmonic constant Forecast；Forecast tidal level value is subtracted using tidal station actual measurement tidal level, obtains the remaining water level value of the tidal station, then based on all tidal stations Remaining water level value, using interpolation algorithm, obtain surveying remaining water level model in area.
- 4. tide rectification method according to claim 1, it is characterised in that the tidal station that the step is distributed according to area is surveyed The tidal observation data of observation, which build instantaneous residual water level correction model, also to be included：The remainder error in ocean numerical models simulation error and remaining water level extraction process in itself is extracted, using interpolation algorithm, Build instantaneous residual water level correction model.
- 5. tide rectification method according to claim 1, it is characterised in that the instantaneous remaining water level correction model expression formula For R (t)=h (t)-MSL-T (t)-Δ, wherein, R (t) represents Yu Shuiwei sequences, and h (t) represents measured water level sequence, and MSL is represented Mean sea level, T (t) represent astronomical tidal level, and Δ represents observation error.
- 6. tide rectification method according to claim 1, it is characterised in that the instantaneous residual water level correction model expression Formula is Ei(t)=hi(t)-Ti(t)-Ri(t), wherein, Ei(t) the remaining tide station modified value corresponding to tidal station i is represented；hi(t) Represent tidal station i actual measurement tidal level values, Ti(t) tidal station i simulation astronomical tide place value, R are representedi(t) tidal station i remaining water is represented Place value.
- 7. tide rectification method according to claim 1, it is characterised in that the instantaneous water level model expression is Tj(t) =hj(t)+Rj(t)+Cj(t), wherein, Tj(t) it is tide rectification value at depth measurement point j required by t, hj(t) it is at depth measurement point j The astronomical tidal level analogue value of t, Rj(t) be t at depth measurement point j remaining water level correction value, Cj(t) when being t at depth measurement point j The residual water level correction value at quarter.
- 8. a kind of storage medium, it is characterised in that be stored with a plurality of instruction, the instruction is suitable to by processor loading and right of execution Profit requires the step of tide rectification method based on structure instantaneous water level model described in 1-7 any one.
- 9. a kind of tide rectification device based on structure instantaneous water level model, it is characterised in that including processor, be adapted for carrying out each Instruction；And storage device, suitable for storing a plurality of instruction, the instruction is suitable to be loaded by processor and perform the claims The step of described in 1-7 any one based on the tide rectification method for building instantaneous water level model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711009326.0A CN107883931B (en) | 2017-10-25 | 2017-10-25 | Tidal level correction method and device based on instantaneous water level model construction and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711009326.0A CN107883931B (en) | 2017-10-25 | 2017-10-25 | Tidal level correction method and device based on instantaneous water level model construction and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107883931A true CN107883931A (en) | 2018-04-06 |
| CN107883931B CN107883931B (en) | 2020-04-28 |
Family
ID=61782318
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711009326.0A Active CN107883931B (en) | 2017-10-25 | 2017-10-25 | Tidal level correction method and device based on instantaneous water level model construction and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107883931B (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109033494A (en) * | 2018-06-01 | 2018-12-18 | 上海达华测绘有限公司 | A kind of coastal remote region tidal level projectional technique |
| CN111368406A (en) * | 2020-02-26 | 2020-07-03 | 天津大学 | Continuous depth datum plane construction method |
| CN111399012A (en) * | 2020-03-20 | 2020-07-10 | 中国矿业大学 | Method for monitoring reservoir water level by utilizing GNSS tri-frequency phase combined data |
| CN111723523A (en) * | 2020-06-22 | 2020-09-29 | 中山大学 | Estuary residual water level prediction method based on cascade neural network |
| CN111753461A (en) * | 2020-05-12 | 2020-10-09 | 中山大学 | Tidal water level correction method, target residual water level acquisition method, device and equipment |
| CN111899336A (en) * | 2020-07-30 | 2020-11-06 | 连云港港口工程设计研究院有限公司 | Method for obtaining minimum tide surface model of unknown sea area theory |
| CN112161670A (en) * | 2020-09-22 | 2021-01-01 | 长江水利委员会水文局长江上游水文水资源勘测局 | River channel water level correction method considering time sequence and spatial structure |
| CN113326470A (en) * | 2021-04-11 | 2021-08-31 | 桂林理工大学 | Remote sensing water depth inversion tidal height correction method |
| CN113324627A (en) * | 2021-08-04 | 2021-08-31 | 自然资源部第一海洋研究所 | Buoy-based sea level observation and prediction method |
| CN114693002A (en) * | 2022-05-23 | 2022-07-01 | 中国海洋大学 | Tide level prediction method, device, electronic equipment and computer storage medium |
| CN114969664A (en) * | 2022-06-01 | 2022-08-30 | 广州市城市规划勘测设计研究院 | Water level correction method, device, equipment and medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH03130624A (en) * | 1989-10-16 | 1991-06-04 | Marine Instr Co Ltd | Abnormal-wave removing circuit of wave height meter |
| CN102693355A (en) * | 2012-03-22 | 2012-09-26 | 国家海洋局第一海洋研究所 | High-precision water level calculation technology for rare tide station data control |
| CN104764445A (en) * | 2015-04-20 | 2015-07-08 | 中测新图(北京)遥感技术有限责任公司 | Method and device for determining coordinates of underwater object point |
| CN106679761A (en) * | 2017-01-06 | 2017-05-17 | 国家海洋信息中心 | Method and system for forecasting comprehensive water levels of sea areas in real time |
-
2017
- 2017-10-25 CN CN201711009326.0A patent/CN107883931B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH03130624A (en) * | 1989-10-16 | 1991-06-04 | Marine Instr Co Ltd | Abnormal-wave removing circuit of wave height meter |
| CN102693355A (en) * | 2012-03-22 | 2012-09-26 | 国家海洋局第一海洋研究所 | High-precision water level calculation technology for rare tide station data control |
| CN104764445A (en) * | 2015-04-20 | 2015-07-08 | 中测新图(北京)遥感技术有限责任公司 | Method and device for determining coordinates of underwater object point |
| CN106679761A (en) * | 2017-01-06 | 2017-05-17 | 国家海洋信息中心 | Method and system for forecasting comprehensive water levels of sea areas in real time |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109033494A (en) * | 2018-06-01 | 2018-12-18 | 上海达华测绘有限公司 | A kind of coastal remote region tidal level projectional technique |
| CN109033494B (en) * | 2018-06-01 | 2023-11-03 | 上海达华测绘有限公司 | Coastal remote area tide level calculation method |
| CN111368406A (en) * | 2020-02-26 | 2020-07-03 | 天津大学 | Continuous depth datum plane construction method |
| CN111399012A (en) * | 2020-03-20 | 2020-07-10 | 中国矿业大学 | Method for monitoring reservoir water level by utilizing GNSS tri-frequency phase combined data |
| CN111753461A (en) * | 2020-05-12 | 2020-10-09 | 中山大学 | Tidal water level correction method, target residual water level acquisition method, device and equipment |
| CN111723523A (en) * | 2020-06-22 | 2020-09-29 | 中山大学 | Estuary residual water level prediction method based on cascade neural network |
| CN111723523B (en) * | 2020-06-22 | 2024-04-30 | 中山大学 | Estuary surplus water level prediction method based on cascade neural network |
| CN111899336B (en) * | 2020-07-30 | 2021-10-19 | 连云港港口工程设计研究院有限公司 | Method for obtaining minimum tide surface model of unknown sea area theory |
| CN111899336A (en) * | 2020-07-30 | 2020-11-06 | 连云港港口工程设计研究院有限公司 | Method for obtaining minimum tide surface model of unknown sea area theory |
| CN112161670A (en) * | 2020-09-22 | 2021-01-01 | 长江水利委员会水文局长江上游水文水资源勘测局 | River channel water level correction method considering time sequence and spatial structure |
| CN113326470A (en) * | 2021-04-11 | 2021-08-31 | 桂林理工大学 | Remote sensing water depth inversion tidal height correction method |
| CN113326470B (en) * | 2021-04-11 | 2022-08-16 | 桂林理工大学 | Remote sensing water depth inversion tidal height correction method |
| CN113324627A (en) * | 2021-08-04 | 2021-08-31 | 自然资源部第一海洋研究所 | Buoy-based sea level observation and prediction method |
| CN114693002A (en) * | 2022-05-23 | 2022-07-01 | 中国海洋大学 | Tide level prediction method, device, electronic equipment and computer storage medium |
| CN114969664A (en) * | 2022-06-01 | 2022-08-30 | 广州市城市规划勘测设计研究院 | Water level correction method, device, equipment and medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107883931B (en) | 2020-04-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107883931A (en) | Tide rectification method and device, storage medium based on structure instantaneous water level model | |
| Johnson et al. | Ocean bottom pressure seasonal cycles and decadal trends from GRACE Release‐05: Ocean circulation implications | |
| Zijl et al. | Improved water-level forecasting for the Northwest European Shelf and North Sea through direct modelling of tide, surge and non-linear interaction | |
| Barton et al. | The Navy's Earth System Prediction Capability: A new global coupled atmosphere‐ocean‐sea ice prediction system designed for daily to subseasonal forecasting | |
| Sun et al. | Estimation of sea level variability in the China Sea and its vicinity using the SARIMA and LSTM models | |
| Vinayachandran et al. | Arabian Sea mini warm pool and the monsoon onset vortex | |
| Ishida et al. | Hourly-scale coastal sea level modeling in a changing climate using long short-term memory neural network | |
| Suursaar | Locally calibrated wave hindcasts in the Estonian coastal sea in 1966-2011 | |
| Vivier et al. | Causes of large‐scale sea level variations in the Southern Ocean: Analyses of sea level and a barotropic model | |
| Dong et al. | Toward improved parameterizations of reservoir operation in ungauged basins: A synergistic framework coupling satellite remote sensing, hydrologic modeling, and conceptual operation schemes | |
| Fu et al. | Improving the accuracy of barotropic and internal tides embedded in a high-resolution global ocean circulation model of MITgcm | |
| Zaron | Topographic and frictional controls on tides in the Sea of Okhotsk | |
| Chaudhuri et al. | Response of the Gulf Stream transport to characteristic high and low phases of the North Atlantic Oscillation | |
| Hudson et al. | Remote measurements of tides and river slope using an airborne lidar instrument | |
| Woodworth | Advances in the observation and understanding of changes in sea level and tides | |
| Zhang et al. | Numerical investigation of the control factors driving Zhe-Min Coastal Current | |
| Book et al. | Data assimilation modeling of the barotropic tides in the Korea/Tsushima Strait | |
| Dorfschäfer et al. | A new approach for estimating salinity in the Southwest Atlantic and its application in a data assimilation evaluation experiment | |
| Hirose et al. | Numerical simulation and satellite altimeter data assimilation of the Japan Sea circulation | |
| Grünler et al. | Estimation of tidally influenced estuarine river discharge from space using along‐track InSAR technology: A model‐based feasibility study | |
| Xu et al. | Tide model CST1 of China and its application for the water level reducer of bathymetric data | |
| Yang et al. | Non-linear interactions between tides and storm surges during extreme weather events over the eastern Canadian shelf | |
| McCardell et al. | Heat flux estimates from a synthesis of satellite observations and a hydrodynamic model (with application to Long Island Sound) | |
| Ren et al. | Data Assimilation with highfrequency (HF) radar surface currents at a marine renewable energy test site | |
| Baranes et al. | Tidally-driven interannual variation in extreme sea level probabilities in the Gulf of Maine |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |