首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
The results obtained in this study show that as the dart leader tip passes a given point on the defunct return stroke channel the electric field increases within a fraction of a microsecond to values larger than the critical electric field necessary for the initiation of cold electron runaway in low-density air comprising the channel. These results are in support of the hypothesis that cold runaway electron breakdown may play a role in the emission of X-ray bursts by dart leaders. The calculations also show that the peak power dissipated by a typical dart leader is about 300–500 MW/m and the energy dissipated within the first 10 μs or so is about 500–600 J/m. Furthermore, the minimum resistance and the maximum radius of the core of a typical dart leader are estimated to be about 3 Ω/m and 0.003 m, respectively.  相似文献   

2.
An interpretation of the nature of the sudden ionospheric disturbance in terms of response to X-ray flux enhancement in the band 1–20 Å has been made by many authors. Last decades investigations revealed presence of important qualitative distinctions in spatially temporal pattern of geomagnetic response to solar flares featuring harder radiation spectra (with quanta energies above 100 keV). These distinctions can not be adequately described by classical theory implying ionization growing on E and D ionosperic layers and intensification of Sq-current system. In this respect, solar flare on 4 November 2003 characterizing by existence of two separate (time lag ~45 min) spectral maximums in X-rays range (average quantum energy <100 keV) and in γ-rays range (average quantum energy >100 keV), represents convenient proving ground for study of specifics the geomagnetic response to bursts marked by different hardness. In current article, we show that this flare has a number of unusual features including specific variation of accompanying current system and magnetospheric manifestation that is observed in trapped radiation fluxes and magnetic field on geosynchronous orbit. Possible physical mechanism leading to intensification of magnetospheric–ionospheric current system is discussed.  相似文献   

3.
Narrow bipolar events (NBEs) are a distinct class of intra-cloud lightning discharge. In this paper we present observations of 10 negative and 67 positive such events in East China. Positive NBEs occurred at 7–12 km altitude above mean sea level (MSL) with a mean altitude of 9.5 km, and negative NBEs occurred at 14–16 km altitude. Electrical/channel characteristics of these events were derived from NBE pulse waveforms based on the transmission-line model. On average, the peak current moment and the charge moment change of a NBE event is 15 kA km, and 0.12 C km, respectively. The mean time for the propagation of current front along the channel is 2.2 μs. The upper limit on channel length for NBEs in this study is 510–1060 m, the lower limit on discharge current amplitude is 12.5–43.2 kA, and the minimum charge transfer is 0.1–0.3 C.  相似文献   

4.
We have determined the post-perovskite phase transition boundary in MgSiO3 in a wide temperature range from 1640 to 4380 K at 119–171 GPa on the basis of synchrotron X-ray diffraction measurements in-situ at high-pressure and -temperature in a laser-heated diamond-anvil cell (LHDAC). The results show a considerably high positive Clapeyron slope of + 13.3 ± 1.0 MPa/K and a transition temperature of about 3520 ± 70 K at the core–mantle boundary (CMB) pressure. The thermal structure in D″ layer can be tightly constrained from precisely determined post-perovskite phase transition boundary and the depths of paired seismic discontinuities. These suggest that temperature at the CMB may be around 3700 K, somewhat lower than previously thought. A minimum bound on the global heat flow from the core is estimated to be 6.6 ± 0.5 TW.  相似文献   

5.
Variations of electron fluxes with energies 300–600 keV in the region of quasitrapping are analyzed using data of the low orbiting Coronas-F satellite. Enhancements in the electron fluxes with energies above 300 keV are observed at the polar boundary of the outer radiation belt. Meteor-3M satellite data, OVATION and AP models of the position of the auroral oval are used to determine the position of analyzed increases in the energetic electrons with respect to the position of the auroral oval. There is a significant number of events when these increases were observed at a few consequent orbits crossing the outer radiation belt boundary. Studied increases in relativistic electron fluxes are localized at the latitudes of the auroral oval. Different mechanisms of formation of observed enhancements are discussed. The possibility of the appearance of increases due to formation of local particle traps is analyzed using Tsyganenko geomagnetic field models. The role of the formation of local particle traps at the boundary of the outer radiation belt and its possible influence to the formation of the outer radiation belt is discussed.  相似文献   

