首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
A pair of parallel cold wires separated in either the vertical or lateral direction was used to obtain the three components x, y, z of the temperature derivative in the streamwise, lateral and vertical directions, respectively. The average absolute skewness values of x and z are nonzero and approximately equal, while the skewness of y is approximately zero. These results appear to be consistent with the presence of a large, three-dimensional organised structure in the surface layer. There is an apparent low-frequency contamination in the spectral density of y and z due mainly to small errors in estimating the sensitivity of the cold wires. The temperature derivatives were high-pass filtered, the filter being set to remove possible contributions from the large structure and to minimise low-frequency sensitivity contamination. The filtered rms ratios \~x/\~y and \~x/\~z were in the range 0.7 to 0.9, a result in qualitative agreement with that obtained in the laboratory boundary layer by Sreenivasan et al. (1977). The skewness of filtered x or z is negligible, consistent with local isotropy of small-scale temperature fluctuations and in support of the high wavenumber spectral isotropy discussed in Antonia and Chambers (1978).  相似文献   

2.
This is one of a series of papers on the Askervein Hill Project. It presents results on the variations in mean wind speed at fixed heights (z) above the ground from linear arrays of anemometer posts and towers. Most of the data are for z = 10 m but some are for z = 3 m. Selected and directionally grouped data from the 55 Mean Flow runs are presented together with mean flow data from Askervein '83 Turbulence runs. Comparisons are made between the data and guideline estimates of fractional speed-up ratio at hilltop locations and between the data and MS3DJH/3 model predictions along the tower lines. There is good agreement in most cases.  相似文献   

3.
It is well known that if turbulent mass convection is modelled as diffusion, errors result unless trajectories from the source (ath) to the point of observation (z p ) comprise many statistically-independent segments (Taylor, 1921). We show that this is not guaranteed merely by the Lagrangian timescale () at the source being small (e.g., source at ground), but that a better criterion istmax[(h), (z p )], wheret is a typical travel time toz p .  相似文献   

4.
It is shown that for the purpose of trajectory simulation, the vertical velocityw L (t) of a fluid element, which is moving in a system (such as a forest canopy, or the unstably stratified atmospheric surface layer) whose turbulent velocity scale w is height-dependent, must be chosen from a frequency-distribution which is asymmetric aboutw L = 0. If the gradient w /z varies only slowly with height, correct trajectories may be obtained by adding a bias (where L is the length scale) to a fluctuating velocity chosen from a symmetric distribution with variance w 2(z).  相似文献   

5.
Recently Wilson and Flesch (Boundary-Layer Meteorology, 84, 411-426, 1997) suggested that the average increment d z to the orientation = arctan(w/u) of the Lagrangian velocity-fluctuation vector can be used to distinguish the better Lagrangian stochastic models within the well-mixed class. Here it is demonstrated that the specification of d z constitutes neither a sufficient or universally applicable criterion to distinguish the better Lagrangian stochastic models within the well-mixed class. The hypothesis made by Wilson and Flesch that Lagrangian stochastic models with /PE irrotational are zero-spin models, having d z=0, is proven  相似文献   

6.
The standard deviation of vertical two-point longitudinal velocity fluctuation differences is analyzed experimentally with eleven sets of turbulence measurements obtained at the NASA 150-m ground-winds tower site at Cape Kennedy, Florida. It is concluded that /u *0 is proportional to (fz/u *0)0.22, where the coefficient of proportionality is a function of fz/u *0 and u *0/fL 0. The quantities f and L0 denote the Coriolis parameter and the surface Monin-Obukhov stability length, respectively; u *0 is the surface friction velocity; z is the vertical distance between the two points over which the velocity difference is calculated; and zz is the mean height of the mid-point of the interval z above natural grade. The results of the analysis are valid for 20<-u *0/fL 0<2000.  相似文献   

