Method, system, medium and equipment for simulating wind field of micro-terrain in mountainous area influenced by typhoon

Abstract

The invention provides a typhoon-influenced mountainous area micro-terrain wind field simulation method, system, medium and equipment, and the method comprises the steps: constructing a typhoon boundary database according to the actually measured typhoon wind speed and the time history data of typhoon distribution along the height; according to the actually measured typhoon wind speed in a typhoon boundary database, time history data of typhoon distributed along the height and a typhoon boundary defined by a user, adopting a CFD simulation technology for simulation to obtain the inflow wind speed of each point on an inflow boundary surface under the typhoon boundary defined by the user; constructing a micro-terrain rigid body boundary according to the measured data of the micro-terrain area of the mountainous area; and according to the inflow wind speed of each point on the inflow boundary surface under the microtopography rigid body boundary and the typhoon boundary, simulating the wind field at the monitoring point in the mountainous area microtopography area by adopting a CFD analogue simulation technology. The wind field in the coastal mountainous terrain can be accurately simulated, and basic wind load input parameters are provided for design of power transmission lines in the coastal mountainous area.

Description

Technical Field

[0001] The present invention belongs to the technical field of disaster prevention and mitigation of power transmission lines, and specifically relates to a method, system, medium and equipment for simulating micro-topography wind fields in mountainous areas affected by typhoons. Background Art

[0002] With the improvement of economic development level, a large number of expensive three-dimensional ultrasonic equipment, Doppler laser radar, etc. have been widely used in wind field measurement. A large amount of typhoon wind field measured data has also been collected at a few points arranged along the coast, providing a data source for the analysis of typhoon wind field characteristic parameters. Through the analysis of these measured data, it is found that typhoons show stronger non-stationarity compared with conventional winds. The statistical characteristics of non-stationary random processes, such as mean, variance, autocorrelation function, etc., change with time. Under general climatic conditions, it can be assumed that the statistical characteristics of wind speed are constant, which conforms to the assumption of stationary random process. However, under typhoon conditions, the statistical characteristics of wind speed will change with time and space, so the assumption of stationarity is no longer met. If the study is still conducted under the framework of stationary random process, the results will inevitably be unreliable. If the assumption of stationary random process is strictly followed, nearly half of the data will be discarded, which is a great waste.

[0003] Especially in typhoons in coastal micro-topography mountainous areas, the wind load calculation of mountain wind fields is currently generally specified by correction coefficients based on wind pressure or wind speed. However, this correction coefficient is mainly based on a relatively simple two-dimensional typical mountain model, or by obtaining near-surface wind field characteristics, wind speed acceleration ratio and other information through a scaled wind tunnel test model. This processing method requires more manpower and material resources, a longer time period, and the errors in the obtained wind field characteristic parameters, especially the near-surface characteristic parameters, are large. Therefore, in coastal mountainous terrain, the real state under the action of typhoons is difficult to be accurately simulated. Summary of the invention

[0004] In order to overcome the above-mentioned deficiencies of the prior art, the present invention proposes a method for simulating the wind field in a mountainous area affected by a typhoon, comprising:

[0005] A typhoon boundary database is constructed based on the measured typhoon wind speed and the time-course data of the typhoon's height distribution;

[0006] According to the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary, the CFD simulation technology is used to perform simulation to obtain the inflow wind speed at each point on the inflow boundary surface under the user-defined typhoon boundary;

[0007] Based on the measured data of the micro-topography area in the mountainous area, the rigid boundary of the micro-topography is constructed;

[0008] According to the inflow wind speed at each point on the inflow boundary surface under the micro-topography rigid boundary and the typhoon boundary, the wind field at the monitoring point in the mountain micro-topography area is simulated by using the CFD simulation technology.

[0009] Preferably, the step of constructing a micro-topography rigid boundary based on measured data of a micro-topography region in a mountainous area includes:

[0010] A plurality of boundary points of the selected mountain micro-topography area are obtained to obtain measured data of each boundary point, wherein the measured data includes surface elevation, longitude coordinate data and latitude coordinate data;

[0011] The micro-topography rigid body boundary is constructed based on the surface elevation, longitude coordinate data and latitude coordinate data of each boundary point.

[0012] Preferably, the step of constructing a typhoon boundary database based on the measured typhoon wind speed and the time course data of the typhoon's height distribution includes:

[0013] A statistical typhoon boundary database is constructed based on the time-varying average wind expression function, typhoon wind spectrum expression function, typhoon horizontal and vertical correlation expression function, wind profile coefficient and empirical value of turbulence intensity at set height;

[0014] Based on the measured typhoon wind speed, the time history data of the typhoon distribution along the height, and the typhoon horizontal and vertical correlation expression function, a measured typical typhoon boundary database is constructed;

[0015] Based on the typhoon wind speed and the time history data of the typhoon along the height distribution input by the user, and the typhoon horizontal and vertical correlation expression function, a user-defined typhoon boundary database is constructed;

[0016] A typhoon boundary database is determined based on the statistical typhoon boundary database, the measured typical typhoon boundary database and the user-defined typhoon boundary database.

[0017] Preferably, the simulation is performed using CFD simulation technology based on the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary to obtain the inflow wind speed at each point on the inflow boundary surface under the user-defined typhoon boundary, including:

[0018] According to the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary, the user-defined typhoon boundary is calculated using CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the typhoon boundary database;

[0019] According to the measured typical typhoon boundary database, the user-defined typhoon boundary is calculated using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the measured typical typhoon boundary database;

[0020] According to the user-defined typhoon boundary database, the user-defined typhoon boundary is calculated using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the user-defined typhoon boundary database.

[0021] Preferably, the user-defined typhoon boundary includes one or more of the following information:

[0022] The average wind speed of the typhoon at a set height and within a set time, the wind direction range and angular interval of the typhoon to be analyzed, the empirical values ​​of the typhoon's wind profile coefficient and the height turbulence intensity, the fitting values ​​of the typhoon's actual wind direction, wind profile coefficient and the turbulence intensity at a set height, the selected expression of the time-varying average wind expression function, the expression of the typhoon wind spectrum expression function, the expression of the typhoon's horizontal and vertical correlation expression function, and the wind field grid size.

[0023] Preferably, after simulating the wind field at the monitoring point in the mountain micro-topography area by using the CFD simulation technology according to the inflow wind speed at each point on the inflow boundary surface under the micro-topography rigid boundary and the typhoon boundary, the method further comprises:

[0024] Based on the analyzed wind speed and wind speed condition data query conditions within the typhoon wind direction range, the wind speed and wind speed condition data query conditions of the typical typhoon and the user-defined typhoon, the wind field at the monitoring point is analyzed and extracted in the wind field simulation results;

[0025] In the analyzed and extracted wind field, determine the wind speed whose height is within a set height interval and is a multiple of the set spacing, wherein the set height interval is a height range within which the wind field is distributed;

[0026] Based on the determined wind speed at a height within a set height range and a multiple of a set spacing, the wind load on the power transmission lines in the mountainous area within the area where the typhoon landing point is located is calculated.

[0027] Preferably, the expression of the time-varying average wind expression function is as follows:

[0028]

[0029] In the above formula, U(t) is the time-varying average wind speed of the typhoon at time t, h 1k(t) is the kth screening component of the sample to be screened at time t, where the sample to be screened is an empirical value based on the measured typhoon height turbulence intensity; let U1(t) be the non-steady wind speed, c1(t) = h 1k (t), U1(t)-c1(t)=r(t), r(t) is defined as the new sample to be screened, r k (t) is the kth screening weight of the new sample to be screened at time t.

