首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 70 毫秒
1.
The dispersion and deposition of particulate organic matter from a fish cage located in an idealized curved channel with a 90° bend are studied for different horizontal grid resolutions. The model system consists of a three-dimensional, random-walk particle tracking model coupled to a terrain-following ocean model. The particle tracking model is a Lagrangian particle tracking simulator which uses the local flow field, simulated by the ocean model, for advection of the particles and random walk to simulate the turbulent diffusion. The sinking of particles is modeled by imposing an individual particle settling velocity. As the homogeneous water flows through the bend in the channel, the results show that a cross-channel secondary circulation is developed. The motion of this flow is similar to a helical motion where the water in the upper layers moves towards the outer bank and towards the inner bank in the lower layers. The intensity of the secondary circulation will depend on the viscosity scheme and increases as the horizontal grid resolution decreases which significantly affects the distribution of the particles on the seabed. The presence of the secondary circulation leads to that most of the particles that settle, settle close to the inner bank of the channel.  相似文献   

2.
Flocculation has an important impact on particle trapping in estuarine turbidity maximum (ETM) through associated increases in particle settling velocity. To quantify the importance of the flocculation processes, a size-resolved flocculation model is implemented into an ocean circulation model to simulate fine-grained particle trapping in an ETM. The model resolves the particle size from robust small flocs, about 30 μm, to very large flocs, over 1000 μm. An idealized two-dimensional model study is performed to simulate along-channel variations of suspended sediment concentrations driven by gravitational circulation and tidal currents. The results indicate that the flocculation processes play a key role in generating strong tidal asymmetrical variations in suspended sediment concentration and particle trapping. Comparison with observations suggests that the flocculation model produces realistic characteristics of an ETM.  相似文献   

3.
海潮对卫星重力场恢复的影响   总被引:7,自引:2,他引:7       下载免费PDF全文
本文讨论了海潮对卫星重力测量的影响问题. 首先介绍了海潮对卫星重力测量影响的基本理论;采用FES02和TPXO6海潮模型计算了海潮负荷对卫星重力结果前60阶的影响;并用两个模型之间的差异作为海潮模型精度的估计量,据此计算了海潮模型误差对卫星重力结果的影响. 与GRACE恢复的重力场精度的比较说明:海潮对重力场40阶以下的影响都超过了目前重力场恢复精度;尽管由于卫星测高技术的发展,海潮模型的精度有了很大的提高,但目前的全球海潮模型用于GRACE重力场恢复的前12阶的改正还是不够精确. 另外,我们也利用中国东海和南海潮汐资料以及FES02海潮模型讨论了中国近海潮汐效应对GRACE观测的影响. 结果说明该影响与海潮模型的误差相当. 这反映了当前海潮模型的不确定度,因此通过结合全球验潮站资料有望提高海潮对卫星重力测量的改正精度.  相似文献   

4.
《Continental Shelf Research》2005,25(9):1071-1080
Visible-band satellite pictures of the Irish Sea reveal the presence of isolated areas of enhanced turbidity which are geographically fixed and present all year, although they are more strongly marked in winter. The positions of these maxima coincide with areas of fast tidal currents: some of the dissipated tidal energy is used to raise fine sediments into suspension, producing the turbidity. However, there is no obvious source of fine sediment at the maxima, and without a source they would be expected to diffuse away, down the turbidity gradient. Their continued presence is therefore a puzzle, and it has proved difficult to reproduce these features in numerical models. Recent observations, presented here, suggest a possible mechanism by which isolated turbidity maxima may be maintained. These measurements show that the gradients of concentration on the sides of the turbidity maximum are in opposite senses for different particle sizes, and that there is a flux of fine particles out of the maximum and one of larger particles into it. A mechanism which can explain this observation, and which can also explain the continued presence of turbidity maxima isolated from a local source is as follows: the high turbulent energy levels at the centre of the patch tear flocculated particles apart. These fine particles then diffuse away down the gradient of fine particles and to areas of lower energy. Here they aggregate to form larger particles which diffuse back down the gradient of large particles towards the centre of the turbidity maximum, where they are torn up and the cycle continues. The source of material for the turbidity maxima is therefore larger flocculated particles in the surrounding water. This idea is tested quantitatively with an analytical solution to the steady-state diffusion equation incorporating aggregation and dis-aggregation of particles as simple functions of turbulent energy. The solution shows that isolated maxima of fine suspended particles can be maintained at regions of high turbulence without the necessity to invoke a local (non-sustainable) source. Within the maximum, strong vertical mixing lifts the slow settling fine particles to the surface to produce an isolated surface turbidity maximum as observed. We conclude that for models to successfully produce turbidity maxima in the presence of diffusion they should incorporate at least two particle size classes and aggregation and dis-aggregation of particles according to the local level of turbulence.  相似文献   