6.
Radiation with energies up to about 250 keV associated with the dart leader phase of rocket-triggered lightning were reported by Dwyer et al. (2004). The mechanism of X-ray generation by dart leaders, however, is unknown at present. Recently, Cooray et al., in pressPlease provide complete bibliographic details for Ref. cooray e al., in press and Saleh et al., in press if available., Cooray et al., 2009a developed physical concepts and mathematical techniques necessary to calculate the electric field associated with the tip of dart leaders. We have utilized the results of these calculations together with the energy dependent frictional force on electrons, as presented by Moss et al. (2006), to evaluate the maximum energy an electron will receive in accelerating in the dart-leader-tip electric field. The main assumptions made in performing the calculations are: (a) the dart leader channel is straight and vertical; (b) the path of the electrons are straight inside the channel; and (c) the decay of the channel temperature is uniform along the length of the dart leader. In the calculation, we have taken into account the fact that the electric field is changing both in space and time and that the gas in the defunct return stroke channel is at atmospheric pressure and at elevated temperature (i.e. reduced gas density). The results of the calculation show that for a given dart leader current there is a critical defunct-return-stroke-channel temperature above which the cold electron runaway becomes feasible. For a typical dart leader, this temperature is around 2500 K. This critical temperature decreases with increase in dart leader current. Since the temperature of the defunct return stroke channel may lie in the range of 2000–4000 K, the results show that the electric field at the tip of dart leaders is capable of accelerating electrons to MeV energy levels.  相似文献   

7.
We have considered spatial distributions of positive lightning discharges in the east of Siberia for the summer seasons of 2003–2007 and properties of their electromagnetic signals with the ELF “slow tail”, which, as is known, can be accompanied by sprites. There are two main regions of positive discharges located in the south and west of Yakutsk. Two other “centers” (the northeastern and the eastern) are located in high-mountainous regions. In these regions the positive discharges intensity can exceed the negative discharge intensity.The electromagnetic signals in the ELF range (usually in the form of two half-cycles) were observed after the VLF atmospherics were recorded in the high-latitude regions. The delay of ELF pulses relative to the corresponding atmospherics was 0–7 ms. The long (up to 350 ms) events of quasi-periodic ELF oscillations with the period of about 7 ms (which corresponds to the quasi-period of ELF pulses) were revealed.  相似文献   

8.
The suggestion that the polar cap can completely disappear under certain northward IMF conditions is still controversial. We know that the size of the polar cap is strongly controlled by the interplanetary magnetic field (IMF). Under a southward IMF, the polar cap is usually large and filled with weak diffuse polar rain electrons. The polar cap shrinks under a northward IMF. Here we use the global auroral images and coincident particle measurements on May 15, 2005 to show that the discrete arcs (due to precipitation of both electrons and ions) expanded from the dayside oval to the nightside oval and filled the whole polar ionosphere after a long (8 h) and strong (~5–30 nT) northward IMF Bz, The observations suggested that the polar cap disappeared under a closed magnetosphere.  相似文献   

9.
Equivalent dose (De) values were measured by using medium aliquots of different grain size quartz fractions of five lakeshore sediments from the arid region of north China. There are two different relationships between De values and grain sizes of these five samples. The first relationship is that the De values obtained from various grain sizes are in agreement within 1 delta errors. The second relationship is that De values are similar to each other for fractions between 125 and 300 μm, while the De value of the 63–90 μm fraction is 40~55% smaller than others. For example, the De values obtained for sample #3 are 20.15 ± 1.19 Gy, 19.80 ± 0.83 Gy and 20.93 ± 1.06 Gy for fractions of 90–125, 125–150 and 250–300 μm respectively, but are 10.79 ± 0.84 Gy for the 63–90 μm fraction. The second relationship can't be interpreted by previous studies of both dosimetry and heterogeneous bleaching. It is deduced for sample #2, #3 and #6 that fine particles (<90 μm) intruded after the dominant sedimentation. Comparison of OSL ages from different grain size fractions of sample #2 with a radiocarbon age from the same lithologic layer supports that fractions coarser than 125 μm yield more reliable burial ages, while the fraction finer than 90 μm yields underestimated ages for some lakeshore sediments from this arid region.  相似文献   

