CN110580392A - Polynomial spectrum fitting method for representing near-island reef shallow water wave energy characteristics - Google Patents
Polynomial spectrum fitting method for representing near-island reef shallow water wave energy characteristics Download PDFInfo
- Publication number
- CN110580392A CN110580392A CN201910845506.5A CN201910845506A CN110580392A CN 110580392 A CN110580392 A CN 110580392A CN 201910845506 A CN201910845506 A CN 201910845506A CN 110580392 A CN110580392 A CN 110580392A
- Authority
- CN
- China
- Prior art keywords
- peak
- spectrum
- wave
- wave spectrum
- polynomial
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 127
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000002902 bimodal effect Effects 0.000 claims abstract description 12
- 230000003595 spectral effect Effects 0.000 claims description 10
- 238000012216 screening Methods 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 abstract description 7
- 238000010276 construction Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Measurement Of Radiation (AREA)
Abstract
A polynomial spectrum fitting method for representing the energy characteristics of shallow water waves near an island reef belongs to the field of ocean engineering. The method provides corresponding polynomial spectrum characterization forms aiming at a single-peak spectrum and a double-peak spectrum respectively according to the shallow water wave spectrum characteristics. The proposed polynomial wave spectrum characterization form can express the near-island wave spectrum energy distribution condition in the measured data. Compared with the common wave spectrum form, the polynomial form provided by the invention fits the high-frequency part of the wave spectrum more in line with the characteristics of a shallow water wave spectrum, and can complete the fitting work of a bimodal spectrum.
Description
Technical Field
the invention relates to a polynomial wave spectrum fitting method for representing the energy characteristics of shallow water waves near an island reef, and belongs to the field of ocean engineering.
Background
The distribution form of the wave energy has important significance for ocean scientific research and ocean engineering construction. For a short-term wave environment within 20 minutes to 6 hours, the wave height history is generally regarded as a steady random process, and the energy spectrum density form relative to the wave frequency is obtained according to the wave height history transformation. Common spectrum types used for characterizing wave energy characteristics are JONSWAP, PM, Newman, ITTC two-parameter spectra and the like. These spectra are generally applicable to the unimodal wave spectra in deep water regions, and have been validated on a large body of measured data.
the energy characteristics of shallow water waves are different from those of deep water, and the energy characteristics of shallow water waves cannot be well expressed by utilizing the existing spectral shape or the correction form of the spectral shape. Meanwhile, a large number of bimodal spectra exist in the island-like reef wave environment, the common spectra are designed for a single-peak spectrum, and the energy distribution condition of the bimodal spectra is not fully considered.
Aiming at the island-near shallow water wave environment which is different from the deep water wave environment, the conventional common wave spectrum form cannot well express the internal energy distribution, so that a spectrum form characterization method aiming at a single-peak spectrum and a double-peak spectrum in the shallow water wave environment is required to be provided.
Disclosure of Invention
in order to solve the problems in the prior art, the invention provides a polynomial wave spectrum fitting method for representing the energy characteristics of shallow water waves near an island reef. The method provides corresponding polynomial fitting forms for the shallow water monomodal spectrum and the bimodal spectrum of the juxtasonic reef respectively. Based on wave spectrum data obtained by ocean actual measurement, the method utilizes the proposed polynomial representation form to fit actual measurement spectrum types according to a least square method on the basis of determining spectrum type classification. Fitting of a single peak spectrum and a double peak spectrum in the actually measured wave spectrum in the shallow water environment near the island is achieved, and common parameters of a target sea area are given.
In order to achieve the purpose of the invention, the technical scheme adopted by the invention patent is as follows: a polynomial spectrum fitting method for characterizing the energy characteristics of shallow water waves near an island reef comprises the following steps:
a. And acquiring the island-near shallow water wave spectrum of the target sea area in different time periods by using the wave rider type ocean buoy.
b. And determining the positions of the main peak value and the secondary peak value of the actually measured wave spectrum.
c. And determining the actually measured spectrum to be in a single-peak distribution form or a double-peak distribution form according to the relation between the main peak value and the secondary peak value.
d. And determining a corresponding polynomial wave spectrum fitting form according to the actually measured wave spectrum form.
e. and fitting the actually measured spectrum shape by using the proposed polynomial wave spectrum according to a least square method.
f. obtaining parameters to be fitted in the polynomial wave spectrum, completing the expression of the wave energy form of the target sea area at a specific time, and determining common parameters.
The method comprises the following specific steps:
In step 101, actually measuring energy spectrum data S (f) of waves in a target sea area by using a wave knight type ocean buoy;
In step 102, a main peak value and a secondary peak value of the measured wave spectrum are determined; wherein the main peak position peak1At the maximum of the wave spectral energy, i.e. S (peak)1)=max[S(f)](ii) a Peak position peak2Is represented by the formula [ S (f). ((f-peak))1)2]Determining the position of the maximum value of (2);
In step 103, selecting a dual-peak distribution form in the actually measured wave spectrum according to the position relationship between the main peak and the secondary peak, wherein the screening standard is as follows:
a) The main peak value and the secondary peak value are both in the range of 0.05-0.25 Hz;
b) The interval between the main peak and the secondary peak is greater than or equal to 0.05Hz, i.e. | peak1-peak2|≥0.05Hz;
c) The wave spectral energy density of the secondary peak is greater than 30% of the primary peak, i.e. S (peak)2)≥0.3S(peak1);
d) The valleys of the wave spectrum between the major and minor peaks are less than 2/3 for the minor peak, i.e.
When the actually measured wave spectrum simultaneously meets the conditions, determining that the wave spectrum is a wave spectrum in a double-peak value form, otherwise, determining that the wave spectrum is a wave spectrum in a single-peak value form;
In step 104, respectively proposing polynomial wave spectrum forms of a single-peak distribution form and a double-peak distribution form of the island-near reef shallow water wave spectrum;
The polynomial form for the single-peak wave spectrum is:Wherein A, C, a and b are parameters to be fitted;
The polynomial form for the bimodal waveform spectrum is:wherein A, B, C, D, a and B are parameters to be fitted;
In step 105, fitting a typical measured wave spectrum by using a least square method according to the corresponding polynomial wave spectrum expression form;
In step 106, the common parameters in the polynomial of the unimodal wave spectrum and the polynomial of the bimodal wave spectrum of the target sea area are obtained.
The invention has the following beneficial effects:
(1) the method can be used for representing a single-peak spectrum form in the island-like reef shallow water environment, expanding the method for expressing the wave energy by the wave spectrum from a deep water area to a shallow water area, and determining common parameters.
(2) The method can express a double-peak spectrum form which often appears in the island-like reef shallow water environment, and expand a single-peak energy distribution form to a double-peak form, thereby providing theoretical basis and guarantee for marine scientific research and marine engineering construction.
(3) Compared with the common wave spectrum form, the polynomial form provided by the invention is more suitable for the high-frequency part of the wave spectrum to be fitted with the characteristics of the shallow water wave spectrum.
Drawings
FIG. 1 is a flow chart of a polynomial spectrum fitting method for characterizing shallow water wave energy of a juxtapose reef.
FIG. 2 is a typical pattern of a measured single peak wave spectrum for a wave rider buoy.
FIG. 3 is a typical pattern of a measured bimodal wave spectrum for a wave rider buoy.
FIG. 4 is a typical fitting result of a polynomial of the unimodal spectral type.
FIG. 5 is a typical fitting result of a bimodal spectral pattern polynomial.
FIG. 6 shows the comparison of the polynomial spectrum pattern with the measured spectrum pattern under the common parameters.
Detailed Description
the technical scheme of the invention is further explained in detail by combining the attached drawings:
The method is suitable for the characterization method of the near-island reef shallow water wave energy distribution form, and is realized by polynomial fitting based on wave spectrum data obtained by actual measurement.
Fig. 1 is a flow chart of a wave spectrum polynomial characterization method suitable for shallow water near an island, and this embodiment takes a wave spectrum measured by a wave knight-type ocean buoy as target information, which illustrates that the method mainly includes the following steps: in step 101, energy spectrum data s (f) of the wave is actually obtained in the target sea area by using a wave knight type ocean buoy. In this embodiment, the measured wave spectrum s (f) is used as the target information.
In step 102, the primary and secondary peaks are first determined for the measured wave spectrum. Wherein the main peak position peak1At the maximum of the wave spectral energy, i.e. S (peak)1)=max[S(f)]. Peak position peak2is represented by the formula [ S (f). ((f-peak))1)2]is determined by the maximum position of (2).
In step 103, selecting a dual-peak distribution form in the actually measured wave spectrum according to the position relationship between the main peak and the secondary peak, wherein the screening standard is as follows:
e) The major peak and the minor peak are both in the range of 0.05-0.25 Hz.
f) The interval between the main peak and the secondary peak is greater than or equal to 0.05Hz, i.e. | peak1-peak2|≥0.05Hz。
g) the wave spectral energy density of the secondary peak is greater than 30% of the primary peak, i.e. S (peak)2)≥0.3S(peak1)。
h) The valleys of the wave spectrum between the major and minor peaks are less than 2/3 for the minor peak, i.e.
And when the actually measured wave spectrum simultaneously meets the conditions, determining that the wave spectrum is in a double-peak form, otherwise, determining that the wave spectrum is in a single-peak form. A typical observed single peak wave spectrum pattern is shown in FIG. 2, and a typical observed double peak wave spectrum pattern is shown in FIG. 3.
in step 104, polynomial wave spectrum forms for the island-near reef shallow water wave spectrum unimodal distribution form and the double-peaked distribution form are proposed respectively. The polynomial form for the single-peak wave spectrum is:Wherein A, C, a and b are parameters to be fitted; the polynomial form for the bimodal waveform spectrum is:
wherein A, B, C, D, a and B are parameters to be fitted. In step 105, according to the corresponding polynomial wave spectrum expression form, fitting a typical measured wave spectrum by using a least square method, wherein a fitting result for a single-peak wave spectrum form is shown in fig. 4, and a fitting result for a double-peak wave spectrum form is shown in fig. 5.
In step 106, the common parameters for obtaining the single-peak polynomial wave spectrum form of the target sea area are: a is 0.008 and C is 2.5 × 10-9Where a is 7 and b is 9.5. Common parameters for the bimodal polynomial wave spectrum form are: a is 0.008 and B is 1.4X 10-9,C=2.5×10-9,D=0.8×10-9,a=7,b=9.5。
the common parameters are used for expressing the form of the target sea wave spectrum, and the comparison result with the measured data is shown in fig. 6. As can be seen from the figure, the wave spectrum obtained in the form of the single-peak polynomial wave spectrum is basically consistent with the wave spectrum obtained from actually measured data in waveform, and the fitting accuracy is better than that of common PM spectrum and Newman spectrum, so that the method can be used for preparing fitting single-peak spectrum and double-peak spectrum. Compared with the common wave spectrum form, the polynomial form provided by the invention is more suitable for the characteristic of the shallow water wave spectrum by fitting the high-frequency part of the wave spectrum.
Based on the near-island reef shallow water wave energy distribution form, a polynomial wave spectrum characterization method is provided based on the wave spectrum actually measured by the wave rider type ocean buoy. The method can meet the wave characterization of the unimodal wave energy distribution and the double-peaked wave energy distribution aiming at the shallow water wave environment, realizes the determination of the common parameters of the target sea area unimodal polynomial, provides reference for ocean scientific research and provides theoretical guarantee for ocean engineering construction.
Claims (1)
1. a polynomial spectrum fitting method for characterizing shallow water wave energy characteristics of a juxtasonic reef is characterized by comprising the following steps of:
(1) Actually measuring and obtaining energy spectrum data S (f) of waves in a target sea area by utilizing a wave knight type ocean buoy;
(2) Firstly, determining a main peak value and a secondary peak value of the actually measured wave spectrum; wherein the main peak position peak1at the maximum of the wave spectral energy, i.e. S (peak)1)=max[S(f)](ii) a Peak position peak2Is represented by the formula [ S (f). ((f-peak))1)2]determining the position of the maximum value of (2);
(3) Selecting a double-peak distribution form in the actually measured wave spectrum according to the position relation of the main peak and the secondary peak, wherein the screening standard is as follows:
a) The main peak value and the secondary peak value are both in the range of 0.05-0.25 Hz;
b) The interval between the main peak and the secondary peak is greater than or equal to 0.05Hz, i.e. | peak1-peak2|≥0.05Hz;
c) The wave spectral energy density of the secondary peak is greater than 30% of the primary peak, i.e. S (peak)2)≥0.3S(peak1);
d) the valleys of the wave spectrum between the major and minor peaks are less than 2/3 for the minor peak, i.e.
When the actually measured wave spectrum simultaneously meets the conditions, determining that the wave spectrum is a wave spectrum in a double-peak value form, otherwise, determining that the wave spectrum is a wave spectrum in a single-peak value form;
(4) respectively proposing a polynomial wave spectrum form aiming at a single-peak distribution form and a double-peak distribution form of the island-near reef shallow water wave spectrum;
The polynomial form for the single-peak wave spectrum is:Wherein A, C, a and b are parameters to be fitted;
the polynomial form for the bimodal waveform spectrum is:Wherein A, B, C, D, a and B are parameters to be fitted;
(5) fitting the typical actually measured wave spectrum by using a least square method according to the corresponding polynomial wave spectrum expression form;
(6) And obtaining the common parameters in the polynomial of the unimodal wave spectrum and the polynomial of the bimodal wave spectrum of the target sea area.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910845506.5A CN110580392B (en) | 2019-09-06 | 2019-09-06 | Polynomial spectrum fitting method for representing near-island reef shallow water wave energy characteristics |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910845506.5A CN110580392B (en) | 2019-09-06 | 2019-09-06 | Polynomial spectrum fitting method for representing near-island reef shallow water wave energy characteristics |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110580392A true CN110580392A (en) | 2019-12-17 |
CN110580392B CN110580392B (en) | 2021-10-22 |
Family
ID=68812702
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910845506.5A Active CN110580392B (en) | 2019-09-06 | 2019-09-06 | Polynomial spectrum fitting method for representing near-island reef shallow water wave energy characteristics |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110580392B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114970380A (en) * | 2021-12-13 | 2022-08-30 | 昆明理工大学 | Island factor-based sea wave simulation method for offshore island |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050228625A1 (en) * | 2004-04-12 | 2005-10-13 | Lachman Lawrence M | Method and system for modeling ocean waves |
CN106021933A (en) * | 2016-05-24 | 2016-10-12 | 国家***第海洋研究所 | Wave and turbulence nonlinear transactional analysis (TA) method |
CN106338274A (en) * | 2016-08-25 | 2017-01-18 | 厦门斯坦道科学仪器股份有限公司 | Wave measurement apparatus for correcting wave characteristic parameters based on multi-parameters, and wave measurement method thereof |
CN108318103A (en) * | 2018-04-03 | 2018-07-24 | 大连理工大学 | A kind of wave height measuring device and method can be used for practical marine site |
CN109782287A (en) * | 2019-03-12 | 2019-05-21 | 国家***北海预报中心 | A method of it is surged wave height by high-frequency ground wave radar inverting |
-
2019
- 2019-09-06 CN CN201910845506.5A patent/CN110580392B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050228625A1 (en) * | 2004-04-12 | 2005-10-13 | Lachman Lawrence M | Method and system for modeling ocean waves |
CN106021933A (en) * | 2016-05-24 | 2016-10-12 | 国家***第海洋研究所 | Wave and turbulence nonlinear transactional analysis (TA) method |
CN106338274A (en) * | 2016-08-25 | 2017-01-18 | 厦门斯坦道科学仪器股份有限公司 | Wave measurement apparatus for correcting wave characteristic parameters based on multi-parameters, and wave measurement method thereof |
CN108318103A (en) * | 2018-04-03 | 2018-07-24 | 大连理工大学 | A kind of wave height measuring device and method can be used for practical marine site |
CN109782287A (en) * | 2019-03-12 | 2019-05-21 | 国家***北海预报中心 | A method of it is surged wave height by high-frequency ground wave radar inverting |
Non-Patent Citations (1)
Title |
---|
邹丽 等: "理论风浪谱修正方法及数值验证", 《哈尔滨工程大学学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114970380A (en) * | 2021-12-13 | 2022-08-30 | 昆明理工大学 | Island factor-based sea wave simulation method for offshore island |
Also Published As
Publication number | Publication date |
---|---|
CN110580392B (en) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110580392B (en) | Polynomial spectrum fitting method for representing near-island reef shallow water wave energy characteristics | |
CN109143193B (en) | High-frequency ground wave radar wave field inversion method based on second-order peak | |
CN110177934A (en) | The noise analysis and control of wind turbine | |
CN103076604B (en) | Method for measuring distance of low-frequency underwater sound pulse signal on basis of frequency dispersion features | |
CN112763180B (en) | Method for quickly predicting aerodynamic noise of pantograph model of high-speed train in acoustic wind tunnel | |
CN105301572B (en) | Multi-peak low sidelobe phase-coded signal and mismatched filter co-design method | |
CN114091377B (en) | Method for calculating wake flow wind speed of dynamic double-Gaussian wind turbine based on spatial variation | |
CN107340499A (en) | The sane low-sidelobe beam forming method rebuild based on covariance matrix | |
CN103971027A (en) | Rice automatic mapping method based on ratio index of water body and vegetation index changes | |
CN105783843A (en) | Sensor noise signal inhibition method in three-point-method rotation error separation process | |
Traverso et al. | Spectral analysis of the underwater acoustic noise radiated by ships with controllable pitch propellers | |
CN110109087B (en) | Sonar irregular detection range display method and system | |
CN103559373A (en) | Method for optimizing technological parameters of cold ring rolling on basis of fuzzy grey relation | |
CN116680882A (en) | Ocean depth benchmark construction method and device based on tide splitting recombination and tide checking constraint | |
CN114252509B (en) | Three-section type locking landslide precursor identification method based on acoustic emission signals | |
CN115408860A (en) | Abnormal value monitoring and correcting method for wind speed-power data of wind power plant | |
CN104062635A (en) | Sea clutter image sea wave principal wave period retrieval method based on minimum interval thought | |
CN110837085B (en) | Fluctuation index calculation method for underwater target discrimination | |
CN112162235B (en) | Smooth segmented stochastic resonance enhanced acoustic vector signal orientation method | |
CN113049220B (en) | Method for identifying wave and surge component in mixed waves | |
Johansson et al. | Full scale measurement of underwater radiated noise from a coastal tanker | |
CN112434467A (en) | Low-target-strength shape design method for underwater rotating body target | |
Jacobsen et al. | Nearshore versus offshore: comparative cost and competitive advantages | |
CN106143850A (en) | Propeller blade section and application process thereof | |
CN118520690A (en) | Calculation method for workability of offshore construction of engineering ship |
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 |