首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The location, shape and stacking pattern of deep‐marine clastic sediments on drifting stage passive continental margins are strongly influenced by the slope and basin floor topography. The tectonic control on sediment routes and dispersal patterns, however, is less understood on rift margins, particularly the impact of subaqueous transfer zones or relay ramps. In this study, an area of the Palaeocene marine syn‐rift succession in the Vøring Basin is mapped in detail to unravel the relationship between fault geometries and sedimentary infill patterns. Using root‐mean‐square (RMS) amplitudes and deposit thicknesses interpreted from seismic data, sedimentary elements in the Fenris Graben and the Gjallar Ridge are related to the fault patterns and the overall basin geometry. Older deposits are found to be aligned parallel to the basin axis, with the greatest sediment thicknesses on the hanging walls and adjacent to rotated faults. The main sediment supply is interpreted to be sourced from the Vøring Marginal High and Greenland, presumably containing a significant proportion of coarser grained material and comprising numerous local depocentres. With continued rifting and decreased fault activity, finer grained deposition draped the previous basin infill and smoothed the basin floor topography. Deposits close to the foot of relay ramps along the Gjallar Ridge, however, suggest that the high may have acted as a local sediment source leading to local depocentres. Transfer zones played a significant role in sediment transport during the early rifting phase, and were able to maintain some influence into the late rifting and early drifting stage. Identification of early‐ and late‐stage transfer zones may therefore help in locating coarser grained depocentres and potential hydrocarbon reservoirs.  相似文献   

2.
Although the trajectory and geometry of clinoforms in different types of basins have been described in many studies, few studies discuss the influence of halokinesis on clinoforms in salt-related basins. In this study, we analyse the Lower Cretaceous clinoforms in the Tiddlybanken Basin, Norwegian Barents Sea to evaluate the impact of salt mobilization on the geometry and trajectory of clinoforms as well as its implications on sediment partitioning. To accomplish this objective, we use a multidisciplinary approach consisting of seismic and well-interpretation, 3D structural restoration, and forward stratigraphic modelling. The results show that salt mobilization affects prograding clinoforms by: (a) causing lateral variations in progradation rates, resulting in complex palaeogeography, (b) increasing slope angles, which affect the equilibrium of the clinoform profile and can trigger slope-readjustment processes and (c) producing lateral and temporal variations in accommodation space, leading to different clinoform trajectories, stacking patterns and reservoir distribution along the basin. Forward stratigraphic modelling shows that in salt-related basins and other tectonically active basins, the isolated use of conventional methods for clinoform analysis might lead to potential interpretation pitfalls such as misinterpretation of trajectories and overestimation of foreset angles, which can have negative consequences for exploration models.  相似文献   

3.
《Basin Research》2018,30(2):279-301
Spatio‐temporal analysis of basins formed along sheared margins has received much less attention than those formed along orthogonally extended margins. Knowledge about the structural evolution of such basins is important for petroleum exploration but there has been a lack of studies that document these based on 3D seismic reflection data. In this study, we demonstrate how partitioning of strain during deformation of the central and southern part of the Sørvestsnaget Basin along the Senja Shear Margin, Norwegian Barents Sea, facilitated coeval shortening and extension. This is achieved through quantitative analysis of syn‐kinematic growth strata using 3D seismic data. Our results show that during Cenozoic extensional faulting, folds and thrusts developed coevally and orthogonal to sub‐orthogonal to normal faults. We attribute this strain partitioning to be a result of the right‐lateral oblique plate motions along the margin. Rotation of fold hinge‐lines and indications of hinge‐parallel extension indicate that the dominating deformation mechanism in the central and southern Sørvestsnaget Basin during opening along the Senja Shear Margin was transtensional. We also argue that interpretation of shortening structures attributed to inversion along the margin should consider that partitioning of strain may result in shortening structures that are coeval with extensional faults and not a result of a separate compressional phase.  相似文献   