5.
Seasonal circulation of the Bohai Sea (BS) in 1992 was investigated using Lagrangian particle tracking method. The hydrography of the BS was simulated based on an unstructured grid, finite-volume, three-dimensional primitive equation ocean model. With the use of the unstructured triangular grid, the model can easily fit the irregular coastal boundary of the BS. The simulated tides, tidal current, and thermohaline field agreed well with the observations. The transport of particles has three-dimensional structure in the BS. Compared with central Bohai and Bohai Strait, the differences of particles’ transportation between surface and bottom layer in three bays are small. The circulation in the summer is stronger than that in the winter, with the average residual velocity in the surface layer being about 3.7 cm/s during the summer while only 1.8 cm/s during the winter. Using the same model, several well-designed numerical experiments were performed to investigate the effect of oceanic tide, river discharge, wind stress, and thermal stratification on the circulation. It is shown that winds play an important role in the circulation of the BS during both the winter and the summer. Density circulation is important during the summer; however, it is negligible during the winter. River runoff only affects the area around the river mouth. Compared with wind and thermohaline effect, the contribution of tides is small during the summer, and the circulation under only M2 tidal constituent could not reflect the actual circulation of the BS.  相似文献   

6.
The MARS-3D model in conjunction with the particle tracking module Ichthyop is used to study circulation and tracer dynamics under a variety of forcing conditions in the eastern English Channel, and in the Boulogne-sur-Mer harbour (referred to hereafter as BLH). Results of hydrodynamic modelling are validated against the tidal gauge data, VHF radar surface velocities and ADCP measurements. Lagrangian tracking experiments are performed with passive particles to study tracer dispersal along the northern French coast, with special emphasis on the BLH. Simulations revealed an anticyclonic eddy generated in the harbour at rising tide. Tracers, released during flood tide at the Liane river mouth, move northward with powerful clockwise rotating current. After the high water, the current direction changes to westward, and tracers leave the harbour through the open boundary. During ebb tide, currents convergence along the western open boundary but no eddy is formed, surface currents inside the harbour are much weaker and the tracer excursion length is small. After the current reversal at low water, particles are advected shoreward resulting in a significant increase of the residence time of tracers released during ebb tide. The effect of wind on particle dispersion was found to be particularly strong. Under strong SW wind, the residence time of particles released during flood tide increases from 1.5 to 6 days. For release during ebb tide, SW wind weakens the southward tidally induced drift and thus the residence time decreases. Similar effects are observed when the freshwater inflow to the harbour is increased from 2 to 10 m3/s during the ebb tide flow. For flood tide conditions, the effect of freshwater inflow is less significant. We also demonstrate an example of innovative coastal management targeted at the reduction of the residence time of the pathogenic material accidentally released in the harbour.  相似文献   

7.
Based on pressure tide-gauge observations, sea-level records are derived for ten sites along the coast of West Greenland. The ocean tidal signal is extracted by a harmonic tidal analysis. The accuracy of the determined tidal constants is discussed in detail. The tides account for 85% of the observed sea-level standard deviation. The tide gauge records reveal significant shallow-water tidal effects, in particular compound and overtide amplitudes reaching 5 cm. The propagation of the tidal waves into the fjords depends strongly on local conditions and is in some cases accompanied by an amplification of the tidal amplitudes. The observed tidal signals are compared to the predictions of the global ocean tide model FES2004. At the outer coast, a good agreement is found. Inside the fjords, however, the model performs worse and tide gauge observations may still be indispensable when accurate tidal signals are required.  相似文献   

8.
Jouanneau  Nicolas  Sentchev  Alexei  Dumas  Franck 《Ocean Dynamics》2013,63(11):1321-1340

