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

Abstract

The present application provides a method, apparatus, device, and storage medium for determining the buoyancy production term of atmospheric turbulent kinetic energy. The method for determining the buoyancy production term of atmospheric turbulent kinetic energy obtains the wind speeds at several moments through a wind profiler radar within a specified time length, determines the dissipation rate and the pulsating components corresponding to the wind speed in different coordinate axes directions according to the wind speeds at several moments, determines the turbulent trend term, the shear production term, and the turbulent transport term according to the pulsating components corresponding to different coordinate axes directions. Using a wind profiler radar to collect wind speeds is beneficial to improving the spatial resolution, data continuity, and expanding the detection range of data, determining more accurate dissipation rate, turbulent trend term, shear production term, and turbulent transport term, so as to be able to more accurately deduce the buoyancy production term of atmospheric turbulent kinetic energy according to the dissipation rate, turbulent trend term, shear production term, and turbulent transport term, ensuring the accuracy of the buoyancy production term.

Description

Technical Field

[0001] This application relates to the field of meteorological observation, and particularly to a method, device, equipment and storage medium for determining the buoyancy production term of atmospheric turbulent kinetic energy. Background Art

[0002] Currently, in meteorological monitoring, the buoyancy production term of turbulent kinetic energy, as a key physical quantity describing the generation and dissipation process of turbulence, is crucial for revealing the essence and characteristics of turbulence. The monitoring accuracy of the buoyancy production term of turbulent kinetic energy will affect the accuracy of meteorological monitoring. At present, on the one hand, data is usually collected by means of meteorological towers or drones, with low spatial resolution, poor data continuity and limited detection range; on the other hand, the buoyancy production term is directly calculated by collecting data such as temperature. Due to the above-mentioned disadvantages of the data, it is easy to affect the accuracy of the buoyancy production term. Summary of the Invention

[0003] Embodiments of this application provide a method, device, equipment and storage medium for determining the buoyancy production term of atmospheric turbulent kinetic energy to solve at least one problem existing in the related art. The technical solutions are as follows:

[0004] In a first aspect, embodiments of this application provide a method for determining the buoyancy production term of atmospheric turbulent kinetic energy, including:

[0005] Obtaining wind speeds at several moments through a wind profiler radar within a specified time length;

[0006] Determining the dissipation rate and the pulsation components corresponding to the wind speed in different coordinate axis directions according to the wind speeds at several moments;

[0007] Determining the turbulent tendency term, the shear production term and the turbulent transport term according to the pulsation components corresponding to different coordinate axis directions;

[0008] Determining the buoyancy production term of atmospheric turbulent kinetic energy according to the dissipation rate, the turbulent tendency term, the shear production term and the turbulent transport term.

[0009] In an implementation manner, the determining the dissipation rate according to the wind speeds at several moments includes:

[0010] Determining the average wind speed according to the wind speeds at several moments;

[0011] Processing the wind speeds at several moments through fast Fourier transform to determine the turbulence spectrum, the first frequency and the second frequency;

[0012] Determining the first wave number according to the first frequency and the average wind speed, and determining the second wave number according to the second frequency and the average wind speed;

[0013] Determine a first spectral value according to the turbulence spectrum and the first wave number, determine a second spectral value according to the turbulence spectrum and the second wave number, and determine a dissipation rate according to the first spectral value, the second spectral value, the first wave number, and the second wave number.

[0014] In one embodiment, determining the pulsating components corresponding to the wind speed in different coordinate axes according to the wind speed at several moments includes:

[0015] Determine a first component of the average wind speed in a first coordinate axis, a second component in a second coordinate axis, and a third component in a third coordinate axis perpendicular thereto according to the wind speed at several moments;

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

[0017] In one embodiment, determining a turbulence trend term according to the pulsating components corresponding to different coordinate axes includes:

[0018] Determine a first sum value of the square of the first pulsating component, the square of the second pulsating component, and the square of the third pulsating component;

