CN106777949A - A kind of long-term trend Forecasting Methodology of the wave wave direction based on analyze data again - Google Patents
A kind of long-term trend Forecasting Methodology of the wave wave direction based on analyze data again Download PDFInfo
- Publication number
- CN106777949A CN106777949A CN201611121425.3A CN201611121425A CN106777949A CN 106777949 A CN106777949 A CN 106777949A CN 201611121425 A CN201611121425 A CN 201611121425A CN 106777949 A CN106777949 A CN 106777949A
- Authority
- CN
- China
- Prior art keywords
- wave
- wave direction
- sea
- matrix
- value
- 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
- 230000007774 longterm Effects 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000006243 chemical reaction Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 38
- 239000013598 vector Substances 0.000 claims description 16
- 230000003111 delayed effect Effects 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 claims description 4
- 108010014173 Factor X Proteins 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 238000012950 reanalysis Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 3
- 230000003628 erosive effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000006641 stabilisation Effects 0.000 description 3
- 238000011105 stabilization Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- 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
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of long-term trend Forecasting Methodology of the wave wave direction based on analyze data again, it is characterised in that step includes:First, time weather forecast data when collection ERA Interim are each;2nd, each lattice point coordinate is obtained;3rd, sea-level pressure gradient G X, GY anomaly value and standard deviation are calculated;4th, GX, GY anomaly value principal component are analyzed;5th, sea area data are carried out with Box Cox conversion;6th, the predictive factor of wave wave direction is calculated;7th, the standard deviation of wave direction and predictive factor is calculated;8th, predictive factor brings forecast model into;9th, wave direction lagged value brings model into;Tenth, GX, GY prediction on the basis of EOF;11, optimum choice predictive factor;12, model prediction wave wave direction;13, assessment prediction level;14, wave wave direction long-term trend are calculated;15, wave direction long-term trend figure is drawn.The long-term trend of secondary wave direction when the present invention can forecast many, and accuracy rate is high.
Description
Technical field
The invention belongs to ocean wave parameter forecasting technique field, more particularly to a kind of wave wave direction based on analyze data again
Long-term trend Forecasting Methodology.
Background technology
Wave has very important influence to the production and living of people, such as coastal port construction, waterway engineering all with
Wave has substantial connection.Wave wave direction is exactly the folder of an important parameter for reflecting wave feature, seashore dike line and wave wave direction
Angle, directly influence Coastal erosion and develop feature, therefore analysis prediction wave wave direction trend can help to it is more scientific
Harbour and coast protection works are built, is had important practical significance.Traditional observation method such as buoy etc., although can be accurate
Acquisition wave wave direction change information, but they can only obtain change of the wave in fixing point, and coverage rate also has very much
Limit, the buoy for being difficult to obtain the continuous wave of the sea more than 20 years in China Seas at present observes data.With satellite remote sensing
The maturation of technology, satellite data is gradually employed, though there is wider coverage about the satellite data of wave wave direction, at most
The simply data of nearly 20 years, this just seriously constrains the reliability to wave wave direction long-term trend research.How existing skill is overcome
The deficiency of art is one of problem demanding prompt solution in current ocean wave parameter forecasting technique field.
The content of the invention
Goal of the invention:A kind of the long-term of wave wave direction based on analyze data again is provided to overcome the deficiencies in the prior art
Trend forecasting method, the present invention is entered using Box-Cox conversion using the analyze data source again of global advanced stabilization to initial data
Row amendment, then according to meteorological datas such as revised sea-level pressure gradient, wave wave direction, using principal component analytical method and length
Phase wave direction trend formula, the long-term trend of secondary wave wave direction when calculating and predicting each, with very strong operability.
Technical scheme:A kind of long-term trend Forecasting Methodology of the wave wave direction based on analyze data again for proposing, its feature
It is to comprise the following specific steps that:
Step one, the ERA-Interim for collecting the pre- measured center of European Study of Meso Scale Weather based on mesh point mode analyzes number again
According to collection 20~30 years section it is each when time weather forecast data, wherein when each time weather forecast data to refer to 4~8 small
The sea-level pressure SLP and wave wave direction data of Shi Yici;
Step 2, the coordinate of time weather forecast data institute style point, with the coordinate as foundation, carries when obtaining collected each
Take with it is described each when time weather forecast data institute style point corresponding sea-level pressure the gradient matrix GX and GY of coordinate, such as
(1), shown in (2) formula, wave wave direction matrix θ, as shown in (3) formula, including m spatial point, each spatial point contains n sight
Survey data:
Wherein, GXmnSea-level pressure gradient when being the n-th of m-th spatial point time longitude coordinate direction value, GYmn
Sea-level pressure gradient when being the n-th of m-th spatial point time latitude coordinate direction value, θmnIt is m-th the n-th of spatial point
When time wave direction, m is the number of spatial point, when n is observation time.
Step 3, calculating be based on mesh point mode ERA-Interim it is each when time sea-level pressure gradient matrix GX
With the average M of GYXAnd MY, then subtract average M with the original value of sea-level pressure gradient matrix GX and GYXAnd MY, obtain based on lattice
Dot pattern it is each when time sea-level pressure gradient matrix GX and GY anomaly value PXAnd PY, and calculate sea-level pressure gradient
Matrix GX and GY anomaly value PXAnd PYStandard deviation SXAnd SY, as shown in (4), (5) formula:
In above-mentioned (4), (5) formula:Wherein, n represents secondary during observation, and i represents empty
Between point, i=1 ... m, j represent observation data, j=1 ... n;
Step 4, to sea-level pressure gradient matrix GX and GY anomaly value PXAnd PYBe respectively EOF analysis, obtain it is different into
Divide and each composition is to the contribution rate of population variance, retain preceding 30 EOF and principal component;Wherein:
To PXCovariance calculating is carried out, real symmetric matrix L is obtainedm×m, wherein:
The transposition of T representing matrixs;
Then covariance matrix L is soughtm×mCharacteristic vector V and characteristic value Λ, as shown in (6) formula, to meet LV=Λ V, its
In
Wherein, λ1≥λ2≥,...,≥λm(6),
Matrix V is orthogonal matrix, and the jth column element of matrix V is exactly eigenvalue λjCorresponding characteristic vector, wherein j values
From 1 to m;
According to real symmetric matrix Lm×mCharacteristic vector V and characteristic value Λ, calculate each characteristic vector variance contribution ratio and
The accumulative variance contribution ratio of preceding several characteristic vectors;L is ranked up according to characteristic value order from big to small, is made number one
Be EOF1, by that analogy;
Step 5, to collected according to step one and step 2 based on lattice point it is each when time wave wave direction matrix θ and sea
Plane barometric gradient matrix GX and GY carry out Box-Cox conversion, the wave wave direction tr θ after being convertedtWith sea-level pressure ladder
Degree trGXt、trGYt;
Step 6, to corresponding tr θ on each lattice pointt, with k-th principal component PCk,tWith k-th master of delayed 4 hours
Composition PCk,t-428 PC when calculating its coefficient correlation, and taking coefficient correlation highestk,tOr PCk,t-4As the pre- of wave wave direction
Survey the factor, wherein PCk,tOr PCk,t-4Represent principal component, k represents ordinal number, when t is represented time, t-4 represent delayed 4 hours when
It is secondary;
Step 7, calculates the standard deviation S of wave wave directionθ1With 30 predictive factor Xk,tStandard deviation SXk, preserve standby
With;
Step 8, brings the predictive factor that step 6 is obtained into forecast model, and i-th model and i-th are compared with F statistics
+ 1 model predicts the outcome, so as to select optimal predictive factor;
Step 9, the wave wave direction that will be late by a step also brings model into, as one of predictive factor, under integrated forecasting for the moment
The wave wave direction of secondary each lattice point, Optimized model parameter obtains final mask;Wherein model is as shown in (7) formula:
θ in above-mentioned (7) formulatIt is the wave wave direction by conversion on each mesh point, a is constant term, bkCorrespond to
Xk,tCoefficient, θt-pIt is the wave wave direction of delayed p, cpCorrespond to θt-pCoefficient, p is with the related parameter of predictand
Lag coefficient, Xk,tIt is k-th predictor based on SLP, utCan be represented with M ranks autoregression model, if M=0, utJust
It is white noise;
Step 10, on the basis of preceding 30 EOF that step 4 is obtained to it is each when time SLP gradient fields be predicted, obtain
To k-th principal component PCk,t;
Step 11, the S saved backup with step 7XkWeigh 30 predictive factor X of selectionk,t;
Step 12, the predictive factor that step 8 and step 11 are obtained is brought into the final mask of step 9, predicts mesh
Secondary wave wave direction when each in the timestamp phase, the wave direction value that will be predicted reverts to the value before Box-Cox conversion, saves as lattice point mould
Formula file;
Step 13, using the evaluation index assessment prediction level such as RMSE;
Step 14, with the wave wave direction of step 12 prediction as foundation, wave wave direction is calculated with trend computing formula
Long-term trend, finally give the long-term trend result of wave wave direction;
Step 15, according to the result of step 12, corresponds to corresponding lattice point coordinate, draws out wave wave direction and becomes for a long time
Gesture figure.
Beneficial effect:One is there is presently no the technical method for calculating wave wave direction long-term trend.The present invention is using the whole world
The analyze data source again of advanced stabilization, method is set up with decades even analyzing again across century-old wave wave direction data
On the basis of data, so as to solve the reliability of wave wave direction secular trend analysis;Two is that the present invention is converted using Box-Cox
Initial data is modified, then according to meteorological datas such as revised sea-level pressure gradient, wave wave direction, using EOF points
Analysis method and long-term wave direction trend formula, the long-term trend of secondary wave direction when calculating and predicting each;Three is that the present invention can help to
How the coastal building of arrangement or the trend of sea wall of science, and research mitigate erosion of the wave to seashore and waterfront structure
There is important scientific value, also can be used as the important evidence of the long term evolution rule of research littoral zone.
Brief description of the drawings
Fig. 1 is a kind of flow of the long-term trend Forecasting Methodology of wave wave direction based on analyze data again proposed by the present invention
Block diagram.
Fig. 2 is to be painted using a kind of long-term trend Forecasting Methodology of the wave wave direction based on analyze data again proposed by the present invention
The long-term trend result schematic diagram of the Pacific Ocean Partial Sea Area autumn wave wave direction of system.
Specific embodiment
Specific embodiment of the invention is described in further detail with reference to the accompanying drawings and examples.
Now by taking the Partial Sea Area of the Pacific Ocean as an example, using a kind of wave wave direction based on analyze data again proposed by the present invention
Long-term trend Forecasting Methodology carrys out the long-term trend of Study on Predicting Wave wave direction, and with reference to Fig. 1, its specific steps includes as follows:
Step one, collects the 1981- of the ERA-Interim reanalysis datasets of Chinese certain sea region based on mesh point mode
Time sea-level pressure SLP and wave wave direction data during each during 2000, the data break be every 6 hours once;
Step 2, obtains the coordinate of collected 6 hours data institute style points once, with the coordinate as foundation, extract with
The coordinate of time weather forecast data institute style point corresponding sea-level pressure gradient matrix GX and GY when described each, such as (1),
(2) shown in formula, wave wave direction matrix θ, as shown in (3) formula, including m spatial point, each spatial point contains n observation number
According to:
Step 3, secondary sea-level pressure gradient G X's and GY is equal when ERA-Interim of the calculating based on mesh point mode is each
Value MXAnd MY, then subtract average M with original value GX and GYXAnd MY, obtain based on mesh point mode it is each when time GX and GY anomaly
Value PXAnd PY, and calculate GX and GY anomaly values PXAnd PYStandard deviation SXAnd SY, as shown in (4), (5) formula:
In above-mentioned (4), (5) formula:Wherein, n represents secondary during observation, and i represents i-th
Individual spatial point, i=1 ... m, j represent j-th observation data, j=1 ... n.
Step 4, to GX and GY anomaly values PXAnd PYEOF analyses are done respectively, obtain heterogeneity and each composition to population variance
Contribution rate, retain preceding 30 EOF and principal component;Wherein:
To PXCovariance calculating is carried out, real symmetric matrix L is obtainedm×m, wherein:
The transposition of T representing matrixs;
Then covariance matrix L is soughtm×mCharacteristic vector V and characteristic value Λ, as shown in (6) formula, to meet LV=L Λ, its
In:
Wherein, λ1≥λ2≥,...,≥λm(6),
Matrix V is orthogonal matrix, and the jth column element of matrix V is exactly eigenvalue λjCorresponding characteristic vector, wherein j values
From 1 to m;
According to real symmetric matrix Lm×mCharacteristic vector V and characteristic value Λ, calculate each characteristic vector variance contribution ratio and
The accumulative variance contribution ratio of preceding several characteristic vectors, bigger corresponding characteristic vector and the time coefficient of representing of variance contribution is in data
Middle development law is more notable;L is ranked up according to characteristic value order from big to small, that make number one is EOF1, with this
Analogize;
Step 5, to collected according to step one and step 2 based on lattice point it is each when time wave wave direction matrix θ and sea
Plane barometric gradient GX and GY carry out Box-Cox conversion, the wave wave direction tr θ after being convertedtWith sea-level pressure gradient
trGXt、trGYt;
Step 6, to corresponding tr θ on each lattice pointt, with k-th principal component PCk,tWith k-th master of delayed 4 hours
Composition PCk,t-428 k-th principal component PC when calculating its coefficient correlation, and taking coefficient correlation highestk,tOr delayed 4 hours
K-th principal component PCk,t-4As the predictive factor of wave wave direction;
Step 7, calculates the standard deviation S θ of wave wave directionlWith 30 predictive factor Xk,tStandard deviation SXk, preserve standby
With;
Step 8, brings the predictive factor that step 6 is obtained into forecast model, and i-th model and i-th are compared with F statistics
+ 1 model predicts the outcome, so as to select optimal predictive factor;
Step 9, the wave wave direction that will be late by a step also brings model into, as one of predictive factor, under integrated forecasting for the moment
The wave wave direction of secondary each lattice point, Optimized model parameter obtains final mask;Wherein model is as shown in (7) formula:
θ in above-mentioned (7) formulatIt is the wave wave direction by conversion on each mesh point, a is constant term, bkCorrespond to
Xk,tCoefficient, θt-pIt is the wave wave direction of delayed p, cpCorrespond to θt-pCoefficient, p is with the related parameter of predictand
Lag coefficient, Xk,tIt is k-th predictive factor based on SLP, utCan be represented with M ranks autoregression model, if M=0, utJust
It is white noise;
Step 10, to 6 hours 2001- once on the basis of preceding 30 EOF of the 1981-2000 that step 4 is obtained
The SLP gradient fields of 2010 are predicted, and obtain PCk,t;
Step 11, the S saved backup with step 7XkWeigh 30 predictive factor X of selectionk,t;
Step 12, the predictive factor that step 8 and step 11 are obtained is brought into the final mask of step 9, prediction
Secondary wave wave direction when 2001-2010 is each, the wave direction value that will be predicted reverts to the value before Box-Cox conversion, saves as lattice
Dot pattern file;
Step 13, using the evaluation index assessment prediction level such as RMSE;
RMSE (root-mean-square error) is root-mean-square error, also known as standard error, and it is defined asI=1,2,3 ... n.In definite measured number of times, RMSE is represented with following formula:In formula, n is measurement
Number of times;diIt is one group of measured value and the deviation of average value.
Step 14, calculates wave wave direction long-term trend, and the wave wave direction predicted with step 12 uses trend as foundation
Computing formula is calculated, and finally gives the long-term trend result of wave wave direction;
Step 15, according to the result of step 12, corresponds to corresponding lattice point coordinate, draws out wave wave direction and becomes for a long time
Gesture figure.
Fig. 2 is to be painted using a kind of long-term trend Forecasting Methodology of the wave wave direction based on analyze data again proposed by the present invention
The short-term trend result schematic diagram of the Pacific waters summer wave wave direction of system, wherein abscissa are (radian/year).Fig. 2 can have
Effect instructs the sea wall of coastal area to arrange, operable for safeguarding that shore stabilization, prevention Coastal erosion have important scientific value
Property is strong.
All explanations not related to belong to techniques known in specific embodiment of the invention, refer to known skill
Art is carried out.The present invention can be corroded the prediction of the long-term trend of wave wave direction and prevention sea wall through validation trial
Play good directive function.Above specific embodiment and embodiment are to proposed by the present invention a kind of based on analyze data again
Wave wave direction long-term trend Forecasting Methodology technological thought specific support, it is impossible to protection scope of the present invention is limited with this,
It is every according to technological thought proposed by the present invention, any equivalent variations done on the basis of the technical program or equivalent change
It is dynamic, still fall within the scope of technical solution of the present invention protection.
Claims (3)
1. a kind of long-term trend Forecasting Methodology of the wave wave direction based on analyze data again, it is characterised in that including in detail below
Step:
Step one, collects the ERA-Interim reanalysis datasets of the pre- measured center of European Study of Meso Scale Weather based on mesh point mode
20~30 years section it is each when time weather forecast data, wherein time weather forecast data are referred to 4~8 hours one when each
Secondary sea-level pressure SLP and wave wave direction data;
Step 2, the coordinate of time weather forecast data institute style point when obtaining collected each, with the coordinate as foundation, extract with
The coordinate of time weather forecast data institute style point corresponding sea-level pressure gradient matrix GX and GY when described each, such as (1),
(2) shown in formula, wave wave direction matrix θ, as shown in (3) formula, including m spatial point, each spatial point contains n observation number
According to:
Wherein, GXmnSea-level pressure gradient when being the n-th of m-th spatial point time longitude coordinate direction value, GYmnIt is m
Individual spatial point n-th when time sea-level pressure gradient latitude coordinate direction value, θmnWhen being the n-th of m-th spatial point time
Wave direction, m is the number of spatial point, when n is observation time;
Step 3, secondary sea-level pressure gradient matrix GX's and GY is equal when ERA-Interim of the calculating based on mesh point mode is each
Value MXAnd MY, then subtract average M with the original value GX and GY of sea-level pressure gradient matrixXAnd MY, obtain based on mesh point mode
The anomaly value P of secondary sea-level pressure gradient matrix GX and GY when eachXAnd PY, and calculate sea-level pressure gradient matrix GX and
GY anomaly values PXAnd PYStandard deviation SXAnd SY, as shown in (4), (5) formula:
In above-mentioned (4), (5) formula:Wherein, n represents secondary during observation, i representation spaces point,
I=1 ... m, j represent j-th observation data, j=1 ... n;
Step 4, to sea-level pressure gradient matrix GX and GY anomaly value PXAnd PYBe respectively EOF analysis, obtain heterogeneity and
Each composition retains preceding 30 EOF and principal component to the contribution rate of population variance;Wherein:
To PXCovariance calculating is carried out, real symmetric matrix L is obtainedm×m, wherein:
The transposition of T representing matrixs;
Then covariance matrix L is soughtm×mCharacteristic vector V and characteristic value Λ, as shown in (6) formula, to meet LV=Λ V, wherein,
Matrix V is orthogonal matrix, and the jth column element of matrix V is exactly eigenvalue λjCorresponding characteristic vector, wherein j values from 1 to
m;
According to real symmetric matrix Lm×mCharacteristic vector V and characteristic value Λ, calculate the variance contribution ratio of each characteristic vector and former
The accumulative variance contribution ratio of individual characteristic vector;L is ranked up according to characteristic value order from big to small, what is made number one is
EOF1, by that analogy;
Step 5, to collected according to step one and step 2 based on lattice point it is each when time wave wave direction data and sea level gas
Pressure gradient matrix GX and GY carry out Box-Cox conversion, the wave wave direction tr θ after being convertedtWith sea-level pressure gradient trGXt、
trGYt;
Step 6, to corresponding tr θ on each lattice pointt, with k-th principal component PCK, tWith k-th principal component of delayed 4 hours
PCK, t-428 PC when calculating its coefficient correlation, and taking coefficient correlation highestK, tOr PCK, t-4As wave wave direction prediction because
Son;
Step 7, calculates the standard deviation of wave wave directionWith 30 predictive factor Xk,tStandard deviation SXk, save backup;
Step 8, brings the predictive factor that step 6 is obtained into forecast model, and i-th model and i+1 are compared with F statistics
Model predicts the outcome, so as to selecting optimal predictive factor;
Step 9, the wave wave direction that will be late by a step also brings model into, as one of predictive factor, a period of time time under integrated forecasting
The wave wave direction of each lattice point, Optimized model parameter obtains final mask;Wherein model is as shown in (7) formula:
θ in above-mentioned (7) formulatIt is the wave wave direction by conversion on each mesh point, a is constant term, bkCorrespond to Xk,t's
Coefficient, θt-pIt is the wave wave direction of delayed p, cpCorrespond to θt-pCoefficient, p is the delayed system with the related parameter of predictand
Number, Xk,tIt is k-th predictive factor based on SLP, utCan be represented with M ranks autoregression model, if M=0, utIt is exactly white noise
Sound;
Step 10, on the basis of preceding 30 EOF that step 4 is obtained to it is each when time SLP gradient fields be predicted, obtain
PCK, t;
Step 11, the S saved backup with step 7XkWeigh 30 predictive factor X of selectionk,t;
Step 12, the predictive factor that step 8 and step 11 are obtained is brought the final mask of step 9 into, during prediction target
Secondary wave wave direction when each in phase, the wave wave direction value that will be predicted reverts to the value before Box-Cox conversion, saves as lattice point mould
Formula file;
Step 13, using the evaluation index assessment prediction level such as RMSE;
Step 14, with the wave wave direction of step 12 prediction as foundation, the long-term of wave wave direction is calculated with trend computing formula
Trend, finally gives the long-term trend result of wave wave direction;
Step 15, according to the result of step 12, corresponds to corresponding lattice point coordinate, draws out wave wave direction long-term trend
Figure.
2. the long-term trend Forecasting Methodology of a kind of wave wave direction based on analyze data again according to claim 1, it is special
Levy and be, time weather forecast data refer to 6 hours sea-level pressure SLP and wave ripple once when each described in step one
To.
3. the long-term trend Forecasting Methodology of a kind of wave wave direction based on analyze data again according to claim 1 and 2, its
It is characterised by, the evaluation index such as RMSE described in step 13 refers to root-mean-square error, and it is defined asI=1,2,
3,…n;In definite measured number of times, RMSE is represented with following formula:In formula, n is pendulous frequency;diIt is one group of measurement
The deviation of value and average value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611121425.3A CN106777949B (en) | 2016-12-08 | 2016-12-08 | A kind of long-term trend prediction technique based on the wave wave direction for analyzing data again |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611121425.3A CN106777949B (en) | 2016-12-08 | 2016-12-08 | A kind of long-term trend prediction technique based on the wave wave direction for analyzing data again |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106777949A true CN106777949A (en) | 2017-05-31 |
CN106777949B CN106777949B (en) | 2019-09-10 |
Family
ID=58877319
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611121425.3A Expired - Fee Related CN106777949B (en) | 2016-12-08 | 2016-12-08 | A kind of long-term trend prediction technique based on the wave wave direction for analyzing data again |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106777949B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763160A (en) * | 2018-05-28 | 2018-11-06 | 河海大学 | Method and its device based on 20CR data prediction wave significant wave heights |
CN110633864A (en) * | 2019-09-20 | 2019-12-31 | 国网湖南省电力有限公司 | Wind speed numerical prediction correction method and system based on range deviation |
CN114819053A (en) * | 2022-03-04 | 2022-07-29 | 中国人民解放军国防科技大学 | Average wave direction forecast deviation correction method based on space-time convolution LSTM |
CN116453310A (en) * | 2023-06-14 | 2023-07-18 | 昆山工业大数据创新中心有限公司 | Method and device for alarming temperature abnormality of copper bar |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104036123A (en) * | 2014-05-29 | 2014-09-10 | 河海大学 | Short time trend predicating method for sea wave significant wave height based on reanalysis data |
CN104050514A (en) * | 2014-05-29 | 2014-09-17 | 河海大学 | Sea wave significant wave height long-term trend prediction method based on reanalysis data |
-
2016
- 2016-12-08 CN CN201611121425.3A patent/CN106777949B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104036123A (en) * | 2014-05-29 | 2014-09-10 | 河海大学 | Short time trend predicating method for sea wave significant wave height based on reanalysis data |
CN104050514A (en) * | 2014-05-29 | 2014-09-17 | 河海大学 | Sea wave significant wave height long-term trend prediction method based on reanalysis data |
Non-Patent Citations (2)
Title |
---|
FOLLEY M ET.: "Analysis of the nearshore wave energy resource", 《RENEWABLE ENERGY》 * |
万勇等: "基于ERA-Interim再分析数据的OE-W01区块波浪能资源评估", 《资源科学》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763160A (en) * | 2018-05-28 | 2018-11-06 | 河海大学 | Method and its device based on 20CR data prediction wave significant wave heights |
CN110633864A (en) * | 2019-09-20 | 2019-12-31 | 国网湖南省电力有限公司 | Wind speed numerical prediction correction method and system based on range deviation |
CN110633864B (en) * | 2019-09-20 | 2023-10-10 | 国网湖南省电力有限公司 | Wind speed numerical forecast correction method and system based on distance flat deviation |
CN114819053A (en) * | 2022-03-04 | 2022-07-29 | 中国人民解放军国防科技大学 | Average wave direction forecast deviation correction method based on space-time convolution LSTM |
CN114819053B (en) * | 2022-03-04 | 2024-05-03 | 中国人民解放军国防科技大学 | Average wave direction forecast deviation correcting method based on space-time convolution LSTM |
CN116453310A (en) * | 2023-06-14 | 2023-07-18 | 昆山工业大数据创新中心有限公司 | Method and device for alarming temperature abnormality of copper bar |
CN116453310B (en) * | 2023-06-14 | 2023-09-05 | 昆山工业大数据创新中心有限公司 | Method and device for alarming temperature abnormality of copper bar |
Also Published As
Publication number | Publication date |
---|---|
CN106777949B (en) | 2019-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104050514B (en) | A kind of long-term trend Forecasting Methodology of the wave significant wave height based on analyze data again | |
Fang et al. | Employing extreme value theory to establish nutrient criteria in bay waters: a case study of Xiangshan Bay | |
Zhao et al. | Three-dimensional hydrodynamic and water quality model for TMDL development of Lake Fuxian, China | |
Han et al. | A convolutional neural network using surface data to predict subsurface temperatures in the Pacific Ocean | |
Drushka et al. | In situ observations of Madden–Julian oscillation mixed layer dynamics in the Indian and western Pacific Oceans | |
CN106777949B (en) | A kind of long-term trend prediction technique based on the wave wave direction for analyzing data again | |
Placke et al. | Long-term mean circulation of the Baltic Sea as represented by various ocean circulation models | |
CN109543356A (en) | Consider the ocean interior temperature-salinity structure remote sensing inversion method of Space atmosphere | |
CN104021308A (en) | ERA (ECMWF (European Centre for Medium-range Weather Forecasts) Reanalysis)-Interim and ERA40-based method for predicting significant wave height of sea waves | |
Abadie et al. | Wave climate and longshore drift on the South Aquitaine coast | |
CN108268969A (en) | Regional Economic Development form analysis and Forecasting Methodology and device based on remotely-sensed data | |
CN104036123B (en) | Short time trend predicating method for sea wave significant wave height based on reanalysis data | |
KR101264873B1 (en) | Wind power density prediction method using stepewise regression procedure | |
CN115540828A (en) | Internal wave forecasting method based on wall sensor | |
CN108549961A (en) | A method of wave significant wave height is estimated based on CMIP5 | |
Li | Evaluation method of ecological vulnerability of scenic spots based on entropy weight TOPSIS model | |
Chaudhuri et al. | Response of the Gulf Stream transport to characteristic high and low phases of the North Atlantic Oscillation | |
Deng et al. | A method for assessing the coastline recession due to the sea level rise by assuming stationary wind-wave climate | |
CN117114176A (en) | Land utilization change prediction method and system based on data analysis and machine learning | |
Wang et al. | Xihe: A data-driven model for global ocean eddy-resolving forecasting | |
Wang et al. | Long-term variation of storm surge-associated waves in the Bohai Sea | |
CN116562448A (en) | Quick prediction method for flood inundation of flood storage area based on coupling of grid analysis and artificial intelligence technology | |
CN104021434A (en) | Method for forecasting sea wave significant wave height based on ERA-Interim | |
CN106779206A (en) | A kind of short-term trend Forecasting Methodology of the wave wave direction based on analyze data again | |
Bakış et al. | Analysis and comparison of spatial rainfall distribution applying different interpolation methods in Porsuk river basin, Turkey |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190910 |
|
CF01 | Termination of patent right due to non-payment of annual fee |