The MARS-3D model in conjunction with the particle tracking module Ichthyop is used to study circulation and tracer dynamics under a variety of forcing conditions in the eastern English Channel, and in the Boulogne-sur-Mer harbour (referred to hereafter as BLH). Results of hydrodynamic modelling are validated against the tidal gauge data, VHF radar surface velocities and ADCP measurements. Lagrangian tracking experiments are performed with passive particles to study tracer dispersal along the northern French coast, with special emphasis on the BLH. Simulations revealed an anticyclonic eddy generated in the harbour at rising tide. Tracers, released during flood tide at the Liane river mouth, move northward with powerful clockwise rotating current. After the high water, the current direction changes to westward, and tracers leave the harbour through the open boundary. During ebb tide, currents convergence along the western open boundary but no eddy is formed, surface currents inside the harbour are much weaker and the tracer excursion length is small. After the current reversal at low water, particles are advected shoreward resulting in a significant increase of the residence time of tracers released during ebb tide. The effect of wind on particle dispersion was found to be particularly strong. Under strong SW wind, the residence time of particles released during flood tide increases from 1.5 to 6 days. For release during ebb tide, SW wind weakens the southward tidally induced drift and thus the residence time decreases. Similar effects are observed when the freshwater inflow to the harbour is increased from 2 to 10 m3/s during the ebb tide flow. For flood tide conditions, the effect of freshwater inflow is less significant. We also demonstrate an example of innovative coastal management targeted at the reduction of the residence time of the pathogenic material accidentally released in the harbour.

  相似文献   

9.
Non-linear tidal dynamics are investigated in a network that consists of a semi-enclosed main channel and a secondary channel at an arbitrary position. The water motion, governed by the one-dimensional shallow water equations, is forced by an incoming tidal wave. Solutions are obtained with the method of characteristics. The overall aim is to quantify and understand the spatial structure of different tidal harmonics (the principal tide and its non-linear overtides) and of tidal asymmetry for both the vertical and the horizontal tide in the main channel for different locations of the secondary channel. This is of practical interest in the context of possible construction of secondary channels to reduce tidal range in estuaries. Moreover, tidal asymmetry is an important factor in driving net sediment transport. Analysis of the different tidal harmonics shows that their characteristics are similar to those obtained with an earlier linear model. In particular, amplitudes of the harmonics are reduced landward of the secondary channel if the latter is positioned less than a quarter wavelength of the respective tidal wave away from the landward boundary. Thus, the distortions of the tide due to the presence of the secondary channel are generated locally and afterwards propagate through the network. Tidal asymmetry is quantified by examining tidal range, flood-to-ebb ratio and the duration of the falling tide and the duration between maximum flood and maximum ebb. A spatial non-uniform reduction in tidal range is observed that shows very localised increase and decrease depending on the position of the secondary channel. The changes in the velocity characteristics induce changes in net sediment transport. It turns out that the direction of the peak current, derived from the flood-to-ebb ratio, is not sensitive to the position of the secondary channel, whereas the duration between flood and ebb can change from more to less than half the tidal cycle. However, the changes in the velocity asymmetries are confined to a small region.  相似文献   

10.
中国大陆精密重力潮汐改正模型   总被引:12,自引:4,他引:8       下载免费PDF全文
利用理论和实验重力固体潮模型,充分考虑全球海潮和中国近海潮汐的负荷效应,建立了中国大陆的精密重力潮汐改正模型.结果表明,采用不同的固体潮模型会对重力潮汐结果产生相对变化幅度小于0.06%的差异;在沿海地区海潮负荷的影响约为整个潮汐的4%,而中部地区约为1%,其中中国近海潮汐模型的影响约占整个海潮负荷的10%,内插或外推潮波的负荷约占海潮负荷的3%.通过比较实测的重力数据表明,本文给出的重力潮汐改正模型的精度远远优于0.5×10-8 m·s-2,说明了本文构建的模型的实用性,可为中国大陆高精度重力测量提供有效参考和精密的改正模型.  相似文献   

