首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The study of coda waves has recently attracted increasing attention from seismologists. This is due to the fact that it is viewed as a new means by which the stress accumulation stage preceding a large earthquake can be measured, since the scattering paths nearly uniformly cover a fairly large region around the focus and observation stations, compared with the direct ray paths. To date, we have had many reports on the temporal variation of the relation between coda duration and amplitude magnitude, and that of the coda attenuationQ c –1 which is estimated from coda amplitude decay. Some of these have shown a precursor-like behavior; however, others seem to have shown a coseismic change. We have critically reviewed these reports, and discussed what these observational facts tell us about the change in the heterogeneous crust. We found significant temporal variations, not only in the mean but also in the scatter ofQ c –1 , associated with the mainshock occurrence. The formation of new cracks, the reopening and growing of existing cracks, the interaction of these cracks, and the pore water movement through these cracks might correspond to such variations. In addition, we may expect an inhomogeneous distribution of crack clusters in a fairly large region, compared with the aftershock region. The gradual appearance of such crack clusters seems to be the most plausible mechanism by which coda decay gradients are caused to largely scatter in the stress accumulation stage.  相似文献   

2.
Double-torsion specimens of two granitic rocks were prepared in several directions with reference to microcracks fabric. Even for the same rock and at the same stress levels, the observed crack velocities in two granitic rocks were dependent on both the propagation direction and the opening direction. The maximum difference by several orders of magnitude was found for both rocks. The highest crack velocity was observed when the subcritical crack was parallel to most of the preexisting cracks. The maximum critical stress intensity factor was about twice as high as the minimum one in different directions. An analysis for a thin plate having anisotropic elasticity under torsional load showed that the observed difference in the crack velocity and the critical stress intensity factor was not an error due to conventional equations derived on the assumption of isotropic elasticity but the true material's property. As the preferred orientation of microcracks has been pointed out for many granitic rocks, we can conclude that the anisotropic nature of the fracture resistance of the two granitic rocks used in this study was not exceptional. A region of a transport-limited velocity was not found for rocks, even at the velocity of 10–2 m/s, that is almost equal to the theoretical limit of the stress corrosion cracking.  相似文献   

3.
Summary The problem of expressing analytically the magnetic torque, acting on the electrically conducting part of the Earth's mantle, is treated as a function of the system of convection on the surface of the core. The changes of velocities in the system of convection are estimated for decadic changes of the Earth's rotation and for the perturbation of the Earth's rotation in 1897. As regards the decadic changes of the Earth's rotation a change of velocity in the system of convection at the surface of the core of the order of 10–4 m/s corresponds, and as regards the perturbation of the Earth's rotation in 1897 (10–3 s/year) a change of velocity of 10–3 m/s reduced to the whole surface of the core corresponds, and 10–2 m/s corresponds for the region of the focus of the world geomagnetic anomaly (dimension of this region is 106 m).  相似文献   

4.
Fracture and frictional sliding are considered as phenomena involving brittle failure. Brittle failure is preceded by the formation of small (subcritical) cracks. In non-water-saturated rock, the distribution, shape and size of these suberitical cracks determine the change in the physical properties prior to failure. A model is proposed which suggests that the spatial and temporal distribution, shape and size of subcritical cracks within a stressed rock depend upon the rate of deformation and the volatile content.As a rock is stressed beyond about 50 percent of its ultimate failure stress, dilatancy is initiated. With increasing stress a broad zone of cracks develops within the dilatant region. The seismic velocities through this zone decrease markedly and the cracks grow more numerous., changing in size and shape. Before brittle failure of the rock occurs, the subcritical cracks interact, leading to a concentration of the zone. During the stage when the zone narrows, the seismic velocities in crease in the surrounding volume due to local rotation of stresses and consequent closure of some cracks. In most laboratory experiments the stage during which the velocity increases and the now intense deformation zone becomes narrow is very short and difficult to observe experimentally. At very low strain rates and with volatiles present, the crack growth and subsequent interaction lead to the narrowing of the intense deformation zone and therefore to an observable increase in velocity.The above is based upon an interpretation of a number of experiments. Using optical holography we have observed the development and subsequent intensification of a deformation zone. Ultrasonic velocity measurements showed a distinct anomaly (decrease followed by an increase) before failure. The anomaly was only detectable at our lowest experimental strain rates (3×10–8/sec).  相似文献   

5.
In a series of experiments at 0.5–1.3 GPa and 1050–1200°C we have monitored the transport, via crack propagation, of CO2 into well-annealed olivine and quartz aggregates. The objectives were to determine (1) the extent and rate of fluid penetration; (2) the effect of varying both P-T conditions and microstructure; and (3) the fluid penetration pathways. Experiments on CO2 penetration into dunite annealed in the absence of MgO indicate rapid and pervasive fluid transport on a grain-dimension scale, but a limited penetration distance ( 1 mm). Additional experiments on dunite annealed in the presence of MgO (either dispersed or present at both ends), however, resulted in CO2 penetration that was both pervasive on the scale of individual grains and almost always completely through the 5 mm long samples. The abundance of fine (10 μm) grains in the MgO-free dunite, in contrast to the much larger grain sizes of the samples annealed with MgO present, suggests the difference in fluid penetration behavior may arise because the strength variation in dunite scales with the grain size. Effects arising from changes in olivine point defect chemistry, however, are an additional possibility. The response of synthetic quartzite to CO2 overpressure is distinct from that of dunite: Quartzite experiences rapid and complete penetration of CO2, via a macroscopically visible system of transgranular fractures, over the range of P-T conditions investigated.The small amount of porosity ( 2–3%) present in most rock samples fabricated for this study, lacks three-dimensional connectivity, thus precluding any enhanced fluid penetration via porous flow. Pores could possibly enhance fluid penetration as the result of a small reduction in resistance to fracture, but the probable abundance of strength-controlling flaws in natural rocks is likely to produce similar behavior.The results of our experiments on olivine and olivine + MgO suggest that the transport of pressurized CO2 in very olivine-rich mantle environments will be pervasive on the scale of individual grains and its extent may be dependent on rock microstructure and/or crystal chemical effects. Such pervasive fluid transport, perhaps associated with magma decarbonation, may have interesting implications for both magma transport and local LREE enrichment of adjacent mantle wall-rock. The ease with which quartzite is penetrated by CO2 at the conditions of our experiments underscores the possible role of decarbonation reactions in crustal permeability-enhancement processes.  相似文献   

6.
This study presents the results of experimental compaction while measuring ultrasonic velocities of sands with different grain size, shape, sorting and mineralogy. Uniaxial mechanical compaction tests up to a maximum of 50 MPa effective stress were performed on 29 dry sand aggregates derived from eight different sands to measure the rock properties. A good agreement was found between the Gassmann saturated bulk moduli of dry and brine saturated tests of selected sands. Sand samples with poor sorting showed low initial porosity while sands with high grain angularity had high initial porosity. The sand compaction tests showed that at a given stress well‐sorted, coarse‐grained sands were more compressible and had higher velocities (Vp and Vs) than fine‐grained sands when the mineralogy was similar. This can be attributed to grain crushing, where coarser grains lead to high compressibility and large grain‐to‐grain contact areas result in high velocities. At medium to high stresses the angular coarse to medium grained sands (both sorted sands and un‐sorted whole sands) showed high compaction and velocities (Vp and Vs). The small grain‐to‐grain contact areas promote higher deformation at grain contacts, more crushing and increased porosity loss resulting in high velocities. Compaction and velocities (Vp and Vs) increased with decreasing sorting in sands. However, at the same porosity, the velocities in whole sands were slightly lower than in the well‐sorted sands indicating the presence of loose smaller grains in‐between the framework grains. Quartz‐poor sands (containing less than 55% quartz) showed higher velocities (Vp and Vs) compared to that of quartz‐rich sands. This could be the result of sintering and enlargement of grain contacts of ductile mineral grains in the quartz‐poor sands increasing the effective bulk and shear stiffness. Tests both from wet measurements and Gassmann brine substitution showed a decreasing Vp/Vs ratio with increasing effective stress. The quartz‐rich sands separated out towards the higher side of the Vp/Vs range. The Gassmann brine substituted Vp and Vs plotted against effective stress provide a measure of the expected velocity range to be found in these and similar sands during mechanical compaction. Deviations of actual well log data from experimental data may indicate uplift, the presence of hydrocarbon, overpressure and/or cementation. Data from this study may help to model velocity‐depth trends and to improve the characterization of reservoir sands from well log data in a low temperature (<80–100o C) zone where compaction of sands is mostly mechanical.  相似文献   

7.
The frictional properties of a crushed granite gouge and of gouges rich in montmorillonite, illite, and serpentine minerals have been investigated at temperatures as high as 600°C, confining pressures as high as 2.5 kbar, a pore pressure of 30 bar, and sliding velocities of 4.8 and 4.8×10–2 m/sec. The gouges showed nearly identical strength behaviors at the two sliding velocities; all four gouges, however, showed a greater tendency to stick-slip movement and somewhat higher stress drops in the experiments at 4.8×10–2 m/sec. Varying the sliding velocity also had an effect on the mineral assemblages and deformation textures developed in the heated gouges. The principal mineralogical difference was that at 400°C and 1 kbar confining pressure a serpentine breakdown reaction occurred in the experiments at 4.8×10–2 m/sec but not in those at 4.8 m/sec. The textures developed in the gouge layers were in part functions of the gouge type and the temperature, but changes in the sliding velocity affected, among other features, the degree of mineral deformation and the orientation of some fractures.  相似文献   

8.
Recently obtained data on oxygen diffusion in feldspars, quartz, and hornblende permit the prediction of the apparent18O16O temperatures that would be measured in a rock that consisted only of those three minerals, and cooled slowly from high temperature. The computed temperatures would be based on the differences in the18O16O ratios between coexisting pairs of minerals. The present calculation takes into account the diffusion rates for oxygen as a function of temperature, the cooling rate of the rock, the mineral grain sizes, and the mode of the rock. For mineral grains 1 mm in radius, and a cooling rate of 10°C/m.y., the minimum difference in apparent temperature between quartz-feldspar and feldspar-hornblende pairs will be 115°C, despite the assumption of a normal, uneventful, slow cooling history to room temperature. Further, the apparent quartz-hornblende temperature will range over 30°C (590–620°C) depending on the mode of the rock. For a cooling rate of 1000°C/m.y., the apparent difference in temperature can be as much as 400°C. Consequently, consistency in temperatures obtained by oxygen isotope analysis should not be expected in most high-grade metamorphic rocks or igneous rocks which are cooled slowly. Departures from the pattern of temperatures obtained in this model would imply a very rapid quench from high temperature, or a complex history for the rock. For some minerals, including hornblende, the relation between temperature and the equilibrium fractionation of oxygen isotopes between coexisting phases has been derived from observed relations in natural specimens. The choice of the specimens used for such calibrations needs to be re-evaluated in light of these findings. This may result in a change in the equilibrium equation constants.An example from the literature, the San Jose tonalite, Baja California, Mexico, was modelled and yieldsδ18O concentrations in the minerals that correspond closely with the measured values. This suggests that the model used is appropriate, that the rock has had a simple thermal history, and that it cooled at 100–200°C/m.y. over the temperature range 800–500°C. The set of paleotemperatures obtained for a rock will, in general, yield neither the mineral closure temperatures nor the formation or crystallization temperatures. On the other hand, the cooling rate of the rock may be derived from the data. This, in turn, may have important tectonic implications with regard to denudation and uplift rates.  相似文献   

9.
We theoretically study the scattering ofP, SV andSH waves by a zonal distribution of cracks, which simulates a fault fracture zone. An investigation is conducted how the geometrical properties of the crack distribution and the frictional characteristics of the crack surface are reflected in the attenuation and dispersion of incident waves, as well as in the amplitudes of the transmitted and reflected waves from the zone. If the crack distribution within the fault zone changes temporally during the preparation process of the expected earthquake, it will be important for earthquake prediction to monitor it, utilizing the scattering-induced wave phenomena.We consider the two-dimensional problem. Aligned cracks with the same length are assumed to be randomly distributed in a zone with a finite width, on which elastic waves are assumed to be incident. The distribution of cracks is assumed to be homogeneous and sparse. The crack surface is assumed to be stress-free, or to undergo viscous friction; the latter case simulates fluid-filled cracks. The opening displacement of the crack is assumed to be negligibly small. The idea of the mean wave formalism is employed in the analysis, and Foldy's approximation is assumed.When the crack surface is stress-free, it is commonly observed for every wave mode (P, SV andSH) that the attenuation coefficientQ –1 peaks aroundka1, the phase velocity is almost independent ofk in the rangeka<1 and it increases monotonically withk in the rangeka>1, wherek is the intrinsicS wavenumber anda is the half length of the crack. The effect of the friction is to shift the peak ofQ –1 and the corner of the phase velocity curve to the low wavenumber range. The high wavenumber asymptote ofQ –1 is proportional tok –1 independently of model parameters and the wave modes. If the seismological observation thatQ –1 ofS waves has a peak at around 0.5 Hz in the earth's crust is combined with our results, the upper limit of crack size within the crust is estimated about 4 km. The information regarding the transmitted and reflected waves, such as the high wavenumber limit of the amplitude of the transmitted wave etc., allows estimation of the strength of the friction.  相似文献   

10.
The velocities with which grains were observed to emerge from a sand bed after an intersaltation collision at u* = 40 cm s?1 are reported for four bed attitudes, from horizontal bed to adverse bed slope 15°. The principal effect of bed angle is to alter the magnitude and direction of the ricochet velocity. However, emergent velocities of dislodged grains are consistent with reptation path lengths comparable to the length of the upwind face of ripples in the corresponding wind. Calculations of the loss of forward momentum at collision, using the data for the range of bed attitudes studied suggest that creep is most vigorous on the sloping upwind face of the ripple and diminishes at the crest. As a result, the crest would be expected to accumulate the coarse material which moves predominantly by creep. The saltations originating in ricochet from the sloping back of the ripple are more vigorous and more concentrated in plan than are those originating at the crest. However, the saltation path length is at least an order of magnitude greater than the ripple wavelength and the probability distribution of path lengths is quite dispersed. Consequently it is very unlikely that these spatial patterns of ricochet are preserved sufficiently distinctly in the saltation cloud and subsequent collision distribution to be the agent of ripple development. This study therefore supports a view of moving grain interaction with the bed in which saltation provides the power to mobilize grains but ripple growth is associated with reptation and particularly with a pattern of impact which develops with the bed relief. Creep is more active on upwind facing slopes than at the crest, which therefore is a zone of net creep grain deposition.  相似文献   

11.
Load relaxation and cross-head displacement rate-change experiments have been used to establish log10 stress intensity factor (K) versus log10 crack velocity (v) diagrams for double torsion specimens, of synthetic quartz cracked on thea plane in liquid water and moist air.For crack propagation normal toz and normal tor at 20°C,K Ic (the critical stress intensity factor) was found to be 0.852±0.045 MN·m–3/2 and 1.002±0.048 MN·m–3/2, respectively.Subcritical crack growth at velocities from 10–3 m·s–1 to 10–9 m·s–1 at temperatures from 20°C to 80°C is believed to be facilitated by chemical reaction between the siloxane bonds of the quartz and the water or water vapour of the environment (stress corrosion). The slopes, of isotherms in theK-v diagrams are dependent upon crystallographic orientation. The isotherms have a slope of 12±0.6 for cracking normal tor and 19.9±1.7 for cracking normal toz. The activation enthalpy for crack propagation in the former orientation in liquid water at temperatures from 20°C to 80°C is 52.5±3.8 kJ·mole–1.A discussion is presented of the characteristics of theK-v diagrams for quartz.  相似文献   

12.
Pyroclastic flow emplacement is strongly influenced by eruption column height. A surface along which kinetic energy is zero theoretically connects the loci of eruption column collapse with all coeval ignimbrite termini. This surface is reconstructed as a two-dimensional energy line for the 1912 Katmai pyroclastic flow in the Valley of Ten Thousand Smokes from mapped flow termini and the runup of the ignimbrite onto obstructions and through passes. Extrapolation of the energy line to the vicinity of the source vent at Novarupta suggests the eruption column which generated the ignimbrite eruption was approximately 425 m high. The 1912 pyroclastic flow travelled about 25 km downvalley. Empirical velocity data calculated from runup elevations and surveyed centrifugal superelevations indicate initial velocities near Novarupta were greater than 79–88 m s–1. The flow progressively decelerated and was travelling only 2–8 m s–1 when it crossed a moraine 16 km downvalley. The constant slope of the energy line away from Novarupta suggests the flow was systematically slowed by internal and basal friction. Using a simple physical model to calculate flow velocities and a constant kinetic friction coefficient (Heim coefficient) of 0.04 derived from the reconstructed energy line, the flow is estimated to have decelerated at an average rate of –0.16 m s–2 and to have taken approximately 9.5 minutes to travel 25 km down the Valley of Ten Thousand Smokes. The shear strength of the flowing ignimbrite at the moraine was approximately 0.5 kPa, and its Bingham viscosity when it crossed the moraine was 3.5 × 103 P. If the flow was Newtonian, its viscosity was 4.2 × 103 P. Reynolds and Froude numbers at the moraine were only 41–62 and 0.84–1.04, respectively, indicating laminar, subcritical flow.  相似文献   

13.
The distribution of electric charge on the marine aerosol was determined near the sea surface of the Indian Ocean and the North Atlantic during the final stage of the Snellius II-expedition. Mean values for small ion concentrationsn +=455 cm–3 andn –=340 cm–3 were found over the Atlantic, whilen +=310 andn –=250 cm–3 were the mean values over the Indian Ocean. The ration +/n– increased from 1.2 to 2.0 with decreasing wind velocity. At wind velocities below 5 m/s 75% of the net space charge near the ocean surface was found to be carried by small ions.  相似文献   

14.
The systematic time differences observed in the onset of postsunset VHF scintillations recorded simultaneously at Ujjain (Geogr. lat. 23.2○N, Geogr. long. 75.6○E) and Bhopal (Geogr. lat. 23.2○N, Geogr. long. 77.6○E), situated at the peak of the anomaly crest in the Indian region, have been analysed to determine the zonal drifts of scintillation-producing irregularities. The method is based on the assumption that the horizontal movement of irregularities does not change while crossing the F-region cross-over points of these stations. The calculated velocities of irregularities indicate an eastward drift decreasing from about 180 m s−1 to 55 m s−1 during the course of night. In the premidnight period, the drifts are reduced under the magnetically disturbed conditions. The average east-west extension of irregularites is found to be in the range of 200–500 km.  相似文献   

15.
Experimental investigation of inertial flow processes in porous media   总被引：3，自引：0，他引：3  
  相似文献   

16.
The effect of randomly distributed cracks on the attenuation and dispersion ofSH waves is theoretically studied. If earthquake ruptures are caused by sudden coalescence of preexisting cracks, it will be crucial for earthquake prediction to monitor the temporal variation of the crack distribution. Our aim is to investigate how the property of crack distribution is reflected in the attenuation and dispersion of elastic waves.We introduce the stochastic property, in the mathematical analysis, for the distributions of crack location, crack size and crack orientation. The crack size distribution is assumed to be described by a power law probability density (p(a) a – fora minaa max according to recent seismological and experimental knowledge, wherea is a half crack length and the range 13 is assumed. The distribution of crack location is assumed to be homogeneous for the sake of mathematical simplicity, and a low crack density is assumed. The stochastic property of each crack is assumed to be independent of that of the other cracks. We assume two models, that is, the aligned crack model and the randomly oriented crack model, for the distribution of crack orientation. All cracks are assumed to be aligned in the former model. The orientation of each crack is assumed to be random in the latter model, and the homogeneous distribution is assumed for the crack orientation. The idea of the mean wave formalism is employed in the analysis, and Foldy's approximation is assumed.We observe the following features common to both the aligned crack model and the randomly oriented crack model. The attenuation coefficientQ –1 decays in proportion tok –1 in the high frequency range and its growth is proportional tok 2 in the low frequency range, wherek is the intrinsic wave number. This asymptotic behavior is parameter-independent, too. The attenuation coefficientQ –1 has a broader peak as increases and/ora min/a max decreases. The nondimensional peak wave numberk p a max at whichQ –1 takes the peak value is almost independent ofa min/a max for =1 and 2 while it considerably depends ona min/a max for =3. The phase velocity is almost independent ofk in the rangeka max<1 and increases monotonically ask increases in the rangeka max>1. While the magnitude ofQ –1 and the phase velocity considerably depend on the orientation of the crack in the aligned crack model, the above feature does not depend on the crack orientation.The accumulation of seismological measurements suggests thatQ –1 ofS waves has a peak at around 0.5 Hz. If this observation is combined with our theoretical results onk p a max, the probable range ofa max of the crack distribution in the earth can be estimated for =1 or 2. If we assume 4 km/sec as theS wave velocity of the matrix medium,a max is estimated to range from 2 to 5 km. We cannot estimatea max in a narrow range for =3.  相似文献   

17.
Analysis of the etch-pit size distributions (PSDs) observed on potassium feldspar and hornblende grains in a soil catena in loess (age = 12,500 y) reveals natural mineral etching rates. Rates estimated for hornblende (6 to 9×10–15mol/m2s) are based on consistent crystallographically controlled etch pits, while rates estimated for potassium feldspar (2×10–15mol/m2s) are based on irregularly shaped pits. Although little difference in etching rate is observed between soil horizons, the highest etching rates generally occur in the upper B horizons where pH values are lowest. Decreasing soil drainage correlates with an increase in pit density,n°, probably due to increased grain wetting, while decreased drainage correlates with a decrease in pit growth rate (G), probably due to increased dissolved solute concentrations. The PSD model predicts that etching rate is a function ofn° and ofG 4. Etching rates calculated for potassium feldspar do not vary with drainage, while those of hornblende decrease with decreasing drainage. Estimated etching rates are lower than bulk dissolution rates measured in the laboratory.  相似文献   

18.
Ground tilt measurements demonstrate that Askja is in a state of unrest, and that in the period 1988–1991 a maximum 48±3 rad tilt occurred down towards the centre of the caldera. This is consistent with 126 mm of deflation at the centre of the caldera with a 2.5–3.0 km depth to the source of deformation. The volume of the subsidence bowl is 6.2x106 m3. When combined with high precision microgravity measurements, the overall change in sub-surface mass may be quantified. After correction for the observed elevation change using the free air gradient of gravity measured for each station, the total decrease in mass is estimated to be less than 109 kg. A small residual ground inflation and net gravity increase in the southeastern part of the caldera may be caused by dyke intrusion in this region. The minimum dimensions of such an intrusion or complex of intrusions are 1 m width, up to 100 m deep and up to several hundred metres thick.  相似文献   

19.
A microprobe study of Ol–An nodules (allivalites) as sampled from the basaltic andesites of Malyi Semyachik Volcano in Kamchatka revealed small inclusions of a mafic mineral in crystals of olivine and anorthite and in the intergranular mass. The mineral was later identified as spinel. A microprobe analysis of the spinel showed that the grains that were found in minerals are not zonal, with spinel inclusions of different chemical compositions occurring in a single crystal. In contrast to this situation, the spinel crystals in the intergranular basis are zonal and exhibit a reaction rim. The spinel compositions form a long and well-pronounced trend that occurs in the solvus zone of solid solutions that has not been adequately described in the literature. The existence of the resulting spinel trend in the solvus zone can be explained by an early capture of spinel grains by growing silicate crystals and by their rapid cooling after the eruption, which chilled the metastable solution. Spinel grains occasionally exhibited initial phases in the decomposition of a solid solution owing to their long residence in a cooling lava flow at temperatures below the solvus. It is supposed that these spinels were formed somewhat earlier or simultaneously with the crystallization of the silicate phase. This diversity in the compositions of non-zonal spinel can be explained by thermo diffusion that led to homogeneous compositions in the initially zonal spinel crystals after they had been included in the host mineral crystals.  相似文献   

20.
An experimental study was carried out on a granitic mylonite (La Bresse, France) to analyze the influence of pore microstructure on transport properties. Different crack networks were obtained by a controlled thermal treatment. Microstructures were analyzed by means of gas adsorption and mercury porosimetry. Transport properties have been investigated by measuring gas permeability and electrical conductivity. The dependence of permeability on confining pressure shows an exponential decrease, characteristic of a porosity made of cracks. Correlations between measured parameters have been analyzed by comparing them with relations deduced from theoretical models. Linking the formation factor to the porosity leads to a rather low tortuosity value (about 2.4), characterizing a medium with a well connected porosity. Correlation between permeabilityk and formation factorF leads to a power-law relationk F –n wheren2.9, which is consistent with a crack model describing the behavior of the thermally treated rock.  相似文献   

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

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