10.
The southernmost segment of the Andes of southern Patagonia and Tierra del Fuego forms a ~ 700 km long orogenic re-entrant with an interlimb angle of ~ 90° known as Patagonian orocline. No reliable paleomagnetic evidence has been gathered so far to assess whether this great orogenic bend is a primary arc formed over an articulated paleomargin, or is due to bending of a previously less curved (or rectilinear) chain. Here we report on an extensive paleomagnetic and anisotropy of magnetic susceptibility (AMS) study carried out on 22 sites (298 oriented cores), predominantly sampled in Eocene marine clays from the external Magallanes belt of Tierra del Fuego. Five sites (out of six giving reliable paleomagnetic results) containing magnetite and subordinate iron sulphides yield a positive fold test at the 99% significance level, and document no significant rotation since ~ 50 Ma. Thus, the Patagonian orocline is either a primary bend, or an orocline formed after Cretaceous–earliest Tertiary rotations. Our data imply that the opening of the Drake Passage between South America and Antarctica (probably causing the onset of Antarctica glaciation and global climate cooling), was definitely not related to the formation of the Patagonian orocline, but was likely the sole consequence of the 32 ± 2 Ma Scotia plate spreading. Well-defined magnetic lineations gathered at 18 sites from the Magallanes belt are sub-parallel to (mostly E–W) local fold axes, while they trend randomly at two sites from the Magallanes foreland. Our and previous AMS data consistently show that the Fuegian Andes were characterized by a N–S compression and northward displacing fold–thrust sheets during Eocene–early Miocene times (50–20 Ma), an unexpected kinematics considering coeval South America–Antarctica relative motion. Both paleomagnetic and AMS data suggest no significant influence from the E–W left-lateral Magallanes–Fagnano strike–slip fault system (MFFS), running a few kilometres south of our sampling sites. We thus speculate that strike–slip fault offset in the Fuegian Andes may range in the lower bound values (~ 20 km) among those proposed so far. In any case our data exclude any influence of strike–slip tectonics on the genesis of the great orogenic bend called Patagonian orocline.  相似文献   

11.
SHRIMP zircon U–Pb dating, mineral chemical, element geochemical and Sr–Nd–Pb–Hf isotopic data have been determined for the Yulong monzogranite-porphyry in the eastern Tibet, China. The Yulong porphyry was emplaced into Triassic strata at about 39 Ma. The rocks are weakly peraluminous and show shoshonitic affinity, i.e., alkalis-rich, high K2O contents with high K2O / Na2O ratios, enrichment in LREE and LILE. They also show some affinities with the adakite, e.g., high SiO2 and Al2O3, and low MgO contents, depleted in Y and Yb, and enrichment in Sr with high Sr / Y and La / Yb ratios, and no Eu anomalies. The Yulong porphyry has radiogenic 87Sr / 86Sr (0.7063–0.7070) and unradiogenic 143Nd / 144Nd (εNd =  2.0 to − 3.0) ratios. The Pb isotopic compositions of feldspar phenocrysts separated from the Yulong porphyry show a narrow range of 206Pb / 204Pb ratios (18.71–18.82) and unusually radiogenic 207Pb / 204Pb (15.65–15.67) and 208Pb / 204Pb (38.87–39.00) ratios. In situ Hf isotopic composition of zircons that have been SHRIMP U–Pb dated is characterized by clearly positive initial εHf values, ranging from + 3.1 to + 5.9, most between + 4 and + 5. Phenocryst clinopyroxene geothermometry of the Yulong porphyry indicates that the primary magmas had anomalously high temperature (> 1200 °C). The source depth for the Yulong porphyry is at least 100 km inferred by the metasomatic volatile phase (phlogopite–carbonate) relations. Detailed geochemical and Sr–Nd–Pb–Hf isotopic compositions not only rule out fractional crystallization or assimilation-fractional crystallization processes, but also deny the possibility of partial melting of subducted oceanic crust or basaltic lower crust. Instead, low degree (1–5%) partial melting of a metasomatized lithosphere (phlogopite–garnet clinopyroxenite) is compatible with the data. This example gives a case study that granite can be derived directly by partial melting of an enriched lithospheric mantle, which is important to understand the source and origin of diverse granites.  相似文献   

