首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
文章使用2019年7月5日—20日在珠江河口伶仃洋定点连续观测的海表面叶绿素a质量浓度、海表面气温、气压、风速、风向、海表温度、盐度、流速、流向、遥感降雨量数据和中等分辨率成像光谱仪可见光波段影像, 利用小波分析和集成经验模态分解方法分析了观测期间内伶仃洋海表面叶绿素a的时间变化特征及其影响因子。分析结果表明, 观测期间海水表层叶绿素a质量浓度的变化范围为0.44~1.75µg·L-1, 平均值为0.80µg·L-1, 其变化周期主要为6h、12h和24h。其与相对应周期的潮流存在明显的相位关系, 并且在降雨后两者的相位关系发生了转换。7月5日—12日, 叶绿素a与潮流基本呈反相位关系, 涨急时叶绿素a质量浓度低, 落急时叶绿素a质量浓度较高, 浓度相差约为0.3µg·L-1。珠江流域在7月8日—13日发生了一次强降雨过程, 降雨前后海水表层叶绿素a质量浓度在6h、12h和24h周期波段的振幅由0.02~0.09µg·L-1增加到0.15µg·L-1左右。同时, 降雨对珠江河口的叶绿素a质量浓度造成了一个持续80h的增加过程, 浓度增加了0.3µg·L-1。发生降雨后, 7月13日—20日期间潮流滞后于叶绿素a约6h, 水位最高时叶绿素a质量浓度最低, 水位最低时叶绿素a质量浓度最高。由以上结果可以看出, 降雨不仅引起了河口区叶绿素a质量浓度的增加, 还造成了叶绿素a和潮流间相位关系的转换。  相似文献   

2.
Experiments were conducted in an oscillatory water tunnel to investigate what effects temporal inertia has on the resistance of a granular medium. The flow law governing stationary porous media flow is reasonably well accepted and understood but the effects of unsteady flow have commonly been neglected. The present research was designed to assess the magnitude of the acceleration effects on media of uniformly packed spheres of equal diameter and on one sample of randomly placed stone.Oscillatory flow tests were made in a large oscillatory flume tunnel with periods varying from 3 to 12 seconds. The influence of properties of the medium (grain size and porosity) were tested by using spheres with two different diameters and packing each size sphere in different geometric arrangements. Tests made on a stone sample provided a qualitative assessment of the effects of more random material properties.For the experiments described in this paper, the Forchheimer unsteady-stationary flow law described the oscillatory measurements well when velocities and energy losses were maximum. Empirical coefficients determined from steady-stationary flow were generally found to apply to the unsteady flow, however some evidence of dependency on the period of oscillation was noted.  相似文献   

3.
In this paper, the hydrodynamic coefficients of a horizontal semi-immersed cylinder in steady current and oscillatory flow combining with constant current are obtained via forced oscillation experiments in a towing tank. Three non-dimensional parameters (Re, KC and Fr) are introduced to investigate their effects on the hydrodynamic coefficients. The experimental results show that overtopping is evident and dominates when the Reynolds number exceeds 5×105 in the experiment. Under steady current condition, overtopping increases the drag coefficient significantly at high Reynolds numbers. Under oscillatory flow with constant current condition, the added mass coefficient can even reach a maximum value about 3.5 due to overtopping while the influence of overtopping on the drag coefficient is minor.  相似文献   

4.
Non-linear loads on a fixed body due to waves and a current are investigated. Potential theory is used to describe the flow, and a three-dimensional (3D) boundary element method (BEM), combined with a time-stepping procedure, is used to solve the problem. The exact free-surface boundary conditions are expanded about the still-water level by Taylor series so that the solution is evaluated on a time-invariant geometry. A formulation correct to second order in the wave steepness and to first order in the current speed is used. Numerical results are obtained for the first-order and the second-order oscillatory forces and for the second-order mean force on a fixed vertical circular cylinder in waves and a current. The second-order oscillatory forces on the body in waves and current are new results, while the remaining force components are verified by comparison with established numerical and analytical models. It is shown that the current can have a significant influence on the forces, and especially on the amplitude of the second-order oscillatory component.  相似文献   

5.
The drag and lift force are measured on circular cylinders fitted with end plates in a wind tunnel. The gap between the cylinder and the wall, G, the thickness of the turbulent boundary layer along the wall, δ, and the Reynolds number, Re, are varied in the following ranges: 0 < G/D < 2, 0.12 < δ/D < 0.97 and 4.8 × 104 Re 3 × 105. The lift and drag coefficients are presented in terms of a new variable G/δ.

