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 PDF

Info

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
Application number
CN201910845506.5A
Other languages
Chinese (zh)
Other versions
CN110580392B (en
Inventor
邹丽
姜胜超
张九鸣
裴玉国
王振
宗智
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201910845506.5A priority Critical patent/CN110580392B/en
Publication of CN110580392A publication Critical patent/CN110580392A/en
Application granted granted Critical
Publication of CN110580392B publication Critical patent/CN110580392B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Polynomial spectrum fitting method for representing near-island reef shallow water wave energy characteristics
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.
CN201910845506.5A 2019-09-06 2019-09-06 Polynomial spectrum fitting method for representing near-island reef shallow water wave energy characteristics Active CN110580392B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
邹丽 等: "理论风浪谱修正方法及数值验证", 《哈尔滨工程大学学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
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