CN107526910A - A kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation - Google Patents
A kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation Download PDFInfo
- Publication number
- CN107526910A CN107526910A CN201610453592.1A CN201610453592A CN107526910A CN 107526910 A CN107526910 A CN 107526910A CN 201610453592 A CN201610453592 A CN 201610453592A CN 107526910 A CN107526910 A CN 107526910A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- component
- wind field
- mfrac
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A kind of wind field diagnostic method in being evaluated the invention discloses nuclear facilities Accident Off-site Consequence, this method include：The simulation region scope of nuclear facilities accident is determined centered on nuclear facilities, obtains meteorological data and geographical data in the range of simulation region, and three-dimensional initial wind field is set up according to the meteorological data and geographical data got；The first successive step of three-dimensional initial wind field is carried out according to numerical weather forecast data and geographical data afterwards, data are adjusted again to the wind field after first successive step according to making weather observations, and obtain final dimensional wind.Can be on the basis of the meteorological data and environmental information data got using this method, the local meteorological field of complicated local orographic condition can more finely be reflected by constructing, and the impact evaluation for ambient air radioactive substance after the release of follow-up nuclear facilities accident provides meteorological field data.
Description
Technical field
The present invention relates to nuclear facilities accident consequence analysis technical field, and in particular to a kind of nuclear facilities Accident Off-site Consequence is commented
Data preprocessing method in valency.
Background technology
Nuclear facilities or other nuclear facilities substantial deviation accidental conditions, accident occurs when, the release of radioactive substance
Or due control may be lost, reached unacceptable level, it will cause serious consequence.Therefore, to nuclear facilities
The strick precaution control and analysis of accident harm, have important strategic importance.
Nuclear facilities damage sequence evaluation system be by by region meteorological data and radioactive pollutant migration model and after
The evaluation system that fruit evaluation becomes one, for the damage sequence estimated, evaluated and be shown in particular range, the system is
Nuclear facilities or other nuclear facilities Emergency Preparedness and the important component of emergency protective action decision-making.Nuclear facilities damage sequence is evaluated
System, it is necessary first to the data such as the classification source item different with content, meteorology, environmental information are pre-processed, can be to above-mentioned
Each data realize that effective processing has vital influence on the analysis result of follow-up damage sequence, and the present invention is exactly to be directed to
The problem and propose the wind field diagnostic method in a kind of evaluation of nuclear facilities Accident Off-site Consequence.
The content of the invention
For defect present in prior art, it is an object of the invention to provide one kind can construct local meteorological field
Nuclear facilities Accident Off-site Consequence evaluation in wind field diagnostic method.
To achieve the above object, the technical solution adopted by the present invention is as follows：
A kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation, comprises the following steps：
Step 1: centered on nuclear facilities, the simulation region scope of nuclear facilities accident is determined, obtains the meteorological data of simulation region
And geographical data, and establish according to the meteorological data and geographic position data of acquisition the three-dimensional initial wind field of simulation region；Institute
Stating meteorological data includes make weather observations data and the numerical weather forecast data of meteorological department of weather station of simulation region；Institute
Stating geographical data includes latitude and longitude information, underlying surface elevation, vegetation, water body distribution characteristics and the land use of simulation region；
Three-dimensional initial wind field is multi-layer net wind field；
Step 2: according to the numerical weather forecast data of the meteorological department and geographical data to described three-dimensional
Initial wind field carries out just successive step, obtains the wind field after first successive step；The just successive step is included according to landform kinematics effect
Three-dimensional initial wind field three-dimensional initial wind field is adjusted, three-dimensional initial wind field three-dimensional initial wind field carried out according to slope stream effect
Adjust and the blocking effect of wind field is adjusted to three-dimensional initial wind field three-dimensional initial wind field according to landform thermokinetics；
Step 3: being handled according to the data of making weather observations the wind field after first successive step, final wind field is obtained；Institute
State the adjustment that processing includes carrying out the wind field after first successive step successively interpolation processing, smoothing processing and vertical velocity component.
Further, the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation as described above, in step 2,
It is described three-dimensional initial wind field three-dimensional initial wind field is adjusted according to landform kinematics effect including：Imitated according to landform kinematics
Reply vertical direction wind speed component is adjusted and horizontal direction wind speed component is adjusted according to landform kinematics effect；
It is according to the mode that landform kinematics effect is adjusted to vertical direction wind speed component：
2.1) the vertical component w of wind speed in cartesian coordinate system is calculated, calculation formula is：
Wherein, V is simulation region zone leveling wind speed, htFor Terrain Elevation,For htGradient, represent difference in height everywhere,
K is the exponential damping coefficient relevant with atmospheric stability, and z is the vertical coordinate under cartesian coordinate system, and N is Brunt-Visa
Frequency is drawn, | V | it is the absolute value of zone leveling wind speed, g is acceleration of gravity, and θ is environment position temperature；
2.2) the vertical component W of wind speed under landform tracking coordinate system is calculated, calculation formula is：
Wherein, u, v be respectively under cartesian coordinate system wind speed in the component of x-axis horizontal direction and in y-axis horizontal direction
Component,For gradient of the Terrain Elevation in x-axis direction,For gradient of the Terrain Elevation in y-axis direction；X-axis horizontal direction is
East-west, y-axis horizontal direction is south-north direction；
It is according to the mode that landform kinematics effect is modified to horizontal direction wind field component：Using deviation minimum side
Method is adjusted to Horizontal Winds component using iterative method, until the three-dimensional deviation after adjustment is less than setting value.
Further, the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation as described above, in step 2,
It is according to the mode that slope stream effect is adjusted to three-dimensional initial wind field：
1) slope stream wind speed caused by the stream effect of slope is calculated；Slope stream wind speed S calculation formula is：
S={ [Qhgx sinα/[(ρcpT)(CD+k)]}1/3[1-exp(-x/Le)]1/3
Le=h/ (CD+k)
Wherein, QhFor sensible heat flux, g is acceleration of gravity, and x flows the distance apart from top of the slope for slope, and α is the slope flowing water straight angle
Degree, ρ are atmospheric density, cpFor heat capacity of air, T is air themperature, CDFor ground drag coefficient, k is the entrainment system at the top of the fluid layer of slope
Number, LeFor equilibrium length yardstick, h is slope stream height；
2) wind speed is flowed according to slope and slope flow horizontal angle [alpha] calculates the slope stream wind speed u of x-axis horizontal directionsWith y-axis level side
To slope stream wind speed vs；
3) slope stream wind speed is added to three-dimensional initial wind field, the horizontal wind speed of three-dimensional initial wind field is adjusted, adjusted
Formula be：
u1'=u1+us
v1'=v1+vs
Wherein, u1' and v1' it is respectively that the x-axis horizontal direction after effect is modified to three-dimensional initial wind field is flowed according to slope
The wind speed component of wind speed component and y-axis horizontal direction, u1And v1The wind speed of x-axis horizontal direction point in respectively three-dimensional initial wind field
The wind speed component of amount and y-axis horizontal direction.
Further, the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation as described above, in step 1),
When slope stream is descending stream,
CD=K=4 × 10-2
H=0.05 Δs Z
Sin α=minimum (sin α, Δ Z/x)
When slope stream is mountain upslope flow, (CD+ k)~1, x sin α=Δ Z.
Wherein, Δ Z is slope stream height.
Further, the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation as described above, in step 2,
It is described the blocking effect of wind field is adjusted to three-dimensional initial wind field according to landform thermokinetics including：
The local Froude number Fr of each mesh point in three-dimensional initial wind field is calculated, calculation formula is：
Δht=(hmax)ij-(z)ijk
Wherein, V is the wind speed of mesh point, and N is that frequency, Δ h are drawn in Brunt-VisatFor effective barrier height, (hmax)ijFor
Maximum Terrain Elevation in mesh point (i, j) radius of influence, (z)ijkFor the height of mesh point (i, j) in high dead level k；
Judge whether the local Froude number Fr of mesh point is less than the critical Froude number of setting, if it is not, then not at the beginning of three-dimensional
Beginning wind field is adjusted, if so, whether have upward component, if so, then by the wind direction of mesh point if then judging the wind field of mesh point
The direction tangent with landform is adjusted to, if it is not, not being adjusted to three-dimensional initial wind field then.
Further, the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation as described above, in step 3,
According to making weather observations, data are to the interpolation formula of the wind field progress interpolation processing after first successive step：
Wherein, (u, v) '2For the initial wind speed component of the horizontal direction of mesh point in the wind field after interpolation processing, (u, v)1
For the wind speed component of the horizontal direction of mesh point in the wind field after the just successive step, k is the identification number of weather station,
(uobs,vobs)kFor the observation wind speed component of the horizontal direction of k-th of weather station, R is the wind field after default just successive step
Weight number；RkFor the distance of k-th of weather station to mesh point；
According to making weather observations, data are smoothed to wind field after interpolation processing, and smoothing processing formula is：
(ui,j)2"=0.5ui,j+0.125(ui-1,j+ui+1,j+ui,j-1+ui,j+1)
(vi,j)2"=0.5vi,j+0.125(vi-1,j+vi+1,j+vi,j-1+vi,j+1)
Its, (ui,j)2″、(vi,j)2" be respectively smoothing processing after in the horizontal direction at mesh point (i, j) place x-axis direction
Wind speed component and the wind speed component in y-axis direction；ui,j、vi,jAt mesh point (i, j) place after interpolation processing respectively before smoothing processing
Horizontal direction in the wind speed component in x-axis direction and the wind speed component in y-axis direction；
According to making weather observations, data include to the mode of the wind field progress vertical velocity adjustment after smoothing processing：According to not pressing
Contracting mass-conservation equation calculates vertical velocity component from the Horizontal Winds component after smoothing processing, the vertical velocity calculated
Component is final vertical velocity component, and the Horizontal Winds after smoothing processing are final Horizontal Winds, or,
Vertical velocity point is being calculated from the Horizontal Winds component after smoothing processing according to not compression quality conservation equation
After amount, the vertical velocity calculated is adjusted, the vertical velocity component after adjustment is final vertical velocity component.
Further, the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation as described above, according to not pressing
Contracting mass-conservation equation calculates vertical velocity component W from the Horizontal Winds component after smoothing processing1(z) formula is：
Wherein, u ", v " is respectively the wind speed component in x-axis direction and the wind in y-axis direction in the horizontal direction after smoothing processing
Fast component.
Further, the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation as described above, according to not
After compression quality conservation equation calculates vertical velocity component from the Horizontal Winds component after smoothing processing, hung down to what is calculated
Straight wind speed is adjusted, the vertical velocity component W after being adjusted2(z) calculation formula is：
W2(z)=W1(z)-(z/ztop)W1(z=ztop)
Wherein, z be cartesian coordinate system under vertical coordinate, ztopFor the overhead height of simulation region, W1(z=ztop) it is root
The vertical speed at the top of simulation region calculated according to not compression quality conservation equation from the Horizontal Winds component after smoothing processing
Component；
Calculate the vertical velocity component W after adjustment2(z) after, two horizontal directions are adjusted using minimum deviation method for uncertain
Wind speed component, minimum deviation method for uncertain keep vertical component be W2(z) in the case of constant, using alternative manner by two
The deviation adjusting of the wind speed component of horizontal direction is to the deviation threshold less than setting, two water after being adjusted using deviation method
Square to horizontal component be final Horizontal Winds component.
Further, the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation as described above, calculates tune
Vertical velocity component W after whole2(z) after, the adjustment of the wind speed component of two horizontal directions is adjusted using minimum deviation method for uncertain
Mode be：
Calculate the wind speed component deviation D of the horizontal wind direction at each mesh point (i, j, k) place in grid cellijk, calculate public
Formula is：
Wherein, zk+1/2-zk-1/2The floor height in vertical direction is represented, Δ x and Δ y are illustrated respectively in x-axis and y-axis direction
The length of side of grid cell；
The horizontal wind speed component of surrounding's grid cell of mesh point (i, j, k) is adjusted, makes the deviation D of mesh pointijkIt is less than
The deviation threshold ε of setting, adjustment formula are：
(unew)i+1,j,k=ui+1,j,k+uadj
(unew)i-1,j,k=ui-1,j,k-uadj
(vnew)i+1,j,k=vi,j+1,k+vadj
(vnew)i,j-1,k=vi,j-1,k-vadj
Wherein, uadjAnd vadjThe adjustment amount of wind speed component respectively in x-axis and y-axis direction；
Iteration uses above-mentioned adjustment mode, until the deviation of all grid cells is respectively less than the deviation threshold ε set.
Further, the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation as described above, to level
Before the wind field in direction carries out space interpolation, the data of each weather station of earth's surface can be extrapolated to higher level, including：
Power-law equation extrapolation, theory of similarity extrapolation, the custom dimensions factor are extrapolated.
The beneficial effects of the present invention are：Method provided by the present invention, obtaining weather forecast initial fields or environment gas
On the basis of observation data, heat exchange, the effect of mechanical flow-disturbing are considered, construction is estimated finer local meteorological field, set for core
The estimation of ambient air radioactive substance concentration provides meteorological field data after applying nuclear accident release.
Brief description of the drawings
Fig. 1 is a kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation in the specific embodiment of the invention
Flow chart.
Embodiment
With reference to Figure of description, the present invention is described in further detail with embodiment.
In present embodiment, the Diagnosis Field pattern established using method provided by the present invention is referred to as DIAMET meteorologies
Pattern.
Fig. 1 shows the wind field in a kind of nuclear facilities Accident Off-site Consequence evaluation provided in the specific embodiment of the invention
The flow chart of diagnostic method, in Fig. 1 it can be seen from this method mainly include following four step：
Step S100：Centered on nuclear facilities, the simulation region scope of nuclear facilities accident is determined, is obtained in the range of simulation region
Meteorological data and geographical data；
Certain uncertainty be present in the acquisition degree and mode of the potentially useful meteorological data of nuclear facilities.Present embodiment
In, the meteorological data mainly includes make weather observations data and the meteorological department's number of the fixation weather station in the range of simulation region
It is worth data of weather forecast.It is meteorological that the data of making weather observations include the more key elements of wind, temperature, radiation, precipitation etc. that observation station observes
Data.
In actual applications, the meteorological data got can be made a distinction according to different demands, for example, according to
The scope for obtaining meteorological data is distinguished, and can be divided into mesoscale scope (centered on nuclear facilities radius 101~102Within scope)
With small range scale (centered on nuclear facilities radius 101Within) meteorological data；When nuclear facilities accident occurs, set according to core
Agent thus occur the moment, due to being influenceed by regional condition etc., from meteorological data obtain timing separation, can be divided into accident occur before
The data at moment, the data at some moment after accident occurs.
Obtaining for the numerical weather forecast data of meteorological part can be directly from the operation of meteorological numerical forecast department such as state
The Meteorology Forecast System extraction meteorological field such as T639, GFS of the Global Scale of family's level and the WRF/MM5 of regional scale of local level,
After the conversion of necessary form, the numerical weather forecast data of respective regions are obtained.
Because nuclear facilities wind-field model is a local model, horizontal scale is 102Km or so, it is desirable to provide as high as possible
Horizontal resolution meteorological field, by taking WRF patterns as an example, the assimilation of data is observed using three-dimensional variation (3DVAR) method
Calculate, the horizontal resolution that forecast system is finely assimilated in region can bring up to 15km.The layering of system in vertical direction is 35
Layer, top layer air pressure is 50hPa.System uses Euler's mass coordinate, using Runge-Kutta three rank Time integral schemes.System
Physical parameter scheme is respectively used by system experiment：RRTm long-wave radiations scheme, Dudhia shortwave radiations scheme, MYJ are near
Ground layered scheme, land surface emissivity are from heat diffusion scheme, MYJ PBL schemes, Betts-Miller-Janjic cumulus parameters
Change scheme, microphysical processes scheme are Lin schemes.Time effect forecast extended to 72 hours from current 48 hours, during data output
Between interval from business once bring up within 3 hours 30 minutes once.Twice a day forecast (universal time：When 00 and when 12).
Geographical data is indivisible important component in nuclear facilities incident response system, is mainly set according to core
Apply the situation of surrounding, obtain latitude and longitude information in the range of accident simulation region, underlying surface elevation, vegetation, water body, population, food
The characteristics such as chain, user's vocational distribution, land use, it is contemplated that the convenience of whole system design management, should additionally wrap
Include the release of source item and the conversion of respective regions demographic data.
Step S200：The three-dimensional initial wind of simulation region scope is established according to the meteorological data of acquisition and geographic position data
；
In present embodiment, the three-dimensional initial wind field is multi-layer net wind field, DIAMET Meteorological Models using one by
X, U, V, W expression is respectively adopted in the multi-layer net system that tri- directions of Y, Z are formed, the corresponding wind speed component on three directions, should
X and Y-axis are set in system and represents east-west trend and south-north trend respectively, Z represents vertical direction, this just with it is usually used
On corresponding horizontal plane wind field component --- u components (X-axis) and v components (Y-axis) are consistent.Mercator projection grid
(UTM) and generally use a kind of grid system.
After multi-layer net system is established, wind field value can be adjusted by input variable, can be with Lambert projection
Grid is consistent, so as to offset the influence brought due to earth curvature.DIAMET Meteorological Models use the normal latitude that user specifies
With calculate conic constants with reference to longitude, and according to reference to longitude come calculating observation station in the upward distance in east-west.Herein
On the basis of, the value of observation wind field and prediction wind field is adjusted further according to these values, can be consistent with Lambert mapping.Gas
Used as default standard latitude used by module and with reference to longitude with the MM4-FDDA databases of U.S. EPA consistent.
In practical application, initial wind field can be that a dimensional wind or one are made up of zone leveling wind field
Constant wind field.Zone leveling wind field can be by carrying out that vertical direction is average to aerological sounding meteorological data and temporal interpolation
Method obtains, and simply can also be specified by outside.If zone leveling wind field is obtained using the method calculated, then not only
Need to specify the atmosphere that calculate average wind value, and to point out the aerological sounding station for calculating average wind field value；Or adopt
With all station datas, 1/r is used2Interpolation (anti-square distance interpolation) produces the three-dimensional initial wind field of a spatial variations.Three
The foundation of dimension initial wind field establishes mode using existing three-dimensional initial wind field.
It can be used in present embodiment and meteorologic model or the grid of CSUMM mesoscale models generation are predicted by WRF/MM5
Wind field is as three-dimensional initial wind field.Operation is allowed to receive mesh spacing bigger than diagnostic model wind field and different vertical grid point
The forecast model wind field of cloth.This method allows multiple characteristics of wind field to be introduced in Diagnosis Field result, such as in ground gas
The lake surface vane stream that can not be measured as observation station and high-altitude backflow.
Step S300：According to the numerical weather forecast data of the meteorological department and geographical data to described three-dimensional first
Beginning wind field carries out just successive step, obtains the wind field after first successive step；
In present embodiment, the just successive step includes adjusting three-dimensional initial wind field according to landform kinematics effect
It is whole, according to slope flow effect three-dimensional initial wind field is adjusted and according to landform thermokinetics to the blocking effect of wind field to three-dimensional
Initial wind field is adjusted.The concrete mode of three kinds of adjustment is as follows：
1. landform kinematics effect
It is described three-dimensional initial wind field is adjusted according to landform kinematics effect including：According to landform kinematics effect pair
Vertical direction wind speed component is adjusted and horizontal direction wind speed component is adjusted according to landform kinematics effect；
It is according to the mode that landform kinematics effect is adjusted to vertical direction wind speed component：
2.1) the vertical component w of wind speed in cartesian coordinate system is calculated, calculation formula is：
Wherein, V is simulation region zone leveling wind speed, htFor Terrain Elevation,For htGradient, represent difference in height everywhere,
K is the exponential damping coefficient (exponential damping coefficient with the increase of atmospheric stability increase) relevant with atmospheric stability, and z is flute
Vertical coordinate under karr coordinate system, N are that frequency is drawn in Brunt-Visa, | V | it is the absolute value of zone leveling wind speed, g is gravity
Acceleration, θ are environment position temperature；
2.2) the vertical component W of wind speed under landform tracking coordinate system is calculated, calculation formula is：
Wherein, u, v be respectively under cartesian coordinate system wind speed in the component of x-axis horizontal direction and in y-axis horizontal direction
Component,For gradient of the Terrain Elevation in x-axis direction,For gradient of the Terrain Elevation in y-axis direction；X-axis horizontal direction is
East-west, y-axis horizontal direction is south-north direction；
Influence caused by wind field can enter landform kinematics effect by using to initial predicted wind field in the horizontal direction
Row minimum deviation method for uncertain is estimated.The mode being modified according to landform kinematics effect to horizontal direction wind field component
For：Horizontal Winds component is adjusted using iterative method using minimum deviation method for uncertain, until the three-dimensional deviation after adjustment is small
In setting value.
2. effect is flowed on slope
DIAMET Meteorological Models estimate the slope fluxion magnitude of complicated landform using empirical method.The direction of slope stream is typically recognized
To be consistent with the direction of draining.Slope flow vector is added to preliminary conjecture grid wind field, produces the wind field of amendment.
In present embodiment, it is according to the concrete mode that slope stream effect is adjusted to three-dimensional initial wind field：
1) slope stream wind speed caused by the stream effect of slope is calculated；Slope stream wind speed S calculation is as follows：
S=Se [1-exp (- x/Le)]1/3 (1)
Se=[hg (Δ θ/θ) sin α/(CD+k)]1/3 (2)
Le=h/ (CD+k) (3)
Wherein, Se is slope mobile equilibrium speed, LeFor equilibrium length yardstick, x is distance (x this at of the slope stream apart from top of the slope
It is unrelated with horizontal direction x directions), h is slope stream height, and g is acceleration of gravity, and Δ θ is temperature defect in position caused by environment, and θ is ring
Border position temperature, α are slope flow horizontal angle, CDFor ground drag coefficient, k is the entrainment coefficient at the top of the fluid layer of slope；
Formula (2) is substituted into formula (1), obtained：
It need to consider heat balance when calculating slope flow velocity and spending, relationship below be present
D (h Δ θ)/dt=Qhθ/(ρcpT) (5)
Wherein, t represents time, QhFor sensible heat flux, ρ is atmospheric density, cpFor heat capacity of air, T is air themperature；
Equation below (flowing direction along slope, it is assumed that d/dt=d/dx) is obtained to formula (5) integration：
H Δs θ=Qhθx/(ρcpT) (6)
(6) substitution (4) formula is obtained into slope flow velocity degree calculation formula is：
S={ [Qhgx sinα/[(ρcpT)(CD+k)]}1/3[1-exp(-x/Le)]1/3 (7)
Le=h/ (CD+k)
Formula (7) is wave current velocity expression, and reality need to calculate respectively when calculating according to mountain upslope flow and descending stream, for
Descending stream：In formula (7), there is related parameter as follows：
CD=K=4 × 10-2
Le=h/ (CD+ k) in h=0.05 Δs Z
Sin α=minimum (sin α, Δ Z/x)
Δ Z represents slope stream height (discrepancy in elevation of top of the slope to hillside fields), and other specification is by mode computation, such formula (7)
Calculate descending flow velocity degree；
For mountain upslope flow, in formula (7), there is related parameter as follows：
(CD+ k)~1
X sin α=Δ Z
Due to the exp (- x/L in formula (7)e)≤≤ 1, formula (7) can be written as following formula：
2) wind speed is flowed according to slope and slope flow horizontal angle [alpha] calculates the slope stream wind speed u of x-axis horizontal directionsWith y-axis level side
To slope stream wind speed vs；
3) slope stream wind speed is added to three-dimensional initial wind field, the horizontal wind speed of three-dimensional initial wind field is adjusted, adjusted
Formula be：
u1'=u1+us
v1'=v1+vs
Wherein, u1' and v1' it is respectively that the x-axis horizontal direction after effect is modified to three-dimensional initial wind field is flowed according to slope
The wind speed component of wind speed component and y-axis horizontal direction, u1And v1The wind speed of x-axis horizontal direction point in respectively three-dimensional initial wind field
The wind speed component of amount and y-axis horizontal direction.
3. blocking effect
In present embodiment, landform heat power has carried out parameter to the blocking effect of wind field with local Froude number forms formula
Change.The local Froude number Fr of each mesh point in three-dimensional initial wind field is calculated, calculation formula is：
Δht=(hmax)ij-(z)ijk
Wherein, V is the wind speed of mesh point, and N is that frequency (Brunt is drawn in Brunt-VisaFrequency), Δ htTo have
Imitate barrier height, (hmax)ijThe maximum Terrain Elevation of (TERRAD), (z) in mesh point (i, j) radius of influenceijkFor high dead level k
The height of middle mesh point (i, j)；
Froude numbers are calculated to each mesh point, judge whether the local Froude number Fr of mesh point is less than setting afterwards
Critical Froude number, if it is not, be not adjusted to three-dimensional initial wind field then, if so, then judge mesh point wind field whether have to
On component, if so, the wind direction of mesh point then is adjusted into the direction tangent with landform, if it is not, then not to three-dimensional initial wind field
It is adjusted.
Step S400：It is modified, is obtained most according to the three-dimensional grid wind field after successive step at the beginning of the data of making weather observations
Whole dimensional wind.
The step is using diagnostic mode processing wind field, will observe data and imported into what S300 was just formed after successive step
The process of grid wind field.The step for include following four sub-steps：
Interpolation
Surface wind observes data vertical extrapolation
Smoothing processing
O ' the Brien adjustment of vertical velocity
Deviation minimizes
1. interpolation processing
Observation data in weather station are imported into step 1 wind field using apart from inverse ratio weighting method.According to data of making weather observations
To after first successive step wind field carry out interpolation processing interpolation formula be：
Wherein, (u, v) '2For the initial wind speed component of the horizontal direction of mesh point in the wind field after interpolation processing, (u, v)1
For the wind speed component of the horizontal direction of mesh point in the wind field after the just successive step, k is the identification number of weather station,
(uobs,vobs)kFor the observation wind speed component of the horizontal direction of k-th of weather station, R is the wind field after default just successive step
Weighting parameters (weight number), specified by user；RkFor the distance of k-th of weather station to mesh point.
This interpolation method allows to carry out big weight assignment to data of making weather observations in the region near from observation station, and in no gas
As the region of observation data, the wind field after first successive step is occupied an leading position in weight assignment.The method of weighting, each vertical
Used in layer.Determine that earth's surface meteorological measuring is used for minimum regardless of the option of application earth's surface observation Data Extrapolation
The calculating of layer wind field layer is only most suitable.If exceed maximum effect radius with a distance from the mesh point of observation station from one, then
This observation data can not have into row interpolation.
2. Surface wind observes data vertical extrapolation
Before space interpolation is carried out to the wind field of horizontal direction, the data of each weather station of earth's surface can be to higher
Layer is extrapolated.Specifically include：Power-law equation extrapolation, theory of similarity extrapolation, the custom dimensions factor are extrapolated etc..
If surface weather station is closer apart from aerological sounding weather station, and aerological sounding weather station can provide effective spy
Empty data, then avoid the need for carrying out vertical extrapolation to the data of surface weather station.
● power-law equation is extrapolated：
uz=um·(z/zm)P (10)
Wherein, Z be DIAMET grid bodies central point height, m；zmFor observed altitude, m；uzTo be extrapolated at height z
Wind speed u components, m/s；umFor the wind speed u components of observation, m/s；P is power law index.
Similar formula can be used for the v components of calculation of wind speed.
Carried according to Douglas and kessler1988 in DWM, P=0.143 is applied to top, P=
0.286 is applied to water surface.The mean sea level of waters grid is labeled as zero.
● the custom dimensions factor is extrapolated
A series of scale factors are given, more than DIAMET earth's surfaces each layer is corresponding one, all can be used to NZ layers from the 2nd layer
Lower formula calculates wind field：
ui=u1·FEXTRPi (11)
Wherein, i represents the DIAMET wind fields number of plies (i=2,3 ..., NZ)；u1Represent the u components of the Wind Fields of Layer 1；uiTable
Show the u components of Layer i wind fields；FEXTRPiRepresent i-th layer of scale factor that user gives.Same formula can be used for
Calculate the v components of wind field.
● the theory of similarity is extrapolated
The method of the third extrapolation is based on van Ulden and Holtslg (1985) research, and this Extrapolation method uses phase
The influence of earth's surface wind speed and direction is expanded to high dead level with data have been observed like theory.So, more than earth's surface or mixed layer
The wind speed and direction of 200m each layer will all change.Here is the calculating of van Ulden and Holtslag (1985) Extrapolation method
Formula, wind shift are provided with the change of height by equation (12)：
D (z)/D (h)=d1[1-exp(-d2z/h)] (12)
Herein, D (z) represents torsion angle the z of floor height center at, the torsion angle at D (h) expression reference altitudes h, and experience is normal
Measure d1=1.58, d2=1.0.
It is based on the similar reason of not peaceful-difficult to understand cloth Hough of surface layer that van Ulden calculate with wind profile that Holtslag is proposed
Opinion.Due to relevant with stability, it is necessary to utilize formula 14 based on height and Mo-length difficult to understand) and the stable letter of formula (15) calculating
Number.By stability function and measurement height, the roughness length of grid is updated to equation (12) where layer centre-height and weather station
Wind speed at computation layer centre-height.Wind speed and direction after correction can finally convert back u and v components, remain follow-up insert
Used in value program.After windup-degree is calculated, to increase the angle, (wind angle moves the wind direction in the Northern Hemisphere clockwise
It is dynamic), and be to subtract the angle (moving counterclockwise) by the wind direction in the Southern Hemisphere.
Formula (1.2-22) gives theory of similarity equation and carrys out calculation of wind speed profile：
Wherein, U (z) refers to the wind speed at DIAMET layers center；U(z1) refer to the wind speed that airspeedometer is highly located；Z0Refer to coarse length
Degree；Z1Refer to airspeedometer height；ΨMRefer to stability function.
Formula (1.2-23) gives the expression formula of stability function under unstable state：
X=(1-16z/L)1/4 (15)
Formula (16) is the expression formula of stability function at steady state：
ψM=-17 [1-exp (- 0.29z/L)] (16)
3. smoothing processing
The final wind field of observation data generation is added in wind field after first successive step to be smoothed, to reduce in wind field
Caused discontinuity.Smoothing processing formula in DIAMET Meteorological Models is as follows：
Smoothing processing formula is：
(ui,j)2"=0.5ui,j+0.125(ui-1,j+ui+1,j+ui,j-1+ui,j+1)
(vi,j)2"=0.5vi,j+0.125(vi-1,j+vi+1,j+vi,j-1+vi,j+1)
Its, (ui,j)2″、(vi,j)2" be respectively smoothing processing after in the horizontal direction at mesh point (i, j) place x-axis direction
Wind speed component and the wind speed component in y-axis direction；ui,j、vi,jAt mesh point (i, j) place after interpolation processing respectively before smoothing processing
Horizontal direction in x-axis direction wind speed component and y-axis direction wind speed component (according to the calculation formula meter of above-mentioned interpolation processing
Air quantity after the first successive step calculated)；
4. the calculating of vertical speed
In DIAMET Meteorological Models, there are two kinds of options for being used for calculating vertical speed.The first, passes through not compression quality
Conservation equation directly calculates vertical speed from smoothed Horizontal Winds component, and the vertical velocity component calculated is final
Vertical velocity component, the Horizontal Winds after smoothing processing are final Horizontal Winds.Second method is not compressed in basis
After mass-conservation equation calculates vertical velocity component from the Horizontal Winds component after smoothing processing, by adjusting vertical speed
Profile and make the value at the top of model area be zero, then, wind field horizontal component is adjusted to again and new vertical velocity field is kept
Quality compatibility.
The first, vertically-supplying air is calculated according to not compression quality conservation equation from the Horizontal Winds component after smoothing processing
Fast component W1(z) formula is：
Wherein, u ", v " is respectively the wind speed component in x-axis direction and the wind in y-axis direction in the horizontal direction after smoothing processing
Fast component.
If selected for first method, then the compatible vertical speed of quality can be as final vertical speed (i.e. w=
w1), moreover, making again to be adjusted wind field horizontal component in this way.Final Horizontal Winds are exactly by smoothly
The wind field obtained after reason formula is smooth.
Godden and Lurmann (1993) is it is proposed that this computational methods there may come a time when that clathrum can be caused top hangs down
Straight speed is unpractical excessive.In order to avoid the appearance of this problem, second method can be used.
Second, calculated according to not compression quality conservation equation from the Horizontal Winds component after smoothing processing vertically
After wind speed component, the vertical velocity calculated is adjusted, the vertical velocity component W after being adjusted2(z) calculation formula
For：
W2(z)=W1(z)-(z/ztop)W1(z=ztop)
Wherein, z be cartesian coordinate system under vertical coordinate, ztopFor the overhead height of simulation region, W1(z=ztop) it is root
The vertical speed at the top of simulation region calculated according to not compression quality conservation equation from the Horizontal Winds component after smoothing processing
Component；
Calculate the vertical velocity component W after adjustment2(z) after, two horizontal directions are adjusted using minimum deviation method for uncertain
Wind speed component, minimum deviation method for uncertain keep vertical component be W2(z) in the case of constant, using alternative manner by two
The deviation adjusting of the wind speed component of horizontal direction is to the deviation threshold less than setting, two water after being adjusted using deviation method
Square to horizontal component be final Horizontal Winds component.
In present embodiment, the vertical velocity component W after being adjusted2(z) after, two are adjusted using minimum deviation method for uncertain
The adjustment of the wind speed component of individual horizontal direction, ensure that vertical direction wind speed is constant, pass through iteration adjustment horizontal wind speed component so that
The deviation of each mesh point is less than the maximum deflection difference value of user's setting, that is, meets
Wherein, u and v is Horizontal Winds component；W is that the vertical speed that landform is tracked in coordinate system is W2(z)；ε is maximum
Tolerance value.
The concrete mode that deviation minimizes the following is：
For mesh point (i, j, k), the wind speed component deviation D of the horizontal wind direction at calculating mesh point (i, j, k) placeijk, meter
Calculating formula is：
Wherein, zk+1/2Middle point coordinates between vertical direction kth layer and kth+1, zk-1/2For vertical direction kth layer with
Middle point coordinates between kth -1, zk+1/2-zk-1/2The floor height in vertical direction is represented, Δ x and Δ y are illustrated respectively in x-axis and y-axis
The length of side of grid cell on direction；The subscript of w, u, v appeared in formula represent three-dimensional coordinate, such as wi,j,k+1/2Represent flat
Sliding processing recoil is designated as the vertical direction wind speed at (i, j, k+1/2) place.
Afterwards, the horizontal wind speed component of surrounding's grid cell of mesh point (i, j, k) is adjusted, makes the deviation of mesh point
DijkMake deviation DijkLess than deviation threshold ε, adjustment formula is：
(unew)i+1,j,k=ui+1,j,k+uadj
(unew)i-1,j,k=ui-1,j,k-uadj
(vnew)i+1,j,k=vi,j+1,k+vadj
(vnew)i,j-1,k=vi,j-1,k-vadj
Wherein, uadjAnd vadjThe adjustment amount of wind speed component respectively in x-axis and y-axis direction；unew、vnewAfter representing adjustment
Wind speed component in x-axis and y-axis direction, unewAnd vnewSubscript represent three-dimensional coordinate.
Every time after the deviation of a mesh point is eliminated, deviation can be then being produced around it at point, iteration uses above-mentioned
Adjustment mode, until the deviation of all grid cells is respectively less than the deviation threshold ε set.
Calculate more than, by first successive step of the three-dimensional initial wind field through landform kinematics, slope stream and blockage effect, and insert
Value, smooth etc. obtain final three-dimensional meteorological wind field after Adjustable calculation again.
Obviously, those skilled in the art can carry out the essence of various changes and modification without departing from the present invention to the present invention
God and scope.So, if these modifications and variations of the present invention belong to the scope of the claims in the present invention and its equivalent technology
Within, then the present invention is also intended to comprising including these changes and modification.
Claims (10)
1. a kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation, comprises the following steps：
Step 1: centered on nuclear facilities, the simulation region scope of nuclear facilities accident is determined, obtains the meteorological data and ground of simulation region
Environmental data is managed, and the three-dimensional initial wind field of simulation region is established according to the meteorological data and geographic position data of acquisition；The gas
Image data includes make weather observations data and the numerical weather forecast data of meteorological department of the weather station of simulation region；Describedly
Reason environmental data includes latitude and longitude information, underlying surface elevation, vegetation, water body distribution characteristics and the land use of simulation region；It is three-dimensional
Initial wind field is multi-layer net wind field；
Step 2: according to the numerical weather forecast data of the meteorological department and geographical data to the three-dimensional initial wind field
Just successive step is carried out, obtains the wind field after first successive step；The just successive step is included according to landform kinematics effect at the beginning of three-dimensional
Beginning wind field is adjusted, three-dimensional initial wind field is adjusted and according to resistance of the landform thermokinetics to wind field according to slope stream effect
Plug effect is adjusted to three-dimensional initial wind field；
Step 3: being handled according to the data of making weather observations the wind field after first successive step, final wind field is obtained；The place
Reason includes carrying out the wind field after first successive step successively the adjustment of interpolation processing, smoothing processing and vertical velocity component.
2. the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation according to claim 1, in step 2,
It is described three-dimensional initial wind field is adjusted according to landform kinematics effect including：According to landform kinematics effect to vertical direction
Wind speed component is adjusted and horizontal direction wind speed component is adjusted according to landform kinematics effect；
It is according to the mode that landform kinematics effect is adjusted to vertical direction wind speed component：
2.1) the vertical component w of wind speed in cartesian coordinate system is calculated, calculation formula is：
<mrow>
<mi>w</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>&CenterDot;</mo>
<mo>&dtri;</mo>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>k</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mfrac>
<mi>N</mi>
<mrow>
<mo>|</mo>
<mi>V</mi>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>N</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mi>g</mi>
<mi>&theta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&theta;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>z</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mrow>
Wherein, V is simulated domain mean wind speed, htFor Terrain Elevation,For htGradient, represent difference in height everywhere, k be with it is big
The relevant exponential damping coefficient of gas stability, z are the vertical coordinate under cartesian coordinate system, and N is that frequency is drawn in Brunt-Visa, |
V | it is the absolute value of zone leveling wind speed, g is acceleration of gravity, and θ is environment position temperature；
2.2) the vertical component W of wind speed under landform tracking coordinate system is calculated, calculation formula is：
<mrow>
<mi>W</mi>
<mo>=</mo>
<mi>w</mi>
<mo>-</mo>
<mi>u</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>v</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
</mrow>
Wherein, u, v are respectively component of the wind speed in x-axis horizontal direction and the component in y-axis horizontal direction under cartesian coordinate system,For gradient of the Terrain Elevation in x-axis direction,For gradient of the Terrain Elevation in y-axis direction；X-axis horizontal direction is thing
Trend, y-axis horizontal direction is south-north direction；
It is according to the mode that landform kinematics effect is modified to horizontal direction wind field component：Using minimum deviation method for uncertain pair
Horizontal Winds component is adjusted using iterative method, until the three-dimensional deviation after adjustment is less than setting value.
3. the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation according to claim 1, in step 2,
It is according to the mode that slope stream effect is adjusted to three-dimensional initial wind field：
1) slope stream wind speed caused by the stream effect of slope is calculated；Slope stream wind speed S calculation formula is：
S={ [Qh g x sinα/[(ρcp T)(CD+k)]}1/3[1-exp(-x/Le)]1/3
Le=h/ (CD+k)
Wherein, QhFor sensible heat flux, g is acceleration of gravity, and x flows the distance apart from top of the slope for slope, and α is slope flow horizontal angle, and ρ is
Atmospheric density, cpFor heat capacity of air, T is air themperature, CDFor ground drag coefficient, k is the entrainment coefficient at the top of the fluid layer of slope, Le
For equilibrium length yardstick, h is slope stream height；
2) wind speed is flowed according to slope and slope flow horizontal angle [alpha] calculates the slope stream wind speed u of x-axis horizontal directionsWith y-axis horizontal direction
Slope stream wind speed vs；
3) slope stream wind speed is added to three-dimensional initial wind field, the horizontal wind speed of three-dimensional initial wind field is adjusted, the public affairs of adjustment
Formula is：
uZ'=u1+us
v1'=v1+vs
Wherein, u1' and v1' it is respectively the wind speed that the x-axis horizontal direction after effect is modified to three-dimensional initial wind field is flowed according to slope
The wind speed component of component and y-axis horizontal direction, u1And v1The wind speed component and y of x-axis horizontal direction in respectively three-dimensional initial wind field
The wind speed component of axle horizontal direction.
4. the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation according to claim 3, in step 1),
When slope stream is descending stream,
CD=K=4 × 10-2
H=0.05 Δs Z
Sin α=minimum (sin α, Δ Z/x)
When slope stream is mountain upslope flow, (CD+ k)~1, x sin α=Δ Z；
Wherein, Δ Z is slope stream height.
5. the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation according to claim 1, in step 2,
It is described the blocking effect of wind field is adjusted to three-dimensional initial wind field according to landform thermokinetics including：
The local Froude number Fr of each mesh point in three-dimensional initial wind field is calculated, calculation formula is：
<mrow>
<mi>F</mi>
<mi>r</mi>
<mo>=</mo>
<mfrac>
<mi>V</mi>
<mrow>
<msub>
<mi>N&Delta;h</mi>
<mi>t</mi>
</msub>
</mrow>
</mfrac>
</mrow>
Δht=(hmax)ij-(z)ijk
Wherein, V is the wind speed of mesh point, and N is that frequency, Δ h are drawn in Brunt-VisatFor effective barrier height, (hmax)ijFor grid
Maximum Terrain Elevation in point (i, j) radius of influence, (z)ijkFor the height of mesh point (i, j) in high dead level k；
Judge whether the local Froude number Fr of mesh point is less than the critical Froude number of setting, if it is not, then not to three-dimensional initial wind
Field is adjusted, if so, whether have upward component, if so, then adjusting the wind direction of mesh point if then judging the wind field of mesh point
For the direction tangent with landform, if it is not, not being adjusted to three-dimensional initial wind field then.
6. the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation according to claim 1, its feature exist
In in step 3, according to making weather observations, data are to the interpolation formula of the wind field progress interpolation processing after first successive step：
<mrow>
<msub>
<msup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>,</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>&prime;</mo>
</msup>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>,</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mn>1</mn>
</msub>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>k</mi>
</munder>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>k</mi>
</msub>
<mrow>
<msup>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mrow>
<mfrac>
<mn>1</mn>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>k</mi>
</munder>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</mfrac>
</mrow>
Wherein, (u, v) '2For the initial wind speed component of the horizontal direction of mesh point in the wind field after interpolation processing, (u, v)1For institute
State the wind speed component of the horizontal direction of mesh point in the wind field after first successive step, k is the identification number of weather station, (uobs,
vobs)kFor the observation wind speed component of the horizontal direction of k-th of weather station, R is the power of the wind field after default just successive step
Tuple；RkFor the distance of k-th of weather station to mesh point；
According to making weather observations, data are smoothed to wind field after interpolation processing, and smoothing processing formula is：
(ui,j)2"=0.5ui,j+0.125(ui-1,j+ui+1,j+ui,j-1+ui,j+1)
(vi,j)2"=0.5vi,j+0.125(vi-1,j+vi+1,j+vi,j-1+vi,j+1)
Its, (ui,j)2″、(vi,j)2" be respectively smoothing processing after in the horizontal direction at mesh point (i, j) place x-axis direction wind speed
Component and the wind speed component in y-axis direction；ui,j、vi,jIn the water at mesh point (i, j) place after interpolation processing respectively before smoothing processing
Square to the wind speed component in middle x-axis direction and the wind speed component in y-axis direction；
According to making weather observations, data include to the mode of the wind field progress vertical velocity adjustment after smoothing processing：According to not compressing matter
Amount conservation equation calculates vertical velocity component from the Horizontal Winds component after smoothing processing, the vertical velocity component calculated
For final vertical velocity component, the Horizontal Winds after smoothing processing are final Horizontal Winds, or,
After according to compression quality conservation equation does not calculate vertical velocity component from the Horizontal Winds component after smoothing processing,
The vertical velocity calculated is adjusted, the vertical velocity component after adjustment is final vertical velocity component.
7. the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation according to claim 6, its feature exist
In calculating vertical velocity component W from the Horizontal Winds component after smoothing processing according to not compression quality conservation equation1(z)
Formula be：
<mrow>
<mfrac>
<mrow>
<msup>
<mi>du</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mi>dv</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
</mrow>
<mrow>
<mi>d</mi>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>dW</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>z</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
</mrow>
Wherein, u ", v " is respectively the wind speed component in x-axis direction and the wind speed point in y-axis direction in the horizontal direction after smoothing processing
Amount.
8. the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation according to claim 7, its feature exist
In, after according to compression quality conservation equation does not calculate vertical velocity component from the Horizontal Winds component after smoothing processing,
The vertical velocity calculated is adjusted, the vertical velocity component W after being adjusted2(z) calculation formula is：
W2(z)=W1(z)-(z/ztop)W1(z=ztop)
Wherein, z be cartesian coordinate system under vertical coordinate, ztopFor the overhead height of simulation region, W1(z=ztop) it is according to not
The vertical velocity component at the top of simulation region that compression quality conservation equation calculates from the Horizontal Winds component after smoothing processing；
Calculate the vertical velocity component W after adjustment2(z) after, the wind speed of two horizontal directions is adjusted using minimum deviation method for uncertain
Component, minimum deviation method for uncertain are W keeping vertical component2(z) in the case of constant, using alternative manner by two level sides
To wind speed component deviation adjusting to less than setting deviation threshold, using deviation method adjust after two horizontal directions
Horizontal component be final Horizontal Winds component.
9. the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation according to claim 8, its feature exist
In, calculate adjustment after vertical velocity component W2(z) after, the wind speed of two horizontal directions is adjusted using minimum deviation method for uncertain
The mode of the adjustment of component is：
Calculate the wind speed component deviation D of the horizontal wind direction at each mesh point (i, j, k) place in grid cellijk, calculation formula
For：
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>z</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mi>&Delta;</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mi>&Delta;</mi>
<mi>y</mi>
</mrow>
</mfrac>
</mrow>
Wherein, zk+1/2-zk-1/2The floor height in vertical direction is represented, Δ x and Δ y are illustrated respectively in grid list on x-axis and y-axis direction
The length of side of member；
The horizontal wind speed component of surrounding's grid cell of mesh point (i, j, k) is adjusted, makes the deviation D of mesh pointijkLess than setting
Deviation threshold ε, adjustment formula be：
(unew)i+1,j,k=ui+1,j,k+uadj
(unew)i-1,j,k=ui-1,j,k-uadj
(vnew)i+1,j,k=vi,j+1,k+vadj
(vnew)i,j-1,k=vi,j-1,k-vadj
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>a</mi>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mi>&Delta;</mi>
<mi>x</mi>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>a</mi>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mi>&Delta;</mi>
<mi>y</mi>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
Wherein, uadjAnd vadjThe adjustment amount of wind speed component respectively in x-axis and y-axis direction；
Iteration uses above-mentioned adjustment mode, until the deviation of all grid cells is respectively less than the deviation threshold ε set.
10. the wind field diagnostic method in a kind of nuclear facilities Accident Off-site Consequence evaluation according to claim 6, its feature exist
In before space interpolation is carried out to the wind field of horizontal direction, the data of each weather station of earth's surface can be entered to higher level
Row extrapolation, including：Power-law equation extrapolation, theory of similarity extrapolation, the custom dimensions factor are extrapolated.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610453592.1A CN107526910A (en) | 2016-06-21 | 2016-06-21 | A kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610453592.1A CN107526910A (en) | 2016-06-21 | 2016-06-21 | A kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107526910A true CN107526910A (en) | 2017-12-29 |
Family
ID=60735168
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610453592.1A Pending CN107526910A (en) | 2016-06-21 | 2016-06-21 | A kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107526910A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111582547A (en) * | 2020-04-09 | 2020-08-25 | 中国科学院国家空间科学中心 | Method for acquiring wind field distribution at different places by using wind field data set |
| CN112037290A (en) * | 2020-08-28 | 2020-12-04 | 中科三清科技有限公司 | Data processing method, device and equipment |
| CN116842877A (en) * | 2023-09-04 | 2023-10-03 | 江苏省环境监测中心 | Small-scale three-dimensional wind field reconstruction algorithm based on multi-source data comprehensive utilization |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101726748A (en) * | 2008-10-27 | 2010-06-09 | 北京广利核系统工程有限公司 | Method for displaying evaluation data of nuclear radiation aftermath |
| CN101834001A (en) * | 2010-04-12 | 2010-09-15 | 中国人民解放军防化指挥工程学院 | Optimal control method aiming at nuclear accident loss and emergency action cost |
| CN102628876A (en) * | 2012-02-13 | 2012-08-08 | 甘肃省电力公司风电技术中心 | Ultra-short term prediction method comprising real-time upstream and downstream effect monitoring |
| CN103348132A (en) * | 2011-02-07 | 2013-10-09 | 维斯塔斯风力系统有限公司 | Wind turbine having heliplatform arrangement and method of using same |
| CN103956196A (en) * | 2014-03-31 | 2014-07-30 | 中国核电工程有限公司 | Liquid film evaporation cooling plate of liquid water collecting and cooling apparatus |
| CN204294017U (en) * | 2014-11-13 | 2015-04-29 | 大亚湾核电运营管理有限责任公司 | A kind of nuclear island air cleaner |
| CN104750953A (en) * | 2013-12-26 | 2015-07-01 | 中国辐射防护研究院 | Method for simulating middle-and-small-scale airborne substance atmosphere transportation ensemble diffusion |
| CN104900826A (en) * | 2015-04-21 | 2015-09-09 | 天津瑞能电气有限公司 | Wind power generation wind-electricity turbine pitch system battery case and temperature control method thereof |
| CN104950129A (en) * | 2014-03-31 | 2015-09-30 | 江苏核电有限公司 | Nuclear power station meteorological information emergency early-warning complete set equipment and method |
-
2016
- 2016-06-21 CN CN201610453592.1A patent/CN107526910A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101726748A (en) * | 2008-10-27 | 2010-06-09 | 北京广利核系统工程有限公司 | Method for displaying evaluation data of nuclear radiation aftermath |
| CN101834001A (en) * | 2010-04-12 | 2010-09-15 | 中国人民解放军防化指挥工程学院 | Optimal control method aiming at nuclear accident loss and emergency action cost |
| CN103348132A (en) * | 2011-02-07 | 2013-10-09 | 维斯塔斯风力系统有限公司 | Wind turbine having heliplatform arrangement and method of using same |
| CN102628876A (en) * | 2012-02-13 | 2012-08-08 | 甘肃省电力公司风电技术中心 | Ultra-short term prediction method comprising real-time upstream and downstream effect monitoring |
| CN104750953A (en) * | 2013-12-26 | 2015-07-01 | 中国辐射防护研究院 | Method for simulating middle-and-small-scale airborne substance atmosphere transportation ensemble diffusion |
| CN103956196A (en) * | 2014-03-31 | 2014-07-30 | 中国核电工程有限公司 | Liquid film evaporation cooling plate of liquid water collecting and cooling apparatus |
| CN104950129A (en) * | 2014-03-31 | 2015-09-30 | 江苏核电有限公司 | Nuclear power station meteorological information emergency early-warning complete set equipment and method |
| CN204294017U (en) * | 2014-11-13 | 2015-04-29 | 大亚湾核电运营管理有限责任公司 | A kind of nuclear island air cleaner |
| CN104900826A (en) * | 2015-04-21 | 2015-09-09 | 天津瑞能电气有限公司 | Wind power generation wind-electricity turbine pitch system battery case and temperature control method thereof |
Non-Patent Citations (6)
| Title |
|---|
| JAY WU .ETC: "Three-dimensional dose evaluation system using real-time wind field information for nuclear accidents in Taiwan", 《NUCLEAR INSTRUMENTS AND METHODS IN PHYSICS RESEARCH A：ACCELERATORS,SPECTROMETERS,DETECTORS AND ASSOCIATED EQUIPMENT》 * |
| JOSEPH S.SCIRE .ETC: "《A User"s Guide for the CALMET Meteorological Model（Version 5）》", 31 January 2000 * |
| 乔清党等: ""核事故应急气载放射性物质扩散的精细化模拟方法研究"", 《核安全》 * |
| 伯鑫等: "大气扩散CALPUFF模型技术综述", 《环境监测管理与技术》 * |
| 刘爱华: "核设施气载核素浓度及辐射剂量计算研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 胡啸峰等: "核事故下源项排放率对放射性物质大气扩散及沉降的影响", 《科学技术与工程》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111582547A (en) * | 2020-04-09 | 2020-08-25 | 中国科学院国家空间科学中心 | Method for acquiring wind field distribution at different places by using wind field data set |
| CN112037290A (en) * | 2020-08-28 | 2020-12-04 | 中科三清科技有限公司 | Data processing method, device and equipment |
| CN112037290B (en) * | 2020-08-28 | 2021-05-11 | 中科三清科技有限公司 | Data processing method, device and equipment |
| CN116842877A (en) * | 2023-09-04 | 2023-10-03 | 江苏省环境监测中心 | Small-scale three-dimensional wind field reconstruction algorithm based on multi-source data comprehensive utilization |
| CN116842877B (en) * | 2023-09-04 | 2023-11-28 | 江苏省环境监测中心 | Small-scale three-dimensional wind field reconstruction algorithm based on multi-source data comprehensive utilization |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| He et al. | High-resolution dataset of urban canopy parameters for Beijing and its application to the integrated WRF/Urban modelling system | |
| Huang et al. | Impact of different reanalysis data on WRF dynamical downscaling over China | |
| CN107678075A (en) | A kind of urban heat land effect monitoring method and system based on domestic satellite | |
| Järvinen et al. | Estimation of ECHAM5 climate model closure parameters with adaptive MCMC | |
| CN107390298A (en) | A kind of analogy method and device of Complex Mountain underlying surface strong wind | |
| Kang et al. | Development of an observation processing package for data assimilation in KIAPS | |
| CN102354348A (en) | Watershed scale soil moisture remote sensing data assimilation method | |
| Challa et al. | A simulation study of mesoscale coastal circulations in Mississippi Gulf coast | |
| Li et al. | Numerical simulation study of the effect of buildings and complex terrain on the low-level winds at an airport in typhoon situation | |
| CN110175375A (en) | A kind of earth's surface Calculation method for solar radiation based on deep learning | |
| Hamdi et al. | Evaluating the performance of SURFEXv5 as a new land surface scheme for the ALADINcy36 and ALARO-0 models | |
| Chen et al. | Refined urban canopy parameters and their impacts on simulation of urbanization-induced climate change | |
| Krishnamurti | An experiment in numerical prediction in equatorial latitudes | |
| Rauscher et al. | Domain choice in an experimental nested modeling prediction system for South America | |
| Moussiopoulos et al. | A refined diagnostic wind model | |
| KR101504100B1 (en) | Physics parameterization package construction method of numerical weather prediction model using cubed-sphere horizontal grid system based on lorenz vertical discretization method and hardware device implementing the same | |
| CN112925870B (en) | Population spatialization method and system | |
| CN107526910A (en) | A kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation | |
| CN106054283A (en) | Method and device for inverting upper troposphere and lower stratosphere wind fields | |
| Hadi | Diagnosis of the best method for wind speed extrapolation | |
| CN115994496B (en) | Urban park high resolution atmospheric CO 2 Numerical simulation method of concentration three-dimensional field | |
| CN113985490B (en) | Method and device for carrying out surface gravity simulation by utilizing terrain and crust density data | |
| Wang et al. | Evaluation of five planetary boundary layer schemes in WRF over China's largest semi-fixed desert | |
| Wang et al. | A new horizontal mixing-length formulation for numerical simulations of tropical cyclones | |
| Qian et al. | The effect of grid spacing and domain size on the quality of ensemble regional climate downscaling over South Asia during the northeasterly monsoon |
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 |