首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
汶川-映秀MS8.0地震的介质破裂与深部物质运移的动力机制   总被引:1,自引:0,他引:1  
2008年5月22日14时28分在龙门山造山带发生了震惊世界的汶川—映秀MS8.0大地震。这次大地震形成了迄今为止在三维空间分布最为复杂、在水平方向(NE)延伸最长的逆冲型破裂带。笔者分析和讨论了大地震发生后在地表引发的破裂现象和震源深处与其周边地域介质在受力作用下的破裂与"破裂链"的逐步形成和辐射效应。基于该大地震"孕育"、发生和发展的深部介质和构造环境与深层动力过程的研究提出:汶川—映秀8.0级地震的发震断层是龙门山断裂系地表以不同角度西倾的3条断裂带向深部汇聚,青藏高原东北缘松潘—甘孜地域下地壳与上地幔盖层物质向东南减薄,在受到四川盆地"刚性"壳、幔介质阻隔下而沿龙门山断裂系抬升,二者在震源深处,即(15±5)km处强烈碰撞而形成的一条NE向的深部汇聚断裂带。汇聚断裂带即为MS8.0地震和一系列余震的发震断裂带,而不是简单地、形式地将地表的某一条或某两条断裂带视为发震断裂带,因为震源是一个体积。大地震发生后的地表破裂、次生灾害等,即浅层过程的调查和分析对恢复生产、重建家园有着重要意义。显然,在地表和深井中若能长期进行破裂效应的观测,以捕捉地震强烈活动地区震源深处介质与结构发生的初始微破裂和其形成"破裂链...  相似文献   

2.
龙门山断裂带深部构造和物性分布的分段特征   总被引:8,自引:0,他引:8  
根据龙门山断裂带周边的固定数字地震台网和流动地震观测获得的宽频带地震记录,用多种地震学方法研究该地区的地壳上地幔结构。深部结构研究表明,龙门山断裂带物性分布具有显著的分段特征。用远震接收函数H-k叠加方法计算了各个台站的地壳厚度和波速比。地壳厚度总体变化是,地壳从东向西增厚,最小厚度为37.8 km,最大厚度是68.1 km。从东南向西北横跨龙门山断裂带的地壳急剧增厚,从41.5 km增厚至52.5 km。但是,龙门山断裂带两侧地壳厚度的差异在断裂带的南段和北段是不同的。在南段,地壳厚度急剧变化的分界线在中央断裂附近;在中段,分界线在后山断裂附近;在北段,则断裂带两侧地壳厚度差异很小。泊松比的空间分布是,松潘—甘孜地体北部和西秦岭造山带具有低泊松比(ν<0.26),扬子地台具有低—中泊松比(ν<0.27),松潘—甘孜地体南部、三江褶皱带和四川盆地具有中—高泊松比(0.26<ν<0.29)。除龙门山断裂带南段及其附近,大部分地区均不具有超高的泊松比(ν>0.30)。龙门山断裂带南段地壳具有高泊松比(ν>0.30),而北段地壳则为中—低泊松比。高泊松比可以看成是铁镁质组分增加和/或部分熔融的证据,表明那里的下地壳部分熔融是可能的。松潘—甘孜地体东南部地区的下地壳处于富含流体或温度较高的部分熔融状态,它有助于青藏高原的下地壳物质向东运动。青藏高原东部中、上地壳向东运动受刚性强度较大的扬子地台的阻挡,沿龙门山断裂带产生应变能积累。当应变达到临界值,发生急剧的摩擦滑动,释放积累的应变能,产生汶川Ms8.0地震。汶川地震在龙门山断裂带不同地段,表现出不同的破裂特征和余震分布,可能与断层带的分段深部构造差异有关。  相似文献   

