首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Large asymmetric bedforms known as dunes commonly dominate the bed of sand rivers. Due to the turbulence generation over their stoss and lee sides, dunes are of central importance in predicting hydraulic roughness and water levels. During floods in steep alluvial rivers, dunes are observed to grow rapidly as flow strength increases, undergoing an unstable transition regime, after which they are washed out in what is called upper stage plane bed. This transition of dunes to upper stage plane bed is associated with high transport of bed sediment in suspension and large decrease in bedform roughness. In the present study, we aim to improve the prediction of dune development and dune transition to upper stage plane bed by introducing the transport of suspended sediment in an existing dune evolution model. In addition, flume experiments are carried out to investigate dune development under bed load and suspended load dominated transport regimes, and to get insight in the time scales related to the transition of dunes to upper stage plane bed. Simulations with the extended model including the transport of suspended sediment show significant improvement in the prediction of equilibrium dune parameters (e.g. dune height, dune length, dune steepness, dune migration rate, dune lee side slope) both under bed load dominant and suspended load dominant transport regimes. The chosen modeling approach also allows us to model the transition of dunes to upper stage plane bed which was not possible with the original dune evolution model. The extended model predicts change in the dune shapes as was observed in the flume experiments with decreasing dune heights and dune lee slopes. Furthermore, the time scale of dune transition to upper stage plane bed was quite well predicted by the extended model. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

2.
Salinity difference between terrestrial river discharge and oceanic tidal water plays a role in modifying the local flow field and, as a consequence, estuarine morphodynamics. Although widely recognized, recent numerical studies exploring the long-term morphological evolution of river-influenced estuaries with two-dimensional, depth-averaged models have mostly neglected salinity. Using a three-dimensional morphodynamic model, we aim to gain more insight into the effect of salinity on the morphodynamics of fluvio-deltaic systems. Model results indicate that the resultant estuarine morphology established after 600 years differs remarkably when a salinity gradient is included. A fan-shaped river-mouth delta exhibits less seaward expansion and is cut through by narrower channels when salinity is included. The inclusion of salinity tends to generate estuarine circulation, which favours landward sediment transport and hence limits the growth of the delta while enhancing the development of intertidal areas. The formation of deltaic channel–shoal patterns resulting from morphodynamic evolution tends to strengthen salinity stratification, which is characterized by an increased gradient Richardson number. The direction of the depth-averaged residual sediment transport over a tide may be opposite to the direction of residual velocity, indicating the significant influence of baroclinic effects on the net sediment transport direction (and hence morphological change). The effect of salinity on morphological evolution becomes less profound when the strength of tidal or fluvial forcing is dominant over the other. The effects of sediment type and flocculation, which are particularly important when salinity gradients are present, are also discussed. Overall, this study highlights that neglecting salinity to simulate long-term estuarine morphodynamics requires more careful justification, particularly when the environment is characterized by fine sediment types (favouring suspended transport), and relatively large river discharge and estuarine depth (favouring baroclinic effects). © 2020 John Wiley & Sons, Ltd.  相似文献   

3.
The morphodynamics of shallow, vertically well-mixed estuaries, characterised by tidal flats and deeper channels, have been investigated. This paper examines what contributes to flood/ebb-dominant sediment transport in localised regions through a 2D model study (using the TELEMAC modelling system). The Dyfi Estuary in Wales, UK has been used as a case study and, together with idealised estuary shapes, shows that shallow water depths lead to flood dominance in the inner estuary whilst tidal flats and deep channels cause ebb dominance in the outer estuary. For medium sands and with an artificially ‘flattened’ bathymetry (i.e. no tidal flats), the net sediment transport switches from ebb-dominant to flood-dominant where the parameter a/h (local tidal amplitude ÷ local tidally averaged water depth) exceeds 1.2. Sea level rise will reduce this critical value of a/h and also reduce the ebb-directed sediment transport significantly, leading to a flood-dominated estuarine system. A similar pattern, albeit with greater transport, was simulated with tidal flats included and also with a reduced grain size. This suggests that analogous classifications for flood/ebb asymmetry of the tide in estuaries as a whole may not represent the local sediment transport in sufficient detail. Through the Dyfi simulations, the above criterion involving a/h is shown to be complicated further by augmented flow past a spit at the estuary mouth which gives rise to a self-maintaining scour hole. Simulations of one year of bed evolution in an idealised flat-bottomed estuary, including tidal flow past a spit, recreate the flood/ebb dominance on either side of the spit and the formation of a scour hole in between. The erosion rate at the centre of the hole is reduced as the hole deepens, suggesting the establishment of a self-maintaining equilibrium state.  相似文献   

