首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 677 毫秒
1.
The critical and transition conditions of a gaseous explosion at constant pressure and varying pressure is investigated. The study covered the solutions of the problem in pressure-time plane, temperature-time plane, as well as pressure-temperature plane. The explosion characteristics of combustible gases using the definition of criticality as an inflection point in critical trajectory. as well as other definitions of criticality are presented. The analytical study showed that the critical parameters for constant pressure explosion are lower than those for constant volumes where θ1tr is lowered from 0.5 to 0.29 and θatr is also lowered from 0.25 to 0.17. The sensitivity of the critical pressure to the ambient temperature and the boundary between ignition and non-ignition regions are also offered. It was also found that the critical parameters in the temperature–pressure plane were different from those obtained in the temperature-time plane or the pressure–time plane.  相似文献   

2.
Analytical and numerical solutions are used to determine the critical conditions for thermal explosion of autocatalytic reaction. The solutions covers both the reaction governed by the Arrhenius kinetics equation and the Frank-Kamenetskii approximation for that equation. The definition of criticality as the point at which d2θ/2=0, d3θ/3=0 and /〉0 is used here. The study is dealt with low and high exothermicity (B) of the reaction and their effects on the critical parameters. The numerical solutions cover the whole reaction from start at β = 0 up to the end at β = 1.0. All trajectories from subcritical, critical to supercritical are offered. The effects of different parameters such as B, ψ and θa (ambient temperature) on the critical conditions are presented. The results showed that the lower the autocatalytic factor (β0) is, the pronounced autocatalytic reaction explosion. The analytical solution offered analytical expressions for the critical condition and the different limits of the solutions are clarified. It was found that the numerical results confirm the analytical solution.  相似文献   

3.
The coupling of gas explosion flame and shock wave is analyzed. In the gas explosion process, shock wave is affected by the flame directly, and shock wave also induces the flame. Inhibiting explosion can be achieved by the interference between the flame and shock wave propagation. If the coupling effects can be damaged, the adverse effects caused by the explosion should be mitigated and controlled. According to the structure characteristics of foam ceramics, the coupling effects mechanism of ceramic foam on gas explosion flame and shock wave is researched. When the explosion goes through the structure of foam ceramics, the flame can be quenched and the shock wave be attenuated. After the flame is quenched, the supply of precursor shock wave energy is cut off. Due to lack of energy supply, the destructive effects of blast wave will be reduced effectively. Coupling effects of the flame and shock wave can be damaged by the special structure of foam ceramics. Studies suggest that a certain function to represent the structure characteristics of foam ceramics must exist. For a certain material of foam ceramics, the sure porosity δ and the pore diameter d also can be get, which is the key to research and develop foam ceramic suppression technology of gas explosion.  相似文献   

4.
Flameproof enclosures having internal electrical components are generally used in classified hazardous areas such as underground coalmines, refineries and places where explosive gas atmosphere may be formed. Flameproof enclosure can withstand the pressure developed during an internal explosion of an explosive mixture due to electrical arc, spark or hot surface of internal electrical components. The internal electrical component of a flameproof enclosure can form ignition source and also work as an obstacle in the explosion wave propagation. The ignition source position and obstacle in a flameproof enclosure have significant effect on explosion pressure development and rate of explosion pressure rise. To study this effect three cylindrical flameproof enclosures with different diameters and heights are chosen to perform the experiment. The explosive mixture used for the experiment is stoichiometric composition of methane in air at normal atmospheric pressure and temperature.It is observed that the development of maximum explosion pressure (Pmax) and maximum rate of explosion pressure rise (dp/dt)ex in a cylindrical flameproof enclosure are influenced by the position of ignition source, presence of internal metal or non-metal obstacles (component). The severity index, KG is also calculated for the cylindrical enclosures and found that it is influenced by position of ignition source as well as blockage ratios (BR) of the obstacles in the enclosures.  相似文献   

