首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 976 毫秒
1.
Quasi arithmetic and Archimedean functionals are used to build new classes of spectral densities for processes defined on any d-dimensional lattice \mathbbZd{\mathbb{Z}^d} and random fields defined on the d-dimensional Euclidean space \mathbbRd{\mathbb{R}^d}, given simple margins. We discuss the mathematical features of the proposed constructions, and show rigorously as well as through examples, that these new classes of spectra generalize celebrated classes introduced in the literature. Additionally, we obtain permissible spectral densities as linear combinations of quasi arithmetic or Archimedean functionals, whose associated correlation functions may attain negative values or oscillate between positive and negative ones. We finally show that these new classes of spectral densities can be used for nonseparable processes that are not necessarily diagonally symmetric.  相似文献   

2.
An extension of the multidimensional Born inversion technique for acoustic waves is described. In earlier work, a perturbation in reference sound velocity was determined by assuming that the reference velocity was constant. In this extension, we allow the reference velocity to be a function of the depth variable z. The output of this method is a high-frequency bandlimited reflectivity function of the subsurface. The reflectivity function is an array of bandlimited singular functions scaled by the normal reflection strength. Each singular function is a Dirac delta function of a scalar argument which measures distance normal to a reflecting interface. Thus, the reflectivity function is an indicator map of subsurface reflectors equivalent to the map produced by migration. In addition to the assumption of small perturbation, the method requires that the reflection data reside in the high frequency regime in a well-defined sense. The method is based on the derivation of an integral equation for the perturbation in sound velocity from a known reference velocity. When the reference velocity is constant, the integral equation admits an analytic solution as a multifold integral of the reflection data. Further high frequency asymptotic analysis simplifies this integral considerably and leads to an extremely efficient numerical algorithm for computing the reflectivity function. The development of a computer code to implement this constant-reference-velocity solution is published elsewhere. For a reference velocity c(z) we can no longer invert the integral equation exactly. However, we can write down an asymptotic high-frequency approximation for the kernal of the integral equation and an asymptotic solution for the perturbation. The computer implementation of this result is designed along the same lines as the code for constant background velocity. In tests the total processing time for this algorithm with depth-dependent background velocity is usually considerably less than that required by a standard Kirchhoff migration algorithm. The method is implemented as a migration technique and compared with alternative migration algorithms on the flanks of the salt dome.  相似文献   

3.
The modified tracer gas technique is used to determine the reaeration coefficient in six different water bodies of the Itajaí River basin, three with rural land use and three in urban areas. Propane was used as the tracer gas and Rhodamine WT as the conservative tracer, providing information on dilution, mixing and dispersion. Liquefied petroleum gas was used instead of high purity propane, aimed at reducing the costs associated with the field trials. Reaeration‐rate coefficients observed in the field ranged from 25.8 to 367.7 d?1. Two data sets could be observed where smaller streams had substantially larger coefficients of between 133.1 and 367.7 d?1, while the larger streams had values ??ranging from 25.8 to 54.5 d?1. Five empirical equations were evaluated by comparing the values ??obtained in the field. The equations proposed by Tsivolgou and Wallace and Tsivoglou and Neal showed greater adherence to the values ??determined in the tests. Reaeration‐rate coefficients obtained in the field were correlated with the hydrodynamic characteristics of the watercourses, thus establishing a mathematical function through which to obtain estimates for future evaluations. The R2 value obtained using this equation was 0.959, indicating a high correlation between the calculated values ??and those estimated in the field. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

4.
Short periodic oscillations with the periods from 10 up to 110 days of the hemispheric components of effective atmospheric angular momentum (EAAM) excitation function and their correlation with polar motion excitation function have been analyzed. The EAAM data of the Japan Meteorological Agency (JMA) computed for the two hemispheres and the very long baseline interferometry (VLBI) polar motion NGS 92 R01 data (NGS 1992), determined by the National Geodetic Survey were applied. The distinct oscillations with periods of about 28, 35–55 and 60–80 days were detected in the y-component of both polar motion excitation function and northern EAAM excitation functions containing wind and pressure, with and without inverted barometric correction terms. The y-component of the polar motion excitation function is significanly correlated (correlation coefficient equal to 0.55-0.75) with the y-components of the northern EAAM excitation functions mentioned above, which are mostly induced by the atmospheric circulation over lands. No meaningful correlation between polar motion excitation function and the southern EAAM excitation functions was found. The x-components of the EAAM and polar motion excitation functions are not significantly correlated. The strong short periodic variation of the length of day (LOD) and y in the early 1988 seems to be caused by the above-mentioned 35–55 days oscillations of the northern hemisphere atmosphere. This variation can be related to the rapid passing from the El Niño to the La Niña phenomenon or from the minimum to the maximum in the Southern Oscillation Index in 1987-1989.  相似文献   