4.
Geomorphological characteristics of tidal basins control hydrodynamics and sediment transport potential within such basins, for example, by adjusting the balance in tidal asymmetry. In this study we examine the effects of entrance geometry on tidal velocity asymmetry, slack water asymmetry, bed shear stress patterns and hypsometric profile shapes by comparison of six shallow meso-tidal basins of Tauranga Harbour, New Zealand. Numerical model results show how tidal distortion increases with distance from a basin entrance. A simple ratio between basin width and entrance width defines levels of basin dilation. Sub-basins with a constricted geometry and deep entrance channels are associated with small bed shear stress values and high rates of flood-directed tidal velocity asymmetry in the sheltered basin centres, indicating a large potential for sediment deposition of larger particles. Moreover, slack water asymmetry within these basins is weakly ebb-directed, indicating a small potential for transport of fine sediments out of the basins. The constricted depositional basins are characterized by convex hypsometric profiles with elevated intertidal regions. Unconstricted geometries are associated with larger bed shear stress values and more ebb-directed tidal velocity asymmetry within basin centres, suggesting limited potential for overall sediment deposition. The slack tide duration asymmetry is weakly flood-dominant indicating that limited input of fine sediment into the basins is possible. The comparatively high-energy conditions within these exposed basins are associated with a less convex hypsometric intertidal profile. The ability to estimate tidal asymmetries is advantageous when developing management strategies related to ecosystem functioning, navigability or coastal protection in specific geomorphic settings. © 2019 John Wiley & Sons, Ltd.  相似文献   

5.
Dune bedforms and salt‐wedge intrusions are common features in many estuaries with sand beds, and yet little is known about the interactions between the two. Flow visualization with an echosounder and velocity measurements with an acoustic Doppler current profiler over areas of flat‐bed and sand dunes in the highly‐stratified Fraser River estuary, Canada, were used to examine the effect of dunes on interfacial mixing. As the salt‐wedge migrates upstream over the flat‐bed, mixing is restricted to the lower portion of the water column. However, as the salt‐wedge migrates into the dune field from the flat bed, there is a dramatic change in the flow, and large internal in‐phase waves develop over each of the larger dunes, with water from the salt‐wedge reaching the surface of the estuary. The friction Richardson number shows that bed friction is more important in interfacial mixing over the dunes than over the flat‐bed, and a plot of internal Froude Number versus obstacle (dune) height shows that the salt‐wedge flow over the dunes is mainly supercritical. Such bedforms can be expected to cause similar effects in interfacial mixing in other estuaries and sediment‐laden density currents, and may thus be influential in fluid mixing and sediment transport. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

6.
Tidal propagation in estuaries is affected by friction and fresh water discharge, besides changes in the depth and morphology of the channel. Main distortions imply variations in the mean water level and asymmetry. Tidal asymmetry can be important as a mechanism for sediment accumulation and turbidity maximum formation in estuaries, while mean water level changes can affect navigation depths. Data from several gauges stations from the Amazon estuary and the adjacent coast were analyzed and a 2DH hydrodynamic model was configured in a domain covering the continental shelf up to the last section of the river where the tidal signature is observed. Based on data, theoretical and numerical results, the various influences in the generation of estuarine harmonics are presented, including that of fresh water discharge. It is shown that the main overtide, M4, derived from the most important astronomic component in the Amazon estuary, M2, is responsible for the tidal wave asymmetry. This harmonic has its maximum amplitude at the mouth, where minimum depths are found, and then decreases while tide propagates inside the estuary. Also, the numerical results show that the discharge does not affect water level asymmetry; however, the Amazon river discharge plays an important role in the behavior of the horizontal tide. The main compound tide in Amazon estuary, Msf, generated from the combination of the M2 and S2, can be strong enough to provoke neap low waters lower than spring ones. The results show this component increasing while going upstream in the estuary, reaching a maximum and then slightly decaying.  相似文献   