11.
Sediment particles are often colonized by biofilm in a natural aquatic ecological system, especially in eutrophic water body. A series of laboratory experiments on particle size gradation, drag coefficient and settling velocity were conducted after natural sediment was colonized by biofilm for 5, 10, 15 and 20 days. Particle image acquisition, particle tracking techniques of Particle Image Velocimetry and Particle Tracking Velocimetry were utilized to analyze the changes of these properties. The experimental results indicate that the size gradation, the drag force exerted on bio-particles, and the settling velocity of bio-particles underwent significant change due to the growth of biofilm onto the sediment surface. The study proposes a characteristic particle size formula and a bio-particle settling velocity formula based on the regression of experiment results, that the settling velocity is only 50% to 60% as the single particle which has the same diameter and density. However, biofilm growth causes large particle which the settling velocities are approximately 10 times larger than that of primary particles. These results may be specifically used in the low energy reservoir or lake environment.  相似文献   

12.
First, we investigated some aspects of tsunami–tide interactions based on idealized numerical experiments. Theoretically, by changing total ocean depth, tidal elevations influence the speed and magnitude of tsunami waves in shallow regions with dominating tidal signals. We tested this assumption by employing a simple 1-D model that describes propagation of tidal waves in a channel with gradually increasing depth and the interaction of the tidal waves with tsunamis generated at the channel's open boundary. Important conclusions from these studies are that computed elevations by simulating the tsunami and the tide together differ significantly from linear superposing of the sea surface heights obtained when simulating the tide and the tsunami separately, and that maximum tsunami–tide interaction depends on tidal amplitude and phase. The major cause of this tsunami–tide interaction is tidally induced ocean depth that changes the conditions of tsunami propagation, amplification, and dissipation. Interactions occur by means of momentum advection, bottom friction, and variable water flux due to changing total depth and velocity. We found the major cause of tsunami–tide interactions to be changing depth. Secondly, we investigate tsunami–tide interactions in Cook Inlet, Alaska, employing a high-resolution 2-D numerical model. Cook Inlet has high tides and a history of strong tsunamis and is a potential candidate for tsunami impacts in the future. In agreement with previous findings, we find that the impacts of tsunamis depend on basin bathymetries and coastline configurations, and they can, in particular, depend on tsunami–tide interactions. In regions with strong tides and tsunamis, these interactions can result in either intensification or damping of cumulative tsunami and tide impacts, depending on mean basin depth, which is regulated by tides. Thus, it is not possible to predict the effect of tsunami–tide interaction in regions with strong tides without making preliminary investigations of the area. One approach to reduce uncertainties in tsunami impact in regions with high tides is to simulate tsunamis together with tidal forcing.  相似文献   

13.
Satellite altimetry and GRACE observations carry both the signature of ocean tides and have in general complementary potential to resolve tidal constituents. It is therefore straightforward to perform a combined estimation of a global ocean tide model based on these two data sources. The present paper develops and applies a three step procedure for generating such a combined ocean tide model. First, the processing of multi-mission altimetry data is described along with the harmonic analysis applied to derive initially a pure empirical ocean tide model. Then the capability of GRACE to sense particular tidal constituents is elaborated and an approach to estimate tidal constituents from GRACE is outlined. In a third step a combination strategy with optimal stochastic data treatment is developed and applied to the altimetry-only tide model EOT08a and four years of GRACE observations, leading to the combined model EOT08ag. The differential contributions of GRACE to EOT08ag remain small and are mainly concentrated to the Arctic Ocean, an area with little or poor altimetry data. In comparison with other tide models, EOT08ag is validated by K-band range residuals, the impact on gravity field modelling and on precise orbit determination and by variance reduction of crossover differences. None of these comparison exhibits a significant improvement over the altimetry-only tide model except for a few areas above 60°N. Overall the improvements of the combination remain small and appear to stay below the current GRACE baseline accuracy.  相似文献   