3.
汶川8.0级特大地震震源断裂特征及其动力学分析   总被引:1,自引:1,他引:0       下载免费PDF全文
笔者根据地壳表层地质构造和地震地质研究与地震测深和大地电磁测深成果,运用现代构造解析理论与方法,论证了汶川8.0级特大地震的深部构造环境,探讨了汶川8.0级特大地震震源区的地震断裂和震源断裂基本特征与相互关系,及其形成的地球动力学问题。笔者认为龙门山碰撞造山带深处隐伏壳幔韧性剪切带可能是汶川8.0级特大地震主震区的震源断裂,而地壳表层发育的映秀断裂带、北川断裂带和彭县—灌县断裂带等可能是汶川8.0级特大地震主震区的地震断裂,该区震源断裂与地震断裂既有显著区别,又有密切联系。研究表明,在印度板块与太平洋板块和菲律宾海板块对欧亚板块俯冲碰撞的动力学作用下,形成上扬子地块向青藏高原东缘碰撞—楔入以及青藏高原东缘向东仰冲,深部向东俯冲的动力学态势,造成龙门山碰撞造山带切割莫霍界面的壳幔韧性剪切带向中上地壳扩展,应力高度集中与能量快速释放破裂,从而引起汶川8.0级特大地震的发生,以及大型同震地表破裂带的形成。探索震源断裂与地震断裂区别与联系,对进一步研究地震机制与发震动力学以及防震减灾有重要意义。  相似文献   

4.
2008年5月12日在青藏高原东缘龙门山断裂带中段发生汶川8.0级特大地震。大震发生时释放应力并对震源区及外围构造应力场产生影响,受汶川地震断层破裂方式和强度空间差异性的影响,震后龙门山断裂带地壳应力场也应表现差异特征,至今鲜有针对该科学问题深入的分析和讨论。经过系统收集、梳理汶川地震后沿龙门山断裂带水压致裂地应力测量数据与2008年汶川地震中强余震序列震源机制解资料,对汶川地震后龙门山断裂带中上地壳构造应力场进行厘定,通过与震前构造应力场对比,深入探讨了汶川8.0级地震对龙门山断裂带地壳应力场的影响,进而对汶川震后应力调整过程及青藏高原东缘龙门山地区深部构造变形模式进行研究,研究结果表明:受汶川8.0级地震的影响,震后龙门山断裂带地壳构造应力场空间分布具有差异性,近地表至上地壳15 km深度范围,映秀—青川段最大主应力方向为北西西向、地应力状态为逆走滑型,青川东北部最大主应力方向偏转至北东东向、应力状态转变为走滑型;15~25km深度范围,龙门山断裂带最大主应力方向仍为北西—北西西向、应力状态以逆冲型为主。汶川8.0级地震后,龙门山断裂带中地壳北西西向逆冲挤压的构造应力特征进一步支持了青藏高原东缘龙门山地区东西两侧刚性块体碰撞挤压、逆冲推覆的动力学模式。  相似文献   

5.
青藏高原东缘地壳上地幔结构及其动力学意义   总被引:3,自引:0,他引:3       下载免费PDF全文
张忠杰  陈赘  田小波 《地质科学》2009,44(4):1136-1150
本文综述了我们在青藏高原东缘实施的垂直切过龙门山断裂带宽频带地震探测的研究成果,揭示了研究区复杂的地壳上地幔结构,结果表明松潘-甘孜地块与四川盆地西缘莫霍面深度为58 km与40 km±,在龙门山断裂带下方存在约15 km的莫霍面错断; 松潘-甘孜与龙门山断裂带域地壳纵横波速度比Vp/Vs比值远大于173,预示着粘性下地壳流或基性/超基性物质的存在。探讨了研究区强烈的盆山之间以及深部不同层圈之间的相互作用,推断四川盆地对青藏高原东缘软流圈驱动的物质东向逃逸阻挡作用可能深达整个上地幔。  相似文献   

6.
2008年5月12日在低地形变速率的龙门山断裂带上突发汶川强震,引发人们对该地震孕震机制的思考.本文根据GPS观测资料确定边界条件,通过三维粘弹性数值模拟探讨了汶川地震的孕震机理,计算了该区域岩石圈的应力增加速率和积累过程,以及汶川地震同震应力变化与震后应力松弛,在此基础上估算了汶川8.0级大地震的复发周期.数值模拟结果表明:印度板块对欧亚板块的推挤造成青藏高原的物质东流,高原中、下地壳物质在龙门山断裂带处遭到相对坚硬的四川盆地的阻挡之后,部分中、下地壳物质在龙门山断裂带下堆积产生应力集中.两个重要因素为应力集中提供了重要控制作用:其一是青藏高原中、下地壳较低的粘滞系数与四川盆地中、下地壳较高的粘滞系数的差异,其二是从青藏高原到四川盆地的Moho面深度在龙门山断裂带的突变.低应变速率的龙门断裂带岩石圈在数千年时间尺度的应力积累过程中,脆性上地壳的应力随时间近乎线性增长,并且上地壳深部的应力增长率超过浅部,6000年内应力积累最大量达到-21.6MPa,应力增长速率为-0.0036MPa/a;而柔性的中、下地壳以及岩石圈上地幔的应力在增长一段时间之后趋于稳定.在空间上,龙门山断裂带受到的压应力从断层西南向北东方向逐渐减小,而剪应力从西南到北东方向逐渐增大,应力状态有利于地震发生时断层的破裂方式从西南的逆冲运动向北东的逆冲兼走滑运动的方式发展.通过应力积累与地震应力降的计算得到汶川8.0级大地震的复发周期约为5400年.  相似文献   

