Method, device and equipment for determining buoyancy generation item of atmospheric turbulence kinetic energy and storage medium

Abstract

The invention provides a method, device and equipment for determining a buoyancy generation item of atmospheric turbulence kinetic energy and a storage medium, and the method comprises the steps: obtaining the wind speed at a plurality of moments through a wind measurement radar within a specified time length, and obtaining the buoyancy generation item of the atmospheric turbulence kinetic energy according to the wind speed at the plurality of moments; the dissipation rate and pulsation components corresponding to the wind speed in different coordinate axis directions are determined, a turbulence trend item, a shear generation item and a turbulence transmission item are determined according to the pulsation components corresponding to the different coordinate axis directions, and the wind speed is acquired by adopting a wind measurement radar, so that the spatial resolution and the data continuity of data are improved, and the detection range is expanded. And more accurate dissipation rate, turbulence trend term, shear generation term and turbulence transmission term are determined, so that the buoyancy generation term of the atmospheric turbulence kinetic energy can be more accurately deduced according to the dissipation rate, turbulence trend term, shear generation term and turbulence transmission term, and the accuracy of the buoyancy generation term is ensured.

Description

Technical Field

[0001] The present application relates to the field of meteorological observation, and in particular to a method, device, equipment and storage medium for determining the buoyancy generation term of atmospheric turbulent kinetic energy. Background Art

[0002] At present, in meteorological monitoring, the buoyancy generation term of turbulent kinetic energy is a key physical quantity that describes the process of turbulence generation and dissipation. Its research is crucial to revealing the nature and characteristics of turbulence. The monitoring accuracy of the buoyancy generation term of turbulent kinetic energy will affect the accuracy of meteorological monitoring. At present, on the one hand, data is usually collected by observation means such as meteorological towers or drones, and the spatial resolution of the data is low, the data continuity is poor, and the detection range is limited; on the other hand, the buoyancy generation term is calculated directly using collected data such as temperature. Due to the above-mentioned shortcomings in the data, it is easy to affect the accuracy of the buoyancy generation term. Summary of the invention

[0003] The embodiments of the present application provide a method, device, equipment and storage medium for determining the buoyancy generation term of atmospheric turbulent kinetic energy to solve at least one problem existing in the related art. The technical solution is as follows: In a first aspect, an embodiment of the present application provides a method for determining a buoyancy generation term of atmospheric turbulent kinetic energy, comprising: Within a specified time period, the wind speed at several moments is obtained through wind radar; Determine the dissipation rate and the pulsating components of the wind speed corresponding to different coordinate axis directions according to the wind speed at a plurality of moments; According to the pulsating components corresponding to different coordinate axis directions, the turbulence trend term, the shear generation term and the turbulence transport term are determined; The buoyancy generation term of the atmospheric turbulence kinetic energy is determined according to the dissipation rate, the turbulence tendency term, the shear generation term and the turbulence transport term.

[0004] In one embodiment, determining the dissipation rate according to the wind speed at several moments includes: Determining an average wind speed based on the wind speeds at a number of times; Processing the wind speed at a plurality of moments by fast Fourier transform to determine a turbulence spectrum, a first frequency, and a second frequency; Determining a first wave number according to the first frequency and the average wind speed, and determining a second wave number according to the second frequency and the average wind speed; A first spectrum value is determined according to the turbulence spectrum and the first wave number, a second spectrum value is determined according to the turbulence spectrum and the second wave number, and a dissipation rate is determined according to the first spectrum value, the second spectrum value, the first wave number, and the second wave number.

[0005] In one embodiment, determining the pulsating components of the wind speed corresponding to different coordinate axis directions according to the wind speed at a plurality of moments includes: Determine, based on the wind speeds at a plurality of moments, a first component of the average wind speed on the first coordinate axis, a second component of the second coordinate axis, and a third component of the third coordinate axis in the vertical direction; Determine the difference between the first current component of the wind speed on the first coordinate axis at the current moment and the first component, and obtain the first pulsating component corresponding to the wind speed in the direction of the first coordinate axis; determine the difference between the second current component of the wind speed on the second coordinate axis at the current moment and the second component, and obtain the second pulsating component corresponding to the wind speed in the direction of the second coordinate axis; determine the difference between the third current component of the wind speed on the third coordinate axis at the current moment and the third component, and obtain the third pulsating component corresponding to the wind speed in the direction of the third coordinate axis.

