首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
This paper presents a Bayesian methodology for updating the seismic hazard curves. The methodology is based on the comparison of predictive exceedance rates of a fixed acceleration level (given by the seismic hazard curves) and the observed exceedance rates in some selected sites. The application of the methodology needs, firstly, the definition of a prior probabilistic seismic hazard assessment based in a logic tree. Each main branch corresponds to a probabilistic model of calculus of seismic hazard. The method considers that, initially (or a priori), the weights of all branches of the logic tree are equivalent. Secondly, the method needs to compile the observations in the region. They are introduced in a database containing the recorded acceleration data (during the instrumental period). Nevertheless, the instrumental period in stable zones (as France) shows only very low acceleration levels recorded during a short observation period. Then, a method to enlarge the REX (number of observations) is presented taking into account the historical data and defining “synthetic” accelerations in the sites of observation. The synthetic REX allows to expand the period of observation and to increase the acceleration thresholds used in the Bayesian updating process. The application of the Bayesian approach leads to a new and more objective definition of the weights of each branch of the logic tree and, therefore, to new seismic hazard curves (mean and centiles). The Bayesian approach doesn’t change the probabilistic models (seismic hazard curves). It only modifies the weights of each branch of the logic tree.  相似文献   

2.
-- We investigate the impact of different rupture and attenuation models for the Cascadia subduction zone by simulating seismic hazard models for the Pacific Northwest of the U.S. at 2% probability of exceedance in 50 years. We calculate the sensitivity of hazard (probabilistic ground motions) to the source parameters and the attenuation relations for both intraslab and interface earthquakes and present these in the framework of the standard USGS hazard model that includes crustal earthquakes. Our results indicate that allowing the deep intraslab earthquakes to occur anywhere along the subduction zone increases the peak ground acceleration hazard near Portland, Oregon by about 20%. Alternative attenuation relations for deep earthquakes can result in ground motions that differ by a factor of two. The hazard uncertainty for the plate interface and intraslab earthquakes is analyzed through a Monte-Carlo logic tree approach and indicates a seismic hazard exceeding 1 g (0.2 s spectral acceleration) consistent with the U.S. National Seismic Hazard Maps in western Washington, Oregon, and California and an overall coefficient of variation that ranges from 0.1 to 0.4. Sensitivity studies indicate that the paleoseismic chronology and the magnitude of great plate interface earthquakes contribute significantly to the hazard uncertainty estimates for this region. Paleoseismic data indicate that the mean earthquake recurrence interval for great earthquakes is about 500 years and that it has been 300 years since the last great earthquake. We calculate the probability of such a great earthquake along the Cascadia plate interface to be about 14% when considering a time-dependent model and about 10% when considering a time-independent Poisson model during the next 50-year interval.  相似文献   

3.
Microtremor Measurements for the Microzonation of Dinar   总被引：3，自引：0，他引：3  
v--vThe geotechnical site conditions in Dinar town located in western Turkey were investigated after the 1995 Dinar earthquake based on borings, in situ penetration tests, seismic wave velocity measurements, and microtremor records. The variation of damage distribution within the town was evaluated with respect to 23 district damage ratios calculated, based on the detailed damage survey conducted by the General Directorate of Disaster Affairs. Site amplifications were estimated from microtremor spectral ratios and microzonation was performed using a GIS methodology. The results of in situ penetration tests and seismic wave velocity measurements as well as the damage distribution were compared with the amplification zonation obtained from microtremor records. The results indicate the applicability of microtremor spectral ratios for assessing the local site conditions and site amplifications.  相似文献   

4.
This study describes the methodology implemented to establish the ground-motion logic-tree for national probabilistic seismic hazard map of Turkey for shallow active crustal regions. The presented procedure provides quantitative information to guide the hazard experts while establishing the logic tree to capture the epistemic uncertainty in ground-motion characterization. It uses non-data-driven and data-driven testing methods to identify and rank candidate ground-motion prediction equations (GMPEs) under a specific ground-motion database. The candidate GMPEs are subjected to visual inspection and are classified into center, body and range (CBR) spectral estimates for a proper consideration of epistemic uncertainty. The GMPEs classified into CBR are then used in a suite of seismic hazard sensitivity analysis to establish the most suitable GMPE logic-tree whose spectral estimates are not biased by any one of the GMPEs in the logic-tree structure. The sensitivity analysis considers normalized spectral ordinates and is not manipulated by the spectral amplitudes. The proposed procedure is inherited from the relevant studies of the Earthquake Model of the Middle East (EMME; www.efehr.org:8080/jetspeed/portal/emme.psml) regional seismic hazard project. This paper also highlights the similarities and differences in ground-motion characterization between EMME and our approach.  相似文献   

