首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
High energy dynamic compaction (HEDC) is adopted in a coastal reclamation area because the grain size of backfilled soil mostly ranges between 20 cm and 100 cm. The in situ tests for evaluating the effectiveness of HEDC were performed on the backfilled soil ground. The crater depth per drop and the whole test zone elevations before and after HEDC were measured and analyzed. Dynamic penetration tests and spectral analysis of surface wave (SASW) tests were used for investigating the improvement depth. Furthermore, the allowable bearing capacity of HEDC treated ground was determined based on the results of plate-load tests. It was found that HEDC did not cause the ground surface heave during construction, and was more effective than low energy dynamic compaction (LEDC) in terms of applied energy utilization. Based on the test results, the improvement depth of HEDC at this site was not less than 14 m, and there was no obvious weak layer within the range of improvement depth. The allowable bearing capacities were larger than 160 kPa. The investigation results indicate that the HEDC technique is an effective way for improving backfilled coarse-grained soil in coastal reclamation areas. This technique helps to achieve both greater improvement depths and higher ground bearing capacities as compared with LEDC.  相似文献   

2.
A regional ocean circulation model with four-dimensional variational data assimilation scheme is configured to study the ocean state of the Indian Ocean region (65°E–95°E; 5°N–20°N) covering the Arabian Sea (AS) and Bay of Bengal (BoB). The state estimation setup uses 10 km horizontal resolution and 5 m vertical resolution in the upper ocean. The in-situ temperature and salinity, satellite-derived observations of sea surface height, and blended (in-situ and satellite-derived) observations of sea surface temperature alongwith their associated uncertainties are used for data assimilation with the regionally configured ocean model. The ocean state estimation is carried out for 61 days (1 June to 31 July 2013). The assimilated fields are closer to observations compared to other global state estimates. The mixed layer depth (MLD) of the region shows deepening during the period of assimilation with AS showing higher MLD compared to the BoB. An empirical forecast equation is derived for the prediction of MLD using the air–sea forcing variables as predictors. The surface and sub-surface (50 m) heat and salt budget tendencies of the region are also investigated. It is found that at the sub-surface, only the advection and diffusion temperature and salt tendencies are important.  相似文献   

3.
Several cold vents are observed at the northern Cascadia margin offshore Vancouver Island in a 10 km2 region around Integrated Ocean Drilling Program Expedition 311 Site U1328. All vents are linked to fault systems that provide pathways for upward migrating fluids and at three vents methane plumes were detected acoustically in the water column. Downhole temperature measurements at Site U1328 revealed a geothermal gradient of 0.056 ± 0.004°C/m. With the measured in situ pore-water salinities the base of methane hydrate stability is predicted at 218–245 meters below seafloor. Heat-probe measurements conducted across Site U1328 and other nearby vents showed an average thermal gradient of 0.054 ± 0.004°C/m. Assuming that the bottom-simulating reflector (BSR) marks the base of the gas hydrate stability zone variations in BSR depths were used to investigate the linkages between the base of the gas hydrate stability zone and fluid migration. Variations in BSR depth can be attributed to lithology-related velocity changes or variations of in situ pore-fluid compositions. Prominent BSR depressions and reduced heat flow are seen below topographic highs, but only a portion of the heat flow reduction can be due to topography-linked cooling. More than half of the reduction may be due to thrust faulting or to pore-water freshening. Distinct changes in BSR depth below seafloor are observed at all cold vents studied and some portion of the observed decrease in the BSR depth was attributed to fault-related upwelling of warmer fluids. The observed decrease in BSR depth below seafloor underneath the vents ranges between 7 and 24 m (equivalent to temperature shifts of 0.07–0.15°C).  相似文献   