5.
Explosions will, in most cases, generate blast waves. While simple models (e.g., Multi Energy Method) are useful for simple explosion geometries, most practical explosions are far from trivial and require detailed analyses. For a reliable estimate of the blast from a gas explosion it is necessary to know the explosion strength. The source explosion may not be symmetric; the pressure waves will be reflected or deflected when hitting objects, or even worse, the blast waves may propagate inside buildings or tunnels with a very low rate of decay. The use of computational fluid dynamics (CFD) explosion models for near and far field blast wave predictions has many advantages. These include more precise estimates of the energy and resulting pressure of the blast wave, as well as the ability to evaluate non-symmetrical effects caused by realistic geometries, gas cloud variations and ignition locations. This is essential when evaluating the likelihood of a given leak source as cause of an explosion or equally when evaluating the potential risk associated with a given leak source for a consequence analysis.In addition, unlike simple methods, CFD explosion models can also evaluate detailed dynamic effects in the near and far field, which include time dependent pressure loads as well as reflection and focusing of the blast waves. This is particularly valuable when assessing actual near-field blast damage during an explosion investigation or potential near-field damage during a risk analysis for a facility. One main challenge in applying CFD, however, is that these models require more information about the actual facility, including geometry details and process information. Collecting the necessary geometry and process data may be quite time consuming. This paper will show some blast prediction validation examples for the CFD model FLACS. It will also provide examples of how directional effects or interaction with objects can significantly influence the dynamics of the blast wave. Finally, the challenge of obtaining useful predictions with insufficient details regarding the geometry will also be addressed.  相似文献   

6.
The number of explosive attacks on civilian buildings has recently increased and the pattern of damage inflicted on structures when an explosion takes place at altitude remains quite difficult to predict. The primary aim of the work reported here was to enhance the understanding of how blast waves from an explosion at altitude interact with the ground and with a structure. Small-scale experiments were conducted using a propane–oxygen stoichiometric mixture as explosive. This approach is original because it models high-explosive detonation in terms of gaseous charge explosion using TNT equivalents. Several non-dimensional laws are expressed and validated by experiments. These relationships allow determination of the propagation of a blast wave and its interaction with a structure as a function of the position of the explosive charge when the explosion occurs at altitude. Then, from knowledge of the blast loading, using Hopkinson's scaling law and TNT equivalents, we can predict the interaction of blast waves with the ground and a structure on a real scale. Simulations were performed using the Autodyn code, and good correlation with the experimental results was obtained.  相似文献   

7.
The dust explosion behaviors induced by two different combustion mechanisms (homogeneous and heterogeneous mechanisms) were comparatively investigated, based on the experiments under different dust concentrations, particle sizes and initial pressures in Siwek 20-L chamber. Based on the thermo-gravimetric analysis (TGA), sweet potato dust and magnesium dust were selected as the representative dusts with homogeneous and heterogeneous combustion mechanisms, respectively. Experiments find that these two dusts have different behaviors in the explosion kinetics due to different combustion mechanisms. For sweet potato dust, the explosion pressure pmax, the pressure rise rate (dp/dt)max and the combustion fraction η exhibit similar variation trends as dust concentration increases and they all reach to the maximum values at the worst-case concentration; while for magnesium dust, the variation of (dp/dt)max is somewhat different from that of pmax, that is, the (dp/dt)max will achieve the maximum at the concentration higher than the worst-case and keep stabilized with further increase of dust concentration. As the particle size decreases, the (dp/dt)max for sweet potato dust will increasingly rise and gradually approach to a stabilized value, but for magnesium dust, the increase of (dp/dt)max becomes pronounced only in the range of smaller particle sizes. To account the effect of initial pressure on pmax under different combustion mechanisms, a dimensionless pressure PR was introduced to denote the relative intensity of explosion. It is found that, for sweet potato dust, the increased initial pressure will promote the explosion process (or with high PR) for the dust cloud with high concentration due to the augmented oxygen concentration, but for the dust cloud with low concentration, the increased initial pressure will suppress the explosion process due to the increased resistance in devolatilization. For magnesium dust, the rise of initial pressure will generally promote the explosion process even for the dust cloud with low concentration; however, in the case of small particle size, the promotion of increased initial pressure to the explosion process is not so pronounced.  相似文献   

