首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 250 毫秒
1.
多波束声呐记录的海底后向散射片段(Snippet)数据处理成角度响应曲线和地理编码(Mosaic)图像可以 帮助识别海底底质类型和反映地貌形态,这一过程包括辐射校正、角度响应改正(AVG)和几何地理编码,但不同的多波束系统硬件在辐射校正和角度响应改正方法上存在差异且传统处理方法忽略了声呐系统本身的指向性模型随时间变化的事实。以声呐方程为基础,针对Kongsberg EM 多波束系统提出了一套完整的Snippet数据处理流程,并分析了各步骤中存在的可变性,给出了每一步的处理建议,最后将此方法应用于EM2040浅水多波束实测数据,并验证了该方法的有效性和实用性。  相似文献   

2.
The relationship between acoustic backscatter, sediment characteristics and benthic habitat is examined using high-resolution sidescan sonar data collected at the Loch Linnhe artificial reef site on the west coast of Scotland. The site is typical for the continental shelf of NW Europe, with a mix of seabed environments from muddy to coarse, stony substrata on a 10–100 m length scale. A sidescan sonar mosaic was produced and classified according to derived backscatter parameters (mean, median and standard deviation of the backscatter values) using an unsupervised classification procedure. The accuracy of the final classified map was assessed by comparison with a ground-truthing survey in which the biological habitat was derived from underwater video footage. The sidescan correctly predicted seabed surface characteristics of observed biological habitat with 78% accuracy. A second, and more challenging test of the acoustic data to correctly predict biological habitat was made by comparing it with data from 21 grab sampling stations. These stations were divided into three groups using multivariate statistical techniques based on their backscatter properties. Benthic assemblage structure was found to be significantly distinct between the high and low, and the medium and low backscatter stations. There was a low to moderate but significant correlation between the multivariate patterns of acoustic backscatter, benthic assemblage structure, and particle size distribution. The work shows that even in areas with subtle and gradational changes in substratum, the sidescan was able to predict biological community with an acceptable accuracy.  相似文献   

3.
徐韦  程和琴  黄知  郑树伟  陈钢 《海洋学报》2019,41(1):172-182
近海海底地形探测与沉积物精确分类对涉海工程建设、生物栖息地反演以及海底资源勘查与开发具有重要的现实意义。以澳洲Joseph Bonaparte湾为例,利用多波束测深技术获取了该海湾约880 km2水域的水深数据与背向散射强度信号,结合同步采集的54个海底表层沉积物样品,通过随机决策树模型对该海域海底表层沉积物进行了分类研究。结果表明:(1)利用随机决策树模型分析该海域沉积物类型与背向散射强度的关系时,当模型内部参数设置:树的总数为200,最小分裂节点为2,每棵树的最大分裂级数为5时,可提高预测准确率;(2)该参数设置下,利用13°和37°入射角的背向散射强度预测该海域沉积物类型时,准确率最高,其值为83.3%,且在研究海域,砂质砾和砾质砂分布在背向散射强度较强的深槽或海沟等地区,而砾质泥质砂和含砾泥质砂主要分布在背向散射强度较弱的浅水海域。分析还发现,当水深数据作为预测海底表层沉积物类型的特征变量时,有可能降低最终预测结果的准确率。  相似文献   

4.
The sediment backscatter strength measured by multibeam echosounders is a key feature for seafloor mapping either qualitative (image mosaics) or quantitative (extraction of classifying features). An important phenomenon, often underestimated, is the dependence of the backscatter level on the azimuth angle imposed by the survey line directions: strong level differences at varying azimuth can be observed in case of organized roughness of the seabed, usually caused by tide currents over sandy sediments. This paper presents a number of experimental results obtained from shallow-water cruises using a 300-kHz multibeam echosounder and specially dedicated to the study of this azimuthal effect, with a specific configuration of the survey strategy involving a systematic coverage of reference areas following “compass rose” patterns. The results show for some areas a very strong dependence of the backscatter level, up to about 10-dB differences at intermediate oblique angles, although the presence of these ripples cannot be observed directly—neither from the bathymetry data nor from the sonar image, due to the insufficient resolution capability of the sonar. An elementary modeling of backscattering from rippled interfaces explains and comforts these observations. The consequences of this backscatter dependence upon survey azimuth on the current strategies of backscatter data acquisition and exploitation are discussed.  相似文献   

5.
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.  相似文献   

6.
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.  相似文献   

