首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
This paper addresses inhomogeneities in the short-period S-wave attenuation field in the lithosphere beneath Altai. A technique based on the analysis of the amplitude ratios of Sn and Pn waves is used. High S-wave attenuation areas are identified in the West Altai, which are related to the source zones of recent large earthquakes, viz., the 1990 Zaisan earthquake and the 2003 Chuya earthquake. Associated with the Chuya earthquake, a large ringlike seismogenic structure had been formed since 1976. It is also found that ringlike seismogenic structures are confined to high S-wave attenuation areas unrelated to large historical earthquakes. It is supposed that processes paving the way for strong earthquakes are taking place in these areas. The magnitudes of probable earthquakes are estimated using the earlier derived correlation dependences of the sizes of ringlike seismogenic structures and the threshold values of magnitudes on the energy of principal earthquakes with prevailing focal mechanisms taken into consideration. The sources of some earthquakes are likely to occur near to the planned gas pipeline route from Western Siberia to China, which should be taken into account. The relationship of anomalies in the S-wave attenuation field and the ringlike seismogenic structures to a high content of deep-seated fluids in the lithosphere is discussed.  相似文献   

2.
The features of the seismic regime before the strongest earthquakes of Taiwan in the late 20th (Chi-Chi on September 21, 1999, Mw = 7.6) and the early 21st century (March 31, 2002, Mw = 7.4) are analyzed. Based on 1990–1999 and 1994–2002 data, respectively, retrospective analysis of three seismic regime parameters are studied: the total annual number of earthquakes NΣ in the range of ML = 2.5–5.5 and Mw = 3.0–7.0; the total annual quantity of released seismic energy ΣE, J; and angular coefficient b of earthquake recurrence graphs. Two explicit subperiods are revealed in the course of the seismic regime: quiescence in 1990–1996 before the Chi-Chi earthquake and in 1994–1997 before the March 2002 earthquake; in 1997–1999 and 1998–2002, respectively, seismic activation is observed. Due to the predominance of weak earthquakes during the Chi-Chi earthquake preparation, factor b appeared relatively higher (–1.16 on average); in contrast, before the March 2002 earthquake, due to the occurrence of foreshocks with Mw = 6.8–7.0, the factor b values appeared relatively lower (–0.55 and–0.74 for the quiescence and activation subperiods, respectively). Despite the fundamental difference in the seismotectonic situation between the domains where two mainshocks occurred and significantly difference energy ranges of the initial seismic events, the analysis results are similar for both earthquakes. In both cases, the mainshock occurred at the peak of released energy, which can be considered a coincidence. Solid verification of this positive tendency requires the accumulation of seismological statistics.  相似文献   

3.
The estimated characteristics of the atmospheric boundary layer, obtained by the simulation of wind wave fields using three versions of the WAM numerical model are compared with the well-known empirical dependences of drag coefficient C d on wind speed U 10 and wave age A, as well as with the dependence of dimensionless roughness height z n on inverse wave age u*/с р. Calculations carried out for several years in the areas of the Pacific and Indian oceans, based on the ERA-interim and CFSR wind reanalyses have shown good agreement between the model and empirical dependences C d (U 10) and C d (A). The range of estimated variability for z n (u*/с р ) has been found to be significantly less than empirical. It has been also found that estimated values of wind speed U 10W (t) are overestimated from 5 to 10% in all versions of WAM models compared with the input wind reanalysis U 10R (t) at the moments of appearance maximum values of wind U 10R (t). The reasons for the established features of the WAM model and their dependence on the model version are discussed.  相似文献   

4.
The species composition, phytoplankton abundance, and relative yield of the variable fluorescence (F v /F m ) were determined in the mesotrophic Nhatrang Bay in October–November of 2004. The species diversity (250 taxonomic units) and heterogeneity of the phytoplankton structure were high. With respect to the number of species and their abundance, diatoms prevailed. In selected parts of the bay, dinoflagellates dominated. The mean biomass in the water column under 1 m2 (B t ) varied from 2.3 to 64.4 mg C/m3 being 31.0 mg C/m3 on average. The values of B t were the lowest at the stations nearest to the river mouth. Seaward, B t increased. The values of B t increased with depth at some stations and decreased at others. In the surface sea layers, the biomass was lower than that in the underlying waters. The values of F v /F m ranged from 0.10 to 0.64 (at a mean value of 0.49). The lowest values of F v /F m were observed in the area close to the seaport. Over the greater part of the bay, the values of F v /F m were higher than 0.47. Such values are indicative of the relatively high potential photosynthetic activity of the phytoplankton. The abundance and species diversity were higher than those in the dry season (March–April).  相似文献   