It is found that the lift coefficient is governed by the gap to diameter ratio G/D while the drag coefficient is dominated by the ratio of gap to thickness of the boundary layer, G/δ.  相似文献   


6.
Computer simulations of steady plus low-amplitude oscillatory flow about a circular cylinder are reported at a fixed Reynolds number of 150 based on the steady component. The conventional Keleugan–Carpenter number based on the oscillatory component is fixed at π/5. The oscillation frequency is varied so as to study a wide spectrum of flows where inertial forces dominate at one end and viscous drag forces at the other as a function of the modified Keleugan–Carpenter number. The hydrodynamic force on the cylinder in-line with the flow direction is represented by Morison's equation and an extended version with three terms. The drag and inertia coefficients in Morison's equation are determined by least-squares fits to data directly computed from integration of skin friction and pressure distributions around the periphery of the cylinder. The root-mean-square value of the residue of reconstructed minus directly-computed forces varies between 2 and 41% depending on the flow parameters. Comparable results can be obtained with a semi-theoretical approach using inviscid inertia and quasi-steady viscous drag terms. Physical explanations for the variation of the force coefficients are provided and implications for pertinent flow–structure interactions are discussed.  相似文献   

7.
Kinematics of extreme waves in deep water   总被引：2，自引：0，他引：2  
The velocity profiles under crest of a total of 62 different steep wave events in deep water are measured in laboratory using particle image velocimetry. The waves take place in the leading unsteady part of a wave train, focusing wave fields and random wave series. Complementary fully nonlinear theoretical/numerical wave computations are performed. The experimental velocities have been put on a nondimensional form in the following way: from the wave record (at a fixed point) the (local) trough-to-trough period, TTT and the maximal elevation above mean water level, ηm of an individual steep wave event are identified. The local wavenumber, k and an estimate of the wave slope, ε are evaluated from ω2/(gk)=1+ε2, where ω=2π/TTT and g denotes the acceleration of gravity. A reference fluid velocity, is then defined. Deep water waves with a fluid velocity up to 75% of the estimated wave speed are measured. The corresponding kηm is 0.62. A strong collapse of the nondimensional experimental velocity profiles is found. This is also true with the fully nonlinear computations of transient waves. There is excellent agreement between the present measurements and previously published Laser Doppler Anemometry data. A surprising result, obtained by comparison, is that the nondimensional experimental velocities fit with the exponential profile, i.e. eky, y the vertical coordinate, with y=0 in the mean water level.  相似文献   

8.
Wedge entry into initially calm water   总被引：2，自引：0，他引：2  
This paper presents results of calculations based on the Cauchy's theorem method of Vinje and Brevig1 for the two-dimensional entry of wedges of various angles into initially calm water. The problem has a long history which is briefly reviewed in the introduction, and significant progress has been made with both linear theories (valid for low entry speed) and with theories which treat the free surface conditions exactly but with the assumptions of zero gravity and constant speed of entry. This simplifies the problem to one which is self-similar in dimensionless space variables ξ = x/vt and η = y/vt and this has a number of consequences. For wedges with half-angles up to about 45° and with high entry speeds, the numerical approach, which includes gravity, validates these assumptions and the agreement between both free surface displacements and pressure distributions on the wetted wedge surface is excellent except in the region of the jet of fluid which rises up the side of the wedge. Because the potential flow initial value problem is singular at the intersection of the free surface and wedge surface, exact numerical resolution of the jet is not possible. Nevertheless, the rest of the fluid motion is insensitive to the treatment of the jet, which itself may be calculated quite realistically. Of particular interest (but little practical relevance) is the pressure on the upper part of the wedge surface (in the jet region) which according to self-similar theories is very small but positive, but which is calculated to be small but negative by the numerical scheme. This effect, which is enhanced when gravity is included, is insensitive to the numerical resolution of the jet and suggests that the jet may separate from the wedge surface, the new intersection point being where the pressure vanishes on the wedge surface. A modified numerical scheme allows this to happen and the results are in qualitative agreement with the experiments by Greenhow and Lin.2

The numerical method presented here is extremely versatile and a number of other effects may be explored. Examples of transient motion, non-constant speed of entry, oblique entry and complete penetration of the surface so that a cavity is formed behind the wedge are presented.  相似文献   


