首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
In-water and remote measurements of ocean color   总被引：2，自引：0，他引：2  
Spectral measurements of downwelling irradiance, E d(), above the surface, and of upwelling irradiance just below the surface, E u(), allow computation of spectral values of the diffuse reflectance R() = E u()/E d(); this yields full information about the true color and brightness of the ocean. Typical results are presented and interpreted for waters very different in turbidity and phytoplankton content. Conversely, the possibility of infering the water content from R() data at selected wavelengths is examined in terms of the respective number of equations and unknowns. The necessary use of assumptions and of empirical laws is emphasized.The magnitude of the useful signal emerging from the water, and the magnitude of the additional signals due to specular reflexion at the interface and to atmospheric scattering, are compared on the basis of spectroradiometric measurements performed within and above the sea, from different altitudes. These unwanted signals are dominant, causing a drastic change in the spectral composition of the light received by a remote sensor. The evaluation of the atmospheric effect must be very precise in order to recover the marine signal with a sufficient accuracy for a meaningful application of any kind of algorithm.This work was supported in part by the Centre National d'Exploitation des Océans (under contract CNEXO 77/1695) and in part by the Centre National de la Recherche Scientifique (RCP 247 & ERA 278).  相似文献   

2.
Zusammenfassung Bisher wurde der Umkehreffekt vonGötz nur mittels Photozelle oder Zählrohr untersucht. Die Anwendung des Multipliers als Meßinstrument bietet wegen dessen wesentlich größerer Empfindlichkeit viele Vorteile. So kann damit noch bei Zenitdistanzen der Sonnez>90° bei kleineren Wellenlängen als bisher und bei engerem Spalt des Spektrometers (bessere Monochromasie) gemessen werden.Resultate: 1. Die Verschiebung der Minima der Umkehrkurven,i n /i m =f(z), erfolgtgleichsinnig mit abnehmender Wellenlänge n nach kleineren Zenitdistanzen (entgegen den Angaben vonSchein undStoll, aber in Übereinstimmung mit der Theorie nachGötz undDobson). 2. Zum ersten Mal beziehen sich hier die Angaben aufeffektive Wellenlängen, die ihrerseits noch eine Funktion vonz sind. 3. Aus den Messungen bei den kleinsten, bisher noch nicht zugänglichen Wellenlängen folgt, daß deren Streuniveaus bei tiefen Sonnenständen oberhalb der Schichten größten Ozongehaltes liegen müssen.
Summary Until now the Umkehr-effect of Götz (inversion of intensity-ratio in zenith-scattered sunlight) has been investigated only by means of photocells or Geiger counters. The application of the multiplier as measuring instrument is of greater advantage because of its much higher sensitivity. Measurements can be made with it even with the sun's zenith distancesz>90° at shorter wave-lengths and with a narrower slit of the spectrometer (better monochromasy).Results: (1) The shift of the minima of the Umkehr-curvesi n /i m =f(z), erfolgtgleichsinnig mit abnehmender Wellenlänge n shows the same sense with decreasing wave-length n to shorter zenith distances (in opposition to the assertions ofSchein andStoll, but in agreement with the theory ofGötz andDobson). (2) For the first time the statements refer here to effective wave-lengths which in their turn are a function ofz. (3) From measurements of the shortest wave-lengths which until now have not yet been accessible, it can be concluded that at low sun heights the scattering level is situated above the layers with the greatest ozone content.
Résumé Jusqu'à présent le «Umkehreffect» deGötz, c. à d. l'effet d'inversion des relations d'intensité dans le rayonnement solaire zénithal diffus, n'a été étudié qu'au moyen de cellules photo-électriques ou de compteurs Geiger. En utilisant le multiplicateur comme instrument de mesure on obtient l'avantage d'une bien plus grande sensibilité. Ainsi il est possible de mesurer encore à des distances zénithales du soleil dez>90° par longueurs d'ondes plus courtes que jusqu'ici et avec une fente du spectromètre plus étroite, ce qui nous garantit une meilleure monochromasie.Résultats: 1) Le déplacement des minima des courbes d'inversioni n /i m =f(z), erfolgtgleichsinnig mit abnehmender Wellenlänge n s'effectue dans le même sens que la diminution des longuers d'ondes n , par distances zénithales plus petites (contrairement aux assertions deSchein et deStoll, mais en concordance avec la théorie selonGötz etDobson). — 2) Pour la première fois les données se rapportent ici à des longueurs d'ondes effectives qui, de leur côté, sont encore fonctions dez. — 3) Il ressort des mesures des plus courtes longueurs d'ondes, qui, jusqu'ici, restaient hors d'atteinte, que pour des hauteurs du soleil basses, leur niveau de fiffusion se trouve au-dessus des couches ayant le plus grand contenu d'ozone.