[0006] In one embodiment, determining the turbulence trend term according to the pulsation components corresponding to different coordinate axis directions includes: determining a first sum of the square of the first pulsation component, the square of the second pulsation component, and the square of the third pulsation component; The turbulent kinetic energy is determined according to half of the first sum value, and the turbulent kinetic energy is time-derived to obtain a turbulent trend term.

[0007] In one embodiment, determining the shear generation term according to the pulsation components corresponding to different coordinate axis directions includes: Determine a first derivative of the first current component in the direction of the third coordinate axis and a second derivative of the second current component in the direction of the third coordinate axis; determining a first covariance of the first pulsation component and the third pulsation component and a second covariance of the second pulsation component and the third pulsation component, and determining a first product of the first covariance and the first derivative and a second product of the second covariance and the second derivative; A shear generation term is determined according to a difference between the inverse of the first product and the second product.

[0008] In one embodiment, determining the turbulent transport term according to the pulsating components corresponding to different coordinate axis directions includes: determining a third product of the square of the first pulsation component and the third pulsation component, a fourth product of the square of the second pulsation component and the third pulsation component, and the cube of the third pulsation component; determining a second sum of the third product, the fourth product and the cube of the third pulsation component, and determining a third derivative of the second sum in the direction of the third coordinate axis; A turbulent transport term is determined based on half of the inverse of the third derivative.

[0009] In one embodiment, determining the buoyancy generation term of the atmospheric turbulence kinetic energy according to the dissipation rate, the turbulence trend term, the shear generation term, and the turbulence transport term comprises: determining a third sum of the turbulence tendency term and the dissipation rate; The shear generation term and the turbulence transport term are simultaneously subtracted from the third sum to obtain the buoyancy generation term of the atmospheric turbulence kinetic energy.

[0010] In a second aspect, an embodiment of the present application provides a device for determining a buoyancy generation term of atmospheric turbulent kinetic energy, comprising: An acquisition module is used to obtain the wind speed at several moments through a wind measuring radar within a specified time length; A first determination module is used to determine the dissipation rate and the pulsating components corresponding to the wind speed in different coordinate axis directions according to the wind speed at a plurality of moments; A second determination module is used to determine the turbulence trend term, the shear generation term and the turbulence transport term according to the pulsation components corresponding to different coordinate axis directions; The third determination module is used to determine the buoyancy generation term of the atmospheric turbulence kinetic energy according to the dissipation rate, the turbulence trend term, the shear generation term and the turbulence transport term.

[0011] In a third aspect, an embodiment of the present application provides an electronic device, comprising: a processor and a memory, wherein the memory stores instructions, and the instructions are loaded and executed by the processor to implement a method in any one of the above-mentioned embodiments.

[0012] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium, which stores a computer program. When the computer program is executed, the method in any one of the above-mentioned embodiments is implemented.

[0013] The beneficial effects of the above technical solution include at least: By obtaining the wind speed at several moments within a specified time length through a wind measuring radar, the dissipation rate and the pulsating components of the wind speed in different coordinate axis directions are determined according to the wind speed at several moments, and the turbulence trend term, shear generation term and turbulence transport term are determined according to the pulsating components corresponding to different coordinate axis directions. The use of wind measuring radar to collect wind speed is beneficial to improving the spatial resolution and continuity of the data and expanding the detection range, and determining more accurate dissipation rate, turbulence trend term, shear generation term and turbulence transport term, so that the buoyancy generation term of the atmospheric turbulent kinetic energy can be more accurately deduced based on the dissipation rate, turbulence trend term, shear generation term and turbulence transport term, thereby ensuring the accuracy of the buoyancy generation term.

[0014] The above summary is for illustrative purposes only and is not intended to be limiting in any way. In addition to the illustrative aspects, embodiments and features described above, further aspects, embodiments and features of the present application will be readily apparent by reference to the accompanying drawings and the following detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] In the accompanying drawings, unless otherwise specified, the same reference numerals throughout the multiple drawings represent the same or similar parts or elements. These drawings are not necessarily drawn to scale. It should be understood that these drawings only depict some embodiments disclosed in the present application and should not be regarded as limiting the scope of the present application.