[0030] Preferably, the expression of the typhoon wind spectrum expression function is as follows:

[0031] S(ω, t) = |A(ω, t)| 2 S(ω)

[0032] In the above formula: S(ω, t) is the typhoon wind spectrum at time t, A(ω, t) is the modulation parameter of the typhoon wind speed at time t, S(ω) is the power spectral density of the typhoon reaching the zero-mean stationary random process, t is the typhoon operation time, and ω is the angular frequency of the typhoon at time t.

[0033] Preferably, the modulation parameter of the typhoon wind speed at time t satisfies the following formula:

[0034]

[0035] In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, h is the height in the typhoon, is the average wind speed at height h in the △t time period, ω is the angular frequency of the typhoon at time t, t is the time of typhoon operation, and the △t time period is within the interval of typhoon operation time within 0-t.

[0036] Preferably, the typhoon wind spectrum at time t satisfies the following formula:

[0037]

[0038] In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, U * (△t) is the friction speed of the typhoon in the △t time period, and the △t time period is within the interval of the typhoon operation time within the time interval of 0-t; ω is the angular frequency of the typhoon at time t, and h is the height in the typhoon; is the average wind speed at height h during the time period △t.

[0039] Preferably, the expression of the typhoon horizontal and vertical correlation expression function is as follows:

[0040]

[0041] In the above formula: r ij (ω) is the spatial characteristics of the typhoon; C x , C y , C z They represent the attenuation coefficients of the typhoon from point i to point j, left and right, up and down, and front and back, respectively; is the average wind speed at the i-th point in the typhoon, is the average wind speed at the jth point in the typhoon; (x i ,y i 、z i ) is the three-dimensional coordinate of point i in the typhoon, (x j ,y j 、z j ) is the three-dimensional coordinate of point j located in the typhoon.

[0042] The present invention also proposes a system for simulating wind fields in mountainous areas affected by typhoons, the system comprising:

[0043] A typhoon boundary database construction module is used to construct a typhoon boundary database based on the measured typhoon wind speed and the time course data of the typhoon along the height distribution;

[0044] A simulation calculation module is used to perform simulation using CFD simulation technology based on the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary to obtain the inflow wind speed at each point on the inflow boundary surface under the user-defined typhoon boundary;

[0045] Micro-topography rigid boundary construction module, used to construct micro-topography rigid boundary based on measured data of micro-topography area in mountainous area;

[0046] The wind field simulation module is used to simulate the wind field at the monitoring point in the mountain micro-topography area by using the CFD simulation technology according to the inflow wind speed at each point on the inflow boundary surface under the micro-topography rigid boundary and the typhoon boundary.

[0047] Preferably, the micro-topography rigid body boundary construction module includes:

[0048] A boundary data acquisition submodule is used to select a plurality of boundary points of the mountain micro-topography area and acquire measured data of each boundary point, wherein the measured data includes surface elevation, longitude coordinate data and latitude coordinate data;

[0049] The boundary construction submodule is used to construct a micro-topography rigid body boundary based on the surface elevation, longitude coordinate data and latitude coordinate data of each boundary point.

[0050] Preferably, the typhoon boundary database construction module includes:

[0051] A statistical typhoon boundary database construction submodule is used to construct a statistical typhoon boundary database based on the time-varying average wind expression function, the typhoon wind spectrum expression function, the typhoon horizontal and vertical correlation expression function, the wind profile coefficient, and the empirical value of the set height turbulence intensity;

[0052] A submodule for constructing a measured typical typhoon boundary database, for constructing a measured typical typhoon boundary database based on the measured typhoon wind speed, the time course data of the typhoon along the height distribution, and the typhoon horizontal and vertical correlation expression function;

[0053] A user-defined typhoon boundary database construction submodule is used to construct a user-defined typhoon boundary database based on the typhoon wind speed and the time history data of the typhoon along the height distribution input by the user, and the typhoon horizontal and vertical correlation expression function;

[0054] The typhoon boundary database determination submodule is used to determine the typhoon boundary database based on the statistical typhoon boundary database, the measured typical typhoon boundary database and the user-defined typhoon boundary database.

[0055] Preferably, the simulation calculation module includes:

[0056] The first submodule of simulation calculation is used to calculate the user-defined typhoon boundary using CFD simulation technology according to the measured typhoon wind speed in the typhoon boundary database, the time course data of the typhoon along the height distribution, and the user-defined typhoon boundary, and obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the typhoon boundary database;

[0057] The second submodule of simulation calculation is used to calculate the user-defined typhoon boundary according to the measured typical typhoon boundary database by using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the measured typical typhoon boundary database;

[0058] The third submodule of simulation calculation is used to calculate the user-defined typhoon boundary according to the user-defined typhoon boundary database using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the user-defined typhoon boundary database.

[0059] Preferably, the user-defined typhoon boundary includes one or more of the following information:

[0060] The average wind speed of the typhoon at a set height and within a set time, the wind direction range and angular interval of the typhoon to be analyzed, the empirical values ​​of the typhoon's wind profile coefficient and the height turbulence intensity, the fitting values ​​of the typhoon's actual wind direction, wind profile coefficient and the turbulence intensity at a set height, the selected expression of the time-varying average wind expression function, the expression of the typhoon wind spectrum expression function, the expression of the typhoon's horizontal and vertical correlation expression function, and the wind field grid size.

[0061] Preferably, the system further comprises:

[0062] A wind field analysis and extraction module, for analyzing and extracting the wind field at the monitoring point in the wind field simulation results based on the wind speed and wind speed working condition data query conditions within the analyzed typhoon wind direction range, the wind speed and wind speed working condition data query conditions of the typical typhoon and the user-defined typhoon;

[0063] A wind speed determination module is used to determine the wind speed within a set height interval and a multiple of a set spacing in the analyzed and extracted wind field, wherein the set height interval is a height range within which the wind field is distributed;

[0064] The transmission line wind load calculation module is used to calculate the transmission line wind load in the mountainous area within the area where the typhoon landing point is located based on the wind speed determined to be within the set height range and a multiple of the set spacing.

[0065] Preferably, the expression of the time-varying average wind expression function is as follows:

[0066]

[0067] In the above formula, U(t) is the time-varying average wind speed of the typhoon at time t, h 1k (t) is the kth screening component of the sample to be screened at time t, where the sample to be screened is an empirical value based on the measured typhoon height turbulence intensity; let U1(t) be the non-steady wind speed, c1(t) = h 1k (t), U1(t)-c1(t)=r(t), r(t) is defined as the new sample to be screened, r k (t) is the kth screening weight of the new sample to be screened at time t.

[0068] Preferably, the expression of the typhoon wind spectrum expression function is as follows:

[0069] S(ω, t) = |A(ω, t) 2 S(ω)

[0070] In the above formula: S(ω, t) is the typhoon wind spectrum at time t, A(ω, t) is the modulation parameter of the typhoon wind speed at time t, S(ω) is the power spectral density of the typhoon reaching the zero-mean stationary random process, t is the typhoon operation time, and ω is the angular frequency of the typhoon at time t.

[0071] Preferably, the modulation parameter of the typhoon wind speed at time t satisfies the following formula:

[0072]

[0073] In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, h is the height in the typhoon, is the average wind speed at height h in the △t time period, ω is the angular frequency of the typhoon at time t, t is the time of typhoon operation, and the △t time period is within the interval of typhoon operation time within 0-t.

[0074] Preferably, the typhoon wind spectrum at time t satisfies the following formula:

[0075]