9.
为了探究北部湾养殖区域香港牡蛎(Crassostrea hongkongensis)体内的异养细菌和弧菌数量及其耐药概况变化, 对不同养殖场牡蛎体内的异养细菌进行分离培养, 并统计其数量, 通过药敏纸片扩散等方法研究了细菌的耐药状况。结果显示: 牡蛎在高死亡率养殖环境中体内的异养细菌[(8.6±0.4)×106CFU·g-1]和弧菌[(9.5±0.4)×105CFU·g-1]数量较高, 在中死亡率环境中体内的异养细菌[(6.9±0.2)×106CFU·g-1]和弧菌[(4.5±0.6)×105CFU·g-1]数量次之, 在低死亡率养殖环境中体内的异养细菌[(3.3±0.1)×106CFU·g-1]和弧菌[(2.5±0.6)×105CFU·g-1]数量最低。耐药细菌主要为革兰氏阴性菌, 对β-内酰胺类(青霉素)、糖肽类(万古霉素)的耐药率较高, 对四环素类(四环素、多西环素)的耐药率次之, 对氨基糖苷类(链霉素、庆大霉素、妥布霉素、新霉素)、大环内酯类(红霉素)、喹诺酮类(诺氟沙星、环丙沙星、氧氟沙星、恩诺沙星)的耐药率较低。在高死亡率环境中牡蛎体内的多重耐药菌占79.7%, 其耐药谱型(48种)较广; 在中度死亡率环境中牡蛎体内的多重耐药菌占66.2%, 其耐药谱型为30种; 在低死亡率环境中牡蛎体内的多重耐药菌占58.4%, 其耐药谱型为17种。本文探究了牡蛎死亡率与其体内异养细菌数量和细菌耐药性的关系, 结果显示牡蛎在高死亡率环境中体内的耐药细菌数量多、耐药谱型较广, 低死亡率环境中牡蛎体内的耐药细菌数量较少, 异养细菌数量与牡蛎死亡率呈正相关关系, 两者相关系数为0.996。  相似文献   

10.
海洋中光后向散射系数的变化包含了浮游植物生物量的信息, 可应用于卫星遥感和光学剖面观测平台获取海洋中大时空尺度-高分辨率剖面的浮游植物生物量变化特征。本文选取了琼东上升流影响下生物—光学变异性较为显著的海域, 基于2013年航次实测数据, 建立了颗粒物后向散射系数(bbp)与叶绿素a浓度(Chl a)间的区域性关系模型。模型假定颗粒物后向散射系数由不随叶绿素浓度变化的固定背景值, 以及较大粒级(>2μm)和pico级(微微型, <2μm)两类浮游植物的后向散射贡献累加所得。采集的数据集进行了模型检验, 结果表明, 模型能很好地模拟琼东海域水体的bbp与Chl a间的变化趋势, 性能优于常用的幂函数关系模型, 尤其在低叶绿素浓度范围, 很好地解决幂函数显著低估的现象; 琼东海域的bbp和Chl a关系存在显著的水层变化, 底层后向散射固定背景值显著高于上层水体背景值, 表明底层受上升流的影响, 水体中不随Chl a共变的颗粒物浓度增大, 其后向散射相应增强; 叶绿素最大层的后向散射固定背景值显著低于上层其他水体的固定背景值, 后向散射固定背景值的贡献百分比约为21%~35%; 随着叶绿素浓度增大, 较大粒级的浮游植物对颗粒物后向散射系数的贡献也显著增大, 可达到50%以上, pico级浮游植物贡献稳定在40%附近。本研究的结果将为琼东海域浮游植物生物量的光学遥感、生物地球化学过程研究提供更为精确的区域性模型和基础支撑数据。  相似文献   

11.
In this article a method to estimate the mean upcrossing intensity, μ(u), of a stochastic process is proposed. It is assumed that the stochastic process is a sum of a Gaussian process and a second-order correction term. The method is based on the two-dimensional Saddlepoint approximation. The accuracy of the method is tested on processes having analytical solutions for μ(u). Numerical examples are given where the stochastic process represents (i) the horizontal response of a floating offshore structure in a Gaussian sea, and (ii) the response of a structure subjected to a Gaussian wind velocity process. In addition, the estimates are compared to empirical upcrossing intensities of simulated responses. For case (ii), the obtained μ(u) estimates are compared to estimates obtained by numerical integration.  相似文献   

12.
The authors estimated the effects of a cooling water discharge from a coastal power plant on underwater light levels using 6 years of integrated hourly irradiance measurements which were taken at the sea bottom in water depths of 10–14m. They applied a before-after-control-impact design (BAC1) to the irradiance recordings. The basic BAC1 design estimates the power plant effect as the difference between the means of impact-control differences observed in separate periods after and before the power plant began operation. This design subtracts out constant spatial differences and uniform temporal changes, leaving only the time-by-location interaction which is the power plant effect. They modified the basic BACI analysis in several ways to deal with high natural variability, serial correlation in the data, and the effects of current direction. The daily irradiance data were divided into two sets of days in which the impact stations were downcurrent or upcurrent from the power plant discharge. For downcurrent days, the authors found significant (p < 0·0.016) reductions in irradiance ranging from 20 to 28% at three impact stations. For the upcurrent days, they found insignificant increases of about 14%.  相似文献   