4.
The nature of the Lower Crustal Body (LCB) underneath the western part of the Vøring margin (NE Atlantic) is studied with three scenarios of its extension history: (a) The LCB is Caledonian crust. (b) Half the LCB is Caledonian crust and the other half is emplaced as magmatic underplating in Late Palaeocene. (c) The entire LCB is emplaced as magmatic underplating. The extension of the margin transect is obtained with a procedure that accounts for the extension and thinning of the sedimentary basins. This procedure has been extended to include magmatic underplating. The lithosphere is modelled with deposition of sediments and four rift phases since the Early Devonian until today. The forward modelling is mass conservative and the present‐day thicknesses of the formations, crust, LCB and magmatic underplate are reproduced. The state of the lithosphere and the sedimentary basins are shown and compared at the beginning and at the end of the rift phases. It is concluded that the scenario with the LCB as only underplating requires an unrealistic amount of extension. A scenario where underplating accounts for maximum half the LCB is more likely. Two different interpretations for the Moho underneath the Utgard High are tested: one with a shallow base‐crust and another with a deeper and flatter base‐crust. Tectonic modelling of the two versions favours the latter interpretation. The modelling shows that the Late Jurassic rift phase was much more prominent than the Late Cretaceous and Palaeocene rift phase for all cases of underplating. A strong Late Jurassic rift phase is consistent with the accumulation space needed for the thick Cretaceous formations. There are no observations of magmatism from the Late Jurassic, although this rift phase is stronger than the Cretaceous and Palaeocene rift phase.  相似文献   

5.
《Basin Research》2018,30(5):942-964
Mass wasting is triggered on many continental slopes by a number of mechanisms, including seismic shaking, high sedimentation rates, the presence of weak geological units and gas hydrate dissociation. In this study, the morphology of a Late Miocene–Early Pliocene mass‐transport complex (MTC) on the Utgard High is unravelled and discussed in relation to possible trigger mechanisms. The approach used here includes 3D seismic interpretation and the analysis of variance attribute maps. The interpreted MTC is located on the crest and flanks of the Utgard High and is composed of three mass‐transport deposits with seismic characters varying from transparent and chaotic seismic facies at the base to slightly deformed layers composed of mounds and rafted blocks in the middle and chaotic to transparent reflections at the top. Lithologically, the MTC consists predominantly of claystone with high gamma ray and low density and resistivity values, demonstrating that the associated mounds represent remobilized ooze sediments. A vertical stack of six magmatic sills emplaced from 55.6 to 56.3 Ma into the Upper Cretaceous shales is interpreted at depths of 3,000–5,500 ms two‐way travel time (TWTT). In association with these magmatic sills are several hydrothermal vent complexes that interacted with the top MTC horizon, signifying that episodic and secondary fluid‐venting events might be the principal mechanism facilitating mass wasting in the study area. In addition, the remobilization of ooze sediments into mounds is hypothesized to be dependent on fluids and clayey layers. As a corollary of this work, the importance of relict and recurrent episodes of fluid flow in the Vøring Basin and their influence on the geotechnical integrity of the overburden and later mass wasting is established.  相似文献   

6.
An extensive, reprocessed two‐dimensional (2D) seismic data set was utilized together with available well data to study the Tiddlybanken Basin in the southeastern Norwegian Barents Sea, which is revealed to be an excellent example of base salt rift structures, evaporite accumulations and evolution of salt structures. Late Devonian–early Carboniferous NE‐SW regional extensional stress affected the study area and gave rise to three half‐grabens that are separated by a NW‐SE to NNW‐SSE trending horst and an affiliated interference transfer zone. The arcuate nature of the horst is believed to be the effect of pre‐existing Timanian basement grain, whereas the interference zone formed due to the combined effect of a Timanian (basement) lineament and the geometrical arrangement of the opposing master faults. The interference transfer zone acted as a physical barrier, controlling the facies distribution and sedimentary thickness of three‐layered evaporitic sequences (LES). During the late Triassic, the northwestern part of a salt wall was developed due to passive diapirism and its evolution was influenced by halite lithology between the three‐LES. The central and southeastern parts of the salt wall did not progress beyond the pedestal stage due to lack of halite in the deepest evaporitic sequence. During the Triassic–Jurassic transition, far‐field stresses from the Novaya Zemlya fold‐and‐thrust belt reactivated the pre‐salt Carboniferous rift structures. The reactivation led to the development of the Signalhorn Dome, rejuvenated the northwestern part of the salt wall and affected the sedimentation rates in the southeastern broad basin. The salt wall together with the Signalhorn Dome and the Carboniferous pre‐salt structures were again reactivated during post‐Early Cretaceous, in response to regional compressional stresses. During this main tectonic inversion phase, the northwestern and southeastern parts of the salt wall were rejuvenated; however, salt reactivation was minimized towards the interference transfer zone beneath the centre of the salt wall.  相似文献   