7.
In tidal estuaries, quantifying freshwater discharge is still a difficult problem that has not yet been overcome due to the inherent difficulty in measuring and analysing the tidal discharge, especially during periods of low river flow. Because observations are often made in the stations further upstream, where the ratio of river to tidal discharge is large, it remains difficult to determine the discharge rate in the saline region. Freshwater discharge estimation is even more difficult in a branched estuary system having multiple diversion channels that connect with each other at a junction. To date, several methods have been developed for estimating freshwater discharge in estuaries. The most widely used are analytical and conceptual models that employ salinity as the principal trace and numerical simulations. However, these methods are very time consuming and costly as they require large sets of observations before the computations can take place. This paper presents a simple approach to investigating the discharge distribution over branched channels by considering the energy loss due to friction. We develop an analytical model that can obtain the discharge rate quantitatively at a junction where the main flow bifurcates into two branches. The model uses the bed roughness, tidal water level, and cross‐sectional profile under tidally averaged conditions as input data. Two selected estuarine systems in the Hiroshima delta in Japan and the Mekong delta in Vietnam have been investigated. Computations of the newly developed model show good agreement with earlier published results computed by sophisticated analytical and numerical models.  相似文献   

8.
This paper presents an approach to incorporate time‐dependent dune evolution in the determination of bed roughness coefficients applied in hydraulic models. Dune roughness is calculated by using the process‐based dune evolution model of Paarlberg et al. ( 2009 ) and the empirical dune roughness predictor of Van Rijn ( 1984 ). The approach is illustrated by applying it to a river of simple geometry in the 1‐D hydraulic model SOBEK for two different flood wave shapes. Calculated dune heights clearly show a dependency on rate of change in discharge with time: dunes grow to larger heights for a flood wave with a smaller rate of change. Bed roughness coefficients computed using the new approach can be up to 10% higher than roughness coefficients based on calibration, with the largest differences at low flows. As a result of this larger bed roughness, computed water depths can be up to 15% larger at low flow. The new approach helps to reduce uncertainties in bed roughness coefficients of flow models, especially for river systems with strong variations in discharge with time. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

9.
Contemporary hydrodynamics and morphological change are examined in a shallow microtidal estuary, located on a wave-dominated coast (Port Stephens, NSW, Australia). Process-based numerical modelling is undertaken by combining modules for hydrodynamics, waves, sediment transport and bathymetry updates. Model results suggest that the complex estuarine bathymetry and geometry give rise to spatial variations in the tidal currents and a marked asymmetry between ebb and flood flows. Sediment transport paths correspond with tidal asymmetry patterns. The SE storms significantly enhance the quantities of sediment transport, while locally generated waves by the westerly strong winds also are capable of causing sediment entrainment and contribute to the delta morphological change. The wave/wind-induced currents are not uniform with flow over shoals driven in the same direction as waves/winds while a reverse flow occurring in the adjacent channel. The conceptual sediment transport model developed in this study shows flood-directed transport occurs on the flood ramp while ebb-directed net transport occurs in the tidal channels and at the estuary entrance. Accretion of the intertidal sand shoals and deepening of tidal channels, as revealed by the model, suggest that sediment-infilling becomes advanced, which may lead to an ebb-dominated estuary. It is likely that a switch from flood- to ebb-dominance occurs during the estuary evolution, and the present-day estuary acts as a sediment source rather than sediment sink to the coastal system. This is conflictive to the expectation drawn from the estuarine morphology; however, it is consistent with previous research suggesting that, in an infilling estuary, an increase in build-up of intertidal flats/shoals can eventually shift an estuary towards ebb dominance. Thus, field data are needed to validate the result presented here, and further study is required to investigate a variety of estuaries in the Australian area.  相似文献   