5.
The tsunami warning system in the Russian Far East employs the medium-period magnitude MS (BB) by Vaniek–Soloviev. However, its use may lead to inadequacies and underestimates for the tsunamigenic potential of an earthquake. Specifically, this can happen in the case of a so-called tsunami–earthquake. This kind of earthquakes with a nonstandard spectrum was revealed by H. Kanamori in 1972. This problem can be overcome by using a magnitude scale that deals with longer period seismic waves. This study develops a technique for determining the magnitudes at regional distances (from 70 to 4500 km) using the amplitudes of surface seismic waves of periods of 40 and 80 s. At distances of 70–250 km, the amplitude of the joint group of shear and surface waves is used. For the new magnitudes designated M S(40) and M S(80), experimental calibration curves are constructed using more than 1250 three-component records at 12 stations of the region. The magnitudes are calibrated so as to produce an unbiased estimate of the moment magnitude M w in the critical range 7.5–8.8. The rms error of the single-station estimate M w is around 0.27. At distances below 250 km and M w ≥ 8.3, the estimate of M w obtained by the proposed technique becomes saturated at the level of M w ~ 8.3, which is acceptable for operative analysis because no missed alarms arise. The technique can be used in operational tsunami warning based on seismological data. This can markedly decrease the number of false alarms.  相似文献   

6.
Some characteristics of seismicity in Southern California are studied. It is found that ring-shaped seismicity structures with threshold magnitudes Mth of 4.1, 4.1, and 3.8 formed prior to three large (M w > 7.0) earthquakes in 1992, 1999, and 2010, respectively. The sizes of these structures are several times smaller than for intracontinental strike-slip events with similar magnitudes. Two ring-shaped structures are identified in areas east of the city of Los Angeles, where relatively large earthquakes have not occurred for at least 150 years. The magnitudes of large events which can occur in the areas of these structures are estimated on the basis of the previously obtained correlation dependence of ring sizes on magnitudes of the strike-slip earthquakes. Large events with magnitudes of M w = 6.9 ± 0.2 and M w = 8.6 ± 0.2 can occur in the area to the east of the city of Los Angeles and in the rupture zone of the 1857 great Fort Tejon earthquake, respectively. We believe that ring-structure formation, similarly to the other regions, is connected with deep-seated fluid migration.  相似文献   

7.
Chlorophyll-a concentration (C chl) variations in the cross section within and outside the Peter the Great Bay shelf during different stages of the winter–spring phytoplankton bloom in 2003–2005 has been considered based on a ship (obtained during the R/V Akademik M.A. Lavrent’ev voyage of February 26 to March 9, 2003) and MODIS-Aqua spectroradiometer and the SeaWiFS color-scanner satellite data. A comparison of the C chl variability obtained from the ship and satellite data indicates that these data are inconsistent. According to satellite data obtained at the MUMM atmospheric correction, the C chl variability is distorted less than the NIR-correction data. Studying the variations in the coefficients of light absorption by the detritus and yellow substance (a dg) and light backscattering by suspended particles (b bp), C chl, chlorophyll-a fluorescence (F chl) according to the satellite data allow us to state that the variations in the discrepancy between the satellite and ship C chl values are mainly caused by the variations in the content of the detritus and yellow substance in water. Based on the satellite data, it has been revealed that the a dg values increase with increasing wind mixing after the phytoplankton bloom (about 2–5 km areas where the a dg, C chl, F chl, and bbp values abruptly increased in 2005, apparently due to eddy formation). It has been indicated that the F chl characteristic, which is close to C chl, increases when the favorable conditions for the phytoplankton bloom deteriorate. Therefore, this characteristic cannot be used to identify C chl under the indicated conditions.  相似文献   