[0016] Figure 1 A schematic flow chart of the steps of a method for determining the buoyancy generation term of atmospheric turbulence kinetic energy according to an embodiment of the present application; Figure 2 A structural block diagram of a device for determining a buoyancy generation term of atmospheric turbulent kinetic energy according to an embodiment of the present application; Figure 3 This is a structural block diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0017] In the following, only some exemplary embodiments are briefly described. As those skilled in the art will appreciate, the described embodiments may be modified in various ways without departing from the spirit or scope of the present application. Therefore, the drawings and descriptions are considered to be exemplary and non-restrictive in nature.

[0018] Reference Figure 1 , a flow chart of a method for determining the buoyancy generation term of atmospheric turbulent kinetic energy according to an embodiment of the present application is shown. The method for determining the buoyancy generation term of atmospheric turbulent kinetic energy may at least include steps S100-S400: S100: obtaining wind speeds at a number of moments by means of a wind measuring radar within a specified time period.

[0019] Optionally, a time window of a specified time length N can be set, and the wind measuring radar is installed in a place where the vertical wind speed profile can be measured. Within the specified time length N, the wind measuring radar uses a scanning method such as VAD scanning or DBS scanning to achieve wind speed measurement, and obtains the wind speed at several moments, specifically the three-dimensional wind speed W ( , ), is the vertical height, t is the time, that is, W ( , ) is the height at time t The wind speed at different vertical heights at different times can be obtained. The wind speed at each moment includes the component of the first coordinate axis (horizontally eastward), the component of the second coordinate axis (horizontally northward), and the component of the third coordinate axis (vertically upward / height direction). It should be noted that the wind measuring radar can be a Doppler radar or a wind measuring laser radar or other radar equipment that can be used to measure the wind field.

[0020] S200, determining the dissipation rate and the pulsating components of the wind speed corresponding to different coordinate axis directions according to the wind speed at a plurality of moments.

[0021] S300, determining a turbulence trend term, a shear generation term, and a turbulence transport term according to pulsation components corresponding to different coordinate axis directions.

[0022] S400, determining the buoyancy generation term of the atmospheric turbulence kinetic energy according to the dissipation rate, the turbulence trend term, the shear generation term, and the turbulence transport term.

[0023] The technical solution of the embodiment of the present application is to obtain the wind speed at several moments within a specified time length through a wind measuring radar, determine the dissipation rate and the pulsating components of the wind speed corresponding to different coordinate axis directions according to the wind speed at several moments, and determine the turbulence trend term, shear generation term and turbulence transport term according to the pulsating components corresponding to different coordinate axis directions. The use of wind measuring radar to collect wind speed is beneficial to improving the spatial resolution and continuity of data and expanding the detection range, and determining more accurate dissipation rate, turbulence trend term, shear generation term and turbulence transport term, so that the buoyancy generation term of atmospheric turbulent kinetic energy can be more accurately deduced based on the dissipation rate, turbulence trend term, shear generation term and turbulence transport term, thereby ensuring the accuracy of the buoyancy generation term.

[0024] In one embodiment, the dissipation rate is determined according to the wind speed at several moments in step S200, including steps S210-S240: S210. Determine an average wind speed based on wind speeds at a number of moments.

[0025] Optionally, the wind speeds at several moments are summed up, and then the sum is divided by the number of moments to obtain the average wind speed U.

[0026] S220, processing the wind speeds at several moments by fast Fourier transform to determine the turbulence spectrum, the first frequency, and the second frequency.

[0027] Optionally, the wind speed at several moments is processed by fast Fourier transform to determine the turbulence spectrum , first frequency and the second frequency .in, is the wave number (m -1 ), the first frequency and the second frequency The frequencies may be any two of the frequencies determined after the fast Fourier transform processing. The method and process of the fast Fourier transform processing refer to the existing method and will not be described in detail.

[0028] S230. Determine a first wave number according to the first frequency and the average wind speed, and determine a second wave number according to the second frequency and the average wind speed.

[0029] Optionally, the wave number determination formula is:

[0030] Specifically, the first frequency and the second frequency Substitute into the frequency , so that the corresponding first wave number can be obtained And the second wave .

[0031] S240, determining a first spectrum value according to the turbulence spectrum and the first wave number, determining a second spectrum value according to the turbulence spectrum and the second wave number, and determining a dissipation rate according to the first spectrum value, the second spectrum value, the first wave number, and the second wave number.

[0032] Specifically, according to the turbulence spectrum obtained , respectively, different first wave numbers And the second wave Substitute into, so as to determine the first spectrum value and the second spectrum value Then, based on the first spectrum value , the second spectrum value , the first wave number And the second wave , determine the dissipation rate , the specific formula is:

[0033] in, is the Kolmogorov constant (0.55), It is -5 / 3.

[0034] In one embodiment, step S200 determines the pulsating components of the wind speed corresponding to different coordinate axis directions according to the wind speed at several moments, including steps S250-S260, and the execution order of S250 and S210 is not limited: S250. Determine, according to wind speeds at a plurality of moments, a first component of the average wind speed on the first coordinate axis, a second component of the second coordinate axis, and a third component of the third coordinate axis in the vertical direction.

[0035] Optionally, based on the average wind speed U determined in step S210, a first component of the average wind speed U on the first coordinate axis, a second component of the second coordinate axis, and a third component of the third coordinate axis in the vertical direction may be determined.

[0036] S260. Determine the difference between the first current component and the first component of the wind speed on the first coordinate axis at the current moment, and obtain the first pulsating component corresponding to the wind speed in the direction of the first coordinate axis; determine the difference between the second current component and the second component of the wind speed on the second coordinate axis at the current moment, and obtain the second pulsating component corresponding to the wind speed in the direction of the second coordinate axis; determine the difference between the third current component and the third component of the wind speed on the third coordinate axis at the current moment, and obtain the third pulsating component corresponding to the wind speed in the direction of the third coordinate axis.

[0037] Optionally, the current moment is determined in each moment, such as the last moment, and then the component of the wind speed on the first coordinate axis (recorded as the first current component) at the current moment (for example, t moment, t=10) can be determined. ), the second component of the wind speed at the current moment on the second coordinate axis (recorded as the second current component ) and the third component of the wind speed at the current moment on the third coordinate axis (referred to as the third current component ). Then, the difference between the first current component and the first component is determined respectively to obtain the first pulsating component of the wind speed corresponding to the first coordinate axis direction , determine the difference between the second current component and the second component, and obtain the second pulsating component of the wind speed corresponding to the second coordinate axis direction , determine the difference between the third current component and the third component, and obtain the third pulsating component corresponding to the wind speed in the direction of the third coordinate axis .

[0038] In one embodiment, determining the turbulence trend item according to the pulsation components corresponding to different coordinate axis directions in step S300 includes steps S310-S320: S310: Determine a first sum of the square of the first pulsation component, the square of the second pulsation component, and the square of the third pulsation component.

[0039] Optionally, the square of the first pulsation component is determined , the square of the second pulsation component and the square of the third pulsation component The first sum of .

[0040] S320. Determine the turbulent kinetic energy based on half of the first sum, and perform time differentiation on the turbulent kinetic energy to obtain a turbulent trend term.

[0041] Optionally, the turbulent kinetic energy E is determined based on half of the first sum, i.e. , and then the turbulent kinetic energy E is derived with time to obtain the turbulent trend term .

[0042] In one embodiment, the step S300 determines the shear generation item according to the pulsation components corresponding to different coordinate axis directions, including steps S330-S350: S330, determining a first derivative of the first current component in the direction of the third coordinate axis and a second derivative of the second current component in the direction of the third coordinate axis.

[0043] Optionally, determine the first derivative of the first current component in the direction of the third coordinate axis (z direction) , and the second derivative of the second current component in the direction of the third coordinate axis .

[0044] S340, determining a first covariance between the first pulsation component and the third pulsation component and a second covariance between the second pulsation component and the third pulsation component, and determining a first product of the first covariance and the first derivative and a second product of the second covariance and the second derivative.

[0045] Optionally, the first pulsating component may be determined based on the first component of the average wind speed in the first coordinate axis, the second component of the second coordinate axis, and the third component of the third coordinate axis in the vertical direction. With the third pulsation component The first covariance and the second pulsation component With the third pulsation component The second covariance , and determine the first covariance The first product with the first derivative and the second covariance The second product with the second derivative .

[0046] S350: Determine a shear generation item according to a difference between the inverse of the first product and the second product.

[0047] Optionally, cut the generated items The calculation formula is:

[0048] In one embodiment, the turbulent transport term is determined according to the pulsating components corresponding to different coordinate axis directions in step S300, including steps S360-S380: S360. Determine a third product of the square of the first pulsation component and the third pulsation component, a fourth product of the square of the second pulsation component and the third pulsation component, and the cube of the third pulsation component.