4.
2D and 3D seismic reflection and well log data from Andaman deep water basin are analyzed to investigate geophysical evidence related to gas hydrate accumulation and saturation. Analysis of seismic data reveals the presence of a bottom simulating reflector (BSR) in the area showing all the characteristics of a classical BSR associated with gas hydrate accumulation. Double BSRs are also observed on some seismic sections of area (Area B) that suggest substantial changes in pressure–temperature (P–T) conditions in the past. The manifestation of changes in P–T conditions can also be marked by the varying gas hydrate stability zone thickness (200–650 m) in the area. The 3D seismic data of Area B located in the ponded fill, west of Alcock Rise has been pre-stack depth migrated. A significant velocity inversion across the BSR (1,950–1,650 m/s) has been observed on the velocity model obtained from pre-stack depth migration. The areas with low velocity of the order of 1,450 m/s below the BSR and high amplitudes indicate presence of dissociated or free gas beneath the hydrate layer. The amplitude variation with offset analysis of BSR depicts increase in amplitude with offset, a similar trend as observed for the BSR associated with the gas hydrate accumulations. The presence of gas hydrate shown by logging results from a drilled well for hydrocarbon exploration in Area B, where gas hydrate deposit was predicted from seismic evidence, validate our findings. The base of the hydrate layer derived from the resistivity and acoustic transit-time logs is in agreement with the depth of hydrate layer interpreted from the pre-stack depth migrated seismic section. The resistivity and acoustic transit-time logs indicate 30-m-thick hydrate layer at the depth interval of 1,865–1,895 m with 30 % hydrate saturation. The total hydrate bound gas in Area B is estimated to be 1.8 × 1010 m3, which is comparable (by volume) to the reserves in major conventional gas fields.  相似文献   

5.
Remote sensing bathymetry inversion can quickly obtain water depth data of large areas, but this process relies on a large number of in-situ depth data points. USV-based (Unmanned Surface Vehicle) technique can obtain the bathymetry data of shallow water where ordinary ships are inaccessible, but this technique is inefficient and generally only data along survey line can be collected. The combination of USV and high-resolution remote sensing provides a new solution for water depth surveying and mapping around an island. This paper focuses on the key techniques, using USV sounding data and GeoEye-1 multispectral remote sensing images covering the region of Wuzhizhou island in the experiment. The results show that the MAE (Mean Absolute Error) of USV sounding is 0.25 m, while the MRE (Mean Relative Error) is 1.41%, and the MRE of remote sensing bathymetry aided by USV sounding can be controlled within 20%. Errors are mainly from areas shallower than 5 m, and are also affected by the USV sounding position accuracy. It shows that it is feasible to combine the USV sounding and high-resolution remote sensing bathymetry, and this technique has broad application prospects in the field of bathymetry in large shallow areas.  相似文献   

6.
This article presents a case history of determination of effective depth of prefabricated vertical drains (PVDs) under embankment loading on a very soft clay deposit in central China, near Jiujiang, Jiangxi Province. The height of the embankment was 5.3 m and construction time was about one year. The PVDs were installed to a depth of 8.5 m at a spacing of 1.5 m in a triangular pattern. Field observations and the finite element method (FEM) were employed to analyze the performance of the soft deposit during embankment construction. The influential depth of the embankment loading was evaluated based on settlement, excess pore pressure, and stress increase in subsoil, both from the observed data and FEM analysis. The effective PVD depth was determined in the following ways: (1) the depth of 5% subsoil settlement of surface settlement; (2) vertical stress increase in subsoil of 25% in-situ stress; and (3) consolidation time/PVD depth relation by FEM. Based on the analysis, the effective depth of PVDs was determined to be between 10 and 12.8 m for this field case.  相似文献   

7.
Coral reef calcareous sediment, a special category of rock-soil material, has representative geological structure and environmental characteristics. It is widely distributed in shallow areas of the tropical ocean; therefore the exploration instruments and technologies for engineering geology studies of coral reef calcareous sediments are very different from those used in land or deep-sea. Obtaining undisturbed cores from the Holocene unconsolidated stratum of coral reefs has been a key problem in the field of marine geology and environment surveying. The authors have designed a novel floating drilling platform equipped with a drilling machine, and successfully achieved undisturbed cores from both reef flats (with water depth 0.5 to 2 m) and a enclosed lagoon (with water depth 2–12 m) of Yongshu Reef (9°32–9°42 N, 112°52–113°04 E), southern South China Sea. Based on the detailed observation on the cores and the analysis from engineering geology, Yongshu Reef was split into reef three engineering geological zones: leading edge, reef flat (including outer reef flat, middle reef flat and inner reef flat) and lagoon. The sediments are classified in the stratum as fine sand, medium sand, coarse sand, gravel and weak-cemented reef limestone.  相似文献   