5.
Isotope signatures in precipitation from the Global Network for Isotopes in Precipitation around the Mediterranean basin and literature data are compared with isotopic data from a large karstic aquifer in southeast Spain to explain the origin and type of the precipitation events dominating recharge. Analysis of the deuterium excess d at the scale of the Mediterranean basin and at the regional scale allows us to understand the isotopic context of the study area: Campo de Dalias and the Sierra de Gador (Almería province). The origin of precipitation can be determined from its d value. The d value changes as a function of the initial evaporation condition. It depends on the relative humidity and temperature during the evaporation producing the water vapour of the clouds. The water vapour, which dominates the study area, is generated in two areas: the Atlantic Ocean (d = 10‰) and the western Mediterranean basin (d = 15‰). With increasing precipitation volume, the western Mediterranean character dominates. These heavier storms contribute mainly to recharge, as illustrated by the d value of 13·6‰ in deep groundwater of the Campo de Dalias. Weighted d values increase with the volume of precipitation, giving a significant relationship for the southern and eastern coasts of the Iberian Peninsula. This selectivity of d to monthly precipitation was used to estimate the return period of precipitation leading to aquifer recharge at 0·9–4·9 years. Moderate rainfall, which occurs more frequently, still represents ～60–90% of the total precipitation. One of the challenges to meet ever‐growing water demands is to increase recharge from moderate events yielding intermediate quantities per event, but forming the bulk of the annual precipitation. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

6.
In this study, closed form equations as functions of the isolator, bridge and ground motion properties are formulated to calculate the optimum characteristic strength, Qd and post-elastic stiffness, kd, of the isolator to minimize the maximum isolator displacement (MID) and force (MIF) for seismic isolated bridges (SIBs). For this purpose, first, sensitivity analyses are conducted to identify the bridge, isolator and ground motion parameters that affect the optimum values of Qd and kd. Next, for the identified parameters, nonlinear time history analyses of typical SIBs are conducted to determine the optimum values of Qd and kd for a wide range of values of the parameters. Next, nonlinear regression analyses of the available data are conducted to obtain closed form equations for the optimum values of Qd and kd, to minimize the MID and MIF. The equations are then simplified for various site soil conditions. It is observed that the optimum Qd and kd are highly dependent on the site soil condition. Furthermore, the optimum Qd is found to be a linear function of the peak ground acceleration.  相似文献   

7.
The Canadian Land Surface Scheme (CLASS) was modified to correct an underestimation of the winter albedo in evergreen needleleaf forests. Default values for the visible and near‐infrared albedo of a canopy with intercepted snow, αVIS,cs and αNIR,cs, respectively, were too small, and the fraction of the canopy covered with snow, fsnow, increased too slowly with interception, producing a damped albedo response. A new model for fsnow is based on zI*, the effective depth of newly intercepted snow required to increase the canopy albedo to its maximum, which corresponds in the model with fsnow = 1. Snow unloading rates were extracted from visual assessments of photographs and modelled based on relationships with meteorological variables, replacing the time‐based method employed in CLASS. These parameterizations were tested in CLASS version 3.6 at boreal black spruce and jack pine forests in Saskatchewan, Canada, a subalpine Norway spruce and silver fir forest at Alptal, Switzerland, and a boreal maritime forest at Hitsujigaoka, Japan. Model configurations were assessed based on the index of agreement, d, relating simulated and observed daily albedo. The new model employs αVIS,cs = 0.27, αNIR,cs = 0.38 and zI* = 3 cm. The best single‐variable snow unloading algorithm, determined by the average cross‐site d, was based on wind speed. Two model configurations employing ensemble averages of the unloading rate as a function of total incoming radiation and wind speed, and air temperature and wind speed, respectively, produced larger minimum cross‐site d values but a smaller average. The default configuration of CLASS 3.6 produced a cross‐site average d from October to April of 0.58. The best model employing a single parameter (wind speed at the canopy top) for modelling the unloading rate produced an average d of 0.86, while the two‐parameter ensemble‐average unloading models produced a minimum d of 0.81 and an average d of 0.84. © 2015 Her Majesty the Queen in Right of Canada. Hydrological Processes published by John Wiley & Sons, Ltd.  相似文献   

