首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 187 毫秒
1.
通过断裂超微变形构造分析,对柴达木盆地周边断裂构造发育的应力-应变环境进行了较为深入的探讨。研究结果表明:新高泉断裂超微构造以密度较低的平直位错线为典型特征,反映了相对较低的应力-应变环境;大头羊逆冲断裂位错密度高,位错形态复杂,呈弯曲状延伸,局部发育位错环,应力-应变环境较高;绿梁山逆冲断裂的超微构造以亚颗粒的普遍发育为显著特征,表明应变环境已达恢复阶段;纳赤台逆冲断裂的超微构造较为复杂,位错线、位错壁和亚颗粒都有不同程度的发育,表明变形的复杂性。根据超微构造的特征,这对不同断裂的古构造应力和应变速率进行了计算。   相似文献   

2.
平南幔源包体中橄榄石的显微构造研究及其意义   总被引:1,自引:0,他引:1  
平南玄武岩中的尖晶石二辉橄榄岩包体的平衡温度为930~980℃, 平衡深度为59~74km, 包体中橄榄石的扭折带滑移系多为(010) [100], 但也有(001) [100]的滑移系类型; 斜方辉石的滑移系为(100) [001], 它们均为高温低应变速率下的滑移系, 说明该区的上地幔主要是在高温低应变速率条件下经历了塑性变形作用.橄榄石位错组态多样, 有自由位错、位错壁、位错弓弯、缠结、{110}滑移带, 反映了上地幔的塑性变形特征.根据位错壁的大小估算, 上地幔差异流动应力为24.5~42.1MPa, 流动速率为2.93×10-17~8.36×10-16s-1, 有效粘度为1.72×1023~2.80×1024 Pa·s, 特征与中国东部新生代上地幔较为一致, 反映同处于拉张环境.   相似文献   

3.
头道桥韧性剪切带是研究区一条较典型的、发育比较完整的韧性剪切带,形成糜棱岩、碎裂岩系列.镜下观察区内糜棱岩具多种显微变形结构,主要表现为粒间位移、晶内变形、恢复作用、重结晶作用,从变形环境分析该韧性剪切带具分带性.其变形温度在350~500℃,压力在300~500 MPa,发育最大深度约为22 km.估算差异应力和应变速率分别在39.96~96.48 MPa和(2.17~23.60)×10-12 s-1之间,表明该区糜棱岩是缓慢变形的结果.  相似文献   

4.
劝农山地区位于长春市东南部,处于佳-伊断裂和西拉木伦河缝合带交汇处.详细野外调查发现,该区曾遭受强烈韧性剪切变形,剪切带内岩石普遍糜棱岩化,主要由下二叠统范家屯组(P1f)钙质糜棱岩与侵入其中的燕山期花岗质糜棱岩组成,变形程度处于初糜棱岩至糜棱岩之间,多具有糜棱结构.岩石应变类型主要为压扁型应变,偏一般压缩,为L=S型构造岩,指示其形成于挤压型剪切带的构造环境.多种宏微观韧性剪切变形标志,指示明显的左行剪切运动.电子探针方解石-白云石地质温度计、方解石和石英EBSD组构特征、方解石e双晶形态以及石英长石变形行为等均显示岩石具有低温塑性流变特点,变形环境不超过绿片岩相.剪切带内应变速率偏高,应变集中带应变速率最大,在10-6.95~10-8.89之间,远离强变形带应变速率在10-9.25~10-12.17之间,糜棱岩化作用过程中差异应力下限应大致为51.27~65.46 MPa,代表剪切带糜棱岩化作用为低温中等强度应变,在稍快的应变速率条件下形成.压溶扩散和双晶滑移为劝农山韧性剪切带变形初期的主要变形机制,随着递进变形,逐渐以双晶滑移和晶内滑移为主,递进变形晚期,局部强变形域内发生了粒间滑移.劝农山韧性剪切带形成与早白垩世中晚期伊泽纳崎板块NNW向高斜度斜向俯冲于欧亚大陆之下有关,是佳-伊断裂带左旋走滑事件的局部表现.   相似文献   

5.
开展断层脆塑性转化带的变形机制、断层带流体-岩石相互作用、断层愈合作用等研究,对理解间震期、同震加载、震后滑动阶段断层的变形机制转化、强震孕育和发生具有重要意义。笔者采用Carrara大理岩,在温度300~700℃、围压300 MPa和600~800 MPa、应变速率1×10-4/s~1×10-5/s和1×10-7/s~2.5×10-6/s、水含量0.005%~0.01%和0.1%~0.5%条件下,开展了轴向压缩变形实验与裂隙愈合实验。通过偏光显微镜、扫描电镜与能谱分析,研究了实验变形样品的微观结构与变形机制,讨论了水、温度、围压、应变速率对脆塑性转化和变形机制的影响,以及裂隙愈合对断层强度和流体压力变化的制约作用。实验结果表明:(1)Carrara大理岩在低温(300~400℃)、低应变速率和高含水条件下发生了压溶,其中,在低温低应变速率(1×10-7/s)条件下为压溶蠕变,在低温中等应变速率(5×10-7/s)条件下为压溶+碎裂流动。(2)在低温(400℃)、中-高应变速率和低含水条件下发生了位错滑移(双晶滑移、机械双晶)与碎裂流动,局部伴有压溶作用。(3)在中温(500℃)、各应变速率和各含水条件下发生了位错滑移(双晶滑移、机械双晶)与动态重结晶作用。(4)在高温(600~700℃)条件下,动态重结晶作用成为主要变形机制。(5)在压溶和动态重结晶作用下,在脆性变形阶段产生的裂隙与孔隙被愈合。断层强度恢复程度受裂隙和孔隙愈合程度控制。温度、水和应变速率对大理岩脆塑性转化和变形机制的影响非常显著,在相同温度与应变速率条件下,水降低了样品强度,促进了压溶和塑性变形。增加应变和水含量,能够显著促进裂缝和孔隙愈合。根据实验结果推测:在快速变形的同震和震后滑动阶段,断层脆塑性转化带以碎裂变形为主;在缓慢变形的间震期,断层脆塑性转化带以压溶和动态重结晶为主。在塑性变形作用下,同震滑动产生的裂隙被愈合,不仅恢复了断层带强度,而且为断层带内部形成高压流体创造了条件。  相似文献   

6.
桐柏杂岩北界剪切带的构造变形特征研究   总被引:1,自引:0,他引:1  
桐柏杂岩北界剪切带是桐柏—大别造山带发育的一条大型的左行走滑韧性剪切带,发育糜棱岩、超糜棱岩以及强直片麻岩等强变形岩石。通过对剪切带内岩石的矿物共生组合、剪切带的位移量计算、石英的分维数计算、石英颗粒的差异应力以及应变速率的估算等多种方法分析结果显示:桐柏杂岩北界剪切带的位移量为38.3~41.1km;动态重结晶石英颗粒边界统计数据的相关系数R均大于0.96,动态重结晶颗粒边界具有统计意义上的自相似性;分维数为1.207~1.302之间,从剪切带的核部向边部,颗粒粒径和分维值的变化是逐渐变大的。石英的分维值反映出石英的动态重结晶与剪切带密切相关,变质变形由剪切带核部向边部逐渐减弱。桐柏杂岩北界剪切带的差异应力值为23.583~46.983 MPa之间,应变速率为1.70×10~(-11)~9.74×10~(-11)s~(-1)。通过对比剪切带内动态重结晶石英颗粒的粒径、差异应力和应变速率可以发现,差异应力与应变速率表现出一定的正相关性,都随着从剪切带的核部到边部慢慢变小;而石英的粒径则与差异应力和应变速率表现出负相关性,从剪切带核部到边部,随着石英粒径的变大,差异应力和应变速率变小。差异应力越大,导致石英颗粒的受力强度越大,变形越强,表现出来的就是应变速率越强。从显微构造变形角度来研究认识该剪切带,对揭示桐柏—大别造山带的构造演化以及构造模型的建立有着重要的构造意义。  相似文献   