14.
The Sable Gully is a broad deep underwater canyon located to the east of Sable Island on the edge of the Scotian Shelf. Being the home of many marine species including the endangered Northern Bottlenose Whale, the Gully was designated as a marine protected area (MPA) in 2004. Better understanding of physical environmental conditions over this MPA is needed for sustainable ecosystem management. In this study, a multi-nested ocean circulation model and a particle tracking model are used to examine the three-dimensional (3D) circulation and movement of particles carried passively by the flow over the Sable Gully. The 3D circulation model is driven by tides, wind, and surface heat/freshwater fluxes. The model performance is assessed by comparing the results with the previous numerical tidal results and current meter observations made in the Gully. The simulated tidal circulation over the Gully and adjacent waters is relatively strong on shallow banks and relatively weak on the continental slope. Below the depth of the Gully rim ( ~ 200 m), the tidal currents are constrained by the thalweg of the Gully and amplified toward the Gully head. The simulated subtidal circulation in the Gully has a complex spatial structure and significant seasonal variability. The simulated time-dependent 3D flow fields are then used in a particle tracking model to study the particle movements, downstream and upstream areas, and residence time of the Gully. Based on the movements of particles released at the depth of the Gully rim and tracked forward in time, the e-folding residence time is estimated to be about 7 and 13 days in February and August 2006, respectively. The Gully flanks are identified as high retention areas with the typical residence time of 10 and 20 days in February and August 2006, respectively. Tracking particles with and without tides reveals that tidal circulation reduces the value of residence time in the Gully, particularly along the Gully flanks.  相似文献   

15.
Hydrographic conditions, and particularly current speeds, have a strong influence on the management of fish cage culture. These hydrodynamic conditions can be used to predict particle movement within the water column and the results used to optimise environmental conditions for effective site selection, setting of environmental quality standards, waste dispersion, and potential disease transfer. To this end, a 3D hydrodynamic model, MOHID, has been coupled to a particle tracking model to study the effects of mean current speed, quiescent water periods and bulk water circulation in Mulroy Bay, Co. Donegal Ireland, an Irish fjard (shallow fjordic system) important to the aquaculture industry. A Lagangrian method simulated the instantaneous release of “particles” emulating discharge from finfish cages to show the behaviour of waste in terms of water circulation and water exchange. The 3D spatial models were used to identify areas of mixed and stratified water using a version of the Simpson-Hunter criteria, and to use this in conjunction with models of current flow for appropriate site selection for salmon aquaculture.The modelled outcomes for stratification were in good agreement with the direct measurements of water column stratification based on observed density profiles. Calculations of the Simpson-Hunter tidal parameter indicated that most of Mulroy Bay was potentially stratified with a well mixed region over the shallow channels where the water is faster flowing. The fjard was characterised by areas of both very low and high mean current speeds, with some areas having long periods of quiescent water. The residual current and the particle tracking animations created through the models revealed an anticlockwise eddy that may influence waste dispersion and potential for disease transfer, among salmon cages and which ensures that the retention time of waste substances from cages is extended.The hydrodynamic model results were incorporated into the ArcViewTM GIS system for visualisation and interrogation of results and to allow effective holistic environmental management and site selection for fish cage aquaculture.  相似文献   

16.
Traditionally, ocean tides have been modeled in frequency domain with a forcing from selected tidal constituents. It is a natural approach; however, it implicitly neglects non-linearities of ocean dynamics. An alternative approach is time-domain modeling with a forcing given by the full lunisolar potential, i.e., all tidal waves are a priori included. This approach has been applied in several ocean tide models; however, some challenging tasks still remain, for example, assimilation of satellite altimetry data. In this paper, we introduce the assimilative scheme applicable in a time-domain model, which is an alternative to existing techniques used in assimilative ocean tide models. We present results from DEBOT, a global barotropic ocean tide model, which has two modes: DEBOT-h, a purely hydrodynamical mode, and DEBOT-a, an assimilative mode. The accuracy of DEBOT in both modes is assessed through a series of tests against tide gauge data which demonstrate that DEBOT is comparable to state-of-the-art global ocean tide models for major tidal constituents. Furthermore, as signals of all tidal frequencies are included in DEBOT, we also discuss modeling of minor tidal constituents and non-linear compound tides. Our modeling approach can be useful for those applications where the frequency domain approach is not suitable.  相似文献   

