首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 187 毫秒
1.
A discussion is presented about the mechanisms that govern the spatial and seasonal variability in sand-wave height and migration speed in the 4 km wide Marsdiep tidal inlet, the Netherlands. Since 1998, current velocities and water depths have been recorded with an ADCP that is mounted under the ferry ‘Schulpengat’. In this paper, the current measurements were used to explain the sand-wave observations presented in Buijsman and Ridderinkhof [this issue. Long-term evolution of sand waves in the Marsdiep inlet. I: high-resolution observations. Continental Shelf Research, doi:10.1016/j.csr.2007.10.011]. Across nearly the entire inlet, the sand waves migrate in the flood direction. In the flood-dominated southern part of the inlet, the ‘measured’ (i.e. based on sand-wave shape and migration speed) and predicted bedload transport agree in direction, magnitude, and trends, whereas in the ebb-dominated northern part the predicted bedload and suspended load transport is opposite to the sand-wave migration. In the southern part, 55% of the bedload transport is due to tidal asymmetries and 45% due to residual currents. In addition to the well-known tidal asymmetries, asymmetries that arise from the interaction of M2M2 and its overtides with S2S2 and its compound tides are also important. It is hypothesised that in the northern part of the inlet the advection of suspended sand and lag effects govern the sand-wave migration. The relative importance of suspended load transport also explains why the sand waves have smaller lee-slope angles, are smaller, more rounded, and more three-dimensional in the northern half of the inlet. The sand waves in this part of the inlet feature the largest seasonal variability in height and migration speed. This seasonal variability may be attributed to the tides or a seasonal fluctuation in fall velocity. In both cases sediment transport is enhanced in winter, increasing sand-wave migration and decreasing sand-wave height. The influence of storms and estuarine circulation on the sand-wave variability is negligible.  相似文献   

2.
We deployed bottom-mounted quadrapod equipped with acoustic Doppler current profiler (ADCP), acoustic Doppler velocimeter (ADV), and optical backscatter sensor (OBS) over two semidiurnal tidal cycles along the western coast of the Yellow Sea, China. In combination with shipboard profiling of CTD and LISST-100, we resolved the temporal and spatial distributions of tidal currents, turbulent kinetic energy (TKE), suspended sediment concentration (SSC) and particle size distributions. During the observations, tidal-induced bottom shear stress was the main stirring factor. However, weak tidal flow during the ebb phase was accompanied by two large SSC and median size events. The interactions of seiche-induced oscillations with weak ebb flow induced multiple flow reversals and provided a source of turbulence production, which stripped up the benthic fluff layers (only several millimeters) around the Jiaozhou Bay mouth. Several different methods for inferring mean suspended sediment settling velocity agreed well under peak currents, including estimates using LISST-based Stokes’ settling law, and ADCP-based Rouse profiles, ADV-based inertial-dissipation balance and Reynolds flux. Suspended particles in the study site can be roughly classified into two types according to settling behavior: a smaller, denser class consistent with silt and clay and a larger, less dense class consistent with loosely aggregated flocs. In the present work, we prove that acoustic approaches are robust in simultaneously and non-intrusively estimating hydrodynamics, SSC and settling velocities, which is especially applicable for studying sediment dynamics in tidal environments with moderate concentration levels.  相似文献   

3.
The echo-amplitude of a 23-m-deep bottom-mounted acoustic Doppler current profiler (ADCP) shows regular spikes up to 30 dB above background level when a ship passes nearby, due to deep penetration of bubble clouds. This is evidenced from regularly occurring spikes in echo-data that are simultaneous with ferry crossings in a narrow sea-strait. The bubbles can nearly reach the bottom and are comparable in magnitude to near-bottom scattering off suspended material in vigorous tidal currents exceeding 1 m s−1 in magnitude. The bubble clouds mask the sea surface from the echo-amplitude, which hampers the use of an ADCP for estimating atmospheric parameters and near-surface currents, under such conditions. The echo-spikes associated with the ferry are confirmed with coinciding dips in bottom pressure up to 1200 N m−2 and with deviations up to 10° in the ADCP's heading due to pressure waves and magnetic field disturbances from under the ferry and from its rear, respectively.  相似文献   