8.
Land surface soil moisture (SSM) is an important variable for hydrological, ecological, and meteorological applications. A multi‐linear model has recently been proposed to determine the SSM content from the combined diurnal evolution of both land surface temperature (LST) and net surface shortwave radiation (NSSR) with the parameters TN (the LST mid‐morning rising rate divided by the NSSR rising rate during the same period) and td (the time of daily maximum temperature). However, in addition to the problem that all the coefficients of the multi‐linear model depend on the atmospheric conditions, the model also suffers from the problems of the nonlinearity of TN as a function of the SSM content and the uncertainty of determining the td from the diurnal evolution of the LST. To address these problems, a modified multi‐linear model was developed using the logarithm of TN and normalizing td by the mid‐morning temperature difference instead of using the TN and td. Except for the constant term, the coefficients of all other variables in the modified multi‐linear model proved to be independent of the atmospheric conditions. Using the relevant simulation data, results from the modified multi‐linear model show that the SSM content can be determined with a root mean square error (RMSE) of 0.030m3/m3, provided that the constant term is known or estimated day to day. The validation of the model was conducted using the field measurements at the Langfang site in 2008 in China. A higher correlation is achieved (coefficient of determination: R2 = 0.624, RMSE = 0.107m3/m3) between the measured SSM content and the SSM content estimated using the modified multi‐linear model with the coefficients determined from the simulation data. Another experiment is also conducted to estimate the SSM content using the modified model with the constant term calibrated each day by one‐spot measurements at the site. The estimation result has a relatively larger error (RMSE = 0.125m3/m3). Additionally, the uncertainty of the determination of the coefficients is analysed using the field measurements, and the results indicate that the SSM content obtained using the modified model accurately characterizes the surface soil moisture condition. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

9.
It has generally been assumed that diffusive sediment transport on soil‐mantled hillslopes is linearly dependent on hillslope gradient. Fieldwork was done near Santa Barbara, California, to develop a sediment transport equation for bioturbation by the pocket gopher (Thomomys bottae) and to determine whether it supports linear diffusion. The route taken by the sediment is divided into two parts, a subsurface path followed by a surface path. The first is the transport of soil through the burrow to the burrow opening. The second is the discharge of sediment from the burrow opening onto the hillslope surface. The total volumetric sediment flux, as a function of hillslope gradient, is found to be: qs (cm3 cm−1 a−1) = 176(dz/dx)3 − 189(dz/dx)2 + 68(dz/dx) + 34(dz/dx)0·4. This result does not support the use of linear diffusion for hillslopes where gopher bioturbation is the dominant mode of sediment transport. A one‐dimensional hillslope evolution program was used to evolve hillslope profiles according to non‐linear and linear diffusion and to compare them to a typical hillslope. The non‐linear case more closely resembles the actual profile with a convex cap at the divide leading into a straight midslope section. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