[0049] Optionally, the square of the first pulsation component is determined With the third pulsation component The third product of , the square of the second pulsation component With the third pulsation component The fourth product of and the cube of the third pulsation component .

[0050] S370, determining a second sum of the third product, the fourth product and the cube of the third pulsating component, and determining a third derivative of the second sum in the direction of the third coordinate axis.

[0051] Optionally, a second sum of the third product, the fourth product and the cube of the third pulsation component is determined: , and determine the third derivative of the second sum in the direction of the third coordinate axis .

[0052] S380. Determine a turbulent transport term according to half of the inverse of the third derivative.

[0053] Optionally, the turbulent transport term The final calculation formula is:

[0054] Among them, there is no limit to the order in which S360, S330 and S310 are executed.

[0055] In one implementation, step S400 specifically includes: Determining the turbulence trend term and dissipation rate The third sum is obtained by subtracting the shear generation term and the turbulent transport term from the third sum , we get the buoyancy generation term of atmospheric turbulent kinetic energy .

[0056] It should be noted that the turbulent kinetic energy budget equation is:

[0057] in, The buoyancy generation term for the atmospheric turbulence kinetic energy that needs to be solved; It is a pressure transmission term. Since the static pressure pulsation of the atmosphere is very small (0.01-0.05 hPa), it is extremely difficult to measure. Some studies have tried to use microbarometers and vertical displacement sonar anemometers to calculate the pressure-velocity covariance, but the results obtained are very scattered and there is no clear trend towards stability, indicating that it can be ignored in actual operation. In practical applications, It can be considered that = 0, so the final turbulent kinetic energy balance equation can be determined as Then, based on the residual method, the buoyancy generation term of the atmospheric turbulent kinetic energy that needs to be solved can be determined :

[0058] That is, by determining the turbulence tendency term and dissipation rate The third sum is obtained by subtracting the shear generation term and the turbulent transport term from the third sum , we get the buoyancy generation term of atmospheric turbulent kinetic energy .

[0059] In the embodiment of the present application, a wind measuring radar is used to obtain three-dimensional wind speed. Based on the characteristics of the wind measuring radar with high temporal resolution, three-dimensional spatial resolution, continuous acquisition capability and long detection height, continuous, high temporal resolution, high spatial resolution and high detection range of more accurate wind speed data can be obtained, so that the determined turbulence trend item , dissipation rate , turbulent transport term and shear generation More accurately, the buoyancy generation term of atmospheric turbulence kinetic energy is then determined using the residual method , ultimately increasing the buoyancy generation factor accuracy.

[0060] Reference Figure 2 , shows a structural block diagram of a device for determining a buoyancy generation term of atmospheric turbulent kinetic energy according to an embodiment of the present application, the device may include: An acquisition module is used to obtain the wind speed at several moments through a wind measuring radar within a specified time length; A first determination module is used to determine the dissipation rate and the pulsating components of the wind speed corresponding to different coordinate axis directions according to the wind speed at a plurality of moments; A second determination module is used to determine the turbulence trend term, the shear generation term and the turbulence transport term according to the pulsation components corresponding to different coordinate axis directions; The third determination module is used to determine the buoyancy generation term of the atmospheric turbulence kinetic energy according to the dissipation rate, the turbulence trend term, the shear generation term and the turbulence transport term.

[0061] The functions of each module in the device of the embodiment of the present application can be found in the corresponding description in the above method, and will not be repeated here.

[0062] Reference Figure 3, shows a structural block diagram of an electronic device according to an embodiment of the present application, the electronic device includes: a memory 310 and a processor 320, the memory 310 stores instructions that can be run on the processor 320, the processor 320 loads and executes the instructions to implement the method for determining the buoyancy generation term of atmospheric turbulent kinetic energy in the above embodiment. The number of the memory 310 and the processor 320 can be one or more.

[0063] In one embodiment, the electronic device further includes a communication interface 330 for communicating with external devices and performing data exchange transmission. If the memory 310, the processor 320 and the communication interface 330 are implemented independently, the memory 310, the processor 320 and the communication interface 330 can be interconnected through a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 3 Only one thick line is used in the diagram, but this does not mean that there is only one bus or only one type of bus.

[0064] Optionally, in a specific implementation, if the memory 310, the processor 320 and the communication interface 330 are integrated on a chip, the memory 310, the processor 320 and the communication interface 330 can communicate with each other through an internal interface.