4.
In tidal environments, the response of suspended sediment concentration (SSC) to the current velocity is not instantaneous, the SSC lagging behind the velocity (phase lag), and the amplitude of SSC variation decreasing with height above the bed (amplitude attenuation). In order to quantitatively describe this phenomenon, a one-dimensional vertical advection–diffusion equation of SSC is derived analytically for uniform unsteady tidal flow by defining a concentration boundary condition using a constant vertical eddy diffusivity and sediment settling velocity. The solution, in simple and straightforward terms, shows that the vertical phase lag increases linearly with the height above the bed, while the amplitude of the SSC variation decreases exponentially with the height. The relationship between the SSC and the normalized current velocity can be represented by an ellipse or a line, depending on the phase lag. The lag of sediment movement or “diffusion/settling lag” is the mechanism generating the phase lag effect. Field observations used for validation show that the theoretically predicted and the observed curves of the vertical SSC phase lag and amplitude attenuation show reasonable agreement. The procedure proposed in this paper substantially simplifies the modeling of suspended matter transport in tidal flows.  相似文献   

5.
A tidal bore is a water discontinuity at the leading edge of a ood tide wave in estuaries with a large tidal range and funneling topography. New measurements were done in the Garonne River tidal bore on 14 15 November 2016, at a site previously investigated between 2010 and 2015. The data focused on long, continuous, high-frequency records of instantaneous velocity and suspended sediment con- centration (SSC) estimate for several hours during the late ebb, tidal bore passage and ood tide. The bore passage drastically modi ed the ow eld, with very intense turbulent and sediment mixing. This was evidenced with large and rapid uctuations of both velocity and Reynolds stress, as well as large SSCs during the ood tide. Granulometry data indicated larger grain sizes of suspended sediment in water samples compared to sediment bed material, with a broader distribution, shortly after the tidal bore. The tidal bore induced a sudden suspended sediment ux reversal and a large increase in suspended sedi- ment ux magnitude. The time-variations of turbulent velocity and suspended sediment properties indicated large uctuations throughout the entire data set. The ratio of integral time scales of SSC to velocity in the x-direction was on average TE,SSC/TE,x 0.16 during the late ebb tide, compared to TE,SSC/ TE,x 0.09 during the late ood tide. The results imply different time scales between turbulent velocities and suspended sediment concentrations.  相似文献   

6.
The effect of horizontal grid resolution on the horizontal relative dispersion of particle pairs has been investigated on a short time scale, i.e. one tidal M 2 cycle. Of particular interest is the tidal effect on dispersion and transports in coastal waters where small-scale flow features are important. A three-dimensional ocean model has been applied to simulate the tidal flow through the Moskstraumen Maelstrom outside Lofoten in northern Norway, well known for its strong current and whirlpools (Gjevik et al., Nature 388(6645):837–838, 1997; Moe et al., Cont Shelf Res 22(3):485–504, 2002). Simulations with spatial resolution down to 50 m have been carried out. Lagrangian tracers were passively advected with the flow, and Lyapunov exponents and power law exponents have been calculated to analyse the separation statistics. It is found that the relative dispersion of particles on a short time scale (12–24 h) is very sensitive to the grid size and that the spatial variability is also very large, ranging from 0 to 100 km2 over a distance of 100 m. This means that models for prediction of transport and dispersion of oil spills, fish eggs, sea lice etc. using a single diffusion coefficient will be of limited value, unless the models actually resolves the small-scale eddies of the tidal current.  相似文献   

7.
In a tidal channel with adjacent tidal flats, along–channel momentum is dissipated on the flats during rising tides. This leads to a sink of along–channel momentum. Using a perturbative method, it is shown that the momentum sink slightly reduces the M2 amplitude of both the sea surface elevation and current velocity and favours flood dominant tides. These changes in tidal characteristics (phase and amplitude of sea surface elevations and currents) are noticeable if widths of tidal flats are at least of the same order as the channel width, and amplitudes and gradients of along–channel velocity are large. The M2 amplitudes are reduced because stagnant water flows from the flats into the channel, thereby slowing down the current. The M4 amplitudes and phases change because the momentum sink acts as an advective term during the fall of the tide, such a term generates flood dominant currents. For a prototype embayment that resembles the Marsdiep–Vlie double–inlet system of the Western Wadden Sea, it is found that for both the sea surface elevation and current velocity, including the momentum sink, lead to a decrease of approximately 2% in M2 amplitudes and an increase of approximately 25% in M4 amplitudes. As a result, the net import of coarse sediment is increased by approximately 35%, while the transport of fine sediment is hardly influenced by the momentum sink. For the Marsdiep–Vlie system, the M2 sea surface amplitude obtained from the idealised model is similar to that computed with a realistic three–dimensional numerical model whilst the comparison with regard to M4 improves if momentum sink is accounted for.  相似文献   