7.
Summary A dispersion model is proposed to predict the continuous vertical variation of the dispersion parameters y and z in case of hot pollutant release to the atmosphere. In such a case, the plume rises far above the ground and is subject to varying levels of turbulence. The framework in this paper can be divided into three approaches: (1) determination of the eddy diffusivitiesK y (z, y ) andK z (z, z ) as functions of height above ground and plume dimensions, (2) determination of both the plume rise and its vertical velocity using a modified version of Brigg's formula, and (3) numerical solution of actual problems with buoyant plumes at each time step. The model results have been applied to a case of pollutant release from fire destruction of a chemical storehouse roof.With 15 Figures  相似文献   

8.
Summary The standard equations for the theory of atmospheric tides are solved here by an integral representation on the continuous spectrum of free oscillations. The model profile of back-ground temperature is that of the U.S. Standard Atmosphere in the lower and middle atmosphere, and in the lower thermosphere, above which an isothermal top extends to arbitrarily great heights. The top is warm enough to bring both the Lamb and the Pekeris modes into the continuous spectrum.Computations are made for semidiurnal lunar tidal pressure at sea level at the equator, and the contributions are partitioned according to vertical as well as horizontal structure. Almost all the response is taken up by the Lamb and Pekeris modes of the slowest westward-propagating gravity wave. At sea level, the Lamb-mode response is direct and is relatively insensitive to details of the temperature profile. The Pekeris mode at sea level has an indirect response-in competition with the Lamb mode-and, as has been known since the time of its discovery, it is quite sensitive to the temperature profile, in particular to stratopause temperature. In the standard atmosphere the Lamb mode contributes about +0.078 mb to tidal surface pressure at the equator and the Pekeris mode about –0.048 mb.The aim of this investigation is to illustrate some consequences of representing the tide in terms of the structures of free oscillations. To simplify that task as much as possible, all modifying influences were omitted, such as background wind and ocean or earth tide. Perhaps the main defect of this paper's implementation of the free-oscillation spectrum is that, in contrast to the conventional expansion in the structures of forced oscillations, it does not include dissipation, either implicity or explicity, and thus does not satisfy causality. Dissipation could be added implicity by means of an impedance condition, for example, which would cause up-going energy flux to exceed downgoing flux at the base of the isothermal top layer. To achieve complete causality, however, the dissipation must be modeled explicity. Nevertheless, since the Lamb and Pekeris modes are strongly trapped in the lower and middle atmosphere, where dissipation is rather weak (except possibly in the surface boundary layer), more realistic modeling is not likely to change the broad features of the present results.Symbols a earth's mean radius; expansion coefficient in (5.3) - b recursion variable in (7.4); proximity to resonance in (9.2) - c sound speed in (2.2); specific heatc p in (2.2) - f Coriolis parameter 2sin in (2.2) - g standard surface gravity - h equivalent depth - i ; discretization index in (7.3) - j index for horizontal structure - k index for horizontal structure; upward unit vectork in (2.2) - m wave number in longitude - n spherical-harmonic degree; number of grid layers in a model layer - p tidal pressure perturbation; background pressurep 0 - q heating function (energy per mass per time) - r tidal state vector in (2.1) - s tidal entropy perturbation; background entropys 0 - t time - u tidal horizontal velocityu - w tidal vertical component of velocity - x excitation vector defined in (2.3); vertical coordinate lnp */p 0 [except in (3.8), where it is lnp /p 0] - y vertical-structure function in (7.1) - z geopotential height - A constant defined in (6.2) - C spherical-harmonic expansion coefficient in (3.6) - D vertical cross section defined in (5.6) and (5.9) - E eigenstate vector - F vertical-structure function for eigenstate pressure in (3.2) [re-defined with WKB scaling in (7.2)] - G vertical-structure function for eigenstate vertical velocity in (3.2) [re-defined with WKB scaling in (7.2)] - H pressure-scale height - I mode intensity defined in (8.1) - K quadratic form defined in (4.4) - L quadratic form defined in (4.4); horizontal-structure magnification factor defined in (5.11) - M vertical-structure magnification factor defined in (4.6) - P eigenstate pressure in (3.2); tidal pressure in (6.2) - R tidal state vector in (5.1) - S eigenstate entropy in (3.2); spherical surface area, in differential dS - T background molecular-scale (NOAA, 1976) absolute temperatureT 0 - U eigenstate horizontal velocityU in (3.2); coefficient in (7.3) - V horizontal-structure functionV for eigenstate horizontal velocity in (3.2); recursion variable in (7.3) - W eigenstate vertical velocity in (3.2) - X excitation vector in (5.1) - Y surface spherical harmonic in (3.7) - Z Hough function defined in (3.6) - +dH/dz - (1––)/2 - Kronecker delta; Dirac delta; correction operator in (7.6) - equilibrium tide elevation - (square-root of Hough-function eigenvalue) - ratio of specific gas constant to specific heat for air=2/7 - longitude - - - background density 0 - eigenstate frequency in (3.1) - proxy for heating functionq =c P/t - latitude - tide frequency - operator for the limitz - horizontal-structure function for eigenstate pressure in (3.2) - Hough function defined in (6.2) - earth's rotation speed - horizontal gradient operator - ()0 background variable - ()* surface value of background variable - () value at base of isothermal top layer - Õ state vector with zerow-component - , energy product defined in (2.4) - | | energy norm - ()* complex conjugate With 10 Figures  相似文献   