Mit 9 Textabbildungen  相似文献   

3.
STAR (System for Transfer of Atmospheric Radiation) was developed to calculate accurately and efficiently the irradiance, the actinic flux, and the radiance in the troposphere. Additionally a very efficient calculation scheme to computer photolysis frequencies for 21 different gases was evolved. STAR includes representative data bases for atmospheric constituents, especially aerosol particles. With this model package a sensitivity study of the influence of different parameter on photolysis frequencies in particular of O3 to Singlet D oxygen atoms, of NO2, and of HCHO was performed. The results show the quantitative effects of the influence of the solar zenith angle, the ozone concentration and vertical profile, the aerosol particles, the surface albedo, the temperature, the pressure, the concentration of NO2, and different types of clouds on the photolysis frequencies.Notation I A(, ) actinic flux - I H(, ) irradiance - L(, , , ) radiance - wavelength - azimuth angle - cosine of zenith angle - s cosine of solar zenith angle - optical depth - s scattering coefficient - c extinction coefficient - o single scattering albedo - p mix mixed phase function - g mix mixed asymmetry factor - J gas photolysis frequency  相似文献   

4.
A two-dimensional mesoscale soil-atmosphere model is used to simulate the triggering of atmospheric convection by horizontally varying soil water content. The variation is periodic with a wavelength between 4 and 40 km, which is considered a realistic scale for the variation of land surface characteristics. Three stages of convection can be clearly discerned: a short initial stage when convection sets in and where the size of the conective cells is determined by , a mature stage with well developed cells whose size is still determined by , and a decay/transformation stage, characterized by the formation of narrow regions of strong updrafts and wide regions of moderate downdrafts, independent of . Parameters relevant for the transition are given, and the importance of the feedback between soil and atmosphere is demonstrated. The dependence of convective parameters, e.g., height of the convective layer, vertical velocity and fluxes of heat and moisture on is investigated. The calculations of the mature stage are compared with the predictions of a linear model.  相似文献   

5.
The relation between the turbulence Reynolds numberR and a Reynolds numberz* based on the friction velocity and height from the ground is established using direct measurements of the r.m.s. longitudinal velocity and turbulent energy dissipation in the atmospheric surface layer. Measurements of the relative magnitude of components of the turbulent kinetic energy budget in the stability range 0 >z/L 0.4 indicate that local balance between production and dissipation is maintained. Approximate expressions, in terms of readily measured micrometeorological quantities, are proposed for the Taylor microscale and the Kolmogorov length scale .  相似文献   

6.
Equilibrium evaporation beneath a growing convective boundary layer   总被引：1，自引：1，他引：0  
Expressions for the equilibrium surface Bowen ratio ( s ) and equilibrium evaporation are derived for a growing convective boundary layer (CBL) in terms of the Bowen ratio at the top of the mixed layer i and the entrainment parameter A R . If AR is put equal to zero, the solution for s becomes-that previously obtained for the zero entrainment or closed box model. The Priestley-Taylor parameter is also calculated and plotted in terms ofA R and i . Realistic combinations of the atmospheric parameters give values of in the range 1.1 to 1.4.  相似文献   