[0076] In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, U * (△t) is the friction speed of the typhoon in the △t time period, and the △t time period is within the interval of the typhoon operation time within the time interval of 0-t; ω is the angular frequency of the typhoon at time t, and h is the height in the typhoon; is the average wind speed at height h during the time period △t.

[0077] Preferably, the expression of the typhoon horizontal and vertical correlation expression function is as follows:

[0078]

[0079] In the above formula: r ij (ω) is the spatial characteristics of the typhoon; C x , C y , C z They represent the attenuation coefficients of the typhoon from point i to point j, left and right, up and down, and front and back, respectively; is the average wind speed at the i-th point in the typhoon, is the average wind speed at the jth point in the typhoon; (x i ,y i 、z i ) is the three-dimensional coordinate of point i in the typhoon, (x j ,y j 、z j) is the three-dimensional coordinate of point j located in the typhoon.

[0080] The present invention application also proposes an electronic device, comprising: at least one processor and a memory; the memory and the processor are connected via a bus;

[0081] The memory is used to store one or more programs;

[0082] When the one or more programs are executed by the at least one processor, the method for simulating the wind field in the mountainous area affected by a typhoon is implemented.

[0083] The present invention application also proposes a readable storage medium on which an execution program is stored. When the execution program is executed, the method for simulating the wind field in the mountainous area affected by the typhoon is implemented.

[0084] Compared with the closest prior art, the present invention has the following beneficial effects:

[0085] The present invention provides a method, system, medium and equipment for simulating wind fields in micro-topography of mountainous areas affected by typhoons. In the present invention, a typhoon boundary database is constructed according to the measured typhoon wind speed and the time-course data of the typhoon distribution along the height; according to the measured typhoon wind speed and the time-course data of the typhoon distribution along the height in the typhoon boundary database and the user-defined typhoon boundary, CFD simulation technology is used for simulation to obtain the inflow wind speed at each point on the inflow boundary surface under the user-defined typhoon boundary; according to the measured data of the micro-topography area in the mountainous area, a micro-topography rigid body boundary is constructed; according to the micro-topography rigid body boundary and the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary, CFD simulation technology is used to simulate the wind field at the monitoring point in the micro-topography area in the mountainous area. Comprehensively considering the influence of the measured typhoon and the micro-topography of the mountainous area on the wind field simulation analysis results, the wind field under the action of the typhoon in the coastal mountainous terrain can be more accurately simulated. In addition, in the present invention, by calculating the wind load of the mountainous area power transmission line near the coastal typhoon landing point, basic wind load input parameters can also be provided for the design of the coastal mountainous area power transmission line. Due to the introduction of simulation technology, the solution of the present invention can obtain wind field characteristic parameters with less manpower and material resources and a shorter time period; since the simulation is combined with the measured data of the mountainous micro-topography area, the error of the near-surface characteristic parameters is smaller. BRIEF DESCRIPTION OF THE DRAWINGS

[0086] Figure 1 The present invention provides a micro-topography wind field simulation analysis flow chart considering the impact of typhoon inflow;

[0087] Figure 2 The basic parameter setting management window of typhoon inflow boundary provided by the present invention application;

[0088] Figure 3 The micro-topography rigid body boundary generation and management interface provided by the present invention application;

[0089] Figure 4 The user specifies and adds monitoring points interface provided by the present invention application;

[0090] Figure 5 The present invention application provides a management interface for outputting wind speed simulation analysis results at a user-specified monitoring point;

[0091] Figure 6 A schematic diagram of the operation of an electronic device provided by the present invention;

[0092] Figure 7 A schematic flow chart of a method for simulating a micro-topography wind field in a mountainous area affected by a typhoon provided in the present invention;

[0093] Figure 8 This is a structural diagram of a system for simulating micro-topography wind fields in mountainous areas affected by typhoons provided in the present invention. DETAILED DESCRIPTION

[0094] The specific implementation methods of the present invention are further described in detail below with reference to the accompanying drawings.

[0095] Embodiment 1:

[0096] like Figure 1 , Figure 7 A method for simulating wind fields in mountainous areas affected by typhoons is shown in the figure, including the following contents:

[0097] Step 1: Construct a typhoon boundary database based on the measured typhoon wind speed and the time-course data of the typhoon's height distribution;

[0098] Step 2: According to the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary, a CFD (Computational Fluid Dynamics) simulation technology is used to perform simulation to obtain the inflow wind speed at each point on the inflow boundary surface under the user-defined typhoon boundary;

[0099] Step 3: Construct the micro-topography rigid boundary based on the measured data of the mountain micro-topography area;

[0100] Step 4: Based on the inflow wind speed at each point on the inflow boundary surface under the micro-topography rigid boundary and the typhoon boundary, the CFD simulation technology is used to simulate the wind field at the monitoring point in the mountain micro-topography area.

[0101] For example, the above step 1 may include the following process:

[0102] A statistical typhoon boundary database is constructed based on the time-varying average wind expression function, typhoon wind spectrum expression function, typhoon horizontal and vertical correlation expression function, wind profile coefficient and empirical value of turbulence intensity at set height;

[0103] The expression of the time-varying average wind expression function is as follows:

[0104]

[0105] In the above formula, U(t) is the time-varying average wind speed of the typhoon at time t, h 1k (t) is the kth screening component of the sample to be screened at time t, where the sample to be screened is an empirical value based on the measured typhoon height turbulence intensity; let U1(t) be the non-steady wind speed, c1(t) = h 1k (t), U1(t)-c1(t)=r(t), r(t) is defined as the new sample to be screened, r k (t) is the kth screening weight of the new sample to be screened at time t.

[0106] The expression of the typhoon wind spectrum expression function is as follows:

[0107] S(ω, t) = |A(ω, t) 2 S(ω)

[0108] In the above formula: S(ω, t) is the typhoon wind spectrum at time t, A(ω, t) is the modulation parameter of the typhoon wind speed at time t, S(ω) is the power spectral density of the typhoon reaching the zero-mean stationary random process, t is the typhoon operation time, and ω is the angular frequency of the typhoon at time t.

[0109] Wherein, the modulation parameter of the typhoon wind speed at time t satisfies the following formula:

[0110]

[0111] In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, h is the height in the typhoon, is the average wind speed at height h in the △t time period, ω is the angular frequency of the typhoon at time t, t is the time of typhoon operation, and the △t time period is within the interval of typhoon operation time within 0-t.

[0112] Wherein, the typhoon wind spectrum at time t satisfies the following formula:

[0113]

[0114] In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, U * (△t) is the friction speed of the typhoon in the △t time period, and the △t time period is within the interval of the typhoon operation time within the time interval of 0-t; ω is the angular frequency of the typhoon at time t, and h is the height in the typhoon; is the average wind speed at height h during the time period △t.

[0115] The expression of the typhoon horizontal and vertical correlation expression function is as follows:

[0116]

[0117] In the above formula: r ij (ω) is the spatial characteristics of the typhoon; C x , C y , C z They represent the attenuation coefficients of the typhoon from point i to point j, left and right, up and down, and front and back, respectively; is the average wind speed at the i-th point in the typhoon, is the average wind speed at the jth point in the typhoon; (x i ,y i 、z i ) is the three-dimensional coordinate of point i in the typhoon, (x j ,y j 、z j ) is the three-dimensional coordinate of point j located in the typhoon.

[0118] As mentioned above, based on the time-varying average wind expression function, the typhoon wind spectrum expression function, the expression functions of horizontal correlation and vertical correlation, and the wind profile coefficient, a set of statistically significant statistical typhoon boundaries can be established according to the empirical value of the turbulence intensity at a height of 10 m, and the statistical typhoon boundaries can be used as the statistical typhoon boundary database;

