首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Tests of models for high-frequency seafloor backscatter   总被引:3,自引:0,他引:3  
The interaction of high-frequency sound with the seafloor is inherently a stochastic process. Inversion techniques must, therefore employ good stochastic models for bottom acoustic scattering. An assortment of physical models for bottom backscattering strength is tested by comparison with scattering strength data obtained at 40 kHz at three shallow water sites spanning a range of sediment types from fine silt to coarse sand. These acoustic data are accompanied by sediment physical property data obtained by core sample analysis and in situ probes. In addition, stereo photography was used to measure the power spectrum of bottom relief on centimeter scales. These physical data provided the inputs needed to test the backscatter models, which treat scattering from both the rough sediment-water interface and the sediment volume. For the three sites considered here, the perturbation model for scattering from a slightly rough fluid seafloor performs well. Volume scattering is predicted to be weak except at a site having a layer of methane bubbles  相似文献   

2.
多波束反向散射强度数据处理研究   总被引:8,自引:5,他引:8  
在探讨多波束测深系统反向散射强度与海底底质类型的关系基础上,研究影响反向散射强度的各种因素,主要分析了海底地形起伏、中央波束区反射信号对反向散射强度的影响,并给出了消除这些影响的方法;将处理后的“纯”反向散射强度数据镶嵌生成海底声像图,为海底底质类型划分以及地貌解译提供了基础数据和辅助判读依据.  相似文献   

3.
High-resolution multi-beam/single-beam bathymetric data and seismic profiling data from the latest surveys are used to map and interpret the detailed seafloor geomorphology of the western region of the North Yellow Sea (NYS), China. The mapping area covers 156 410 km2, and incorporates a flat shelf plain, subaqueous accumulation shoals, tidal scouring troughs, and tidal sand ridge groups. Offshore areas with water depths less than 50 m in the western region of the NYS are mainly covered by thick, loose sediments, forming wide spread accumulation geomorphological features; these include the Liaodong Peninsula subaqueous accumulation system containing shoals and rugged scouring troughs, and the large mud wedge of the Shandong Peninsula. In the central part of the NYS, there is a relatively flat residual shelf plain with coarser sediment deposits. This flat shelf plain has a water depth larger than 50 m and a thin layer of sediment, on which there is a large pockmark field caused by seafloor seepage. These geomorphological structures indicate that modern sedimentary processes are the main driving force controlling the sculpture of the current seafloor surface landform. Extensive strong tidal current systems and abundant sediment sources provide the critical external forces and essential conditions for the formation of seafloor geomorphology. The tectonic basement controls the macroscopic morphological shape of the NYS, but is reflected very little in the seafloor geomorphic elements. Our results provide a detailed seafloor geomorphological map of the western region of the NYS, an area that has not previously mapped and also provide a scientific framework for further research into offshore seafloor geomorphology, shelf sedimentary processes, and submarine engineering construction in this region.  相似文献   

4.
Hydrographic quality bathymetry and quantitative acoustic backscatter data are now being acquired in shallow water on a routine basis using high frequency multibeam sonars. The data provided by these systems produce hitherto unobtainable information about geomorphology and seafloor geologic processes in the coastal zone and on the continental shelf.Before one can use the multibeam data for hydrography or quantitative acoustic backscatter studies, however, it is essential to be able to correct for systematic errors in the data. For bathymetric data, artifacts common to deep-water systems (roll, refraction, positioning) need to be corrected. In addition, the potentially far greater effects of tides, heave, vessel lift/squat, antenna motion and internal time delays become of increasing importance in shallower water. Such artifacts now cause greater errors in hydrographic data quality than bottom detection. Many of these artifacts are a result of imperfect motion sensing, however, new methods such as differential GPS hold great potential for resolving such limitations. For backscatter data, while the system response is well characterised, significant post processing is required to remove residual effects of imaging geometry, gain adjustments and water column effects. With the removal of these system artifacts and the establishment of a calibrated test site in intertidal regions (where the seabed may be intimately examined by eye) one can build up a sediment classification scheme for routine regional seafloor identification.When properly processed, high frequency multibeam sonar data can provide a view of seafloor geology and geomorphology at resolutions of as little as a few decimetres. Specific applications include quantitative estimation of sediment transport rates in large-scale sediment waves, volume effects of iceberg scouring, extent and style of seafloor mass-wasting and delineation of structural trends in bedrock. In addition, the imagery potentially provides a means of quantitative classification of seafloor lithology, allowing sedimentologists the ability to examine spatial distributions of seabed sediment type without resorting to subjective estimation or prohibitively expensive bottom-sampling programs. Using Simrad EM100 and EM1000 sonars as an example, this paper illustrates the nature and scale of possible artifacts, the necessary post-processing steps and shows specific applications of these sonars.  相似文献   