7.
李瑞红  刘育  李海林  郑小礼  赵海  孙政 《岩石学报》2014,30(9):2546-2558
新城金矿床是典型的"焦家式"破碎带蚀变岩型金矿,矿体形态和规模都严格受到断裂破碎带控制,是探讨复杂构造-流体耦合成矿系统控矿构造变形环境研究的理想选区。断裂破碎带中构造岩既是构造变形行为的载体,也是相应变形环境的受体。论文在新城金矿详细露头构造解析的基础上,系统采集该矿床控矿断裂破碎带定向构造岩样品,进行显微构造和EBSD组构分析。研究区构造岩显微构造特征主要表现为韧性变形和脆性变形。韧性变形有波状消光、带状消光、亚晶粒、动态重结晶、核幔构造、丝带构造、碎(残)斑系、扭折带、变形纹、机械双晶、蠕英结构、云母鱼等;脆性变形有书斜构造和显微裂隙等。长石(残)斑系、扭折带、变形纹、蠕英结构和石英颗粒边界迁移动态重结晶、丝带构造等矿物变形特征表明断裂带成矿前以高温韧性变形为主;石英波状消光、亚晶粒、亚颗粒旋转和膨凸动态重结晶、方解石机械双晶、长石显微裂隙充填物等矿物变形反映成矿期兼有中低温韧性变形和脆性变形;压剪性穿晶裂隙则反映出成矿后主要是低温脆性变形。根据差应力、应变测量和EBSD组构分析,将新城金矿床控矿构造变形环境可以分为3个构造期:成矿前在NW-SE向挤压作用下发生韧-脆性左行剪切变形,600~700℃,差应力61.37~111.09MPa,应变测量轴比a/c为2.295~3.978,动态重结晶石英颗粒边界分维值为1.466~1.599,反映矿区为高温中高压高应变带变形环境,应变速率较大;成矿期为NW-SE向逐渐NEE-SWW向转变的挤压作用,发生压剪性脆性变形,200~500℃,差应力65.91~135.68MPa,应变测量轴比a/c为1.403~2.204,动态重结晶石英颗粒边界分维值为1.321~1.378,反映矿区成矿期为中低温中高压低应变带变形环境,反应速率较小;成矿后在NWW-SEE向挤压作用下发生压剪变形,150~300℃,反映低温低压脆性变形环境。  相似文献   

8.
西昆仑山库地变质橄榄岩的显微构造特征   总被引:4,自引:1,他引:4       下载免费PDF全文
边千韬  林传勇 《地质科学》1995,30(3):275-282
西昆仑山库地变质橄榄岩块是肢解了的震旦纪一早古生代蛇绿岩碎片。变质橄榄岩以残斑结构为主,呈现粗粗结构向残班结构的过渡和残斑结构向板状等粒结构的过渡。此岩石有明显的变形结构,橄榄石变形时的滑移系为{0kl}[100].橄揽石发育有代表高温(800℃以上)的位错构造,也有代表低温的位错构造。表明此变质橄榄岩经受了上地幔环境下发生的高温塑性流变和侵位过程中及侵位后所经受的低温高应变速率条件下的变形。上地幔环境下变形时的差异应力约为80MPa,应变速率约为1.96×10-12s-1-7.03×10-10s-1.此变质橄榄岩代表由大陆裂谷发展起来的洋盆下的地幔,产出的构造环境为洋脊-裂谷带。  相似文献   

9.
江南造山带东段皖赣相邻区璜茅-五城-屯溪韧性剪切带野外现象、室内显微构造特征表明其经历了较为强烈的韧性剪切变形。岩石经历了NNE向右行剪切改造,变质变形程度为初糜棱岩-糜棱岩阶段,少量达超糜棱岩阶段。结合石英普遍发育边界迁移重结晶、长石多以塑性拉长和少量膨凸重结晶的变形行为特征、石英C轴组构显示出中-高温柱面和菱面滑移为主、重结晶石英颗粒大小-频数及分形维数(1.1646~1.2007),Kruhl温度计测算,揭示了璜茅-五城-屯溪韧性剪切带以中高温位错蠕变变形为主,变形和糜棱岩化环境在450~600℃的中深地壳范围,相当于中高绿片岩相-低角闪岩相。同时,岩石变形的古差异应力为20.76~30.04 MPa,估算所得应变速率主要介于10~(-14)~10~(-11) s~(-1),指示了中高温、低应变速率的韧性变形条件。已有的变形年龄显示,该韧性剪切带为加里东构造事件产物,主要的动力来源与江山-绍兴断裂带发育的陆内俯冲作用主应力在远离造山带的转变相关。  相似文献   