8.
In a series of observing system simulations, we test whether the Atlantic meridional overturning circulation (AMOC) can be observed based on the existing Line W deep western boundary array. We simulate a Line W array, which is extended to the surface and to the east to cover the basin to the Bermuda Rise. In the analyzed ocean circulation model ORCA025, such an extended Line W array captures the main characteristics of the western boundary current. Potential trans-basin observing systems for the AMOC are tested by combining the extended Line W array with a mid-ocean transport estimate obtained from thermal wind “measurements” and Ekman transport to the total AMOC (similarly to Hirschi et al., Geophys Res Lett 30(7):1413, 2003). First, we close Line W zonally supplementing the western boundary array with several “moorings” in the basin (Line W-32°N). Second, we supplement the western boundary array with a combination of observations at Bermuda and the eastern part of the RAPID array at 26°N (Line W-B-RAPID). Both, a small number of density profiles across the basin and also only sampling the eastern and western boundary, capture the variability of the AMOC at Line W-32°N and Line W-B-RAPID. In the analyzed model, the AMOC variability at both Line W-32°N and Line W-B-RAPID is dominated by the western boundary current variability. Away from the western boundary, the mid-ocean transport (east of Bermuda) shows no significant relation between the two Line W-based sections and 26°N. Hence, a Line W-based AMOC estimate could yield an estimate of the meridional transport that is independent of the 26°N RAPID estimate. The model-based observing system simulations presented here provide support for the use of Line W as a cornerstone for a trans-basin AMOC observing system.  相似文献   

9.
We present semi-analytical solutions for suspended sediment concentration (SSC) and residual sediment transport in a simple mathematical model of a short tidal embayment. These solutions allow us to investigate in some detail the characteristic tidal and semi-tidal variation of SSC and the processes leading to residual sediment transport, including settling and scour lags, the roles of ‘local’ and ‘advective’ contributions, and the presence of internally or externally generated overtides. By interpreting the transport mechanisms in terms of the classic conceptual models of settling lag we clarify how these models may be expressed in mathematical terms. Our results suggest that settling lag is usually a more important process than scour lag, and that a local model which neglects advection may predict the direction of net sediment transport incorrectly. Finally, we discuss our results in the context of other transport processes and morphodynamic feedback.  相似文献   

10.
The potential for rapid coastline modification in the face of sea-level rise or other stressors is alarming, since coasts are often densely populated and support valuable infrastructure. In addition to coastal submergence, nutrient-related water pollution is a growing concern for coastal wetlands. Previous studies found that the Suspended Sediment Concentration (SSC) of coastal wetlands acts as a first-order control of their sustainability, but SSC dynamics are poorly understood. Our study focuses on the Virginia Coast Reserve (VCR) Long Term Ecological Research (LTER) site, a shallow multiple tidal inlet system in the USA. We apply numerical modelling (Delft3D-SWAN) and subsequent analyses to determine SSC dynamics within the VCR. In particular, we consider two important controls on SSC in the system: vegetation (seagrass and salt marsh) and offshore waves. Our results show that vegetation colonies and increased wave energy lengthen water residence time. The reduction in the tidal prism decreases SSC export from the bay via tidal inlets, leading to increased sediment retention in the bay. We found that alongshore currents can enhance lagoon SSC by importing fine sediments from an adjacent inlet along the coastline. Our numerical experiments on vegetation seasonality can improve the understanding of wave climate impact on coastal bay sediment budget. Offshore waves increase sediment export from coastal bays, particularly during winter seasons with low vegetation density. Therefore, our study can help managers and stakeholders to understand how to implement restoration strategies for the VCR. © 2020 John Wiley & Sons, Ltd.  相似文献   