8.
To explore the effects of temperature changes on dinoflagellate bloom succession in the coastal waters of the East China Sea, changes in the growth, photophysiology, and Rubisco gene expression of Prorocentrum donghaiense and Karenia mikimotoi, two harmful algal species, were investigated at different temperatures (16 to 28°C). The maximal specific growth rate and the maximal mRNA expression of Rubisco gene in P. donghaiense and K. mikimotoi occurred at 20 and 24°C, respectively. The photosynthetic activity of P. donghaiense was generally stable, but K. mikimotoi photosynthesis increased when temperatures rose from 16 to 28°C. The effective photochemical efficiency (F q /F m ) and the maximal relative electron transfer rate (rETRmax) of K. mikimotoi increased significantly with increasing temperature, and the lowest and highest values occurred at 16 and 28°C, respectively. It seems that P. donghaiense has higher photosynthetic capacity than K. mikimotoi due to its higher F q /F m , rETRmax, and photosynthetic efficiency (α). However, K. mikimotoi has a higher growth rate than P. donghaiense. These results suggest that the photosynthetic activity and genetic responses of dinoflagellates are species-dependent. It is likely that temperature changes affect species composition during blooms, leading to the observed patterns of bloom succession.  相似文献   

9.
引入拖曳系数参数化的海冰自由漂流模拟研究   总被引:2,自引:0,他引:2  
Many interesting characteristics of sea ice drift depend on the atmospheric drag coefficient(C_a) and oceanic drag coefficient(C_w).Parameterizations of drag coefficients rather than constant values provide us a way to look insight into the dependence of these characteristics on sea ice conditions.In the present study,the parameterized ice drag coefficients are included into a free-drift sea ice dynamic model,and the wind factor α and the deflection angle θ between sea ice drift and wind velocity as well as the ratio of C_a to C_w are studied to investigate their dependence on the impact factors such as local drag coefficients,floe and ridge geometry.The results reveal that in an idealized steady ocean,C_a/C_w increases obviously with the increasing ice concentration for small ice floes in the marginal ice zone,while it remains at a steady level(0.2-0.25) for large floes in the central ice zone.The wind factor α increases rapidly at first and approaches a steady level of 0.018 when A is greater than 20%.And the deflection angle θ drops rapidly from an initial value of approximate 80° and decreases slowly as A is greater than20%without a steady level like α.The values of these parameters agree well with the previously reported observations in Arctic.The ridging intensity is an important parameter to determine the dominant contribution of the ratio of skin friction drag coefficient(C_s' /C_s) and the ratio of ridge form drag coefficient(C_r'/C_r) to the value of C_a/C_w,α,and θ,because of the dominance of ridge form drag for large ridging intensity and skin friction for small ridging intensity among the total drag forces.Parameterization of sea ice drag coefficients has the potential to be embedded into ice dynamic models to better account for the variability of sea ice in the transient Arctic Ocean.  相似文献   

10.
The paper presents a new model of electronic-vibrational kinetics of the products of ozone and molecular oxygen photodissociation in the terrestrial middle atmosphere. The model includes 45 excited states of the oxygen molecules O2(b 1, Σ g + ,v= 0−2), O2 (a 1Δ g , v= 0−5), and O2(X 3Σ g , v= 1−35) and of the metastable atom O (1 D) and over 100 aeronomic reactions. The model takes into account the dependence of quantum yields of the production of O2(a 1Δ g , v= 0−5) in a singlet channel of ozone photolysis in the Hartley band on the wavelength of photolytic emission. Taking account of the electronic-vibrational kinetics is important in retrieval of the vertical profiles of ozone concentration from measured intensities of the Atm and IR Atm emissions of the oxygen bands above 65 km and leads to an increase in the ozone concentration retrieved from the 1.27-μm emission, in contrast to the previous model of pure electronic kinetics. Sensitivity analysis of the new model is made for variations in the concentrations of atmospheric constituents ([O2], [N2], [O(3P)], [O3], [CO2]), the gas temperature, rate constants of the reactions, and quantum yields of the reaction products. A group of reactions that most strongly affect the uncertainty of ozone retrieval from measured intensities of atmospheric emissions of molecular oxygen O2(b 1Σ g + , v) and O2(a 1Δ g , v) has been determined. Original Russian Text ? V.A. Yankovsky, V.A. Kuleshov, R.O. Manuilova, A.O. Semenov, 2007, published in Izvestiya AN. Fizika Atmosfery i Okeana, 2007, Vol. 43, No. 4, pp. 557–569.  相似文献   