5.
A representation of seismic hazard is proposed for Italy based on the zone-free approach developed by Woo (BSSA 86(2):353–362, 1996a), which is based on a kernel estimation method governed by concepts of fractal geometry and self-organized seismicity, not requiring the definition of seismogenic zoning. The purpose is to assess the influence of seismogenic zoning on the results obtained for the probabilistic seismic hazard analysis (PSHA) of Italy using the standard Cornell’s method. The hazard has been estimated for outcropping rock site conditions in terms of maps and uniform hazard spectra for a selected site, with 10 % probability of exceedance in 50 years. Both spectral acceleration and spectral displacement have been considered as ground motion parameters. Differences in the results of PSHA between the two methods are compared and discussed. The analysis shows that, in areas such as Italy, characterized by a reliable earthquake catalog and in which faults are generally not easily identifiable, a zone-free approach can be considered a valuable tool to address epistemic uncertainty within a logic tree framework.  相似文献   

6.
任梦依  刘哲 《地震学报》2022,44(6):1035-1048
基于广义帕累托分布构建地震活动性模型,因其输入参数取值难以避免不确定性,导致依据该模型所得的地震危险性估计结果具有不确定性。鉴于此,本文选取青藏高原东北缘为研究区,提出了基于全域敏感性分析的地震危险性估计的不确定性分析流程和方法。首先,利用地震活动性广义帕累托模型,进行研究区地震危险性估计;然后,选取地震记录的起始时间和震级阈值作为地震活动性模型的输入参数,采用具有全域敏感性分析功能的E-FAST方法,对上述两个参数的不确定性以及两参数之间的相互作用对地震危险性估计不确定性的影响进行定量分析。结果表明：地震危险性估计结果（不同重现期的震级重现水平、震级上限及相应的置信区间）对两个输入参数中的震级阈值更为敏感;不同重现期的地震危险性估计结果对震级阈值的敏感程度不同;对不同的重现期而言,在影响地震危险性估计结果的不确定性上,两个输入参数之间存在非线性效应,且非线性效应程度不同。本文提出的不确定性分析流程和方法,可以推广应用于基于其它类型地震活动性模型的地震危险性估计不确定性分析。  相似文献   

7.
The Seismic Hazard Harmonization in Europe (SHARE) project, which began in June 2009, aims at establishing new standards for probabilistic seismic hazard assessment in the Euro-Mediterranean region. In this context, a logic tree for ground-motion prediction in Europe has been constructed. Ground-motion prediction equations (GMPEs) and weights have been determined so that the logic tree captures epistemic uncertainty in ground-motion prediction for six different tectonic regimes in Europe. Here we present the strategy that we adopted to build such a logic tree. This strategy has the particularity of combining two complementary and independent approaches: expert judgment and data testing. A set of six experts was asked to weight pre-selected GMPEs while the ability of these GMPEs to predict available data was evaluated with the method of Scherbaum et al. (Bull Seismol Soc Am 99:3234?C3247, 2009). Results of both approaches were taken into account to commonly select the smallest set of GMPEs to capture the uncertainty in ground-motion prediction in Europe. For stable continental regions, two models, both from eastern North America, have been selected for shields, and three GMPEs from active shallow crustal regions have been added for continental crust. For subduction zones, four models, all non-European, have been chosen. Finally, for active shallow crustal regions, we selected four models, each of them from a different host region but only two of them were kept for long periods. In most cases, a common agreement has been also reached for the weights. In case of divergence, a sensitivity analysis of the weights on the seismic hazard has been conducted, showing that once the GMPEs have been selected, the associated set of weights has a smaller influence on the hazard.  相似文献   