12.
The Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), a limb sounding Fourier transform interferometer in the mid-infrared band, on board the polar-orbiting ENVIronmental SATellite (ENVISAT) was launched in March 2002 by the European Space Agency (ESA). For the MIPAS data validation, three balloon flights with MIPAS-B, a balloon version of MIPAS, were carried out on the 24th and 25th of September 2002 from Aire sur l’Adour (France 44°N, 0°E), on the 20th and 21th of March, and on the 2nd and 3rd of July 2003 from Esrange, Kiruna (Sweden 68°N, 21°E). The MIPAS operational data version 4.61 for the temperature vertical profiles were compared with the correlative MIPAS-B measurements via the coincident comparison and the trajectory comparison approaches, respectively. The precision of the MIPAS temperature was estimated to be 1.04–2.48 K in the region 123–8 hPa (15–33 km) and beyond the expectation. The systematic difference shows that the absolute accuracy of the MIPAS temperature is within the total combined errors between 356 and 3 hPa (8–39 km). The agreements in the middle stratosphere are better than in the lower stratosphere and upper troposphere. The maximum difference of 1.9 K (0.9%) was found at 228 hPa (11 km). The biases of MIPAS measurements are ?0.39 K (?0.19%) and +0.25 K (+0.12%) (“+” positive bias, “?” negative bias) corresponding to the altitude regions 143–31 hPa (14–24 km) and 31–5 hPa (24–37 km), respectively.  相似文献   

13.
The propagation features of nighttime whistlers to low-latitude station, Suva (−18.2°, 178.3°, geomag. lat. −22.1°, geomag. long. 253.5°, L=1.15), Fiji, from preliminary observations made during the period from September 2003–2005, are reported. The observations of ELF–VLF signals commenced in September 2003 using the VLF set-up of World Wide Lightning Location Network at our station. The whistlers were observed during the severe magnetic storm of 20–22 November 2003 and moderate magnetic storm of 17–19 July 2005. A whistler with dispersion D=12.7 s1/2 occurred on 22 November at 00:11 h LT. On 20 July at 01:00 h LT, a short whistler with dispersion D=20.9 s1/2 and two whistler events having two-component whistlers with D=15.8, 16.7 s1/2 and 16.7, 17.3 s1/2 were observed. Non-ducted pro-longitudinal mode of the whistler propagation supported by negative latitudinal electron density gradients in the ionosphere that are enhanced by magnetic storms, seems most likely mode of propagation for the whistlers with dispersion of 12.7–17.3 s1/2 to this low-latitude station.  相似文献   

14.
The characteristics (annual and diurnal cycle, polarity, multiplicity and first stroke peak current) of ∼4.3×106 cloud-to-ground flashes recorded in the Iberian Peninsula during the first decade of measurements of the lightning detection network installed in Spain are analyzed. The mean monthly variation shows maximum lightning activity between May and September, while minimum values are observed in January and February. The mean diurnal cycle shows maximum values at 1700 LT and minimum values at around 1000 LT. The average maximum flash density (not corrected for detection efficiency) is 2.1 flashes km−2 year−1. Maximum lightning activity is associated with mountainous areas. The effect of the Mediterranean Sea is also seen. The percentage of positive flashes is 9%, although this changes over the year from 6.5% in June to 22.6% in January. The average multiplicity is found to be 2.0 for the negative flashes and 1.1 for the positive flashes, and the percentages of single-stroke flashes are 53.6% and 89%, respectively. The monthly distribution of multiplicity for negative flashes peaks in the summer and minimum values are found in the winter. The multiplicity of the positive flashes does not seem to be function of the month. The median (mean) first stroke peak currents are found to be 23.5 kA (27.3 kA) for the negative flashes and 35.3 kA (47.1 kA) for the positive flashes. For both polarities, the peak current is higher in the summer than in the winter. The percentage of positive flashes and the mean peak currents for negative flashes are higher over the sea areas than over land.  相似文献   