9.
In a recent paper, the author introduced a new viscous boundary layer, called the mesolayer, in turbulent shear flow. Its importance stems from its location between the inner and outer regions which are controlled by the law of the wall and Reynolds number similarity, respectively. This intrusion prevents the classical overlap assumption which appears to be fundamental in the derivation of the classical logarithmic behavior. The mesolayer has a thickness proportional to Taylor's microscale . This, and the analogy between the energy equation for the spectrum function of isotropic turbulence and the momentum equation for shear flow, suggest the existence of a similar region in wavenumber space with wavenumber k ~ -1. This mesoregion separates the inner region k ~ k s(where k s-1 and is the Kolmogorov length) and the outer region k k e(where k e -1 and l is the energy-containing eddy size) and again invalidates the overlap assumption which appears to be fundamental in the derivation of the classical k -5/3-behavior of the inertial subrange.Incorporation of the mesoregion into the argument leads to a new theory with k -5/3-behavior in two regions (-1 k k s) and (k e k -1) although with two different coefficients of proportionality (Kolmogorov constants). This leads to a wandering of the spectrum curve about the classical k -5/3 line similar to a wandering in turbulent shear flow about the logarithmic curve. This is clearly indicated by the data for the variation of the Kolmogorov constant.Other data support the new theory. In particular, the location of the point k mwhere the curve of the nonlinear energy-transfer function goes through zero shows agreement with the theory, i.e., k m-1.  相似文献   

10.
Summary A fixed platform (Fig.3), installed 100 m from the shoreline in 3 m water depth, was instrumented with velocity, temperature and wave-height sensors. 132 data (10 minutes averages) were analysed to calculate the wind stress; from these, 99 data were used to investigate the vertical distribution of the wind stress; all data are presented with Table 1.It was postulated that the total stress, t being constant with height, is made up additively of two components, the wave-supporting stress, w , and the turbulent stress, c ; see Eq. 1. The vertical distribution of these two components is shown schematically in Fig. 1.The total stress, t , evaluated outside the zone of wave influence, is given in the classical way with Fig. 4. The wave-supporting stress, w (z), was evaluated from the data according to a relation proposed by Kitaigorodskii et al. (1984); it is given with Fig. 5. A height-dependency is clearly evident. The turbulent stress c (z), was evaluated with data of the velocity gradient; it is given with Fig. 6. A height-dependency is not evident.The field data from the lake of Geneva give evidence that the additive relation of Eq. 1 seems to be justified.With 6 Figures  相似文献   

11.
In the summer of 1988/89 flights were carried out in the Coorong coastal area of South Australia to investigate sea-breeze fronts. The flights yielded data sets of the structure of the fronts in the cross-frontal direction with a spatial resolution of approximately 3 m. The study is focused on the budgets of sensible and latent heat in the vicinity of the front and on frontogenesis/frontolysis processes which are closely related to budget considerations.The frontogenesis relationships and the budgets were established on a 2 km length scale by low-pass filtering of the space series. As the wind components were measured with high accuracy, all processes which determine frontogenesis could be evaluated and are displayed in x,z-cross-sections: these are the confluence, shear and diabatic effects, all of which play a role in q/x-, q/z-, /x- as well as /z-frontogenesis. A detailed analysis is given for two different states of frontal development. The presented results shed much light on the governing physical processes in the frontal region with strong emphasis on the effects of confluence-generated updrafts, on shear instabilities causing bulges and clefts in the frontal surface as well as producing the elevated frontal head, and on processes related to differential heating and moistening.  相似文献   