[0065] An embodiment of the present application provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the method for determining the buoyancy generation term of atmospheric turbulent kinetic energy provided in the above embodiment.

[0066] An embodiment of the present application also provides a chip, which includes a processor for calling and executing instructions stored in the memory from the memory, so that a communication device equipped with the chip executes the method provided by the embodiment of the present application.

[0067] An embodiment of the present application also provides a chip, including: an input interface, an output interface, a processor and a memory, wherein the input interface, the output interface, the processor and the memory are connected via an internal connection path, and the processor is used to execute the code in the memory. When the code is executed, the processor is used to execute the method provided in the embodiment of the application.

[0068] It should be understood that the processor may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc. It is worth noting that the processor may be a processor supporting the advanced RISC machines (ARM) architecture.

[0069] Further, optionally, the above-mentioned memory may include a read-only memory and a random access memory, and may also include a non-volatile random access memory. The memory may be a volatile memory or a non-volatile memory, or may include both volatile and non-volatile memories. Among them, the non-volatile memory may include a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or a flash memory. The volatile memory may include a random access memory (RAM), which is used as an external cache. By way of example but not limitation, many forms of RAM are available. For example, static RAM (SRAM), dynamic random access memory (DRAM), synchronous DRAM (SDRAM), double data rate synchronous dynamic random access memory (DDR SDRAM), enhanced synchronous dynamic random access memory (ESDRAM), synchronous link dynamic random access memory (SLDRAM) and direct rambus RAM (DR RAM).

[0070] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware or any combination thereof. When implemented using software, it can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the process or function according to the present application is generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions can be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium.

[0071] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example" or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. Moreover, the described specific features, structures, materials or characteristics may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine different embodiments or examples described in this specification and the features of different embodiments or examples, unless they are contradictory.

[0072] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Therefore, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of this application, the meaning of "plurality" is two or more, unless otherwise clearly and specifically defined.

[0073] Any process or method description in the flow chart or otherwise described herein can be understood to represent a module, fragment or portion of a code including one or more executable instructions for implementing the steps of a specific logical function or process. And the scope of the preferred embodiment of the present application includes other implementations, in which the functions may not be performed in the order shown or discussed, including in a substantially simultaneous manner or in a reverse order according to the functions involved.

[0074] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, which can be embodied in any computer-readable medium for use by an instruction execution system, apparatus or device (such as a computer-based system, a system including a processor or other system that can fetch instructions from an instruction execution system, apparatus or device and execute instructions), or used in combination with these instruction execution systems, apparatuses or devices.

[0075] It should be understood that the various parts of the present application can be implemented with hardware, software, firmware or a combination thereof. In the above embodiments, multiple steps or methods can be implemented with software or firmware stored in a memory and executed by a suitable instruction execution system. All or part of the steps of the above embodiment method can be completed by instructing the relevant hardware through a program, which can be stored in a computer-readable storage medium, and when the program is executed, it includes one of the steps of the method embodiment or a combination thereof.

[0076] In addition, each functional unit in each embodiment of the present application can be integrated into a processing module, or each unit can exist physically separately, or two or more units can be integrated into one module. The above-mentioned integrated module can be implemented in the form of hardware or in the form of a software functional module. If the above-mentioned integrated module is implemented in the form of a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium. The storage medium can be a read-only memory, a disk or an optical disk, etc.

[0077] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any technician familiar with the technical field can easily think of various changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

Claims

1. A method for determining the buoyancy generation term of atmospheric turbulent kinetic energy, characterized in that: include: Within a specified time period, the wind speed at several moments is obtained through wind radar; Determine the dissipation rate and the pulsating components of the wind speed corresponding to different coordinate axis directions according to the wind speed at a plurality of moments; According to the pulsating components corresponding to different coordinate axis directions, the turbulence trend term, the shear generation term and the turbulence transport term are determined; The buoyancy generation term of the atmospheric turbulence kinetic energy is determined according to the dissipation rate, the turbulence tendency term, the shear generation term and the turbulence transport term.

2. The method for determining the buoyancy generation term of atmospheric turbulence kinetic energy according to claim 1, characterized in that: Determining the dissipation rate according to the wind speed at several moments includes: Determining an average wind speed based on the wind speeds at a number of times; Processing the wind speed at a plurality of moments by fast Fourier transform to determine a turbulence spectrum, a first frequency, and a second frequency; Determining a first wave number according to the first frequency and the average wind speed, and determining a second wave number according to the second frequency and the average wind speed; A first spectrum value is determined according to the turbulence spectrum and the first wave number, a second spectrum value is determined according to the turbulence spectrum and the second wave number, and a dissipation rate is determined according to the first spectrum value, the second spectrum value, the first wave number, and the second wave number.