11.
The possibility of estimating the parameters of surface pulsed sources from data on acoustic waves recorded in the atmosphere is studied. Experimental values are given for peak pressure P + of recorded acoustic signals, wave-profile area S + in their positive phase, and length t + of this phase, and the approximations of these parameters are obtained within wide ranges of source energy 10–3 < E < 1010 kg TNT and scaled distances 1 < R/E 1/3 < 4 × 104 m/kg1/3. Conventional methods of estimating the acoustic energy E according to data obtained from acoustic measurements in the atmosphere are analyzed, and ways to improve their accuracy are proposed. The influence of the type of explosions on the parameters P +, S +, and t + of acoustic signals at long distances R/E 1/3 > 500 m/kg1/3 from explosions is shown.  相似文献   

12.
The results of joint analysis of temperature variations near mesopause from long-term measurements at the Zvenigorod Scientific Station of the Obukhov Institute of Atmospheric Physics RAS in 1960–2015 and variations of surface air temperature characterizing global climate change. Together with variations of temperature at the mesopause T ms from measurements of the hydroxyl emissions we analyzed the temperature variations near mesopause T m reduced to the same level of solar activity. The observed strong decrease in temperature near mesopause during last decades, particularly in winter, with its tendency to slow down since the 1980’s is was detected against the background of general increase in the surface air temperature of the Northern Hemisphere T NHs and the Earth as a whole. It was revealed a sharp drop in winter temperature near mesopause in 1970s. and its synchronicity with the shift in climatic features at the surface associated with changes in formation of El Nino events and their impact on the global climate. The general significant negative correlation of temperature variations near mesopause and T NHs detected from 56-year observational data was not accompanied by any significant coherence between the most long-period temperature variations from the cross-wavelet analysis. To assess the possible manifestation of this coherence the results of numerical simulations with a global climate model were used. According to model simulations for the 20–21 centuries taking into account anthropogenic forcings for significant coherence between long-term variations T m and T NHs the prolonged observations are required for temperature near mesopause–about a century or more.  相似文献   

13.
The atmospheric effect on the characteristics of infrasonic signals from explosions has been studied. New methods have been proposed to remotely estimate the energy of explosions using the data of infrasonic wave registration. One method is based on the law of conservation of acoustic pulse I, which is equal to the product of the wave profile area S/2 of the studied infrasonic signal and the distance to the source EI [kt] = 1.38 × 10–10 (I [kg/s])1.482. The second method is based on the relationship between the explosion energy and the dominant period T of the recorded signal, EТ [kt] =1.02 × (Т [s]2/σ)3/2, where σ is a dimensionless distance used for determining the degree of manifestation of nonlinear effects in the propagation of sound along ray trajectories. When compared to the conventional EW (Whitaker’s) relation, the advantage of the EI relation is that it can be used for pulsed sources located at an arbitrary height over the land surface and having an arbitrary form of the initial-pulse profile and for any type of infrasonic arrivals. A distinctive feature of the expression for EТ is that the atmospheric effect on the characteristics of recorded infrasonic signals is explicitly taken into account. These methods have been tested using infrasonic data recorded at a distance of 322 km from the sources (30 explosions caused by a fire that occurred at the Pugachevo armory in Udmurtia on June 2, 2011). For the same explosion, empirical relations have been found between energy values obtained by different methods: EI = 1.107 × E W , E Т = 2.201 × E I .  相似文献   