7.
Cold seeps on the Hikurangi Margin off New Zealand exhibit various seabed morphologies producing different intensity patterns in sidescan backscatter images. Acoustic backscatter characteristics of 25 investigated seep sites fall into four distinct types characterised by variations in backscatter intensity, distribution and inferred structural heights. The types reflect different carbonate morphologies including up to 20-m-high structures (type 1), low-relief crusts (type 2), scattered blocks (type 3) and carbonate-free sites (type 4). Each seep corresponds to a single type; intermediates were not observed. This correlates well with published data on seep fauna at each site, with the four types representing three different faunal habitats of successive stages of seep development. Backscatter signatures in sidescan sonar images of cold seeps may therefore serve as a convenient proxy for variations in faunal habitats.  相似文献   

8.
济州岛南部海域海底声呐图像分析与声学底质分类   总被引:2,自引:2,他引:0  
东海北部外陆架靠近济州岛南部海域,是黄海槽向冲绳海槽延伸的部分,属于黑潮分支黄海暖流的通道入口,分布着脊槽相间的海底底形,对其海底声呐图像的处理分析及声学底质分类的分析研究,有助于了解该通道海底底形表层纹理特征及沉积物分布规律。基于在济州岛南部海域获取的多波束声呐数据,应用图像处理技术和方法,对数据进行了处理,获得了海底声呐影像图,并对其表层纹理特征进行了描述和分析;同时,基于多波束反向散射强度数据,结合19组海底地质取样数据,建立研究区海底反向散射强度与沉积物粒度特征之间的统计关系模型,并以改进的学习向量量化神经网络方法,实现对海底粉砂质砂、黏土质砂以及砂-粉砂-黏土3种底质类型的快速自动分类识别。  相似文献   

9.
Acoustic seabed classification is a useful tool for monitoring marine benthic habitats over broad-scales (>1 km2) and meso-scales (10 m2–1 km2). Its utility in this context was evaluated using two approaches: by describing natural changes in the temporal distribution of marine biotopes across the broad-scale (4 km2), and by attempting to detect specific experimentally-induced changes to kelp-dominated biotopes across the meso-scale (100 m2). For the first approach, acoustic backscatter mosaics were constructed using sidescan sonar and multibeam echosounder data collected from Church Bay (Rathlin Island, Northern Ireland) in 1999, 2008 and 2009. The mosaics were manually segmented into acoustic facies, which were ground-truthed using a drop-video camera. Biotopes were classified from the video by multivariate exploratory analysis and cross-tabulated with the acoustic facies, showing a positive correlation. These results were integrated with bathymetric data to map the distribution of seven unique biotopes in Church Bay. Kappa analysis showed the biotope distribution was highly similar between the biotope maps, possibly due to the stability of bedforms shaped by the tidal regime around Rathlin Island. The greatest biotope change in this approach was represented by seasonal and annual changes in the growth of the seagrass, Zostera marina. In the second approach, sidescan sonar data were collected before and after the removal of 100 m2 of kelp from three sites. Comparison of the data revealed no differences between the high-resolution backscatter imagery. It is concluded that acoustic seabed classification can be used to monitor change over broad- and meso-scales but not necessarily for all biotopes; its success depends on the type of acoustic system employed and the biological characteristics of the target biotope.  相似文献   

10.
The increased use of backscatter measurements in time series for environmental monitoring necessitates the comparability of individual results. With the current lack of pre-calibrated multibeam echosounder systems for absolute backscatter measurement, a pragmatic solution is the use of natural reference areas for ensuring regular assessment of the backscatter measurement repeatability. This method mainly relies on the assumption of a sufficiently stable reference area regarding its backscatter signature. The aptitude of a natural area to provide a stable and uniform backscatter response must be carefully considered and demonstrated by a sufficiently long time-series of measurements. Furthermore, this approach requires a strict control of the acquisition and processing parameters. If all these conditions are met, stability check and relative calibration of a system are possible by comparison with the averaged backscatter values for the area. Based on a common multibeam echosounder and sampling campaign completed by available bathymetric and backscatter time series, the suitability as a backscatter reference area of three different candidates was evaluated. Two among them, Carré Renard and Kwinte, prove to be excellent choices, while the third one, Western Solent, lacks sufficient data over time, but remains a valuable candidate. The case studies and the available backscatter data on these areas prove the applicability of this method. The expansion of the number of commonly used reference areas and the growth of the number of multibeam echosounder controlled thereon could greatly contribute to the further development of quantitative applications based on multibeam echosounder backscatter measurements.  相似文献   