3. The method for determining the buoyancy generation term of atmospheric turbulent kinetic energy according to claim 1, characterized in that: Determining the pulsating components corresponding to the wind speed in different coordinate axis directions according to the wind speed at a plurality of moments includes: Determine, based on the wind speeds at a plurality of moments, a first component of the average wind speed on the first coordinate axis, a second component of the second coordinate axis, and a third component of the third coordinate axis in the vertical direction; Determine the difference between the first current component of the wind speed on the first coordinate axis at the current moment and the first component, and obtain the first pulsating component corresponding to the wind speed in the direction of the first coordinate axis; determine the difference between the second current component of the wind speed on the second coordinate axis at the current moment and the second component, and obtain the second pulsating component corresponding to the wind speed in the direction of the second coordinate axis; determine the difference between the third current component of the wind speed on the third coordinate axis at the current moment and the third component, and obtain the third pulsating component corresponding to the wind speed in the direction of the third coordinate axis.

4. The method for determining the buoyancy generation term of atmospheric turbulent kinetic energy according to claim 3, characterized in that: According to the pulsation components corresponding to different coordinate axis directions, the turbulence trend items are determined to include: determining a first sum of the square of the first pulsation component, the square of the second pulsation component, and the square of the third pulsation component; The turbulent kinetic energy is determined according to half of the first sum value, and the turbulent kinetic energy is time-derived to obtain a turbulent trend term.

5. The method for determining the buoyancy generation term of atmospheric turbulent kinetic energy according to claim 3, characterized in that: According to the pulsation components corresponding to different coordinate axis directions, the shear generation items are determined to include: Determine a first derivative of the first current component in the direction of the third coordinate axis and a second derivative of the second current component in the direction of the third coordinate axis; determining a first covariance of the first pulsation component and the third pulsation component and a second covariance of the second pulsation component and the third pulsation component, and determining a first product of the first covariance and the first derivative and a second product of the second covariance and the second derivative; A shear generation term is determined according to a difference between the inverse of the first product and the second product.

6. The method for determining the buoyancy generation term of atmospheric turbulent kinetic energy according to claim 3, characterized in that: According to the pulsating components corresponding to different coordinate axis directions, the turbulent transport items are determined to include: determining a third product of the square of the first pulsation component and the third pulsation component, a fourth product of the square of the second pulsation component and the third pulsation component, and the cube of the third pulsation component; determining a second sum of the third product, the fourth product and the cube of the third pulsation component, and determining a third derivative of the second sum in the direction of the third coordinate axis; A turbulent transport term is determined based on half of the inverse of the third derivative.

7. The method for determining the buoyancy generation term of atmospheric turbulent kinetic energy according to claim 1, characterized in that: Determining the buoyancy generation term of atmospheric turbulence kinetic energy according to the dissipation rate, the turbulence tendency term, the shear generation term, and the turbulence transport term comprises: determining a third sum of the turbulence tendency term and the dissipation rate; The shear generation term and the turbulence transport term are simultaneously subtracted from the third sum to obtain the buoyancy generation term of the atmospheric turbulence kinetic energy.

8. A device for determining the buoyancy generation term of atmospheric turbulent kinetic energy, characterized in that: include: An acquisition module is used to obtain the wind speed at several moments through a wind measuring radar within a specified time length; A first determination module is used to determine the dissipation rate and the pulsating components corresponding to the wind speed in different coordinate axis directions according to the wind speed at a plurality of moments; A second determination module is used to determine the turbulence trend term, the shear generation term and the turbulence transport term according to the pulsation components corresponding to different coordinate axis directions; The third determination module is used to determine the buoyancy generation term of the atmospheric turbulence kinetic energy according to the dissipation rate, the turbulence trend term, the shear generation term and the turbulence transport term.

9. An electronic device, characterized in that: include: A processor and a memory, wherein the memory stores instructions, and the instructions are loaded and executed by the processor to implement the method according to any one of claims 1 to 7.

10. A computer-readable storage medium, wherein a computer program is stored in the computer-readable storage medium, and when the computer program is executed, the method according to any one of claims 1 to 7 is implemented.

Images