7.
粘弹性数值模拟龙门山断裂带应力积累及大震复发周期   总被引:4,自引:0,他引:4  
柳畅  朱伯靖  石耀霖 《地质学报》2012,86(1):157-169
2008年5月12日在低地形变速率的龙门山断裂带上突发汶川强震,引发人们对该地震孕震机制的思考。本文根据GPS观测资料确定边界条件,通过三维粘弹性数值模拟探讨了汶川地震的孕震机理,计算了该区域岩石圈的应力增加速率和积累过程,以及汶川地震同震应力变化与震后应力松弛,在此基础上估算了汶川8.0级大地震的复发周期。数值模拟结果表明:印度板块对欧亚板块的推挤造成青藏高原的物质东流,高原中、下地壳物质在龙门山断裂带处遭到相对坚硬的四川盆地的阻挡之后,部分中、下地壳物质在龙门山断裂带下堆积产生应力集中。两个重要因素为应力集中提供了重要控制作用:其一是青藏高原中、下地壳较低的粘滞系数与四川盆地中、下地壳较高的粘滞系数的差异,其二是从青藏高原到四川盆地的Moho面深度在龙门山断裂带的突变。低应变速率的龙门断裂带岩石圈在数千年时间尺度的应力积累过程中,脆性上地壳的应力随时间近乎线性增长,并且上地壳深部的应力增长率超过浅部,6000年内应力积累最大量达到-21.6MPa,应力增长速率为-0.0036MPa/a;而柔性的中、下地壳以及岩石圈上地幔的应力在增长一段时间之后趋于稳定。在空间上,龙门山断裂带受到的压应力从断层西南向北东方向逐渐减小,而剪应力从西南到北东方向逐渐增大,应力状态有利于地震发生时断层的破裂方式从西南的逆冲运动向北东的逆冲兼走滑运动的方式发展。通过应力积累与地震应力降的计算得到汶川8.0级大地震的复发周期约为5400年。  相似文献   

8.
王芃  张忠杰  张晰  韩颜颜  王敏玲  侯爵  徐涛 《岩石学报》2014,30(4):1179-1187
龙门山是我国东西构造、地貌分界线的重要组成部分。其两侧的岩石圈结构差异,是形成龙门山造山带的主要原因之一,并对龙门山的构造演化起着持续影响。为了解龙门山两侧壳幔结构差异,本文从重力角度探讨跨龙门山地区的地壳密度结构。我们使用EGM2008模型的重力异常数据,以最新的阿坝-遂宁人工源地震剖面速度模型为基础,得到了龙门山造山带中段及其邻区的精细地壳密度结构。密度结构显示松潘-甘孜地区和四川盆地分别具有软弱和坚硬的下地壳。根据本文所得到的地壳密度结构模型,我们认为龙门山的隆升主要受印度洋板块与欧亚大陆板块的陆-陆碰撞作用影响,强烈的挤压作用使青藏高原物质向东运移,东移物质在青藏高原东缘龙门山地区受到坚硬的四川盆地的阻挡转而向上运移,造成了龙门山的隆升。  相似文献   