[0119] Based on the measured typhoon wind speed, the time history data of the typhoon distribution along the height, and the typhoon horizontal and vertical correlation expression function, a measured typical typhoon boundary database is constructed;

[0120] As mentioned above, based on the measured typhoon wind speed and the time-series data of its distribution along the height, as typical typhoon samples, the expression functions of horizontal correlation and vertical correlation are jointly called to establish multiple groups of typical typhoon boundaries with typical significance, and the typical typhoon boundaries with typical significance are used as the measured typical typhoon boundary database.

[0121] Based on the typhoon wind speed and the time history data of the typhoon along the height distribution input by the user, and the typhoon horizontal and vertical correlation expression function, a user-defined typhoon boundary database is constructed;

[0122] As mentioned above, based on the measured typhoon wind speed and its time course data along the height distribution input by the user, as typical typhoon samples, the expression functions of horizontal correlation and vertical correlation are jointly called to establish a user-defined typhoon boundary database;

[0123] A typhoon boundary database is determined based on the statistical typhoon boundary database, the measured typical typhoon boundary database and the user-defined typhoon boundary database.

[0124] For example, the above step 2 may include the following process:

[0125] According to the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary, the user-defined typhoon boundary is calculated using CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the typhoon boundary database;

[0126] As mentioned above, based on the above typhoon boundary database, set the time-varying average wind expression function, typhoon wind spectrum expression function, horizontal correlation and vertical correlation expression function selection function and window, set the 10-minute average wind speed of the typhoon at a height of 10m, set the input function and window of the typhoon wind direction range and angle interval that need to be analyzed, set the empirical value input function and window of the wind profile coefficient and the turbulence intensity at a height of 10m, and set the user-defined input grid subdivision function. Based on these parameters, the initial value calculation and initial value assignment of the inflow wind speed at each point on the inflow boundary surface under the statistical typhoon boundary are realized in the fluid simulation analysis software to realize the typhoon inflow boundary. CFD simulation technology can be realized by but not limited to fluid simulation analysis software.

[0127] According to the measured typical typhoon boundary database, the user-defined typhoon boundary is calculated using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the measured typical typhoon boundary database;

[0128] The above-mentioned method obtains the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the measured typical typhoon boundary database, specifically: setting the typical typhoon wind speed time series data reading function and window, setting the time-varying average wind expression function, typhoon wind spectrum expression function, horizontal correlation and vertical correlation expression function selection function and window, setting the typhoon actual wind direction input function and window, setting the wind profile coefficient and 10m height turbulence intensity fitting value display window, setting the user-defined input grid subdivision function. Based on these parameters, the initial value calculation and initial value assignment of the inflow wind speed at each point on the inflow boundary surface under the typical typhoon boundary are realized in the fluid simulation analysis software to realize the typhoon inflow boundary.

[0129] According to the user-defined typhoon boundary database, the user-defined typhoon boundary is calculated using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the user-defined typhoon boundary database.

[0130] The above-mentioned method obtains the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the user-defined typhoon boundary database, specifically: setting the user-defined typhoon typhoon wind speed time series data reading function and window, setting the selection function and window of the expression function of horizontal correlation and vertical correlation, setting the typhoon actual wind direction input function and window, setting the wind profile coefficient and the 10m height turbulence intensity fitting value display window, and setting the user-defined input grid subdivision function. Based on these parameters, the initial value calculation and initial value assignment of the inflow wind speed at each point on the inflow boundary surface under the user-defined typhoon boundary are realized in the fluid simulation analysis software to realize the typhoon inflow boundary.

[0131] The user-defined typhoon boundary includes one or more of the following information, as follows:

[0132] The average wind speed of the typhoon at a set height and within a set time, the wind direction range and angular interval of the typhoon to be analyzed, the empirical values ​​of the typhoon's wind profile coefficient and the height turbulence intensity, the fitting values ​​of the typhoon's actual wind direction, wind profile coefficient and the turbulence intensity at a set height, the selected expression of the time-varying average wind expression function, the expression of the typhoon wind spectrum expression function, the expression of the typhoon's horizontal and vertical correlation expression function, and the wind field grid size.

[0133] By describing the differences in typhoon inflow wind fields, it is possible to consider the strong shear, high turbulence and non-stationary characteristics of the typhoon wind field in the inflow boundary; preparations are made for combining the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary with the micro-topography rigid boundary in step 3, facilitating the combination of the two and realizing the simulation of micro-topography wind fields in mountainous areas.

[0134] For example, the above step 3 may include the following process:

[0135] A plurality of boundary points of the selected mountain micro-topography area are obtained to obtain measured data of each boundary point, wherein the measured data includes surface elevation, longitude coordinate data and latitude coordinate data;

[0136] The micro-topography rigid body boundary is constructed based on the surface elevation, longitude coordinate data and latitude coordinate data of each boundary point.

[0137] The virtual simulation of the micro-topography area in the mountainous area is realized, and the simulated micro-topography rigid boundary can be used to construct the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary.

[0138] After the above step 4, the following steps may also be included:

[0139] Based on the analyzed wind speed and wind speed condition data query conditions within the typhoon wind direction range, the wind speed and wind speed condition data query conditions of the typical typhoon and the user-defined typhoon, the wind field at the monitoring point is analyzed and extracted in the wind field simulation results;

[0140] As mentioned above, when analyzing and extracting the wind field at the monitoring point, the 10m height wind speed and wind speed condition data query function and window for each working condition within the typhoon wind direction range that has been analyzed are set, and the 10m height wind speed and wind speed condition data query function and window for typical typhoon and user-defined typhoon conditions are set, so as to complete the analysis and extraction of the wind field characteristic parameters of the micro-topography area under typhoon inflow conditions.

[0141] In the analyzed and extracted wind field, determine the wind speed whose height is within a set height interval and is a multiple of the set spacing, wherein the set height interval is a height range within which the wind field is distributed;

[0142] During the analysis and extraction, a general fluid dynamics simulation solver is used to solve the wind speed at a height of 10m at the user-specified monitoring point in the micro-topography area under the typhoon inflow boundary conditions, as well as the wind speed distribution along the height at the point. The fluid simulation analysis software can be a fluid dynamics simulation solver, thereby simulating various parameters of the typhoon.

[0143] Based on the determined wind speed at a height within a set height range and a multiple of a set spacing, the wind load on the power transmission lines in the mountainous area within the area where the typhoon landing point is located is calculated.

[0144] The simulated wind field is used to extract wind speeds within the set height range and in multiples of the set spacing, for example, between 10 meters and 400 meters and in multiples of 10 meters, to calculate the wind load on the power transmission lines in the mountainous area where the typhoon landfall is located.

[0145] By establishing a rigid boundary that can reflect the actual terrain changes in mountainous areas, the influence of surface micro-topography on the near-ground wind field is comprehensively considered under typhoon inflow conditions; by extracting the wind speed at a height of 10m between 10m and 400m in the micro-topography area under the typhoon inflow boundary and the wind speed change along the height, the wind load of the mountainous power transmission line near the coastal typhoon landing point that can reflect the comprehensive influence of typhoon characteristics and micro-topography is finally formed. This method can be used to accurately analyze the wind load of transmission lines in coastal mountainous areas and improve the safety of line wind resistance design.