10.
Tides are often considered to be the dominant hydrodynamic process within mesotidal estuaries although waves can also have a large influence on intertidal erosion rates. Here, we use a combination of hydrodynamic measurements and sediment deposition records to determine the conditions under which observed waves are ‘morphologically significant’, in which case they influence tidal and suspended sediment flux asymmetry and subsequently infilling over geomorphological timescales. Morphological significant conditions were evaluated using data from contrasting arms in a dendritic mesotidal estuary, in which the orientation of the arms relative to the prevailing wind results in a marked difference in wave conditions, deposition rates and morphology. By defining the morphological significance of waves as a product of the magnitude of bed shear stress and frequency of occurrence, even small (but frequently occurring) winds are shown to be capable of generating waves that are morphologically significant given sufficient fetch. In the arm in which fetch length is restricted, only stronger but rare storm events can influence sediment flux and therefore tides are more morphologically significant over longer timescales. Water depth within this mesotidal estuary is shown to be a critical parameter in controlling morphological significance; the rapid attenuation of short period waves with depth results in contrasting patterns of erosion occurring during neaps and accretion during springs. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

11.
This paper discusses a model which simulates dune development resulting from aeolian saltation transport. The model was developed for application to coastal foredunes, but is also applicable to sandy deserts with transverse dunes. Sediment transport is calculated using published deterministic and empirical relationships, describing the influence of meteorological conditions, topography, sediment characteristics and vegetation. A so-called adaptation length is incorporated to calculate the development of transport equilibrium along the profile. Changes in topography are derived from the predicted transport, using the continuity equation. Vegetation height is incorporated in the model as a dynamic variable. Vegetation can be buried during transport events, which results in important changes in the sediment transport rates. The sediment transport model is dynamically linked to a second-order closure air flow model, which predicts friction velocities over the profile, influenced by topography and surface roughness. Modelling results are shown for (a) the growth and migration of bare, initially sine-shaped dunes, and (b) dune building on a partly vegetated and initially flat surface. Results show that the bare symmetrical dunes change into asymmetric shapes with a slipface on the lee side. This result could only be achieved in combination with the secondorder closure model for the calculation of air flow. The simulations with the partly vegetated surfaces reveal that the resulting dune morphology strongly depends on the value of the adaptation length parameter and on the vegetation height. The latter result implies that the dynamical interaction between aeolian activity and vegetation (reaction to burial, growth rates) is highly relevant in dune geomorphology and deserves much attention in future studies. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

12.
Tidal sand dune dynamics is observed for two tidal cycles in the Arcachon tidal inlet, southwest France. An array of instruments is deployed to measure bathymetric and current variations along dune profiles. Based on the measurements, dune crest horizontal and vertical displacements are quantified and show important dynamics in phase with tidal currents. We observed superimposed ripples on the dune stoss side and front, migrating and changing polarity as tidal currents reverse. A 2D RANS numerical model is used to simulate the morphodynamic evolution of a flat non-cohesive sand bed submitted to a tidal current. The model reproduces the bed evolution until a field of sand bedforms is obtained that are comparable with observed superimposed ripples in terms of geometrical dimensions and dynamics. The model is then applied to simulate the dynamics of a field of large sand dunes of similar size as the dunes observed in situ. In both cases, simulation results compare well with measurements qualitatively and quantitatively. This research allows for a better understanding of tidal sand dune and superimposed ripple morphodynamics and opens new perspectives for the use of numerical models to predict their evolution.  相似文献   

13.
Dune mobility and vegetation cover in the Southwest Kalahari desert   总被引：1，自引：0，他引：1  
As part of a wider project investigating the palaeoenvironmental significance of partially vegetated linear dunes in the southwest Kalahari, data collected in the latter part of 1992 concerning dune movement and vegetation cover suggest that sediment transport is occurring on some dune surfaces, and that the majority of surface activity occurs on the crests and upper slopes of the dunes. The data suggest that the limiting variables on surface sediment movement vary on different parts of a dune. On interdunes and lower dune slopes the primary limiting variable is available wind energy, while on dune crests and upper slopes it is vegetation cover. Ground cover by litter has much greater importance in protecting the surface sediment from erosion than rooted vegetation. From individual data points, no evidence is found to support a threshold vegetation cover below which sediment movement occurs. Rather, a gradient of activity is suggested whereby a reduction in vegetation cover increases the potential for sediment movement and surface change. However, dunes with differing amounts of mean vegetation cover display differing degrees of surface activity, and at this scale, a vegetation cover threshold in the region of 14 per cent may be recognized.  相似文献   