8.
瓦斯对煤尘爆炸特性影响的实验研究   总被引:2,自引:3,他引:2  
瓦斯的存在对煤尘爆炸特性的理论计算和数值仿真的结果与实际数据有一定差距,因此,通过不同浓度瓦斯与煤尘共存条件下爆炸实验研究,得出了矿井瓦斯对煤尘的最低着火温度、最小点火能量、爆炸下限浓度、最大爆炸压力和最大爆炸压力上升速度等爆炸特性影响的规律即瓦斯对煤尘最低着火温度影响不大;瓦斯可使煤尘的最小点火能量减小,尤其是对难于点燃的煤尘;混合物的爆炸下限浓度随瓦斯浓度的增加而降低;混合物的最大爆炸压力上升速度由于瓦斯的存在而增强,而最大爆炸压力几乎没有变化。同时研究了瓦斯对无爆炸性煤尘的影响。实验研究的结论对于现场防止煤尘爆炸的发生具有指导意义。  相似文献   

9.
To further elucidate the influence mechanism of side vents on the dynamic characteristics of gas explosions in tubes is helpful to design more reasonable vent layouts. In this paper, 9.5% methane-air explosion experiments were conducted in a tube with two side-vented ducts, and the effects of vent layouts and vent areas on the dynamic characteristics of explosion overpressure and flame propagation speed were investigated. The results demonstrate that under the same condition with a single vent area of 100 mm × 100 mm, when only the end vent is open, the maximum explosion overpressure and the maximum flame propagation speed are the highest among the five vent layouts. When the side vents 1 and 2 and the end vent are open, the maximum explosion overpressure is the lowest, and an unusual discovery is that the flame front changes into a hemispherical shape, finger shape, quasi-plane shape, tulip shape and wrinkled structure. When only side vent 1 is open, a unique Helmholtz oscillation occurs, and a new discovery is that there is a consistent oscillation relationship among the overpressure, flame propagation speed and flame structure. Helmholtz oscillation occurs only when a single vent area is 100 mm × 100 mm–60 mm × 60 mm, and the oscillation degree decreases with decreasing vent area. During the vent failure stage, the maximum explosion overpressure is generated, the flame front begins to appear irregular shape, and the flame propagation speed shows a prominent characteristic peak. After the vent failure stage, the driving effect of the end vent on the flame is higher than that of the side vent on the flame. Furthermore, the correlation equations of the mathematical relationships among the maximum explosion overpressure Pred, the static activation pressure Pstat and the vent coefficient Kv under four vent layouts are established, respectively.  相似文献   

10.
Elongated congestion patterns are common at chemical processing and petroleum refining facilities due to the arrangement of processing units. The accidental vapor cloud explosion (VCE) which occurred at the Buncefield, UK facility involved an elongated congested volume formed by the trees and undergrowth along the site boundary. Although elongated congested volumes are common, there have been few evaluations reported for the blast loads produced by elongated VCEs. Standard VCE blast load prediction techniques do not directly consider the impact of this congested volume geometry versus a more compact geometry.This paper discusses an evaluation performed to characterize the blast loads from elongated VCEs and to identify some significant differences in the resulting blast wave shape versus those predicted by well-known VCE blast load methodologies (e.g., BST and TNO MEM). The standard blast curves are based on an assumption that the portion of the flammable gas cloud participating in the VCE is hemispherical and located at grade level. The results of this evaluation showed that the blast wave shape for an elongated VCE in the near-field along the long-axis direction is similar to that for an acoustic wave generated in hemispherical VCEs with a low flame speed. Like an acoustic wave, an elongated VCE blast wave has a very quick transition from the positive phase peak pressure to the negative phase peak pressure, relative to the positive phase duration. The magnitude of the applied negative pressure on a building face depends strongly on the transition time between the positive and negative phase peak pressures, and this applied negative phase can be important to structural response under certain conditions. The main purpose of this evaluation was to extend previous work in order to investigate how an elongated VCE geometry impacts the resultant blast wave shape in the near-field. The influence of the normalized flame travel distance and the flame speed on the blast wave shape was examined. Deflagration and deflagration-to-detonation transition regimes were also identified for unconfined elongated VCEs as a function of the normalized flame travel distance and flame speed attained at a specified flame travel distance.  相似文献   