[0019] Determine the turbulent kinetic energy according to half of the first sum value, and perform a time derivative on the turbulent kinetic energy to obtain a turbulence trend term.

[0020] In one embodiment, determining a shear production term according to the pulsating components corresponding to different coordinate axes includes:

[0021] Determine a first derivative of the first current component in the third coordinate axis direction and a second derivative of the second current component in the third coordinate axis direction;

[0022] Determine a first covariance between the first pulsating component and the third pulsating component and a second covariance between the second pulsating component and the third pulsating component, and determine a first product of the first covariance and the first derivative and a second product of the second covariance and the second derivative;

[0023] Determine a shear production term according to a difference between the negative value of the first product and the second product.

[0024] In one embodiment, determining the turbulent transport term according to the pulsation components corresponding to different coordinate axis directions includes:

[0025] Determine the third product of the square of the first pulsation component and the third pulsation component, the fourth product of the square of the second pulsation component and the third pulsation component, and the cube of the third pulsation component;

[0026] Determine the second sum value of the third product, the fourth product and the cube of the third pulsation component, and determine the third derivative of the second sum value in the direction of the third coordinate axis;

[0027] Determine the turbulent transport term according to half of the opposite number of the third derivative.

[0028] In one embodiment, determining the buoyancy generation term of the atmospheric turbulent kinetic energy according to the dissipation rate, the turbulent trend term, the shear production term and the turbulent transport term includes:

[0029] Determine the third sum value of the turbulent trend term and the dissipation rate;

[0030] Subtract the shear production term and the turbulent transport term from the third sum value at the same time to obtain the buoyancy generation term of the atmospheric turbulent kinetic energy.

[0031] In a second aspect, an embodiment of the present application provides a device for determining the buoyancy generation term of atmospheric turbulent kinetic energy, including:

[0032] An acquisition module, configured to obtain wind speeds at several moments through a wind profiler radar within a specified time length;

[0033] A first determination module, configured to determine the dissipation rate and the pulsation components corresponding to the wind speed in different coordinate axis directions according to the wind speeds at several moments;

[0034] A second determination module, configured to determine the turbulent trend term, the shear production term and the turbulent transport term according to the pulsation components corresponding to different coordinate axis directions;

[0035] A third determination module, configured to determine the buoyancy generation term of the atmospheric turbulent kinetic energy according to the dissipation rate, the turbulent trend term, the shear production term and the turbulent transport term.

[0036] In a third aspect, an embodiment of the present application provides an electronic device, including: a processor and a memory, where instructions are stored in the memory, and the instructions are loaded and executed by the processor to implement the method in any one of the above aspects.

[0037] Fourthly, an embodiment of the present application provides a computer-readable storage medium storing a computer program, which when executed implements the method in any of the above aspects.

[0038] The beneficial effects in the above technical solutions at least include:

[0039] By obtaining the wind speeds at several moments through a wind measurement radar within a specified time length, determining the dissipation rate and the pulsation components corresponding to the wind speed in different coordinate axis directions according to the wind speeds at several moments, and determining the turbulent trend term, the shear production term, and the turbulent transport term according to the pulsation components corresponding to different coordinate axis directions, using a wind measurement radar to collect wind speeds is beneficial to improving the spatial resolution, data continuity of the data, and expanding the detection range, determining more accurate dissipation rate, turbulent trend term, shear production term, and turbulent transport term, so that the buoyancy production term of the atmospheric turbulent kinetic energy can be more accurately deduced according to the dissipation rate, turbulent trend term, shear production term, and turbulent transport term, ensuring the accuracy of the buoyancy production term.

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

[0041] In the drawings, unless otherwise specified, the same reference numerals throughout the several views denote the same or similar components 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.

[0042] Figure 1 It is a schematic flowchart of the steps of a method for determining the buoyancy production term of the atmospheric turbulent kinetic energy according to an embodiment of the present application;

[0043] Figure 2 It is a structural block diagram of a device for determining the buoyancy production term of the atmospheric turbulent kinetic energy according to an embodiment of the present application;