15.
Explosion deep seismic sounding data sections of high quality had been obtained with RV Meteor in the Reykjanes Iceland Seismic Project (RRISP77 [Angenheister, G., Gebrande, H., Miller, H., Goldflam, P., Weigel, W., Jacoby, W.R., Pálmason, G., Björnsson, S., Einarsson, P., Pavlenkova, N.I., Zverev, S., Litvinenko, I.V., Loncarecic, B., Solomon, S., 1980. Reykjanes Ridge Iceland Seismic Experiment (RRISP 77). J. Geophys. 47, 228–238]) which close an information gap near 62°N. Preliminary results were presented by Weigel [Weigel, W., 1980. Aufbau des Reykjanes Rückens nach refraktionsseismischen Messungen. In: Weigel, W. (Ed.), Reykjanes Rücken, Island, Norwegischer Kontinentalrand. Abschlusskolloquium, Hamburg zur Meteor-Expedition, vol. 45. DFG, Bonn, pp. 53–61], and here we report on the data and results of interpretation. Clear refracted phases to 90 km distance permit crustal and uppermost mantle structure to be modelled by ray tracing. The apparent P-wave velocities are around 4.5, 6–6.5, 7–7.6 and 8.2–8.7 km/s, but no wide-angle reflections have been clearly seen. Accompanying sparker reflection data reveal thin sediment ponds in the axial zone and up to 400 m thick sediments at 10 Ma crustal age. Ray tracing reveals the following model below the sediments: (1) a distinct, 1–2 km thick upper crust (layer 2A) with Vp increasing with age (to 10 Ma) from <3.4 to 4.9 km/s and with a vertical gradient of 0.1–0.2 km/s/km, (2) a lower crust or layer 3 beginning at depths of 2 (axis) to 4 km (10 Ma age) below sea level with 6.1–6.8 km/s and similar vertical gradients as above, (3) the lower crust bottoms at 5.2–9.5 km depth below sea level (0–10 Ma) with a marked discontinuity, underneath which (4) Vp rises from about 7.5–7.8 km/s (0–10 Ma) with a positive vertical gradient of, again, 0.1–0.2 km/s/km such that 8 km/s would be reached at 12 km and deeper near the axis. Our preferred interpretation is that the mantle begins at the distinct discontinuity (“Moho”), but a deeper “Moho” of Vp  8 km/s cannot be excluded. From Iceland southward to 60°N several experiments show a decrease of crustal thickness from 14 to 8 km. Velocity trends with age across the ridge reflect cooling and filling of cracks, and thickness trends probably suggest volcanic productivity variations as previously suggested.Gravity inversion concentrates on a profile across the ridge with the above seismic a priori information; with 0.2–0.5 km depth uncertainty it leads to a good fit (±2.5 mGal where seismic data exist). Best fitting densities are (in kg/m3) for sediments, 2180; upper crust, 2450–2570; lower crust, 2850–2940; mantle lithosphere, 3215–3240 with a deficit for an asthenospheric wedge of no more than −100 kg/m3. The morphological ridges and troughs superimposed on the SE ridge flank are partly correlated, partly anti-correlated with the Bouguer anomaly and suggest that variable crustal density variations accompany the morphology variations.  相似文献   

16.
The study on topsoil contamination due to heavy metals was carried out by using the Magnetic susceptibility (MS) measurements in Izmit industrial city, northern Turkey. We attempted to investigate correlations between the concentration of selected heavy metals and the MS from 41 sample sites around Izmit Gulf. These investigations let us quantify and standardize the MS method, which may have consequences for long term monitoring of anthropogenic pollution, especially in urban areas. The MS surfer contour map based on the topsoil measurements was compiled with a randomly ranged distance density. The soil samples collected throughout the industrial areas, the parks, road sides and residential areas were also analyzed by Atomic Absorption Spectrometer. Heavy metals Cu, Ni, Cr and Pb show strong correlations with MS, while Zn and Co show a weak correlation with MS. Moreover, the Tomlinson pollution load index (PLI) shows insignificant correlation with the MS.The MS was examined vertically (0–30 cm) with respect to anthropogenic and/or lithogenic influences at the fourteen sample sites. The maximum values were mostly observed in depths of 2–5 cm and the MS values on the depth profiles vary between 10 × 10? 8 m3 kg? 1 and 203 × 10? 8 m3 kg? 1. The study revealed that MS is an inexpensive, fast and non-destructive method for the detection and mapping of contaminated soils.  相似文献   