9.
汶川地震断裂带科学钻探1号井(WFSD-1)的ASR三维地应力测试结果表明,龙门山前陆逆冲带与其下伏的龙门山前陆盆地和上覆的松潘-甘孜地块的构造及地应力状态存在有重大差异。从整体上看,在汶川地震中,龙门山前陆逆冲带表现为在强烈的区域性挤压背景下,深部物质沿壳内拆离层自SW向NE方向的"层状"流动,在地壳上部转化为沿映秀-北川断裂(YBF)的快速垂向挤出,而其西侧的松潘-甘孜地块作自SE往NW方向的重力滑覆,东侧的龙门山前陆盆地则表现为自NE往SW方向的走滑或右行旋转。晚新生代以来,扬子地块相对于青藏高原东缘的龙门山造山带并无明显的或大尺度的陆内俯冲作用发生。龙门山前陆逆冲带深部高温低粘度物质垂直向上的、快速的流动和挤出,直接导致了"5.12"汶川地震的发生,而松潘-甘孜地块E向扩展导致龙门山前陆带的强烈挤压和陆壳增厚及深部应力和地震能量的积聚则是诱导深部位移场发生突变和物质快速垂向挤出的主因,E向扩展是深部地震能量积聚和快速垂向挤出作用的必要条件,而非地震发生的直接原因。ASR地应力测试得出的主压应力方向完全平行于GPS同震速度场的位移方向,似乎表明ASR测试获得的原地应力场或许真实地反映了或最接近于地震过程中的构造应力状态。  相似文献   

10.
龙门山断裂带是中国大陆地壳中著名的造山带和地震带。为了研究龙门山逆冲推覆构造的深部结构及其与青藏高原物质东向逃逸之间的关系,布置了跨越西秦岭造山带、松潘—甘孜褶皱带、四川盆地的4条大地电磁测深剖面,每条剖面有宽频大地电磁测深点20个,共计80个宽频大地电磁测深数据以及16个长周期数据;进行了二维、三维反演,得到了电性模型。研究结果显示:在30 km以浅,龙门山断裂带基本上与电性梯度变化带相重合,说明龙门山断裂带是深达地壳中部的深大断裂带;松潘-—甘孜褶皱带在20 km以浅呈现以高阻为主的较为复杂的电性结构分布特征,这与其复杂的地表构造有关;在20 km以深,接近龙门山断裂带附近呈现为高阻,推测此高阻可能是对四川盆地的基底的反映,表明松潘—甘孜褶皱带以龙门山断裂为界推覆至四川盆地之上;推断青藏高原物质东向逃逸,有可能从松潘—甘孜褶皱带深部沿礼县—宕昌一带向东北方向逃逸,这一运移主要是发生在18 km以下、30 km以上的中下地壳范围内,物质逃逸的主要形式是部分熔融。  相似文献   

11.
Komatiites are mantle-derived ultramafic volcanic rocks. Komatiites have been discovered in several States of India, notably in Karnataka. Studies on the distribution of trace-elements in the komatiites of India are very few. This paper proposes a simple, accurate, precise, rapid, and non-destructive wavelength-dispersive x-ray fluorescence (WDXRF) spectrometric technique for determining Sc, V, Cr, Co, Ni, Cu, Zn, Rb, Sr, Y, Zr, Nb, Ba, Pb, and Th in komatiites, and discusses the accuracy, precision, limits of detection, x-ray spectral-line interferences, inter-element effects, speed, advantages, and limitations of the technique. The accuracy of the technique is excellent (within 3%) for Sc, V, Cr, Co, Ni, Cu, Zn, Rb, Sr, Zr, Nb, Ba, Pb, and Th and very good (within 4%) for Y. The precision is also excellent (within 3%) for Sc, V, Cr, Co, Ni, Cu, Zn, Rb, Sr, Y, Zr, Nb, Ba, Pb, and Th. The limits of detection are: 1 ppm for Sc and V; 2 ppm for Cr, Co, and Ni; 3 ppm for Cu, Zn, Rb, and Sr; 4 ppm for Y and Zr; 6 ppm for Nb; 10 ppm for Ba; 13 ppm for Pb; and 14 ppm for Th. The time taken for determining Sc, V, Cr, Co, Ni, Cu, Zn, Rb, Sr, Y, Zr, Nb, Ba, Pb, and Th in a batch of 24 samples of komatiites, for a replication of four analyses per sample, by one operator, using a manual WDXRF spectrometer, is only 60 hours.  相似文献   