[0044] Figure 3 It is a structural block diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0045] In the following, only some exemplary embodiments are briefly described. As those skilled in the art can recognize, the described embodiments can be modified in various different ways without departing from the spirit or scope of the present application. Therefore, the drawings and the description are considered to be exemplary in nature and not restrictive.

[0046] Reference Figure 1 , a flowchart 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:

[0047] S100. Within a specified time length, obtain the wind speeds at several moments through a wind profiler radar.

[0048] Optionally, a time window with a specified time length N can be set. Install the wind profiler radar at a place where the vertical wind speed profile can be measured. Within the specified time length N, use scanning methods such as VAD scanning or DBS scanning that can achieve wind speed measurement through the wind profiler radar to obtain the wind speeds at several moments. Specifically, the three-dimensional wind speed W( , ), is the vertical height, t is the moment, that is, W( , ) is the wind speed at height at moment t. The wind speeds at different vertical heights at different moments can be obtained. Among them, the wind speed at each moment includes the component along the first coordinate axis (horizontal eastward), the component along the second coordinate axis (horizontal northward), and the component along the third coordinate axis (vertical upward / height direction). It should be noted that the wind profiler radar can be a Doppler radar or a wind profiling lidar and other radar devices that can be used to measure the wind field.

[0049] S200. According to the wind speeds at several moments, determine the dissipation rate and the pulsating components corresponding to the wind speed in different coordinate axis directions.

[0050] S300. According to the pulsating components corresponding to different coordinate axis directions, determine the turbulent tendency term, the shear production term, and the turbulent transport term.

[0051] S400. According to the dissipation rate, the turbulent tendency term, the shear production term, and the turbulent transport term, determine the buoyancy generation term of atmospheric turbulent kinetic energy.

[0052] The technical solution of the embodiment of the present application, by obtaining the wind speeds at several moments through a wind profiler radar within a specified time length, determining the dissipation rate and the pulsating components corresponding to the wind speed in different coordinate axis directions according to the wind speeds at several moments, determining the turbulent tendency term, the shear production term, and the turbulent transport term according to the pulsating components corresponding to different coordinate axis directions, using a wind profiler radar to collect wind speeds is beneficial to improving the spatial resolution, data continuity of the data, and expanding the detection range, determining more accurate dissipation rate, turbulent tendency term, shear production term, and turbulent transport term, so as to be able to more accurately deduce the buoyancy generation term of atmospheric turbulent kinetic energy according to the dissipation rate, turbulent tendency term, shear production term, and turbulent transport term, and ensure the accuracy of the buoyancy generation term.

[0053] In one embodiment, determining the dissipation rate according to the wind speeds at several moments in step S200 includes steps S210 - S240:

[0054] S210. Determine the average wind speed according to the wind speeds at several moments.

[0055] Optionally, sum the wind speeds at several moments, and then divide the sum result by the number of moments to obtain the average wind speed U.

[0056] S220. Process the wind speeds at several moments through fast Fourier transform to determine the turbulence spectrum, the first frequency, and the second frequency.

[0057] Optionally, process the wind speeds at several moments through fast Fourier transform to determine the turbulence spectrum , the first frequency and the second frequency . Wherein, is the wave number (m -1 ), the first frequency and the second frequency can be any two frequencies among the frequencies determined after the fast Fourier transform processing. The method process of the fast Fourier transform processing refers to the existing method and will not be elaborated here.

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

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

[0060]

[0061] Specifically, substitute the first frequency and the second frequency into the frequency respectively, so as to obtain the corresponding first wave number and the second wave number .

[0062] S240. Determine the first spectral value according to the turbulence spectrum and the first wave number, determine the second spectral value according to the turbulence spectrum and the second wave number, and determine the dissipation rate according to the first spectral value, the second spectral value, the first wave number, and the second wave number.