[0146] In summary, the method proposed in the present invention solves the problem that the typhoon boundary is not considered or considered very generally in the micro-topography wind field simulation analysis by using the statistical typhoon boundary database, the measured typical typhoon database and the user-defined typhoon boundary data during the wind field simulation analysis. It can be compatible with more diverse typhoon samples, better reflect the variability of the typhoon wind field, and comprehensively consider the impact of typhoon and micro-topography on the wind load of the transmission line. This method is based on fluid dynamics simulation technology, and by describing the differences in the typhoon inflow wind field, it realizes the consideration of the strong shear, high turbulence and non-stationary characteristics of the typhoon wind field in the inflow boundary.

[0147] Embodiment 2:

[0148] The embodiment of the present invention provides a specific embodiment, in which a wind field simulation analysis is performed for a wind field in a mountain micro-topography where typhoons are frequent and the impact is more prominent. In order to reflect the impact of the typhoon wind field specificity on the wind field simulation analysis results, it is necessary to set the inlet boundary of the wind field simulation analysis to be able to reflect the typhoon average wind speed profile, turbulence intensity profile and time-varying characteristics. Then, the wind field under the influence of the rigid boundary of the mountain micro-topography is simulated and solved, so as to obtain the wind field simulation analysis results that can simultaneously consider the typhoon wind field characteristics and the micro-topography, and also provide basic wind load input parameters for the design of power transmission lines in coastal mountainous areas. By analyzing the wind speed data under typhoon conditions in the framework of non-stationary random processes, the information in the data can be fully utilized to obtain more accurate and reliable analysis results. This method can be used to accurately analyze wind loads on power transmission lines in coastal mountainous areas and improve the safety of wind-resistant design of lines.

[0149] The specific embodiment of the method for simulating the wind field in mountainous area micro-topography considering the influence of typhoon characteristics is as follows:

[0150] (1) Establishing a typhoon boundary database

[0151] 1) Based on the time-varying average wind expression function, typhoon wind spectrum expression function, horizontal correlation and vertical correlation expression functions, wind profile coefficient, and empirical value of turbulence intensity at 10m height, a set of statistically significant typhoon boundary databases are established.

[0152] ① The time-varying average wind expression function is obtained by iterative calculation through the following process:

[0153] The screening method for setting the kth screening component of the sample to be screened at time t is as follows: taking the empirical value of the measured turbulence intensity of the typhoon as the sample to be screened, and the sample to be screened at time t is screened for the first time to obtain the first screening component h1(t);

[0154] Usually, the component h1(t) obtained after the first screening is not a function yet and needs to be screened multiple times. The component h1(t) obtained after the first screening is screened repeatedly, and the screening is iterated k times. After k screenings, h 1k (t) is a function, definition: c1(t) = h 1k (t);

[0155] c1(t) contains the shortest period component of the non-steady wind speed U1(t). The remainder r(t) after separating c1(t) from U1(t) is:

[0156] U1(t)-c1(t)=r(t)

[0157] Because r(t) still contains components with a longer period, r(t) is taken as a new sample to be screened and subjected to the same screening process until the residue becomes very small or the residue is a monotonic function, and r is obtained. k (t);

[0158] Finally, the time-varying average wind expression function is expressed as: the kth screening component h of the sample to be screened 1k (t) and the kth screening component r of the new sample to be screened k The sum of (t) is the time-varying average wind speed history;

[0159] The expression of the time-varying average wind expression function is as follows:

[0160]

[0161] In the above formula, U(t) is the time-varying average wind speed history of the typhoon at time t;

[0162] ② The typhoon wind spectrum of the non-stationary random typhoon process is set as follows:

[0163] First, the modulation parameter of the typhoon wind speed at time t satisfies the following formula:

[0164]

[0165] In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, h is the height in the typhoon, is the average wind speed at height h in the △t time period, ω is the angular frequency of the typhoon at time t, t is the time of typhoon operation, and the △t time period is within the interval of typhoon operation time within 0-t.

[0166] Secondly, the typhoon wind spectrum at time t satisfies the following formula:

[0167]

[0168] In the above formula: U * (△t) is the friction speed of the typhoon in the △t time period, and the △t time period is within the interval of the typhoon operation time within the time interval of 0-t; ω is the angular frequency of the typhoon at time t, and h is the height in the typhoon; is the average wind speed at height h during the time period △t.

[0169] ③ Horizontal and vertical correlation expression functions, the three-dimensional expression of the spatial characteristics of typhoon is:

[0170]

[0171] Where: r ij (ω) is the spatial characteristics of the typhoon; C x , C y , C z They represent the attenuation coefficients of the typhoon from point i to point j, left and right, up and down, and front and back, respectively; is the average wind speed at the i-th point in the typhoon, is the average wind speed at the jth point in the typhoon; (x i ,y i 、z i ) is the three-dimensional coordinate of point i in the typhoon, (x j ,y j 、z j ) is the three-dimensional coordinate of point j located in the typhoon.

[0172] ④ The wind profile coefficient can be referred to the research results of relevant literature, and a safer value of 0.12 is taken for the wind profile coefficient.

[0173] ⑤ The turbulence intensity at 10m can refer to the research results of relevant literature. The turbulence intensity at 10m is taken as 1.24 times I 10 , I 10 Enter parameters for the user.

[0174] ⑥ Based on the above formula, the external input parameters and model functions of the statistical typhoon boundary with statistical significance are formed, providing a basis for further forming the statistical typhoon boundary in the software.

[0175] 2) Based on the measured typhoon wind speed and its time-course data along the height distribution, as typical typhoon samples, the expression function (1) of horizontal correlation and vertical correlation is jointly called, and through interpolation and extension, according to the boundary grid size defined by the user, multiple groups of typical typhoon boundary databases with typical significance are established.

[0176] 3) Based on the measured typhoon wind speed and its time-course data along the height distribution input by the user, as a typical typhoon sample, the expression function (1) of horizontal correlation and vertical correlation is jointly called, and through interpolation and extension, according to the user-defined boundary grid size, a user-defined typhoon boundary database is established.

[0177] (2) Realization of typhoon inflow boundary

[0178] ① Based on the statistical typhoon boundary database, set the selection function and window of the expression function of time-varying mean wind function, typhoon wind spectrum function, horizontal correlation and vertical correlation.

[0179] ② Set the 10-minute average wind speed of the typhoon at a height of 10m, and set the input function and window of the typhoon wind direction range and angle interval that need to be analyzed.

[0180] ③ Set the empirical value input function and window of wind profile coefficient and turbulence intensity at 10m height;

[0181] ④Set user-defined input grid subdivision function.

[0182] ⑤ Based on these parameters, the initial value calculation and initial value assignment of the inflow wind speed at each point on the inflow boundary surface under the statistical typhoon boundary are realized in the fluid simulation analysis software to realize the typhoon inflow boundary. The interface of the fluid simulation analysis software is as follows Figure 2 shown.

[0183] (3) Wind field simulation in micro-topography area under typhoon inflow conditions;

[0184] 1) Input the surface elevation and longitude and latitude coordinate data of the mountain micro-topography area selected by simulation analysis, establish the micro-topography rigid body boundary, and the management interface is as follows: Figure 3 As shown;

[0185] 2) With the help of a general fluid dynamics simulation solver, the wind speed at a height of 10m at the user-specified monitoring point in the micro-topography area under the typhoon inflow boundary conditions, as well as the wind speed distribution along the height at the point, are solved. The management interface of the fluid simulation analysis software is as follows: Figure 4 As shown;

[0186] 3) Set up the 10m height wind speed and wind speed condition data query function and window for each working condition within the typhoon wind direction range that has been analyzed, set up the 10m height wind speed and wind speed condition data query function and window for typical typhoon and user-defined typhoon conditions, complete the analysis and extraction of wind field characteristic parameters in the micro-topography area under typhoon inflow conditions, and the management interface of the fluid simulation analysis software is as follows: Figure 5 shown.