12.
From measured one-dimensional spectra of velocity and temperature variance, the universal functions of the Monin-Obukhov similarity theory are calculated for the range –2 z/L + 2. The calculations show good agreement with observations with the exception of a range –1 z/L 0 in which the function m , i.e., the nondimensional mean shear, is overestimated. This overestimation is shown to be caused by neglecting the spectral divergence of a vertical transport of turbulent kinetic energy. The integral of the spectral divergence over the entire wave number space is suggested to be negligibly small in comparison with production and dissipation of turbulent kinetic energy.Notation a,b,c contants (see Equations (–4)) - Ci constants i=u, v, w, (see Equation (5) - kme,kmT peak wave numbers of 3-d moel spectra of turbulent kinetic energy and of temperature variance, respectively - kmi peak wave numbers of 1-d spectra of velocity components i=u, v, w and of temperature fluctuations i= - ksb, kc characteristics wave numbers of energy-feeding by mechanical effects being modified by mean buoyancy, and of convective energy feeding, respectively - L Monin-Obukhov length - % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Gabeivayaaraaaaa!3C5B!\[{\rm{\bar T}}\] difference of mean temperature and mean potential temperature - T* Monin-Obukhov temperature scale - velocity of mean flow in positive x-direction - u* friction velocity - u, v, w components of velocity fluctuations - z height above ground - von Kármanán constant - temperature fluctuation - m nondimensional mean shear - H nondimensional mean temperature gradient - nondimensional rate of lolecular dissipation of turbulent kinetic energy - D nondimensional divergence of vertical transports of turbulent linetic energy  相似文献   

13.
Careful micrometeorological measurements on an empty parking lot allowed determination of the surface fluxes of sensible heatH and of momentum by applying profile equations derived from Monin-Obukhov similarity theory with two sets of the stability correction function for momentum m and sensible heat h . These fluxes were compared with reference values ofH independently determined by means of an eddy correlation technique. In general, better agreement was found betweenH values derived from profiles with the stability functions of Brutsaert (1992) and referenceH values, than when the Businger-Dyer functions were used to deriveH. The disagreement in the latter comparison was especially serious under strongly unstable conditions, with the value ofy=–z/L (wherez is the height andL is the Obukhov length) larger than 10. A closer look at the procedure for calculatingH from the profiles revealed that the large differences between theH values derived with these two different versions of the stability correction functions were caused by the small differences of the h values, and not by the larger differences of the m values. This result stems from the strong sensitivity of the resultingH values on the choice of h .  相似文献   

14.
A Field Study of the Mean Pressure About a Windbreak   总被引:3,自引:0,他引:3  
To provide additional field data for assessingwindbreak flow models, mean ground-level pressurehas been measured upstream and downstream from along porous fence (height H = 1.25 m, resistancecoefficient k r = 2.4). Measurements were madeduring periods of near-neutral stability and near-normallyincident flow, with the fence standing on bare soil(roughness length, z 0 0.8 cm;H/z 0 160), or within a plant canopy. The mean pressure field,measured far from the ends of the fence, was foundto be quite insensitive to mean wind direction( , zero for perpendicular flow), for| | less than about 25°.In the absence of a canopy, during each measurementperiod the minimum pressure occurred at the closestsampling location to leeward of the windbreak, thepressure-gradient in most cases beingmaximally-adverse in the immediate lee, and decayingwith increasing downwind distance (x). On one day ofmeasurements, however, the pressure gradient over2 x/H 6 (H = windbreak height) resembled theleeward plateau identified by Wang and Taklein their numerical studies. Perhaps thisoccasional feature was only due to instrumenterror. Nevertheless a plateau of sorts wasindicated in similar measurements by Judd andPrendergast (with H = 1.92 m, z 0 1.2 cm;H/z 0 160, k r 3). Therefore,existence of a leeward pressure plateau behind athin fence cannot be definitely ruled out.When the windbreak was placed in a canopy, minimumsurface pressure was displaced downwind. Thisagrees with the wind-tunnel study of Judd, Raupach and Finnigan,and is consistent with a simple simulation reported here.  相似文献   

15.
From measurements in the atmospheric surface layer over a paddy field, the Kolmogorov constants for CO2 and longitudinal wind velocity were obtained. In this study, the nondimensional dissipation rate nc = (1–16 v )-1/2 for CO2 variance and = (1–16 v )-1/4 v for turbulent energy were used, assuming the equality of the local production term and the local dissipation term, and neglecting the divergence flux term in the budget equation. The value of the constant for CO2 was consistent with recent determinations for temperature and humidity. The constant for longitudinal wind velocity showed good agreement with other recent observations.  相似文献   

16.
Frequency spectra of atmospheric turbulenceS (f) in the inertial subrange are considered in the free convection regime over the sea surface in a case of motionless instrument measurements (Eulerian frequency spectra). The frequency spectra formulaef * S (f)/ 2 =c (f */f)5/3 for wind velocity (=1–3), temperature (=t) and humidity (=e) fluctuations are derived on the basis of similarity theory and the –5/3 law. These relations also can be derived from a consideration of convective large-scale advection of small eddies. The frequency scalef * = (N 1 2/)1/2 (H/z 2)1/3 is the lower bound of the inertial subrange and it is of order 10–2 Hz.The spectra formulae are compared with direct measurements of atmospheric turbulence from the fixed research tower in the coastal zone of the Black Sea in calm weather. It is shown that these formulae are realized at least over two to three decades of the frequency range (approximately from 10–2 to 10 Hz) and values of the numerical coefficients are found. The derived formulae can be used for calculations of sensible and latent heat fluxes by measuring the high-frequency range of spectra at a fixed point at low wind speeds when the conventional inertial dissipation method is not applicable.  相似文献   

17.
Components of the radiation and energy balances were measured over a clear-cut area and a mature, mixed forest during the summer of 1981 at the Petawawa National Forestry Institute, Chalk River, Ontario. The work concentrated on the clear-cut site which supported a canopy layer composed primarily of bracken fern and logging remnants.Forty days of radiation data were collected at the clear-cut site. After the first four weeks of measurements (the green season), most of the ferns quickly died, and their foliage changed appearance from a green to brownish colour (the brown season). The daily mean reflection coefficient of solar radiation determined over the green season was 0.20 and decreased to 0.13 for the brown season. The corresponding value for the forest was 0.13, based on a limited amount of data. The clear-cut site received 11% and 21% less net radiation than the forest on a 24-hr and daylight-hours basis, respectively, as a consequence of the higher reflection coefficient and larger daytime longwave radiation emission.A reversing temperature difference measurement system (RTDMS), incorporating ten-junction thermopiles was employed at each site in order to determine Bowen ratios () via differential psychrometry. Both systems performed well, especially the RTDMS over the forest which was capable of resolving very small differences of temperature, typically less than 0.2 °C over a height of 3 m. The mean hourly Bowen ratio, calculated from values from 0800 to 1600 hr, varied from 0.2 to 1.0 for the forest and from 0.4 to 0.8 for the clear-cut site in the green season.A significant canopy heat storage component of the energy balance, Q S , was found at the clear-cut site. In the early morning, a portion of the available energy was used to heat the biomass materials and air within the canopy layer. The stored heat within the canopy was released later in the day, increasing the available energy total.The daily mean value of the Priestley-Taylor coefficient (Priestley and Taylor, 1972) for the green season at the clear-cut site was 1.14, and individual values tended to increase during wet surface conditions and decrease when the surface dried. The daylight mean value during dry canopy conditions at the forest was 1.05, and much higher values occurred when the canopy was wet. The enhancement of for the wet forest was a result of the evaporation of intercepted rain (which is not limited by stomatal resistance) and the concomitant transfer of sensible heat to the forest.  相似文献   

18.
TheConvectiveDiffusionObserved byRemoteSensors (CONDORS) field experiment conducted at the Boulder Atmospheric Observatory used innovative techniques to obtain three-dimensional mappings of plume concentration fields, /Q, of oil fog detected by lidar and chaff detected by Doppler radar. It included extensive meteorological measurements and, in 1983, tracer gases measured at a single sampling arc. Final results from ten hours of elevated and surface release data are summarized here. Many intercomparisons were made. Oil fog /Q measured 40m above the arc are mostly in good agreement withSF 6 values, except in a few instances with large spacial inhomogeneities over short distances. After a correction scheme was applied to compensate for the effect of its settling speed, chaff dy/Q agreed well with those of oil except in two cases of oil fog hot spots. Mass or frequency distribution vs. azimuth or elevation angle comparisons were made for chaff, oil, and wind, with mostly good agreements. Spacial standard deviations, y and z, of chaff and oil agree overall and are consistent at short range with velocity standard deviations vand w 0.6w* (the convective scale velocity), as measured atz>100m. Surface release y is enhanced up to 60% at smallx, consistent with the Prairie Grass measurements and with larger v and reduced wind speed measured near the surface. Decreased y at small dimensionless average times is also noted. Finally, convectively scaled dy, C y, were plotted versus dimensionlessx andz for oil, chaff, and corrected chaff for each 30–60 min period. Aggregated CONDORSC y fields compare well with laboratory tank and LES numerical simulations; surface-released oil fog compares expecially well with the tank experiments. However, large deviations from the norm occurred in individual averaging periods; these deviations correlated strongly with anomalies in measured distributions.On assignment to the US Environmental Protection Agency, Atmospheric Research and Exposure Assessment Laboratory, RTP, NC.  相似文献   

19.
Analytical solutions for the Ekman layer   总被引:1,自引:0,他引:1  
The PBL equation that governs the transition from the constant-stress surface layer to the geostrophic wind in a neutrally stratified atmosphere for which the eddy viscosityK(z) is assumed to vary smoothly from the surface-layer value U *z (0.4,U *=friction velocity,z=elevation) to the geostrophic asymptoteK GU *d forzd is solved through an expansion in fd/U *1 (f=Coriolis parameter). The resulting solution is separated into Ekman's constant-K solution an inner component that reduces to the classical logarithmic form forzd and isO() relative to the Ekman component forzd. The approximationKU *d is supported by the solution of Nee and Kovasznay's phenomenological transport equation forK(z), which yieldsKU *d exp(–z/d), where is an empirical constant for which observation implies, 1. The parametersA andB in Kazanskii and Monin's similarity relation forG/U * (G=geostrophic velocity) are determined as functions of . The predicted values ofG/U * and the turning angle are in agreement with the observed values for the Leipzig wind profile. The predicted value ofB based on the assumption of asymptotically constantK is 4.5, while that based on the Nee-Kovasznay model is 5.1; these compare with the observed value of 4.7 for the Leipzig profile. A thermal wind correction, an asymptotic solution for arbitraryK(z) and 1, and an exact (unrestricted ) solution forK(z)=U *d[1–exp(–z/d)] are developed in appendices.  相似文献   

20.
Nine profiles of the temperature structure parameter C T 2 and the standard deviation of vertical velocity fluctuations ( w) in the convective boundary layer (CBL) were obtained with a monostatic Doppler sodar during the second intensive field campaign of the First ISLSCP Field Experiment in 1987. The results were analyzed by using local similarity theory. Local similarity curves depend on four parameters: the height of the mixed layer (z i ), the depth of the interfacial layer (), and the temperature fluxes at the top of the mixed layer (Q i ) and the surface (Q o). Values of these parameters were inferred from sodar data by using the similarity curve for C T 2 and observations at three points in its profile. The effects of entrainment processes on the profiles of C T 2 and wnear the top of the CBL appeared to be described well by local similarity theory. Inferred estimates of surface temperature flux, however, were underestimated in comparison to fluxes measured by eddy correlation. The measured values of wappeared to be slightly smaller than estimates based on available parmeterizations. These discrepancies might have been caused by experimental error or, more likely, by the distortion of turbulence structure above the site by flow over the nonuniform terrain at the observation site.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号