5.
A new highly precise source of data has recently become available using multibeam sonar systems in hydrography. Multibeam sonar systems can provide hydrographic quality depth data as well as high-resolution seafloor sonar images. We utilize the seafloor backscatter strength data of each beam from multibeam sonar and the automatic classification technology so that we can get the seafloor type identification maps. In this article, analyzing all kinds of error effects in backscatter strength, data are based on the relationship between backscatter strength and seafloor types. We emphasize particularly analyzing the influences of local bottom slope and near nadir reflection in backscatter strength data. We also give the correction algorithms and results of these two influent factors. After processing the raw backscatter strength data and correcting error effects, we can get processed backscatter strength data which reflect the features of seafloor types only. Applying the processed backscatter strength data and mosaicked seafloor sonar images, we engage in seafloor classification and geomorphy interpretation in future research.  相似文献   

6.
Processing Multibeam Backscatter Data   总被引:1,自引:0,他引:1  
A new highly precise source of data has recently become available using multibeam sonar systems in hydrography. Multibeam sonar systems can provide hydrographic quality depth data as well as high-resolution seafloor sonar images. We utilize the seafloor backscatter strength data of each beam from multibeam sonar and the automatic classification technology so that we can get the seafloor type identification maps. In this article, analyzing all kinds of error effects in backscatter strength, data are based on the relationship between backscatter strength and seafloor types. We emphasize particularly analyzing the influences of local bottom slope and near nadir reflection in backscatter strength data. We also give the correction algorithms and results of these two influent factors. After processing the raw backscatter strength data and correcting error effects, we can get processed backscatter strength data which reflect the features of seafloor types only. Applying the processed backscatter strength data and mosaicked seafloor sonar images, we engage in seafloor classification and geomorphy interpretation in future research.  相似文献   

7.
Recent developments in seafloor imaging and mapping techniques greatly improved our capability of identifying marine geohazards affecting continental margins. Geomorphic features can be detected in great detail by high-resolution multibeam imaging and regarded as geohazard indicators; the most common include slide scars and deposits, canyon headscarps and steep erosional flanks, fault-related seafloor unevenness, mud volcanoes, pockmarks, gravity flow deposits, erosional scours and bedforms indicating sediment mobility at diverse temporal/spatial scale. These processes are widespread on Italian continental margins and are potential indicators of geohazard for human settlements and infrastructures in the offshore and coastal zones. The national Project MaGIC (Marine Geohazards along the Italian Coasts) aims at documenting potential geohazards based on the acquisition of high-resolution multibeam bathymetry and on the production of maps of the geohazard-related geomorphic features for most of the Italian continental margins. With reference to this issue, we discuss some of the most frequent problems dealing with reconnaissance, interpretation and cartographic representation of geohazard-related geomorphic features at a regional scale.  相似文献   

8.
Here we apply quantitative technique to describe the seafloor seepages based on the multi-beam backscatter and bathymetric investigations to characterize the pockmark morphology. The variable seafloor backscatter strength for coarser seafloor sediments are related to the diagenesis derived from biodegraded seepages. In this regard, box counting method is used to estimate ‘fractal dimension’ for backscatter imagery data of 398 blocks. These blocks are further sub-grouped into six classes depending on the spread of pockmark related seepages. The study area lies 102 km west off Marmagao along the central west coast of India which contains pre-dominantly (70%) gas-charged sediments. Comparison between the estimated self-similar fractals reveals that there is approximately 97% correlation between the box (Dbox) and information (Dinfo) dimensions. Box dimension–derived fractal dimension values, suggest that the seepages are more along the fault trace in deeper waters, in comparison to sparsely distributed shallow water seepages. Besides, this poor seepage is confined within the smooth to moderately rough seafloor. It is established that the high backscatter strength along the upper slope of the pockmark region having higher fractal dimensions reflects multifractal behavior of seepage distribution. Entire area indicates patchy seepage patterns as supported by estimated fractal values showing intermittent fluctuations, which emphasizes non-linear behavior. Estimated self organizing criticality (SOC) parameters for six representative blocks reveal that the nature of pockmark, fault trace, sediment nature coupled with slumping of pockmark’s wall, sediment movement due to bottom currents are controlling the dynamic balance in the area seepage system. Further, our study emphasizing the multifractal behavior of seepage blocks, clearly depicts the drift in the seepage pattern.  相似文献   