14.
海水石油类物质吸收系数遥感化提取算法研究   总被引:2,自引:0,他引:2  
Establishing the remote sensing algorithm of retrieving the absorption coefficient of seawater petroleum substances is an efficient way to improve the accuracy of retrieving a seawater petroleum concentration using a remote sensing technology. A remote sensing reflectance is a basic physical parameter in water color remote sensing. Apply it to directly retrieve the absorption coefficient of seawater petroleum substances is of potential advantage. The absorption coefficient of waters containing petroleum [ACWCP, a_o(λ)], consists of the absorption coefficient of pure water [ACPW, a_w(λ)], plankton [ACP, a_(ph)(λ)], colored scraps [ACCS, a_(d,g)(λ)], and petroleum substance [ACPS, a_(oil)(λ)]. Among those, ACCS consists of the absorption coefficient of nonalgal particle [ACNP, a_d(λ)] and colored dissolved organic matter [ACCDOM, a_g(λ)]. For waters containing petroleum, the retrieved ACCS using the existing method is a combination absorption coefficient of ACNP,ACCDOM and ACPA [CAC, a_(d,g,oil)(λ)]. Therefore, the principle question is how to extract ACPS from CAC.Through the analysis of the three proportion tests conducted between the year of 2013 and 2015 and the corresponding remote sensing data, an algorithm of retrieving the absorption coefficient of petroleum substances is proposed based on remote sensing reflectance. First of all, ACPS and CAC are retrieved from the reflectance using the quasi-analytical algorithm(QAA), with some parameter modified. Secondly, given the fact that the backscatter coefficient [BC, b_(bp)(555)] of total particles at 555 nm can be obtained completely from the reflectance, the relation between BC and ACNP in petroleum contaminated water can be established. As a result, ACNP can be calculated. Then, combining the remote sensing retrieving algorithm of a_g(440), the method of achieving the spectral slope of the absorption coefficient can be established, from which ACCDOM,can be calculated. Finally, ACPS can be computed as the residual. The accuracy of ACPS based on this algorithm is 86% compared with the in situ measurements.  相似文献   

15.
The current study was carried out from May 2014 to April 2015 to estimate the stock status of P. viridis in Marudu Bay. The gonad development was monitored by histological examination, while the population parameters including asymptotic length (L), growth coefficient (K), mortality rate (Z, F and M), exploitation level (E) and recruitment of P. viridis were estimated using the lengthfrequency data. Results of the current study demonstrated that P. viridis in Marudu Bay spawned throughout the year with two major peaks, one in April to May and another one in October to December. The recruitment pattern was continuous with the peak in May to June 2014, which corresponded to the first spawning peak in April. However, no significant recruitment was observed from the second spawning peak due to the difference in spawning timing between male and female populations. The estimated asymptotic length (L), growth coefficient (K), total mortality (Z), natural mortality (M), fishing mortality (F) and growth performance (ф) of P. viridis in Marudu Bay were estimate to be 117 mm, 0.97 yr-1, 4.39 yr-1, 1.23 yr-1, 3.16 yr-1 and 4.123, respectively. The exponent b of the lengthweight relationship was 2.4 and exploitation level (E) was 0.72. The high mortality, low condition indices and negative allometric of P. viridis in Marudu Bay is caused by a lack of suitable food in the surrounding water.  相似文献   

16.
Age, growth and mortality of Mactra chinensis were investigated during the period from October 2012 to September 2013 in Busan, South Korea. The monthly variation of the marginal index (MI) of the shell and chondrophore showed that the ring of this species was formed once a year during July. We estimated the age of M. chinensis by reading the external rings on the shell and the growth bands of the chondrophore to compare growth parameters between the two growth characters. The age of this species ranged from 0 to 8 years (shell-based age reading) and from 0 to 10 years (chondrophore-based age reading). Based on external rings and growth bands of chondrophore for the same period, the von Bertalanffy growth functions were expressed by the equation, L t = 101.53[1-exp {-0.15(t + 0.75)}] and L t = 90.03[1-exp {-0.20(t + 0.50)}], respectively. The likelihood test showed that there was a significant difference in L (P < 0.001), K (P < 0.001), to (P < 0.001) estimated from non-linear regression between the two growth characters.  相似文献   

17.
C波段紧缩极化SAR海冰探测能力评估   总被引:1,自引:0,他引:1  
The C-band synthetic aperture radar(SAR) data from the Bohai Sea of China, the Labrador Sea in the Arctic and the Weddell Sea in the Antarctic are used to analyze and discuss the sea ice full polarimetric information reconstruction ability under compact polarimetric modes. The type of compact polarimetric mode which has the highest reconstructed accuracy is analyzed, along with the performance impact of the reconstructed pseudo quad-pol SAR data on the sea ice detection and sea ice classification. According to the assessment and analysis, it is recommended to adopt the CTLR mode for reconstructing the polarimetric parameters σ_(HH)~0,σ_(VV)~0, H and α,while for reconstructing the polarimetric parameters σ_(HV)~0, ρ_(H-V), λ_1 and λ_2, it is recommended to use the π/4 mode.Moreover, it is recommended to use the π/4 mode in studying the action effects between the electromagnetic waves and sea ice, but it is recommended to use the CTLR mode for studying the sea ice classification.  相似文献   