14.
Both natural changes (e.g., tidal forcing from the ocean and global sea level rise) and human-induced changes (e.g., dredging for navigation, sand excavation, and land reclamation) exert considerable influences on the long-term evolution of tidal regimes in estuaries. Evaluating the impacts of these factors on tidal-regime shifts is particularly important for the protection and management of estuarine environments. In this study, an analytical approach is developed to investigate the impacts of estuarine morphological alterations (mean water depth and width convergence length) on tidal hydrodynamics in Lingdingyang Bay, Southeast China. Based on the observed tidal levels from two tidal gauging stations along the channel, tidal wave celerity and tidal damping/amplification rate of different tidal constituents are computed using tidal amplitude and phase of tidal constituents extracted from a standard harmonic analysis. We show that the minimum mean water depth for the whole estuary occurred in 2006, whereas a shift in tidal wave celerity for the M2 tide component occurred in 2009. As such, the study period (1990–2016) could be separated into pre-human (1990–2009) and post-human (2010–2016) phases. Our results show that the damping/amplification rate and celerity of the M2 tide have increased by 31% (from 7 to 9.2 m−1) and 28% (from 7 to 9 m·s−1) respectively, as a consequence of the substantial impacts of human interventions. The proposed analytical method is subsequently applied to analyse the historical development of tidal hydrodynamics and regime shifts induced by human interventions, thus linking the evolution of estuarine morphology to the dominant tidal hydrodynamics along the channel. The observed tidal regime shift is primarily caused by channel deepening, which substantially enlarged the estuary and reduced effective bottom friction resulting in faster celerity and stronger wave amplification. Our proposed method for quantifying the impacts of human interventions on tidal regime shifts can inform evidence-based guidelines for evaluating hydraulic responses to future engineering activities.  相似文献   

15.
Submarine dune dynamics are controlled by tidal currents and wind forces. According to the relative influence of these forces and the nature of dune sediment, different bedform behaviors can be observed. The footprint of the different hydrodynamic agents is recorded into the internal architecture of dunes. This paper is concerned with bedforms that compose the thick sediment wedge located in the eastern English Channel, off the Bay of Somme. This sedimentary archive constitutes an interesting feature to achieve a better understanding of seabed sediment dynamics and its timeline building stages. The dynamics of large submarine dunes, which are organized in fields, are studied thanks to bathymetric and seismic data over the periods 1937–1993 and 1993–2007. Dune morphology presents low lee and stoss side slopes (on average 8° and 3°, respectively) and dune migration rate is not very high. Dune movements are in the direction of residual tidal currents, i.e. toward the east, with mean migration rates around 0·8 to 5 ± 0·25 m yr?1 and up to 6·6 ± 0·7 m yr?1, respectively, at multi‐decennial and decennial time scales. The dune internal architecture is complex with superimposed eastward prograding units, displaying locally opposite progradation. Second‐order discontinuities (dip of 0·5°–4° perpendicular to dune crests) constitute dune master bedding. By counting the number of second‐order reflectors between 1937–1993 and 1993–2007, the formation periodicity of these bounding surfaces is estimated to range from 4 to 18 years. These time intervals coincide with the long‐term tidal cyclicities and also with the inter‐annual to decennial variability of storm activity in northern Europe. Two theories were made to interpret the dune internal structures: the second‐order surfaces are interpreted either as the depositional surfaces corresponding to the marks of weak energy periods (weak tidal and storm action), or as erosive surfaces due to an opposite direction of dune migration provoked temporarily by exceptional storms from the northeast. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