8.
A spectacular hummocky topography was discovered offshore of the south-eastern slope of the Nisyros island volcano in the eastern sector of the Aegean volcanic arc in 2000–2001, using multibeam bathymetric mapping and seismic profiling, and interpreted as part of a volcanic debris avalanche originating onland. During E/V Nautilus cruise NA011 in 2010, a detailed side-scan sonar and ROV exploration aimed at evaluating the surface morphology of this avalanche field. Combining the new data with selected older datasets reveals that the debris avalanche is characterized by numerous (at least 78) variously sized and shaped hummocks. Some of these are distinctly round, either scattered or aligned in groups, whereas others are elongated in the form of ridges. This is consistent with existing models accounting for variations in the longitudinal and lateral velocity ratio of landslides. Maximum dimensions reach 60 m in height above the sea bottom, 220 m in length and 230 m in width. The structures outline a large tongue-shaped, submarine hummock field of about 22.2 km2, approx. 4.8 km wide and 4.6 km long and with an estimated volume of 0.277 km3. Due to its characteristic shape, the collapsed volcanic flank is interpreted to represent a singular failing event, involving a rapid and virtually instantaneous downslope movement of the slide mass into the sea. Indeed, the H/L (height of 280 m vs. run-out of 7 km) ratio for the Nisyros slide is 0.04; plotted against volume, this falls within the theoretical bounds as well as measured values typical of submarine landslides. The timing of the event is probably related to the extrusion of Nikia lavas and their subsequent failure and formation of a main scarp observed at about 120 m depth on an 8-km-long seismic profile and a map of slope angle distribution, at the depth where the palaeo-coastline was located 40 ka ago. An inferred age of ca. 40 ka for the avalanche awaits confirmation based on dating of core material.  相似文献   

9.
A central theme of the ongoing GEOTRACES program is to improve the understanding of processes occurring at ocean interfaces with continents, sediments, and ocean crust. In this context, we studied the distributions of Al, Mn, Fe, Co, Ni, Cu, Zn, Cd, and Pb around the Juan de Fuca Ridge (JdFR) in total dissolvable (td), dissolved (d), and labile particulate (lp) fractions, which represent a fraction in unfiltered samples, filtered samples through an AcroPak capsule filter, and the difference between td and d, respectively. Al and Fe were dominated by lp-species, while Ni, Zn, and Cd were dominated by d-species with undetectable amounts of lp-species. Major findings in this study are as follows: (1) The continental margin (CM) provided large sources of Al, Mn, Fe, and Co from the surface to ~2000 m in depth. The supply from CM caused high surface concentrations of dMn and dCo, a subsurface (100–300 m depth) maximum of dCo, and intermediate (500–2000 m depth) maxima of lpAl and lpFe. The supply of dFe from CM was ~10 times that from the high-temperature hydrothermal activity at station BD21, which is located at ~3 km from the Middle Valley venting site and ~ 200 km from Vancouver Island. (2) DPb was maximum at the top layer of North Pacific Intermediate Water, probably owing to isopycnal transport of anthropogenic Pb via advection of subducted surface waters. Although dCo and dPb had different sources in the upper water, they showed a strong linearity below 300 m (r 2 = 0.95, n = 38), indicating concurrent scavenging. (3) A high-temperature hydrothermal plume occurred at a depth of 2300 m at BD21, accounting for maxima of dAl, dMn, dFe, lpCu, and lpPb and a minimum of dCu. (4) Strong bottom maxima of lpAl, lpMn, lpFe, lpCo, and lpPb occurred above the abyssal plain at the western foot of the JdFR, indicating resuspension of sediments. However, bottom maxima of d-species were apparent only for dAl and dCu.  相似文献   