10.
Spatial and temporal variations of the isotopic composition of precipitation over Thailand were investigated. The local meteoric water line for Thailand deviates slightly from the global meteoric water line, with lower slopes (7.62 ± 0.07, 7.59 ± 0.08) and intercepts (6.42 ± 0.39, 6.22 ± 0.42) using ordinary and precipitation weighted methods. Differences in spatial and temporal δ18O distributions between the tropical monsoon and tropical savanna climate zones were found due to differing moisture source contributions and seasonal precipitation patterns. The temporal data reveals that the northeast monsoon rains originate from isotopically-enriched local moisture with isotope values of −9.36 to −0.09‰ (mean − 3.73 ± 0.42‰), whereas the southwest monsoon clouds had a more significant rainout effect from Rayleigh distillation, with isotope values of −9.56 to −1.78‰ (mean − 5.40 ± 0.38‰). The precipitation amount at each site was negatively correlated with δ18O (−0.24 to −3.20‰ per 100 mm, R2 = 0.1–0.9). Furthermore, δ18O was negatively correlated with geography (latitude, altitude) for the southwest monsoon periods, as expected based on other observed correlations. However, an inverse correlation was seen in the northeast monsoon due to differing moisture transportation as part of the continental effect. The correlation coefficient (R) was higher in the southwest monsoon (−0.84 for latitude effect, −0.64 for altitude effect) than the northeast monsoon (0.67 for latitude effect, 0.35 for altitude effect). The spatial pattern of isotopic composition reflects the southwest monsoon more clearly than the northeast monsoon, but the two monsoons also have a cancelling impact on orographic patterns. An agreement of the δ18O and deuterium excess (d-excess) was a negative correlation and found to reflect precipitation sources and re-evaporation processes. The d-excess was slightly higher for the northeast monsoon, bringing moisture from the Pacific Ocean and travelling across the continent before reaching the observed stations. By contrast, the d-excess was relatively lower for the Indian Ocean's moisture in the southwest monsoon.  相似文献   

11.
In this study, the removal of zinc(II) ion from an aqueous solution by pistachio shells (PS) is investigated. The dynamic behavior of the adsorption is examined on the effects of pH, adsorbent dosage, and contact time. The adsorption rates are determined quantitatively and simulated by the Lagergren first order, pseudo‐second order, Elovich, and intra‐particle diffusion kinetic models. The adsorption kinetic models are also tested for validity. The thermodynamic parameters, which are also deduced from adsorption experiments, are very useful in elucidating the nature of adsorption. The experimental results reveal that the optimum pH value and the contact time for the adsorption of Zn2+ onto PS are found as 6 and 10 min, respectively. According to these parameters, adsorption process follows the pseudo‐second order kinetic model with high correlation coefficients (R2 = 0.999). The obtained results demonstrate that PS is a reasonably effective adsorbent for the removal of Zn2+ from aqueous leachate of hazardous waste.  相似文献   

12.
The dispersion of inertial particles continuously emitted from a point source is analytically investigated in the limit of small but finite inertia. Our focus is on the evolution equation of the particle joint probability density function p(x,?v,?t), x and v being the particle position and velocity, respectively. For arbitrary inertia, position and velocity variables are coupled, with the result that p(x,?v,?t) can be determined by solving a partial differential equation in a 2d-dimensional space, d being the physical-space dimensionality. For small (but nevertheless finite) inertia, (x,?v)-variables decouple and the determination of p(x,?v,?t) is reduced to solve a system of two standard forced advection–diffusion equations in the space variables x. The latter equations are derived here from first principles, i.e., from the well-known Lagrangian evolution equations for position and particle velocity.  相似文献   

13.
Average steady source flow in heterogeneous porous formations is modelled by regarding the hydraulic conductivity K(x) as a stationary random space function (RSF). As a consequence, the flow variables become RSFs as well, and we are interested into calculating their moments. This problem has been intensively studied in the case of a Neumann type boundary condition at the source. However, there are many applications (such as well-type flows) for which the required boundary condition is that of Dirichlet. In order to fulfill such a requirement the strength of the source must be proportional to K(x), and therefore the source itself results a RSF. To solve flows driven by sources whose strength is spatially variable, we have used a perturbation procedure similar to that developed by Indelman and Abramovich (Water Resour Res 30:3385–3393, 1994) to analyze flows generated by sources of deterministic strength. Due to the linearity of the mathematical problem, we have focused on the explicit derivation of the mean head distribution G d (x) generated by a unit pulse. Such a distribution represents the fundamental solution to the average flow equations, and it is termed as mean Green function. The function G d (x) is derived here at the second order of approximation in the variance σ2 of the fluctuation (where K A is the mean value of K(x)), for arbitrary correlation function ρ(x), and any dimensionality d of the flow domain. We represent G d (x) as product between the homogeneous Green function G d (0)(x) valid in a domain with constant K A , and a distortion term Ψ d (x) = 1 + σ2ψ d (x) which modifies G d (0)(x) to account for the medium heterogeneity. In the case of isotropic formations ψ d (x) is expressed via one quadrature. This quadrature can be analytically calculated after adopting specific (e.g.. exponential and Gaussian) shape for ρ(x). These general results are subsequently used to investigate flow toward a partially-penetrating well in a semi-infinite domain. Indeed, we construct a σ2-order approximation to the mean as well as variance of the head by replacing the well with a singular segment. It is shown how the well-length combined with the medium heterogeneity affects the head distribution. We have introduced the concept of equivalent conductivity K eq(r,z). The main result is the relationship where the characteristic function ψ(w)(r,z) adjusts the homogeneous conductivity K A to account for the impact of the heterogeneity. In this way, a procedure can be developed to identify the aquifer hydraulic properties by means of field-scale head measurements. Finally, in the case of a fully penetrating well we have expressed the equivalent conductivity in analytical form, and we have shown that (being the effective conductivity for mean uniform flow), in agreement with the numerical simulations of Firmani et al. (Water Resour Res 42:W03422, 2006).  相似文献   