8.
Located on the Alpine-Himalayan earthquake belt, Iran is one of the seismically active regions of the world. Northern Iran, south of Caspian Basin, a hazardous subduction zone, is a densely populated and developing area of the country. Historical and instrumental documented seismicity indicates the occurrence of severe earthquakes leading to many deaths and large losses in the region. With growth of seismological and tectonic data, updated seismic hazard assessment is a worthwhile issue in emergency management programs and long-term developing plans in urban and rural areas of this region. In the present study, being armed with up-to-date information required for seismic hazard assessment including geological data and active tectonic setting for thorough investigation of the active and potential seismogenic sources, and historical and instrumental events for compiling the earthquake catalogue, probabilistic seismic hazard assessment is carried out for the region using three recent ground motion prediction equations. The logic tree method is utilized to capture epistemic uncertainty of the seismic hazard assessment in delineation of the seismic sources and selection of attenuation relations. The results are compared to a recent practice in code-prescribed seismic hazard of the region and are discussed in detail to explore their variation in each branch of logic tree approach. Also, seismic hazard maps of peak ground acceleration in rock site for 475- and 2,475-year return periods are provided for the region.  相似文献   

9.
-- The study addresses the evaluation of earthquake hazard parameters such as maximum regional magnitude (Mmax) and the slope of Gutenberg-Richter law # (where b=# log e) for the Hellenic Wadati-Benioff zone and the overriding lithospheric plate in the area of Crete and its surroundings. The seismicity of the area is divided in a cellular (1.0° 2 1.0°) manner allowing analysis of the localized earthquake hazard parameters and graphical representation of their spatial variation. Our approach incorporates the recently updated earthquake catalogue for Greece and the adjacent areas, the consideration of the morphology of the deep seismically active structures in the studied area and use of a probabilistic procedure for estimating the earthquake hazard parameters.¶One of the main inconsistencies in the earthquake hazard assessment is the estimation of the maximum magnitude and the related uncertaint y. The Bayesian approach, applied in the present, is a straightforward technique for evaluating the earthquake hazard parameters and is based on the following assumptions: Poissonian character of seismic events flow, a frequency-magnitude law of Gutenberg-Richter's type with cutoff maximal value for estimated parameter and a seismic catalogue, having a rather sizeable number of events (i.e., 50 events at least per cell). For five cells in which the number of events is less than 50, an effort is made to produce synthetic data. The re-assessed parameters obtained from the synthetic data show no significant difference and the real data (of the five cells) are finally taken into account although the estimated uncertainty is high.¶For four random cells we constructed hazard curves showing the probabilities that a certain magnitude M will be exceeded in one year and the return periods (in years) that are expected for a given magnitude. These are particularly useful for the mapping of earthquake hazard in regions of either low or high seismic activity, as is Crete and the adjacent area.¶The obtained results show that the W and E parts of both subducting and overriding plates differ in the spatial distribution of all the estimated earthquake hazard parameters. The Mmax distribution indicates strong coupling between the western portions of the interacting plates (Mmax > 6.3) to the south of 36°N. The smaller values of Mmax (Mmax < 6.3) estimated in the SE part of the studied area indicate weak coupling between the eastern portions of the subducting and overriding plates.¶Values of b > 1.0 are found to the south and east of Crete for the Wadati-Benioff zone, and over the central part of the island and the area to the northeast of it (cell 11) for the continental wedge, which suggests nonuniform stress field and/or heterogeneous material.  相似文献   

10.
v--v Continuous seismic threshold monitoring is a technique that has been developed over the past several years to assess the upper magnitude limit of possible seismic events that might have occurred in a geographical target area. The method provides continuous time monitoring at a given confidence level, and can be applied in a site-specific, regional or global context.¶In this paper (Part 1) and a companion paper (Part 2) we address the problem of optimizing the site-specific approach in order to achieve the highest possible automatic monitoring capability of particularly interesting areas. The present paper addresses the application of the method to cases where a regional monitoring network is available. We have in particular analyzed events from the region around the Novaya Zemlya nuclear test site to develop a set of optimized processing parameters for the arrays SPITS, ARCES, FINES, and NORES. From analysis of the calibration events we have derived values for beam-forming steering delays, filter bands, short-term average (STA) lengths, phase travel times (P and S waves), and amplitude-magnitude relationships for each array. By using these parameters for threshold monitoring of the Novaya Zemlya testing area, we obtain a monitoring capability varying between mb 2.0 and 2.5 during normal noise conditions.¶The advantage of using a network, rather than a single station or array, for monitoring purposes becomes particularly evident during intervals with high global seismic activity (aftershock sequences), high seismic noise levels (wind, water waves, ice cracks) or station outages. For the time period November-December 1997, all time intervals with network magnitude thresholds exceeding mb 2.5 were visually analyzed, and we found that all of these threshold peaks could be explained by teleseismic, regional, or local signals from events outside the Novaya Zemlya testing area. We could therefore conclude within the confidence level provided by the method, that no seismic event of magnitude exceeding 2.5 occurred at the Novaya Zemlya test site during this two-month time interval.¶As an example of particular interest in a monitoring context, we apply optimized threshold processing of the SPITS array for a time interval around 16 August 1997 mb 3.5 event in the Kara Sea. We show that this processing enables us to detect a second, smaller event from the same site (mb 2.6), occurring about 4 hours later. This second event was not defined automatically by standard processing.  相似文献   