17.
Environmental management of marine fish culture in Hong Kong   总被引:4,自引:0,他引:4  
Marine fish farming is an important commercial practice in Hong Kong. Marine fish farms located in eutrophic coastal waters often face the threat of severe dissolved oxygen depletion associated with algal blooms and red tides. On the other hand, mariculture activities also contribute to pollution. The sustainable management of mariculture requires proper siting of the fish farms and stocking density control. Both of these are related to the carrying capacity of the water body concerned, which is mainly governed by its flushing characteristics. A simple method to determine the carrying capacity of a fish farm has been developed by using three-dimensional (3D) hydrodynamic modelling and its effective coupling with a diagenetic water quality model. A systematic methodology using numerical tracer experiments has been developed to compute the tidal flushing in a fish farm. The flushing time is determined from the results of a numerical tracer experiment using robust 3D hydrodynamic and mass transport models. A unit tracer concentration is initially prescribed inside the region of interest and zero elsewhere; the subsequent mass transport and the mass removal process are then tracked. The fish farms are usually situated in well-sheltered shallow embayments and may not connect directly to the open water. It is found that it is necessary to define both "local" and "system-wide" flushing times to represent the effectiveness of the mass exchange with the surrounding water body and the open sea respectively. A diagenetic water quality model simulating the sediment-water-pollutant interaction is employed to address the response of the water column and the benthic layer to pollution discharges. With the flushing rate reliably computed, the carrying capacity of the fish farm can be determined in terms of key water quality parameters: chlorophyll-a, dissolved oxygen, organic nitrogen and potential lowest dissolved oxygen level on a day of negligible photosynthetic production. The predictions are well-supported by field data.  相似文献   

18.
The seasonal cycle of the main lunar tidal constituent M 2 is studied globally by an analysis of a high-resolution ocean circulation and tide model (STORMTIDE) simulation, of 19 years of satellite altimeter data, and of multiyear tide-gauge records. The barotropic seasonal tidal variability is dominant in coastal and polar regions with relative changes of the tidal amplitude of 5–10 %. A comparison with the observations shows that the ocean circulation and tide model captures the seasonal pattern of the M 2 tide reasonably well. There are two main processes leading to the seasonal variability in the barotropic tide: First, seasonal changes in stratification on the continental shelf affect the vertical profile of eddy viscosity and, in turn, the vertical current profile. Second, the frictional effect between sea-ice and the surface ocean layer leads to seasonally varying tidal transport. We estimate from the model simulation that the M 2 tidal energy dissipation at the sea surface varies seasonally in the Arctic (ocean regions north of 60°N) between 2 and 34 GW, whereas in the Southern Ocean, it varies between 0.5 and 2 GW. The M 2 internal tide is mainly affected by stratification, and the induced modified phase speed of the internal waves leads to amplitude differences in the surface tide signal of 0.005–0.0150 m. The seasonal signals of the M 2 surface tide are large compared to the accuracy demands of satellite altimetry and gravity observations and emphasize the importance to consider seasonal tidal variability in the correction processes of satellite data.  相似文献   

19.
Although there have been studies on the tide in convergent bay (or estuary), the tide change in terms of phase speed, amplitude, and phase difference between elevation and tidal current from a coastal ocean to a convergent bay has not been clearly shown so far. This study systematically examines the change of tidal wave characteristics from the eastern Yellow Sea to the Asan Bay, a strongly convergent bay on the west coast of Korea, using observations and an analytical model. As the tidal wave propagates from the eastern Yellow Sea into the Asan Bay, the phase speed, amplitude, and phase difference between elevation and tidal current increase along the channel. Such a phenomenon represents a unique example of tide change from a coastal ocean to a convergent bay, indicating dominance of convergence over friction in the Asan Bay. Both analytically computed tidal amplitude and travelling time compare well with observations. In the Asan Bay, the influence of the reflected wave is only felt in the upper one fifth of the bay and is almost unperceivable in the rest of the bay. The analytical analyses presented in this paper are particularly useful for understanding the relative importance of channel convergence, bottom friction, and reflected wave on the tidal characteristics change along the channel and the proposed method could be applicable to other estuaries.  相似文献   

20.
A Lagrangian model is used to evaluate source regions of particles collected in the sediment traps at the DYFAMED(Dynamique des Flux Atmosphériques en Méditerranée) station by tracking particles backwards from March 1 to August 31, 2001. The analysis suggests that source regions depend on the flow fields, the settling speed of the particles, and the deployment depths of the traps. Monthly variation is observed in the distribution patterns of source regions, which is caused by the currents. The source regions are located around the traps and up to hundreds of kilometers away. As the settling speed increases with the particle diameters, the distance to the source regions decreases. The vertical flux can be approximately estimated in 1D for the particles with diameters larger than 500 μm. Furthermore, traps moored at various depths at the DYFAMED can collect particles that originated from different regions in the Ligurian Sea.  相似文献   

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

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