9.
Application of quantitative angular backscatter modelling to manganese nodule-bearing areas of the Central Indian Ocean Basin (CIOB) has been initiated at NIO during the year 1998. Studies were aimed to establish the suitability of seafloor backscattering in delineating seafloor parameters characteristic of nodule-rich sediments. In this paper, processed Hydrosweep multi-beam backscatter data from 45 spot locations in the CIOB (where nodule samples are available) were analysed to estimate seafloor and sediment volume roughness parameters. The application of a composite roughness model to a nodule-bearing region (6,600 km2) of the CIOB, to determine seafloor interface roughness parameters from a multi-beam backscatter dataset, shows only four power law sets. The results attest 80% of the nodule-bearing seafloor to be smooth in terms of interface roughness parameters at micro-topographic level. The sediment volume roughness parameters are dominant only in 29% of the smooth interface roughness sites. This indicates that 51% of the seafloor area possesses negligible (interface and volume) roughness. A critical analysis using pseudo-side-scan records from 12 selected locations in the study area affirms the combined importance of the seafloor interface and sediment volume roughness parameters for precise determination of manganese nodule abundance.  相似文献   

10.
Selection of a set of dominant echo features to classify seafloor sediments using a multilayer perceptron neural network is investigated at two acoustic frequencies (33 and 210 kHz). Several sets of inputs with different combinations of two, three, four, five, and six echo features are exploited with three-layer neural networks. The performances of the networks are analyzed to assess the most discriminating set of echo features for classification of seafloor sediments. The results of the overall average performances reveal that backscatter strength and time spread are the two most important echo features at 33 kHz, whereas backscatter strength has higher discriminating characteristics at 210 kHz for seafloor sediment classification. In addition, a set of four echo features consisting of backscatter strength, time-spread, statistical skewness, and Hausdroff dimension gives the highest success at both the acoustic frequencies.  相似文献   

11.
Inhomogeneous substrate analysis using EM300 backscatter imagery   总被引:2,自引:0,他引:2  
Backscatter reflectivity from multibeam echo-sounders provides a powerful tool to efficiently characterize seafloor substrates. A comprehensive EM300 bathymetric and backscatter survey has been completed of Cook Strait, in central New Zealand. This paper presents a detailed analysis of the realtime corrections applied to the raw EM300 multibeam data and additional corrections required to compute angular variations of the backscatter strength. The corrections, including the local absorption coefficient, the influence of seafloor topography and sound refraction in the water column, are determined for different Cook Strait seafloor substrates. Modifying MB-System software code, we extracted the backscatter signal parameters in order to quantify the raw backscatter strength and apply additional processing. Profiles of backscatter strength versus incidence angle were computed for a variety of sites characterized by flat seafloor and homogeneous substrates, and for which ground-truth data were available. For each homogeneous site, different but characteristic backscatter profiles are observed that can be interpreted in terms of sediment facies. To analyze heterogeneous substrates, we present a statistical technique, based on a 3-dimensional distribution of (incidence angle, backscatter strength) couples that preserves geological information of the substrate components. This analysis, using backscatter data acquired on a submarine volcano, north of New Zealand, clearly differentiates soft sediments and lava flows within a heterogeneous substrate.  相似文献   

