Asymmetric baffle float wave energy device broadband wave attenuation method and system
Patent Information
- Application Number
- CN202610839647.6
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2026-06-11
- Publication Date
- 2026-09-01
- Estimated Expiration
- 2046-06-11
AI Technical Summary
[0004]本申请提供了一种非对称挡板浮子波浪能装置宽频消波方法及系统,解决了现有非对称挡板浮子波浪能装置因无法区分背浪侧透射波物理成因而导致消波控制方向错误的问题,以及阻尼增大模式触发后散射透射分量补偿性增大缺乏专项抑制手段的问题,在宽频带入射波条件下稳定维持了较低的背浪侧透射系数
[0010]Using the relationship between the component ratio and a preset threshold as the switching criterion, when the component ratio is not lower than the preset threshold, the damping coefficient of the power output device is increased based on the positive difference value; when the component ratio is lower than the preset threshold, the damping coefficient of the power output device is decreased based on the positive difference value. This forms a dual-mode reverse switching control law with completely opposite directions. The technical contribution of this control law is that increasing the damping can suppress the buoy's heave amplitude from the radiation source, thereby reducing the radiation transmission component; decreasing the damping can amplify the buoy's heave amplitude, enhance the dynamic interception of the incident wave by the asymmetric baffle, thereby reducing the scattered transmission component. The two intervention directions are automatically switched in real time driven by the component ratio, overcoming the fundamental defect of the existing unidirectional fixed damping adjustment strategy, which cannot identify the cause of wave suppression failure and thus leads to incorrect control direction. Furthermore, using the increment of the scattered transmission amplitude relative to the reference time as the criterion, when the increment exceeds the scattering compensation threshold, the back wave side guide plate is deployed to the preset working angle. The guide plate is used to compensate for the increase in scattering components caused by the decrease in float heave amplitude under the special compensation damping mode of the refracted and deflected scattered transmission waves. After the scattering transmission amplitude returns to the compensation range, the guide plate is retracted, forming a functional division closed loop between active PTO control and passive guide structure. The two types of transmission components are controlled synchronously in a wide frequency band, realizing stable broadband wave suppression of the asymmetric baffle float wave energy device under different incident wave frequency conditions.
Smart Images