12.
最新的流行病学研究表明,空气中较高浓度的悬浮细颗粒可能对人类的健康有不利的影响。根据该项研究显示,由于心脏病、慢性呼吸问题和肺功能指标恶化而导致死亡率的升高与细尘粒子有关。这些研究结果已经促使欧盟于1999年4月出台了限制空气中二氧化硫、二氧化氮、氧化氮、铅和颗粒物含量的法案(1999/30/EC),对各项指标包括对可吸入PM10颗粒的浓度提出了新的限制性指标。PM10颗粒是指可以通过预分级器分离采集的气体动力学直径小于10μm的细颗粒。目前研究的兴趣重点逐步偏向PM2.5这些更细微颗粒物,PM2.5这种颗粒物对健康有明显的不利影响。在欧盟指令2008/50/EC中,对PM10和PM2.5都提  相似文献   

13.
This paper reports the first results of a study of 11 isotope systems (3He/4He, 40Ar/36Ar, 34S/32S, 65Cu/63Cu, 62Ni/60Ni, 87Sr/86Sr, 143Nd/144Nd, 206–208Pb/204Pb, Hf–Nd, U–Pb, and Re–Os) in the rocks and ores of the Cu–Ni–PGE deposits of the Norilsk ore district. Almost all the results were obtained at the Center of Isotopic Research of the Karpinskii All-Russia Research Institute of Geology. The use of a number of independent genetic isotopic signatures and comprehensive isotopic knowledge provided a methodic basis for the interpretation of approximately 5000 isotopic analyses of various elements. The presence of materials from two sources, crust and mantle, was detected in the composition of the rocks and ores. The contribution of the crustal source is especially significant in the paleofluids (gas–liquid microinclusions) of the ore-forming medium. Crustal solutions were probably a transport medium during ore formation. Air argon is dominant in the ores, which indicates a connection between the paleofluids and the atmosphere. This suggests intense groundwater circulation during the crystallization of ore minerals. The age of the rocks and ores of the Norilsk deposits was determined. The stage of orebody formation is restricted to a narrow age interval of 250 ± 10 Ma. An isotopic criterion was proposed for the ore-bearing potential of mafic intrusions in the Norilsk–Taimyr region. It includes several interrelated isotopic ratios of various elements: He, Ar, S, and others.  相似文献   

14.
本文介绍样品经四硼酸锂熔融制成玻璃小饼。采用Lachance模式和理论a系数来校正元素间的效应,由3080E型X-射线荧光谱仪和DF-350B数据处理系统完成硅酸盐中十三个项目的测定。  相似文献   

15.
本文拟定了一种以熔融法制备样片,用X射线荧光光谱测定硅酸盐类样品中Si、Fe、Al、Tj、Mn、Ca、Mg、K、Na、P等元素的分析方法。在对不同靶材X光管和分光晶体实验对比的基础上,选择了最佳的测量条件。该法具有快速、准确,测量范围广,检测限低,价格便宜等优点。经过近百个各种类型标样或管理样品的分析对比表明,本法不仅适用于硅酸盐类岩石样品的分析,还适用于铁矿、铝土矿、碳酸盐类岩矿样品以及水泥、耐火材料等样品的分析。  相似文献   

16.
Mining induced subsidence can significantly affect mining costs where major surface facilities and natural environment need to be protected. Overburden grout injection is a technology used to control coal mine subsidence by injecting the mine waste material extracted from the coal back into the inter-burden rock during longwall mining. The flowing slurry is here categorised as a nonlinear viscous cohesive (Bingham plastic) fluid. During longwall mining the grout slurry is pumped into the separated beds of the rock mass through a central vertical borehole, which is drilled deep into the inter-burden rock strata above the coal seam. However, a blockage can occur in the injection system when the slurry velocity falls below a certain critical threshold velocity, indicating a material phase change from cohesive-viscous to cohesive-frictional. In situ field injection tests through boreholes have been simulated at a smaller scale at the CSIRO laboratory in Brisbane by pumping the slurry through a radial disk (gap = 4 mm) from its centre. Laboratory experiments indicate a general, nonlinear, cohesive, viscous, frictional model for shear behaviour of the slurry, in which the material shear parameters are functions of the disk radial distance. Complete dimensional and dimensionless analytical solutions have been developed based on an approach related to Bingham–Herschel–Bulkley fluid mechanics. The derived formulae include relations for minimum pump pressure, local pressure and pressure gradient, wall shear stress, volume rate, velocity and velocity gradient. The theoretical results match the experimental measurements. The experiments covered slurries with maximum particle sizes of 0.5 to 2 mm with about 50% being larger than 100 µm. The viscosities at the various solids concentrations were measured with a standard torsion viscometer. This study differs from the previous research in several distinct aspects, namely, consideration of the variable shear parameters rather than fixed values, inclusion of total nonlinear behaviour, and implementation of a friction function to mimic behaviour of the deposited and consolidating stiff slurry, which can cause a significant pressure rise as a result of the increased shear resistance.  相似文献   