12.
While the average seafloor backscatter strength within a narrow range of grazing angles can be used as a first-order classification tool, this technique often fails to distinguish seafloors of known differing geological character. In order to resolve such ambiguities, it is necessary to examine the variation in backscatter strength as a function of grazing angle. For this purpose, a series of multiply overlapping GLORIA sidescan sonar images (6.5 kHz) have been obtained in water depths ranging from 1000 to 2500 m. To constrain the placement of acoustic backscatter measurements and to measure the true impinging angle of the incident wave, the corresponding seafloor was simultaneously surveyed using the Seabeam multibeam system. As a result of the multiple overlap, the angular response of seafloor backscatter strength may be derived for regions much smaller than the swath width. By using the derived angular response of seafloor backscatter strength in regions for which sediment samples exist, an empirical seafloor classification scheme is proposed based on the shape, variance, and magnitude of the angular response. Because of the observed variability in the shape of the angular response with differing seafloor types, routine normalization of single-pass swath data to an equivalent single grazing angle image cannot be achieved. As a result, for the case of single-pass surveys, confident seafloor classification may only be possible for regions approaching the scale of the swath width  相似文献   

13.
A brief review of the published evidence of current deposits around Italy is the occasion to test the robustness of matching bottom current velocity models and seafloor morphologies to identify contourite drifts not yet documented. We present the result of the regional hydrodynamic model MARS3D in the Northern Tyrrhenian and Ligurian Sea with horizontal resolution of 1.2 km and 60 levels with focus on bottom current: data are integrated over summer and winter 2013 as representative of low and high intensity current conditions.The Eastern Ligurian margin is impacted by the Levantine Intermediate Water (LIW) with modeled mean velocity of bottom current up to 20 cm s−1 in winter 2013 and calculated bottom shear stress exceeding 0.2 N m−2 in water depth of 400–800 m. By crossing this information with seafloor morphology and geometry of seismic reflections, we identify a sediment drift formerly overlooked at ca 1000 m water depth. The Portofino separated mounded drift has a maximum thickness of at least 150 m and occurs in an area of mean current velocity minimum. Independent evidence to support the interpretation include bottom current modelling, seafloor morphology, seismic reflection geometry and sediment core facies. The adjacent areas impacted by stronger bottom currents present features likely resulted from bottom current erosion such as a marine terrace and elongated pockmarks.Compared to former interpretation of seafloor morphology in the study area, our results have an impact on the assessment of marine geohazards: submarine landslides offshore Portofino are small in size and coexist with sediment erosion and preferential accumulation features (sediment drifts) originated by current-dominated sedimentary processes. Furthermore, our results propel a more general discussion about contourite identification in the Italian seas and possible implications.  相似文献   

14.
This paper examines the potential for remote classification of seafloor terrains using a combination of quantitative acoustic backscatter measurements and high resolution bathymetry derived from two classes of sonar systems currently used by the marine research community: multibeam echo-sounders and bathymetric sidescans sonar systems. The high-resolution bathymetry is important, not only to determine the topography of the area surveyed, but to provide accurate bottom slope corrections needed to convert the arrival angles of the seafloor echoes received by the sonars into true angles of incidence. An angular dependence of seafloor acoustic backscatter can then be derived for each region surveyed, making it possible to construct maps of acoustic backscattering strength in geographic coordinates over the areas of interest. Such maps, when combined with the high-resolution bathymetric maps normally compiled from the data output by the above sonar systems, could be very effective tools to quantify bottom types on a regional basis, and to develop automatic seafloor classification routines.  相似文献   

15.
Using automated supervised segmentation of multibeam backscatter data to delineate seafloor substrates is a relatively novel technique. Low-frequency multibeam echosounders (MBES), such as the 12-kHz EM120, present particular difficulties since the signal can penetrate several metres into the seafloor, depending on substrate type. We present a case study illustrating how a non-targeted dataset may be used to derive information from multibeam backscatter data regarding distribution of substrate types. The results allow us to assess limitations associated with low frequency MBES where sub-bottom layering is present, and test the accuracy of automated supervised segmentation performed using SonarScope® software. This is done through comparison of predicted and observed substrate from backscatter facies-derived classes and substrate data, reinforced using quantitative statistical analysis based on a confusion matrix. We use sediment samples, video transects and sub-bottom profiles acquired on the Chatham Rise, east of New Zealand. Inferences on the substrate types are made using the Generic Seafloor Acoustic Backscatter (GSAB) model, and the extents of the backscatter classes are delineated by automated supervised segmentation. Correlating substrate data to backscatter classes revealed that backscatter amplitude may correspond to lithologies up to 4 m below the seafloor. Our results emphasise several issues related to substrate characterisation using backscatter classification, primarily because the GSAB model does not only relate to grain size and roughness properties of substrate, but also accounts for other parameters that influence backscatter. Better understanding these limitations allows us to derive first-order interpretations of sediment properties from automated supervised segmentation.  相似文献   