14.
A new method is presented for the determination of the parameters of the pulse contained in pulse‐like records. The Mavroeidis and Papageorgiou (M&P) wavelet is used for the mathematical representation of the pulse, but the proposed methodology could be modified to cover other types of wavelets as well. First, the period of the pulse is determined from the peak of the Sd × Sv product spectrum, a new concept that is introduced herein and is defined as the product of the velocity and the displacement response spectra. The remaining parameters of the M&P wavelet are derived from the targeted response spectrum of the ground motion applying a new relationship that is established between the cumulative absolute displacement (CAD) of a wavelet and its peak spectral amplitude. The method follows a well‐defined procedure that can be easily implemented in a computer code for the automatic determination of the pulse parameters of a given ground motion. As an application, the pulses contained in 91 NGA records that have been classified as pulse‐like by Baker are determined. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

15.
青藏高原地区分布的湖泊数量众多、面积较大、分布范围广泛.受制于恶劣的自然条件,对该地区湖泊的光学吸收特性以及光合有效辐射衰减系数(Kd(PAR))的研究鲜有成果.本文依据2014和2015年间采集的13个典型高海拔、湖泊面积较大的湖库的现场实测数据和实验室测定数据,分析了采样湖库各个采样点的Kd(PAR)特征以及有色可溶性有机物(CDOM)、藻类颗粒物吸收及非色素颗粒物吸收特性,计算并分析了Kd(PAR)与透明度以及光学活性物质的关系.研究结果表明:青藏高原地区各湖库平均各项颗粒物吸收系数均较低,总颗粒物吸收系数在400~700 nm波段内不超过0.14 m~(-1)、CDOM吸收系数在355 nm波长处最高,为1.23 m~(-1)、最低接近于0、藻类颗粒物吸收特性不明显;实验数据完整的巴木错、格仁错和班公错的主导吸收组分各异,其中巴木错为CDOM吸收主导,格仁错与班公错为非色素颗粒物吸收主导;青藏高原采样湖泊总体Kd(PAR)平均值较小,仅为0.26 m~(-1),样点最大值出现在可鲁克湖(1.17 m~(-1)),最小值出现在普莫雍错(0.10 m~(-1));在采样湖泊中Kd(PAR)与透明度呈显著相关;Kd(PAR)与CDOM的相关性最强,叶绿素a浓度次之,与总悬浮颗粒物浓度的相关性最不显著.  相似文献   

16.
Summary A method has been derived for computing the gravity effect of a segment of an infinite homogeneous elliptical cylinder. The initial data was represented by formulae expressing the components of the field of gravity of a homogeneous two-dimensional body by means of line integrals. The method is based on the integration of theln R function over the boundary of the cross-section of the attracting body, R being the distance from a fixed point in which the gravity effect is determined. The problem was solved in confocal co-ordinates.  相似文献   