7.
The Porcupine Basin is a Mesozoic failed rift located in the North Atlantic margin, SW of Ireland, in which a postrift phase of extensional faulting and reactivation of synrift faults occurred during the Mid–Late Eocene. Fault zones are known to act as either conduits or barriers for fluid flow and to contribute to overpressure. Yet, little is known about the distribution of fluids and their relation to the tectono‐stratigraphic architecture of the Porcupine Basin. One way to tackle this aspect is by assessing seismic (Vp) and petrophysical (e.g., porosity) properties of the basin stratigraphy. Here, we use for the first time in the Porcupine Basin 10‐km‐long‐streamer data to perform traveltime tomography of first arrivals and retrieve the 2D Vp structure of the postrift sequence along a ~130‐km‐long EW profile across the northern Porcupine Basin. A new Vp–density relationship is derived from the exploration wells tied to the seismic line to estimate density and bulk porosity of the Cenozoic postrift sequence from the tomographic result. The Vp model covers the shallowest 4 km of the basin and reveals a steeper vertical velocity gradient in the centre of the basin than in the flanks. This variation together with a relatively thick Neogene and Quaternary sediment accumulation in the centre of the basin suggests higher overburden pressure and compaction compared to the margins, implying fluid flow towards the edges of the basin driven by differential compaction. The Vp model also reveals two prominent subvertical low‐velocity bodies on the western margin of the basin. The tomographic model in combination with the time‐migrated seismic section shows that whereas the first anomaly spatially coincides with the western basin‐bounding fault, the second body occurs within the hangingwall of the fault, where no major faulting is observed. Porosity estimates suggest that this latter anomaly indicates pore overpressure of sandier Early–Mid Eocene units. Lithological well control together with fault displacement analysis suggests that the western basin‐bounding fault can act as a hydraulic barrier for fluids migrating from the centre of the basin towards its flanks, favouring fluid compartmentalization and overpressure of sandier units of its hangingwall.  相似文献   

8.
The Crotone Basin was generated in the late Cenozoic as a forearc basin of the Ionian arc‐trench system. New data are gained through detailed field mapping, high‐resolution stratigraphic analysis of a key area and examination of offshore well data and seismic reflection profiles. Major unconformities divide the basin fill into major sequences, which reveal a three‐stage internal organization thought to reflect geodynamic events of the Calabrian arc and backarc area closely. The first stage is characterized by extensional block faulting and uplift followed by rapid drowning during high subsidence and transtension in the basin along a major NNW‐ to NW‐striking fault system. This stage is interpreted to reflect resumption of rollback after an episode of slab tearing triggered by transitory docking of continental lithosphere in the trench. The initial uplift is inferred to reflect decoupling and rebound after the transitory coupling phase. The second stage is characterized by increased subsidence and continued extension/transtension. This trend presumably reflects a decreasing rate of rollback resulting from a tendency towards viscous coupling after acceleration of slab downwelling. The third stage is characterized by short‐lived transpression along major shear zones and local inversion of former basins. This is inferred to reflect entrance into the trench of buoyant continental lithosphere, resulting in significant deceleration of slab rollback and consequently a break in, or slowing of, backarc extension, and predominance of the effects of compression related to Africa–Europe convergence. Overall, the above evolution resulted in the formation of a progressively narrower and rapidly retreating slab, inducing extreme rates of backarc extension, and may have played a critical role in determining the intermittent nature of the backarc rifting.  相似文献   