[0063] Specifically, according to the obtained turbulence spectrum , substitute different first wave numbers and second wave numbers respectively, so as to determine the first spectral value and the second spectral value . Then, based on the first spectral value , the second spectral value , the first wave number and the second wave number , the dissipation rate is determined. The specific formula is:

[0064]

[0065] where is the Kolmogorov constant (0.55), is -5 / 3.

[0066] In one embodiment, in step S200, according to the wind speeds at several moments, the pulsating components corresponding to the wind speed in different coordinate axes are determined, including steps S250 - S260. The execution order of S250 and S210 is not limited:

[0067] S250. According to the wind speeds at several moments, determine the first component of the average wind speed in the first coordinate axis, the second component in the second coordinate axis, and the third component in the third coordinate axis in the vertical direction.

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

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

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

[0071] In one implementation, in step S300, according to the pulsating components corresponding to different coordinate axis directions, determining the turbulence trend term includes steps S310 - S320:

[0072] S310. Determine the first sum value of the square of the first pulsating component, the square of the second pulsating component, and the square of the third pulsating component.

[0073] Optionally, determine the square of the first pulsating component , the square of the second pulsating component and the square of the third pulsating component of the first sum value .

[0074] S320. According to half of the first sum value, determine the turbulent kinetic energy, and perform time differentiation on the turbulent kinetic energy to obtain the turbulence trend term.

[0075] Optionally, according to half of the first sum value, determine the turbulent kinetic energy E, that is , and then perform time differentiation on the turbulent kinetic energy E to obtain the turbulence trend term .

[0076] In one implementation, in step S300, according to the pulsating components corresponding to different coordinate axis directions, determining the shear production term includes steps S330 - S350:

[0077] S330. Determine the first derivative of the first current component in the third coordinate axis direction and the second derivative of the second current component in the third coordinate axis direction.

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

[0079] S340. Determine the first covariance between the first pulsating component and the third pulsating component and the second covariance between the second pulsating component and the third pulsating component, and determine the first product of the first covariance and the first derivative and the second product of the second covariance and the second derivative.

[0080] Optionally, based on the first component of the mean wind speed in the first coordinate axis, the second component in the second coordinate axis, and the third component in the third coordinate axis in the vertical direction, determine the first pulsating component and the third pulsating component of the first covariance and a second pulsation component and a third pulsation component of the second covariance , and determine a first covariance of the first product of the first derivative and a second covariance of the second product of the second derivative .

[0081] S350. Determine a shear generation term according to the difference between the negative of the first product and the second product.

[0082] Optionally, the calculation formula of the shear generation term is as follows:

[0083] In one implementation, in step S300, according to the pulsation components corresponding to different coordinate axis directions, determine a turbulent transport term, including steps S360 - S380:

[0084] S360. Determine the third product of the square of the first pulsation component and the third pulsation component, the fourth product of the square of the second pulsation component and the third pulsation component, and the cube of the third pulsation component.

[0085] Optionally, determine the square of the first pulsation component and the third pulsation component of the third product , the square of the second pulsation component and the third pulsation component of the fourth product and the cube of the third pulsation component .

[0086] S370. Determine the second sum of the third product, the fourth product and the cube of the third pulsation component, and determine the third derivative of the second sum in the third coordinate axis direction.

[0087] Optionally, determine the second sum of the third product, the fourth product and the cube of the third pulsation component , and determine the third derivative of the second sum in the third coordinate axis direction .

[0088] S380. Determine the turbulent transport term according to half of the negative of the third derivative.

[0089] Optionally, the final calculation formula of the turbulent transport term is as follows:

[0090]

[0091] Among them, there is no specified order of execution between S360, S330, and S310.

[0092] In one embodiment, step S400 specifically includes:

[0093] Determine the turbulent trend term and the dissipation rate of the third sum value, and subtract the shear production term and the turbulent transport term from the third sum value at the same time , to obtain the buoyancy production term of the atmospheric turbulent kinetic energy .

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

[0095]