10.
Sea floor morphology plays an important role in many scientific disciplines such as ecology, hydrology and sedimentology since geomorphic features can act as physical controls for e.g. species distribution, oceanographically flow-path estimations or sedimentation processes. In this study, we provide a terrain analysis of the Weddell Sea based on the 500 m × 500 m resolution bathymetry data provided by the mapping project IBCSO. Seventeen seabed classes are recognized at the sea floor based on a fine and broad scale Benthic Positioning Index calculation highlighting the diversity of the glacially carved shelf. Beside the morphology, slope, aspect, terrain rugosity and hillshade were calculated and supplied to the data archive PANGAEA. Applying zonal statistics to the geomorphic features identified unambiguously the shelf edge of the Weddell Sea with a width of 45–70 km and a mean depth of about 1200 m ranging from 270 m to 4300 m. A complex morphology of troughs, flat ridges, pinnacles, steep slopes, seamounts, outcrops, and narrow ridges, structures with approx. 5–7 km width, build an approx. 40–70 km long swath along the shelf edge. The study shows where scarps and depressions control the connection between shelf and abyssal and where high and low declination within the scarps e.g. occur. For evaluation purpose, 428 grain size samples were added to the seabed class map. The mean values of mud, sand and gravel of those samples falling into a single seabed class was calculated, respectively, and assigned to a sediment texture class according to a common sediment classification scheme.  相似文献   

11.
Ocean General Circulation Model (OGCM) simulations from 1970–2007 are used to study the upper ocean heat content variability in the Tropical Indian Ocean (TIO). Model computed heat contents up to 50 m (denoted by HC50 m hereafter) representing upper ocean heat content and 300 m (HC300 m) representing heat content up to thermocline depth are first compared with heat contents computed from observations of two buoys in the TIO. It is found that there is good agreement between the model and observations. Fourier analysis of heat content is carried out in different regions of TIO. The amplitudes of semi-annual variability for HC50 m and HC300 m are observed to be greater than those for the annual variability in the Bay of Bengal, while in the Arabian Sea there is a mixed result. Heat content tendency is known to be governed by net surface heat flux and horizontal as well as vertical heat transports. For understanding the relative importance of these processes, a detailed analysis of these terms in the tendency equation is carried out. Rossby wave is observed in the annual mode of heat transport while equatorial jet and Kelvin waves are observed in the semi-annual mode of heart transport. Finally, the correlation between heat content and sea surface temperature (SST) and sea level anomaly (SLA), taken one at a time, is computed. It is found that the correlation improves significantly when both these quantities are together taken into account.  相似文献   

12.
To better understand the vertical distribution of phytoplankton in the tropical and subtropical North Pacific, we used fast repetition rate fluorometry to investigate the photo-physiological condition of the phytoplankton assemblage in this region between February and March 2007. Along 155°E, between the equator and 24°N, the peak of fluorescence (F m), an indication of the deep chlorophyll maximum (DCM), was deeper than the top of the nitracline and occurred at the 2.4 ± 1.3 % (mean ± SD) light depth (relative to 0 m). The photochemical efficiency (F v/F m) and effective absorption cross-section of photosystem II (σPSII) were low at the surface but increased rapidly at depths between the top of the nitracline (40–138 m) and the DCM (70–158 m), an indication that the photo-physiological condition of the phytoplankton improved below the top of the nitracline. The depth of the maximal F v/F m [Z(F v/F m max)] was 18–32 m deeper than the DCM and corresponded to the 0.8 ± 0.2 % light depth. The values of F v/F m at the Z(F v/F m max) were 20 % higher than those at the DCM and averaged 0.48 ± 0.01. These results suggest that the phytoplankton assemblage beneath the DCM had a high potential photosynthetic performance capacity and was growing by using the very low ambient light in this region.  相似文献   

13.
We have developed new systems capable of profiling to >1000 m for measuring in situ pH and fugacity of CO2 (fCO2) in the ocean using spectrophotometric analysis (pH and CO2 profilers). The in situ pH is determined by detecting the color change of the pH indicator (m-cresol purple). It can withstand ambient pressure to 1000 m depth. The CO2 profiler analyzed in situ fCO2 by detecting the change of pH in an inner solution, equilibrated with the seawater through a gas permeable membrane. It can be operated to 2500 m depth. We used an amorphous fluoropolymer tubing form of AF-2400 for the gas permeable membrane due to its high gas permeability coefficients. The inner solution was a mixture of 2 μM bromocresol purple (BCP) and 5 μM sodium hydroxide. This system gave us a response time of 1 minute, which is twice as fast as previous systems. The precisions of pH and CO2 profilers were within 0.002 and 2.5% respectively. We have used these profilers to study the North Pacific, obtaining good agreement with the difference between the data from profilers and a discrete bottle of 0.002 ± 0.005 pH (SE, n = 25) and −0.4 ± 3 μatm (SE, n = 31).  相似文献   

14.
Abstract