11.
Historical seismic data and seismogenic information are quite scarce for the low seismicity region, and modeling the parameters uncertainties based on probabilistic model is suspicious. The convex set theory-based seismic hazard analysis approach is proposed. The uncertainties of b value, the annual occurrence rate v and the upper bound magnitude Mu are described by the envelop bound convex model and the ellipsoidal bound convex model. Convex analysis method and China probabilistic seismic hazard analysis methodology are combined to perform a bound seismic hazard analysis for Ningbo city, China. The seismic intensity interval obtained using the bound seismic hazard analysis is compared with that calculated using China probabilistic seismic hazard analysis methodology. The sensitivity analysis indicates that the interval of seismic intensity is most sensitive to the annual occurrence rate v. Furthermore, the different convex models have little effect on the interval of seismic intensity.  相似文献   

12.
The basic seismic load parameters for the upcoming national design regulation for DIN EN 1998-1/NA result from the reassessment of the seismic hazard supported by the German Institution for Civil Engineering (DIBt). This 2016 version of the national seismic hazard assessment for Germany is based on a comprehensive involvement of all accessible uncertainties in models and parameters and includes the provision of a rational framework for integrating ranges of epistemic uncertainties and aleatory variabilities in a comprehensive and transparent way. The developed seismic hazard model incorporates significant improvements over previous versions. It is based on updated and extended databases, it includes robust methods to evolve sets of models representing epistemic uncertainties, and a selection of the latest generation of ground motion prediction equations. The new earthquake model is presented here, which consists of a logic tree with 4040 end branches and essential innovations employed for a realistic approach. The output specifications were designed according to the user oriented needs as suggested by two review teams supervising the entire project. Seismic load parameters, for rock conditions of \(v_{S30}\) = 800 m/s, are calculated for three hazard levels (10, 5 and 2% probability of occurrence or exceedance within 50 years) and delivered in the form of uniform hazard spectra, within the spectral period range 0.02–3 s, and seismic hazard maps for peak ground acceleration, spectral response accelerations and for macroseismic intensities. Results are supplied as the mean, the median and the 84th percentile. A broad analysis of resulting uncertainties of calculated seismic load parameters is included. The stability of the hazard maps with respect to previous versions and the cross-border comparison is emphasized.  相似文献   

13.
Modern earthquake loss models make use of earthquake catalogs relevant to the seismic hazard assessment upon seismicity and seismotectonic analysis. The main objective of this paper is to investigate a recently compiled catalog (National Institute of Meteorology or INM catalog: 412-2011) and to generate seismic hazard maps through classical probabilistic seismic hazard assessment (PSHA) and smoothed-gridded seismicity models for Tunisia. It is now established with the local earthquake bulletin that the recent seismicity of Tunisia is sparse and moderate. Therefore, efforts must be undertaken to elaborate a robust hazard analysis for risk assessment and seismic design purposes. These recommendations follow the recently published reports by the World Bank that describe the seismic risk in Tunis City as being beyond a tolerable level with an MSK intensity level of VII. Some attempts were made during the past two decades to assess the seismic hazard for Tunisia and they have mostly failed to properly investigate the historical and instrumental seismicity catalog. This limitation also exists for the key aspect of epistemic and random uncertainties impact on the final seismic hazard assessment. This study also investigates new ground motion prediction equations suitable for use in Tunisia. The methodology applied herein uses, for the first time in PSHA of Tunisia, seismicity parameters integrated in logic tree framework to capture epistemic uncertainties through three different seismic source models. It also makes use of the recently released version of OpenQuake engine; an open-source tool for seismic hazard and risk assessment developed in the framework of the Global Earthquake Model.  相似文献   