13.
《Coastal Engineering》2006,53(11):897-913
For the general purposes of morphodynamic computations in coastal zones, simple formula-based models are usually employed to evaluate sediment transport. Sediment transport rates are computed as a function of the bottom shear stress or the near bed flow velocity and it is generally assumed that the sediment particles react immediately to changes in flow conditions. It has been recognized, through recent laboratory experiments in both rippled and plane bed sheet flow conditions that sediment reacts to the flow in a complex manner, involving non-steady processes resulting from memory and settling/entrainment delay effects. These processes may be important in the cross-shore direction, where sediment transport is mainly caused by the oscillatory motions induced by surface short gravity waves.The aim of the present work is to develop a semi-unsteady, practical model, to predict the total (bed load and suspended load) sediment transport rates in wave or combined wave-current flow conditions that are characteristic of the coastal zone. The unsteady effects are reproduced indirectly by taking into account the delayed settling of sediment particles. The net sediment transport rates are computed from the total bottom shear stress and the model takes into account the velocity and acceleration asymmetries of the waves as they propagate towards the shore.A comparison has been carried out between the computed net sediment transport rates with a large data set of experimental results for different flow conditions (wave-current flows, purely oscillatory flow, skewed waves and steady currents) in different regimes (plane bed and rippled bed) with fine, medium and coarse uniform sand. The numerical results obtained are reasonably accurate within a factor of 2. Based on this analysis, the limits and validity of the present formulation are discussed.  相似文献   

14.
A Deep Draft Semi-submersible (DDS) under certain flow conditions could be subjected to Vortex-Induced Motions (VIM), which significantly influences the loads on and life fatigue of the moorings and the risers. To investigate the VIM of a DDS with four rectangular section columns in waves coupled with a uniform current, a numerical study using the computational fluid dynamics (CFD) method was conducted. The issues of the VIM of multi-column floaters can be conveniently converted to the issues of oscillating cylinders in fluid cross flows. This paper looks into the CFD numerical simulation of infinite cylinders having rectangular sections in a two-dimensional sinusoidal time-dependent flow field coupled with a uniform current. The resulted hydrodynamic forces and motion responses in different oscillatory flows plus currents both aligned in the same direction for the incidence of 135° of the DDS relative to the flow are compared with the ones in current only cases. The results show that the VIM response of this geometric arrangement of a DDS with four rectangular columns in a current combined with oscillatory flows is more evident than that in the current only case. The oscillatory flows and waves have the significant influence on the VIM response, forces and trajectory, in-plane motions of the DDS.  相似文献   

15.
Boundary roughness and bedforms in the surf zone   总被引：1，自引：0，他引：1  
Hydrodynamical models of the nearshore system frequently assume that a single friction coefficient is sufficient to represent flow conditions at a point in the surf zone. Furthermore, models attempting to relate bed configuration to surf zone flows have relied primarily upon the wave orbital velocity as an indicator of potential bedforms, and thus as the control on boundary roughness. The data presented here point out potential errors arising from either of these approaches. The results of a field experiment conducted at Wendake Beach, Ontario, show that at a single location in an active surf zone, the Darcy-Weisbach friction coefficient, f, varied by approximately 250% (in this case between 0.016 and 0.041).

It is also shown that existing bedform models, based upon primary wave motions alone, do not accurately predict conditions at this study site. For a relatively constant wave orbital velocity and velocity asymmetry, it is found that changes in bed roughness, as a result of bedform development, are reflected mainly in the vertical profile of the long-shore current velocity. A sequence of bedforms, from oscillatory ripples through flat bed, is inferred from the data, and found to be supported by diver observations and preserved primary sedimentary structures.  相似文献   


16.
The dimensions of sand ripples in full-scale oscillatory flows   总被引：1，自引：0，他引：1  
New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D50 ≥ 0.30 mm and 3D ripples occur when D50 ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.  相似文献   

17.
在三亚海域设置3个站位, 分别于各站位连续采集12个月表层水样。利用流式细胞仪进行浮游病毒及浮游细菌丰度的测定, 并对两者及其与环境因子之间的相关性进行研究, 同时对不同站位之间进行差异性分析。结果表明, 调查海域浮游病毒丰度(平均7.63×106viruses·mL-1)高于浮游细菌丰度(平均1.52×106cells·mL-1)。浮游病毒及浮游细菌丰度在三亚河口最高, 且有明显的季节变化, 冬春季高于夏秋季, 并且其与鹿回头半岛西侧、小东海之间均有极显著差异(P<0.01), 不同类群病毒代表的宿主类群也有所不同。调查海域总体水平浮游病毒丰度与浮游细菌丰度显著正相关(r=0.800, P<0.01); 叶绿素(Chl a)和氮盐(NO- 2、NO- 2+NO- 3、NH+ 4)是影响两者的关键因子。  相似文献   