17.
18.
岩石密度和超高压岩石折返速率   总被引:1,自引:3,他引:1  
在常温常压条件下对中国大陆科学钻CCSD主孔100-3000米岩心样品进行了密度测量,建立了密度连续剖面,并界定了不同超高压岩石的密度值。通过对比高温高压物性实验资料,岩石密度随着退变程度增强而降低,榴辉岩密度变化序列为3.52g/cm3→3.49g/cm3→3.07g/cm3→2.93g/cm3。超高压长英质岩石密度变化序列为3.00g/cm→2.80g/cm3→2.65g/cm3。上述实验资料是讨论不同折返阶段岩石所受浮力的基础,为研究折返速率大小提供了基本参数。本文通过折返板块运动平衡时,上浮力与粘滞力平衡这一关系式,定量研究了大陆俯冲板块的折返速率,认为密度差产生上浮力从而引起折返,温度对板块折返速率的影响最为显著;密度差大小、折返角度、折返板块大小对折返速率也有直接的影响。定量模拟分析表明,在温度高于850℃时,板块的折返速率可以超过100mm/a;当温度降至700℃时,折返速率则低于1.5mm/a。作者认为在折返早期,温度较高,板块快速折返至60-70km榴辉岩相深度;随着传导散热,温度降低,板块以较慢的速率折返至中下地壳。折返速率的估算表明,浮力是板块折返第一阶段(从>100km深部折返至<40km的中下地壳)的主要驱动力。  相似文献   

19.
《Applied Geochemistry》2001,16(2):137-159
Five hundred and ninety-eight samples of terrestrial moss (Hylocomium splendens and Pleurozium schreberi) collected from a 188,000 km2 area of the central Barents region (NE Norway, N Finland, NW Russia) were analysed by ICP-AES and ICP-MS. Analytical results for Al, B, Ba, Ca, K, La, Mg, Mn, Na, P, Rb, Si, Sr, Th, U and Y concentrations are reported here. Graphical methods of data analysis, such as geochemical maps, cumulative frequency diagrams, boxplots and scatterplots, are used to interpret the origin of the patterns for these elements. None of the elements reported here are emitted in significant amounts from the smelting industry on the Kola Peninsula. Despite the conventional view that moss chemistry reflects atmospheric element input, the nature of the underlying mineral substrate (regolith or bedrock) is found to have a considerable influence on moss composition for several elements. This influence of the chemistry of the mineral substrate can take place in a variety of ways. (1) It can be completely natural, reflecting the ability of higher plants to take up elements from deep soil horizons and shed them with litterfall onto the surface. (2) It can result from naturally increased soil dust input where vegetation is scarce due to harsh climatic conditions for instance. Alternatively, substrate influence can be enhanced by human activity, such as open-cast mining, creation of ‘technogenic deserts’, or handling, transport and storage of ore and ore products, all of which magnify the natural elemental flux from bedrock to ground vegetation. Seaspray is another natural process affecting moss composition in the area (Mg, Na), and this is most visible in the Norwegian part of the study area. Presence or absence of some plant species, e.g., lichens, seems to influence moss chemistry. This is shown by the low concentrations of B or K in moss on the Finnish and Norwegian side of the (fenced) border with Russia, contrasting with high concentrations on the other side (intensive reindeer husbandry west of the border has selectively depleted the lichen population).  相似文献   

20.
We have developed a set of four synthetic standards for the rare earth and high field strength elements designed for use in the determination of those elements in silicates. The base material is a glass nominally at the eutectic of the MgO-Al2O3-SiO2 ternary system. The nominal doping level was 2 wt% of the elements as oxides. To avoid problems associated with peak interference, the elements were doped in four separate glasses.  相似文献   

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

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