10.
Carrara大理岩高温高压变形实验研究   总被引:1,自引:0,他引:1  
利用高精度的Paterson高温高压流变仪对Carrara大理岩在高温(873~1173K)高压(~300MPa)以及约10-6~10-3s-1应变速率下进行了三轴压缩变形实验。结果表明,在等应变速率条件下,其强度随着温度的升高而降低;在等温和等压条件下,其强度随着应变速率的增加先快速增加而后缓慢增加。在应变速率对差应力的双对数投图中,我们发现随着温度的升高拟合直线的斜率减小,并且在873K和高应变速率时973K温度下Carrara大理岩的流变本构方程服从指数律变化关系;而在高温(1073K和1173K)和973K低应变速率条件下Carrara大理岩的应力指数n为5.3~7.7,且服从幂次律变化关系。因此,Carrara大理岩在本研究的实验条件下主要有两种变形机制,一种是用指数律表示的高应力变形机制;另一种是用幂次律表示的中等应力变形机制。  相似文献   

11.
污染场地健康风险评价的实例研究   总被引:35,自引:1,他引:34  
以常州市A厂有机污染调查为依据,采用篇中讨论的污染场地健康风险评价方法综合评价由土壤、空气、地下水污染带来的人群健康风险。综合分析表明,场地目标污染物为4种苯系物,可能对人体产生肝肾病变、鼻上皮组织恶化和运动协调能力降低等非致癌效应;污染物暴露途径主要有3条,分别暴露于厂区人群和饮用厂区下游民井井水的居民。评价结果显示,厂区人群同时遭受皮肤接触污染土壤和呼吸污染空气带来的非致癌危害。4种污染物的多途径同种污染物累计非致癌危害指数分别为甲苯,2·82×10-2;乙苯,4·94×10-2;间对二甲苯,1·21×10-7;异丙苯,1·05×10-7。两途径多污染物非致癌总危害指数分别为4·26×10-5和7·76×10-2,非致癌综合危害指数为7·76×10-2;厂区下游居民饮用遭受污染的地下水,其多污染物非致癌总危害指数为1·55×10-5。各危害指数均小于1,非致癌健康风险较小。  相似文献   

12.
出露于辽西白塔地区的侏罗纪中酸性火山岩及火山碎屑岩组合,主要是由粗安岩、流纹质岩屑晶屑熔结凝灰岩、粗面英安岩和流纹岩组成。通过系统的LA-ICPMS锆石的U-Pb同位素年代学研究,研究区火山岩形成年代为165~149 Ma。岩石地球化学特征显示,酸性火山岩属于钙碱性火山岩系,中性火山岩分别属于碱性火山岩系和中-高钾钙碱性火山岩系,中性岩具有埃达克质岩的地球化学特征,SiO_2≥56.20%、Al_2O_3≥15.09%、MgO≤2.89%、Sr≥561×10~(-6)、Y≤13.15×10~(-6)、Yb≤1.47×10~(-6)和无明显Eu异常。研究区火山岩普遍具有富集大离子亲石元素(Rb、K、Ba、Sr、Pb)和轻稀土元素,亏损高场强元素(Ti、Nb、Ta、P)的特征。其原始岩浆来源于基性下地壳物质的部分熔融,形成于古太平洋板块俯冲影响下板内挤压造山的构造环境。  相似文献   

13.
陈立华  王焰  易凯  赖河涛 《水文》2016,36(6):89-96
依据钦州市58a平均降雨量和3条入海河流控制站径流量长序列资料,采用滑动平均、线性回归、Spearman、M-K及R/S法综合分析降雨径流的趋势性及突变特征。结果表明,钦州市降雨量总体呈弱增加趋势,增加率为0.742mm/10a,而茅岭江、钦江、大风江流域径流量存在总体减少趋势,减少率分别为0.2×108m~3/10a、1.1×108m~3/10a、0.4×108m~3/10a。运用复Molet小波分析多时间尺度周期性,降雨量序列存在5个时间尺度,22a和15a时间尺度分别为序列第一、三主周期;径流量序列存在3个时间尺度,其中22a和8a时间尺度分别为径流量序列第一、二主周期。  相似文献   

