CN-122018027-A - Convection triggering gravity wave three-dimensional identification method and device

CN122018027ACN 122018027 ACN122018027 ACN 122018027ACN-122018027-A

Abstract

The invention provides a three-dimensional recognition method and device for gravity waves triggered by convection, the method comprises the steps of generating three-dimensional lattice point data by fusing multi-source weather observation data, recognizing a moving path of a convection system according to the three-dimensional lattice point data, sampling vertical profiles of the three-dimensional lattice point data along the moving path within a preset statistical radius range to obtain a time sequence vertical profile set containing power and micro-physical parameters, carrying out spectrum filtering and wave number spectrum analysis of physical constraints on middle beta-scale vertical wind shear time sequence vertical profiles in the time sequence vertical profile set, extracting gravity wave signals, calculating spectrum slope matching degree of the gravity wave signals, gravity wave vertical coherence between adjacent height layers and frequency range matching degree of the gravity wave signals, and obtaining comprehensive confidence degree of the gravity wave signals according to the spectrum slope matching degree, the vertical consistency and the frequency range matching degree. The invention realizes objective, quantitative and three-dimensional accurate recognition of gravitational waves.

Inventors

HUANG XUANXUAN
HUANG JUAN
JIA YAN
JIANG SHUJIE
LI RUNQIU
LI WENJUAN
ZHOU LINGLI
ZHAO LU
CHEN LIE
WANG LIYING
SU GUIYANG
ZHOU KAI
LUO RAN

Assignees

浙江省气象台

Dates

Publication Date: 20260512
Application Date: 20260410

Claims (10)

1. A method for three-dimensional recognition of gravity waves triggered by convection, comprising: Generating three-dimensional grid point data by fusing multi-source meteorological observation data, and identifying a moving path of a streaming system according to the three-dimensional grid point data; performing vertical profile sampling on the three-dimensional lattice point data in a preset statistical radius range along the moving path to obtain a time sequence vertical profile set containing power and micro-physical parameters; Performing spectrum filtering and wave number spectrum analysis of physical constraint on the middle beta-scale vertical wind-cut time sequence vertical profile in the time sequence vertical profile set, and extracting gravitational wave signals; And calculating the spectrum slope matching degree of the gravity wave signal, the gravity wave vertical coherence between adjacent height layers and the frequency range matching degree of the gravity wave signal, and obtaining the comprehensive confidence degree of the gravity wave signal according to the spectrum slope matching degree, the vertical coherence and the frequency range matching degree.
2. The method of claim 1, wherein the step of sampling the vertical profile of the three-dimensional lattice point data along the movement path within a predetermined statistical radius to obtain a time-series vertical profile set including dynamic and micro-physical parameters comprises: For differential propagation phase shift rate in the three-dimensional lattice data A field searching for vertical accumulation within a preset statistical radius range centered on the position at each moment in the movement path The largest horizontal position is used as the convection microphysics centroid; for a vertical wind shear VS field in the three-dimensional grid point data, at each moment in the movement path Searching a horizontal position with the maximum vertical cumulative absolute VS as the optimal point of power disturbance in a preset statistical radius range with the position as the center; extracting a polarization amount vertical profile at the convection microphysical centroid at each moment, and extracting a dynamic parameter vertical profile at the optimal point of dynamic disturbance at each moment; and generating a time sequence vertical profile set containing dynamic and micro physical parameters according to the polarization amount vertical profile and the dynamic parameter vertical profile.
3. The method of claim 1, wherein the performing physically constrained spectral filtering and wavenumber spectral analysis on the mid- β scale vertical wind shear temporal vertical profile in the temporal vertical profile set extracts a gravitational wave signal, comprising: Dividing the middle beta-scale vertical wind-cutting time sequence vertical profile into a plurality of mutually overlapped sequences by adopting a sliding window method, and multiplying the sequences on a specific height layer in each sliding window point by a window function to obtain a windowed sequence; performing discrete Fourier transform on the windowed sequence to obtain a complex frequency spectrum and a power spectrum density of the gravitational wave signal; Determining physical constraints to be met by the frequency of the gravitational wave signal according to the key atmospheric parameters; Determining a band-pass filter according to the physical constraint, and extracting a filtered gravitational wave spectrum from the complex spectrum by adopting the band-pass filter; performing inverse Fourier transform on the filtered gravitational wave spectrum to reconstruct a time domain gravitational wave disturbance signal; and calculating the filtered power spectrum density according to the filtered gravitational wave spectrum.
4. A method of three-dimensional recognition of a convection-triggered gravitational wave as claimed in claim 3, wherein the physical constraints to be met by the frequency of the gravitational wave signal are: ; ; Wherein, the For the frequency of the gravitational wave signal, And The minimum and maximum frequency thresholds of the gravitational wave signal respectively, Is the coriolis frequency of the signal and, Is the buoyancy frequency near the top of the troposphere.
5. The method for three-dimensional recognition of gravity waves triggered by convection according to claim 4, wherein the boundaries of the band-pass filter adopt gaussian transitions to smooth the truncation effect, and the formula is: ; Wherein, the As a response function of the band-pass filter, Is the nyquist frequency.
6. The method of claim 4, wherein calculating a spectral slope match of the gravitational wave signal comprises: In the frequency range of the gravitational wave signal Performing linear least square fitting on the filtered power spectrum density in a double-logarithmic coordinate system to obtain a fitted spectrum slope; and calculating the matching degree of the fitting spectrum slope and the gravity wave linear theory classical value.
7. A method of three-dimensional recognition of gravity waves triggered by convection according to claim 3, characterized in that calculating the vertical coherence of gravity waves between adjacent height layers comprises: Carrying out standardization processing on the time sequences of adjacent vertical layers in the time domain gravity wave disturbance signals to obtain standardized time sequences; Calculating zero hysteresis correlation coefficient, symbol consistency rate and waveform similarity for the standardized time sequences of adjacent vertical layers; And carrying out weighted fusion on the zero-lag correlation coefficient, the symbol consistency rate and the waveform similarity to obtain the gravity wave vertical coherence.
8. The method of claim 4, wherein calculating a frequency range match for the gravitational wave signal comprises: Calculating the gravitational wave signal in the frequency range The ratio of the energy in the frequency range to the total energy in the full frequency band is used as the matching degree of the frequency range.
9. The method of any one of claims 1-8, further comprising, after calculating the degree of matching of the spectral slope of the gravitational wave signal, the vertical coherence of gravitational waves between adjacent height layers, and the degree of matching of the frequency range of the gravitational wave signal: And fusing the spectrum slope matching degree of the gravitational wave signals, the gravitational wave vertical coherence between adjacent height layers and the frequency range matching degree with a micro physical field and a dynamic field in the time sequence vertical profile set to perform three-dimensional visualization.
10. A convection-triggered gravitational wave three-dimensional identification device, comprising: the identification module is used for generating three-dimensional lattice point data by fusing the multi-source meteorological observation data and identifying a moving path of the flow system according to the three-dimensional lattice point data; The sampling module is used for sampling the vertical profile of the three-dimensional lattice point data in a preset statistical radius range along the moving path to obtain a time sequence vertical profile set containing power and micro-physical parameters; The analysis module is used for carrying out spectrum filtering and wave number spectrum analysis of physical constraint on the middle beta-scale vertical wind-cut time sequence vertical profile in the time sequence vertical profile set and extracting gravitational wave signals; The matching module is used for calculating the spectrum slope matching degree of the gravity wave signals, the gravity wave vertical coherence between adjacent height layers and the frequency range matching degree of the gravity wave signals, and obtaining the comprehensive confidence degree of the gravity wave signals according to the spectrum slope matching degree, the vertical coherence and the frequency range matching degree.

Description

Convection triggering gravity wave three-dimensional identification method and device Technical Field The invention relates to the technical field of meteorological monitoring, in particular to a convection triggering gravitational wave three-dimensional identification method and device. Background The convection triggering gravitational wave is a key dynamic process for connecting a small-scale convection system with a medium-scale circulation, accurately identifies the activity of the gravitational wave, and has important precursor indication significance for the short-time forecasting and early warning of strong convection weather (such as thunderstorm strong wind and short-time strong rainfall). However, existing gravitational wave identification techniques rely heavily on subjective experience or spectral analysis of a single physical quantity, with the following prominent bottlenecks: (1) The identification subjectivity lacks objective quantitative criteria, and the automatic and business stable identification is difficult to realize. (2) The analysis dimension is single, and is mostly limited to a two-dimensional plane or a single vertical section, so that the structure, propagation and evolution of gravitational waves in a three-dimensional space cannot be completely revealed. (3) The reliability of the result is difficult to evaluate, and a quantification system is lacked to distinguish a real gravitational wave signal from other disturbances such as environmental noise, turbulence and the like. Disclosure of Invention Aiming at the problems existing in the prior art, the invention provides a convection triggering gravitational wave three-dimensional identification method and device, which can realize objective, quantitative and three-dimensional gravitational wave accurate identification. The invention provides a convection triggering gravitational wave three-dimensional identification method, which comprises the following steps: Generating three-dimensional grid point data by fusing multi-source meteorological observation data, and identifying a moving path of a streaming system according to the three-dimensional grid point data; performing vertical profile sampling on the three-dimensional lattice point data in a preset statistical radius range along the moving path to obtain a time sequence vertical profile set containing power and micro-physical parameters; Performing spectrum filtering and wave number spectrum analysis of physical constraint on the middle beta-scale vertical wind-cut time sequence vertical profile in the time sequence vertical profile set, and extracting gravitational wave signals; And calculating the spectrum slope matching degree of the gravity wave signal, the gravity wave vertical coherence between adjacent height layers and the frequency range matching degree of the gravity wave signal, and obtaining the comprehensive confidence degree of the gravity wave signal according to the spectrum slope matching degree, the vertical coherence and the frequency range matching degree. According to the method for three-dimensionally recognizing the gravity wave triggered by the convection, provided by the invention, vertical profile sampling is carried out on the three-dimensional lattice point data in a preset statistical radius range along the moving path to obtain a time sequence vertical profile set containing power and micro physical parameters, and the method comprises the following steps: For differential propagation phase shift rate in the three-dimensional lattice data A field searching for vertical accumulation within a preset statistical radius range centered on the position at each moment in the movement pathThe largest horizontal position is used as the convection microphysics centroid; for a vertical wind shear VS field in the three-dimensional grid point data, at each moment in the movement path Searching a horizontal position with the maximum vertical cumulative absolute VS as the optimal point of power disturbance in a preset statistical radius range with the position as the center; extracting a polarization amount vertical profile at the convection microphysical centroid at each moment, and extracting a dynamic parameter vertical profile at the optimal point of dynamic disturbance at each moment; and generating a time sequence vertical profile set containing dynamic and micro physical parameters according to the polarization amount vertical profile and the dynamic parameter vertical profile. According to the method for three-dimensionally identifying the gravity wave triggered by convection, which is provided by the invention, spectrum filtering and wave number spectrum analysis of physical constraint are carried out on the middle beta scale vertical wind shear time sequence vertical profile in the time sequence vertical profile set, and the gravity wave signal is extracted, and the method comprises the following steps: Dividing the middle beta-scale vertical wind-c