[0096] Among them, is the buoyancy production term of the atmospheric turbulent kinetic energy to be solved; is the pressure transport term. Since the static pressure pulsation of the atmosphere is very small (0.01 - 0.05 hPa) and extremely difficult to measure, existing research has tried to use a microbarometer and a sonar anemometer with vertical displacement 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. Therefore, the final turbulent kinetic energy budget equation can be determined as , and then based on the residual method, the buoyancy production term of the atmospheric turbulent kinetic energy to be solved can be determined :

[0097]

[0098] That is to say, by determining the turbulent trend term and the dissipation rate of the third sum value, and subtracting the shear production term and the turbulent transport term from the third sum value at the same time , to obtain the buoyancy production term of the atmospheric turbulent kinetic energy .

[0099] In the embodiment of the present application, a wind profiler radar is used to obtain three-dimensional wind speeds. Based on the characteristics of the wind profiler radar having high time resolution, three-dimensional spatial resolution, continuous acquisition ability, and a far detection altitude, more accurate wind speed data with high continuity, high time resolution, high spatial resolution, and a high detection range can be obtained. As a result, the determined turbulent trend term , dissipation rate , turbulent transport term and shear production term are more accurate, and then the residual method is used to determine the buoyancy production term of the atmospheric turbulent kinetic energy , finally improve the buoyancy production term accuracy.

[0100] Refer to Figure 2 , which shows the structural block diagram of the buoyancy production term determination device for atmospheric turbulent kinetic energy according to an embodiment of the present application. The device may include:

[0101] An acquisition module, configured to obtain wind speeds at several moments through a wind measurement radar within a specified time length;

[0102] A first determination module, configured to determine the dissipation rate and the pulsating components corresponding to the wind speed in different coordinate axis directions according to the wind speeds at several moments;

[0103] A second determination module, configured to determine the turbulent trend term, the shear production term, and the turbulent transport term according to the pulsating components corresponding to different coordinate axis directions;

[0104] A third determination module, configured to determine the buoyancy production term of the atmospheric turbulent kinetic energy according to the dissipation rate, the turbulent trend term, the shear production term, and the turbulent transport term.

[0105] For the functions of the modules in the device according to the embodiment of the present application, reference may be made to the corresponding descriptions in the above method, which will not be elaborated here.

[0106] Refer to Figure 3 , which shows the 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. Instructions that can run on the processor 320 are stored in the memory 310. The processor 320 loads and executes the instructions to implement the method for determining the buoyancy production term of the atmospheric turbulent kinetic energy in the above embodiment. Among them, the number of the memory 310 and the processor 320 may be one or more.

[0107] In one implementation, the electronic device further includes a communication interface 330, configured to communicate with external devices and perform data interaction and 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 may be interconnected through a bus and complete communication with each other. The bus may be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 3It is represented only by a thick line, but it does not mean that there is only one bus or one type of bus.

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

[0109] The embodiment of the present application provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, it implements the method for determining the buoyancy generation term of atmospheric turbulence kinetic energy provided in the above embodiment.

[0110] The embodiment of the present application further provides a chip, which includes a processor for calling and running instructions stored in a memory, so that a communication device equipped with the chip executes the method provided in the embodiment of the present application.

[0111] The embodiment of the present application further provides a chip, including: an input interface, an output interface, a processor, and a memory. The input interface, the output interface, the processor, and the memory are connected through an internal connection path. 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.

[0112] It should be understood that the above-mentioned processor may be a central processing unit (CPU), or may also be 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. The general-purpose processor may be a microprocessor or any conventional processor, etc. It is worth noting that the processor may be a processor that supports the advanced RISC machines (ARM) architecture.

[0113] Further, optionally, the above-mentioned memory may include a read-only memory and a random access memory, and may further 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 random access memory (SRAM), dynamic random access memory (DRAM), synchronous dynamic random access memory (SDRAM), double data rate synchronous dynamic random access memory (DDR SDRAM), enhanced synchronous dynamic random access memory (ESDRAM), synchlink dynamic random access memory (SLDRAM), and direct rambus random access memory (DR RAM).