11.
To forestall, control, and mitigate the detrimental effects of aluminium dust, a 20-L near-spherical dust explosion experimental system and an HY16429 type dust-cloud ignition temperature test device were employed to explore the explosion characteristics of micron-sized aluminium powder under different ignition energies, dust particle sizes, and dust cloud concentration (Cdust) values; the minimum ignition temperature (MIT) values of aluminium powder under different dust particle sizes and Cdust were also examined. Flame images at different times were photographed by a high-speed camera. Results revealed that under similar dust-cloud concentrations and with dust particle size increasing from 42.89 to 141.70 μm, the MIT of aluminium powder increased. Under various Cdust values, the MIT of aluminium dust clouds attained peak value when concentrations enhanced. Furthermore, the increase of ignition energy contributed to the increase of the explosion pressure (Pex) and the rate of explosion pressure rise [(dP/dt)ex]. When dust particle size was augmented gradually, the Pex and (dP/dt)ex attenuated. Decreasing particle size lowered both the most violent explosion concentration and explosive limits.  相似文献   

12.
The method of explosion venting is widely used in industrial explosion-proof design due to its simple operation, economical and practical features. A dump vessel vented platform was built. By changing the vacuum level and the gas in the dump vessels and the structural size of linked vessels, the pressure in the explosion vessel and the dump vessel was compared, and the influencing factors of explosion venting investigated. The main conclusions are as follows: In the explosion venting process, the higher the vacuum in the dump vessel, the smaller the pressure peak of the explosion vessel and the dump vessel, and the faster the explosion pressure is lowered. When the dump vessel is under the same vacuum level and the gas in the dump vessel is CO2, the maximum pressure of the explosion vessel and the dump vessel is less than the maximum pressure when the containment medium is air. Under the same vacuum condition, the larger the volume ratio of the dump vessel and the explosion vessel, the smaller the pressure peak of the explosion vessel, the faster the explosion pressure drops, and the volume of the dump vessel reaches or exceeds the explosion vessel. Increasing the volume ratio of the containment vessel to the explosion vessel facilitates protection of the explosion vessel and the containment vessel. Under the same vacuum condition, when the gas explosion in 113 L vessel vents into 22 L vessel, the longer the length of the pipe, the greater the maximum pressure in the spherical vessel. When the gas explosion in 22 L vessel vents into 113 L dump vessel, as the pipeline grows, the maximum pressure in the two vessels decreases, but the reduction is not significant. In practical application, it is recommended to use a vacuum of 0.08Mpa or more for the dump vessel vented, and the containment medium is CO2.In terms of the structural size of the container, it is recommended that the ratio of the receiving container to the explosion container be as large as possible, and the pipe length be as long.  相似文献   

13.
A gas explosion in an underground structure may cause serious damage to the human body and ground buildings and may result in huge economic losses. The pressure of the gas explosion is an important parameter in determining its severity and designating an emergency plan. However, existing empirical and computational fluid dynamics (CFD) methods for pressure prediction are either inaccurate or inefficient when considering multiple influencing factors and their interrelationships. Therefore, for a more efficient and reliable prediction, the present study developed a multifactorial prediction model based on a beetle antennae search (BAS) algorithm improved back propagation (BP) neural network. A total of 317 sets of data which considered factors of geometry, gas, obstacle, vent, and ignition were collected from previous studies. The results showed that the established model can predict pressures accurately by low RMSE (43.4542 and 50.7176) and MAPE (3.9666% and 4.9605%) values and high R2 (0.7696 and 0.7388) values for training and testing datasets, respectively. Meanwhile, the BAS algorithm was applied to improve both the calculation efficiency and the accuracy of the proposed model by enabling a more intelligent hyperparameter tuning method. Furthermore, the permutation importance of input variables was investigated, and the length (L) and the ratio of length and diameter (L/D) of geometry were found to be the most critical factors that affect the explosion pressure level.  相似文献   

14.
The explosion characteristics of anthracite coal dust with/without small amount of CH4 (1.14 vol %) were investigated by using a 20 L spherical explosion apparatus with an emphasis on the roles of oxygen mole fraction and inert gas. Two methods based on overpressure and combustion duration time were used to determine the minimum explosion concentration (MEC) or the lower explosion limit (LEL) of the pure anthracite coal dust and the hybrid coal-methane mixtures, respectively. The experiment results showed that increasing oxygen mole fraction increases the explosion risk of coal dust: with increasing oxygen mole fraction, the explosion pressure (Pex) and the rate of explosion pressure rise ((dp/dt)ex)) increase, while MEC decreases. The explosion risk of anthracite dust was found to be lower after replacing N2 with CO2, suggesting that CO2 has a better inhibition effect on explosion mainly due to its higher specific heat. However, the addition of 1.14% CH4 moderates the inhibition effect of CO2 and the promotion effect of O2 on anthracite dust explosion for some extent, increasing explosion severity and reducing the MEC of anthracite dust. For hybrid anthracite/CH4 mixture explosions, Barknecht's curve was found to be more accurate and conservative than Chatelier's line, but neither are sufficient from the safety considerations. The experimental results provide a certain help for the explosion prevention and suppression in carbonaceous dust industries.  相似文献   

