CN108920877A - A kind of global wave method for numerical simulation based on MIKE21-SW model - Google Patents
A kind of global wave method for numerical simulation based on MIKE21-SW model Download PDFInfo
- Publication number
- CN108920877A CN108920877A CN201810873702.9A CN201810873702A CN108920877A CN 108920877 A CN108920877 A CN 108920877A CN 201810873702 A CN201810873702 A CN 201810873702A CN 108920877 A CN108920877 A CN 108920877A
- Authority
- CN
- China
- Prior art keywords
- wave
- data
- model
- offshore
- typhoon
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000004088 simulation Methods 0.000 title claims abstract description 42
- 238000005259 measurement Methods 0.000 claims abstract description 31
- 238000004519 manufacturing process Methods 0.000 claims abstract description 6
- 238000013480 data collection Methods 0.000 claims abstract description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 35
- 230000003993 interaction Effects 0.000 claims description 21
- VJJPUSNTGOMMGY-NBJJDLTASA-N (8aR,9R)-5-[[(2R,4aR,6R,7R,8R,8aS)-7,8-dihydroxy-2-methyl-4,4a,6,7,8,8a-hexahydropyrano[3,2-d][1,3]dioxin-6-yl]oxy]-9-(4-hydroxy-3,5-dimethoxyphenyl)-5a,6,8a,9-tetrahydro-5H-isobenzofuro[6,5-f][1,3]benzodioxol-8-one Chemical compound COC1=C(O)C(OC)=CC([C@@H]2C3=CC=4OCOC=4C=C3C(O[C@H]3[C@@H]([C@@H](O)[C@@H]4O[C@H](C)OC[C@H]4O3)O)C3[C@@H]2C(OC3)=O)=C1 VJJPUSNTGOMMGY-NBJJDLTASA-N 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 11
- 230000006854 communication Effects 0.000 claims description 9
- 238000000605 extraction Methods 0.000 claims description 9
- 239000013535 sea water Substances 0.000 claims description 9
- 230000002123 temporal effect Effects 0.000 claims description 7
- UZVHFVZFNXBMQJ-UHFFFAOYSA-N butalbital Chemical compound CC(C)CC1(CC=C)C(=O)NC(=O)NC1=O UZVHFVZFNXBMQJ-UHFFFAOYSA-N 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 6
- 230000002035 prolonged effect Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000010835 comparative analysis Methods 0.000 claims description 2
- 230000001429 stepping effect Effects 0.000 claims 2
- 239000006185 dispersion Substances 0.000 claims 1
- 238000012795 verification Methods 0.000 claims 1
- 230000002452 interceptive effect Effects 0.000 abstract description 2
- 238000012821 model calculation Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 8
- 230000008901 benefit Effects 0.000 description 6
- 241000208340 Araliaceae Species 0.000 description 4
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 4
- 235000003140 Panax quinquefolius Nutrition 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 235000008434 ginseng Nutrition 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000032258 transport Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004873 anchoring Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012950 reanalysis Methods 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Computer Security & Cryptography (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Educational Administration (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of global wave method for numerical simulation based on MIKE21-SW model, key step includes：Source data collection, landform production and model parameter setting, the calibration of model parameter, model calculation and the verifying of large-scale model regularly update building for wave time sequence (containing typhoon) and user interface platform, finally obtain a kind of global deep-sea and offshore, it is reliable, high-precision, extract the Wave parameters regularly updated used convenient for user and compose parametric data.The present invention can more accurately carry out the operation and verifying of model, improve the reliability of model simulation results by introducing the retrievable actual measurement landform of offshore and Wave Data；The present invention can preferably simulate the wave of offshore by the grid of encryption offshore；The present invention provides friendly user's interactive platform, the Wave Data for conveniently obtaining the whole world is allowed users to.
Description
Technical field
The present invention relates to natural environment analogue technique field more particularly to a kind of global waves based on MIKE21-SW model
Unrestrained method for numerical simulation.
Background technique
With economic globalization, water industry has obtained quick development, builds port level and significantly improves.Carrying out harbour
It in design and work progress, often needs to obtain the element of wave near offshore and harbour by the method for numerical simulation, and establishes
The wave model in the whole world can be convenient, efficiently obtain the wave situation in engineering sea area, mention for the design and construction of port engineering
Wave boundary condition is provided for wave parameter, and to wave model in offshore and port.
Currently used Numerical modeling of waves method mainly has three categories：First is that the wave number based on Boussinesq equation
It is worth model：The Numerical Wave Model of Boussinesq equation assumes that water particle is horizontal and vertical speed, being capable of simulated waves
Deformation in shallow water, refraction, the nonlinear interaction between diffraction and wave, but the general stability of the model is poor, is only applicable in
In small-scale Numerical modeling of waves；Second is that the Numerical Wave Model based on mild-slope equation：The model assumes that side slope is gentle
In the case of, with the method for perturbation expansion, three-dimensional problem is reduced to two-dimensional problems, is capable of shallow water effect, the refraction of simulated waves
With the physical processes such as diffraction, but the model is also only applicable to small-scale Numerical modeling of waves；Third is that being based on energy-balance equation
Numerical Wave Model：The model can various physical processes in simulated waves communication process, deformation in shallow water, folding such as wave
It penetrates with diffraction etc., due to using energy-balance equation, thus space-time step-length is not influenced by wavelength, without limitation, energy
Enough carry out large scale, prolonged Numerical modeling of waves.
Three kinds of modes are compared, since it is desired that the wave field in the whole world of simulation belongs to large scale, prolonged wave
Numerical simulation, it is higher for the stability requirement of model, so the theory of global wave numerical simulation is typically all to use at present
For energy-balance equation as governing equation, formula is as follows：
S=Sin+Sn1+Sds+Sbot+Ssurf
In formula, λ indicates that longitude coordinate, φ indicate that latitude coordinate, S indicate the summation in related with wave spectrum source and remittance, and N refers to
Wave action density spectrum, t refer to time, SinRefer to that disturbance interacts；Sn1Refer to nonlinear wave-wave interaction；SdsRefer to white
Cap loses item；SbotRefer to that bottom friction loses item；SsurfRefer to wave breaking item.
Currently based on energy-balance equation, existing and more commonly used global wave numerical model mainly has following three
Kind：
(1) NOAA global wave numerical model
NOAA global wave model is by U.S.National Oceanic and atmosphere office (National Oceanic and
Atmospheric Administration, abbreviation NOAA) exploitation global wave model, be based on third generation numerical Simulation of Ocean Wave
III model of method WaveWatch uses CFSR wind field as drive condition, energy using energy-balance equation as governing equation
The physical phenomenons such as refraction, deformation in shallow water, the whitecap of simulated waves.WaveWatch Numerical Wave Model is by Deflt university at first
Based on the Numerical Wave Model that energy-balance equation developed, later by continuously improving and perfect, it is widely used by now
Its third generation model, i.e. WaveWatch III.The model be using physical parameter scheme and accurate Scheme For Numerical Integration, and
The physical mechanism for considering wind driving etc., the numerical simulation of wave is carried out by structured grid.NOAA global wave numerical model
Temporal resolution be 3 hours, spatial resolution most of region be 0.5 ° × 0.5 °, certain respective regions have higher
Grid precision (such as mediterranean region, the northwest Indian Ocean).
The advantages of model, is to have obtained the Numerical modeling of waves data in the whole world, and disadvantage is it under typhoon effect
Wave simulation effect it is not satisfactory, it is not high simultaneously for the wave simulation result precision in offshore region.
(2) ECMWF global wave numerical model
ECMWF global wave numerical model is to develop (the European by European Study of Meso Scale Weather forecasting centre
Centre for Medium-Range Weather Forecasts), using energy-balance equation as its governing equation, make
The wind field of ERA wind field, the analytical wind field data of NCEP-NCARA and the sea NCEP-NCAR surface dynamic reanalysis is used to merge as drive
Dynamic wind field condition.The time resolution of ECMWF global wave numerical model be 6 hours, spatial resolution reach as high as 0.125 ° ×
0.125°。
The advantages of model, is that wave hindcasting data time section is longer, available 1990~2010 years wave hindcastings
Data, disadvantage are that the grid of model is thicker, can not accurately simulation table wind effect under wave condition；And for nearshore wave
Wave simulation effect it is poor.
(3) BMT global wave numerical model
BMT global wave numerical model be by Britain's Ocean Technology Centre (British Maritime Technology,
Abbreviation BMT) exploitation global wave numerical model, be based on III model of third generation numerical Simulation of Ocean Wave method WaveWatch, use
Energy-balance equation uses CFSR wind field as drive condition as governing equation, can the refractions of simulated waves, deformation in shallow water,
The physical phenomenons such as whitecap.The temporal resolution of BMT global wave model is 3 hours, and simulation has obtained 1979~2015 years 37 totally
The global wave numeric data in year.
According to similar, advantage is to have obtained the wave hindcasting in the whole world for the advantage and disadvantage of the model and count off after NOAA global wave
Data, disadvantage is that the wave simulation effect under acting on for typhoon is not satisfactory, simultaneously for the wave simulation in offshore region
As a result precision is not high.
With higher and higher for the precision of wave simulation in the development of Numerical modeling of waves technology and Practical Project,
Therefore urgent need is a kind of now being capable of accurate simulation typhoon and the global wave model of nearshore wave condition simultaneously.
Summary of the invention
It is a kind of for establishing the wave in the whole world the technical problem to be solved by the present invention is in view of the deficiencies of the prior art, provide
Numerical model, with obtain a kind of whole world, deep-sea and offshore, it is reliable, high-precision wave field, based on MIKE21-SW mould
The global wave method for numerical simulation of type.
In order to solve the above technical problems, the present invention adopts the following technical scheme that：A kind of whole world based on MIKE21-SW model
Numerical modeling of waves method, it is characterised in that：It carries out according to the following steps,
The first step, source data collection；
1) collection of C-MAP, ETOP and the obtainable actual measurement terrain data of offshore；
2) collection of wind field data, the wind field provide drive condition for model；
3) collection of the obtainable actual measurement Wave Data of satellite altimeter, buoy and offshore, with the verifying for model；
4) collection of ice sheet data；
Second step, the setting of landform production and Parameters in Mathematical Model；
1) landform makes：The drafting of terrain mesh；
2) setting of Parameters in Mathematical Model：The physical phenomenon for being simulated waves in communication process, including it is wind driving, non-thread
Property wave-wave interaction, the loss of white cap, bottom frictional resistance, wave breaking, ice covering, wave deformation in shallow water and wave refraction, and
Carry out the setting of relevant parameter；
Third step, the calibration of model parameter；
4th step, the operation of model and the verifying of large-scale model；
5th step post-processes obtained Wave Data, obtains the wave time sequence regularly updated containing typhoon
Column；
6th step carries out building for user interface platform, enables users to conveniently and efficiently obtain required Wave Data.
Drafting in the step 1) of second step, about terrain mesh：MIKE21-SW model uses non-structured triangle
Shape grid uses the grid of multinest in grid file manufacturing process, carries out mesh refinement for offshore region, in order to
The accurately wave condition of simulation offshore；
The method for drafting of grid file：The landform data boundary being collected into is imported into mesh generator Mesh first
In Generator, piecemeal then is carried out to the region of simulation, is broadly divided into three pieces：Profundal zone, offshore water and transition region；It is right
The size of grid is respectively set in each region, and offshore area grid is minimum, and profundal zone is maximum, finally combines mesh generator
Generate the non-structured triangular mesh of blockette.
In the step 2) of second step, model correlation is carried out for physical phenomenon of the accurate simulated waves in communication process
The setting of parameter, main physical phenomenon and analogy method are as follows in process of wave propagation：
A) wind drives Sin
Due to sea level upper air drag interaction, wind energy transformation is wave energy, and wind energy input equation is as follows：
Sin(f, θ)=max [alpha, gamma E (f, θ)]
γ=ε β σ x2
In formula, α, β respectively indicate linear and nonlinear rate of growth, and ε refers to the density ratio of air and seawater；σ refers to circle frequency relatively
Rate；C refers to phase velocity；θ, θWRespectively refer to wave direction and wind direction；
B) nonlinear wave-wave interaction Sn1
Sn1Item considers three-wave interaction and four-wave interaction simultaneously, is calculated using DIA method, and calculating formula is as follows：
Sn1(σ, θ)=Sn1+(σ, θ)+Sn1-(σ, θ)
Sn1+(σ, θ)=max [0,2 π αEBcgJ2| sin β | (cE (σ _, θ)) -2c_E (σ _, θ) E (σ, θ)]
Sn1-(σ, θ)=- 2Sn1+(σ, θ)
In formula, αEBRefer to reconciliation parameter；J refers to interaction parameter.
C) white cap loses Sds
Consider that the derivation of energy formula of white cap phenomenon loss is as follows in model：
D) bottom friction loses Sbot
In process of wave propagation, cause the movement of seawater, the friction in seawater and seabed also results in the loss of wave energy, calculates
Formula is as follows：
In formula, CfRefer to the coefficient of friction resistance.
E) wave breaking Ssurf
Wave may be crushed in communication process, region especially smaller in the offshore depth of water, it is necessary to examine
Consider the influence that wave energy is lost in wave breaking, the calculation formula of wave breaking item is as follows：
F) ice covers
For the sea area for having ice sheet to cover, ice sheet can have an impact the propagation of wave.It is covered for simulation ice, it is fixed in model
It is then treated as land, is not involved in calculating by justice ice sheet density when ice sheet density at the position of sea area is greater than definition value；
When ice sheet density at the position of sea area is less than definition value, then participate in calculating.
F) deformation in shallow water and wave refraction
During deformation in shallow water refers to that wave is propagated from deep water to offshore, the changed physical phenomenon of wave height.Wave
Refraction refers to that, as the depth of water shoals, when the isobath oblique of wave ray and sea-floor relief, the direction of wave can change
A kind of physical phenomenon.
In the third step, the calibration of model parameter is broadly divided into two steps：
A) trial operation of model：The trial operation of model, which refers to, carries out model setting using the model parameter of default, then transports
Row model；
B) verifying at model characteristic point：Several actual measurement Wave Data (including wave height, wave represented at position are chosen first
Period and wave direction)；Then the analogue value that corresponding position model is calculated is extracted；Finally measured value and model value are carried out pair
Than analysis, the error of the two is calculated.If the error represented between the measured value and model value at point is larger, model ginseng is modified
Number, re-starts the calibration of model, until error is within the range of permission；If representing between the measured value and model value at point
Error within allowed band, then carry out next step operation.In order to realize the purpose of verifying, need to introduce some judgement ginsengs
Number, such as related coefficient, expression formula are as follows：
In formula, Xi and Yi respectively refer to actual measurement and simulation wave characteristics, including wave height, period, wave direction；WithPoint
Do not refer to that the average value of actual measurement and the analogue value, n are total sample number.
In the 4th step, large-scale model verifying is the wave height that model is calculated, period, wave direction and is collected into
Field data obtained by satellite altimeter, buoy and offshore is verified, if the error between measured value and model value is larger,
Model parameter is then modified, calibration and the operation of model are re-started, until error is within the range of permission；If measured value and mould
Error between offset then carries out next step operation within allowed band.Herein with the verifying in third step model calibration
The main distinction is, that chooses herein is that large-scale data point is verified, and is chosen at several representative points in third step
Data verified；
In the step 1) of the first step, the collection of C-MAP, ETOP and the retrievable actual measurement terrain data of offshore is divided into following
Three steps carry out：
A) extraction of C-MAP terrain data：C-MAP is the topographic database carried in MIKE21 software, is integrated with the whole world
The electronic chart data of landform, for extracting the depth of water of deep water and offshore region and the data on land-based area boundary；
B) extraction of ETOP terrain data：ETOP data provide bathymetric data for large scale numerical model, and ETOP data are adopted
The whole world the ETOPO1 relief model for being 1 ' with resolution ratio includes land landform and ocean bathymetric data；
C) extraction of the retrievable actual measurement terrain data of offshore：The terrain data of offshore actual measurement refer in conjunction with it is specific it is proposed,
Project in construction, by using surveying vessel, using the terrain data for going back to the offshore region that sounding instrument equipment etc. obtains.
In the step 2) of the first step, the collection of wind field data is carried out in three steps：
A) downloading of CFSR wind field data
Select driving wind field of the CFSR wind field as model；The data type needed is selected first：Wind speed, wind direction；Then
The precision of grid is selected, obtainable grid precision includes 0.312 ° × 0.312 °, 0.5 ° × 0.5 °, 1.875 ° × 1.875 °
With 2.5 ° × 2.5 °, the highest wind field data of choice accuracy；It is divided into one hour between the data time of specified downloading, parameter has inputted
It after, that is, is downloaded, file format .grb2, the wind speed and direction data in file comprising needing.
B) downloading of typhoon data
Currently used typhoon data mainly have：American Association typhoon warning center (Joint Typhoon Warning
Center, abbreviation JTWC) best typhoon track data (1950~2017 years) and Japan Meteorological Agency (Japan
Meteorological Agency, abbreviation JMA) best typhoon track observed pile (1951~2017 years) etc..Below with
JTWC typhoon data instance illustrates the method for down loading of typhoon data：Selection needs the regions of data, mainly to have northwest too first
The typhoon data in flat ocean, the Southern Hemisphere and North Indian Ocean；Download typhoon track file and typhoon track, respectively txt format and kml
Format, typhoon track file are mainly the latitude and longitude coordinates for giving a certain field typhoon running track, maximum typhoon radius, maximum
The Typhoon Parameters such as wind speed, center of typhoon air pressure；Typhoon track is mainly loaded into google.earth, in order to intuitively see
Examine the running track of typhoon.
In the step 3) of the first step, the collection of the Wave Data of satellite altimeter, buoy and offshore actual measurement is divided into three steps
It carries out：
A) collection of satellite altimeter Wave Data：The Wave Data of satellite altimeter is the wave number measured by satellite
According to, including wave height, wave period and wave direction, it include the Wave Data that Jason-1, Jason-2 and Saral satellite measure, wherein
The time interval of the Wave Data of the Jason-1 satellite altimeter of collection is in January, 2002~2013 year June, and cycle period is
10 days, track band was 254；The time interval of the Wave Data of the Jason-2 satellite altimeter of collection is 2008~2018 years 7
Month, cycle period is 10 days, and track band is 254；The time interval of the Wave Data of the Saral satellite altimeter of collection is
In March, 2013~2018 year July, cycle period are 35 days, and track band is 1002；Satellite altimeter Wave Data collects process
For：
The satellite type needed is selected first：Jason-1, Jason-2 or Saral；Then the type of data is selected, altogether
There are three types of type, respectively Delayed-time, near-real time and tri- kinds of real time, to obtain prolonged wave
Unrestrained data and select Delayed-time；
B) collection of buoy data
Buoy data select the buoy data of American National buoy data center NDBC, which includes that 100 anchorings are floating
Mark and 60 C-MAN websites, are distributed in the Pacific Ocean and Atlantic Ocean area, provide the Wave Data that temporal resolution is one hour,
Including wave height, wave direction and period；
C) offshore surveys Wave Data：The Wave Data of offshore actual measurement is the Wave Data measured by wave observation station,
Including wave height, wave period and wave direction.
In the step 4) of the first step, the collection of ice sheet data is similar with the collection method of CFSR wind field data, generally
It is the ice sheet data that the whole world is obtained according to disclosed website data.
Compared to the prior art, the advantage of the invention is that：
(1) invention introduces the retrievable actual measurement landform of offshore and wave height data, model can more accurately be carried out
Verifying, improve modeling reliability of structure；
(2) present invention can preferably simulate the wave of offshore by the grid of encryption offshore；
(3) the present invention provides friendly user's interactive platform, the wave for conveniently obtaining the whole world is allowed users to
Unrestrained data；
(4) wave condition under typhoon effect can accurately be simulated.
Detailed description of the invention
Fig. 1 is global wave model realization schematic process flow diagram；
Fig. 2 is that CFSR wind field data download schematic diagram；
Fig. 3 is the grid schematic diagram that Mesh Generator is generated.
Specific embodiment
The present invention is described further With reference to embodiment.
Referring to Fig.1, in the present embodiment, the global wave method for numerical simulation based on MIKE21-SW model is specific to walk
It is rapid as follows：
1, the collection of C-MAP, ETOP and offshore actual measurement terrain data, is divided into following three step and carries out：
1.1, the extraction of C-MAP terrain data：C-MAP is the topographic database carried in MIKE21 software, is integrated with complete
The electronic chart data of ball landform is mainly used for extracting the depth of water in offshore region and the data on land-based area boundary, for the number at deep-sea
According to fewer；
1.2, the extraction of ETOP terrain data：ETOP data are currently used depth of water bank circle data, can be large scale number
Value model, which provides more accurately bathymetric data, ETOP data, many databases, such as ETOPO5, ETOPO2 and ETOPO1,
Current newest ETOP data are the whole world the ETOPO1 relief models that resolution ratio is 1 ', contain land landform and ocean depth of water number
According to.
1.3, the extraction of the retrievable actual measurement terrain data of offshore：The terrain data of offshore actual measurement refers in conjunction with specific quasi-
It builds, project in construction, by using surveying vessel, using the terrain data for going back to the offshore region that sounding instrument equipment etc. obtains.Such number
Generally higher according to precision, disadvantage is that offshore actual measurement terrain data relies on actual items data situation, can not achieve all standing.
2, the collection of wind field data, is mainly carried out in three steps：
2.1, the selection of wind field type
Wind field relatively common at present mainly has CCMP wind field and CFSR wind field, and CCMP wind field by NASA, (navigate by American National
Empty space agency) release, provided the oceanic winds field data since in the July, 1987, almost can with all oceans covering the whole world,
Spatial resolution is 0.25 ° × 0.25 °, the wind speed at temporal resolution 6h, 10m height；CFSR is U.S. NCEP
The third generation atmosphere analytical database again of (National Centers for Environment Prediction) publication,
CFSR has coupled the rear report of atmosphere-ocean-earth's surface-sea ice global simulation system as a result, the wind field data of CFSR contain
Wind speed and direction data at 1979~2010 each mesh points in the whole world above average sea level 10m height.Two kinds of wind fields
Comparatively, the temporal resolution of CFSR wind field is higher, thus select driving wind field of the CFSR wind field as this model.
2.2, the downloading of CFSR wind field data
First according to the disclosed website CFSR, the data type needed is selected：Wind speed, wind direction；Then the essence of grid is selected
Degree, obtainable grid precision have 0.312 ° × 0.312 °, 0.5 ° × 0.5 °, 1.875 ° × 1.875 ° and 2.5 ° × 2.5 °,
The highest wind field data of this choice accuracy；Next be divided into one hour between the data time of specified downloading, parameter input complete with
Afterwards, it can be downloaded, referring to Fig. 2, file format .grb2 contains the wind speed and direction data of needs in file.
2.3, the downloading of typhoon data
Currently used typhoon data mainly have：American Association typhoon warning center (Joint Typhoon Warning
Center, abbreviation JTWC) best typhoon track data (1950~2017 years) and Japan Meteorological Agency (Japan
Meteorological Agency, abbreviation JMA) best typhoon track observed pile (1951~2017 years).Below with
JTWC typhoon data instance illustrates the method for down loading of typhoon data：Selection needs the regions of data, mainly to have northwest too first
The typhoon data in flat ocean, the Southern Hemisphere and North Indian Ocean；Download typhoon track file and typhoon track, respectively txt format and kml
Format, typhoon file are mainly the latitude and longitude coordinates for giving a certain field typhoon operation, maximum typhoon radius, maximum wind velocity, platform
The Typhoon Parameters such as wind central pressure.
3, the collection of the Wave Data of satellite altimeter, buoy and offshore actual measurement is divided into the progress of three steps：
3.1, the collection of satellite altimeter Wave Data：The Wave Data of satellite altimeter refers to the wave measured by satellite
Data, including wave height, wave period and wave direction.The data of currently used satellite altimeter have Jason-1, Jason-2 and Saral
The Wave Data that satellite measures, wherein the time interval of the Wave Data of Jason-1 satellite altimeter is in January, 2002~2013
June in year, cycle period are 10 days or so, and track band is 254；The time interval of the Wave Data of Jason-2 satellite altimeter
For 2008~so far (in July, 2018), cycle period be 10 days or so, track band be 254；Saral satellite altimeter
The time interval of Wave Data is in March, 2013~so far (in July, 2018), and cycle period is 35 days, and track band is 1002.
It is as follows that satellite altimeter Wave Data collects process：
First according to the website of disclosed satellite altimeter, the satellite type needed is selected：Jason-1, Jason-2 or
Saral；Then the type for selecting data, altogether there are three types of type, respectively Delayed-time, near-real time and
Tri- kinds of real time, prolonged Wave Data, is typically chosen Delayed-time in order to obtain.The data of satellite altimeter
Scale is generally bigger, suitable for no buoy or the region of measured data.
3.2, the collection of buoy data
Buoy data more commonly used at present are American National buoy data center NDBC (The National Data
Center buoy data), the data include 100 or so mooring buoys and 60 or so C-MAN (Coastal
Marine Automated Network) website, it is mainly distributed on the Pacific Ocean and Atlantic Ocean area, temporal resolution can be provided
For one hour Wave Data (including wave height, wave direction and period), can by the website of American National buoy data center into
The downloading of row data.
3.3, the retrievable actual measurement Wave Data of offshore：The terrain data of offshore actual measurement is primarily referred to as through wave observation station
The Wave Data measured, including wave height, wave period and wave direction.The advantages of such data, is that precision is generally higher, and disadvantage is
Actual items data situation, generally short-term data data are relied on, space is difficult to realize gamut covering.
4, the drafting of terrain mesh：It is non-structured triangular mesh, grid spirit used by MIKE21-SW model
It is active strong, compared with structured grid (such as：Rectangular mesh), the effect that triangular mesh portrays boundary is more preferable.In grid text
The grid that multinest is used in part manufacturing process carries out mesh refinement for offshore region, in order to accurately simulate offshore
Wave condition.
The method for drafting of grid file：The landform data boundary being collected into is imported into mesh generator Mesh first
In Generator, piecemeal then is carried out to the region of simulation, can be divided mainly into three pieces：Profundal zone, offshore water and transition region；It is right
The size of grid is respectively set in each region, and offshore area grid is minimum；Profundal zone is maximum；Finally combine mesh generator
The non-structured triangular mesh of blockette is generated, the grid schematic diagram that MeshGenerator is generated is as shown in Figure 3.
5, physical phenomenon of the simulated waves in communication process and need to carry out the setting of Parameters in Mathematical Model, passed in wave
Main physical phenomenon and realization equation are as follows during broadcasting：
5.1, disturbance interaction Sin
Due to sea level upper air drag interaction, wind energy transformation is wave energy, and wind energy input equation is as follows：
Sin(f, θ)=max [alpha, gamma E (f, θ)]
γ=ε β σ x2
In formula, α, β respectively indicate linear and nonlinear rate of growth, and ε refers to the density ratio of air and seawater；σ refers to circle frequency relatively
Rate；C refers to phase velocity；θ, θWRespectively refer to wave direction and wind direction.
5.2, nonlinear wave-wave interaction Sn1
Sn1Item considers three-wave interaction and four-wave interaction simultaneously, is calculated using DIA method, and calculating formula is as follows：
Sn1(σ, θ)=Sn1+(σ, θ)+Sn1-(σ, θ)
Sn1+(σ, θ)=mas [0,2 π αEBcgJ2| sin β | (cE (σ -, θ)) -2c_E (σ -, θ) E (σ, θ)]
Sn1-(σ, θ)=- 2Sn1+(σ, θ)
In formula, αEBRefer to reconciliation parameter；J refers to interaction parameter.
5.3, white cap loses item Sds
Consider that the derivation of energy formula of white cap phenomenon loss is as follows in model：
5.4 bottom frictions lose item Sbot
In process of wave propagation, cause the movement of seawater, the friction in seawater and seabed also results in the loss of wave energy, calculates
Formula is as follows：
In formula, CfRefer to the coefficient of friction resistance.
5.5 wave breaking item Ssurf
Wave may be crushed in communication process, region especially smaller in the offshore depth of water, it is necessary to examine
Consider the influence that wave energy is lost in wave breaking, the calculation formula of wave breaking item is as follows：
The covering of 5.6 ice
For the sea area for having ice sheet to cover, ice sheet can have an impact the propagation of wave.It is covered for simulation ice, it is fixed in model
It is then treated as land, is not involved in calculating by justice ice sheet density when ice sheet density at the position of sea area is greater than definition value；
When ice sheet density at the position of sea area is less than definition value, then participate in calculating.
5.7 deformation in shallow water and wave refraction
During deformation in shallow water refers to that wave is propagated from deep water to offshore, the changed physical phenomenon of wave height.Wave
Refraction refers to that, as the depth of water shoals, when the isobath oblique of wave ray and sea-floor relief, the direction of wave can change
A kind of physical phenomenon.
6, the calibration of model parameter
The trial operation of 6.1 models：The trial operation of model, which refers to, carries out model setting using the model parameter of default, then transports
Row model；
Verifying at 6.2 model characteristic points：Choose first it is several represent at position actual measurement Wave Data (including wave height,
Wave period and wave direction)；Then the analogue value that corresponding position model is calculated is extracted；Finally measured value and model value are carried out
Comparative analysis calculates the error of the two.If the error represented between the measured value and model value at point is larger, model ginseng is modified
Number, re-starts the calibration of model, until error is within the range of permission；If representing between the measured value and model value at point
Error within allowed band, then carry out next step operation.In order to realize the purpose of verifying, need to introduce some judgement ginsengs
Number, such as related coefficient, expression formula are as follows：
In formula, Xi and Yi respectively refer to survey the wave characteristics with the analogue value, including wave height, period, wave direction；WithPoint
Do not refer to that the average value of actual measurement and the analogue value, n are total sample number.
7, model calculation and the verifying of large-scale model are the wave height that model is calculated, period, wave direction and are collected into
Satellite altimeter, field data is verified obtained by buoy and offshore, if the error between measured value and model value compared with
Greatly, then model parameter is modified, calibration and the operation of model are re-started, until error is within the range of permission；If measured value
Error between model value then carries out next step operation within allowed band.Herein with the verifying in step 6 model calibration
The main distinction be that chooses herein is that large-scale data point is verified, and step 6 is chosen at several representative points
Data are verified；
8, Wave Data is post-processed, the wave time sequence regularly updated (containing typhoon)
9, user interface platform building is carried out, can conveniently obtain required Wave Data convenient for user.
On the present invention is described in detail, it is described above, be only a preferred embodiment of the present invention, when cannot limit
Determine the scope of the present invention, i.e., it is all according to the made equivalent changes and modifications of the application range, it should still belong in covering scope of the present invention.
Claims (9)
1. a kind of global wave method for numerical simulation based on MIKE21-SW model, it is characterised in that：It carries out according to the following steps,
The first step, source data collection；
1) collection of C-MAP, ETOP and offshore actual measurement terrain data；
2) collection of CFSR wind field data and JTWC/JMA typhoon data；
3) collection of satellite altimeter, buoy and the retrievable actual measurement Wave Data of offshore, with the verifying for model；
4) collection of ice sheet data；
Second step, the setting of landform production and Parameters in Mathematical Model；
1) landform makes：The drafting of terrain mesh；
2) setting of Parameters in Mathematical Model：The physical phenomenon for being simulated waves in communication process, including wind driving, nonlinear wave
Wave interaction, the loss of white cap, bottom frictional resistance, wave breaking, ice covering, wave deformation in shallow water and wave refraction, and carry out
The setting of relevant parameter；
Third step, the calibration of model parameter；
4th step, the operation of model and the verifying of large-scale model；
5th step post-processes obtained Wave Data, obtains the wave time sequence regularly updated containing typhoon；
6th step carries out building for user interface platform, enables users to conveniently and efficiently obtain required Wave Data.
2. the global wave method for numerical simulation according to claim 1 based on MIKE21-SW model, it is characterised in that：
Drafting in the step 1) of second step, about terrain mesh：It is non-structured triangle used by MIKE21-SW model
Grid uses the grid of multinest in grid file manufacturing process, mesh refinement is carried out for offshore region, in order to standard
The wave condition of true simulation offshore；
The method for drafting of grid file：The landform data boundary being collected into is imported into mesh generator Mesh Generator first
In, piecemeal then is carried out to the region of simulation, is divided into three pieces：Profundal zone, offshore water and transition region；Each region is respectively set
The size of grid, offshore area grid is minimum, and profundal zone is maximum, and finally combining mesh generator is the non-knot for producing blockette
The triangular mesh of structure.
3. the global wave method for numerical simulation according to claim 1 based on MIKE21-SW model, it is characterised in that：
In the step 2) of second step, the setting of Parameters in Mathematical Model is carried out for physical phenomenon of the simulated waves in communication process,
Main physical phenomenon and implementation in process of wave propagation is as follows：
A) disturbance interaction Sin
Due to sea level upper air drag interaction, wind energy transformation is wave energy, and wind energy input equation is as follows：
Sin(f, θ)=max [alpha, gamma E (f, θ)]
γ=ε β σ x2
In formula, α, β respectively indicate linear and nonlinear rate of growth, and ε refers to the density ratio of air and seawater；σ refers to opposite circular frequency；c
Refer to phase velocity；θ, θWRespectively refer to wave direction and wind direction；
B) nonlinear wave-wave interaction Sn1
Sn1Item offshore region considers that three-wave interaction and far-reaching sea consider four-wave interaction, is calculated using DIA method, calculating formula
It is as follows：
Sn1(σ, θ)=Sn1+(σ, θ)+Sn1(σ, θ)
Sn1+(σ, θ)=mas [0,2 π αEBcgJ2| sin β | (cE (σ _, θ)) -2c_E (σ _, θ) E (σ, θ)]
Sn1-(σ, θ)=- 2Sn1+(σ, θ)
In formula, αEBRefer to reconciliation parameter；J refers to interaction parameter；
C) white cap loses item Sds
Consider that the derivation of energy formula of white cap phenomenon loss is as follows in model：
D) bottom friction loses item Sbot
In process of wave propagation, cause the movement of seawater, the friction in seawater and seabed also results in the loss of wave energy, calculation formula
It is as follows：
In formula, CfRefer to the coefficient of friction resistance；
E) wave breaking item Ssurf
Wave may be crushed in communication process, region especially smaller in the offshore depth of water, it is necessary to consider wave
The broken influence lost to wave energy of wave, the calculation formula of wave breaking item are as follows：
F) ice covers
For the sea area for having ice sheet to cover, ice sheet can have an impact the propagation of wave；It covers for simulation ice, is defined in model
It is then treated as land, is not involved in calculating by ice sheet density when ice sheet density at the position of sea area is greater than definition value；Work as sea
When ice sheet density at the position of domain is less than definition value, then participate in calculating；
G) deformation in shallow water and wave refraction
During deformation in shallow water refers to that wave is propagated from deep water to offshore, the changed physical phenomenon of wave height；Wave refraction
Refer to as the depth of water shoals, when the isobath oblique of wave ray and sea-floor relief, one kind that the direction of wave can change
Physical phenomenon.
4. the global wave method for numerical simulation according to claim 1 based on MIKE21-SW model, it is characterised in that：
In the third step, the calibration of model parameter is broadly divided into two steps：
A) trial operation of model：The trial operation of model, which refers to, carries out model setting using the model parameter of default, then runs mould
Type；
B) verifying at model characteristic point：Several, which are chosen, in offshore water first represents the actual measurement Wave Data at position, including
Wave height, wave period and wave direction；Then the analogue value that corresponding position model is calculated is extracted；Finally by measured value and model value
It compares and analyzes, calculates the error of the two；If the error represented between the measured value and model value at point is larger, mould is modified
Shape parameter re-starts the calibration of model, until error is within the range of permission；If representing the measured value and model value at point
Between error within allowed band, then carry out next step operation；When carrying out measured data and analogue value contrast verification, in addition to
Carry out sequence data comparative analysis when be also required to introduce some judgement parameters, comprising related coefficient, dispersion, the coefficient of variation and
Deviation etc., wherein the expression formula of related coefficient is as follows：
In formula, Xi and Yi respectively refer to survey the wave characteristics with the analogue value, including wave height, period, wave direction；WithIt respectively refers to
The average value of actual measurement and the analogue value, n is total sample number.
5. the global wave method for numerical simulation according to claim 1 based on MIKE21-SW model, it is characterised in that：
In the 4th step, large-scale model verifying is the wave height that model is calculated, period, wave direction and the satellite altitude being collected into
Meter, buoy and the retrievable field data of offshore are verified, if the error between measured value and model value is larger, modify mould
Shape parameter re-starts calibration and the operation of model, until error is within the range of permission；If between measured value and model value
Error within allowed band, then carry out next step operation；The difference with the verifying in third step model calibration is herein,
That chooses herein is that large-scale data point is verified, and is to choose data that several are represented at point to be tested in third step
Card.
6. the global wave method for numerical simulation according to claim 1 based on MIKE21-SW model, it is characterised in that：
In the step 1) of the first step, the collection of C-MAP, ETOP and the retrievable actual measurement terrain data of offshore are divided into following three stepping
Row：
A) extraction of C-MAP terrain data：C-MAP is the topographic database carried in MIKE21 software, is integrated with global landform
Electronic chart data, for extracting the depth of water of deep water and offshore region and the data on land-based area boundary；
B) extraction of ETOP terrain data：ETOP data provide bathymetric data for large scale numerical model, and ETOP data, which use, divides
The whole world the ETOPO1 relief model that resolution is 1 ' includes land landform and ocean bathymetric data；
C) extraction of the retrievable actual measurement terrain data of offshore：The terrain data of offshore actual measurement refers to that combination is specifically proposed, is building
Project, by using surveying vessel, using the terrain data for going back to the offshore region that sounding instrument equipment etc. obtains.
7. the global wave method for numerical simulation according to claim 1 based on MIKE21-SW model, it is characterised in that：
In the step 2) of the first step, the collection of wind field data is carried out in two steps：
A) downloading of CFSR wind field data
Select driving wind field of the CFSR wind field as model；The data type needed is selected first：Wind speed, wind direction；Then it selects
The precision of grid, obtainable grid precision includes 0.312 ° × 0.312 °, 0.5 ° × 0.5 °, 1.875 ° × 1.875 ° and
2.5 ° × 2.5 °, the highest wind field data of choice accuracy；It is divided into one hour between the data time of specified downloading, parameter input is completed
After, that is, it is downloaded, file format .grb2, the wind speed and direction data in file comprising needing；
B) downloading of typhoon data
Typhoon data used at present have：American Association typhoon warning center JTWC, best typhoon track data be 1950~
2017 and Japan Meteorological Agency JMA, best typhoon track observed pile was 1951~2017 years；Wherein JTWC typhoon data
Method for down loading be：Selection needs the regions of data first, includes the typhoon number of northwest Pacific, the Southern Hemisphere and North Indian Ocean
According to；Downloading typhoon track file and typhoon track, respectively txt format and kml format, typhoon track text give a certain field platform
The Typhoon Parameters such as the latitude and longitude coordinates of wind running track, maximum typhoon radius, maximum wind velocity, center of typhoon air pressure；Typhoon track
It is mainly loaded into google.earth, in order to intuitively observe the running track of typhoon.
8. the global wave method for numerical simulation according to claim 1 based on MIKE21-SW model, it is characterised in that：
In the step 3) of the first step, the collection of satellite altimeter, buoy and the retrievable actual measurement Wave Data of offshore is divided into three steppings
Row：
A) collection of satellite altimeter Wave Data：The Wave parameters of the measurement data of satellite altimeter include wave height, wave period and
Wave direction, more commonly used satellite data includes Jason-1, Jason-2 and Saral, wherein the Jason-1 satellite altitude collected
The time interval of the Wave Data of meter is in January, 2002~2013 year June, and cycle period is 10 days, and track band is 254；It receives
The time interval of the Wave Data of the Jason-2 satellite altimeter of collection is in July, 2008~2018, and cycle period is 10 days, rail
Mark band is 254；The time interval of the Wave Data of the Saral satellite altimeter of collection is in March, 2013~2018 year July,
Cycle period is 35 days, and track band is 1002；Satellite altimeter Wave Data collects process：
The satellite type needed is selected first：Jason-1, Jason-2 or Saral；Then the type of data is selected, one shares three
Tri- kinds of seed type, respectively Delayed-time, near-real time and real time, to obtain prolonged wave number
According to and select Delayed-time；
B) collection of buoy data
Buoy data select the buoy data of American National buoy data center NDBC, the data include 100 mooring buoys and
60 C-MAN websites, are distributed in the Pacific Ocean and Atlantic Ocean area, provide the Wave Data that temporal resolution is one hour, including
Wave height, wave direction and period；
C) offshore is obtainable surveys Wave Data：The Wave Data of offshore actual measurement is the wave number measured by wave observation station
According to, including wave height, wave period and wave direction.
9. the global wave method for numerical simulation according to claim 1 based on MIKE21-SW model, it is characterised in that：
In the step 4) of the first step, the collection of ice sheet data is similar with the collection method of CFSR wind field data, passes through disclosed website
Data come obtain the whole world ice sheet data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810873702.9A CN108920877A (en) | 2018-08-02 | 2018-08-02 | A kind of global wave method for numerical simulation based on MIKE21-SW model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810873702.9A CN108920877A (en) | 2018-08-02 | 2018-08-02 | A kind of global wave method for numerical simulation based on MIKE21-SW model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108920877A true CN108920877A (en) | 2018-11-30 |
Family
ID=64393993
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810873702.9A Pending CN108920877A (en) | 2018-08-02 | 2018-08-02 | A kind of global wave method for numerical simulation based on MIKE21-SW model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108920877A (en) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110728042A (en) * | 2019-09-29 | 2020-01-24 | 河海大学 | Method for rapidly forecasting long-time scale coast wave elements |
| CN111292214A (en) * | 2019-12-18 | 2020-06-16 | 中电投电力工程有限公司 | Method for calculating high-precision Chinese offshore wave characteristic distribution |
| CN112035429A (en) * | 2020-07-13 | 2020-12-04 | 国家海洋信息中心 | Multi-scale isophotic line automatic extraction method based on submarine topography model |
| CN112182759A (en) * | 2020-09-27 | 2021-01-05 | 中交第四航务工程勘察设计院有限公司 | Method for testing wave numerical simulation result based on satellite altimeter data |
| CN112818512A (en) * | 2021-01-04 | 2021-05-18 | 上海勘测设计研究院有限公司 | Storm mode-based maximum single-wave height Hmax calculation method |
| CN113051844A (en) * | 2021-03-12 | 2021-06-29 | 中国科学院烟台海岸带研究所 | Method for identifying and accurately tracing environmental pollution events of ocean and open water areas |
| CN113468773A (en) * | 2021-04-27 | 2021-10-01 | 交通运输部天津水运工程科学研究所 | Global table wind wave numerical simulation method and system |
| CN113553785A (en) * | 2021-07-14 | 2021-10-26 | 海博泰科技(青岛)有限公司 | Open wharf and harbor basin wave forecasting method |
| CN115048806A (en) * | 2022-06-30 | 2022-09-13 | 中交水运规划设计院有限公司 | Wind wave numerical simulation method |
| CN116611270A (en) * | 2023-07-20 | 2023-08-18 | 国家海洋局北海预报中心((国家海洋局青岛海洋预报台)(国家海洋局青岛海洋环境监测中心站)) | Typhoon wave real-time aggregate forecasting method, medium and system |
| CN117195775A (en) * | 2023-09-20 | 2023-12-08 | 上海勘测设计研究院有限公司 | Method, system, medium and device for calculating extreme mixed wave elements of offshore water area |
| CN117909666A (en) * | 2024-03-19 | 2024-04-19 | 青岛哈尔滨工程大学创新发展中心 | Intelligent sea wave correction method and system integrating numerical mode and deep learning |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130268520A1 (en) * | 2012-04-04 | 2013-10-10 | Microsoft Corporation | Incremental Visualization for Structured Data in an Enterprise-level Data Store |
| CN104331599A (en) * | 2014-09-30 | 2015-02-04 | 江苏省交通规划设计院股份有限公司 | Unstructured grid nesting wave numerical simulation method |
| CN105389468A (en) * | 2015-11-06 | 2016-03-09 | 中国海洋大学 | Rip current forecasting method |
| CN105824993A (en) * | 2016-03-10 | 2016-08-03 | 中国海洋大学 | Method for artificial sand beach, sand filling and maintenance engineering construction |
| CN108319772A (en) * | 2018-01-26 | 2018-07-24 | 中国科学院海洋研究所 | A kind of analysis method again of wave long term data |
-
2018
- 2018-08-02 CN CN201810873702.9A patent/CN108920877A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130268520A1 (en) * | 2012-04-04 | 2013-10-10 | Microsoft Corporation | Incremental Visualization for Structured Data in an Enterprise-level Data Store |
| CN104331599A (en) * | 2014-09-30 | 2015-02-04 | 江苏省交通规划设计院股份有限公司 | Unstructured grid nesting wave numerical simulation method |
| CN105389468A (en) * | 2015-11-06 | 2016-03-09 | 中国海洋大学 | Rip current forecasting method |
| CN105824993A (en) * | 2016-03-10 | 2016-08-03 | 中国海洋大学 | Method for artificial sand beach, sand filling and maintenance engineering construction |
| CN108319772A (en) * | 2018-01-26 | 2018-07-24 | 中国科学院海洋研究所 | A kind of analysis method again of wave long term data |
Non-Patent Citations (2)
| Title |
|---|
| 宗芳伊: "近20年南海波浪及波浪能分布、变化研究", 《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑(月刊)》 * |
| 谢鎏晔: "台风风场与波浪场的数值模拟研究", 《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑(月刊)》 * |
Cited By (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110728042B (en) * | 2019-09-29 | 2022-09-16 | 河海大学 | Method for rapidly forecasting long-time scale coast wave elements |
| CN110728042A (en) * | 2019-09-29 | 2020-01-24 | 河海大学 | Method for rapidly forecasting long-time scale coast wave elements |
| CN111292214A (en) * | 2019-12-18 | 2020-06-16 | 中电投电力工程有限公司 | Method for calculating high-precision Chinese offshore wave characteristic distribution |
| CN111292214B (en) * | 2019-12-18 | 2024-05-07 | 上海能源科技发展有限公司 | Method for calculating high-precision China offshore sea wave characteristic distribution |
| CN112035429A (en) * | 2020-07-13 | 2020-12-04 | 国家海洋信息中心 | Multi-scale isophotic line automatic extraction method based on submarine topography model |
| CN112035429B (en) * | 2020-07-13 | 2023-08-22 | 国家海洋信息中心 | Multi-scale isodepth line automatic extraction method based on submarine topography model |
| CN112182759A (en) * | 2020-09-27 | 2021-01-05 | 中交第四航务工程勘察设计院有限公司 | Method for testing wave numerical simulation result based on satellite altimeter data |
| CN112818512A (en) * | 2021-01-04 | 2021-05-18 | 上海勘测设计研究院有限公司 | Storm mode-based maximum single-wave height Hmax calculation method |
| CN112818512B (en) * | 2021-01-04 | 2022-09-30 | 上海勘测设计研究院有限公司 | Storm mode-based maximum single-wave height Hmax calculation method |
| CN113051844A (en) * | 2021-03-12 | 2021-06-29 | 中国科学院烟台海岸带研究所 | Method for identifying and accurately tracing environmental pollution events of ocean and open water areas |
| CN113051844B (en) * | 2021-03-12 | 2022-06-28 | 中国科学院烟台海岸带研究所 | Method for identifying and accurately tracing environmental pollution events of ocean and open water areas |
| CN113468773B (en) * | 2021-04-27 | 2023-02-03 | 交通运输部天津水运工程科学研究所 | Global table storm numerical simulation method and system |
| CN113468773A (en) * | 2021-04-27 | 2021-10-01 | 交通运输部天津水运工程科学研究所 | Global table wind wave numerical simulation method and system |
| CN113553785A (en) * | 2021-07-14 | 2021-10-26 | 海博泰科技(青岛)有限公司 | Open wharf and harbor basin wave forecasting method |
| CN113553785B (en) * | 2021-07-14 | 2023-12-26 | 海博泰科技(青岛)有限公司 | Open type wharf and harbor pool wave forecasting method |
| CN115048806A (en) * | 2022-06-30 | 2022-09-13 | 中交水运规划设计院有限公司 | Wind wave numerical simulation method |
| CN115048806B (en) * | 2022-06-30 | 2024-04-09 | 中交水运规划设计院有限公司 | Wind wave numerical simulation method |
| CN116611270A (en) * | 2023-07-20 | 2023-08-18 | 国家海洋局北海预报中心((国家海洋局青岛海洋预报台)(国家海洋局青岛海洋环境监测中心站)) | Typhoon wave real-time aggregate forecasting method, medium and system |
| CN116611270B (en) * | 2023-07-20 | 2023-10-03 | 国家海洋局北海预报中心((国家海洋局青岛海洋预报台)(国家海洋局青岛海洋环境监测中心站)) | Typhoon wave real-time aggregate forecasting method, medium and system |
| CN117195775A (en) * | 2023-09-20 | 2023-12-08 | 上海勘测设计研究院有限公司 | Method, system, medium and device for calculating extreme mixed wave elements of offshore water area |
| CN117195775B (en) * | 2023-09-20 | 2024-04-05 | 上海勘测设计研究院有限公司 | Method, system, medium and device for calculating extreme mixed wave elements of offshore water area |
| CN117909666A (en) * | 2024-03-19 | 2024-04-19 | 青岛哈尔滨工程大学创新发展中心 | Intelligent sea wave correction method and system integrating numerical mode and deep learning |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108920877A (en) | A kind of global wave method for numerical simulation based on MIKE21-SW model | |
| Garcia et al. | World ocean atlas 2018. Vol. 4: Dissolved inorganic nutrients (phosphate, nitrate and nitrate+ nitrite, silicate) | |
| Garcia et al. | World Ocean Atlas 2018, Volume 3: Dissolved Oxygen, Apparent Oxygen Utilization, and Dissolved Oxygen Saturation. | |
| Ho et al. | Seasonal variability of sea surface height in the South China Sea observed with TOPEX/Poseidon altimeter data | |
| 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 | |
| Allard et al. | The US Navy coupled ocean-wave prediction system | |
| CN110276160A (en) | A kind of region of no relief soil organic matter three-dimensional spatial distribution analogy method | |
| Suursaar | Locally calibrated wave hindcasts in the Estonian coastal sea in 1966-2011 | |
| Rossi et al. | Measurement of sea waves | |
| Harvey et al. | On the upward extension of the polar vortices into the mesosphere | |
| Cheung et al. | Numerical modeling and field evidence of coastal overwash in southern New England from Hurricane Bob and implications for paleotempestology | |
| Nittis et al. | Operational monitoring and forecasting in the Aegean Sea: system limitations and forecasting skill evaluation | |
| Suursaar et al. | Multi-scale analysis of wave conditions and coastal changes in the north-eastern Baltic Sea | |
| Medeiros et al. | A Random Forest model based on lidar and field measurements for parameterizing surface roughness in coastal modeling | |
| Daher et al. | Extraction of tide constituents by harmonic analysis using altimetry satellite data in the Brazilian coast | |
| Hadadpour et al. | Wave hindcasting in Anzali, Caspian Sea: a hybrid approach | |
| Janeković et al. | The Adriatic Sea M2 and K1 tides by 3D model and data assimilation | |
| Nikitin | Storage, processing and visualization data system of drifter observations of surface currents in the World Ocean | |
| Lebedev et al. | Comparative analysis of the North Atlantic surface circulation reproduced by three different methods | |
| Kumar et al. | Numerical simulation of typhoon wind forcing in the Korean seas using a spectral wave model | |
| Myers et al. | VDatum and strategies for national coverage | |
| Roldán et al. | Modified parametric hurricane wind model to improve the asymmetry in the region of maximum winds | |
| Orlić et al. | Report of the Croatian Committee of Geodesy and Geophysics on activities carried out between 2015 and 2018 | |
| Coppini et al. | The Mediterranean Forecasting System–Part 1: Evolution and performance | |
| Anthony et al. | Sea level changes along the coast of Sandakan Town, Sabah, Malaysia: projection and inundation coverage |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181130 |
|
| RJ01 | Rejection of invention patent application after publication |