Figure CN122382932B_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of wave energy utilization technology, and in particular to a broadband wave suppression method and system for an asymmetric baffle float wave energy device. Background Technology
[0002] The integrated design of floating breakwaters and oscillating buoy-type wave energy conversion devices is an important technological direction in the field of nearshore protection and wave energy utilization. In existing technologies, asymmetric triangular baffle buoy wave energy devices replace rectangular symmetrical floats with triangular asymmetric cross-section floats. Utilizing the geometric asymmetry of the floats, they alter the wave radiation and scattering characteristics on the wave-facing and wave-avoiding sides. While maintaining high wave energy capture efficiency, they broaden the effective wave-damping frequency band. The transmission coefficient under shortwave conditions is significantly lower than that of symmetrical rectangular float devices, making it one of the mainstream technological routes that combines power generation and protection functions.
[0003] However, existing asymmetric baffle float wave energy devices have fundamental flaws in wave suppression control. In current technology, the total transmitted wave surface elevation on the back wave side is monitored as a single entity, with the total transmission coefficient as the sole wave suppression evaluation index. The adjustment direction of the damping coefficient of the power output device is fixed in one direction, either always pursuing maximum energy capture power or always pursuing minimum total transmission coefficient. The fundamental flaw of this approach is that the total transmitted wave on the back wave side is physically composed of two components with different origins: scattered transmission components and radiated transmission components. Their responses to damping adjustment are completely opposite—increasing damping suppresses the radiated transmission component but weakens the dynamic shielding effect of the float, leading to an increase in the scattered transmission component; decreasing damping enhances the dynamic shielding effect of the float but amplifies the radiated transmission component. Current technology cannot distinguish the real-time contributions of these two components. Therefore, regardless of the unidirectional fixed damping strategy adopted, when the incident wave frequency deviates from the device's design operating point, control direction errors inevitably occur, leading to rapid degradation of broadband wave suppression performance. Summary of the Invention
[0004] This application provides a broadband wave suppression method and system for an asymmetric baffle float wave energy device, which solves the problem that existing asymmetric baffle float wave energy devices cannot distinguish the physical composition of the back wave side transmitted waves, resulting in incorrect wave suppression control direction, and the lack of specific suppression methods for the compensatory increase of the scattered transmission component after the damping increase mode is triggered. It stably maintains a low back wave side transmission coefficient under broadband incident wave conditions.
[0005] In a first aspect, this application provides a broadband wave attenuation method for an asymmetric baffle float wave energy device, the broadband wave attenuation method for the asymmetric baffle float wave energy device comprising: Step S1: Collect the wave surface elevation on the wave-facing side of the asymmetric baffle float, the total transmitted wave surface elevation on the wave-receiving side, and the float heave velocity to obtain the raw transmitted wave data set. Step S2: The buoy heave velocity in the original transmitted wave data set is phase-shifted by a quarter period to obtain the radiation transmission reference signal. The total transmitted wave surface elevation on the back wave side is subtracted from the radiation transmission reference signal and then the scattering transmission amplitude is obtained by moving root mean square operation. The radiation transmission amplitude is obtained by performing moving root mean square operation on the radiation transmission reference signal. The radiation transmission amplitude is divided by the scattering transmission amplitude to obtain the component ratio. Step S3: Using the relationship between the component ratio and the preset threshold as the switching criterion, when the component ratio is not lower than the preset threshold, the damping coefficient of the power output device is increased based on the positive difference between the component ratio and the preset threshold as the increment; when the component ratio is lower than the preset threshold, the damping coefficient of the power output device is decreased based on the positive difference between the preset threshold and the component ratio as the decrement, thereby obtaining a wave suppression control command. Step S4: Using the increment of the scattering and transmission amplitude relative to the reference time as a criterion, when the increment exceeds the scattering compensation threshold, the back wave side guide plate is deployed to a preset working angle, and the guide plate is retracted after the scattering and transmission amplitude returns to the compensation range.
[0006] Secondly, this application provides a broadband wave suppression system for an asymmetric baffle float wave energy device, the broadband wave suppression system comprising: The acquisition module is used to acquire the wave surface elevation on the wave-facing side of the asymmetric baffle float, the total transmitted wave surface elevation on the wave-receiving side, and the float heave velocity to obtain the raw transmitted wave data set. The calculation module is used to obtain a radiation transmission reference signal by shifting the float heave velocity in the original transmitted wave data group by a quarter-cycle phase, to obtain a scattering transmission amplitude by subtracting the total transmitted wave surface elevation on the back wave side from the radiation transmission reference signal and then performing a sliding root mean square operation on the radiation transmission reference signal, to obtain a radiation transmission amplitude, and to obtain a component ratio by dividing the radiation transmission amplitude by the scattering transmission amplitude. The output module is used to use the relationship between the component ratio and a preset threshold as a switching criterion. When the component ratio is not lower than the preset threshold, the damping coefficient of the power output device is increased based on the positive difference between the component ratio and the preset threshold as an incremental reference. When the component ratio is lower than the preset threshold, the damping coefficient of the power output device is decreased based on the positive difference between the preset threshold and the component ratio as a decremental reference, thereby obtaining a wave suppression control command. The retraction module is used to determine the scattering and transmission amplitude relative to the reference time. When the increment exceeds the scattering compensation threshold, the back wave side guide plate is deployed to a preset working angle, and the guide plate is retracted after the scattering and transmission amplitude returns to the compensation range.
[0007] Thirdly, a broadband wave-damping device for an asymmetric baffle float wave energy device is provided, comprising: a memory and at least one processor, wherein the memory stores instructions; the at least one processor invokes the instructions in the memory to cause the broadband wave-damping device for an asymmetric baffle float wave energy device to perform the aforementioned broadband wave-damping method for an asymmetric baffle float wave energy device.
[0008] Fourthly, a computer-readable storage medium is provided, wherein instructions are stored therein, which, when executed on a computer, cause the computer to perform the aforementioned broadband wave attenuation method of the asymmetric baffle float wave energy device.
[0009] In the technical solution provided in this application, the present invention constructs a raw data set of transmitted waves by arranging wave height sensors on the wave-facing and wave-avoiding sides of an asymmetric baffle float and simultaneously collecting the wave surface elevation on the wave-facing side, the total transmitted wave surface elevation on the wave-avoiding side, and the float's heave velocity. This provides a complete input under a unified time reference for subsequent component separation. Based on this, the float's heave velocity is shifted by a quarter-cycle phase to obtain a radiation transmission reference signal. Then, the difference between the total transmitted wave surface elevation on the wave-avoiding side and the radiation transmission reference signal is calculated, followed by a moving root mean square operation to obtain the scattered transmission amplitude and the radiation transmission amplitude, respectively. These two values are then divided to obtain the component ratio. The core contribution of this component separation algorithm lies in the fact that the total transmitted wave surface elevation on the wave-avoiding side is physically composed of the superposition of scattered transmission and radiation transmission components. Existing technologies treat these two components as a whole, failing to distinguish the physical causes of wave-damping failure. The present invention, however, utilizes the fixed quarter-cycle phase relationship between the radiation transmission wave and the float's heave velocity to, for the first time, separate the two types of components from the aliased signal in real time into independent quantized variables, enabling subsequent control strategies to implement differentiated interventions for different physical causes.
[0010] Using the relationship between the component ratio and a preset threshold as the switching criterion, when the component ratio is not lower than the preset threshold, the damping coefficient of the power output device is increased based on the positive difference value; when the component ratio is lower than the preset threshold, the damping coefficient of the power output device is decreased based on the positive difference value. This forms a dual-mode reverse switching control law with completely opposite directions. The technical contribution of this control law is that increasing the damping can suppress the buoy's heave amplitude from the radiation source, thereby reducing the radiation transmission component; decreasing the damping can amplify the buoy's heave amplitude, enhance the dynamic interception of the incident wave by the asymmetric baffle, thereby reducing the scattered transmission component. The two intervention directions are automatically switched in real time driven by the component ratio, overcoming the fundamental defect of the existing unidirectional fixed damping adjustment strategy, which cannot identify the cause of wave suppression failure and thus leads to incorrect control direction. Furthermore, using the increment of the scattered transmission amplitude relative to the reference time as the criterion, when the increment exceeds the scattering compensation threshold, the back wave side guide plate is deployed to the preset working angle. The guide plate is used to compensate for the increase in scattering components caused by the decrease in float heave amplitude under the special compensation damping mode of the refracted and deflected scattered transmission waves. After the scattering transmission amplitude returns to the compensation range, the guide plate is retracted, forming a functional division closed loop between active PTO control and passive guide structure. The two types of transmission components are controlled synchronously in a wide frequency band, realizing stable broadband wave suppression of the asymmetric baffle float wave energy device under different incident wave frequency conditions. Attached Figure Description
[0011] To more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0012] Figure 1 This is a schematic diagram of an embodiment of the broadband wave attenuation method of the asymmetric baffle float wave energy device in this application. Figure 2 This is a schematic diagram of the dynamic response of scattering and transmission amplitude and the triggering control process of the back wave side guide vane in the embodiments of this application. Detailed Implementation
[0013] This application provides a broadband wave attenuation method and system for an asymmetric baffle float wave energy device. The terms "first," "second," "third," "fourth," etc. (if present) in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments described herein can be implemented in a sequence other than that illustrated or described herein. Furthermore, the terms "comprising" or "having" and any variations thereof are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or device that includes a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or devices.
[0014] For ease of understanding, the specific process of the embodiments of this application is described below. Please refer to [link / reference]. Figure 1 One embodiment of the broadband wave attenuation method for the asymmetric baffle float wave energy device in this application includes: Step S1: Collect the wave surface elevation on the wave-facing side of the asymmetric baffle float, the total transmitted wave surface elevation on the wave-receiving side, and the float heave velocity to obtain the raw transmitted wave data set. Specifically, the raw transmitted wave data set serves as the unified data source for all subsequent calculations in this method. It comprises three types of physical quantities: the wavefront elevation on the wave-facing side, the total transmitted wavefront elevation on the wave-avoiding side, and the buoy's heave velocity. The wavefront elevation on the wave-facing side is collected by a wave height sensor positioned at a depth of 0.05 times the water surface on the wave-facing side of the asymmetric baffle buoy. The total transmitted wavefront elevation on the wave-avoiding side is collected by wave height sensors symmetrically positioned at the same depth on the wave-avoiding side. Both sensors are synchronously triggered using the same sampling frequency to ensure time alignment. The buoy's heave velocity is not directly measured but is obtained by collecting the heave acceleration signal from a vertical accelerometer installed at the buoy's center of mass. The initial value for integration, zero heave velocity under static water equilibrium, is used as the initial value for first-order numerical integration. The initial value is chosen based on the physical fact that the buoy has no velocity under static water conditions. These three types of signals are then timestamped and merged into the raw transmitted wave data set, ensuring that subsequent steps use signals from the same time reference.
[0015] Step S2: The buoy heave velocity in the original transmitted wave data set is phase-shifted by a quarter period to obtain the radiation transmission reference signal. The total transmitted wave surface elevation on the back wave side is subtracted from the radiation transmission reference signal and then the scattering transmission amplitude is obtained by the moving root mean square operation. The radiation transmission amplitude is obtained by performing the moving root mean square operation on the radiation transmission reference signal. The radiation transmission amplitude is divided by the scattering transmission amplitude to obtain the component ratio. Specifically, the radiation transmission reference signal is a quantified proxy signal for the radiation component in the total transmitted wavefront elevation on the back wave side. Its physical basis lies in the following: According to linear potential flow theory, the wavefront elevation radiated towards the back wave side by the heave motion of the asymmetric baffle buoy leads the buoy's heave velocity by a quarter of a period. Therefore, by forward-shifting the buoy's heave velocity time series by a quarter of the incident wave's main period to the corresponding number of sampling points, and then multiplying by the radiation wave transfer coefficient, the radiation transmission reference signal is obtained. The radiation wave transfer coefficient is calculated by dividing the radiation damping coefficient by the product of seawater density, gravitational acceleration, and the buoy's waterline area. The radiation damping coefficient is determined before device installation through a water tank calibration test with the buoy fixed and the power output device damped to zero. The component ratio is a dimensionless number obtained by dividing the radiation transmission amplitude by the scattering transmission amplitude. When the scattering transmission amplitude is lower than the preset minimum value, the quotient of the previous effective moment is used to replace the current calculation result. This protection logic prevents division by zero anomalies in near-still water conditions. The preset minimum value is 0.001 meters, which is lower than the minimum effective wave surface disturbance that can be measured in actual sea conditions, ensuring that the replacement logic is not triggered under normal wave conditions.
[0016] Step S3: Using the relationship between the component ratio and the preset threshold as the switching criterion, when the component ratio is not lower than the preset threshold, the damping coefficient of the power output device is increased based on the positive difference between the component ratio and the preset threshold; when the component ratio is lower than the preset threshold, the damping coefficient of the power output device is decreased based on the positive difference between the preset threshold and the component ratio, thus obtaining the wave suppression control command. Specifically, the preset threshold is the quantitative boundary that distinguishes between radiation transmission-dominated and scattering transmission-dominated operating conditions. It is determined during the calibration phase after device installation and before normal operation: with the float fixed and the power output device damping coefficient set to zero, the component ratio sequence is continuously collected for no less than twenty incident wave cycles, and the arithmetic mean of the component ratio is taken to obtain the baseline mean of the component ratio. This is then multiplied by a preset margin coefficient to obtain the preset threshold. The preset margin coefficient is a fixed value within the range of greater than one and not greater than 1.5. The value is greater than one because the component ratio under calibration conditions reflects the natural ratio of the two types of components under undamped conditions. During normal operation, the intervention of the power output device damping will suppress the radiation component to some extent, and a margin is required to avoid false triggering of the radiation-dominated mode. The value is not greater than 1.5 to prevent the threshold from being too high, which would lead to a sluggish response in the radiation-dominated mode. The core innovation of this step lies in the reverse switching logic of the switching criterion: when the component ratio is not lower than the preset threshold, the damping is increased to suppress the radiation component; when the component ratio is lower than the preset threshold, the damping is decreased to enhance the dynamic shielding of the float. The control responses in the two directions are completely opposite, which is fundamentally different from the unidirectional fixed damping adjustment strategy in the existing technology.
[0017] Step S4: Using the increment of the scattered transmission amplitude relative to the reference time as the criterion, when the increment exceeds the scattering compensation threshold, the back wave side guide plate is deployed to the preset working angle, and the guide plate is retracted after the scattered transmission amplitude returns to the compensation range.
[0018] Specifically, the scattering compensation threshold is a quantitative criterion for triggering the deployment of the deflector on the back wave side. It is calculated as the product of the root mean square (RMS) elevation of the wavefront elevation on the front wave side within the current main cycle of the incident wave and a preset compensation ratio coefficient. The preset compensation ratio coefficient is 0.05, meaning that deployment is triggered when the increase in the scattering transmission amplitude exceeds five percent of the current RMS height of the incident wave. This ratio reflects the minimum detectable change in the scattering transmission component that results in a significant compensatory increase. The preset operating angle is determined during the pre-installation calibration test by measuring the main propagation elevation angle of the scattering transmission wave using a two-dimensional wave pressure sensor array on the back wave side and taking the remainder of the angle. This ensures that the deflector's normal direction is orthogonal to the propagation direction of the scattering transmission wave, achieving maximum refraction deflection. The scattering transmission amplitude at the reference time is the arithmetic mean of the scattering transmission amplitude within the main cycle of the incident wave before the first switch of the wave suppression control command to the radiation-dominant indicator. Using the arithmetic mean rather than the instantaneous value as the reference aims to eliminate interference from signal noise within a single cycle. The preset fallback margin in the deflector retraction criterion is 2% of the scattering and transmission amplitude at the reference time, to prevent the deflector from frequently extending and retracting when the scattering and transmission amplitude fluctuates near the compensation threshold.
[0019] In one specific embodiment, step S1 includes: A wave height sensor is arranged on the wave-facing side and the wave-back side of the asymmetric baffle float. The output signals of the wave height sensor on the wave-facing side and the wave height sensor on the wave-back side are synchronously collected based on a preset sampling frequency to obtain the wave surface elevation time series on the wave-facing side and the total transmitted wave surface elevation time series on the wave-back side. The heave acceleration signal output by the vertical accelerometer at the center of mass of the float is integrated with the initial value of zero heave velocity under the static water level equilibrium state, and a first-order numerical integration is performed to obtain the time series of heave velocity of the float. With the buoy fixed and the damping coefficient of the power output device set to zero, the cross-correlation operation is performed on the wavefront elevation time series on the wave-facing side and the total transmitted wavefront elevation time series on the wave-returning side to obtain the fixed phase difference between the two, which is denoted as the scattering phase calibration value. The time series of wavefront elevation on the wave-facing side, the time series of total transmitted wavefront elevation on the wave-receiving side, the time series of float heave velocity, and the scattering phase calibration values were timestamped and then merged to obtain the original transmitted wave data set.
[0020] Specifically, the preset sampling frequency is set to 100Hz. The actual main frequency range of nearshore waves is 0.05Hz to 5Hz. According to the Nyquist sampling theorem, the sampling frequency should not be less than twice the highest signal frequency to avoid spectral aliasing. The wave front elevation time series mentioned in this method needs to participate in multiple steps of calculation, such as quarter-cycle phase shift, moving root mean square, cross-correlation, and fast Fourier transform. In order to ensure the waveform restoration accuracy and the time resolution of phase shift, 100Hz, which is twenty times the highest frequency of 5Hz, is taken as the preset sampling frequency in engineering. This value is higher than the minimum sampling frequency requirement specified by the Nyquist theorem.
[0021] The wave height sensor on the wave-facing side is positioned at a depth of 0.05 times the water surface on the wave-facing side of the asymmetric baffle buoy, while the wave height sensor on the back side is symmetrically positioned at the same depth on the back side. Both sensors are triggered by the same clock source to ensure strict alignment of the wave surface elevation time series on the wave-facing side with the total transmitted wave surface elevation time series on the back side on the time axis. A vertical accelerometer is installed at the center of mass of the buoy to eliminate interference from the buoy's rotational motion on the vertical acceleration measurement. The first-order numerical integration is discretized using the rectangular integration method, with an integration step size of 0.01 seconds between sampling intervals. The initial integration value is taken as zero heave velocity under static water equilibrium conditions. The physical basis for this initial value is that before the device is started, the buoy is stationary on the water surface, and its vertical velocity is zero. Using this as the initial value avoids drift of the integration constant.
[0022] The scattering phase calibration value is extracted through cross-correlation calculation: the time series of wavefront elevation on the wave-facing side and the time series of total transmitted wavefront elevation on the wave-reverse side are progressively shifted within a data window of no less than 20 incident wave main periods, and the correlation coefficient between the two series at each shift step is calculated; the correlation coefficient is calculated by subtracting the arithmetic mean of the corresponding sampling points of the two series at each shift step to obtain two sets of mean-reduced sequences, and then accumulating the product of the corresponding points of the two sets of mean-reduced sequences within the data window length to obtain the cumulative value. The product value is obtained by summing the squares of the two mean-reduced sequences within the data window length and taking the positive square root to obtain two sets of standardized terms. The summed product value is divided by the product of the two sets of standardized terms to obtain the correlation coefficient at the translation step size. The translation step size corresponding to the maximum value of the correlation coefficient is recorded as the maximum correlation translation step size. The maximum correlation translation step size is multiplied by the sampling interval to obtain the delay time. The delay time is divided by the current incident wave main period and then multiplied by 360 degrees. The resulting angle value is the scattering phase calibration value.
[0023] This calibration is performed under the condition of a fixed float and a power output device with a damping coefficient of zero. The purpose is to eliminate the interference of the radiated wave component generated by the float's motion on phase extraction, ensuring that the scattered phase calibration value only reflects the fixed phase delay of the incident wave reaching the back-wave side after being geometrically scattered by the float. This value is determined by the float's geometric parameters and is a constant under the premise that the device's geometry remains unchanged. Timestamp alignment uses the start time of the wavefront elevation time series acquisition on the wave-facing side as a unified reference. After aligning the total transmitted wavefront elevation time series on the back-wave side, the float heave velocity time series, and the scattered phase calibration value to this reference time, they are merged column-wise into a raw transmitted wave data set. The scattered phase calibration value is appended to the metadata field of the data set in scalar form.
[0024] In one specific embodiment, step S2 includes: Based on the time series of the buoy heave velocity in the original transmitted wave data set, the number of sampling points corresponding to one-quarter of the current incident wave main period is used as the translation step size to perform a forward time series translation on the buoy heave velocity time series. Then, the translated sequence is scaled point by point using the radiation wave transfer coefficient to obtain the radiation transmission reference signal. The radiation wave transfer coefficient is calculated by dividing the radiation damping coefficient measured under the calibration state by the product of the seawater density, gravitational acceleration and the waterline area of the buoy. The total transmitted wave surface elevation time series on the back side of the original transmitted wave data set is subtracted point by point from the radiation transmission reference signal to obtain the scattering transmission time series. The scattering transmission amplitude is obtained by performing a sliding root mean square operation on the scattering transmission time series with the number of sampling points corresponding to the current main period of the incident wave as the window length. Using the number of sampling points corresponding to the current main period of the incident wave as the window length, perform a sliding root mean square operation on the radiation transmission reference signal to obtain the radiation transmission amplitude; Divide the radiation transmission amplitude by the scattering transmission amplitude. When the scattering transmission amplitude is lower than the preset minimum value, replace the current calculation result with the quotient value of the previous effective time to obtain the component ratio.
[0025] Specifically, the physical meaning of the radiation transmission reference signal is: the proxy signal of the wave surface elevation component radiated towards the back wave side by the heave motion generated by the asymmetric baffle float under the excitation of the incident wave. According to the linear potential flow theory, the phase of the wave surface elevation radiated towards the back wave side by the heave motion of the float leads the heave velocity of the float by one-quarter of the main period of the incident wave. Therefore, one-quarter of the number of sampling points obtained by dividing the current main period of the incident wave by the sampling interval of 0.01 seconds is used as the translation step size to perform a forward time series translation on the heave velocity time series of the float. The meaning of forward translation is to move the entire sequence in the positive direction of the time axis by the number of sampling points corresponding to the step size. After translation, the beginning of the sequence is filled with zero values. The radiation wave transmission coefficient is calculated as follows: the radiation damping coefficient measured by the forced oscillation test under calibration conditions is divided by the product of seawater density, gravitational acceleration, and the buoy's waterline area. The seawater density is taken as 1025 kg / m³, the gravitational acceleration as 9.81 m / s², and the buoy's waterline area is directly measured from the buoy's geometric drawings. The radiation damping coefficient is obtained by dividing the in-phase component of the fluid force acting on the buoy after applying forced heave motion of known amplitude and frequency during the calibration test by the amplitude of the heave velocity. The unit of the radiation damping coefficient is Newton-seconds per meter, i.e., kilograms per second; the unit of the seawater density is kilograms per cubic meter, the unit of the gravitational acceleration is meters per second squared, and the unit of the buoy's waterline area is square meters. The unit of their product is kilograms per second squared; the unit of the ratio obtained by dividing the radiation damping coefficient by their product is seconds. Therefore, the dimension of the radiation wave transmission coefficient is seconds. The physical meaning of the radiation wave transfer coefficient is that it is the time conversion coefficient required to convert the buoy heave velocity signal (unit: meters per second) into the radiation transmission wavefront elevation signal (unit: meters). After performing point-by-point multiplication on the buoy heave velocity time series using the radiation wave transfer coefficient, the resulting radiation transmission reference signal, in meters, has the same unit as the total transmission wavefront elevation time series on the back wave side. The two can be subtracted point-by-point in step S2 to obtain the scattering transmission time series. Multiplying the translated buoy heave velocity time series point-by-point by the radiation wave transfer coefficient completes the point-by-point scaling, yielding the radiation transmission reference signal.
[0026] The scattering and transmission time series is the residual sequence obtained by subtracting the total transmitted wavefront elevation time series on the back side from the radiation transmission reference signal point by point. Its physical meaning is the time series of the remaining scattering and transmission wavefront elevation components on the back side after removing the radiation component. The discrete implementation of the moving average root mean square (RMS) operation is as follows: The integer number of sampling points obtained by dividing the current incident wave main period by the sampling interval of 0.01 seconds is used as the window length N. At the current time t, the first N sampling points are taken to form the calculation window. The arithmetic mean of the squared values of each point within the window is calculated, and the positive square root is taken as the current RMS output value. This operation is repeated point by point along the time axis, performing the operation on both the scattering and transmission time series and the radiation transmission reference signal, to obtain the scattering and transmission amplitudes respectively. The current incident wave main period is determined in real time by taking the reciprocal of the frequency corresponding to the peak power spectrum after performing a fast Fourier transform on the wavefront elevation time series on the front side. The window length N is updated once for each incident wave main period. The preset minimum value is 0.001 meters. This value is lower than the minimum effective disturbance amount of the actual measurable wavefront elevation. When the scattering and transmission amplitude is lower than 0.001 meters, it is determined that the current working condition is near still water. The component ratio of the previous effective moment is used to replace the current division result. The previous effective moment is defined as the most recent calculation node where the scattering and transmission amplitude is not lower than 0.001 meters. This substitution logic prevents the component ratio from becoming infinitely large due to the divisor approaching zero under near still water conditions.
[0027] In one specific embodiment, step S3, using the relationship between the component ratio and a preset threshold as the switching criterion, includes: Under the calibration state of fixing the float and setting the damping coefficient of the power output device to zero, continuously collect the component ratios within no less than 20 incident wave cycles, and take the arithmetic mean of the component ratio sequence to obtain the reference mean of the component ratios. Multiply the baseline mean of the component ratios by the preset margin coefficient to obtain the preset threshold. The preset margin coefficient is used as a switching margin parameter to make the preset threshold higher than the baseline mean of the component ratios. The preset margin coefficient is a fixed value within the range of greater than 1 and not greater than 1.5. During normal operation of the device, the component ratio at the end of each incident wave main cycle is used as a node. The component ratio at the current moment is compared with a preset threshold. When the component ratio is not lower than the preset threshold, a radiation dominance indicator is output. When the component ratio is lower than the preset threshold, a scattering dominance indicator is output, and the switching criterion comparison result is obtained.
[0028] Specifically, the calculation of the component ratio benchmark mean is performed under calibration conditions with the buoy fixed and the damping coefficient of the power output device set to zero. This condition is consistent with the calibration conditions of the scattering phase calibration value. The purpose is to eliminate the active suppression of the radiation component by the power output device damping, so that the component ratio reflects the natural proportional relationship between the two types of transmitted components under uncontrolled intervention conditions. The basis for continuously collecting no less than 20 incident wave cycles is that in actual sea conditions, the main period of the incident wave exhibits wave-by-wave fluctuations. The mean of 20 cycles can effectively cover the natural variation range of the wave period, making the benchmark mean statistically representative. If the number of collection cycles is too small, the benchmark mean will be significantly affected by a single wave disturbance, causing the preset threshold to deviate from the actual natural ratio. The arithmetic mean is calculated by summing the component ratios output at the end of each main period of the incident wave during the calibration period and dividing by the total number of collection cycles to obtain the benchmark mean of the component ratio. The preset margin coefficient is used as a switching margin parameter to make the preset threshold higher than the average value of the component ratio benchmark. The preset margin coefficient is set within the range of greater than 1 and not greater than 1.5, with a recommended value of 1.2. The basis for setting the value to greater than 1 is that during normal operation, the intervention of the power output device damping suppresses the radiation component, and the component ratio is naturally lower than the calibration state. A margin needs to be left to avoid falsely triggering the radiation dominance indicator in the early stage of normal operation. The basis for setting the value to not greater than 1.5 is that an excessive margin will cause the radiation dominance mode response to be sluggish, so that even after the radiation transmission component increases significantly, the damping increase action cannot be triggered in time. The preset margin coefficient is set within the range of greater than 1 and not greater than 1.5, with a recommended value of 1.2. The reason for setting a value greater than 1 is that during normal operation, the intervention of the power output device damping suppresses the radiation component, and the component ratio is naturally lower than the calibration state. A margin is required to avoid falsely triggering the radiation dominance indicator in the early stage of normal operation. The reason for setting a value not greater than 1.5 is that an excessive margin will cause the radiation dominance mode response to be sluggish, so that even after the radiation transmission component increases significantly, the damping increase action cannot be triggered in time.
[0029] The radiation dominance flag and the scattering dominance flag are single-bit state variables stored in the wave suppression control state register. The radiation dominance flag is represented by a logic value of 1, and the scattering dominance flag is represented by a logic value of 0. They are mutually exclusive, and only one of them is stored in the register at any given time. The triggering mechanism, with the end time of each incident wave main cycle as the node, is implemented as follows: the current incident wave main cycle is determined by taking the reciprocal of the frequency corresponding to the peak power spectrum after performing a fast Fourier transform on the wavefront elevation time series as described in step S2; the end time of each incident wave main cycle is identified by the zero-crossing detection continuously performed on the wavefront elevation time series. The zero-crossing detection takes the moment when the wavefront elevation crosses from a positive value to a negative value as the end time of the current incident wave main cycle, and triggers a comparison operation between the component ratio and a preset threshold at the end time. The execution logic of the comparison operation is as follows: read the latest component ratio output from step S2 in the current wave suppression control status register, subtract it from the preset threshold, write the radiation dominance flag into the wave suppression control status register when the difference is not less than zero, write the scattering dominance flag into the wave suppression control status register when the difference is less than zero, and write the write operation overwrites the status value of the previous moment to obtain the switching criterion comparison result.
[0030] In one specific embodiment, step S3, when the component ratio is not lower than a preset threshold, increases the damping coefficient of the power output device using the positive difference between the component ratio and the preset threshold as an incremental reference, includes: When the switching criterion comparison result is a radiation-dominant indicator, the component ratio is subtracted from the preset threshold to obtain the radiation excess difference. Multiply the excess radiation difference by the radiation gain coefficient, and then multiply by the reference damping coefficient to obtain the damping increment, where the radiation gain coefficient is a preset positive real number and the reference damping coefficient is the factory calibration value of the device. The damping increment is added to the current damping coefficient of the power output device to obtain the updated damping coefficient in the radiation-dominant mode. The updated damping coefficient is then written into the damping control register of the power output device as a wave-damping control command.
[0031] Specifically, the radiation excess difference is the dimensionless excess of the component ratio exceeding a preset threshold. Its physical meaning is the degree of excess of the current radiated transmission component relative to the scattered transmission component. The larger this value, the higher the proportion of the radiated transmission component in the total transmitted wave on the back wave side, requiring a larger damping increment to suppress the float's heave amplitude and thus reduce the radiated transmission component. The recommended value for the radiation gain coefficient is 0.3. This value is based on the following: when the radiation excess difference is 1, the damping increment is 0.3 times the reference damping coefficient, corresponding to a decrease in the float's heave velocity amplitude approximately equal to the ratio of the total damping increase to the original total damping. Numerical simulations have verified that this increase can bring the radiated transmission amplitude back to within the preset threshold level within one to three incident wave main cycles, without causing excessive suppression of the float's motion leading to a compensatory increase in the scattered transmission component exceeding the scattering compensation threshold. The radiation gain coefficient should not exceed 0.5. Exceeding this value will result in an excessively large damping increment, causing a sharp drop in the float's heave amplitude, and the rate of increase in the scattered transmission component exceeding the deflector deployment response speed, triggering control oscillations. The reference damping coefficient is the optimal linear damping value determined by forced oscillation test under design conditions when the device leaves the factory. This value maximizes the energy harvesting efficiency of the device under design conditions. It serves as the proportional reference for calculating the damping increment, ensuring that the increment amplitude maintains a proportional relationship with the inherent damping characteristics of the device.
[0032] The calculation sequence for the damping increment is as follows: First, the radiation excess difference is multiplied by the radiation gain coefficient to obtain a dimensionless scaling factor. Then, this scaling factor is multiplied by the reference damping coefficient to obtain the damping increment with dimensions of the damping coefficient, measured in Newton-seconds per meter. The initial value of the current power output device's damping coefficient is taken as the reference damping coefficient when the device starts up. Subsequently, at the end of each incident wave main cycle, the updated damping coefficient stored in the power output device's damping control register at the previous moment is read as the current value for calculation. The damping increment is added to the current power output device's damping coefficient to obtain the updated damping coefficient in the radiation-dominant mode. This addition operation is performed immediately within the same calculation node after the switching criterion comparison result is written to the wave suppression control status register. The updated damping coefficient is written to the power output device's damping control register in digital form via the control bus. After the power output device reads this register value, it takes effect at the start of the next incident wave main cycle to ensure that the damping coefficient switching is aligned with the wave cycle boundary and to avoid introducing phase disturbances by changing the damping in the middle of the wave cycle. The damping control command is the updated damping coefficient digital value written into the damping control register of the power output device. Its physical form is the opening command of the hydraulic control valve of the power output device or the set value of the excitation current of the linear generator. The specific form is determined by the type of power output device.
[0033] In one specific embodiment, step S3, when the component ratio is lower than a preset threshold, reduces the damping coefficient of the power output device based on the positive difference between the preset threshold and the component ratio, thereby obtaining a wave suppression control command, including: When the switching criterion comparison result is the scattering-dominant identifier, the difference between the preset threshold and the component ratio is used to obtain the scattering excess difference value; Multiply the scattering excess difference by the scattering gain coefficient, and then multiply it by the reference damping coefficient to obtain the damping reduction. The scattering gain coefficient is a preset positive real number and is independent of the radiation gain coefficient. The damping reduction is subtracted from the current power output device damping coefficient to obtain the updated damping coefficient in the scattering-dominant mode. When the updated damping coefficient in the scattering-dominant mode is lower than the product of the reference damping coefficient and the preset lower limit coefficient, the updated damping coefficient in the scattering-dominant mode is replaced by the product of the reference damping coefficient and the preset lower limit coefficient. The replaced value is written into the power output device damping control register as a wave suppression control command.
[0034] Specifically, the scattering excess difference is a dimensionless excess exceeding the current component ratio by a preset threshold. Physically, it represents the degree of excess of the current scattered transmission component relative to the radiated transmission component. A larger value indicates a higher proportion of the scattered transmission component in the total transmitted wave on the back wave side, requiring a larger damping reduction to amplify the float's heave amplitude and enhance the dynamic interception effect of the asymmetric baffle on the incident wave, thus reducing the scattered transmission component due to the obstruction of the incident wave's diffraction path. The recommended value for the scattering gain coefficient is 0.25, set independently of the radiation gain coefficient of 0.3. This independent setting is based on the fact that the response rate of the scattered transmission component to the reduction in damping is physically asymmetric to the response rate of the radiated transmission component to the increase in damping. The reduction of the scattered transmission component depends on the dynamic shielding effect of the incident wave after the float's heave amplitude increases. This process requires approximately three to five main incident wave cycles to fully establish, and the response rate is slower than the direct suppression of the radiated transmission component. Therefore, the scattering gain coefficient is set slightly smaller than the radiation gain coefficient to avoid excessive damping reduction leading to a sudden increase in the float's heave amplitude and causing structural safety issues. The scattering gain coefficient should not exceed 0.4. If it exceeds this value, the damping reduction in a single cycle will exceed 40% of the reference damping coefficient, the buoy heave amplitude will increase too quickly, the bottom water depth of the asymmetric baffle will exceed the design range, and the geometric asymmetry of the triangular baffle will be destroyed.
[0035] The preset lower limit coefficient is a proportional parameter to ensure that the damping coefficient of the power output device is not lower than the lower limit of safe operation. The recommended value is 0.3, which means that the minimum allowable value of the damping coefficient of the power output device shall not be lower than 30% of the reference damping coefficient. The basis for choosing a value of 0.3 is that when the damping coefficient of the power output device is lower than 30% of the reference damping coefficient, the constraint force of the power output device on the heave motion of the float is insufficient to maintain the stable operation of the float within the design draft range. The float is at risk of the asymmetric baffle being completely submerged due to excessive heave amplitude. After being submerged, the baffle loses its hydrodynamic effect, and the scattered transmission components increase instead of decreasing, causing the control strategy to fail. The lower limit protection logic is executed as follows: after completing the subtraction operation between the damping reduction and the current power output device damping coefficient, the calculation result is immediately compared with the product of the reference damping coefficient and the preset lower limit coefficient. If the subtraction result is lower than the product value, the product value is used to replace the subtraction result. The replaced value is written into the power output device damping control register as the updated damping coefficient in the scattering-dominant mode. The writing timing is the same as in the radiation-dominant mode, which is executed at the end of the current incident wave main cycle. The updated value takes effect at the beginning of the next incident wave main cycle to ensure that the damping coefficient change is strictly aligned with the wave cycle boundary.
[0036] In one specific embodiment, step S4 includes: The moment when the wave suppression control command is first switched to the radiation-dominant flag is recorded as the reference time, and the arithmetic mean of the scattered transmission amplitude within the main period of the incident wave before the reference time is used as the scattered transmission amplitude at the reference time. Taking the end of each incident wave main cycle as a node, the difference between the current scattering and transmission amplitude and the scattering and transmission amplitude at the reference time is used to obtain the scattering and transmission increment. The scattering and transmission increment is compared with the scattering compensation threshold to obtain the guide plate triggering criterion. When the deflector trigger criterion indicates that the scattering transmission increment exceeds the scattering compensation threshold, the back wave side deflector is deployed from the storage position to the preset working angle to obtain the deflector deployment state. With the deflector deployed, the scattering and transmission amplitude is continuously monitored at the end of each incident wave main cycle. When the scattering and transmission amplitude drops to within the sum of the scattering and transmission amplitude at the reference time and the preset fallback margin, the deflector on the back wave side is retracted to the storage position, thus obtaining the deflector retracted state. The preset fallback margin is the product of the scattering and transmission amplitude at the reference time and the preset fallback ratio coefficient. The preset fallback ratio coefficient is a fixed value within the range of greater than 0 and not greater than 0.05.
[0037] Specifically, the reference time is defined as the moment when the value stored in the wave suppression control status register first switches from the scattering-dominant flag to the radiation-dominant flag. This moment strictly corresponds to the writing time of the switching criterion comparison result, with the end of the incident wave main cycle as the trigger node. The calculation window for the scattering and transmission amplitude at the reference time is the complete incident wave main cycle preceding the reference time. The number of sampling points in this cycle is taken as the integer value obtained by dividing the current incident wave main cycle length by the sampling interval of 0.01 seconds. The scattering and transmission amplitude of each sampling point within the window is accumulated and then divided by the total number of sampling points to obtain the scattering and transmission amplitude at the reference time. The scattering compensation threshold is calculated as follows: the root mean square value of the wave surface elevation time series on the wave-facing side within the main period of the incident wave one complete incident wave before the reference time is taken, and multiplied by the preset compensation ratio coefficient of 0.05 to obtain the scattering compensation threshold. The preset compensation ratio coefficient of 0.05 is based on the fact that this value corresponds to an increase in the scattering transmission amplitude exceeding 5% of the root mean square value of the current incident effective wave height. This ratio reflects the minimum detectable change in the scattering transmission component that shows a statistically significant increase. Increments below this ratio are within the normal random fluctuation range of waves and should not trigger the deflector deployment action to avoid frequent mechanical wear.
[0038] The preset working angle is determined by a water tank calibration test before device installation: Under the calibration state with the float fixed and only the back-wave side scattered transmitted wave propagation retained, at least three wave pressure sensors with uniform vertical spacing are arranged on the back-wave side. The phase difference of the scattered transmitted wave received by each sensor is used to infer the main propagation elevation angle of the scattered transmitted wave in the vertical section. The preset working angle is the complementary angle of this elevation angle, so that the normal direction of the guide plate after deployment is orthogonal to the propagation direction of the scattered transmitted wave. This orthogonal relationship causes the scattered transmitted wave to strike the surface of the guide plate at the maximum incident angle and reflect back to the front-wave side according to the solid boundary reflection law, achieving the maximum refraction and deflection effect. The action of the guide plate deploying from the storage position to the preset working angle is driven by a spring reset mechanism. The deployment time is controlled within one-tenth of the current incident wave main period. This time constraint ensures that the deployment action is completed within a single wave cycle and does not cross the boundary of adjacent wave cycles, thus introducing phase interference. The preset fallback margin is the product of the scattered transmission amplitude at the reference time and the preset fallback ratio coefficient. The preset fallback ratio coefficient is a fixed value within the range of greater than 0 and not greater than 0.05. The preset fallback ratio coefficient is used as the lower limit margin parameter for judging the scattered transmission amplitude when the back wave side guide vane switches from the deployed state to the retracted state. The recommended value of the preset fallback ratio coefficient is 0.02, that is, the preset fallback margin is 0.02 times the scattered transmission amplitude at the reference time. The basis for this value is that the scattered transmission amplitude fluctuates periodically during the fallback process after the guide vane is deployed. A margin of 0.02 times can accommodate amplitude jitter within the normal fluctuation range and prevent the scattered transmission amplitude from repeatedly crossing the compensation threshold boundary, causing the guide vane to frequently deploy and retract. The basis for the preset fallback ratio coefficient not exceeding 0.05 is that the upper limit of this ratio does not exceed the ratio coefficient corresponding to the scattering compensation threshold, so as not to affect the trigger sensitivity of the guide vane deployment.
[0039] The deflector retracted state and the deflector deployed state are mutually exclusive. At any given time, the deflector is only in one of the states. After the retraction action is executed, the deflector status bit in the wave suppression control status register is reset to the storage bit flag, waiting for the next deflector trigger criterion over-threshold event.
[0040] Figure 2 This is a schematic diagram illustrating the dynamic response of scattering and transmission amplitude and the triggering control process of the back-wave side guide vane in an embodiment of this application. Figure 2As shown, the horizontal axis represents time, the vertical axis represents the scattered transmission amplitude, the solid line is the time-domain variation curve of the scattered transmission amplitude, the horizontal dashed line is the scattered transmission amplitude at the reference time, the horizontal dotted line is the scattering compensation threshold, and the horizontal dotted line is the upper limit of the retraction criterion. After the wave-damping control command is first switched to the radiation-dominant indicator, the damping coefficient of the power output device increases, causing the float heave amplitude to decrease. The scattered transmission amplitude increases compensatorily and exceeds the scattering compensation threshold, triggering the guide vane trigger criterion. The guide vane on the back wave side unfolds from the storage position to the preset working angle. The area filled with diagonal lines represents the stage of increased scattering compensation, and the area filled with gray represents the stage of guide vane unfolding and wave-damping. Subsequently, the scattered transmission amplitude continues to fall. When it falls to within the sum of the scattered transmission amplitude at the reference time and the preset fall margin, the guide vane retracts to the storage position, and the scattered transmission amplitude stabilizes near the scattered transmission amplitude at the reference time.
[0041] The broadband wave suppression method of the asymmetric baffle float wave energy device in the embodiments of this application has been described above. The broadband wave suppression system of the asymmetric baffle float wave energy device in the embodiments of this application is described below. One embodiment of the broadband wave suppression system of the asymmetric baffle float wave energy device in the embodiments of this application includes: The acquisition module is used to acquire the wave surface elevation on the wave-facing side of the asymmetric baffle float, the total transmitted wave surface elevation on the wave-receiving side, and the float heave velocity to obtain the raw transmitted wave data set. The calculation module is used to obtain a radiation transmission reference signal by shifting the float heave velocity in the original transmitted wave data group by a quarter-cycle phase, to obtain a scattering transmission amplitude by subtracting the total transmitted wave surface elevation on the back wave side from the radiation transmission reference signal and then performing a sliding root mean square operation on the radiation transmission reference signal, to obtain a radiation transmission amplitude, and to obtain a component ratio by dividing the radiation transmission amplitude by the scattering transmission amplitude. The output module is used to use the relationship between the component ratio and a preset threshold as a switching criterion. When the component ratio is not lower than the preset threshold, the damping coefficient of the power output device is increased based on the positive difference between the component ratio and the preset threshold as an incremental reference. When the component ratio is lower than the preset threshold, the damping coefficient of the power output device is decreased based on the positive difference between the preset threshold and the component ratio as a decremental reference, thereby obtaining a wave suppression control command. The retraction module is used to determine the scattering and transmission amplitude relative to the reference time. When the increment exceeds the scattering compensation threshold, the back wave side guide plate is deployed to a preset working angle, and the guide plate is retracted after the scattering and transmission amplitude returns to the compensation range.
[0042] This invention also provides a broadband wave-damping device for an asymmetric baffle float wave energy system, which can be a server. This broadband wave-damping device includes a processor, memory, display screen, input device, network interface, and database connected via a system bus. The processor, designed as a computer, provides computing and control capabilities. The memory of this broadband wave-damping device includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The database of this broadband wave-damping device stores the data corresponding to this embodiment. The network interface of this broadband wave-damping device is used for communication with external terminals via a network connection. When the computer program is executed by the processor, it implements the above-described method.
[0043] The present invention also provides a computer-readable storage medium, which can be a non-volatile computer-readable storage medium or a volatile computer-readable storage medium, wherein the computer-readable storage medium stores instructions that, when the instructions are executed on a computer, cause the computer to perform the steps of the broadband wave attenuation method of the asymmetric baffle float wave energy device.
[0044] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the systems and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.
[0045] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause an asymmetric baffle float wave energy device broadband wave attenuation device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.
[0046] The above embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims
1. A broadband wave attenuation method for an asymmetric baffle float wave energy device, characterized in that, The method includes: Step S1: Collect the wave surface elevation on the wave-facing side of the asymmetric baffle float, the total transmitted wave surface elevation on the wave-receiving side, and the float heave velocity to obtain the raw transmitted wave data set. Step S2: The buoy heave velocity in the original transmitted wave data set is phase-shifted by a quarter period to obtain the radiation transmission reference signal. The total transmitted wave surface elevation on the back wave side is subtracted from the radiation transmission reference signal and then the scattering transmission amplitude is obtained by moving root mean square operation. The radiation transmission amplitude is obtained by performing moving root mean square operation on the radiation transmission reference signal. The radiation transmission amplitude is divided by the scattering transmission amplitude to obtain the component ratio. Step S3: Using the relationship between the component ratio and the preset threshold as the switching criterion, when the component ratio is not lower than the preset threshold, the damping coefficient of the power output device is increased based on the positive difference between the component ratio and the preset threshold as the increment; when the component ratio is lower than the preset threshold, the damping coefficient of the power output device is decreased based on the positive difference between the preset threshold and the component ratio as the decrement, thereby obtaining a wave suppression control command. Step S4: Using the increment of the scattered transmission amplitude relative to the reference time as a criterion, when the increment exceeds the scattering compensation threshold, the back wave side guide plate is deployed to a preset working angle, and the guide plate is retracted after the scattered transmission amplitude returns to the compensation range.
2. The broadband wave attenuation method for the asymmetric baffle float wave energy device according to claim 1, characterized in that, Step S1 includes: A wave height sensor is arranged on the wave-facing side and the wave-back side of the asymmetric baffle float. The output signals of the wave height sensor on the wave-facing side and the wave height sensor on the wave-back side are synchronously collected based on a preset sampling frequency to obtain the wave surface elevation time series on the wave-facing side and the total transmitted wave surface elevation time series on the wave-back side. The heave acceleration signal output by the vertical accelerometer at the center of mass of the float is integrated with the initial value of zero heave velocity under the static water level equilibrium state, and a first-order numerical integration is performed to obtain the time series of heave velocity of the float. With the calibrated state of fixed float and power output device damping coefficient set to zero, the cross-correlation operation is performed on the wave front elevation time series and the total transmitted wave front elevation time series on the back side to obtain the fixed phase difference between the two, which is recorded as the scattering phase calibration value. The time series of wavefront elevation, the time series of total transmitted wavefront elevation, the time series of float heave velocity, and the scattering phase calibration value are time-stamped and merged to obtain the original transmitted wave data set.
3. The broadband wave attenuation method for the asymmetric baffle float wave energy device according to claim 1, characterized in that, Step S2 includes: Based on the time series of the buoy heave velocity in the original transmitted wave data set, the number of sampling points corresponding to one-quarter of the current incident wave main period is used as the translation step size to perform a forward time-series translation on the buoy heave velocity time series. Then, the translated sequence is scaled point by point using the radiation wave transfer coefficient to obtain the radiation transmission reference signal. The radiation wave transfer coefficient is calculated by dividing the radiation damping coefficient measured under calibration conditions by the product of seawater density, gravitational acceleration, and buoy waterline area. The total transmitted wave surface elevation time series on the back wave side in the original transmitted wave data set is subtracted point by point from the radiation transmission reference signal to obtain the scattering transmission time series. The scattering transmission time series is then subjected to a sliding root mean square operation with the number of sampling points corresponding to the current main period of the incident wave as the window length to obtain the scattering transmission amplitude. Using the number of sampling points corresponding to the current main period of the incident wave as the window length, a sliding root mean square operation is performed on the radiation transmission reference signal to obtain the radiation transmission amplitude. Divide the radiation transmission amplitude by the scattering transmission amplitude. When the scattering transmission amplitude is lower than a preset minimum value, replace the current calculation result with the quotient value of the previous effective time to obtain the component ratio.
4. The broadband wave attenuation method for the asymmetric baffle float wave energy device according to claim 3, characterized in that, Step S3, using the relationship between the component ratio and a preset threshold as the switching criterion, includes: Under the calibration state of fixing the float and setting the damping coefficient of the power output device to zero, the component ratios are continuously collected for no less than 20 incident wave cycles, and the arithmetic mean of the component ratio sequence is taken to obtain the component ratio benchmark mean. The preset threshold is obtained by multiplying the mean of the component ratio benchmark by a preset margin coefficient. The preset margin coefficient is used as a switching margin parameter to make the preset threshold higher than the mean of the component ratio benchmark. The preset margin coefficient is a fixed value in the range of greater than 1 and not greater than 1.
5. During normal operation of the device, the component ratio at the end of each incident wave main cycle is used as a node. The component ratio at the current time is compared with the preset threshold. When the component ratio is not lower than the preset threshold, a radiation dominance indicator is output. When the component ratio is lower than the preset threshold, a scattering dominance indicator is output, and a switching criterion comparison result is obtained.
5. The broadband wave attenuation method for the asymmetric baffle float wave energy device according to claim 4, characterized in that, Step S3, when the component ratio is not lower than the preset threshold, involves increasing the power output device damping coefficient using the positive difference between the component ratio and the preset threshold as an incremental benchmark, including: When the switching criterion comparison result is a radiation-dominant indicator, the difference between the component ratio and the preset threshold is used to obtain the radiation excess difference. Multiply the radiation excess difference by the radiation gain coefficient, and then multiply by the reference damping coefficient to obtain the damping increment, wherein the radiation gain coefficient is a preset positive real number, and the reference damping coefficient is the factory calibration value of the device. The damping increment is added to the current power output device damping coefficient to obtain the updated damping coefficient in the radiation-dominant mode. The updated damping coefficient is then written into the power output device damping control register as a wave-damping control command.
6. The broadband wave attenuation method for the asymmetric baffle float wave energy device according to claim 5, characterized in that, In step S3, when the component ratio is lower than the preset threshold, the damping coefficient of the power output device is reduced based on the positive difference between the preset threshold and the component ratio to obtain a wave suppression control command, including: When the switching criterion comparison result is a scattering dominance indicator, the difference between the preset threshold and the component ratio is used to obtain the scattering excess difference value; The scattering excess difference is multiplied by the scattering gain coefficient and then multiplied by the reference damping coefficient to obtain the damping reduction, wherein the scattering gain coefficient is a preset positive real number and is independent of the radiation gain coefficient. The damping reduction is subtracted from the current power output device damping coefficient to obtain the updated damping coefficient in the scattering-dominant mode. When the updated damping coefficient in the scattering-dominant mode is lower than the product of the reference damping coefficient and the preset lower limit coefficient, the updated damping coefficient in the scattering-dominant mode is replaced by the product of the reference damping coefficient and the preset lower limit coefficient. The replaced value is written into the power output device damping control register as the wave suppression control command.
7. The broadband wave attenuation method for the asymmetric baffle float wave energy device according to claim 1, characterized in that, Step S4 includes: The moment when the wave suppression control command first switches to the radiation dominance flag is recorded as the reference time, and the arithmetic mean of the scattered transmission amplitude within the main period of the incident wave preceding the reference time is used as the scattered transmission amplitude at the reference time. Taking the end of each incident wave main cycle as a node, the difference between the current scattering and transmission amplitude and the scattering and transmission amplitude at the reference time is used to obtain the scattering and transmission increment. The scattering and transmission increment is compared with the scattering compensation threshold to obtain the guide plate triggering criterion. When the guide plate trigger criterion indicates that the scattering transmission increment exceeds the scattering compensation threshold, the back wave side guide plate is deployed from the storage position to the preset working angle to obtain the guide plate deployment state; With the deflector deployed, the scattering and transmission amplitude is continuously monitored at the end of each incident wave main cycle. When the scattering and transmission amplitude drops to within the sum of the scattering and transmission amplitude at the reference time and the preset fallback margin, the deflector on the back wave side is retracted to the storage position, thus obtaining the deflector retracted state. The preset fallback margin is the product of the scattering and transmission amplitude at the reference time and the preset fallback ratio coefficient, and the preset fallback ratio coefficient is a fixed value within the range of greater than 0 and not greater than 0.
05.
8. A broadband wave suppression system for an asymmetric baffle float wave energy device, characterized in that, For implementing the broadband wave suppression method of the asymmetric baffle float wave energy device as described in any one of claims 1-7, the broadband wave suppression system of the asymmetric baffle float wave energy device comprises: The acquisition module is used to acquire the wave surface elevation on the wave-facing side of the asymmetric baffle float, the total transmitted wave surface elevation on the wave-receiving side, and the float heave velocity to obtain the raw transmitted wave data set. The calculation module is used to obtain a radiation transmission reference signal by shifting the float heave velocity in the original transmitted wave data group by a quarter-cycle phase, to obtain a scattering transmission amplitude by subtracting the total transmitted wave surface elevation on the back wave side from the radiation transmission reference signal and then performing a sliding root mean square operation on the radiation transmission reference signal, to obtain a radiation transmission amplitude, and to obtain a component ratio by dividing the radiation transmission amplitude by the scattering transmission amplitude. The output module is used to use the relationship between the component ratio and a preset threshold as a switching criterion. When the component ratio is not lower than the preset threshold, the damping coefficient of the power output device is increased based on the positive difference between the component ratio and the preset threshold as an incremental reference. When the component ratio is lower than the preset threshold, the damping coefficient of the power output device is decreased based on the positive difference between the preset threshold and the component ratio as a decremental reference, thereby obtaining a wave suppression control command. The retraction module is used to determine the scattering and transmission amplitude relative to the reference time. When the increment exceeds the scattering compensation threshold, the back wave side guide plate is deployed to a preset working angle, and the guide plate is retracted after the scattering and transmission amplitude returns to the compensation range.
9. A broadband wave suppression device for an asymmetric baffle float wave energy system, characterized in that, It includes a memory and a processor, the memory storing a computer program that can run on the processor, and the processor executing the computer program to implement the broadband wave attenuation method of the asymmetric baffle float wave energy device according to any one of claims 1 to 7.
10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is run by the processor, it causes the processor to execute the broadband wave suppression method of the asymmetric baffle float wave energy device as described in any one of claims 1 to 7.
Citation Information
Patent Citations
Asymmetrical F-type hinged heaving floating type wave suppressor
CN104099895A
Offshore pitching and heaving coupling type wave energy and wind energy integrated power generation device
CN109441703A