16.
Multibeam echosounders (MBES) have become a widely used acoustic remote sensing tool to map and study the seafloor, providing co-located bathymetry and seafloor backscatter. Although the uncertainty associated with MBES-derived bathymetric data has been studied extensively, the question of backscatter uncertainty has been addressed only minimally and hinders the quantitative use of MBES seafloor backscatter. This paper explores approaches to identifying uncertainty sources associated with MBES-derived backscatter measurements. The major sources of uncertainty are catalogued and the magnitudes of their relative contributions to the backscatter uncertainty budget are evaluated. These major uncertainty sources include seafloor insonified area (1–3 dB), absorption coefficient (up to >?6 dB), random fluctuations in echo level (5.5 dB for a Rayleigh distribution), and sonar calibration (device dependent). The magnitudes of these uncertainty sources vary based on how these effects are compensated for during data acquisition and processing. Various cases (no compensation, partial compensation and full compensation) for seafloor insonified area, transmission losses and random fluctuations were modeled to estimate their uncertainties in different scenarios. Uncertainty related to the seafloor insonified area can be reduced significantly by accounting for seafloor slope during backscatter processing while transmission losses can be constrained by collecting full water column absorption coefficient profiles (temperature and salinity profiles). To reduce random fluctuations to below 1 dB, at least 20 samples are recommended to be used while computing mean values. The estimation of uncertainty in backscatter measurements is constrained by the fact that not all instrumental components are characterized and documented sufficiently for commercially available MBES. Further involvement from manufacturers in providing this essential information is critically required.  相似文献   

17.
Both acoustic and sediment surveys were carried out in the Broughton Archipelago, British Columbia, in order to map a former aquaculture site and calibrate acoustic surveys with georeferenced sediment properties. The acoustic surveys included EM3000 Multibeam (including backscatter) and QTC VIEW™ (Series IV) technologies, while the geotechnical survey entailed Van Veen grab sampling of surface sediments and associated analyses. The two acoustic technologies were consistent in their ability to identify distinct regions of seafloor characterized by rock outcrops, consolidated substrates, or gel-mud depositional fields. Both multibeam backscatter data and QTC VIEW™ number-coded classifications were extracted across a range of circular areas located at each georeferenced sampling station (radii: 2, 3, 4, 5, 8, 12, 16, 20 m). Statistical correlations were observed between backscatter and certain geotechnical properties, such as sediment porosity, sediment grain size fractions (<2 μm, silt content), and particulate sulfur concentration. The areal resolution of backscatter extraction was explored in terms of determining a sensitive calibration technique between backscatter and sediment properties. In general the highest r2 values between backscatter and sediment variables were observed across extraction radii between 8 and 20 m. Such groundtruthing techniques could be used to interpolate seafloor characteristics between sampling stations and provide a steering tool for sampling designs associated with benthic monitoring programs.  相似文献   

18.
Acoustic backscatter images of the seafloor obtained with sidescan sonar systems are displayed most often using a flat bottom assumption. Whenever this assumption is not valid, pixels are mapped incorrectly in the image frame, yielding distorted representations of the seafloor. Here, such distortions are corrected by using an appropriate representation of the relief, as measured by the sonar that collected the acoustic backscatter information. In addition, all spatial filtering operations required in the pixel relocation process take the sonar geometry into account. Examples of the process are provided by data collected in the Northeastern Pacific over Fieberling Guyot with the SeaMARC II bathymetric sidescan sonar system and the Sea Beam multibeam echo-sounder. The nearly complete (90%) Sea Beam bathymetry coverage of the Guyot serves as a reference to quantify the distortions found in the backscatter images and to evaluate the accuracy of the corrections performed with SeaMARC II bathymetry. As a byproduct, the processed SeaMARC II bathymetry and the Sea Beam bathymetry adapted to the SeaMARC II sonar geometry exhibit a 35m mean-square difference over the entire area surveyed.On leave at the Naval Research Laboratory, Code 7420, Washington D.C. 20375-5350.  相似文献   