[0187] Embodiment 3:

[0188] like Figure 8 As shown, the present invention discloses a system for simulating wind fields in mountainous areas affected by typhoons, the system comprising:

[0189] A typhoon boundary database construction module is used to construct a typhoon boundary database based on the measured typhoon wind speed and the time course data of the typhoon along the height distribution;

[0190] A simulation calculation module is used to perform simulation using CFD simulation technology based on the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary to obtain the inflow wind speed at each point on the inflow boundary surface under the user-defined typhoon boundary;

[0191] Micro-topography rigid boundary construction module, used to construct micro-topography rigid boundary based on measured data of micro-topography area in mountainous area;

[0192] The wind field simulation module is used to simulate the wind field at the monitoring point in the mountain micro-topography area by using the CFD simulation technology according to the inflow wind speed at each point on the inflow boundary surface under the micro-topography rigid boundary and the typhoon boundary.

[0193] Furthermore, the micro-topography rigid body boundary construction module includes:

[0194] A boundary data acquisition submodule is used to select a plurality of boundary points of the mountain micro-topography area and acquire measured data of each boundary point, wherein the measured data includes surface elevation, longitude coordinate data and latitude coordinate data;

[0195] The boundary construction submodule is used to construct a micro-topography rigid body boundary based on the surface elevation, longitude coordinate data and latitude coordinate data of each boundary point.

[0196] Furthermore, the typhoon boundary database construction module includes:

[0197] A statistical typhoon boundary database construction submodule is used to construct a statistical typhoon boundary database based on the time-varying average wind expression function, the typhoon wind spectrum expression function, the typhoon horizontal and vertical correlation expression function, the wind profile coefficient, and the empirical value of the set height turbulence intensity;

[0198] A submodule for constructing a measured typical typhoon boundary database, for constructing a measured typical typhoon boundary database based on the measured typhoon wind speed, the time course data of the typhoon along the height distribution, and the typhoon horizontal and vertical correlation expression function;

[0199] A user-defined typhoon boundary database construction submodule is used to construct a user-defined typhoon boundary database based on the typhoon wind speed and the time history data of the typhoon along the height distribution input by the user, and the typhoon horizontal and vertical correlation expression function;

[0200] The typhoon boundary database determination submodule is used to determine the typhoon boundary database based on the statistical typhoon boundary database, the measured typical typhoon boundary database and the user-defined typhoon boundary database.

[0201] Furthermore, the simulation calculation module includes:

[0202] The first submodule of simulation calculation is used to calculate the user-defined typhoon boundary using CFD simulation technology according to the measured typhoon wind speed in the typhoon boundary database, the time course data of the typhoon along the height distribution, and the user-defined typhoon boundary, and obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the typhoon boundary database;

[0203] The second submodule of simulation calculation is used to calculate the user-defined typhoon boundary according to the measured typical typhoon boundary database by using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the measured typical typhoon boundary database;

[0204] The third submodule of simulation calculation is used to calculate the user-defined typhoon boundary according to the user-defined typhoon boundary database using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the user-defined typhoon boundary database.

[0205] Furthermore, the user-defined typhoon boundary includes one or more of the following information:

[0206] The average wind speed of the typhoon at a set height and within a set time, the wind direction range and angular interval of the typhoon to be analyzed, the empirical values ​​of the typhoon's wind profile coefficient and the height turbulence intensity, the fitting values ​​of the typhoon's actual wind direction, wind profile coefficient and the turbulence intensity at a set height, the selected expression of the time-varying average wind expression function, the expression of the typhoon wind spectrum expression function, the expression of the typhoon's horizontal and vertical correlation expression function, and the wind field grid size.

[0207] Furthermore, the system further comprises:

[0208] A wind field analysis and extraction module, for analyzing and extracting the wind field at the monitoring point in the wind field simulation results based on the wind speed and wind speed working condition data query conditions within the analyzed typhoon wind direction range, the wind speed and wind speed working condition data query conditions of the typical typhoon and the user-defined typhoon;

[0209] A wind speed determination module is used to determine the wind speed within a set height interval and a multiple of a set spacing in the analyzed and extracted wind field, wherein the set height interval is a height range within which the wind field is distributed;

[0210] The transmission line wind load calculation module is used to calculate the transmission line wind load in the mountainous area within the area where the typhoon landing point is located based on the wind speed determined to be within the set height range and a multiple of the set spacing.

[0211] Furthermore, the expression of the time-varying average wind expression function is as follows:

[0212]

[0213] In the above formula, U(t) is the time-varying average wind speed of the typhoon at time t, h 1k (t) is the kth screening component of the sample to be screened at time t, where the sample to be screened is an empirical value based on the measured typhoon height turbulence intensity; let U1(t) be the non-steady wind speed, c1(t) = h 1k (t), U1(t)-c1(t)=r(t), r(t) is defined as the new sample to be screened, r k (t) is the kth screening weight of the new sample to be screened at time t.

[0214] Furthermore, the expression of the typhoon wind spectrum expression function is as follows:

[0215] S(ω, t) = |A(ω, t)| 2 S(ω)

[0216] In the above formula: S(ω, t) is the typhoon wind spectrum at time t, A(ω, t) is the modulation parameter of the typhoon wind speed at time t, S(ω) is the power spectral density of the typhoon reaching the zero-mean stationary random process, t is the typhoon operation time, and ω is the angular frequency of the typhoon at time t.

[0217] Furthermore, the modulation parameter of the typhoon wind speed at time t satisfies the following formula:

[0218]

[0219] In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, h is the height in the typhoon, is the average wind speed at height h in the △t time period, ω is the angular frequency of the typhoon at time t, t is the time of typhoon operation, and the △t time period is within the interval of typhoon operation time within 0-t.

[0220] Furthermore, the typhoon wind spectrum at time t satisfies the following formula:

[0221]

[0222] In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, U * (△t) is the friction speed of the typhoon in the △t time period, and the △t time period is within the interval of the typhoon operation time within the time interval of 0-t; ω is the angular frequency of the typhoon at time t, and h is the height in the typhoon; is the average wind speed at height h during the time period △t.

[0223] Furthermore, the expression of the typhoon horizontal and vertical correlation expression function is as follows:

[0224]

[0225] In the above formula: r ij (ω) is the spatial characteristics of the typhoon; C x , C y , C z They represent the attenuation coefficients of the typhoon from point i to point j, left and right, up and down, and front and back, respectively; is the average wind speed at the i-th point in the typhoon, is the average wind speed at the jth point in the typhoon; (x i ,y i 、z i ) is the three-dimensional coordinate of point i in the typhoon, (x j ,y j 、z j ) is the three-dimensional coordinate of point j located in the typhoon.

[0226] Embodiment 4:

[0227] like Figure 6 As shown, the present invention also provides an electronic device, which may be a computer device, a single-chip device, an intelligent mobile device, etc. The electronic device in this embodiment may include a processor, a memory, a transceiver component, etc. The memory, the processor, and the transceiver component are connected via a bus; the memory may be used to store an execution program, and an exemplary execution program may include instructions; the processor is used to execute the instructions stored in the memory. The memory may also be used to store data, which may be called and / or modified when the instructions are executed.

[0228] 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 gates or transistor logic devices, discrete hardware components, etc. It is the computing core and control core of the terminal, which is suitable for implementing one or more instructions, specifically suitable for loading and executing one or more instructions in a storage medium to implement the corresponding method flow or corresponding functions, so as to implement the steps of a method for simulating a micro-topography wind field in a mountainous area affected by a typhoon in the above-mentioned embodiment.