9.
10.
Through examination of the scaling relations of faults and the use of seismic stratigraphic techniques, we demonstrate how the temporal and spatial evolution of the fault population in a half-graben basin can be accurately reconstructed. The basin bounded by the ≫62-km-long Strathspey–Brent–Statfjord fault array is located on the western flank of the Late Jurassic northern North Sea rift basin. Along-strike displacement variations, transverse fault-displacement folds and palaeo-fault tips abandoned in the hangingwall all provide evidence that the fault system comprises a hierarchy of linked palaeo-segments. The displacement variations developed while the fault was in a prelinkage, multisegment stage of its growth have not been equilibrated following fault linkage. Using the stratal architecture of synrift sediments, we date the main phase of segment linkage as latest Callovian – middle Oxfordian (10–14 Myr after rift initiation). A dense subpopulation of faults is mapped in the hangingwall to the Strathspey–Brent–Statfjord fault array. The majority of these faults are short, of low displacement and became inactive within 3–4 Myr of the beginning of the extensional event. Subsequently, only the segments of the proto-Strathspey–Brent–Statfjord fault and a conjugate array of antithetic faults located 3.5 km basinward continued to grow to define a graben-like basin geometry. Faults of the antithetic array became inactive ∼11.5 Myr into the rift event, concentrating strain on the linked Strathspey–Brent–Statfjord fault; hence, the basin evolved into a half-graben. As the rift event progressed, strain was localized on a smaller number of active structures with increased rates of displacement. The results of this study suggest that a simple model for the linkage of 2–3 fault segments may not be applicable to a complex multisegment array.  相似文献   

11.
Significant lateral variations in observed temperatures in the Beaufort‐Mackenzie Basin raise the question on the temperature‐controlling factors. Based on the structural configuration of the sediments and underlying crust in the area, we calculate the steady‐state 3D conductive thermal field. Integrated data include the base of the relic permafrost layer representing the 0 °C‐isotherm, public‐domain temperature data (from 227 wells) and thermal conductivity data. For >75% of the wells the predicted temperatures deviate by <10 K from the observed temperatures, which validates the overall model setup and adopted thermal properties. One important trend reproduced by the model is a decrease in temperatures from the western to the eastern basin. While in the west, a maximum temperature of 185 °C is reached at 5000 m below sea level, in the east the maximum temperature is 138 °C. The main cause for this pattern lies in lateral variations in thermal conductivity indicating differences in the shale and sand contents of the different juxtaposed sedimentary units. North‐to‐south temperature trends reveal the superposition of deep and shallow effects. At the southern margin, where the insulating effect of the low‐conductive sediments is missing, temperatures are lowest. Farther north, where the sub‐sedimentary continental crust is thick enough to produce considerable heat and a thick pile of sediments efficiently stores heat, temperatures tend to be highest. Temperatures decrease again towards the northernmost distal parts of the basin, where thinned continental and oceanic crust produce less radiogenic heat. Wells with larger deviations of the purely conductive model from the temperature observations (>15 K at 10% of the wells) and their basin‐wide pattern of misfit tendency (too cold vs. too warm temperature predictions) point to a locally restricted coupling of heat transport to groundwater flow.  相似文献   

12.
We investigate the effects of the cooling of intrusive and extrusive igneous bodies on the temperature history and surface heat flow of the Paraná Basin. The Serra Geral igneous event (130–135 Ma) covered most of this basin with flood basalts. Associated with this event numerous sills and dykes intruded the sediments and basement, and extensive underplating may have occurred in the lower crust and upper mantle beneath the basin. We develop an analytical model of the conductive cooling of tabular intrusive bodies and use it to calculate temperatures within the sediments as a function of time since emplacement. Depending on the thickness of these igneous bodies and the timing of sequential emplacement, the thermal history of a given locus in the basin can range from a simple extended period of higher temperatures to multiple episodes of peak temperatures separated by cooling intervals. The cooling of surface flood basalts, sills and dykes is capable of maintaining temperatures above the normal geothermal gradient temperatures for a few hundred thousand years, while large-scale underplating may influence temperatures for up to 10 million years. We conclude that any residual heat from the cooling of the Serra Geral igneous rocks has long since decayed to insignificant values and that present-day temperatures and heat flow are not affected. However, the burial of the sediments beneath the thick basalt cap caused a permanent temperature increase of up to 50°C in the underlying sediments since the beginning of the Cretaceous.  相似文献   