11.
Two highly active mud volcanoes located in 990–1,265 m water depths were mapped on the northern Egyptian continental slope during the BIONIL expedition of R/V Meteor in October 2006. High-resolution swath bathymetry and backscatter imagery were acquired with an autonomous underwater vehicle (AUV)-mounted multibeam echosounder, operating at a frequency of 200 kHz. Data allowed for the construction of ~1 m pixel bathymetry and backscatter maps. The newly produced maps provide details of the seabed morphology and texture, and insights into the formation of the two mud volcanoes. They also contain key indicators on the distribution of seepage and its tectonic control. The acquisition of high-resolution seafloor bathymetry and acoustic imagery maps with an AUV-mounted multibeam echosounder fills the gap in spatial scale between conventional multibeam data collected from a surface vessel and in situ video observations made from a manned submersible or a remotely operating vehicle.  相似文献   

12.
Reefs are one of the marine habitats listed in Annex I of the European Union's Habitats Directive, which aims to establish a coherent European ecological network of Special Areas of Conservation. EU Member States are required to prepare and propose a national list of sites for evaluation under the scheme, but currently the occurrence of reefs in the United Kingdom's nearshore and offshore areas is not well documented. Here we report on our search for rocky reefs in the central English Channel, which unexpectedly revealed an extensive reef system covering an area of 1100 km2. Prior to our work, it was generally perceived that the seabed in this area comprised mostly gravel, with a few isolated rock outcrops.Our approach to determining the location, extent and character of these reefs incorporated broad, medium and fine-scale analyses over a 3200 km2 area of seabed, using single- and multi-beam acoustic data, ground-truthed by underwater video and stills imagery. A benthic terrain model was developed in ArcGIS to map topographic features at the broad and medium scales. Biotope assignments were made at the fine scale through detailed analysis of video footage obtained from 30 sampling stations. The study area has a complex geological history and lies at the centre of a major bedload-parting zone. Together, these strongly influence the seabed character and the distribution of biotopes. An integrated assessment of the physical and biological features was used to map the study area to level 4 of the EUNIS habitat classification system.Similar physical conditions exist in other areas of the UK continental shelf, raising the prospect of predicting where other rocky reef systems might occur. In the absence of a co-ordinated national seabed survey programme, such predictions, coupled with interpretation of existing single-beam bathymetry data, can help prioritise areas where limited survey resources could be most effectively deployed.  相似文献   

13.
Side-scan sonar, camera observations, and sediment sampling are used to determine the seabed dispersal of mussel debris at three aquaculture farm sites in sheltered to exposed environments in the Marlborough Sounds, New Zealand. At the two sheltered farm sites, side-scan sonagraphs reveal high acoustic backscatter beneath culture sites and low backscatter in the immediate surrounds. This suggests a sharp contrast in seabed material around the edges of the farm. Photography and sediment samples reveal farms are underlain by mounds of shells with biodeposits infilling intershell voids and forming a veneer over entire mounds. In contrast, the surrounding seabed is naturally sedimented soft mud. Sediment from beneath the farms had total organic contents of 8%-19% decreasing sharply to natural levels of 4%-7%, 30 m from the farm's boundaries. The third site is exposed to energetic waves and currents, and exhibits low acoustic backscatter zones beneath the farm. Photographs show the seabed is predominantly clean sand with only scattered clumps of shells and no biodeposits visible. Sediments are typical of the natural setting. This lack of mussel debris indicates that there is sufficient energy to transport and dissipate the shell and biodeposits over a wide area with little impact on the natural sediment.  相似文献   

14.
《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.  相似文献   

15.
Abstract

Concentrations of cadmium (Cd), lead (Pb), copper (Cu), and zinc (Zn) were measured in adult barnacles (Elminius modestus Darwin) from Waitemata and Manukau Harbours in the Auckland area, New Zealand. As in studies on sediments reported in the literature, it was possible to identify areas of likely anthropogenic influence, e.g., around the Auckland Harbour Bridge for Pb, Cu, and Zn. Groups of individuals with highest concentrations for these metals showed 19.8–23.8 mg Pb kg?1, 198–266 mg Cu kg?1, and 4460–6530 mg Zn kg?1 (95% confidence limits, dry weight basis). Cd concentrations found for all barnacles from the Auckland area ranged from 0.8 to 3.1 mg kg?1. Two samples from Omaha Beach, 60 km north of Auckland, were used as a reference. Accordingly, groups of individuals with lowest concentrations for Pb, Cu, and Zn could be allocated to this site using the Student‐Newman‐Keuls Multiple Range Test (0.5–1.3 mg Pb kg?1, 8–10 mg Cu kg?1, and 144–214 mg Zn kg?1 ; 95% confidence limits). Only Cd concentrations were highest at Omaha Beach (8.6–12.1 mg Cd kg?1 ). This result may have arisen from “naturally” increased bio‐availabilities of certain metals in mangrove systems which are reported in the literature. Generally, metal concentrations in barnacles from the Auckland Harbour area and from Omaha Beach were within the wider range for E. modestus as well as other barnacle species reported in the international literature.  相似文献   