[0229] Embodiment 5:

[0230] Based on the same inventive concept, the present invention also provides a readable storage medium, specifically an electronic device readable storage medium (Memory), which is a memory device in the electronic device for storing programs and data. It can be understood that the storage medium here can include both the built-in storage medium in the electronic device and the extended storage medium supported by the electronic device. The storage medium provides a storage space, which stores the operating system of the terminal. In addition, one or more instructions suitable for being loaded and executed by the processor are also stored in the storage space, and these instructions can be one or more execution programs (including program codes). It should be noted that the storage medium here can be a high-speed RAM memory or a non-volatile memory, such as at least one disk memory. The processor loads and executes one or more instructions stored in the storage medium, which can implement the steps of a method for simulating a micro-topography wind field in a mountainous area affected by a typhoon in the above embodiment.

[0231] Those skilled in the art will appreciate that the embodiments of the present invention application can be provided as methods, systems, or computer program products. Therefore, the present invention application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present invention application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0232] The present invention application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present invention application. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0233] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0234] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0235] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention application rather than to limit its protection scope. Although the present invention application is described in detail with reference to the above embodiments, ordinary technicians in the relevant field should understand that after reading the present invention application, those skilled in the art can still make various changes, modifications or equivalent substitutions to the specific implementation methods of the application, but these changes, modifications or equivalent substitutions are all within the protection scope of the claims to be approved.

Claims

1. A method for simulating wind fields in mountainous areas affected by typhoons, characterized in that: include: A typhoon boundary database is constructed based on the measured typhoon wind speed and the time-course data of the typhoon's height distribution; According to the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary, a computational fluid dynamics (CFD) simulation technology is used to perform simulation to obtain the inflow wind speed at each point on the inflow boundary surface under the user-defined typhoon boundary; According to the measured data of micro-topography in mountainous areas, the rigid boundary of micro-topography is constructed; According to the inflow wind speed at each point on the inflow boundary surface under the micro-topography rigid boundary and the typhoon boundary, the wind field at the monitoring point in the mountain micro-topography area is simulated by using the CFD simulation technology.

2. The method according to claim 1, characterized in that The method of constructing a micro-topography rigid body boundary based on measured data of a micro-topography area in a mountainous area includes: A plurality of boundary points of the selected mountain micro-topography area are obtained to obtain measured data of each boundary point, wherein the measured data includes surface elevation, longitude coordinate data and latitude coordinate data; The micro-topography rigid body boundary is constructed based on the surface elevation, longitude coordinate data and latitude coordinate data of each boundary point.

3. The method according to claim 1, characterized in that: The method of constructing a typhoon boundary database based on the measured typhoon wind speed and the time history data of the typhoon along the height distribution includes: A statistical typhoon boundary database is constructed based on the time-varying average wind expression function, typhoon wind spectrum expression function, typhoon horizontal and vertical correlation expression function, wind profile coefficient and empirical value of turbulence intensity at set height; Based on the measured typhoon wind speed, the time history data of the typhoon distribution along the height, and the typhoon horizontal and vertical correlation expression function, a measured typical typhoon boundary database is constructed; Based on the typhoon wind speed and the time history data of the typhoon along the height distribution input by the user, and the typhoon horizontal and vertical correlation expression function, a user-defined typhoon boundary database is constructed; A typhoon boundary database is determined based on the statistical typhoon boundary database, the measured typical typhoon boundary database and the user-defined typhoon boundary database.

4. The method according to claim 1, characterized in that: The method uses CFD simulation technology to simulate the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary to obtain the inflow wind speed at each point on the inflow boundary surface under the user-defined typhoon boundary, including: According to the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary, the user-defined typhoon boundary is calculated using CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the typhoon boundary database; According to the measured typical typhoon boundary database, the user-defined typhoon boundary is calculated using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the measured typical typhoon boundary database; According to the user-defined typhoon boundary database, the user-defined typhoon boundary is calculated using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the user-defined typhoon boundary database.

5. The method according to claim 1 or 4, characterized in that: The user-defined typhoon boundary includes one or more of the following information: The average wind speed of the typhoon at a set height and within a set time, the wind direction range and angular interval of the typhoon to be analyzed, the empirical values ​​of the typhoon's wind profile coefficient and the height turbulence intensity, the fitting values ​​of the typhoon's actual wind direction, wind profile coefficient and the turbulence intensity at a set height, the selected expression of the time-varying average wind expression function, the expression of the typhoon wind spectrum expression function, the expression of the typhoon's horizontal and vertical correlation expression function, and the wind field grid size.

6. The method according to claim 1, 4 or 5, characterized in that: After simulating the wind field at the monitoring point in the mountain micro-topography area by using the CFD simulation technology according to the inflow wind speed at each point on the inflow boundary surface under the micro-topography rigid boundary and the typhoon boundary, the method further includes: Based on the analyzed wind speed and wind speed condition data query conditions within the typhoon wind direction range, the wind speed and wind speed condition data query conditions of the typical typhoon and the user-defined typhoon, the wind field at the monitoring point is analyzed and extracted in the wind field simulation results; In the analyzed and extracted wind field, determine the wind speed whose height is within a set height interval and is a multiple of the set spacing, wherein the set height interval is a height range within which the wind field is distributed; Based on the determined wind speed at a height within a set height range and a multiple of a set spacing, the wind load on the power transmission lines in the mountainous area within the area where the typhoon landing point is located is calculated.

7. The method according to claim 3, characterized in that The expression of the time-varying average wind expression function is as follows: In the above formula, U(t) is the time-varying average wind speed of the typhoon at time t, h 1k (t) is the kth screening component of the sample to be screened at time t, where the sample to be screened is an empirical value based on the measured typhoon height turbulence intensity; let U1(t) be the non-steady wind speed, c1(t) = h 1k (t), U1(t)-c1(t)=r(t), r(t) is defined as the new sample to be screened, r k (t) is the kth screening weight of the new sample to be screened at time t.

8. The method according to claim 3, characterized in that The expression of the typhoon wind spectrum expression function is as follows: S(ω,t)=|A(ω,t) 2 S(ω) In the above formula: S(ω, t) is the typhoon wind spectrum at time t, A(ω, t) is the modulation parameter of the typhoon wind speed at time t, S(ω) is the power spectral density of the typhoon reaching the zero-mean stationary random process, t is the typhoon operation time, and ω is the angular frequency of the typhoon at time t.

9. The method according to claim 8, characterized in that The modulation parameter of the typhoon wind speed at time t satisfies the following formula: In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, h is the height in the typhoon, is the average wind speed at height h in the △t time period, ω is the angular frequency of the typhoon at time t, t is the time of typhoon operation, and the △t time period is within the interval of typhoon operation time within 0-t.

10. The method according to claim 8 or 9, characterized in that: The typhoon wind spectrum at time t satisfies the following formula: In the above formula: is the average fluctuating wind speed of the typhoon between 0-h at time t, U*(△t) is the friction speed of the typhoon in the △t time period, and the △t time period is within the interval of the typhoon operation time between 0-t; ω is the angular frequency of the typhoon at time t, and h is the height in the typhoon; is the average wind speed at height h during the time period △t.