14.
Automated Detection, Extraction, and Measurement of Regional Surface Waves   总被引：3，自引：0，他引：3  
v--vOur goal is to develop and test an effective method to detect, identify, extract, and quantify surface wave signals for weak events observed at regional stations. We describe an automated surface wave detector and extractor designed to work on weak surface wave signals across Eurasia at intermediate periods (8 s-40 s). The method is based on phase-matched filters defined by the Rayleigh wave group travel-time predictions from the broadband group velocity maps presented by Ritzwoller and Levshin (1998) and Ritzwoller et al. (1998) and proceeds in three steps: Signal compression, signal extraction or cleaning, and measurement. First, the dispersed surface wave signals are compressed in time by applying an anti-dispersion or phase-matched filter defined from the group velocity maps. We refer to this as the `compressed signal.' Second, the surface wave is then extracted by filtering `noise' temporally isolated from the time-compressed signal. This filtered signal is then redispersed by applying the inverse of the phase-matched filter. Finally, we adaptively estimate spectral amplitude as well as group and phase velocity on the filtered signal. The method is naturally used as a detector by allowing origin time to slide along the time axis. We describe preliminary results of the application of this method to a set of nuclear explosions and earthquakes that occurred on or near the Chinese Lop Nor test site from 1992 through 1996 and one explosion on the Indian Rajasthan test site that occurred in May of 1998.  相似文献   

15.
Previous comparison studies on seismic isolation have demonstrated its beneficial and detrimental effects on the structural performance of high‐speed rail bridges during earthquakes. Striking a balance between these 2 competing effects requires proper tuning of the controlling design parameters in the design of the seismic isolation system. This results in a challenging problem for practical design in performance‐based engineering, particularly when the uncertainty in seismic loading needs to be explicitly accounted for. This problem can be tackled using a novel probabilistic performance‐based optimum seismic design (PPBOSD) framework, which has been previously proposed as an extension of the performance‐based earthquake engineering methodology. For this purpose, a parametric probabilistic demand hazard analysis is performed over a grid in the seismic isolator parameter space, using high‐throughput cloud‐computing resources, for a California high‐speed rail (CHSR) prototype bridge. The derived probabilistic structural demand hazard results conditional on a seismic hazard level and unconditional, i.e., accounting for all seismic hazard levels, are used to define 2 families of risk features, respectively. Various risk features are explored as functions of the key isolator parameters and are used to construct probabilistic objective and constraint functions in defining well‐posed optimization problems. These optimization problems are solved using a grid‐based, brute‐force approach as an application of the PPBOSD framework, seeking optimum seismic isolator parameters for the CHSR prototype bridge. This research shows the promising use of seismic isolation for CHSR bridges, as well as the potential of the versatile PPBOSD framework in solving probabilistic performance‐based real‐world design problems.  相似文献   

16.
ＴｈｅａｐｐｌｉｃａｔｉｏｎｏｆｓｅｉｓｍｉｃｄａｔａｗｉｔｈｄｉｆｆｅｒｅｎｔｐｒｅｃｉｓｉｏｎｉｎｔｈｅｄｅｔｅｒｍｉｎａｔｉｏｎｏｆｓｅｉｓｍｉｃｉｔｙｐａｒａｍｅｔｅｒｓＸｕｅ－ＳｈｅｎＪＩＮ（金学申）；Ｙｉｎｇ－ＨｕａＤＡＩ（戴英华）；Ｊｕ...  相似文献   

17.
This study presents a time-dependent approach for seismic hazard in Tehran and surrounding areas. Hazard is evaluated by combining background seismic activity, and larger earthquakes may emanate from fault segments. Using available historical and paleoseismological data or empirical relation, the recurrence time and maximum magnitude of characteristic earthquakes for the major faults have been explored. The Brownian passage time (BPT) distribution has been used to calculate equivalent fictitious seismicity rate for major faults in the region. To include ground motion uncertainty, a logic tree and five ground motion prediction equations have been selected based on their applicability in the region. Finally, hazard maps have been presented.  相似文献   