14.
Strain rate during testing, uniaxial or triaxial, has important influence on the measured mechanical properties of rocks. Uniaxial compression tests were performed at nine pre-specified static-to-quasistatic strain rates (ranging from 1 × 10?5 to 1 × 10?1 s?1) on coarse crystal marble. The aim is to gain deep insight into the influence of strain rate on characteristic stresses, deformation properties and conversion of strain energy of such rock. It is found that the strain rate of 5 × 10?3 s?1 is the threshold to delineate the failure modes the tested coarse marble behaves in. At a strain rate less than this threshold, single-plane shear and conjugate X-shaped shear are the main failure modes, while beyond this threshold, extensile and splitting failures are dominant. The stress for crack initiation, the critical stress for dilation, the peak stress, and Young’s modulus are all found to increase with strain rate, with an exception that the above stresses and modulus appear relatively low compared to the strain rate in the range of between 1 × 10?4 and 5 × 10?3 s?1. The pre-peak absorbed strain energy, damage strain energy and elastic strain energy are found to increase with strain rate. In addition, the elastic strain energy stored before peak point favors brittle failure of the specimen, as the more stored elastic energy in the specimen, the stronger the fragmenting.  相似文献   

15.
地下水有机污染人体健康风险评价初探   总被引:47,自引:0,他引:47  
健康风险评价是定量描述污染对人体健康产生危害的重要方法,目前国内主要应用于地表水或污水回用的评价。文中针对地下水中有机污染物,考虑中国人饮水习惯及有机污染物的自然衰减作用,对U·S·EPA推荐模型进行了改进,并以北方某市一典型有机污染区的地下水为例,对地下水中污染物通过食入和皮肤接触两种途径进入人体产生的危害进行了风险计算和评价,分析了其主要风险来源。结果表明,典型区各点的非致癌风险均未超标,但有4个点的致癌风险超过U·S·EPA推荐的可接受风险值(1·0×10-4),其中B408点致癌风险高达1·37×10-3,不宜作为饮用水水源;各个点风险的主要贡献因子均是饮水途径摄入的三氯乙烯和四氯乙烯。但饮用煮沸的水在很大程度上能降低风险,建议不饮用生水。  相似文献   

16.
Current deformation in Pribaikalia, Western and Central Mongolia, and Tuva has been studied from measured horizontal GPS velocities and respective computed strain and rotation rates using 1994–2007 data of the Baikal–Mongolian GPS triangulation network.The GPS velocity field shows two main trends: an NE trend within Jonggaria, the Mongolian Altay, and the Great Lakes Valley and an SE trend in the Hangayn and eastern Gobi Altay mountains, and in the Transbaikalian block of the Amur plate. The velocity magnitudes and vectors are consistent with an SE motion of the Amur plate at a rate of ~2 mm/year.The derived strain pattern includes domains of crustal contraction and extension recognized from the magnitudes of relative strains. Shortening predominates in the Gobi and Mongolian Altay and in the Khamar-Daban Range, where it is at ?2 = (19.2 ± 6.0)×10?9 yr?1 being directed northeastward. Extension domains exist in the Baikal rift and in the Busiyngol–West Hangayn area, where the crust is stretching along NW axes at ?1 = (22.2 ± 3.1) × 10–9 yr–1. The eastern Hangayn dome and the Gobi peneplain on its eastern border show low and unstable strain rates. In central and northern Mongolia (Orhon–Selenge basin), shortening and extension are at similar rates: ?2 = (15.4 ± 5.4)×10?9 yr?1 and ?1 = (18.1 ± 3.1)×10?9 yr?1. The strain pattern changes notably in the area of the Mogod earthquake of 1967.Most of rotation throughout Central Asia is clockwise at a low rate of about Ω = 6×10?9 deg·yr?1. High rates of clockwise rotation are observed in the Hangayn domain (18.1 ± 5.2)×10?9 deg·yr?1, in the Gobi Altay (10.4 ± 7.5)×10?9 deg·yr?1, and in the Orhon–Selenge domain (11.9 ± 5.2)×10?9 deg·yr?1. Counterclockwise rotation is restricted to several domains. One is in western Tuva and northwestern Great Lakes Valley of Mongolia (Ω = 3.7×10?9 deg·yr?1). Two more counterclockwise rotation regions occur on both flanks of the Baikal rift: along the craton edge and in basins of Transbaikalia on the rift eastern border, where rotation rates are as high as (13.0 ± 3.9)×10?9 deg·yr?1, while rotation within the Baikal basin does not exceed the measurement error. Another such domain extends from the eastern Hövsgöl area to the Hangayn northern foothills, with the counterclockwise rotation at a highest rate of (16.3 ± 2.8)×10?9 deg·yr?1.  相似文献   