15.
The effect of pyrolysis and oxidation characteristics on the explosion sensitivity and severity parameters, including the minimum ignition energy MIE, minimum ignition temperature MIT, minimum explosion concentration MEC, maximum explosion pressure Pmax, maximum rate of pressure rise (dP/dt)max and deflagration index Kst, of lauric acid and stearic acid dust clouds was experimentally investigated. A synchronous thermal analyser was used to test the particle thermal characteristics. The functional test apparatuses including the 1.2 L Hartmann-tube apparatus, modified Godbert-Greenwald furnace, and 20 L explosion apparatus were used to test the explosion parameters. The results indicated that the rapid and slow weight loss processes of lauric acid dust followed a one-dimensional diffusion model (D1 model) and a 1.5 order chemical reaction model (F1.5 model), respectively. In addition, the rapid and slow weight loss processes of stearic acid followed a 1.5 order chemical reaction model (F1.5 model) and a three-dimensional diffusion model (D3 model), respectively, and the corresponding average apparent activation energy E and pre-exponential factor A were larger than those of lauric acid. The stearic acid dust explosion had higher values of MIE and MIT, which were mainly dependent on the higher pyrolysis and oxidation temperatures and the larger apparent activation energy E determining the slower rate of chemical bond breakage during pyrolysis and oxidation. In contrast, the lauric acid dust explosion had a higher MEC related to a smaller pre-exponential factor A with a lower amount of released reaction heat and a lower heat release rate during pyrolysis and oxidation. Additionally, due to the competition regime of the higher oxidation reaction heat release and greater consumption of oxygen during explosion, the explosion pressure Pm of the stearic acid dust was larger in low concentration ranges and decayed to an even smaller pressure than with lauric acid when the concentration exceeded 500 g/m3. The rate of explosion pressure rise (dP/dt)m of the stearic acid dust was always larger in the experimental concentration range. The stearic acid dust explosion possessed a higher Pmax, (dP/dt)max and Kst mainly because of a larger pre-exponential factor A related to more active sites participating in the pyrolysis and oxidation reaction. Consequently, the active chemical reaction occurred more violently, and the temperature and overpressure rose faster, indicating a higher explosion hazard class for stearic acid dust.  相似文献   

16.
To reveal the effects of different inert gases on explosion characteristics during low density polyethylene (LDPE) dust explosion and optimize the explosion-proof process, eight N2 (CO2)/air mixed inerting conditions were experimentally studied. Typical inerting conditions with 12 L cylindrical explosive tank were used to study the characteristics on the flame propagation. The thermogravimetric analysis with related theories were used to further explain the mechanism and quantities in low density polyethylene (LDPE) dust explosion with different inert gases. The results showed that the reduction of O2 concentration could effectively delay the progress of flame growth process and weaken the effect of dust combustion reaction. The flame growth process of condition (N2/air (18% O2)) was 2.05 times slower than that of the non-inert condition. The explosion strength was obviously reduced, and the characteristic parameters such as explosion pressure and flame propagation speed were also affected by the decrease of O2 concentration. For LDPE powder, the smaller the median diameter, the greater the explosion intensity and the lower the limiting oxygen content (LOC). The LOC with CO2 was usually higher than that with N2 and the effect of CO2 was significantly better than N2 in inerting.  相似文献   