18.
A vital component of any seismic hazard analysis is a model for predicting the expected distribution of ground motions at a site due to possible earthquake scenarios. The limited nature of the datasets from which such models are derived gives rise to epistemic uncertainty in both the median estimates and the associated aleatory variability of these predictive equations. In order to capture this epistemic uncertainty in a seismic hazard analysis, more than one ground-motion prediction equation must be used, and the tool that is currently employed to combine multiple models is the logic tree. Candidate ground-motion models for a logic tree should be selected in order to obtain the smallest possible suite of equations that can capture the expected range of possible ground motions in the target region. This is achieved by starting from a comprehensive list of available equations and then applying criteria for rejecting those considered inappropriate in terms of quality, derivation or applicability. Once the final list of candidate models is established, adjustments must be applied to achieve parameter compatibility. Additional adjustments can also be applied to remove the effect of systematic differences between host and target regions. These procedures are applied to select and adjust ground-motion models for the analysis of seismic hazard at rock sites in West Central Europe. This region is chosen for illustrative purposes particularly because it highlights the issue of using ground-motion models derived from small magnitude earthquakes in the analysis of hazard due to much larger events. Some of the pitfalls of extrapolating ground-motion models from small to large magnitude earthquakes in low seismicity regions are discussed for the selected target region.  相似文献   

19.
This paper presents a new hazard-consistent ground motion characterization of the Itoiz dam site, located in Northern Spain. Firstly, we propose a methodology with different approximation levels to the expected ground motion at the dam site. Secondly, we apply this methodology taking into account the particular characteristics of the site and of the dam. Hazard calculations were performed following the Probabilistic Seismic Hazard Assessment method using a logic tree, which accounts for different seismic source zonings and different ground-motion attenuation relationships. The study was done in terms of peak ground acceleration and several spectral accelerations of periods coinciding with the fundamental vibration periods of the dam. In order to estimate these ground motions we consider two different dam conditions: when the dam is empty (T?=?0.1?s) and when it is filled with water to its maximum capacity (T?=?0.22?s). Additionally, seismic hazard analysis is done for two return periods: 975?years, related to the project earthquake, and 4,975?years, identified with an extreme event. Soil conditions were also taken into account at the site of the dam. Through the proposed methodology we deal with different forms of characterizing ground motion at the study site. In a first step, we obtain the uniform hazard response spectra for the two return periods. In a second step, a disaggregation analysis is done in order to obtain the controlling earthquakes that can affect the dam. Subsequently, we characterize the ground motion at the dam site in terms of specific response spectra for target motions defined by the expected values SA (T) of T?=?0.1 and 0.22?s for the return periods of 975 and 4,975?years, respectively. Finally, synthetic acceleration time histories for earthquake events matching the controlling parameters are generated using the discrete wave-number method and subsequently analyzed. Because of the short relative distances between the controlling earthquakes and the dam site we considered finite sources in these computations. We conclude that directivity effects should be taken into account as an important variable in this kind of studies for ground motion characteristics.  相似文献   

20.
Semi‐active variable stiffness resettable devices can reduce seismic demands and damages in structures. Despite their advantages, variable stiffness resettable devices are under‐utilized mainly because of the shortage of fundamental research in quantifying the sensitivity of key seismic response parameters, and losses, in structures that use such systems for seismic hazard mitigation. Within this setting, the research summarized herein measures the effectiveness of semi‐active resettable energy dissipating devices in the Single‐Degree‐of‐Freedom domain aiming at quantifying the sensitivity of their seismic response to variation in control parameters and generating the required knowledge to utilize such semi‐active devices in the Multi‐Degree‐of‐Freedom domain. The performance (i.e. maximum relative displacement and peak absolute acceleration demands) of Single‐Degree‐of‐Freedom systems with an array of semi‐active control logics under various dynamic excitation regimes is studied. Two sets of 40 ground motions representing various seismic loading conditions (i.e. pulse‐like and rock‐site ground motions) are used, and an efficient control logic for mitigating these seismic demands is proposed. Numerical results show that proposed control logic enables a decrease of 40–60% for both maximum relative displacement and seismic base shear and 15–25% decrease for peak absolute acceleration. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

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

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