This article examines whether Digital Elevation Model (DEM) resolution affects the accuracy of predicted coastal inundation extent using LISFLOOD-FP, with application to a sandy coastline in New Jersey. DEMs with resolution ranging from 10 to 100 m were created using coastal elevation data from NOAA, using the North American Vertical Datum of 1988. A two-dimensional hydrodynamic flood model was developed in LISFLOOD-FP using each DEM, all of which were calibrated and validated against an observed 24-h tidal cycle and used to simulate a 1.5 m storm surge. While differences in predicted inundated area from all models were within 1.0%, model performance and computational time worsened and decreased with coarser DEM resolution, respectively. This implied that using a structured grid model for modeling coastal flood vulnerability is based on two trade-offs: high DEM resolution coupled with computational intensity, but higher precision in model predictions, and vice versa. Furthermore, water depth predictions from all DEMs were consistent. Using an integrated numerical modeling and GIS approach, a two-scale modeling strategy, where a coarse DEM is used to predict water levels for projection onto a fine DEM was found to be an effective, and computationally efficient approach for obtaining reliable estimates of coastal inundation extent.  相似文献   

15.
Regulations are being developed through the International Seabed Authority (ISBA) for the exploration and mining of cobalt-rich ferromanganese crusts. This paper lays out geologic and geomorphologic criteria that can be used to determine the size and number of exploration and mine-site blocks that will be the focus of much discussion within the ISBA Council deliberations. The surface areas of 155 volcanic edifices in the central equatorial Pacific were measured and used to develop a mine-site model. The mine-site model considers areas above 2,500 m water depth as permissive, and narrows the general area available for exploration and mining to 20% of that permissive area. It is calculated that about eighteen 100 km2 exploration blocks, each composed of five 20 km2 contiguous sub-blocks, would be adequate to identify a 260 km2 20-year-mine site; the mine site would be composed of thirteen of the 20 km2 sub-blocks. In this hypothetical example, the 260 km2 mine site would be spread over four volcanic edifices and comprise 3.7% of the permissive area of the four edifices and 0.01% of the total area of those four edifices. The eighteen 100 km2 exploration blocks would be selected from a limited geographic area. That confinement area is defined as having a long dimension of not more than 1,000 km and an area of not more than 300,000 km2.  相似文献   

16.
Changes in the shore topography (e.g. slope) occur at a scale of hundreds of meters in several locations in the Lusitanian and the Mediterranean Sea provinces. We tested whether differences in the bottom inclination might affect the vertical distribution patterns of two sympatric coastal labrid fishes, the rainbow wrasse Coris julis and the ornate wrasse Thalassoma pavo. Visual censuses were used to determine the distribution and abundance of these labrid species in high (≥30°) and low (≤3°) slope rocky substrates covered by brown macroalgae and at two different depths (shallow, 4–7 m, and deep 14–20 m). Pectoral fin aspect ratio was used as an estimate of swimming performance to potentially explain the patterns observed. Despite the intrinsic biogeographical differences in the overall density of T. pavo and C. julis, on steep coasts the ornate wrasse dominated in shallow waters, whereas the two species coexisted both in shallow and deeper depths on gentle slope coasts. These distribution patterns were consistent across locations, and fin aspect ratio was not a good predictor of between‐habitat use for wrasses. We show that, under specific topographical conditions, the depth segregation pattern seems to be an interactive segregation (likely related to resource competition) rather than a result of selective segregation due to morphological differences in the pectoral fin. Significant ecological changes might occur in locations where the density of T. pavo has recently increased as a result of water warming.  相似文献   

17.
Multibeam bathymetric surveys and single-beam profiles were collected in 2003–2010 from aboard the Ukrainian RV Professor Vodyanitskiy (cruises PV-58 and PV-60, 2003 and 2004), and the German RV Meteor (cruise M-72, legs 1 and 4, 2007) and RV Maria S. Merian (cruise MSM-15, leg 2, 2010) along the continental margin of the NW Black Sea. Integrating published, reprocessed and novel data has revealed the existence of a major continuous channel extending from the Dnepr paleo-delta into greater water depths. It is more than 90 km long, 1.1 km wide and up to 125 m deep. On the upper slope (120–960 m water depth), a number of smaller channels merge into the large, Y-shaped Dnepr Canyon, which then continues obliquely downslope via this submarine channel to at least 1,815 m water depth off the Crimean continental margin, NW Black Sea. The channel could be an important, hitherto unknown link between the shallow oxic and deep anoxic environments of the Black Sea, along which sediment and organic matter could be funneled into the deep-sea basin. This would have far-reaching implications for investigations dealing with marine geology and biology, climate change, as well as oil and natural gas exploitation. The unusual alignment of the channel along the margin of the basin, as well as the location and mode of channel termination in deeper waters deserve future research.  相似文献   