[0114] In the above embodiments, it may be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented using software, it may 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 processes or functions according to the present application are generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions may be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium.

[0115] In the description of this specification, the descriptions referring to terms such as "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. mean that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of this application. Moreover, the specific features, structures, materials, or characteristics described can be combined in a suitable manner in any one or more embodiments or examples. In addition, without contradiction, those skilled in the art can combine and combine the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0116] In addition, the terms "first" and "second" are only used for descriptive purposes and cannot be understood as indicating or implying relative importance or implicitly specifying the quantity of the indicated technical features. Thus, the features defined with "first" and "second" can explicitly or implicitly include at least one of such features. In the description of this application, "a plurality of" means two or more unless otherwise specifically defined.

[0117] Any process or method description represented in the flowchart or described in other ways herein can be understood as representing a module, segment, or part of code including one or more executable instructions for implementing a specific logical function or process. And the scope of the preferred embodiments of this application includes additional implementations, where the functions can be executed in a substantially simultaneous manner or in a reverse order according to the involved functions, rather than in the order shown or discussed.

[0118] The logic and / or steps represented in the flowchart or described in other ways herein, for example, can be considered as a sequenced list of executable instructions for implementing a logical function and can be specifically implemented 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 systems that can fetch and execute instructions from the instruction execution system, apparatus, or device), or in connection with these instruction execution systems, apparatuses, or devices.

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

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

[0121] As described above, the above are only specific embodiments of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art within the technical scope disclosed in the present application can easily think of various changes or substitutions thereof, and these should all be covered by the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims

1. A method for determining the buoyancy generation term of atmospheric turbulent kinetic energy, characterized in that Comprising: Within a specified time length, obtain the wind speeds at several moments by using a VAD scan or a DBS scan with a wind profiler radar, where the wind speeds are three-dimensional wind speeds; Determine the dissipation rate and the pulsation components corresponding to the wind speeds in different coordinate axis directions according to the wind speeds at several moments; Determine the turbulent trend term, the shear production term, and the turbulent transport term according to the pulsation components corresponding to different coordinate axis directions; Determine the buoyancy production term of the atmospheric turbulent kinetic energy according to the dissipation rate, the turbulent trend term, the shear production term, and the turbulent transport term. Specifically: determine the third sum value of the turbulent trend term and the dissipation rate, and subtract the shear production term and the turbulent transport term from the third sum value simultaneously to obtain the buoyancy production term of the atmospheric turbulent kinetic energy; The determining the dissipation rate according to the wind speeds at several moments includes: Determine the average wind speed according to the wind speeds at several moments; Process the wind speeds at several moments through fast Fourier transform to determine the turbulent spectrum, the first frequency, and the second frequency; Determine the first wave number according to the first frequency and the average wind speed, and determine the second wave number according to the second frequency and the average wind speed. The specific formulas are: where U is the average wind speed, substitute the first frequency f1 and the second frequency f0 into the frequency f respectively, and the corresponding wave numbers k are the first wave number k1 and the second wave number k0; Determine the first spectrum value according to the turbulent spectrum and the first wave number, determine the second spectrum value according to the turbulent spectrum and the second wave number, and determine the dissipation rate according to the first spectrum value, the second spectrum value, the first wave number, and the second wave number.

2. The method for determining the buoyancy generation term of atmospheric turbulent kinetic energy according to claim 1, wherein: Determine the pulsation components corresponding to the wind speeds in different coordinate axis directions according to the wind speeds at several moments includes: Determine the first component of the average wind speed in the first coordinate axis, the second component in the second coordinate axis, and the third component in the third coordinate axis perpendicular to the vertical direction according to the wind speeds at several moments; Determine the difference between the first current component of the wind speed at the current moment in the first coordinate axis and the first component to obtain the first pulsation component corresponding to the wind speed in the first coordinate axis direction, determine the difference between the second current component of the wind speed at the current moment in the second coordinate axis and the second component to obtain the second pulsation component corresponding to the wind speed in the second coordinate axis direction, and determine the difference between the third current component of the wind speed at the current moment in the third coordinate axis and the third component to obtain the third pulsation component corresponding to the wind speed in the third coordinate axis direction.