7.
The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O2 of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O2 versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(Hdecomposition–HO 2 reaction), between these two reactions pathways. For (H)[22–0.5(HO 2 reaction)], alkoxy radical decomposition dominates over the reaction with O2 at room temperature and atmospheric pressure of air, while for (H)[25-0.5(HO 2 reaction)], the O2 reaction dominates over decomposition (where the units of H are in kcal mol–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.  相似文献   

8.
Summary During an expedition to the high Andes of Southern Peru in June–July 1977, measurements of direct solar radiation in four spectral bands (0.270–0.530–0.630–0.695–2.900 ) were conducted at six sites in elevations ranging from sea level to 5645 m. These measurements were evaluated in Langley plots to determine total optical depths () and irradiances at the top of the atmosphere. In addition, water vapor optical depths (wv) were calculated from the mean radiosounding over Lima during the expedition, and Rayleigh (ray) and ozone (oz) optical depths were obtained from published tabulations. Subtracting ray, oz, and wv from yielded estimates of aerosol optical depth aer. The components ray and oz decrease from the shorter towards the longer wavelength bands and from the lower towards the higher elevation sites; aer also decreases towards the higher elevations. Particularly pronounced is the decrease of aer and from the lowlands of the Pacific coast to the highlands of the interior, reflecting the effect of a persistent lower-tropospheric inversion and the contrast from the marine boundary layer to the clear atmosphere of the high Andes.With 4 Figures  相似文献   

9.
A Comparative Analysis of Transpiration and Bare Soil Evaporation   总被引：4，自引：0，他引：4  
Transpiration Ev and bare soil evaporation Eb processes are comparatively analysed assuming homogeneous and inhomogeneous areal distributions of volumetric soil moisture content . For a homogeneous areal distribution of we use a deterministic model, while for inhomogeneous distributions a statistical-deterministic diagnostic surface energy balance model is applied. The areal variations of are simulated by Monte-Carlo runs assuming normal distributions of .The numerical experiments are performed for loam. In the experiments we used different parameterizations for vegetation and bare soil surface resistances and strong atmospheric forcing. According to the results theEv()-Eb() differences are great, especially in dry conditions. In spite of this, the available energy flux curves of vegetation Av() and bare soil Ab() surfaces differ much less than the Ev() and Eb() curves. The results suggest that Ev is much more non-linearly related to environmental conditions than Eb. Both Ev and Eb depend on the distribution of , the wetness regime and the parameterization used. With the parameterizations, Eb showed greater variations than Ev. These results are valid when there are no advective effects or mesoscale circulation patterns and the stratification is unstable.  相似文献   

10.
Tropospheric photodissociation rate coefficients (J values) were calculated for NO2, O3, HNO2, CH2O, and CH3CHO using high spectral resolution (0.1 mm wavelength increments), and compared to the J values obtained with numerically degraded resolution (=1, 2, 4, 6, 8, and 10 nm, and several commonly used nonuniform grids). Depending on the molecule, substantial errors can be introduced by the larger increments. Thus for =10 nm, errors are less than 1% for NO2, less than 2% for HNO2, +6.5% to -16% for CH2O, -6.9% to +24% for CH3CHO, and -24% to +110% for O3. The errors for CH2O arise from the fine structure of its absorption spectrum, and are prevalently negative (underestimate of J). The errors for O3, and to a lesser extent for CH3CHO, arise mainly from under-resolving the overlap of the molecular action spectrum and the tropospheric actinic flux in the wavelength region of stratospheric ozone attenuation. The sign of those errors depends on whether the actinic flux is averaged onto the grid before or after the radiative transfer calculation. In all cases studied, grids with 2 nm produced errors no larger than 5%.  相似文献   

11.
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.  相似文献   

12.
Concurrent measurements of the surface energy balance components (net radiation, heat storage, and sensible and latent heat fluxes) were made in three communities (open water, Phragmites australis, Scirpus acutus) in a wetland in north-central Nebraska, U.S.A., during May-October, 1994. The Bowen ratio – energy balance method was used to calculate latent and sensible heat fluxes. This paper presents results from the open water area. The heat stored in water (G) was found to play a major role in the energy exchange over the water surface. During daytime, G consumed 45–60% of R n , the net radiation (seasonally averaged daytime G was about 127 W m–2). At night, G was a significant source of energy (seasonally averaged nighttime G was about -135 Wm–). The diurnal pattern of latent heat flux ( E) did not follow that of R n . On some days, E was near zero during midday periods with large R n . The diurnal variability in E seemed to be significantly affected by temperature inversions formed over the cool water surface. The daily evaporation rate (E) ranged from 2 to 8 mm during the measurement period, and was generally between 70 and 135% of the equilibrium rate.  相似文献   

13.
Drag and drag partition on rough surfaces   总被引：13，自引：0，他引：13  
An analytic treatment of drag and drag partition on rough surfaces is given. The aims are to provide simple predictive expressions for practical applications, and to rationalize existing laboratory and atmospheric data into a single framework. Using dimensional analysis and two physical hypotheses, theoretical predictions are developed for total stress (described by the square root of the canopy drag coefficient), stress partition (described by the ratio S/ of the stress s on the underlying ground surface to total stress ), zero-plane displacement and roughness length. The stress partition prediction is the simple equation S/= 1/(1+), where = CRCS the ratio of element and surface drag coefficients. This prediction agrees very well with data and is free of adjustable constants. Other predictions also agree well with a range of laboratory and atmospheric data.  相似文献   

14.
Functional forms of the universal similarity functions A, B (for wind components parallel and normal to the surface stress), and C (for potential temperature difference) are determined based on the generalized theory of the resistance laws for the Planetary Boundary Layer (PBL). The similarity-profile functions for the surface layer are matched with the velocity and temperature-defect profiles that are assumed to have shapes modified by certain powers of nondimensional height z/h, where h is the PBL height. The powers of the outer-layer profile functions are determined, so that the functions become negligible in the surface layer. To close the temperature defect law, an assumption that the temperature gradient across the top of the PBL is continuous with the stratification of the overlying atmosphere is used. The result of this assumption is that nondimensional momentum and temperature profiles in the PBL can be described in terms of four basic ratios: (1) roughness ratio = /h (2) scale-height ratio =|f|h/u*, (3) ambient stratification parameter =h/*, and (4) stability parameter =h/L, where L is the Monin-Obukhov length, z0 is the surface roughness, is the upper-air stratification, u * is the friction velocity, and * is the temperature scale at the surface. For stable conditions, the scale-height ratio can be related to the atmospheric stability and the upperair stratification, and the generalized similarity and Rossby number similarity theories become identical. Under appropriate boundary conditions, function A is explicitly dependent on the stability parameter , while B is a function of scale-height ratio , which in turn depends on the stability. Function C is shown to be dependent on the stability and the upper-air stratification, due to the closure assumption used for the temperature profile.The suggested functional forms are compared with other empirical approximations by several authors. The general framework used to determine the functional forms needs to be tested against good boundary-layer measurements.  相似文献   

15.
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.  相似文献   

16.
Summary Rainfall in West Africa is examined in relation to monthly mean equivalent potential temperature ( e )at the earth's surface. The study revealed that monthly mean equivalent potential temperature ( e ) and monthly rainfall (R) generally decreased northwards from the equator.A good relationship existed betweenR and e in the northern zone of West Africa (i.e., north of 7.5° N). No definite relationship existed in the southern zone. In the northern zone, the departure of e from its annual mean ( ) first became positive about a month before the onset of the rains. Positive departures from ) generally resulted in more than normal (or average) rainfall in this zone. In general, little or no rainfall occurred in West Africa whenever e was less than 320 K.
Zusammenfassung Der Niederschlag (MonatssummeR) in Westafrika wird in Zusammenhang mit der mittleren monatlichen Äquivalent-temperatur ( e ) an der Erdoberfläche untersucht. Es zeigte sich, daß die Monatswerte beider Elemente im allgemeinen vom Äquator nach Norden abnehmen.ZwischenR und e ergab sich für das nördliche Westafrika (nördlich von 7.5° N) eine gute, für die südliche Zone jedoch keine beweisbare Übereinstimmung. In der nördlichen Zone übertraf e das Jahresmittel erstmals etwa einen Monat vor Beginn der Regenzeit. Positive Abweichungen vom mittleren e hatten immer übernormalen Niederschlag in dieser Zone zur Folge. Dagegen gab es wenig oder keinen Niederschlag in Westafrika, wenn e unter 320 K lag.

With 7 Figures  相似文献   

17.
Résumé Après avoir rappelé quelques notions de calcul tensoriel et les équations de la dynamique atmosphérique en coordonnées généralisées, l'auteur établit l'équation fondamentale de la rotationnelle absolue de l'air. Diverses formes de cette équation sont mentionnées. La première de ces formes est celle d'un bilan qui permet de définir le transport et le taux de production de la rotationnelle absolue de l'air. De la forme particulière de ce transport, il résulte que ce bilan se réduit à un bilan dans un espace à deux dimensions. Deux exemples illustrent les formules générales, le premier, en coordonnées sphériques (, ,r), le second, en coordonnées (, , ) où est un scalaire quelconque (pression atmosphérique, température potentielle, ...). Une autre forme de l'équation fondamentale est celle qui donne le taux d'accroissement individuel d'une composante de la rotationnelle absolue. Cette forme conduit à l'équation d'Ertel. En terminant, l'auteur généralise l'invariant d'Ertel-Rossby.
Summary Having recalled some results of the tensor analysis and the equations of atmospheric motion in general coordinates, the writer establishes the fundamental equation of the absolute vorticity. Different forms of this equation are mentioned. One of these is the equation of balance from which it is possible to deduce the definition of transport and production of vorticity. It is shown that the equation of balance of the absolute vorticity may be reduced to an equation of balance in two-dimensional space. Two examples illustrate the general equations and formulas: the first example, in spherical coordinates (, ,r), the second one in the coordinates (, , ) where is an arbitrary single-valued scalar quantity (pressure, potentiel temperature, ...). Another form of the absolute vorticity equation expresses the individual change of an arbitrary component of the absolute vorticity. This form leads toErtel's equation. Finally, the writer generalizes theErtel-Rossby invariant.
Zusammenfassung Nach einer Rekapitulation einiger Resultate der Tensorrechnung und der atmosphärischen Bewegungsgleichungen in generalisierten Koordinaten wird die absolute Wirbelgleichung aufgestellt und es werden verschiedene Formen dieser Gleichung erwähnt. Eine derselben ist die Bilanzformel, aus der die Definition des Wirbeltransports und der Anteil der Wirbelbildung abgeleitet werden können. Es wird gezeigt, daß die Bilanzgleichung des absoluten Wirbels auf eine Bilanzgleichung im zweidimensionalen Raum reduziert werden kann. Die allgemeinen Formeln werden durch zwei Beispiele erläutert: das erste in Kugelkoordinaten (, ,r), das zweite in den Koordinaten (, , ), wo eine beliebige skalare Größe (luftdruck, potentielle Temperatur usw.) darstellt. Eine andere Form der absoluten Wirbelgleichung gibt die individuelle Zunahme einer beliebigen Komponente der absoluten Wirbelstärke wieder; diese Form führt zurErtelschen Gleichung. Zum Schlusse wird dieErtel-Rossbysche Invariante verallgemeinert.
  相似文献   

18.
A photochemical scheme which includes a detailed treatment of multiple scattering up to solar zenith angles of 96° (developed for use in a GCM) has been used to study partitioning within chemical families. Attention is drawn to the different zenith angle dependence of diffuse radiation for the two spectral regions <310 nm and >310 nm. The effect that this has on the so-called 40 km ozone problem is discussed. The importance of correctly including multiple scattering for polar ozone studies is emphasised.  相似文献   

19.
The sensible heat flux density C and the latent heat flux density E are coupled in the case of a multi-layer model of vegetation. Therefore two linearly independent combinations of C and E, the enthalpy flux density H and the saturation heat flux density J, are introduced. Two electrical analogues, for H and J, are designed. They are equivalent to the resistance scheme for C and E, but uncoupled. Penman's formulas for C and E, which are applicable only to single-layer models, can be expressed equivalently in terms of H and J. This version of Penman's formulas can be extended easily to multi-layer canopies.  相似文献   

20.
Energy partitioning and evaporation were measured over three wetland surfaces in a subarctic coastal marsh during pre-growing and growing periods. These surfaces included an alder/willow woodland, a sedge marsh and a raised backshore sedge meadow. A combination model analysis was used to assess the relative importance of surface resistance and meteorological conditions on the magnitude of the Bowen ratio, , during the growing period.Overall, the three surfaces experienced important site-to-site and seasonal differences in and evaporation, Q E. During the non-foliated period, Q E was largest and was smallest for the open water marsh, while the dry backshore site experienced the smallest Q E and largest . The non-foliated woodland assumed intermediate values of and Q E. After the vegetation covers were established, the woodland assumed the smallest and largest Q E flux. It was also found that at the marsh site increased with the presence of a vegetation cover.Wind direction was always an important factor in determining Q E and at all sites. was substantially larger and Q E was smaller for onshore winds (i.e., originating from James Bay) than for offshore winds. The combination model analysis showed that canopy resistance at all sites was largest during warm offshore winds, which were associated with large saturation deficits. However, the effect of increased canopy resistance on during offshore winds was offset by a large climatological resistance, resulting in small values and large Q E. When winds originated from James Bay, canopy resistance was smaller than for offshore winds, but the climatological resistance also was much smaller, resulting in larger and small Q E. The results have important implications for changes in land cover and climate on the regional water balance.  相似文献   

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

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