18.
Cascading of cold Antarctic shelf water (ASW) initiates compensatory isopycnic upwelling of the warm Circumpolar Deep Water (CDW). The baroclinic/thermoclinic Antarctic slope front (ASF) is formed, and a mesoscale intrusive structure develops on the shelf edge and slope. Mesoscale processes when the ASF peaks are periodically accompanied by local baroclinic instability, which forms a smaller-scale intrusive structure. Therefore, the ASF is naturally subdivided into two layers according to the intrusion scales (vertical δН and horizontal L) and the horizontal parameters of the front (thermoclinity (TL)ρ and baroclinity γρ). Analysis of ASF intrusive layering due to the baroclinic factor supports the following conclusion: the higher the (TL)ρ of the ASF, the greater the intrusion intensity |δθ| (temperature anomaly amplitude), while an increase in γρ of the ASF leads to a decrease in intrusion scales δН and L. Frontal intrusions can be distinguished by a development degree. Regardless of the degree of development, all warm intrusions are characterized by vertical density stratification, while cold intrusions are characterized by density quasihomogeneity. According to field data, the ASF instability process is subdivided into four stages. When theASF is baroclinically unstable, the local baroclinic deformation radius RdL of the front is close in magnitude to the horizontal scale L of the intrusions that form, and their characteristic vertical scale δH is close to the typical vertical scale of front instability.  相似文献   

19.
The age and growth of filefish, Thamnaconus modestus (Günther 1877) in the southern waters of Korea were investigated. Samples were collected with commercial trawl catches during the period from May 2009 to December 2011. Of the 2,626 specimens collected, the sex ratio was not significantly different from 1:1 (P > 0.05). The total length ranged from 11.3 to 42.1 cm. The gonadosomatic index for both sexes was the highest in May to June, indicating that May to June is the main spawning period. The length of females at sexual maturity was 25.92 cm. The length-weight relationship of the filefish was TW = 0.0121TL3.0536 (n = 1,692, r2 = 0.9034, P < 0.001). The age of the sampled individuals was estimated by counting growth rings recorded on the 5th vertebrae; ages ranged from 0 to 9 years. The filefish of the same age displayed a high individual variation in total length. Length-at-age data were fitted by using the Von Bertalanffy growth model. The estimated Von Bertalanffy growth parameters were L = 42.04 cm, k = 0.21 year?1 and t0 = ?1.56 for females, L = 41.20 cm, k = 0.18 year?1 and t0 = ?2.36 for males, and L = 43.16 cm, k = 0.17 year?1 and t0 = ?2.18 for the combination of both male and female. These data can be used as useful biological information for the future fishery management of filefish resources in Korean waters.  相似文献   

20.
Using the three-dimensional nonlinear finite-element thermohydrodynamic model QUODDY-4, we obtain two solutions giving an idea of the role of the nonlinear interaction of tidal harmonics in the formation of their spatial structure. The first solution is induced by the total tide (M 2 + S 2 + K 1 + O 1) at the open boundary and by the total static tide inside the area under consideration; then, the solution obtained is subject to harmonic analysis. The second solution is obtained by specifying such tidal sea level elevations at the open boundary that meet individual tidal harmonics. These two solutions are compared. It is shown that the differences between the solutions for the S 2, K 1, and O 1 waves can be significant, especially near the open boundary between the White and Barents seas. This conclusion remains valid also for the maximum velocities (major semiaxes of the ellipses) of the barotropic (depth-average) tidal current as well as for the average (over the tidal cycle) densities of the total tidal energy and components of the tidal-energy budget. The emergence of this feature indicates that there are resonance modes with frequencies that differ from those of the S 2, K 1, and O 1 harmonics to a lesser extent than the M 2 harmonic frequency. The same conclusion can be made by comparing the values of the amplification factor, which is defined as the ratio between the actual and static tides, for the system of the Barents and White seas as a whole and for Mezen Bay in the White Sea and Czech Bay in the Barents Sea, taken separately.  相似文献   

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

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