16.
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.  相似文献   

17.
Dual-frequency side-looking sonars have the potential to be used as remote sensing tools to characterize subaqueous terrains. In one case study of the carbonate-ooze-coated Blake Plateau off-shore of Georgia, U.S.A., the difference in acoustic attenuation for 50 and 20 mm wavelengths (30 and 72 kHz frequency) permits the discrimination of sub-bottom scatterers from seabed surface textural features to reveal patchy regions where a buried hard ground had been pock-marked by karst-like depressions. In a second study of the Upper Hudson River in New York, U.S.A., related to environmental contaminates, the backscatter response at 15 and 3 mm acoustic wavelengths (100 and 500 kHz frequency) serves as a useful proxy for sediment grain size with coarser detritus distinguished from finer sediments. Sand and gravel regions inferred from the backscatter were confirmed by ground truth sampling.  相似文献   

18.
台湾岛西南海域福尔摩沙海脊冷泉区地形地貌特征分析   总被引:2,自引:0,他引:2  
利用水深数据和ROV近海底影像资料,对福尔摩沙海脊冷泉区的海底地形地貌和冷泉系统的海底表征进行了描述和分析,并讨论了二者之间的响应关系。结果表明,相对于船载多波束数据而言,近海底多波束测深系统所获得的数据能更高精度地反映冷泉区海底地形地貌特征,是研究冷泉系统不可或缺的基础资料。基于ROV近海底观测影像资料,福尔摩沙海脊冷泉系统整体表现为局部被化能自养生物群落覆盖并有流体喷口零星分布的巨大自生碳酸盐岩岩丘,海底表征主要包括形态各异的自生碳酸盐岩结壳或岩体、化能自养生物群落、流体喷口、还原性沉积物等几种形式。研究表明,福尔摩沙海脊冷泉区的地形地貌特征与冷泉系统海底表征具有良好的响应关系,并且该区的地形地貌特征主要受控于出露于海底的自生碳酸盐岩的形态特征及规模。首次揭示了福尔摩沙海脊冷泉区地形地貌特征与其海底表征之间的响应关系,以期为后续的冷泉研究提供必要的背景资料支持。  相似文献   

19.
The 1982 United Nations Convention on the Law of the Sea declares the seabed beyond national jurisdiction and its mineral resources as the “common heritage of mankind” (CHM). This article examines the operationalisation of the CHM principle in the international seabed mining regime, with focus placed on the sharing of benefits derived from mining. The article begins by providing an overview of the CHM principle, before examining four modalities provided for in the Convention, both institutional and substantial, and their role in giving effect to the CHM principle: (1) financial benefits; (2) the “Enterprise”; (3) the parallel system of reserved areas; and, (4) marine scientific research. Finally, overarching issues are discussed and some suggestions on ways forward are presented. The article considers that the deep seabed mining regime is not yet ready to effectively share the benefits derived from the common heritage of mankind. In particular, the future of the Enterprise is uncertain and changes to the so-called parallel system that affect the CHM have received minimal discussion. Moreover, a lack of publicly available research data related to seabed mining is hindering current benefits for humankind. However, work is underway at the International Seabed Authority to establish rules and policies with respect to the sharing of financial benefits. While the CHM principle remains largely untested, approaches that are transparent, inclusive, accountable, and equitable are more likely to be successful.  相似文献   

20.
The purpose of this study was to investigate, characterize and map areas of the seabed of Las Perlas Archipelago (LPA), Republic of Panama using swath-bathymetry acoustic sonar techniques and to assess these methods as tools for feeding information into management zoning policy. In 2007 the LPA was granted conservation protection under national legislation. However, detailed management plans are still pending. Seabed mapping plays a fundamental role in identifying areas which should be prioritized within the management framework. Visual representation of habitat maps provides an effective medium for involving stakeholders in a co-management arena. In this survey, acoustically mapped areas of the seabed were ground-truthed using a combination of benthic grab samples, drop-down video and diver observations. The resulting mapped areas were then incorporated into a Geographic Information System (GIS) for further analysis. The output was a physical characterization of the seabed at three locations selected for being areas of high rugosity (habitat complexity) and, therefore, their potential importance as valuable fish aggregation sites. The rocky reefs and rhodolith beds identified in this survey represent particularly important fish aggregation and nursery habitats which should be considered priorities for protection under the management plans. This survey demonstrated the use of acoustic techniques to spatially resolve topographic features and physical characteristics of the seabed, illustrating their potential value as tools for fisheries management and marine reserve zoning in Las Perlas Archipelago and elsewhere.  相似文献   

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

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