CN112417705A - Island sea area extreme value wave calculation method based on environment contour line - Google Patents
Island sea area extreme value wave calculation method based on environment contour line Download PDFInfo
- Publication number
- CN112417705A CN112417705A CN202011418288.6A CN202011418288A CN112417705A CN 112417705 A CN112417705 A CN 112417705A CN 202011418288 A CN202011418288 A CN 202011418288A CN 112417705 A CN112417705 A CN 112417705A
- Authority
- CN
- China
- Prior art keywords
- wave height
- period
- distribution
- wave
- contour line
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 12
- 230000000737 periodic effect Effects 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 claims 3
- 239000003643 water by type Substances 0.000 claims 1
- 239000000463 material Substances 0.000 description 2
- 235000014653 Carica parviflora Nutrition 0.000 description 1
- 241000243321 Cnidaria Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
Images
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses an island sea area extreme value wave calculation method based on an environment contour line, which comprises the steps of selecting an island sea area target position, obtaining wave height and periodic variation values of a long-period time sequence through field measurement or value postreporting, obtaining a wave height period joint distribution result through statistical processing, and calculating the joint distribution result of the wave height extreme value and the period extreme value. The data result evaluated by the invention is suitable for the work related to extreme value wave element determination in the fields of offshore engineering, ocean engineering and island-like reef engineering, and the calculation method is convenient and effective.
Description
Technical Field
The invention relates to the technical field of ocean engineering, in particular to an island reef sea area extreme value wave calculation method based on an environment contour line.
Background
When waves are transmitted to the island sea area from deep water or a limited water depth area, due to the influence of the island terrain, the waves are broken and bottom materials are rubbed (including permeation), the wave attenuation characteristics are obvious, and seabed bottom materials in different forms, such as silt, gravel or coral reefs, have larger differences, so that the extreme value characteristics of the waves in the island sea area are influenced.
The conventional extreme value wave analysis methods include an extreme value distribution model, a Weibull distribution model, a Pearson III type distribution model and a lognormal distribution model, and the methods fit the wave height probability distribution and calculate the wave height extreme values in different recurrence periods, but the selection of the period extreme values is not described by a more exact model, particularly, the evaluation of the wave height extreme values and the period extreme values is relatively independent, and a key problem is how to establish a model of the wave height and period joint relation and a calculation method for carrying out the wave height extreme value and period extreme value joint distribution by using the model, so that how to quantitatively evaluate the wave extreme value conditions of an island reef sea area, namely the joint distribution of the wave height extreme values and the period extreme values, is a problem which needs to be solved urgently in the fields of offshore engineering, ocean engineering and island reef engineering.
Disclosure of Invention
The invention provides an island sea area extreme value wave calculation method based on an environment contour line aiming at the defects in the background technology.
The invention aims to solve the phenomenon, adopts the following technical scheme, and discloses an island sea area extreme value wave calculation method based on an environment contour line, which comprises the following steps:
s1, selecting object: taking the selected points in the island sea area as objects, obtaining wave height and period change values of the long-period time sequence through field measurement or numerical value postreporting, and obtaining a wave height period combined distribution result through statistical processing;
s2, calculating a wave height probability distribution expression: fitting wave height distribution by adopting a Weibull model, solving fitting empirical parameters by a least square method, and obtaining a probability distribution expression of the wave height;
s3, calculating a conditional probability distribution expression of the period under the given wave height: solving fitting empirical parameters by a least square method by adopting a lognormal distribution model to obtain a conditional probability distribution expression of the period under the given wave height;
s4, calculating the probability density expression of the wave height and period combined distribution: multiplying the probability distribution expression of the wave height and the conditional probability distribution expression of the period based on S1 and S2 to obtain a probability density expression of the wave height and period combined distribution;
s5, determine boundary: setting a wave steepness coefficient, and determining a wave steepness limiting line;
s6, calculating the wave height period joint distribution contour line: obtaining extreme value wave height and period combined distribution contour lines in different reappearance periods by adopting an environment contour line method based on S4; based on the steep limit line of S5, identifying the wave height and period combination value corresponding to the steep limit line in the contour line: and judging that the combination values are unreasonable by the combination values of the large wave height and the small period, and finally obtaining a reasonable wave height period combined distribution contour line.
As a further preferable mode of the present invention, in step S2, the wave height probability obeys weibull distribution, the period conditional probability density function obeys lognormal distribution, the wave height period joint distribution probability is a product of the wave height probability and the period conditional probability, and the wave height probability expression is as follows
In the formula (1), h represents a wave height, α, β, and γ are distribution parameters, and the formula (1) is a three-parameter weibull distribution, and if γ is 0, a two-parameter weibull distribution is adopted.
As a further preferable mode of the present invention, in step S3, the period conditional probability density expression is as follows
The mu and sigma in the formula (2) need to be obtained by actually measuring data or numerical values and then reporting the results, and the fitting formula is as follows
In the formulae (3) and (4), a0,a1,a2,b0,b1And b2Is a fitting parameter obtained by fitting a curve by a nonlinear least square method, and therefore, the joint probability distribution can be expressed as
As a further preferred embodiment of the present invention, in step S5, after the joint probability distribution expression is determined, the wave height probability distribution and the periodic conditional probability distribution are converted into the standard normal U parameter space, respectively
A circle is then established in the U-shaped parameter space that specifies the recurrence period. Taking the observed value recorded every 3 hours as an example, the circle radius of the contour line corresponding to 100 years of the recurrence period is
Finally, the circle is converted into the wave height and the period contour line corresponding to the environment parameter space
The coefficient of steepness is defined as
α＝(h/l)b·tanhkd (11)
(h/l) in the formula (11)bExpressed as the ratio of the height of the breaking wave to the wavelength, k is the wave number and d is the water depth. The steepness factor here takes the value 0.12.
The method calculates the extreme value wave condition based on the environment contour line, obtains the wave heights and the period combined distribution of different recurrence periods by taking the on-site wave measurement result or the wave numerical value post-reporting result as input, and is suitable for the work related to extreme value wave element determination in the fields of near-shore engineering, ocean engineering and near-island reef engineering.
Drawings
FIG. 1 is a schematic view of a station layout according to the present invention;
FIG. 2 is a fitting result of the wave height probability distribution of the present invention;
FIG. 3 is a periodic conditional probability fit for a given wave height according to the present invention;
FIG. 4 is a fitting result of the wave height and period joint distribution of the present invention;
FIG. 5 is a profile of the extreme wave height and period combination distribution for different reconstruction periods of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention provides a technical scheme that: an island sea area extreme value wave calculation method based on an environment contour line comprises the following steps:
s1, selecting object: taking the selected points in the island sea area as objects, obtaining wave height and period change values of the long-period time sequence through field measurement or numerical value postreporting, and obtaining a wave height period combined distribution result through statistical processing;
s2, calculating a wave height probability distribution expression: fitting wave height distribution by adopting a Weibull model, solving fitting empirical parameters by a least square method, and obtaining a probability distribution expression of the wave height;
s3, calculating a conditional probability distribution expression of the period under the given wave height: solving fitting empirical parameters by a least square method by adopting a lognormal distribution model to obtain a conditional probability distribution expression of the period under the given wave height;
s4, calculating the probability density expression of the wave height and period combined distribution: multiplying the probability distribution expression of the wave height and the conditional probability distribution expression of the period based on S1 and S2 to obtain a probability density expression of the wave height and period combined distribution;
s5, determine boundary: setting a wave steepness coefficient, and determining a wave steepness limiting line;
s6, calculating the wave height period joint distribution contour line: obtaining extreme value wave height and period combined distribution contour lines in different reappearance periods by adopting an environment contour line method based on S4; based on the steep limit line of S5, identifying the wave height and period combination value corresponding to the steep limit line in the contour line: and judging that the combination values are unreasonable by the combination values of the large wave height and the small period, and finally obtaining a reasonable wave height period combined distribution contour line.
In step S2, the wave height probability obeys weibull distribution, the period conditional probability density function obeys lognormal distribution, the wave height period joint distribution probability is the product of the wave height probability and the period conditional probability, and the wave height probability expression is as follows
In the formula (1), h represents a wave height, α, β, and γ are distribution parameters, and the formula (1) is a three-parameter weibull distribution, and if γ is 0, a two-parameter weibull distribution is adopted.
In step S3, the periodic conditional probability density expression is as follows
The mu and sigma in the formula (2) need to be obtained by actually measuring data or numerical values and then reporting the results, and the fitting formula is as follows
In the formulae (3) and (4), a0,a1,a2,b0,b1And b2Is a fitting parameter obtained by fitting a curve by a nonlinear least square method, and therefore, the joint probability distribution can be expressed as
In step S5, after the joint probability distribution expression is determined, the wave height probability distribution and the periodic condition probability distribution are respectively converted into a standard normal U parameter space
A circle is then established in the U-shaped parameter space that specifies the recurrence period. Taking the observed value recorded every 3 hours as an example, the circle radius of the contour line corresponding to 100 years of the recurrence period is
Finally, the circle is converted into the wave height and the period contour line corresponding to the environment parameter space
The coefficient of steepness is defined as
α＝(h/l)b·tanhkd (22)
(h/l) in the formula (11)bExpressed as the ratio of the height of the breaking wave to the wavelength, k is the wave number and d is the water depth. The steepness factor here takes the value 0.12.
In conclusion, the wave extreme value condition is calculated based on the environment contour line, the wave height and period combined distribution of different recurrence periods is obtained by taking the wave on-site measurement result or the wave numerical value post-reporting result as input, the estimated data result is suitable for the work related to extreme value wave element determination in the field of offshore engineering, ocean engineering and island-like reef engineering, and the calculation method is convenient and effective.
While there have been shown and described what are at present considered the fundamental principles and essential features of the invention and its advantages, it will be apparent to those skilled in the art that the invention is not limited to the details of the foregoing exemplary embodiments, but is capable of other specific forms without departing from the spirit or essential characteristics thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.
Furthermore, it should be understood that although the present description refers to embodiments, not every embodiment may contain only a single embodiment, and such description is for clarity only, and those skilled in the art should integrate the description, and the embodiments may be combined as appropriate to form other embodiments understood by those skilled in the art.
Claims (4)
1. An island sea area extreme value wave calculation method based on an environment contour line is characterized by comprising the following calculation steps:
s1, selecting object: taking the selected points in the island sea area as objects, obtaining wave height and period change values of the long-period time sequence through field measurement or numerical value postreporting, and obtaining a wave height period combined distribution result through statistical processing;
s2, calculating a wave height probability distribution expression: fitting wave height distribution by adopting a Weibull model, solving fitting empirical parameters by a least square method, and obtaining a probability distribution expression of the wave height;
s3, calculating a conditional probability distribution expression of the period under the given wave height: solving fitting empirical parameters by a least square method by adopting a lognormal distribution model to obtain a conditional probability distribution expression of the period under the given wave height;
s4, calculating the probability density expression of the wave height and period combined distribution: multiplying the probability distribution expression of the wave height and the conditional probability distribution expression of the period based on S1 and S2 to obtain a probability density expression of the wave height and period combined distribution;
s5, determine boundary: setting a wave steepness coefficient, and determining a wave steepness limiting line;
s6, calculating the wave height period joint distribution contour line: obtaining extreme value wave height and period combined distribution contour lines in different reappearance periods by adopting an environment contour line method based on S4; based on the steep limit line of S5, identifying the wave height and period combination value corresponding to the steep limit line in the contour line: and judging that the combination values are unreasonable by the combination values of the large wave height and the small period, and finally obtaining a reasonable wave height period combined distribution contour line.
2. The method for calculating the extreme value wave of the island and reef sea area based on the environmental contour line as claimed in claim 1, wherein in step S2, the wave height probability obeys weibull distribution, the period conditional probability density function obeys log normal distribution, the wave height period joint distribution probability is the product of the wave height probability and the period conditional probability, and the wave height probability expression is as follows
In the formula (1), h represents a wave height, α, β, and γ are distribution parameters, and the formula (1) is a three-parameter weibull distribution, and if γ is 0, a two-parameter weibull distribution is adopted.
3. The method for estimating extreme waves of island waters based on environmental profile according to claim 1, wherein in step S3, the periodic conditional probability density expression is as follows
The mu and sigma in the formula (2) need to be obtained by actually measuring data or numerical values and then reporting the results, and the fitting formula is as follows
In the formulae (3) and (4), a0,a1,a2,b0,b1And b2Is a fitting parameter obtained by fitting a curve by a nonlinear least square method, and therefore, the joint probability distribution can be expressed as
4. The method for calculating the extreme value wave of the island ocean based on the environmental contour line according to claim 1, wherein in step S5, after the joint probability distribution expression is determined, the wave height probability distribution and the periodic conditional probability distribution are respectively converted into the standard normal U parameter space
A circle is then established in the U-shaped parameter space that specifies the recurrence period. Taking the observed value recorded every 3 hours as an example, the circle radius of the contour line corresponding to 100 years of the recurrence period is
Finally, the circle is converted into the wave height and the period contour line corresponding to the environment parameter space
The coefficient of steepness is defined as
α＝(h/l)b·tanh kd (11)
(h/l) in the formula (11)bExpressed as the ratio of the height of the breaking wave to the wavelength, k is the wave number and d is the water depth. The steepness factor here takes the value 0.12.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011418288.6A CN112417705A (en) | 2020-12-07 | 2020-12-07 | Island sea area extreme value wave calculation method based on environment contour line |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011418288.6A CN112417705A (en) | 2020-12-07 | 2020-12-07 | Island sea area extreme value wave calculation method based on environment contour line |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112417705A true CN112417705A (en) | 2021-02-26 |
Family
ID=74775441
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011418288.6A Pending CN112417705A (en) | 2020-12-07 | 2020-12-07 | Island sea area extreme value wave calculation method based on environment contour line |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112417705A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113468726A (en) * | 2021-06-09 | 2021-10-01 | 山东电力工程咨询院有限公司 | Wave period-wave height combined distribution calculation method and system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110954071A (en) * | 2019-12-13 | 2020-04-03 | 哈尔滨工程大学 | Probability statistical analysis method for wave parameters |
| CN111611549A (en) * | 2020-05-25 | 2020-09-01 | 中国海洋大学 | Method for determining typhoon-influenced sea area design wave height based on marine environment parameters |
-
2020
- 2020-12-07 CN CN202011418288.6A patent/CN112417705A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110954071A (en) * | 2019-12-13 | 2020-04-03 | 哈尔滨工程大学 | Probability statistical analysis method for wave parameters |
| CN111611549A (en) * | 2020-05-25 | 2020-09-01 | 中国海洋大学 | Method for determining typhoon-influenced sea area design wave height based on marine environment parameters |
Non-Patent Citations (2)
| Title |
|---|
| 刘小龙 等: "岛礁海域波浪环境轮廓法的应用研究", 《第三十一届全国水动力学研讨会论文集(下册)》, pages 1068 - 1076 * |
| 李玉成 等: "缓坡上不规则波浪的破碎指标", 《水动力学研究与进展(A辑)》, no. 1, pages 21 - 27 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113468726A (en) * | 2021-06-09 | 2021-10-01 | 山东电力工程咨询院有限公司 | Wave period-wave height combined distribution calculation method and system |
| CN113468726B (en) * | 2021-06-09 | 2024-05-03 | 山东电力工程咨询院有限公司 | Wave period-wave height combined distribution calculation method and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sutherland et al. | Evaluating the performance of morphological models | |
| Gori et al. | Tropical cyclone compound flood hazard assessment: From investigating drivers to quantifying extreme water levels | |
| Angnuureng et al. | Shoreline resilience to individual storms and storm clusters on a meso-macrotidal barred beach | |
| Ranasinghe et al. | Estimating coastal recession due to sea level rise: beyond the Bruun rule | |
| Stockdon et al. | A simple model for the spatially-variable coastal response to hurricanes | |
| Stockdon et al. | Empirical parameterization of setup, swash, and runup | |
| Capobianco et al. | Nourishment design and evaluation: applicability of model concepts | |
| Zou et al. | Ensemble prediction of coastal flood risk arising from overtopping by linking meteorological, ocean, coastal and surf zone models | |
| Zhang et al. | Parameterization of nearshore wave breaker index | |
| Rutten et al. | Uncertainty in runup predictions on natural beaches using XBeach nonhydrostatic | |
| Mares-Nasarre et al. | Individual wave overtopping volumes on mound breakwaters in breaking wave conditions and gentle sea bottoms | |
| Chevalier et al. | Impact of cross-reef water fluxes on lagoon dynamics: a simple parameterization for coral lagoon circulation model, with application to the Ouano Lagoon, New Caledonia | |
| CN112417705A (en) | Island sea area extreme value wave calculation method based on environment contour line | |
| Feng et al. | Investigating coastal backwater effects and flooding in the coastal zone using a global river transport model on an unstructured mesh | |
| Li et al. | Tide-modulated wave characteristics and breaking regimes in the intertidal zone of a dissipative beach | |
| Lashley et al. | Comparison of deep-water-parameter-based wave overtopping with wirewall field measurements and social media reports at Crosby (UK) | |
| CN111597506B (en) | Prediction method for near-shore wave breaking parameters and wave height | |
| Tuan et al. | Distribution of wave heights on steep submerged reefs | |
| CN108009127A (en) | The analysis method that a kind of typhoon characterization factor influences sea wave height | |
| Sutherland et al. | Coastal defence vulnerability 2075 | |
| Verdonk | The impact of overtopping on the failure probability of slipped river dikes | |
| Gori et al. | Tropical cyclone compound flood hazard assessment: from investigating dependence to quantifying impacts | |
| Pelikán et al. | Hindcast of wind driven wave heights in water reservoirs. | |
| CN112326194B (en) | Wave energy flux-based measurement and analysis method for calculating actually-measured wave friction coefficient | |
| CN111859707B (en) | Pipeline diameter determining method, urban waterlogging early warning method and device |
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 |