16.
Forecasts of water level during river floods require accurate predictions of the evolution of river dune dimensions, because the hydraulic roughness of the main channel is largely determined by the bed morphology. River dune dimensions are controlled by processes like merging and splitting of dunes. Particularly the process of dune splitting is still poorly understood and – as a result – not yet included in operational dune evolution models. In the current paper, the process of dune splitting is investigated by carrying out laboratory experiments and by means of a sensitivity analysis using a numerical dune evolution model. In the numerical model, we introduced superimposed TRIAS ripples (i.e. triangular asymmetric stoss side‐ripples) on the stoss sides of underlying dunes as soon as these stoss sides exceed a certain critical length. Simulations with the model including dune splitting showed that predictions of equilibrium dune characteristics were significantly improved compared to the model without dune splitting. As dune splitting is implemented in a parameterized way, the computational cost remains low which means that dune evolution can be calculated on the timescale of a flood wave. Subsequently, we used this model to study the mechanism of dune splitting. Literature showed that the initiation of a strong flow separation zone behind a superimposed bedform is one of the main mechanisms behind dune splitting. The flume experiments indicated that besides its height also the lee side slope of the superimposed bedform is an important factor to determine the strength of the flow separation zone and therefore is an important aspect in dune splitting. The sensitivity analysis of the dune evolution model showed that a minimum stoss side length was required to develop a strong flow separation zone. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

17.
Dipper Harbour Creek's lower reaches run through a narrow salt marsh on the Bay of Fundy, New Brunswick, Canada. This 2 km long section of the creek constitutes an extreme example of a tide-dominated estuary exhibiting strong downvalley morphology and sedimentology gradients. Dipper Harbour Creek drains a basin of roughly 8.8 km2, but except during the spring snowmelt freshet, tidal flow so overshadows freshwater flow within the salt marsh reach that the system essentially functions as a tidal creek. To identify and explain the main geomorphic processes controlling the creek system, records were collected in summer 1993 of tidal stage and velocity fluctuations, sand dune migration rates, bed material composition, channel cross-sectional geometry and channel sinuosity. Bed materials become progressively finer upvalley, with deposits of medium to coarse sands concentrated in the highly sinuous central reach of the creek during the summer. Current velocities within the creek are strongly flood-dominant, featuring a consistent low-stage peak in flood velocity, a secondary high-stage flood surge, and a weaker ebb peak occurring around bankfull stage. Under summer low freshwater discharge conditions, the predominant direction of bed sand transport is upvalley. The spring freshet, however, causes a major downvalley shift of sand deposits, suggesting a seasonal cycling of medium to coarse sands within Dipper Harbour Creek.  相似文献   

18.
The fate of mud in an estuary over an entire year was unravelled using complementary, independent, spatially explicit techniques. Sequential ERS-2 SAR and Envisat MERIS-FR data were used to derive synoptic changes in intertidal bottom mud and suspended particulate matter (SPM) in the top of the water column, respectively. These satellite data were combined with in situ measurements and with a high resolution three-dimensional cohesive sediment model, simulating mud transport, resuspension, settling and deposition under the influence of tides, wind, waves and freshwater discharge. The spatial distribution of both bottom mud and SPM as observed by in situ and satellite techniques was largely explained by modelled estuarine circulation, tidal and wind-induced variations in vertical mixing and horizontal advection. The three data sources also showed similar spring-neap and seasonal variations in SPM (all factor 1.5 to 2), but semi-diurnal tidal variations were underestimated by the model. Satellite data revealed that changes in intertidal bottom mud were spatially heterogeneous, but on average mud content doubled during summer, which was confirmed by in situ data. The model did not show such seasonal variation in bed sediment, suggesting that seasonal dynamics are not well explained by the physical factors presently implemented in the model, but may be largely attributed to other (internal) factors, including increased floc size in summer, temporal stabilisation of the sediment by microphytobenthos and a substantially lower roughness of the intertidal bed in summer as observed by the satellite. The effects of such factors on estuarine mud dynamics were evaluated.  相似文献   

19.
A macrotidal estuary with mountain streams(MEMS) is characterized by rapidly rising and falling flood peaks and large tidal ranges and is a typical estuary type with strong flow dynamics that is found worldwide. Understanding the morphodynamic evolution of MEMSs has great significance for river management. The roles of the mountain stream river flood process and macrotides on the morphology evolution still needs further quantitative study. Taking the Oujiang Estuary as an example, the evolution ...  相似文献   

20.
The response of the Yalu River Estuary to human activities was investigated.Changes of sediment dynamics during the past 10 years were explored through hydrodynamic calculation,as well as heavy mineral...  相似文献   

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

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