17.
Permeability, k, and porosity, φ, were measured in mixtures of fine, medium, and coarse sand, where the volume fraction of each of the three components was systematically varied. The k was modeled well by the Kozeny‐Carman equation for three‐component mixtures by using a representative grain size parameter, d, computed by averaging the grain diameters of components recursively, with averaging methods based on whether finer components exist in sufficient volume to fill the pores within coarser components. The φ was modeled well by using linear interpolation with piecewise‐planar models. We explored the use of differing numbers of piecewise‐planar elements in the model, and illustrate the trade‐off between the increased accuracy and the increased data requirements that both come from adding more elements. The k model is a function of both d and φ, but more sensitive to d. The k model gave results consistent with measured values when computed using either measured φ values, or values from any of the φ models.  相似文献   

18.
Summary In adjusting a bilaterally connected polygonal traverse, whether to positionally accurate or inaccurate points of departure, the mean errors are the problem of primary importance. Hitherto, the mean errors of lengths and angles were chosen; the procedure was such that the weights of both quantities were determined from the mean errors chosen, and this fixed the ratio of the mean errors. This ratio did not change as a result of the adjustment, but the absolute values of the mean errors did. Provided the adjustment was carried out on a polygonal traverse with fixed points of departure, this change did not matter. In the case of ellipses of errors in the points of departure, this change is not permissible, because it would include the change of the semi-axes of the ellipses of errors and, therefore, also of the positional rigidity of the points of departure. The contribution of this paper is in the exact method by which it is possible to compute a coefficientc, pertaining to the mean angular error selected (in the case of positionally inaccurate pointsd), which determines the mean errorm=±e √s (orm=±d √s) of the lengthss. The solution is based on the definition of the mean error of a unit weight founded on the work of deformation. In the calculus of observations, the work of deformation has so far been determined as the deformation work of internal forces from the corrections of the individual quantities after the adjustment is concluded. However, it is possible, as was demonstrated in this paper, to express it as the deformation work of external forces, which act during the adjustment in pointO (Fig. 1) on the auxiliary static system in the shape of a console, and which provide it with the necessary deformation. If the external forces are expressed by means of the tensor of the auxiliary system, the equation for the mean error of the unit weight will provide a relation between this error and the coefficientc, ord. If the mean angular error is selected and if its weight is put equal to one, an equation of the fourth degree is obtained in terms ofc (ord), from which it is possible to compute this coefficient. From the external forces, necessary to produce the deformation of the auxiliary system during its adjustment, the corrections of the individual elements of this system may be determined. If we want to determine the ellipse of errors in one of the polygonal points, it is necessary to investigated the shifts of this point, which occur if the point is acted upon by a unit force first in one and then in another direction, perpendicular to the first. Both shifts represent conjugate radii of the deformation ellipse. From this ellipse we proceed to the ellipse of errors. The general solution is supplemented by a numerical example of adjusting a polygonal traverse with positionally inaccurate points and, for sake of comparison, also by the adjustment of this traverse with positionally accurate points.

Anschrift: Husova 5, Praha 1-Staré Město  相似文献   

19.
Abstract

An exceptional case to the model-independent theory of Knobloch (1995) is presented, by investigating a rotating cylindrical annulus of height H and side wall radii r o and r i, with non-slip, perfectly thermally conducting side walls and thermally insulating stress-free ends. Radial heating permits the possibility of either two- or three-dimensional convective solutions being the preferred mode. An analytical solution is obtained for the two-dimensional case and a numerical solution for the three-dimensional solution, which is also applied to the two-dimensional solution. It is shown that both two- and three-dimensional solutions can be realized depending on the aspect ratio, γ = H/d, where d = r o-r i is the thickness of the annulus, the radii ratio λ = r i/r o and the rotation rate of the model. For γ = O(1) and λ = 0.4, the preferred convective solution is three-dimensional when the Taylor number, T < 102 and two-dimensional for T > 102. For small aspect ratios, γ ? 1, the preferred mode is two-dimensional for all rotation rates.  相似文献   

20.
Summary The equilibrium small ion concentration (n) was measured in aerosols for various ionization rates (q) and particle concentrations (Z). The results agree with the formulaq=bZn+n 2. These experiments were carried out with different particle sizes and the results provided a series of values of the combination coefficientb. The size of the particles was determined by the diffusion box method and the diameter (d) ranged from 0.01 to 0.1 microns. The results fitted an empirical formulab=0.29d.  相似文献   

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

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