CN108399507B - Typhoon disaster influence assessment method and device - Google Patents
Typhoon disaster influence assessment method and device Download PDFInfo
- Publication number
- CN108399507B CN108399507B CN201810239480.5A CN201810239480A CN108399507B CN 108399507 B CN108399507 B CN 108399507B CN 201810239480 A CN201810239480 A CN 201810239480A CN 108399507 B CN108399507 B CN 108399507B
- Authority
- CN
- China
- Prior art keywords
- wind speed
- direction angle
- wind
- shear
- wind direction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000001066 destructive effect Effects 0.000 claims abstract description 21
- 238000011156 evaluation Methods 0.000 claims abstract description 18
- 239000013598 vector Substances 0.000 claims description 39
- 238000012545 processing Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 description 11
- 238000004088 simulation Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000005034 decoration Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 208000031968 Cadaver Diseases 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Computer Security & Cryptography (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a method and a device for evaluating typhoon disaster-causing influence. The method comprises the following steps: obtaining lattice point wind field parameters output by simulating a typhoon wind field by using a mesoscale numerical mode; for each lattice point meeting the preset requirement in the same horizontal plane: calculating destructive force parameters of the grid point, wherein the destructive force parameters comprise any one or more of wind speed space shear quantity, wind speed time shear quantity, wind direction angle space shear quantity and wind direction angle time shear quantity; and based on any one or more of the obtained wind speed space shear variable, wind speed time shear variable, wind direction angle space shear variable and wind direction angle time shear variable of each grid point, and combining the wind speed value, thereby evaluating the typhoon disaster-causing influence. The method can comprehensively and accurately reflect the parameter characteristic difference of the destructive force among all parts in the engineering influenced by the typhoon, and ensures the evaluation accuracy of the typhoon disaster-causing influence.
Description
Technical Field
The invention relates to the technical field of wind engineering, in particular to a method and a device for evaluating typhoon disaster-causing influence.
Background
The strong wind disaster is a focus of attention of various industries in the economic society, and is a key point of research in various fields of disaster prevention and reduction. Wind resistance specifications with national/local professional characteristics are established in various professional fields such as building structures, transportation, energy and power, city planning and the like.
In the current wind-resistant design at home and abroad, the maximum value of the wind-resistant parameter of the engineering area is often selected to evaluate the influence of strong wind by point substitution. However, the range of modern large-scale projects, such as large-span and ultra-long bridges, large wind farms and the like, covers dozens or even hundreds of square kilometers, and such large-scale projects are often built on complicated terrain underlying surfaces such as sea-land or island junctions, deep-cut canyons and the like, so that the wind-resistant parameter characteristic difference of the same project site is very obvious. Particularly, under the influence of typhoon, due to the vortex structure of a typhoon wind field, the wind speed and the wind direction and related characteristic quantity thereof change rapidly in space and time, and the situation of wind resistance parameters is more complicated due to the special terrain and the complex underlying surface. In addition, researches show that the asymmetry of the typhoon structure can be increased by the complex terrains and the underlying surfaces, and for example, systems such as induced low pressure, mesoscale leeward vortex, a terrain adduction line, thunderstorm, gravitational wave and the like can be generated when the typhoons are influenced by the terrains in the process of approaching the continents, so that the typhoon structure is changed. The change of typhoon structure often causes complex change of typhoon field, thus the nature and size of the anti-wind parameter are changed.
Therefore, the wind-resistant parameter characteristics of the engineering location are considered by the maximum value of the wind-resistant parameters of the selected engineering area in a point-to-point mode, so that the wind-resistant parameter characteristic difference among all parts in the engineering cannot be reflected comprehensively and accurately, the method is not suitable for actual engineering design, particularly under the influence of typhoon, and the accurate evaluation on the typhoon influence cannot be realized.
Some parametric or semi-empirical half-value wind farm models have also been developed, such as Batts (1980) wind farm models, Shapiro (1983) wind farm models, Yan Meng wind farm models (Meng et al 1995), Vickery (2000) models, CE (the U.S. army Corps of Engineers) wind farm models (1996), etc., to evaluate typhoon effects. However, the models do not consider complex atmospheric physical processes and underlying surface influences, the simulated typhoon intensity and the near-formation wind field have obvious difference from the actual condition, and the accuracy is poor.
Disclosure of Invention
In view of this, the present invention provides a method and an apparatus for evaluating typhoon-induced impact, so as to improve the accuracy of evaluating typhoon-induced impact. The technical scheme is as follows:
based on one aspect of the present invention, the present invention provides a method for evaluating an impact of typhoon, including:
obtaining lattice point wind field parameters output by simulating a typhoon wind field by using a mesoscale numerical mode;
for each lattice point meeting the preset requirement in the same horizontal plane: calculating destructive force parameters of the grid point, wherein the destructive force parameters comprise any one or more of wind speed space shear quantity, wind speed time shear quantity, wind direction angle space shear quantity and wind direction angle time shear quantity; the wind speed space tangent variable is the maximum value of the wind speed vector difference between the grid point and each grid point adjacent to the grid point, the wind direction angle space tangent variable is the maximum value of the wind direction angle vector difference between the grid point and each grid point adjacent to the grid point, the wind speed time tangent variable is the maximum value of the wind speed vector difference between the wind speed of the grid point at the moment t and the wind speed vector differences at the moment t-1 and the moment t +1, and the wind direction angle time tangent variable is the maximum value of the wind direction angle of the grid point at the moment t and the wind direction angle vector differences at the moment t-1 and the moment t + 1;
and evaluating the typhoon disaster-causing influence based on any one or more of the obtained wind speed space shear variable, wind speed time shear variable, wind direction angle space shear variable and wind direction angle time shear variable of each grid point and in combination with the wind speed value.
Optionally, the grid points meeting the preset requirement are grid points with 8 grid points adjacent to each other at the periphery.
Optionally, calculating the wind speed spatial shear amount of the grid point comprises:
using formulas
Calculating grid points
Amount of wind speed space shear
Calculating the wind direction angle space shear amount of the grid point comprises the following steps:
using formulas
Calculating grid points
Wind direction angle space tangent quantity theta (i, j)max；
Calculating the wind speed time shear amount of the grid point comprises the following steps:
using formulas
Calculating grid points
Time shear amount of wind speed
Calculating the wind direction angle time shear amount of the grid point comprises the following steps:
using formulas
Calculating grid points
Time-tangent variable of wind direction angle theta (i, j, t)max；
Wherein i and j are positive integers, and t represents time.
Optionally, the evaluating the typhoon-causing influence based on any one or more of the obtained wind speed space shear variable, wind speed time shear variable, wind direction angle space shear variable, and wind direction angle time shear variable of each grid point, and in combination with the wind speed value, includes:
and determining that the wind speed value is greater than a preset wind speed threshold value, and the disaster-causing influence of typhoon on grid points respectively meeting corresponding preset conditions is large in any one or more of the wind speed space tangent variable, the wind speed time tangent variable, the wind direction angle space tangent variable and the wind direction angle time tangent variable.
Optionally, the mesoscale numerical mode comprises a mesoscale weather forecast mode WRF numerical mode.
Based on another aspect of the present invention, the present invention provides an apparatus for evaluating an impact of typhoon, comprising:
the system comprises an acquisition unit, a processing unit and a control unit, wherein the acquisition unit is used for acquiring grid point wind field parameters output by simulating a typhoon wind field by using a mesoscale numerical mode;
the calculating unit is used for calculating the grid points meeting the preset requirement in the same horizontal plane: calculating destructive force parameters of the grid point, wherein the destructive force parameters comprise any one or more of wind speed space shear quantity, wind speed time shear quantity, wind direction angle space shear quantity and wind direction angle time shear quantity; the wind speed space tangent variable is the maximum value of the wind speed vector difference between the grid point and each grid point adjacent to the grid point, the wind direction angle space tangent variable is the maximum value of the wind direction angle vector difference between the grid point and each grid point adjacent to the grid point, the wind speed time tangent variable is the maximum value of the wind speed vector difference between the wind speed of the grid point at the moment t and the wind speed vector differences at the moment t-1 and the moment t +1, and the wind direction angle time tangent variable is the maximum value of the wind direction angle of the grid point at the moment t and the wind direction angle vector differences at the moment t-1 and the moment t + 1;
and the evaluation unit is used for evaluating the typhoon disaster-causing influence based on any one or more of the obtained wind speed space shear variable, wind speed time shear variable, wind direction angle space shear variable and wind direction angle time shear variable of each grid point and combining the wind speed value.
Optionally, the grid points meeting the preset requirement are grid points with 8 grid points adjacent to each other at the periphery.
Optionally, the computing unit is specifically configured to,
using formulas
Calculating grid points
Amount of wind speed space shear
Using formulas
Calculating grid points
Wind direction angle space tangent quantity theta (i, j)max；
Using formulas
Calculating grid points
Time shear amount of wind speed
Using formulas
Calculating grid points
Time-tangent variable of wind direction angle theta (i, j, t)max；
Wherein i and j are positive integers, and t represents time.
Optionally, the evaluation unit is specifically configured to determine that a wind speed value is greater than a preset wind speed threshold, and that any one or more of the wind speed space tangent variable, the wind speed time tangent variable, the wind direction angle space tangent variable, and the wind direction angle time tangent variable respectively satisfy a grid point of a corresponding preset condition, is subjected to a large typhoon causing effect.
Optionally, the mesoscale numerical mode comprises a mesoscale weather forecast mode WRF numerical mode.
According to the method and the device for evaluating the typhoon disaster-causing influence, the typhoon wind field simulated by the mesoscale numerical mode is utilized, the typhoon intensity is closer to the real typhoon intensity, the simulation precision of a typhoon path, the intensity and the near-ground wind field is ensured, and the evaluation accuracy of the typhoon influence is further ensured. The method can comprehensively and accurately reflect the parameter characteristic difference of the destructive force among all parts in the engineering influenced by typhoon and further ensure the evaluation accuracy of the typhoon induced disaster influence by calculating any one or more of the wind speed space tangent variable, the wind speed time tangent variable, the wind direction angle space shear quantity and the wind direction angle time tangent variable of each lattice point meeting the preset requirement and finally evaluating the typhoon induced disaster influence based on any one or more of the obtained wind speed space tangent variable, the wind speed time tangent variable, the wind direction angle space shear quantity and the wind direction angle time tangent variable of each lattice point and combining the wind speed value.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
Fig. 1 is a flowchart of a method for evaluating typhoon-induced disaster influence according to the present invention;
FIG. 2 shows grid points of the present invention
The wind speed space shear distribution schematic diagram of (1);
FIG. 3 shows grid points of the present invention
A wind speed time shear diagram at time t;
FIG. 4 is a schematic diagram of the time course of wind speed space shear at A, B, C in three typical locations during the typhoon "black lattice ratio";
FIG. 5 is a schematic diagram of the wind direction shear time course changes of A, B, C at three typical locations during the typhoon 'black lattice ratio';
fig. 6 is a schematic structural diagram of an evaluation apparatus for typhoon-induced influence according to the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, the method for evaluating the typhoon-caused influence provided by the present invention includes:
The mesoscale numerical mode, also called mesoscale mode, has gradually improved description on the physical process of atmospheric motion, the influence of the underlying surface, the physical process of the boundary layer and the like, and can better simulate mesoscale structural features in the typhoon. The typhoon wind field simulated by the mesoscale numerical mode has typhoon intensity which is closer to the real typhoon intensity, thereby ensuring the simulation precision of the typhoon path, intensity and the near-ground wind field and further ensuring the evaluation accuracy of typhoon influence.
In The present invention, The adopted mesoscale numerical mode may specifically be a mesoscale WRF (The Weather Research and Forecasting Model, Weather forecast mode) numerical mode. The middle-scale WRF numerical mode adopts a Lambert projection mode, the horizontal direction is an Arakawa-C jumping point grid, and the vertical direction adopts terrain to follow the static air pressure coordinate. The mesoscale WRF numerical mode adopts a 2-layer bidirectional nesting mode, the time step of an outer-layer grid is 36 seconds, the layer top is 50hPa, and the vertical direction is 31 layers. The number of element field grid points obtained by outer grid simulation is 813 multiplied by 548, the horizontal resolution is 3km, and the integral area is 2436km multiplied by 1641 km; the number of element field lattice points obtained by inner layer grid simulation is 1069 multiplied by 991, the horizontal resolution is 1km, and the integration area is 1068km multiplied by 990 km.
In the invention, the output time resolution of the typhoon wind field simulated by the mesoscale WRF numerical mode can be 10min, the lattice point distance can be 1km multiplied by 1km, and the acquired lattice point wind field parameters output by the typhoon wind field simulated by the mesoscale numerical mode can comprise the height of the terrain, the wind speed and the wind direction at the height of 10m following the terrain, and the like.
102, for each lattice point meeting the preset requirement in the same horizontal plane: and calculating destructive force parameters of the grid point, wherein the destructive force parameters comprise any one or more of wind speed space shear quantity, wind speed time shear quantity, wind direction angle space shear quantity and wind direction angle time shear quantity. The wind speed space tangent variable is the maximum value of the wind speed vector difference between the grid point and each grid point adjacent to the grid point, the wind direction angle space tangent variable is the maximum value of the wind direction angle vector difference between the grid point and each grid point adjacent to the grid point, the wind speed time tangent variable is the maximum value of the wind speed vector difference between the wind speed of the grid point at the moment t and the wind speed vector differences at the moment t-1 and the moment t +1, and the wind direction angle time tangent variable is the maximum value of the wind direction angle of the grid point at the moment t and the wind direction angle vector differences at the moment t-1 and the moment t + 1.
In order to ensure the accuracy of evaluating the typhoon-induced influence, the method preferably calculates the destructive force parameters of the grid point and simultaneously comprises four parameters of a wind speed space shear variable, a wind speed time shear variable, a wind direction angle space shear variable and a wind direction angle time shear variable.
In the invention, the grid points meeting the preset requirement in the same horizontal plane are grid points with 8 adjacent grid points around the grid points, namely the grid points comprise 8 adjacent grid points. Take the example shown in FIG. 2, in which the grid points are located at intermediate positions
Comprises 8 adjacent grid points, which are respectively the grid points at the upper positions thereof
Lattice points at the lower position thereof
Grid point at its left position
Grid point at its right position
Grid point at the upper left position
Lattice points at the upper right position thereof
Grid point at the lower left position
Grid point at the lower right position
For grid points with the number of the adjacent grid points less than 8, the grid points which do not meet the preset requirement are determined by the invention and are not processed.
For each lattice point meeting the preset requirement in the same horizontal plane, the method calculates the destructive power parameter of each lattice point meeting the preset requirement.
In the present invention, the destructive power parameter of a grid point refers to a "tangential variable" parameter of space and time of wind speed and wind direction of the grid point, which specifically includes: wind speed space shear variable, wind speed time shear variable, wind direction angle space shear variable and wind direction angle time shear variable.
The invention decomposes the wind damage factor (key factor of structure damage) of the wind to the structure into the wind direction, the rapid change of the wind speed in time (10min) and space (1 multiplied by 1Km), which are called as the time shear of the wind speed, the time shear of the wind direction, the space shear of the wind speed and the space shear of the wind direction.
To quantitatively characterize the size and distribution of this "shear", the present invention uses the lattice point wind field parameters output using the mesoscale numerical model, defined as follows:
the wind speed space tangent variable is the maximum value of the wind speed vector difference between the grid point and each grid point adjacent to the grid point; the wind direction angle space tangent variable is the maximum value of the wind direction angle vector difference between the grid point and each grid point adjacent to the grid point; the wind speed time switching variable is the maximum value of the wind speed vector difference of the grid point at the time t and the wind speed vector difference at the time t-1 and the time t + 1; the wind direction angle time-cut variable is the maximum value of the difference between the wind direction angle of the grid point at the time t and the wind direction angle vectors at the time t-1 and the time t + 1.
Specifically, taking the example shown in FIG. 2, for a grid point
The tangential variable of the wind speed space is a lattice point
8 grid points adjacent to its periphery
The maximum value of the wind speed vector difference between the grid points is shown in the following formula 1:
for lattice points
The wind direction angle space tangent is a lattice point
8 grid points adjacent to its periphery
The maximum value of the wind direction angle vector difference between the grid points is shown in the following formula 2:
wherein i and j are positive integers.
Further, as shown in conjunction with FIG. 3, for a grid point
The time-cut variable of the wind speed is a lattice point
The maximum value of the difference between the wind speed at time t and the wind speed vectors at time t-1 and time t +1 is shown in the following equation 3:
for lattice points
The time tangent of the wind direction angle is a lattice point
The maximum value of the difference between the wind direction angle at the time t and the wind direction angle vectors at the time t-1 and the time t +1 is shown in the following formula 4:
and 103, evaluating the typhoon disaster-causing influence based on any one or more of the obtained wind speed space shear variable, wind speed time shear variable, wind direction angle space shear variable and wind direction angle time shear variable of each grid point and combining the wind speed value.
Based on any one or more of the wind speed space tangent variable, the wind speed time tangent variable, the wind direction angle space tangent variable and the wind direction angle time tangent variable of each grid point in the plurality of grid points obtained in the step 102 of the present invention, the evaluation of the typhoon induced disaster effect is realized by judging any one or more of the wind speed space tangent variable, the wind speed time tangent variable, the wind direction angle space tangent variable and the wind direction angle time tangent variable of each grid point and combining the currently detected wind speed values at each grid point.
In practical application, the wind speed value indicates the wind speed at the grid point, and the larger the value is, the larger the wind speed is, the stronger the typhoon-causing force is. The wind speed space shear variable, the wind speed time shear variable, the wind direction angle space shear variable and the wind direction angle time shear variable are quantities indicating the magnitude of the torque applied to the engineering, the building and the like, and the larger the values of the quantities, the more uneven the force applied to the engineering and the building is, the larger the risk of the generated rotation effect and the accompanying torsional deformation is, namely, the larger the torque is, and the stronger the disaster-causing force of the typhoon is.
Therefore, in the present invention, for a grid point with a wind speed value greater than a preset wind speed threshold, if any one or more of the destructive force parameters of the grid point, that is, any one or more of the wind speed space shear variable, the wind speed time shear variable, the wind direction angle space shear variable, and the wind direction angle time shear variable, is greater than the corresponding preset threshold, it may be determined that the grid point is greatly affected by the typhoon.
It should be noted that, in the present invention, the size of the preset wind speed threshold value can be flexibly set according to actual needs, which is not limited by the present invention; the preset conditions respectively corresponding to the wind speed space shear variable, the wind speed time shear variable, the wind direction angle space shear quantity and the wind direction angle time shear quantity are as follows: the value of the wind speed space cut variable is larger than a first threshold value, the value of the wind speed time cut variable is larger than a second threshold value, the value of the wind direction angle space cut variable is larger than a third threshold value, and the value of the wind direction angle time cut variable is larger than a fourth threshold value.
In order to demonstrate the feasibility of the evaluation method for typhoon-induced influence, the invention uses a mesoscale WRF numerical mode to simulate the typhoon wind condition during the landing period of No. 0814 strong typhoon 'blackroom ratio', and takes two parameters of wind speed space shear quantity and wind direction angle space shear quantity as the consideration factors of the typhoon-induced influence.
Firstly, the WRF numerical mode is used for carrying out fine and accurate numerical simulation on the 'blackroom ratio' of No. 0814 strong typhoon, and the typhoon path, strength and wind field simulated by the mesoscale WRF numerical mode are ensured to be closer to the real typhoon path, strength and wind field. The destructive force parameters are calculated by utilizing the typhoon wind speed and wind direction obtained by simulation, and the result shows that the wind speed space shear of the typhoon is generally distributed in a spiral strip shape around the center of the typhoon, and the large value area is mainly concentrated in two areas: the area near the typhoon eye is the left rear area deviating from the advancing direction of the typhoon; and secondly, a mountain region with large topographic relief. And the method is very close to the real typhoon disaster influence, and the accuracy of the typhoon disaster influence evaluation is ensured.
On the field top of the cloud and fog mountains near the landing path of the 'black lattice ratio', a plurality of lattice points are selected on the windward side (marked as A), the mountain top (marked as B) and the leeward side (marked as C) of the mountains respectively, and the shear parameters of the lattice points are calculated every 10 minutes during the wind speed influence period (23 days 18 hours-24 days 00 hours) of the typhoon process over level 8, as shown in figures 4 and 5, the wind field situation of A, B, C at a certain moment relative to the typhoon in one test example is respectively shown, and the altitude of 3 typical test points is 36.4m, 1296.5m and 47.2m respectively. In combination with the table 1 below.
TABLE 1 mean and standard deviation of spatial shear at each lattice point during "Black lattice ratio" of typhoon
It can be seen that the spatial shear of wind speed is the largest at point B on the top of the hill, followed by point C on the leeward side of the hill, and the spatial shear of wind speed is the smallest at point a on the windward side. From the standard deviation, the maximum standard deviation of the wind speed space tangential variable of the point C in the typhoon process is 0.55, and the minimum standard deviation of the point A is 0.23; the maximum value of the wind direction angle space tangent occurs at point C at 7.42, and other characteristics of the wind direction angle space tangent are similar to the wind speed space shear. Therefore, the shear value and the variation of the windward side point A are minimum no matter the wind speed space shear quantity or the wind direction angle space shear quantity, various parameters are increased even if a typhoon eye wall strong wind area is close, the mountain leeward side point C is most prominent in performance, when the typhoon eye area passes through and the wind direction is changed into south wind, the mountain leeward side is changed into the mountain windward side, the roughness of the windward underlying cushion surface is obviously reduced, and the wind speed and the wind direction space shear are obviously reduced.
It can be understood that the evaluation result of the 'blackroom ratio' of the strong typhoon number 0814 by adopting the evaluation method of the typhoon influence provided by the invention is very close to the actual typhoon disaster-causing influence of the 'blackroom ratio' of the strong typhoon number 0814.
It should be noted that the above verification process is only exemplarily described by taking two parameters, namely, the wind speed space shear amount and the wind direction angle space shear amount, as examples, and in the practical application process, it is preferable that the four parameters, namely, the wind speed space shear amount, the wind speed time shear amount, the wind direction angle space shear amount and the wind direction angle time shear amount, are taken as consideration factors of the typhoon-induced disaster influence at the same time.
Based on the method for evaluating typhoon-causing influence provided by the invention in the previous text, the invention further provides a device for evaluating typhoon-causing influence, as shown in fig. 6, comprising:
the acquiring unit 100 is configured to acquire a lattice point wind field parameter output by simulating a typhoon wind field in a mesoscale numerical mode;
a calculating unit 200, configured to, for each grid point meeting the preset requirement in the same horizontal plane: calculating destructive force parameters of the grid point, wherein the destructive force parameters comprise any one or more of wind speed space shear quantity, wind speed time shear quantity, wind direction angle space shear quantity and wind direction angle time shear quantity; the wind speed space tangent variable is the maximum value of the wind speed vector difference between the grid point and each grid point adjacent to the grid point, the wind direction angle space tangent variable is the maximum value of the wind direction angle vector difference between the grid point and each grid point adjacent to the grid point, the wind speed time tangent variable is the maximum value of the wind speed vector difference between the wind speed of the grid point at the moment t and the wind speed vector differences at the moment t-1 and the moment t +1, and the wind direction angle time tangent variable is the maximum value of the wind direction angle of the grid point at the moment t and the wind direction angle vector differences at the moment t-1 and the moment t + 1;
and the evaluation unit 300 is configured to evaluate the typhoon-caused disaster influence based on any one or more of the obtained wind speed space shear quantity, wind speed time shear quantity, wind direction angle space shear quantity, and wind direction angle time shear quantity of each grid point, and in combination with the wind speed value.
Optionally, the grid points meeting the preset requirement are grid points with 8 grid points adjacent to each other at the periphery. The mesoscale numerical mode may be a mesoscale weather forecast mode WRF numerical mode.
Alternatively, the computing unit 200 in the present invention is specifically adapted to,
using formulas
Calculating grid points
Amount of wind speed space shear
Using formulas
Calculating grid points
Wind direction angle space tangent quantity theta (i, j)max；
Using formulas
Calculating grid points
Time shear amount of wind speed
Using formulas
Calculating grid points
Time-tangent variable of wind direction angle theta (i, j, t)max；
Wherein i and j are positive integers, and t represents time.
The evaluation unit 300 in the present invention is specifically configured to determine that a wind speed value is greater than a preset wind speed threshold, and that any one or more of the wind speed space tangent, the wind speed time tangent, the wind direction angle space tangent, and the wind direction angle time tangent have a large disaster-causing influence on typhoons at grid points that respectively satisfy corresponding preset conditions.
It should be noted that, in the present specification, the embodiments are all described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments may be referred to each other.
Finally, it should also be noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and decorations can be made without departing from the principle of the present invention, and these modifications and decorations should also be regarded as the protection scope of the present invention.
Claims (10)
1. A method for evaluating the influence of typhoon disaster is characterized by comprising the following steps:
obtaining lattice point wind field parameters output by simulating a typhoon wind field by using a mesoscale numerical mode;
for each lattice point meeting the preset requirement in the same horizontal plane: calculating destructive force parameters of the grid point, wherein the destructive force parameters comprise any one or more of wind speed space shear quantity, wind speed time shear quantity, wind direction angle space shear quantity and wind direction angle time shear quantity; the wind speed space tangent variable is the maximum value of the wind speed vector difference between the grid point and each grid point adjacent to the grid point, the wind direction angle space tangent variable is the maximum value of the wind direction angle vector difference between the grid point and each grid point adjacent to the grid point, the wind speed time tangent variable is the maximum value of the wind speed vector difference between the wind speed of the grid point at the moment t and the wind speed vector differences at the moment t-1 and the moment t +1, and the wind direction angle time tangent variable is the maximum value of the wind direction angle of the grid point at the moment t and the wind direction angle vector differences at the moment t-1 and the moment t + 1;
and evaluating the typhoon disaster-causing influence based on any one or more of the obtained wind speed space shear variable, wind speed time shear variable, wind direction angle space shear variable and wind direction angle time shear variable of each grid point and in combination with the wind speed value.
2. The method of claim 1, wherein the grid points satisfying the predetermined requirement are 8 grid points with a number of grid points adjacent to each other.
3. The method of claim 2,
calculating the wind speed space shear amount of the grid point comprises the following steps:
using formulas
Calculating grid points
Amount of wind speed space shear
Calculating the wind direction angle space shear amount of the grid point comprises the following steps:
using formulas
Calculating grid points
Wind direction angle space tangent quantity theta (i, j)max；
Calculating the wind speed time shear amount of the grid point comprises the following steps:
using formulas
Calculating grid points
Time shear amount of wind speed
Calculating the wind direction angle time shear amount of the grid point comprises the following steps:
using formulas
Calculating grid points
Time-tangent variable of wind direction angle theta (i, j, t)max；
Wherein i and j are positive integers, and t represents time.
4. The method according to any one of claims 1 to 3, wherein the evaluating the typhoon-causing influence based on any one or more of the obtained wind speed space shear, wind speed time shear, wind direction angle space shear and wind direction angle time shear of each grid point and in combination with the wind speed value comprises:
and determining that the wind speed value is greater than a preset wind speed threshold value, and the disaster-causing influence of typhoon on grid points respectively meeting corresponding preset conditions is large in any one or more of the wind speed space tangent variable, the wind speed time tangent variable, the wind direction angle space tangent variable and the wind direction angle time tangent variable.
5. The method of any of claims 1-3, wherein the mesoscale numerical mode comprises a mesoscale weather forecast mode, WRF, numerical mode.
6. An evaluation device for typhoon-induced influence, comprising:
the system comprises an acquisition unit, a processing unit and a control unit, wherein the acquisition unit is used for acquiring grid point wind field parameters output by simulating a typhoon wind field by using a mesoscale numerical mode;
the calculating unit is used for calculating the grid points meeting the preset requirement in the same horizontal plane: calculating destructive force parameters of the grid point, wherein the destructive force parameters comprise any one or more of wind speed space shear quantity, wind speed time shear quantity, wind direction angle space shear quantity and wind direction angle time shear quantity; the wind speed space tangent variable is the maximum value of the wind speed vector difference between the grid point and each grid point adjacent to the grid point, the wind direction angle space tangent variable is the maximum value of the wind direction angle vector difference between the grid point and each grid point adjacent to the grid point, the wind speed time tangent variable is the maximum value of the wind speed vector difference between the wind speed of the grid point at the moment t and the wind speed vector differences at the moment t-1 and the moment t +1, and the wind direction angle time tangent variable is the maximum value of the wind direction angle of the grid point at the moment t and the wind direction angle vector differences at the moment t-1 and the moment t + 1;
and the evaluation unit is used for evaluating the typhoon disaster-causing influence based on any one or more of the obtained wind speed space shear variable, wind speed time shear variable, wind direction angle space shear variable and wind direction angle time shear variable of each grid point and combining the wind speed value.
7. The apparatus of claim 6, wherein the grid points satisfying the predetermined requirement are 8 grid points with a number of grid points adjacent to each other.
8. The apparatus according to claim 7, characterized in that the calculation unit is specifically configured to,
using formulas
Calculating grid points
Amount of wind speed space shear
Using formulas
Calculating grid points
Wind direction angle space tangent quantity theta (i, j)max；
Using formulas
Calculating grid points
Time shear amount of wind speed
Using formulas
Calculating grid points
Time-tangent variable of wind direction angle theta (i, j, t)max；
Wherein i and j are positive integers, and t represents time.
9. The apparatus according to any one of claims 6 to 8, wherein the evaluation unit is specifically configured to determine that a wind speed value is greater than a preset wind speed threshold, and that any one or more of the wind speed spatial variation, the wind speed time variation, the wind direction angle spatial variation and the wind direction angle time variation are greatly affected by typhoon at a grid point that respectively satisfies a corresponding preset condition.
10. The apparatus of any of claims 6-8, wherein the mesoscale numerical mode comprises a mesoscale weather forecast mode, WRF, numerical mode.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810239480.5A CN108399507B (en) | 2018-03-22 | 2018-03-22 | Typhoon disaster influence assessment method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810239480.5A CN108399507B (en) | 2018-03-22 | 2018-03-22 | Typhoon disaster influence assessment method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108399507A CN108399507A (en) | 2018-08-14 |
| CN108399507B true CN108399507B (en) | 2022-02-25 |
Family
ID=63091560
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810239480.5A Active CN108399507B (en) | 2018-03-22 | 2018-03-22 | Typhoon disaster influence assessment method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108399507B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113177505B (en) * | 2021-05-14 | 2022-03-04 | 中国科学院大气物理研究所 | Method and device for identifying vortex in meteorological field, storage medium and electronic equipment |
| CN114647018A (en) * | 2022-02-23 | 2022-06-21 | 中国人民解放军国防科技大学 | Squall line strength evaluation method based on cold pool |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103778572A (en) * | 2014-02-24 | 2014-05-07 | 南方电网科学研究院有限责任公司 | Sea wind resource evaluation method based on WRF mode |
| CN104715297A (en) * | 2015-04-10 | 2015-06-17 | 北京安澜尔雅科技有限公司 | Method for predicting wind-induced roof covering loss |
| CN107229825A (en) * | 2017-05-23 | 2017-10-03 | 浙江大学 | A kind of tropical cyclone complete trails analogy method assessed towards calamity source |
| US9816863B1 (en) * | 2014-09-26 | 2017-11-14 | Harris Corporation | Wind field vertical profile estimation using spectral radiance of multiband imagery and temperature profiles |
| CN107526852A (en) * | 2016-06-21 | 2017-12-29 | 中国辐射防护研究院 | A kind of nuclear facilities Accident Off-site Consequence real-time online evaluation method and system |
-
2018
- 2018-03-22 CN CN201810239480.5A patent/CN108399507B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103778572A (en) * | 2014-02-24 | 2014-05-07 | 南方电网科学研究院有限责任公司 | Sea wind resource evaluation method based on WRF mode |
| US9816863B1 (en) * | 2014-09-26 | 2017-11-14 | Harris Corporation | Wind field vertical profile estimation using spectral radiance of multiband imagery and temperature profiles |
| CN104715297A (en) * | 2015-04-10 | 2015-06-17 | 北京安澜尔雅科技有限公司 | Method for predicting wind-induced roof covering loss |
| CN107526852A (en) * | 2016-06-21 | 2017-12-29 | 中国辐射防护研究院 | A kind of nuclear facilities Accident Off-site Consequence real-time online evaluation method and system |
| CN107229825A (en) * | 2017-05-23 | 2017-10-03 | 浙江大学 | A kind of tropical cyclone complete trails analogy method assessed towards calamity source |
Non-Patent Citations (3)
| Title |
|---|
| Evaluation of hurricane ocean vector winds from WindSat;I.S. Adams等;《IEEE Transactions on Geoscience and Remote Sensing》;20060221;第44卷(第3期);656-667 * |
| 不同海岸地形下台风边界层风场精细模拟及其风工程参数特性;薛霖;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20160315(第3期);C038-697 * |
| 浙江海岸台风风暴潮漫堤风险评估研究;卢美;《中国博士学位论文全文数据库 基础科学辑》;20140615(第6期);A010-3 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108399507A (en) | 2018-08-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Saeed et al. | Comparison of six different methods of Weibull distribution for wind power assessment: A case study for a site in the Northern region of Pakistan | |
| EP3761213A1 (en) | A wind field dynamic downscaling method based on aerodynamic parameters of simplified terrain | |
| Beck et al. | Reconstruction of near-surface tornado wind fields from forest damage | |
| Đurišić et al. | A model for vertical wind speed data extrapolation for improving wind resource assessment using WAsP | |
| Kalmikov et al. | Wind power resource assessment in complex urban environments: MIT campus case-study using CFD Analysis | |
| Sheridan | Review of techniques and research for gust forecasting and parameterisation | |
| Tse et al. | A comparative study of typhoon wind profiles derived from field measurements, meso-scale numerical simulations, and wind tunnel physical modeling | |
| Miller et al. | Topographic speed-up effects and observed roof damage on Bermuda following Hurricane Fabian (2003) | |
| da Costa et al. | Computer simulation of atmospheric flows over real forests for wind energy resource evaluation | |
| Larsén et al. | Extreme winds from the NCEP/NCAR reanalysis data | |
| Tabrizi et al. | Extent to which international wind turbine design standard, IEC61400-2 is valid for a rooftop wind installation | |
| CN108399507B (en) | Typhoon disaster influence assessment method and device | |
| Santos et al. | The Alaiz experiment: untangling multi-scale stratified flows over complex terrain | |
| JP4297921B2 (en) | Method for creating contour map of physical quantity and method for estimating physical quantity | |
| Ibrahim et al. | Wind characteristics and gis-based spatial wind mapping study in Malaysia | |
| Konow | Tall wind profiles in heteorogeneous terrain | |
| CN108363882A (en) | A kind of mountain area Transmission Line Design wind speed projectional technique based on power NO emissions reduction pattern | |
| CN110824585B (en) | Method for measuring gale wind-rising mechanism of complex terrain area | |
| KR101800789B1 (en) | Improved Modified Korean Parameter-elevation Regressions on Independent Slopes Model for applying interpolation of wind | |
| Çakmak et al. | Comparative Performance Analysis of Different Wind Fields in the Southern and North-western Coastal Areas of the Black Sea | |
| JP4043779B2 (en) | Basic wind speed map creation method by wind direction and basic wind speed estimation method by wind direction | |
| Tang et al. | Impact of horizontal resolution in CALMET on simulated near-surface wind fields over complex terrain during Super Typhoon Meranti (2016) | |
| Nawri et al. | Evaluation of WRF mesoscale model simulations of surface wind over Iceland | |
| Kim et al. | Long-term wind resource mapping of Korean west-south offshore for the 2.5 GW offshore wind power project | |
| Nemes | Statistical analysis of wind speed profile: a case study from Iasi Region, Romania |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |