首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
《Marine Geology》2005,216(3):169-189
Simultaneous high frequency field measurements of water depth, flow velocity and suspended sediment concentration were made at three fixed locations across the high tide swash and inner surf zones of a dissipative beach. The dominant period of the swash motion was 30–50 s and the results are representative of infragravity swash motion. Suspended sediment concentrations, loads and transport rates in the swash zone were almost one order of magnitude greater than in the inner surf zone. The vertical velocity gradient near the bed and the resulting bed shear stress at the start of the uprush was significantly larger than that at the end of the backwash, despite similar flow velocities. This suggests that the bed friction during the uprush was approximately twice that during the backwash.The suspended sediment profile in the swash zone can be described reasonably well by an exponential shape with a mixing length scale of 0.02–0.03 m. The suspended sediment transport flux measured in the swash zone was related to the bed shear stress through the Shields parameter. If the bed shear stress is derived from the vertical velocity gradient, the proportionality coefficient between shear stress and sediment transport rate is similar for the uprush and the backwash. If the bed shear stress is estimated using the free-stream flow velocity and a constant friction factor, the proportionality factor for the uprush is approximately twice that of the backwash. It is suggested that the uprush is a more efficient transporter of sediment than the backwash, because the larger friction factor during the uprush causes larger bed shear stresses for a given free-stream velocity. This increased transport competency of the uprush is necessary for maintaining the beach, otherwise the comparable strength and greater duration of the backwash would progressively remove sediment from the beach.  相似文献   

2.
Direct bed shear stress measurements in bore-driven swash   总被引:1,自引:0,他引:1  
Direct measurements of bed shear in the swash zone are presented. The data were obtained using a shear plate in medium and large-scale laboratory bore-driven swash and cover a wide range of bed roughness. Data were obtained across the full width of the swash zone and are contrasted with data from the inner surf zone. Estimates of the flow velocities through the full swash cycle were obtained through numerical modelling and calibrated against measured velocity data. The measured stresses and calculated flow velocities were subsequently used to back-calculate instantaneous local skin friction coefficients using the quadratic drag law. The data show rapid temporal variation of the bed shear stress through the leading edge of the uprush, which is typically two–four times greater than the backwash shear stresses at corresponding flow velocity. The measurements indicate strong temporal variation in the skin friction coefficient, particularly in the backwash. The general behaviour of the skin friction coefficient with Reynolds number is consistent with classical theory for certain stages of the swash cycle. A spatial variation in skin friction coefficient is also identified, which is greatest across the surf-swash boundary and likely related to variations in local turbulent intensities. Skin friction coefficients during the uprush are approximately twice those in the backwash at corresponding Reynolds number and cross-shore location. It is suggested that this is a result of the no-slip condition at the tip leading to a continually developing leading edge and boundary layer, into which high velocity fluid and momentum are constantly injected from the flow behind and above the tip region. Finally, the measured stress data are used to determine the asymmetry and cross-shore variation in potential sediment transport predicted by three forms of sediment transport formulae.  相似文献   

3.
《Coastal Engineering》2006,53(4):335-347
This paper investigates cross-shore profile changes of gravel beaches, with particular regard to discussing the tendency for onshore transport and profile steepening in the swash zone. The discussion includes observed morphological changes on a gravel beach from experimental investigations at the Large Wave Flume (GWK) in Hanover, Germany. During the tests all the profile changes occurred in the swash zone, resulting in erosion below the still water line (SWL) and formation of a berm above the SWL. We investigate the profile evolution evaluating the transport rates from a bed load sediment transport formulation coupled with velocities calculated from a set of Boussinesq equations that have been validated for its use in the surf and swash zones [Lynett, P.J., Wu, T.-R., and Liu, L.-F., P., 2002. Modelling wave runup with depth-integrated equations. Coastal Engineering, 46, 89–107; Otta, A.K., and Pedrozo-Acuña, A., 2004. Swash boundary and cross-shore variation of horizontal velocity on a slope. In: J.M. Smith (Editor), Proceedings 29th International Conference on Coastal Engineering. World Scientific, Lisbon, Portugal, pp. 1616–1628]. We discuss the influence of bottom friction on the predicted profiles, using reported friction factors from experimental studies. It is shown that the use of a different friction factor within a realistic range in each phase of the swash (uprush and backwash) improves prediction of the beach profiles, although quantitative agreement between the measured and computed profile evolutions is not satisfactory. Furthermore, if the friction factor and the transport efficiency (C) of the sediment transport formulation are kept the same in the uprush and backwash, accurate representation of profile evolution is not possible. Indeed, the features of the predicted profiles are reversed. However, when the C parameter is set larger during the uprush than during the backwash, the predicted profiles are closer to the observations. Differences between the predicted profiles from setting non-identical C-values and friction factors for the swash phase, are believed to be linked to both the infiltration effects on the flow above the beachface and the more accelerated flow in the uprush.  相似文献   

4.
A new model for the boundary layer development and associated skin friction coefficients and shear stress within the swash zone is presented. The model is developed within a Lagrangian reference frame, following fluid trajectories, and can be applied to both laminar flow and smooth turbulent flow. The model is based on the momentum integral approach for steady, flat-plate boundary layers, with appropriate modifications to account for the unsteady flow regime and flow history. The model results are consistent with previous measurements of bed shear stress and skin friction coefficients within the swash zone. These indicate strong temporal and spatial variation throughout the swash cycle, and a clear distinction between the uprush and backwash phase. This variation has been previously attributed the unsteady flow regime and flow history effects, both of which are accounted for in the new model. Fluid particle trajectories and velocity are computed using the non-linear shallow water wave equations and the boundary layer growth across the entire swash zone is estimated. Predictions of the bed shear stress and skin friction coefficients agree reasonably well with direct bed shear stress measurements reported by Barnes et al. (Barnes, M.P., O’Donaghue, T., Alsina, J.M., Baldock, T.E., 2009. Direct bed shear stress measurements in bore-driven swash. Coastal Engineering 56 (8), 853–867) and, for a given flow velocity, give stresses which are consistent with the bias toward uprush sediment transport which has consistently been observed in measurements. The data and modelling suggest that the backwash boundary layer is initially laminar, which results in the late development of significant bed shear during the backwash, with a transition to a turbulent boundary layer later in the backwash. A new conceptual model for the boundary layer structure at the leading edge of the swash is proposed, which accounts for both the no-slip condition at the bed and the moving wet–dry interface. However, further development of the Lagrangian Boundary Layer Model is required in order to include bore-generated turbulence and to account for variable roughness and mobile beds.  相似文献   

5.
《Marine Geology》2004,203(1-2):109-118
Spatial variations in sediment load in the swash uprush and textural properties of sediment in transport were evaluated to investigate the mechanisms responsible for sediment transport during wave uprush. Four streamer traps were deployed at 2.0-m intervals across the swash zone of a sheltered, microtidal sandy beach at Port Beach, Western Australia, over a 4-day period. During these trapping experiments, offshore significant wave heights were 0.3–0.5 m and wave periods were about 10 s. The average width of the uprush zone was 6.9 m and the average uprush duration was 5.9 s. Cross-shore distributions of sediment load for 70 uprush events reveal a maximum in sediment load landward of the base of the swash (at about 20% of swash width) during single events and a maximum closer to mid-swash (at about 40% of swash width) during multiple events characterized by swash interactions. Settling velocity distributions of trap samples during individual uprush events are similar to distributions found on the beach surface, with the lowest settling velocities (finest sediments) near the base of the swash zone and maximum settling velocities (coarsest sediments) around the mid-swash position. It was found that sediment transport during wave uprush occurs through two distinct mechanisms: (1) sediment entrainment during bore collapse seaward of the base of the swash zone and subsequent advection of this bore-entrained sediment up the beach by wave uprush; and (2) in situ sediment entrainment and transport induced by local shear stresses during wave uprush. Both mechanisms are considered important, but the first mechanism is considered most significant during the early stages of wave uprush when sediment is transported mainly in suspension, while the second mechanism is likely to dominate the mid- to later stages of wave uprush when sediment is transported mainly by sheet flow. The relative importance of the two mechanisms will vary between different beaches with the morphodynamic state of the beach (reflective versus dissipative) expected to play a major role.  相似文献   

6.
Hydrodynamics and sediment transport in the nearshore zone were modeled numerically taking into account turbulent unsteady flow. The flow field was computed using the Reynolds Averaged Navier–Stokes equations with a kε turbulence closure model, while the free surface was tracked using the Volume-Of-Fluid technique. This hydrodynamical model was supplemented with a cross-shore sediment transport formula to calculate profile changes and sediment transport in the surf and swash zones. Based on the numerical solutions, flow characteristics and the effects of breaking waves on sediment transport were studied. The main characteristic of breaking waves, i.e. the instantaneous sediment transport rate, was investigated numerically, as was the spatial distribution of time-averaged sediment transport rates for different grain sizes. The analysis included an evaluation of different values of the wave friction factor and an empirical constant characterizing the uprush and backwash. It was found that the uprush induces a larger instantaneous transport rate than the backwash, indicating that the uprush is more important for sediment transport than the backwash. The results of the present model are in reasonable agreement with other numerical and physical models of nearshore hydrodynamics. The model was found to predict well cross-shore sediment transport and thus it provides a tool for predicting beach morphology change.  相似文献   

7.
《Coastal Engineering》2005,52(1):1-23
We develop solutions for the transport of suspended sediment by a single swash event following the collapse of a bore on a plane beach, and we investigate the morphodynamical role that such transport may play. Although the intrinsic asymmetry between uprush and backwash velocities tends to encourage the export of sediment, we find that swash events may be effective in distributing across the swash zone much or all of the sediment mobilised by bore collapse; additionally, settling lag effects may promote a weak onshore movement of sediment. We quantify both effects in terms of the properties of the sediment and of the swash event, and comment on the relationship between our findings and recent field studies of swash zone sediment transport.  相似文献   

8.
《Coastal Engineering》2005,52(7):633-645
New experimental laboratory data are presented on swash overtopping and sediment overwash on a truncated beach, approximating the conditions at the crest of a beach berm or inter-tidal ridge-runnel. The experiments provide a measure of the uprush sediment transport rate in the swash zone that is unaffected by the difficulties inherent in deploying instrumentation or sediment trapping techniques at laboratory scale. Overtopping flow volumes are compared with an analytical solution for swash flows as well as a simple numerical model, both of which are restricted to individual swash events. The analytical solution underestimates the overtopping volume by an order of magnitude while the model provides good overall agreement with the data and the reason for this difference is discussed. Modelled flow velocities are input to simple sediment transport formulae appropriate to the swash zone in order to predict the overwash sediment transport rates. Calculations performed with traditional expressions for the wave friction factor tend to underestimate the measured transport. Additional sediment transport calculations using standard total load equations are used to derive an optimum constant wave friction factor of fw = 0.024. This is in good agreement with a broad range of published field and laboratory data. However, the influence of long waves and irregular wave run-up on the overtopping and overwash remains to be assessed. The good agreement between modelled and measured sediment transport rates suggests that the model provides accurate predictions of the uprush sediment transport rates in the swash zone, which has application in predicting the growth and height of beach berms.  相似文献   

9.
《Coastal Engineering》2001,42(1):35-52
Measurements were obtained from the swash-zone of a high-energy macrotidal dissipative beach of pore-pressure at four levels below the bed, and cross-shore velocity at a single height above the bed. Time-series from relatively high (Hs≈2.0 m) energy conditions were chosen for analysis from the mid-swash-zone at high tide. Calculation of upwards-directed pore-pressure gradients shows that, in this case, fluidisation of the upper layer of sediment, leading to enhanced backwash transport, is unlikely. An apparent conflict exists in the literature regarding the net effect of infiltration–exfiltration on the sediment transport, through the combined effects of stabilisation–destabilisation and boundary layer modification. This is examined for the data collected using a modified Shields parameter. Inferred instantaneous transport rates integrated over a single swash cycle show a decrease in uprush transport of 10.5% and an increase in backwash transport of 4.5%. Sensitivity tests suggest that there is a critical grain size at which the influence of infiltration–exfiltration changes from offshore to onshore. The exact value of this grain size is highly sensitive to the method used to estimate the friction factor.  相似文献   

10.
Novel laboratory experiments and numerical modelling have been performed to study the advection scales of suspended sediment in the swash zone. An experiment was designed specifically to measure only the sediment picked up seaward of the swash zone and during bore collapse. The advection scales and settling of this sediment were measured during the uprush along a rigid sediment-free beach face by a sediment trap located at varying cross-shore positions. Measurements were made using a number of repeated solitary broken waves or bores. Approximately 25% of the pre-suspended sediment picked up by the bores reaches the mid-swash zone (50% of the horizontal run-up distance), indicating the importance of the sediment advection in the lower swash zone. The pre-suspended sediment is sourced from a region seaward of the shoreline (still water line) which has a width of about 20% of the run-up distance. An Eulerian–Lagrangian numerical model is used to model the advection scales of the suspended sediment. The model resolves the hydrodynamics by solving the non-linear shallow water equations in an Eulerian framework and then solves the advection–diffusion equation for turbulence and suspended sediment in a Lagrangian framework. The model provides good estimates of the measured mass and distribution of sediment advected up the beach face. The results suggest that the correct modelling of turbulence generation prior to and during bore collapse and the advection of the turbulent kinetic energy into the lower swash is important in resolving the contribution of pre-suspended sediment to the net sediment transport in the swash zone.  相似文献   

11.
The result of field experiments, designed to investigate the relative proportions of bedload and suspended load, are described. The ratio of bedload to suspended sediment load in the swash zone is examined in both swash and backwash on four beaches by measuring the amounts collected in a sediment trap. Bedload transport is found to dominate the backwash. The relative proportions of bedload and suspended load change over the tidal cycle, with increasing bedload dominance at low tide. The total amount of sediment transported as swash and backwash is noticeably greater at high tide than at low tide. More sediment is transported on the flood tide than on the ebb.  相似文献   

12.
The morphological changes of multiple intertidal bars (ridges) on a macrotidal beach were examined under low-energy wave conditions during a spring-to-spring tidal cycle. The morphological response was coupled to the tidal water level variations and related residence times for swash processes and surf (breaking waves and bores) over the cross-shore profile. Spring tides induced a large spatial variation in water lines and small residence times for distinct processes. Neap tides narrowed the intertidal area and increased the time for certain processes to work on the sediment at one location. The observed morphological changes could be coupled to the stagnation of processes at a certain bar crest position. The action of surf (breaking waves and bores) played the major role in the onshore migration of the intertidal bars and the simultaneous erosion of the seaward flank. Swash action, responsible for the generation and migration of intertidal bars in microtidal settings, was not the dominant process in causing the observed morphological changes. Intertidal ridges on macrotidal beaches cannot be considered swash bars as suggested by most previous investigations into these morphological features.  相似文献   

13.
This paper describes newly obtained, high-frequency observations of beach face morphological change over numerous tidal cycles on a macrotidal sandy beach made using a large array of ultrasonic altimeters. These measurements enable the net cross-shore sediment fluxes associated with many thousands of individual swash events to be quantified. It is revealed that regardless of the direction of net morphological change on a tidal time scale, measured net fluxes per event are essentially normally distributed, with nearly equal numbers of onshore and offshore-directed events. The majority of swash events cause net cross-shore sediment fluxes smaller than ± 50 kg m− 1 and the mean sediment flux per swash event is only O(± 1 kg m− 1) leading to limited overall morphological change. However, much larger events which deposit or remove hundreds of kilograms of sand per meter width of beach occur at irregular intervals throughout the course of a tide. It was found that swash–swash interactions tend to increase the transport potential of a swash event and the majority of the swash events that cause these larger values of sediment flux include one or more interactions. The majority of the larger sediment fluxes were therefore measured in the lower swash zone, close to the surf/swash boundary where swash–swash interactions are most common. Despite the existence of individual swash events that can cause fluxes of sediment that are comparable to those observed on a tidal time scale, frequent reversals in transport direction act to limit net transport such that the beach face volume remains in a state of dynamic equilibrium and does not rapidly erode or accrete.  相似文献   

14.
T.D. Price  B.G. Ruessink   《Marine Geology》2008,251(1-2):98-109
This paper builds on the work of Masselink [Masselink, G., 1993. Simulating the effects of tides on beach morphodynamics. J. Coast. Res. SI 15, 180–197.] on the use of the residence times of shoaling waves, breaking waves and swash/backwash motions across a cross-shore profile to qualitatively understand temporal beach behaviour. We use a data set of in-situ measurements of wave parameters (height and period) and water depth, and time-exposure video images overlooking our single-barred intertidal measurement array at Egmond aan Zee (Netherlands) to derive boundaries between the shoaling zone, the surf zone and the swash zone. We find that the boundaries are functional dependencies of the local relative wave height on the local wave steepness. This contrasts with the use of constant relative wave heights or water levels in earlier work. We use the obtained boundaries and a standard cross-shore wave transformation model coupled to an inner surf zone bore model to show that large (> 5) relative tide ranges (RTR, defined as the ratio tide range–wave height) indicate shoaling wave processes across almost the entire intertidal profile, with surf processes dominating on the beach face. When the RTR is between 2 and 5, surf processes dominate over the intertidal bar and the lower part of the beach face, while swash has the largest residence times on the upper beach face. Such conditions, associated with surf zone bores propagating across the bar around low tide, were observed to cause the intertidal bar to migrate onshore slowly and the upper beach face to steepen. For RTR values less than about 2, surf zone processes dominate across the intertidal bar, while the dominance of swash processes now extends across most of the beach face. The surf zone processes were now observed to lead to offshore bar migration, while the swash eroded the upper beach face.  相似文献   

15.
The aim of the investigation was to define the mechanisms of sediment transport in the swash zone of microtidal coarse-clastic beaches in the very short term by evaluating the displacement rates of marked pebbles under low-energy wave conditions. Tests were performed at two sites (Marina di Pisa, Ligurian Sea, and Portonovo, central Adriatic Sea) to check the consistency of the data over a range of different grain sizes. Two recovery campaigns were carried out at both sites, one 6 h and the other 24 h after the injection. During the experiments wave action was at a minimum (wave heights never exceeded 0.3 m). The results show that 20% of pebbles ranging in diameter from 30–90 mm moved significantly (more than 0.5 m) already 6 h after the injection, with some tracers being lost (3%). After 24 h, 40% of the pebbles were significantly displaced and 10% were lost. The preferential downslope movement of tracers, which suggests that coarse sediment movement under low-energy conditions is mainly controlled by gravity processes enhanced by steep beachface slopes, represents the novelty of the results reported here. It would appear that swash processes on low-energy beaches cause a significant rate of pebble displacement through the destabilization induced by wave uprush and backwash. Despite the microtidal range, the position of the mean water level plays a major role in changing the beach level at which swash processes can actually trigger pebble movement. The results of this study show that considerable, and mostly seaward-directed, coarse sediment transport takes place even during short fair-weather periods.  相似文献   

16.
本文提出海滩反递变纹层自下而上粒度由小到大,重矿物富集于纹层底部,它是前滨冲流“剪切分选”的产物。原生反递变纹层常被激浪破坏和再造,遇到后期加积海滩过程时,才能保存于海滩层理中。海滩层理的现场观测是研究海滩层理反递变纹层形成机理的重要方法之一。  相似文献   

17.
C.F. Jago  J. Hardisty 《Marine Geology》1984,60(1-4):123-154
The foreshore of Pendine Sands forms the seaward part of an extensive, sandy coastal barrier in a shallow Carmarthen Bay, SW Wales. The sedimentological features of the macrotidal foreshore reflect a tide-induced modification of nearshore wave characteristics. As the tide ebbs, the breaker height may decrease, the surf zone widens and becomes increasingly dissipative, and swash/backwash velocities diminish. A concomitant change from plunging to spilling breakers and increasingly symmetrical swash zone flows are associated with a decreasing beach gradient.

A zero net transport model demonstrates that the beach profile is self-stabilising in the short-term, and periodic levelling has shown that the beach is in long-term equilibrium with prevailing conditions, though this does not preclude a significant dynamic response to changing tides and waves.

The flow regimes of wave-generated currents decline as the tide ebbs, and normal beach processes do not usually affect the lower foreshore. Accordingly, there is an overall seaward-fining of the primary framework component of the sands. In more detail, this framework component displays a slight seaward-coarsening across an upper foreshore dominated by high water swash and surf; a rapid seaward-fining across the mid-foreshore in response to the ebb-attenuating swash zone flow velocities; and a slight seaward-fining across the lower foreshore under the action of nearshore shoaling waves. Bedforms vary from a swash/backwash emplaced flat bed across the upper foreshore to the small ripples of nearshore asymmetric oscillatory flows across the lower foreshore.

The surface sediment veneer is not representative of the subsurface sediments which form in response partly to fairweather conditions, partly to storms. The upper foreshore is characterised by swash/backwash emplaced plane bedding in fine sands frequently disrupted by bubble cavities. The mid-foreshore is composed of coarser-grained shelly traction clogs arranged as landward- and seaward-dipping large-scale cross bedding and/or plane bedding; these are probably storm breaker/surf deposits. The lower foreshore, though partially and sometimes totally bioturbated, shows landward-dipping small-scale cross bedding in very fine sands sorted by nearshore shoaling waves.

Tide- and storm-induced modification of the nearshore flow regimes therefore produces a distinctive shore-normal array of sedimentary facies. Each facies is characterised by diagnostic textural and structural signatures. A prograding sequence of such macrotidal deposits would be similar to, but more extensive than, a comparable microtidal sequence.  相似文献   


18.
Field experiments were conducted on a low-gradient, high-energy sandy beach (Truc Vert, France) and a steep, low-energy gravel beach (Slapton, UK) to examine alongshore-directed currents within the swash zone. At Truc Vert, data were collected over 33 tidal cycles with offshore significant wave heights of 1–4 m and periods of 5–12 s. At Slapton data were collected during 12 tides with wave heights of 0.3–1 m and periods of 4–9 s. The swash motion was predominantly at infragravity frequencies at Truc Vert and incident frequencies at Slapton.  相似文献   

19.
New laboratory experiments have produced detailed measurements of hydrodynamics within swash generated by bore collapse on a steep beach. The experiments are based on a dambreak rig producing a highly repeatable, large-scale swash event, enabling detailed measurements of depths and velocities at a number of locations across the swash zone. Experiments were conducted on two beaches, differentiated by roughness. Results are presented for uprush shoreline motion, flow depths, depth-averaged velocity, velocity profiles and turbulence intensity. Estimates of the time- and spatially-varying bed shear stress are obtained via log-law fitting to the velocity profiles and are compared with the shear plate measurements of Barnes et al. (2009) for similar experimental conditions. Experimental results are compared with model predictions based on a NLSWE model with momentum loss parameterised using the simple quadratic stress law in terms of the depth-averaged velocity. Predicted and measured flow depths and depth-averaged velocities agree reasonably well for much of the swash period, but agreement is not good at the time of bore arrival and towards the end of the backwash. The parameterisation of total momentum loss via the quadratic stress law cannot adequately model the swash bed shear stress at these critical times.  相似文献   

20.
Abstract

On shingle beaches, changes in foreshore elevation and sediment distribution landward of the break point are produced largely by variations in the uprush and backwash of waves. However, very little is known about the forces active in this zone.

A field instrument system which senses and records some of the parameters thought to influence beach erosion and deposition in this zone has been constructed. The equipment is also suitable for the investigation of a number of other shore and nearshore processes including erosion on sandy and rocky shores, and flow processes affecting littoral biological communities.

In the swash zone two sensing heads, a dynamometer and a depth recorder, sense variations in uprush and backwash velocities, energies, discharges, and depths of flow. Both devices are electromechanical and are coupled to a recording unit on land by PVC‐insulated cable. The dynamometer (two force plates mounted back‐to‐back on a compression spring and coupled to variable resistances) has been calibrated, statically and in a flume, to obtain velocity determinations accurate to within 10 cm . sec?1 of true flow speed. Average swash zone velocities lie between 100 and 300 cm . sec?1.

A parallel‐wire resistance gauge mounted an a stilling tube records flow depths. As water level rises and falls in the tube it alters resistance in a control circuit. The land unit, amplifiers and a strip‐chart recorder, receives the output from the dynamometer and flow depth gauge. The recorder is equipped with a trip‐pen so that analysis of wave periods or other variables is possible in the field. With poles at known spacings across the shore and the trip‐pen records, velocity distributions across the swash zone can be obtained. Measurements of velocity made near the bed with the dynamometer can then be related to the local surface velocity profile.

Problems with the instrument system include inability to record velocities at several points simultaneously, and unreliable records of backwash parameters with low breakers on shingle beaches because of the small volume of flow and rapid percolation of water into the beach face.  相似文献   

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

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