18.
The flow around a cylinder, surface mounted or partly buried has been examined at Reynolds number of 1.3×104 and 2.6×104 by using the flow visualization method in the case of steady current. The lengths of the upstream and downstream separation regions were determined and the characteristics of the flow area were discussed for the cylinder having different burial depth to the diameter ratio (G/D). For comparsion, some numerical results obtained from the application of the FLUENT, computational fluid dynamics (CFD) software package have also been presented. Both experimental and numerical results indicated that the lengths of the separation regions near the upstream and downstream of the cylinder decreased with the increasing burial ratio (G/D=0, 0.20, 0.50). The results have significant implications both for the illustration of the potential for using CFD software and highlighting the need for data of two dimensional flow around the buried cylinder.  相似文献   

19.
海表二氧化碳分压(pCO2)是指海洋表层水和大气之间的二氧化碳(CO2)交换处于动态平衡时CO2的含量, 是描述海-气CO2交换的一个主要因子。本文利用2008—2014年覆盖南海大部分海域的海表pCO2观测资料, 结合现场海表温度和海表盐度以及卫星观测的叶绿素a数据, 构建了基于多元线性回归方法的分区域反演模型。模型在水深浅于30m的区域均方根误差为5.3μatm, 其余海区均方根误差为10.8μatm, 与前人基于个别航次的有限区域反演结果的均方根误差相当。利用该模型公式和HYbrid Coordinate Ocean Model(HYCOM)再分析海表温、盐数据及MODIS-Aqua卫星观测的叶绿素a数据进行反演, 得到了时空分辨率为5'×5'的2004—2016年的逐月南海海表pCO2数据。该数据能较好地反映南海海表pCO2在海表温度影响下, 春夏高、秋冬低的季节变化特征, 与前人基于航次观测的研究结果相似, 表明反演模型具有较高的可信度。进一步分析发现, 南海及邻近海域平均海表pCO2具有显著的准十年振荡特征: 2012年附近出现了极小值, 之前表现为降低的趋势, 之后略有升高的趋势。受海表pCO2的影响, 南海海盆平均海-气CO2通量在2012年之前出现了显著降低的趋势, 表明南海释放到大气中的CO2减少, 并在2007年之后的冬季出现了负值(从碳源变为碳汇), 2012年之后变化较为平缓。热带太平洋年代际振荡引起的南海区域海表盐度变化是造成海表pCO2及海-气CO2通量准十年变化的主要原因。分区分析的结果表明, 南海北部海表pCO2变化最为显著, 在南海海表pCO2的季节和准十年变化中都起到非常重要的作用。  相似文献   

20.
层移输沙是海岸带泥沙运动的主要形式之一,其垂向悬沙浓度分布规律的研究一直是海岸工程关心的重点。一般情况下,经典的纯扩散模型被用来描述和解释悬沙浓度的试验数据,该模型认为周期平均悬沙浓度主要由参考浓度、泥沙沉速和泥沙扩散系数确定。泥沙扩散系数可以由泥沙沉速和悬沙浓度的垂向梯度反演得到。既往研究大多直接给出泥沙扩散系数的结果,对于不同反演计算方法间结果差别的研究较少。本研究汇总了已有振荡流层移输沙试验数据,采用曲线拟合方法和直接差分方法计算了相应的泥沙扩散系数,研究表明两种方法得到的计算结果在垂向位置z 0.15 m处差异不大,随着垂向位置的升高,差分方法的计算结果略微大于拟合方法。考虑到拟合方法可以得到连续的泥沙扩散系数垂向分布,本研究推荐使用幂函数形式的曲线拟合方法求解悬移泥沙扩散系数。基于此,对比分析了层移输沙悬沙层泥沙扩散系数随泥沙粒径、振荡流周期、均方根流速和振荡流类型等物理参数的变化规律。在纯振荡流层移输沙条件下,泥沙扩散系数随泥沙粒径的增大而增大,而振荡流周期和均方根流速几乎不影响泥沙扩散系数。在振荡流和定常流共同作用下,泥沙扩散系数受振荡流周期和定常流流速的影响,泥沙扩散系数随着振荡流周期的增大或定常流流速的减小而增大。  相似文献   

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

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