17.
Two identical planar blunt probes of stainless steel material, biased with a bipolar ramp voltage, are used to measure the stratospheric polar conductivities to altitudes of 34 km. One probe (DP) is mounted closer to the gondola, looking downwards and shielded from sunlight, while the other (SP) is mounted looking sideways, away from the gondola. The daytime observations of positive ions in the 29–34 km altitude range with SP, and of negative ions at 34 km with DP, show photoelectric contaminations induced by solar UV radiations in the 190–230 nm band. These contaminations are found to be due to photoemissions from the SP probe steel surface and from the carbon paint that coats the surface of the gondola, respectively. It is found that, a segment of the photocurrent contaminated IV curve, recorded with SP at higher negative probe potentials, is linear, and it can give the ambient positive polar conductivity.  相似文献   

18.
Observational data of the seismic related VHF radio emissions at 41 and 53 MHz obtained at the four stations of Create Island are presented. The epicenter of EQs is located at the distance more than 300 km behind the horizon. It was shown that VHF radiation is generated at the altitudes 1–10 km in the atmosphere over the epicenter of EQs. The theory of generation of electromagnetic radiation by random electric discharges was developed. These discharges are excited by DC electric field enhanced up to the breakdown value in the atmosphere. The field is connected with the electric current flowing in the atmosphere–ionosphere circuit, whose source is generated by convective transport of charged aerosols, which are injected in the atmosphere by soil gases during the enhancement of seismicity. Calculations of the spectrum of electromagnetic radiation are derived, and the theoretical results are confirmed by observation data.  相似文献   

19.
《Marine pollution bulletin》2009,58(6-12):883-888
Mussels (Mytilus galloprovincialis) were transplanted to seven stations around a large shipyard for 126 days to evaluate tributyltin (TBT) contamination. Although the application of TBT-based paints to ships is totally banned in Korea, butyltin compounds were found to accumulate in mussels following transplantation. Concentrations of TBT and total butyltins in transplanted mussels near the shipyard were in the range of 40–350 ng Sn/g and 74–530 ng Sn/g on a dry weight basis, respectively. Obviously, low TBT concentrations (6.0–53 ng Sn/g dw) were determined in mussels at four stations outside the shipyard. A negative gradient of TBT concentrations and TBT portion to total butyltin concentrations were found in both the surface water and transplanted mussels according to distance from the shipyard. In addition, TBT concentrations in surface water and transplanted mussels showed significant correlation (r2 = 0.71; p < 0.001). These results indicate that the shipyard still releases fresh TBT into surrounding waters even after TBT regulation in Korea, and mussel transplantation is useful in evaluating TBT contamination in shipyard area.  相似文献   

20.
Parallel factor analysis of fluorescence excitation emission matrices of surface water samples of a globally large river (Yangtze River, China) watershed identified three classes of fluorescent dissolved organic matter (FDOM) that had ex/em = 280/330 nm, 305/385 nm and 350/450 nm respectively, resembling “peak T”, “peak M” and “peak C” commonly identified in natural water, respectively. Peak T (a tyrosine/tryptophan-like FDOM) did not show correlations to peak M or C which were humic-like substances, while a positive correlation (r = 0.935, p < 0.001) was present between the natural log-transformed maximum fluorescence intensity (Fmax) of peaks T and M indicating a tight link during their production and processing. Fmax values (in Raman unit nm?1) normalized to dissolved organic carbon (DOC) concentration were low, varying in ranges 15.93–85.95, 29.83–83.54 and 19.73–51.05 × 10?5 nm?1 (μmol/L)?1 for peaks T, M and C, respectively, in line with the history of strong photobleaching of the water samples as indicated by fairly high absorption spectral slope ratios (0.75–1.53 with a mean 1.03). Intermediate fluorescence index (FI) (1.46–1.83 with a mean 1.61) and small specific absorption at 254 nm (0.64–1.93 with a mean 1.15 m?1 mg?1 L) of the water samples, indicated the presence of both aquatic microbial DOM (e.g. peak T) and soil DOM (e.g. peak C). Peak C could be substantially removed by UV-A (320–400 nm) irradiation, while peak M was slightly increased when a microbe-containing water was exposed to the same UV-A irradiation. Taken together, peak C was attributed to diffuse soil source while peak M was likely attributed to joint effects of microbial activities and solar irradiation on the chromophores in the sample.  相似文献   

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

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