13.
14.
Sediments from the last interglacial, marine isotope stage 5e (MIS 5e), have been studied for their dinoflagellate cyst content in a core retrieved from the Vøring Plateau, Norwegian Sea. Qualitative and quantitative analyses of the data, and comparison with the surface sample and published Holocene data from the core, reveal distinct differences in hydrological surface conditions between the late Holocene and MIS 5e. A higher number of co-dominant, subordinate species in the last interglacial samples suggests there was a more pronounced seasonality of the surface water at this time. This is supported by the significant presence of Bitectatodinium tepikiense , a species that was virtually absent from the area for most of the Holocene. The seasonality signal is further substantiated by transfer-function reconstructions, which also indicates a stronger stratification of the upper water column during MIS 5e. Moreover, the assemblage data clearly show that optimal, fully marine interglacial conditions prevailed only late in MIS 5e (between ca. 117.5 and 116.5 Kya), which is in contrast with the climatic optimum early in the Holocene. Stable oxygen isotope values from planktic foraminifera for this MIS 5e optimum are comparable with the average Holocene values, but are generally ca. 0.3‰ higher than those of the earlier part of the last interglacial ( sensu stricto ). These higher δ18O values are likely to be the result of the enhanced and prolonged influence of Saalian deglacial meltwater, thus corroborating the existence of a quite differently structured sea surface, as suggested by the dinocyst data.  相似文献   

15.
In this article, we document a large number of focused fluid escape structures using high quality 2D seismic reflection data and multibeam bathymetry data from a poorly known area at the intersection of the northern South China Sea (SCS) and the western SCS. Three types of focused fluid escape systems are identified and described: mud volcanoes, pipes and associated pockmarks. The mud volcanoes occur singly or as clustered groups. The overpressure driving the mud volcanism is argued herein to be related to the generation of thermogenic hydrocarbons. The clustered distributions are related to localized tectonic uplift in the basin. Pipes mainly occur within the Guangle Uplift or accompany the mud volcano clusters. The pipes located within the Guangle Uplift are attributed to carbonate dissolution caused by hydrothermal fluids. Fluids ascended through these structures and were expelled at the palaeo‐seabed or present seabed forming palaeo‐pockmarks and present‐day pockmarks. Some ‘mega‐pockmarks’ show evidence of enlargement due to bottom currents. The marginal basins of the SCS are petroliferous, with attention gradually shifting to the deep‐water area. Our results show that fluid migration must be taken into account when assessing seabed stability. This analysis also improves our understanding the petroleum geology in the study area, and is also useful for predicating where chemosynthetic ecosystems may be located.  相似文献   

16.
Mantle-induced dynamic topography (i.e., subsidence and uplift) has been increasingly recognized as an important process in foreland basin development. However, characterizing and distinguishing the effects (i.e., location, extent and magnitude) of dynamic topography in ancient foreland basins remains challenging because the spatio-temporal footprint of dynamic topography and flexural topography (i.e., generated by topographic loading) can overlap. This study employs 3D flexural backstripping of Upper Cretaceous strata in the central part of the North American Cordilleran foreland basin (CFB) to better quantify the effects of dynamic topography. The extensive stratigraphic database and good age control of the CFB permit the regional application of 3D flexural backstripping in this basin for the first time. Dynamic topography started to influence the development of the CFB during the late Turonian to middle Campanian (90.2–80.2 Ma) and became the dominant subsidence mechanism during the middle to late Campanian (80.2–74.6 Ma). The area influenced by >100 m dynamic subsidence is approximately 400 by 500 km, within which significant (>200 m) dynamic subsidence occurs in an irregular-shaped (i.e., lunate) subregion. The maximum magnitude of dynamic subsidence is 300 ± 100 m based on the 80.2–74.6 Ma tectonic subsidence maps. With the maximum magnitude of dynamic uplift being constrained to be 200–300 m, the gross amount of dynamic topography in the Late Cretaceous CFB is 500–600 m. Although the location of dynamic subsidence revealed by tectonic subsidence maps is generally consistent with isopach map trends, tectonic subsidence maps developed through 3D flexural backstripping provide more accurate constraints of the areal extent, magnitude and rate of dynamic topography (as well as flexural topography) in the CFB through the Late Cretaceous. This improved understanding of dynamic topography in the CFB is critical for refining current geodynamic models of foreland basins and understanding the surface expression of mantle processes.  相似文献   