11. The method according to claim 3, characterized in that The expression of the typhoon horizontal and vertical correlation expression function is as follows: In the above formula: r ij (ω) is the spatial characteristics of the typhoon; C x , C y , C z They represent the attenuation coefficients of the typhoon from point i to point j, left and right, up and down, and front and back, respectively; is the average wind speed at the i-th point in the typhoon, is the average wind speed at the jth point in the typhoon; (x i ,y i 、z i ) is the three-dimensional coordinate of point i in the typhoon, (x j ,y j 、z j ) is the three-dimensional coordinate of point j located in the typhoon.

12. A system for simulating wind fields in mountainous areas affected by typhoons, characterized in that: The system comprises: A typhoon boundary database construction module is used to construct a typhoon boundary database based on the measured typhoon wind speed and the time course data of the typhoon along the height distribution; A simulation calculation module is used to perform simulation using CFD simulation technology based on the measured typhoon wind speed in the typhoon boundary database, the time history data of the typhoon along the height distribution, and the user-defined typhoon boundary to obtain the inflow wind speed at each point on the inflow boundary surface under the user-defined typhoon boundary; Micro-topography rigid boundary construction module, used to construct micro-topography rigid boundary based on measured data of micro-topography area in mountainous area; The wind field simulation module is used to simulate the wind field at the monitoring point in the mountain micro-topography area by using the CFD simulation technology according to the inflow wind speed at each point on the inflow boundary surface under the micro-topography rigid boundary and the typhoon boundary.

13. The system according to claim 12, characterized in that The micro-topography rigid body boundary construction module includes: A boundary data acquisition submodule is used to select a plurality of boundary points of the mountain micro-topography area and acquire measured data of each boundary point, wherein the measured data includes surface elevation, longitude coordinate data and latitude coordinate data; The boundary construction submodule is used to construct a micro-topography rigid body boundary based on the surface elevation, longitude coordinate data and latitude coordinate data of each boundary point.

14. The system according to claim 12, characterized in that The typhoon boundary database construction module includes: A statistical typhoon boundary database construction submodule is used to construct a statistical typhoon boundary database based on the time-varying average wind expression function, the typhoon wind spectrum expression function, the typhoon horizontal and vertical correlation expression function, the wind profile coefficient, and the empirical value of the set height turbulence intensity; A submodule for constructing a measured typical typhoon boundary database, for constructing a measured typical typhoon boundary database based on the measured typhoon wind speed, the time course data of the typhoon along the height distribution, and the typhoon horizontal and vertical correlation expression function; A user-defined typhoon boundary database construction submodule is used to construct a user-defined typhoon boundary database based on the typhoon wind speed and the time history data of the typhoon along the height distribution input by the user, and the typhoon horizontal and vertical correlation expression function; The typhoon boundary database determination submodule is used to determine the typhoon boundary database based on the statistical typhoon boundary database, the measured typical typhoon boundary database and the user-defined typhoon boundary database.

15. The system according to claim 12, characterized in that The simulation calculation module comprises: The first submodule of simulation calculation is used to calculate the user-defined typhoon boundary using CFD simulation technology according to the measured typhoon wind speed in the typhoon boundary database, the time course data of the typhoon along the height distribution, and the user-defined typhoon boundary, and obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the typhoon boundary database; The second submodule of simulation calculation is used to calculate the user-defined typhoon boundary according to the measured typical typhoon boundary database by using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the measured typical typhoon boundary database; The third submodule of simulation calculation is used to calculate the user-defined typhoon boundary according to the user-defined typhoon boundary database using the CFD simulation technology to obtain the inflow wind speed at each point on the inflow boundary surface under the typhoon boundary corresponding to the user-defined typhoon boundary database.

16. The system according to claim 12 or 15, characterized in that The user-defined typhoon boundary includes one or more of the following information: The average wind speed of the typhoon at a set height and within a set time, the wind direction range and angular interval of the typhoon to be analyzed, the empirical values ​​of the typhoon's wind profile coefficient and the height turbulence intensity, the fitting values ​​of the typhoon's actual wind direction, wind profile coefficient and the turbulence intensity at a set height, the selected expression of the time-varying average wind expression function, the expression of the typhoon wind spectrum expression function, the expression of the typhoon's horizontal and vertical correlation expression function, and the wind field grid size.

17. The system according to claim 12, 15 or 16, characterized in that: The system further comprises: A wind field analysis and extraction module, for analyzing and extracting the wind field at the monitoring point in the wind field simulation results based on the wind speed and wind speed working condition data query conditions within the analyzed typhoon wind direction range, the wind speed and wind speed working condition data query conditions of the typical typhoon and the user-defined typhoon; A wind speed determination module is used to determine the wind speed within a set height interval and a multiple of a set spacing in the analyzed and extracted wind field, wherein the set height interval is a height range within which the wind field is distributed; The transmission line wind load calculation module is used to calculate the transmission line wind load in the mountainous area within the area where the typhoon landing point is located based on the wind speed determined to be within the set height range and a multiple of the set spacing.

18. The system according to claim 14, characterized in that The expression of the time-varying average wind expression function is as follows: In the above formula, U(t) is the time-varying average wind speed of the typhoon at time t, h 1k (t) is the kth screening component of the sample to be screened at time t, where the sample to be screened is an empirical value based on the measured typhoon height turbulence intensity; let U1(t) be the non-steady wind speed, c1(t) = h 1k (t), U1(t)-c1(t)=r(t), r(t) is defined as the new sample to be screened, r k (t) is the kth screening weight of the new sample to be screened at time t.

19. The system according to claim 14, characterized in that The expression of the typhoon wind spectrum expression function is as follows: S(ω,t)=|A(ω,t)| 2 S(ω) In the above formula: S(ω, t) is the typhoon wind spectrum at time t, A(ω, t) is the modulation parameter of the typhoon wind speed at time t, S(ω) is the power spectral density of the typhoon reaching the zero-mean stationary random process, t is the typhoon operation time, and ω is the angular frequency of the typhoon at time t.

20. The system according to claim 19, characterized in that The modulation parameter of the typhoon wind speed at time t satisfies the following formula: In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, h is the height in the typhoon, is the average wind speed at height h in the △t time period, ω is the angular frequency of the typhoon at time t, t is the time of typhoon operation, and the △t time period is within the interval of typhoon operation time within 0-t.

21. The system according to claim 19 or 20, characterized in that The typhoon wind spectrum at time t satisfies the following formula: In the above formula: is the average fluctuating wind speed of the typhoon between 0 and h at time t, U * (△t) is the friction speed of the typhoon in the △t time period, and the △t time period is within the interval of the typhoon operation time within the time interval of 0-t; ω is the angular frequency of the typhoon at time t, and h is the height in the typhoon; is the average wind speed at height h during the time period △t.

22. The system according to claim 14, characterized in that The expression of the typhoon horizontal and vertical correlation expression function is as follows: In the above formula: r ij (ω) is the spatial characteristics of the typhoon; C x , C y , C z They represent the attenuation coefficients of the typhoon from point i to point j, left and right, up and down, and front and back, respectively; is the average wind speed at the i-th point in the typhoon, is the average wind speed at the jth point in the typhoon; (x i ,y i 、z i ) is the three-dimensional coordinate of point i in the typhoon, (x j ,y j 、z j ) is the three-dimensional coordinate of point j located in the typhoon.

23. An electronic device, characterized in that: include: at least one processor and memory; The memory and the processor are connected via a bus; The memory is used to store one or more programs; When the one or more programs are executed by the at least one processor, the method for simulating a wind field in a mountainous area affected by a typhoon according to any one of claims 1 to 11 is implemented.

24. A readable storage medium, characterized in that: An execution program is stored thereon, and when the execution program is executed, the method for simulating the wind field in a mountainous area affected by a typhoon in any one of claims 1 to 11 is implemented.

Images