11.
A 2D depth-averaged numerical model is set up to simulate the macro-scale hydrodynamic characteristics, sediment transport patterns and morphological evolution in Hangzhou Bay, a large macro-tidal estuary on the eastern coast of China. By incorporating the shallow water equations, the suspended sediment transport equation and the mass-balance equation for sediment; short-term hydrodynamics, sediment transport and long-term morphological evolution for Hangzhou Bay are simulated and the underlying physical mechanisms are analyzed. The model reproduces the spatial distribution patterns of suspended sediment concentration (SSC) in Hangzhou Bay, characterized by three high SSC zones and two low SSC zones. It also correctly simulates the residual flow, the residual sediment transport and the sediment accumulation patterns in Hangzhou Bay. The model results are in agreement with previous studies based on field measurements. The residual flow and the residual sediment transport are landwards directed in the northern part of the bay and seawards directed in the southern part. Sediment accumulation takes place in most areas of the bay. Harmonic analysis revealed that the tide is flood-dominant in the northern part of the bay and ebb-dominant in the southern part of the bay. The strength of the flood-dominance increases landwards along the northern Hangzhou Bay. In turn sediment transport in Hangzhou Bay is controlled by this tidal asymmetry pattern. In addition, the direction of tidal propagation in the East China Sea, the presence of the archipelago in the southeast and the funnel-shaped geometry of the bay, play important roles for the patterns of sediment transport and sediment accumulation respectively.  相似文献   

12.
To describe the exchange of water and sediment through the Venice Lagoon inlets a 3-D hydrodynamic and sediment transport model has been developed and applied to a domain comprising Venice Lagoon and a part of the Adriatic Sea. The model has been validated for both current velocities and suspended particle concentration against direct observations and from observations empirically derived fluxes from upward-looking acoustic Doppler current profiler probes installed inside each inlet. The model provides estimates of the suspended sediment transport in the lower 3 m of the water column that is not detected by acoustic Doppler current profiler sensors. The bedload model prediction has been validated against measured sand transport rates collected by sand traps deployed in the Lido and Chioggia inlets. Results indicate that, in the Lido inlet, 87% of the total load is in suspension, while the rest moves as bedload.  相似文献   

13.
Suspended sediment concentration (SSC) plays an important role in the estuarine environment.Its spatial or temporal variations in coastal zones and estuaries indicate that sediments are suspended,trans...  相似文献   

14.
The German Wadden Sea (southern North Sea) sediments are composed of both cohesive and non-cohesive deposits. The spatial distribution patterns are mainly driven by wind-induced waves and tidal currents. Transport intensity and duration depend on the hydrodynamic conditions, which vary over time. In this paper, the transport of suspended sediment was investigated on seasonal, tidal and hourly time scales in the back-barrier system of Spiekeroog Island. Long- and short-term data of fair weather periods and two storm events were investigated based on stationary and mobile measurements of currents and waves by Acoustic Doppler Current Profiler (ADCP), in situ particle size and suspended sediment concentration (SSC) measurements by laser in situ scattering and transmissometry (LISST) as well as wind records. The ADCP backscatter intensities were calibrated by means of LISST volume concentration data in order to quantify longer term SSCs and fluxes in the back-barrier system. Values up to 120 mg l−1 were recorded, but concentrations more commonly were below 60 mg l−1. The long-term results confirm former observations of a balanced budget during low-energy (fair weather) conditions in the study area. In general, SSCs were higher during spring tides than during neap tides. The data also clearly show the remobilisation of sediment by tidal current entrainment. The records include two severe storm events, “Britta” (1st November 2006) and “Kyrill” (18th January 2007). The data reveal very complex temporal flow and transport patterns. During both storm events, the export of material was mainly controlled by the interaction of wind, waves and tidal phase. The typical ebb-dominance occurring during fair-weather conditions was temporarily neutralised and even reversed to a flood-dominated situation. During “Kyrill”, the wind and high-waves setup in conjunction with the tidal phase was even able to compress the duration of two successive ebb cycles by over 70%. Although SSCs increased during both storms and higher turbulence lifted particle clouds upwards, an export of suspended matter towards the North Sea was only observed under the conditions taking place during “Britta”. Such fluxes, however, are currently still difficult to quantify because the backscatter intensity during high energy events includes a substantial amount of noise produced by the high turbulence, especially near the water surface.  相似文献   