17.
Explosion characteristics of five alcohol–air (ethanol, 1-butanol, 1-pentanol, 2-pentanol and 3-pentanol) mixtures were experimentally conducted in an isochoric chamber over wide ranges of initial temperature and pressure. The effect of temperature and pressure on the different explosion behaviors among these alcohols with various structures were investigated. Results show that the peak explosion pressure is increased with the decrease of temperature and increase of pressure. Maximum rate of pressure rise is insensitive to the temperature variation while it significantly increases with the increase of initial pressure. Among the 1-, 2-, and 3-pentanol–air mixtures, 1-pentanol has the highest values in peak explosion pressure and maximum rate of pressure rise and 2-pentanol gives the lowest values at the initial pressure of 0.1 MPa. These differences tend to be decreased with the increase of initial pressure. Among the three primary alcohol–air (ethanol, 1-butanol and 1-pentanol) mixtures, a similar explosion behavior is presented at the lean mixture side because of the combined effect of adiabtic temperature and flame propagation speed. At the rich mixture side, 1-pentanol gives the highest values in peak explosion pressure and maximum rate of pressure rise and ethanol gives the lowest values. This phenomenon can be interpretated from the combining influence of heat release and heat loss, since the flame speeds of ethanol-, 1-butanol-, 1-pentanolair mixtures are close at rich mixture side.  相似文献   

18.
The inhibition effect of heptafluoropropane (CF3CHFCF3) on methane explosions under different inhibitor concentrations in a closed vessel was studied. A high-speed camera and a pressure sensor were adopted respectively to record flame propagation characteristics and pressure data. Results indicate that the relationship between flame propagation and pressure rising was correlated. As the equivalent ratio (ϕ)≤1, the pressure presented a trend of rising firstly and then decreasing with increasing CF3CHFCF3 concentration, and it was found that there existed a critical concentration for pressure decrease. As ϕ > 1, the pressure exhibited a decreasing trend. Although the pressure appeared to seemingly increase, the moment that the pressure began to rise (trise) and the moment that the maximum explosion overpressure appeared (tPmax) were obviously delayed. The average rate of pressure rise ((dP/dt)ave) was decreased as the concentration of CF3CHFCF3 increased. It indicates that CF3CHFCF3 can effectively reduce the explosion reaction rate. The critical concentration of CF3CHFCF3 for complete inhibition was determined. Meanwhile, the synergy of CF3CHFCF3-inert gas can improve the inhibition effect. Compared with CF3CHFCF3–N2, the synergy of CF3CHFCF3–CO2 presented a better inhibition effect, and the inhibition effect was increased with increasing inert gas concentration. And the mechanisms of physical and chemical effects on explosion inhibition were analyzed.  相似文献   

19.
Experiments were performed on the influence of pre-ignition turbulence on the course of vented gas and dust explosions. A vertical cylindrical explosion chamber of approximately 100 l volume and a length-to-diameter ratio (l/d) of 4.7 consisting of a steel bottom segment and three glass sections connected by steel flanges was used to perform the experiments. Sixteen small fans evenly distributed within the chamber produced turbulent fluctuations from 0 to 0.45 m/s. A Laser-Doppler-anemometer (LDA) was used to measure the flow and turbulence fields. During the experiments the pressure and in the case of dust explosions the dust concentration were measured. In addition, the flame propagation was observed by a high-speed video camera. A propane/nitrogen/oxygen mixture was used for the gas explosion experiments, while the dust explosions were produced by a cornstarch/air mixture.It turned out that the reduced explosion pressure increased with increasing turbulence intensity. This effect was most pronounced for small vents with low activation pressures, e.g. for bursting disks made from polyethylene foil. In this case, the overpressure at an initial turbulence of 0.45 m/s was twice that for zero initial turbulence.  相似文献   

20.
为研究障碍物截面变化对混合气体爆炸特性参数的影响规律,运用3D动态流体动力学软件AutoReaGas建立障碍物截面不同变化方式和趋势的爆炸模型,模拟分析截面不同变化方式(渐变、突变)和不同变化趋势(由大变小、由小变大)对气体爆炸特性参数(速度、压力场)的影响程度和规律。研究表明:改变障碍物的截面变化方式和趋势都会影响可燃气体的爆炸参数。对于爆炸压力场,截面变化情况下的爆炸超压大于截面不变的情况,截面渐变大于截面突变的情况,截面由大到小变化与截面由小到大变化下的爆炸超压呈现相反变化趋势。对于爆炸速度场,截面不同变化方式和趋势下的速度变化情况均为在开始的一小段时间内先减小然后逐渐增大。  相似文献   

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

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