18.
This paper presents a case study of dynamic compaction (DC) on backfill ground for planned oil tanks located at Nanjing Bay. The ultra-high energy DC level of 18,000 kN · m was applied in the area tamping phase of the dynamic compaction. In combination with ground replacement technique by forcing crushed stones into the underlain cohesive layer, the high-compressibility of this layer was properly tackled. The construction techniques and ground improvement mechanisms are introduced and discussed. Field evaluations before and after the dynamic compaction, including surface wave test, plate load test, and soil sampling, showed that the ultra-high energy dynamic compaction in combination with ground replacement achieved much greater improvement depth and ground bearing capacity than conventional DC. Results show that dynamic compaction of 18,000 kN · m combined with ground replacement is of great merit in dealing with backfill ground with weak cohesive content in harbor development.  相似文献   

19.
In this work, the vertical structure and variability along the western boundary of the Philippines are investigated using direct observations from acoustic Doppler current profiler (ADCP), Doppler volume sampler (DVS) and Aanderaa Seaguard instruments, which are mounted on a subsurface mooring deployed at 8°N, 127°3′E. In climatology, the southward Mindanao Current (MC) and northward Mindanao Undercurrent (MUC) play a dominant role in the upper layer. The mean currents at 1200 and 3500 m flow northward, whereas those at 2500 and 5600 m flow equatorward. The power spectral density (PSD) shows that an intraseasonal signal of 60–80 days is common from the sea surface to the bottom. The semiannual signals are strongest in the MUC layer, and the amplitude then decreases with depth to 3500 m. The seasonal variability at 2500 and 5600 m is similar between the two depths, suggesting a southward current in winter and northward flow in autumn. The current at 3500 m exhibits a northward direction in spring and southward flow in winter. In addition, the linear correlations between mooring data and altimetry products indicate that the variations in surface meridional currents along the western boundary of the Pacific Ocean can reach the bottom via low-frequency processes. The vertical-mode decomposition for observations indicates that the first four modes can effectively capture the original data. The relative contributions of different modes exhibit seasonal variability. The first baroclinic mode plays a dominant role in spring and autumn. In winter and summer, its contribution decreases and becomes comparable to that of the other modes.  相似文献   

20.
In this study, CPT and FVT tests are carried out on young, normally consolidated clay (Busan clay, which is an unusually thick deposit) at several locations in the Nakdong River deltaic plane in Busan, Korea. Two correlation techniques are applied for estimating the cone factors for which the uncorrected and corrected field vane strengths are used as the reference strength. First, due to the scattered data over the entire depth, a direct correlation between the tested results (a method generally used in practice), is applied for each depositional environment (facies) that is detected from the CPT profiles and the geological analysis. Secondly, a special attempt is made to determine the cone factors based on the fact that the ratios of the field vane strength and the (q t  ? σ v0) and (u 2 ? u 0) values to the effective overburden pressure respectively, are approximately constant for each facies (the indirect correlation). It appears that the values obtained from the two techniques are almost identical. However, the direct correlation is sensitively affected by the fabric effects especially in some of the uppermost and lowermost layers of the clay and hence the arithmetic average values of each facies sometimes result in lower correlation coefficients. It is therefore inferred that the indirect correlation is more proficient for estimating the cone factors at each facies and can be used for practical applications. From the results of the indirect correlation, the cone factors of Busan clay are N kt  = 7.13–13.06, N kt,corr  = 9.45–17.07, N Δu  = 6.10–9.88 and N Δu,corr  = 8.41–12.24, the corrected values of which are considerably smaller or larger than the general ranges recommended by Lunne et al. (1997a). However, good correlations between N kt and I P , and between N Δu (or N kt ) and B q , are not found for the clay.  相似文献   

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

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