15.
Using in situ, continuous, high frequency (8–16 Hz) measurements of velocity, suspended sediment concentration (SSC), and salinity, we investigate the factors affecting near-bed sediment flux during and after a meteorological event (cold front) on an intertidal flat in central San Francisco Bay. Hydrodynamic forcing occurs over many frequency bands including wind wave, ocean swell, seiching (500–1000 s), tidal, and infra-tidal frequencies, and varies greatly over the time scale of hours and days. Sediment fluxes occur primarily due to variations in flow and SSC at three different scales: residual (tidally averaged), tidal, and seiching. During the meteorological event, sediment fluxes are dominated by increases in tidally averaged SSC and flow. Runoff and wind-induced circulation contribute to an order of magnitude increase in tidally averaged offshore flow, while waves and seiching motions from wind forcing cause an order of magnitude increase in tidally averaged SSC. Sediment fluxes during calm periods are dominated by asymmetries in SSC over a tidal cycle. Freshwater forcing produces sharp salinity fronts which trap sediment and sweep by the sensors over short (∼30 min) time scales, and occur primarily during the flood. The resulting flood dominance in SSC is magnified or reversed by variations in wind forcing between the flood and ebb. Long-term records show that more than half of wind events (sustained speeds of greater than 5 m/s) occur for 3 h or less, suggesting that asymmetric wind forcing over a tidal cycle commonly occurs. Seiching associated with wind and its variation produces onshore sediment transport. Overall, the changing hydrodynamic and meteorological forcing influence sediment flux at both short (minutes) and long (days) time scales.  相似文献   

16.
During a one‐year period temporal and spatial variations in suspended sediment concentration (SSC) and deposition were studied on a salt and freshwater tidal marsh in the Scheldt estuary (Belgium, SW Netherlands) using automatic water sampling stations and sediment traps. Temporal variations were found to be controlled by tidal inundation. The initial SSC, measured above the marsh surface at the beginning of inundation events, increases linearly with inundation height at high tide. In accordance with this an exponential relationship is observed between inundation time and sedimentation rates, measured over 25 spring–neap cycles. In addition both SSC and sedimentation rates are higher during winter than during summer for the same inundation height or time. Although spatial differences in vegetation characteristics are large between and within the studied salt and freshwater marsh, they do not affect the spatial sedimentation pattern. Sedimentation rates however strongly decrease with increasing (1) surface elevation, (2) distance from the nearest creek or marsh edge and (3) distance from the marsh edge measured along the nearest creek. Based on these three morphometric parameters, the spatio‐temporal sedimentation pattern can be modelled very well using a single multiple regression model for both the salt and freshwater marsh. A method is presented to compute two‐dimensional sedimentation patterns, based on spatial implementation of this regression model. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

17.
Sediment rating curves are commonly used to estimate the suspended sediment load in rivers and streams under the assumption of a constant relation between discharge (Q) and suspended sediment concentrations (SSC) over time. However, temporal variation in the sediment supply of a watershed results in shifts in this relation by increasing variability and by introducing nonlinearities in the form of hysteresis or a path‐dependent relation. In this study, we used a mixed‐effects linear model to estimate an average SSC–Q relation for different periods of time within the hydrologic cycle while accounting for seasonality and hysteresis. We tested the performance of the mixed‐effects model against the standard rating curve, represented by a generalized least squares regression, by comparing observed and predicted sediment loads for a test case on the Chilliwack River, British Columbia, Canada. In our analyses, the mixed‐effects model reflected more accurate patterns of interpolated SSC from Q data than the rating curve, especially for the low‐flow summer months when the SSC–Q relation is less clear. Akaike information criterion scores were lower for the mixed‐effects model than for the standard model, and the mixed‐effects model explained nearly twice as much variance as the standard model (52% vs 27%). The improved performance was achieved by accounting for variability in the SSC–Q relation within each month and across years for the same month using fixed and random effects, respectively, a characteristic disregarded in the sediment rating curve. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