19.
This study applies three classification methods exploiting the angular dependence of acoustic seafloor backscatter along with high resolution sub-bottom profiling for seafloor sediment characterization in the Eckernförde Bay, Baltic Sea Germany. This area is well suited for acoustic backscatter studies due to its shallowness, its smooth bathymetry and the presence of a wide range of sediment types. Backscatter data were acquired using a Seabeam1180 (180 kHz) multibeam echosounder and sub-bottom profiler data were recorded using a SES-2000 parametric sonar transmitting 6 and 12 kHz. The high density of seafloor soundings allowed extracting backscatter layers for five beam angles over a large part of the surveyed area. A Bayesian probability method was employed for sediment classification based on the backscatter variability at a single incidence angle, whereas Maximum Likelihood Classification (MLC) and Principal Components Analysis (PCA) were applied to the multi-angle layers. The Bayesian approach was used for identifying the optimum number of acoustic classes because cluster validation is carried out prior to class assignment and class outputs are ordinal categorical values. The method is based on the principle that backscatter values from a single incidence angle express a normal distribution for a particular sediment type. The resulting Bayesian classes were well correlated to median grain sizes and the percentage of coarse material. The MLC method uses angular response information from five layers of training areas extracted from the Bayesian classification map. The subsequent PCA analysis is based on the transformation of these five layers into two principal components that comprise most of the data variability. These principal components were clustered in five classes after running an external cluster validation test. In general both methods MLC and PCA, separated the various sediment types effectively, showing good agreement (kappa >0.7) with the Bayesian approach which also correlates well with ground truth data (r2?>?0.7). In addition, sub-bottom data were used in conjunction with the Bayesian classification results to characterize acoustic classes with respect to their geological and stratigraphic interpretation. The joined interpretation of seafloor and sub-seafloor data sets proved to be an efficient approach for a better understanding of seafloor backscatter patchiness and to discriminate acoustically similar classes in different geological/bathymetric settings.  相似文献   

20.
《Marine Geology》2005,216(4):275-296
Recent chirp seismic reflection data combined with multibeam bathymetry, backscatter, and analysis of grab samples and short cores provide evidence of significant recent erosion on the outer New Jersey shelf. The timing of erosion is constrained by two factors: (1) truncation at the seafloor of what is interpreted to be the transgressive ravinement surface at the base of the surficial sand sheet, and (2) truncation of apparently moribund sand ridges along erosional swales oriented parallel to the primary direction of modern bottom flow and oblique to the strike of the sand ridges. These observations place the erosion in a marine setting, post-dating the passage of the shoreface ravinement and the evolution of sand ridges that form initially in the near shore environment. Also truncated by marine erosion are shallowly buried, fluvial channel systems, formed during the Last Glacial Maximum and filled during the transgression, and a regional reflector “R” that is > ∼ 40 kyr. Depths of erosion range from a few meters to > 10 m. The seafloor within eroded areas is often marked by “ribbon” morphology, seen primarily in the backscatter data as areas of alternating high and low backscatter elongated in the direction of primary bottom flow. Ribbons are more occasionally observed in the bathymetry; where observed, crests exhibit low backscatter and troughs exhibit high backscatter. Sampling reveals that the high backscatter areas of the ribbons consist of a trimodal admixture of mud, sand and shell hash, with a bimodal distribution of abraded and unabraded sand grains and microfauna. The shell hash is interpreted to be an erosional lag, while the muds and unabraded grains are, in this non-depositional environment, evidence of recent erosion at the seafloor of previously undisturbed strata. The lower-backscatter areas of the ribbon morphology were found to be a well-sorted medium sand unit only a few 10's of cm thick overlying the shelly/muddy/sandy material. Concentrations of well-rounded gravels and cobbles were also found in eroded areas with very high backscatter, and at least one of these appears to be derived from the base of an eroded fluvial channel. Seafloor reworking over the transgressive evolution of the shelf appears to have switched from sand ridge evolution, which is documented to ∼ 40 m water depth, to more strictly erosional modification at greater water depths. We suggest that this change may be related to the reduction with water depth in the effectiveness of sediment resuspension by waves. Resuspension is a critical factor in the grain size sorting during transport by bottom currents over large bedforms like sand ridges. Otherwise, we speculate, displacement of sand by unidirectional currents will erode the seafloor.  相似文献   

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

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