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 PDF

Info

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
Application number
CN201610453592.1A
Other languages
Chinese (zh)
Inventor
徐向军
赵多新
姚仁太
郝宏伟
黄杰
张俊芳
吕明华
黄莎
胡继民
张芳
高卫华
杨彪
韩旻晨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute for Radiation Protection
Original Assignee
China Institute for Radiation Protection
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Institute for Radiation Protection filed Critical China Institute for Radiation Protection
Priority to CN201610453592.1A priority Critical patent/CN107526910A/en
Publication of CN107526910A publication Critical patent/CN107526910A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information 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

A kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation
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>&amp;CenterDot;</mo> <mo>&amp;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>&amp;lsqb;</mo> <mrow> <mo>(</mo> <mfrac> <mi>g</mi> <mi>&amp;theta;</mi> </mfrac> <mo>)</mo> </mrow> <mfrac> <mrow> <mi>d</mi> <mi>&amp;theta;</mi> </mrow> <mrow> <mi>d</mi> <mi>z</mi> </mrow> </mfrac> <mo>&amp;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>&amp;part;</mo> <msub> <mi>h</mi> <mi>t</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>-</mo> <mi>v</mi> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>h</mi> <mi>t</mi> </msub> </mrow> <mrow> <mo>&amp;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&amp;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>&amp;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>&amp;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>&amp;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>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> </mrow> <mrow> <mi>d</mi> <mi>x</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>dv</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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.
CN201610453592.1A 2016-06-21 2016-06-21 A kind of wind field diagnostic method in nuclear facilities Accident Off-site Consequence evaluation Pending CN107526910A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (9)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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&#39;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