17.
埃达克岩是一套钙碱系列岩浆岩,其主要地球化学特征为SiO2 ≥ 56%,Al2O3 ≥ 15%,MgO通常小于3%,低重稀土元素和Y(Y ≤ 18×10-6,Yb ≤ 1.9×10-6),高Sr(多数大于400×10-6),无Eu异常或有轻微的负Eu异常,其成因主要有两类,一种为由俯冲板片的熔融形成,另一种由底侵玄武质下地壳熔融形成,埃达克岩具有重要的成矿意义。  相似文献   

18.
为确定临江荒沟山铅锌矿床成矿流体特征及矿床形成机制,笔者对矿区主成矿阶段(多金属硫化物阶段)中的铅锌矿石进行了流体包裹体岩相学、显微测温及稳定同位素研究。结果显示,包裹体均为气液两相包裹体,属NaCl-H_2O体系,且成矿流体具有中温(209℃~276℃)、中盐度(1.98~16.34wt%)、低密度(0.79~0.96 g/cm~3)的特点;成矿压力约为16.8~33.0 MPa,成矿深度为1.68~3.30 km。氢氧同位素结果显示δ~(18)O_(H_2O)和δD值分别为3.3×10~(-3)~9.9×10~(-3),-115×10~(-3)~-74.8×10~(-3),结合碳氧同位素δ~(18)O-δ~(13)C图解说明成矿流体以幔源流体为主,可能有大气降水混入。硫同位素δ~(34)S值为2.6×10~(-3)~18.9×10~(-3),具有富集重硫的特点;铅同位素~(206)Pb/~(204)Pb,~(207)Pb/~(204)Pb,~(208)Pb/~(204)Pb值分别为15.43~15.57,15.25~15.33,34.88~35.07,说明成矿物质来源于珍珠门组地层,但不排除有深源物质混入。结合区域地质,认为该矿床的成因类型为中温热液脉型铅锌矿,形成于环太平洋体系向华北板块俯冲的构造背景下,形成时代为燕山晚期。  相似文献   

19.
The carbonate aquifers of Lora and Mingo form part of the hydrogeological unit of Sierra de Estepa (SE Spain). By means of time series analysis and a 1D numerical groundwater model, groundwater exploitation was quantified and the mean annual recharge in both systems was estimated (2001–2004). During this period, the Lora and Mingo aquifers received an average groundwater recharge of 0.29 × 106 m3/year and 0.14 × 106 m3/year, respectively, whereas an average of 0.34 × 106 m3/year and 0.21 × 106 m3/year, respectively, was extracted. These conditions led to a conspicuous lowering of the water table in both systems. In addition, the analysis of the evolution of the main hydrogeochemical parameters of the groundwater showed that the increased pumping rates produced an increase in total dissolved solids, and chloride and sodium ions in both aquifers. In the case of the Lora aquifer, the only ion that presented decreased levels was nitrate. The results show that groundwater pumping in both aquifers should not exceed the mean annual recharge of 0.29 × 106 m3/year and 0.14 × 106 m3/year in the Lora and Mingo aquifers, respectively. Nevertheless, it would be advisable to reduce pumping rates to below these values in order to restore piezometric levels and improve groundwater quality for different uses in the future.  相似文献   

20.
滇西主高黎贡韧性剪切带糜棱岩形成时限的初探   总被引:5,自引:2,他引:3  
构造形成时限的研究,可以为相关地质事件提供重要的制约条件。本次工作利用应变速率法进行分析,以主高黎贡韧性剪切带长英质糜棱岩为实例,通过石英自由位错法测得其差异应力130.69MPa,进而得出应变速率为2.906×10-13s-1;利用质量等比线法分析得出其体积因子为0.7314;进而可以求出构造形成时限为3.41万年。这一时间段的确定,为特提斯构造系形成时限的研究,给出了一项可供参考的数据,同时也表明应用应变速率法研究构造形成时限是可行的。  相似文献   

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

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