18.
A numerical shelf circulation model was developed for the Scotian Shelf, using a nested-grid setup consisting of a three-dimensional baroclinic inner model embedded inside a two-dimensional barotropic outer model. The shelf circulation model is based on the Princeton Ocean Model and driven by three-hourly atmospheric forcing provided by a numerical weather forecast model and by tidal forcing specified at the inner model's open boundaries based on pre-calculated tidal harmonic constants. The outer model simulates the depth-mean circulation forced by wind and atmospheric pressure fields over the northwest Atlantic Ocean with a horizontal resolution of 1/12°. The inner model simulates the three-dimensional circulation over the Gulf of St. Lawrence, the Scotian Shelf, and the adjacent slope with a horizontal resolution of 1/16°. The performance of the shelf circulation model is assessed by comparing model results with oceanographic observations made along the Atlantic coast of Nova Scotia and in the vicinity of Sable Island (on the Scotian Shelf) during two periods: October 2000–March 2001 and April–June 2002. Analysis of model results on Sable Island Bank indicates that tidal currents account for as much as ∼80% of the total variance of near-bottom currents, and currents driven by local winds account for ∼30% of the variance of the non-tidal near-bottom currents. Shelf waves generated remotely by winds and propagating into the region also play an important role in the near-bottom circulation on the bank.  相似文献   

19.
Few hyperpycnal flows have ever been observed in marine environments although they are believed to play a critical role in sediment dispersal within estuarine and deltaic depositional systems. The paper describes hyperpycnal flows observed in situ off the Huanghe (Yellow River) mouth, their relationship to tidal cycles, and the mechanisms that drive them. Simultaneous observations at six mooring stations during a cruise off the Huanghe mouth in the flood season of 1995 suggest that hyperpycnal flows observed at the river mouth are initiated by high concentrations of sediment input from river and modulated by tides. Hyperpycnal flows started near the end of ebb tides, when near‐bottom suspended sediment concentration (SSC) increased rapidly and salinity decreased drastically (an inverse salt wedge). The median grain size of suspended particles within the hyperpycnal layer increased, causing strong stratification of the suspended sediments in the water column. Towards the end of flood tides, the hyperpycnal flow attenuated due to frictions at the upper and lower boundaries of the flow and tidal mixing, which collapsed the stratification of the water column. Both sediment concentration and median grain size of suspended particles within the bottom layer significantly decreased. The coarser sediment particles were deposited and the hyperpycnal flows stopped. The intra‐tidal behaviors of hyperpycnal flows are closely associated with the variations of SSC, salinity, and stratification of the water column. As nearly 90% of riverine sediment is delivered to the sea during the flood seasons when hyperpycnal flows are active, hyperpycnal flows at the Huanghe mouth and the river's high sediment loads have caused rapid accretion of the Huanghe delta. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

20.
The macro-tidal Gulf of Kachchh, covering nearly 7000 km(2), is located about 150 km south of the Indus River mouth. In spite of semi-arid climate and lack of major rivers flowing into it, the Gulf is highly turbid with suspended sediment concentrations (SSC) during October-November 2002 ranging between 0.5 and 674 mgl(-1). Highly turbid waters are observed towards the northern portion of the mouth of the Gulf, at the head of the Gulf and adjacent to the numerous shoals present within the Gulf. Perennial high SSC in the Gulf is due to resuspension of sediments by strong tidal currents, shallow bathymetry and presence of fine-grained sediments on the sea floor. Numerical model studies show that there is a dynamic barrier in the central Gulf, which prevents the exchange of water and suspended sediments between the outer and inner Gulf. This dynamic barrier associated with strong east-west tidal currents restricts the turbid waters mainly to the northern Gulf, resulting in relatively clear waters (SSC<10 mgl(-1)) in the southern and central portions of the Gulf. Laser particle size distribution, clay mineralogy and geochemistry of the suspended matter show that the main source of sediments to the Gulf of Kachchh is the Indus River. Although the Indus discharge has been severely curtailed in the recent decades due to construction of numerous dams and barrages, the Gulf of Kachchh continues to receive resuspended sediments from the numerous meso and macro-tidal creeks of the Indus delta. The sediments at the head of the Gulf appear to be a mixture of sediments derived from the Indus as well as the numerous seasonal rivers draining the Rann of Kachchh.  相似文献   

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

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