17.
We have analysed the fundamental mode of Love and Rayleigh waves generated by 12 earthquakes located in the mid-Atlantic ridge and Jan Mayen fracture zone. Using the multiple filter analysis technique, we isolated the Rayleigh and Love wave group velocities for periods between 10 and 50  s. The surface wave propagation paths were divided into five groups, and average group velocities calculated for each group. The average group velocities were inverted and produced shear wave velocity models that correspond to a quasi-continental oceanic structure in the Greenland–Norwegian Sea region. Although resolution is poor at shallow depth, we obtained crustal thickness values of about 18  km in the Norwegian Sea area and 9  km in the region between Svalbard and Iceland. The abnormally thick crust in the Norwegian Sea area is ascribed to magmatic underplating and the thermal blanketing effect of sedimentary layers. Maximum crustal shear velocities vary between 3.5 and 3.9  km  s−1 for most paths. An average lithospheric thickness of 60  km was observed, which is lower than expected for oceanic-type structure of similar age. We also observed low shear wave velocities in the lower crust and upper mantle. We suggest that high heat flow extending to depths of about 30  km beneath the surface can account for the thin lithosphere and observed low velocities. Anisotropy coefficients of 1–5 per cent in the shallow layers and >7 per cent in the upper mantle point to the existence of polarization anisotropy in the region.  相似文献   

18.
19.
In this study, we integrate 3D seismic reflection, wireline log, biostratigraphic and core data from the Egersund Basin, Norwegian North Sea to determine the impact of syn‐depositional salt movement and associated growth faulting on the sedimentology and stratigraphic architecture of the Middle‐to‐Upper Jurassic, net‐transgressive, syn‐rift succession. Borehole data indicate that Middle‐to‐Upper Jurassic strata consist of low‐energy, wave‐dominated offshore and shoreface deposits and coal‐bearing coastal‐plain deposits. These deposits are arranged in four parasequences that are aggradationally to retrogradationally stacked to form a net‐transgressive succession that is up to 150‐m thick, at least 20 km in depositional strike (SW‐NE) extent, and >70 km in depositional dip (NW‐SE) extent. In this rift‐margin location, changes in thickness but not facies are noted across active salt structures. Abrupt facies changes, from shoreface sandstones to offshore mudstones, only occur across large displacement, basement‐involved normal faults. Comparisons to other tectonically active salt‐influenced basins suggest that facies changes across syn‐depositional salt structures are observed only where expansion indices are >2. Subsidence between salt walls resulted in local preservation of coastal‐plain deposits that cap shoreface parasequences, which were locally removed by transgressive erosion in adjacent areas of lower subsidence. The depositional dip that characterizes the Egersund Basin is unusual and likely resulted from its marginal location within the evolving North Sea rift and an extra‐basinal sediment supply from the Norwegian mainland.  相似文献   

20.
We propose and test a conceptual framework for evaluating the relative timing of different types of sedimentary indicators of tectonism in alluvial foreland basin settings. We take the first occurrence of a detrital grain from a newly exposed source‐area lithology to provide the best indicator of the onset of tectonic uplift in the source area. Source‐area unroofing may lag behind initial uplift because of the type, thickness and structure of rocks in the uplifted mountain block, drainage patterns and climate. However, once exposed, advective transport disperses grains quickly throughout fluvial systems. Because of increased subsidence rate from thrust belt loading, an increase in sedimentation rate begins coincident with tectonic load emplacement within the flexural half‐width of the basin. However, farther out into the basin increased sedimentation rates lag behind the composition signal because of time lags associated with propagation of the thrust load and attendant sediment loads into the basin. The progradation of syntectonic gravel lags behind all of these signals as a direct function of the relative proportion of gravel fraction within transported sediment and rates and geometry of subsidence, which selectively traps the coarsest grain‐size fractions in the most proximal parts of the basin. We demonstrate this signal attenuation in the syntectonic Horta–Gandesa alluvial system (late Eocene–Oligocene), exposed along the southeast margin of the Ebro Basin, Spain. The results demonstrate that: (1) the time spans between the compositional signal and the progradation of the gravel front can be geologically significant, on the order of more than a million years within as little as 20 km of the thrust front; and (2) time lags between the signals increase with distance away from the deformation front. No lag time was observed between the first appearance of a new clast composition and the arrival of gravel front when the thrust front was within a few tens of metres from the depositional site. In contrast, the time lag was 0.5–1 Myr when the thrust front was about 5–6 km away and it increased to >1 Myr when the deformation front was about 8 km away. At the most extreme position, when the thrust front was 15–20 km away, the gravel front never reached the study area.  相似文献   

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

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