3. The method for determining the buoyancy generation term of atmospheric turbulent kinetic energy according to claim 2, wherein: Determine the turbulent trend term according to the pulsation components corresponding to different coordinate axis directions includes: Determine the first sum value of the square of the first pulsation component, the square of the second pulsation component, and the square of the third pulsation component; Determine the turbulent kinetic energy according to half of the first sum value, and perform time differentiation on the turbulent kinetic energy to obtain the turbulent trend term.

4. The method for determining the buoyancy generation term of atmospheric turbulent kinetic energy according to claim 2, wherein: Determine the shear production term according to the pulsation components corresponding to different coordinate axis directions includes: Determine the first derivative of the first current component in the third coordinate axis direction and the second derivative of the second current component in the third coordinate axis direction; Determine the first covariance between the first pulsation component and the third pulsation component and the second covariance between the second pulsation component and the third pulsation component, and determine the first product of the first covariance and the first derivative and the second product of the second covariance and the second derivative; Determine the shear production term according to the difference between the negative value of the first product and the second product.

5. The method for determining the buoyancy generation term of atmospheric turbulent kinetic energy according to claim 2, wherein: According to the pulsation components corresponding to different coordinate axes directions, determine that the turbulent transport term includes: Determine the third product of the square of the first pulsation component and the third pulsation component, the fourth product of the square of the second pulsation component and the third pulsation component, and the cube of the third pulsation component; Determine the second sum value of the third product, the fourth product and the cube of the third pulsation component, and determine the third derivative of the second sum value in the third coordinate axis direction; Determine the turbulent transport term according to half of the negative value of the third derivative.

6. An apparatus for determining the buoyancy generation term of atmospheric turbulent kinetic energy, characterized in that, Include: An acquisition module, configured to obtain the wind speeds at several moments within a specified time length by using a wind profiler radar in the VAD scanning or DBS scanning mode, and the wind speeds are three-dimensional wind speeds; A first determination module, configured to determine the dissipation rate and the pulsation components corresponding to the wind speeds in different coordinate axis directions according to the wind speeds at several moments; A second determination module, configured to determine the turbulent tendency term, the shear production term and the turbulent transport term according to the pulsation components corresponding to different coordinate axis directions; A third determination module, configured to determine the buoyancy production term of the atmospheric turbulent kinetic energy according to the dissipation rate, the turbulent tendency term, the shear production term and the turbulent transport term. Specifically: determine the third sum value of the turbulent tendency term and the dissipation rate, and subtract the shear production term and the turbulent transport term from the third sum value at the same time to obtain the buoyancy production term of the atmospheric turbulent kinetic energy; The determining the dissipation rate according to the wind speeds at several moments includes: Determine the average wind speed according to the wind speeds at several moments; Process the wind speeds at several moments through fast Fourier transform to determine the turbulent spectrum, the first frequency and the second frequency; Determine the first wave number according to the first frequency and the average wind speed, and determine the second wave number according to the second frequency and the average wind speed. The specific formula is: Where U is the average wind speed, substitute the first frequency f1 and the second frequency f0 into the frequency f respectively, and the corresponding wave numbers k are the first wave number k1 and the second wave number k0; Determine the first spectrum value according to the turbulent spectrum and the first wave number, determine the second spectrum value according to the turbulent spectrum and the second wave number, and determine the dissipation rate according to the first spectrum value, the second spectrum value, the first wave number and the second wave number.

7. An electronic device, characterized in that, Include: A processor and a memory, instructions are stored in the memory, and the instructions are loaded and executed by the processor to implement the method according to any one of claims 1-5.

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

Images