首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 375 毫秒
1.
龙门山断裂带多参数深部结构成像与地震成因研究   总被引:6,自引:4,他引:2       下载免费PDF全文
通过反演大量的纵、横波地震数据,获得了沿龙门山断裂带及周边区域的深部三维精细结构,结合前人二维大地电磁探测研究成果,提出龙门山断裂带地壳形变与深部速度结构和导电率不均匀性有关,探讨了2008年汶川和2013年芦山地震的诱发和产生与流体侵入及地壳形变的密切关系.本研究发现,2008年汶川地震发生在高速度、高泊松比和低电导率的区域,2013年芦山地震则位于高速度、低泊松比和低电导率的发震层.在上地壳中,四川前陆盆地的低速、低泊松比和低阻异常与松潘一甘孜地块的高速、高泊松比和高阻异常形成了鲜明的对比.在龙门山断裂带下方的两个低速和低阻块体,将龙门山断裂带分成南、中和北三段.我们的研究认为,这两个异常体与来自松潘甘孜地块的下地壳和(或)上地幔的局部熔融或流体侵入到龙门山断裂带的脆弱区有关.基于对汶川和芦山地震的余震分布特征及震源区的地震波速度、泊松比及电阻率参数分析,揭示了龙门山断裂带深部剧烈的地壳形变与流体应力积累对2008年汶川和2013年芦山地震的触发及其地震破裂过程具有重要的控制作用.  相似文献   

2.
通过反演由大量的纵、横波地震数据组成的综合数据集,获得了南北地震带地壳的多参数三维精细结构,探讨和分析了南北地震带的高地震活动性和强震频发的原因.成像结果表明,尽管1976年松潘一平武地震(M7.2)与2008年汶川地震(M8.0)以及2013年芦山地震(M7.0)均发生在高速、低泊松比异常区域,并且在其震源的下方均有一低速、高泊松比异常区域.我们认为,上述三个地震的触发与流体侵入导致的地壳形变之间有密切的联系.1955年炉霍地震(M7.4)和1973年康定地震(M7.1)均发生在鲜水河断裂带上,其震源中心区域表现为低速、高泊松比异常,可以解释为下地壳中的流体沿断层面上涌.在震源区的周边区域兼有高速、高泊松比异常,低速、高泊松比异常以及高速、低泊松比异常,可能分别与含流体的岩石、沿断裂带发育的变质岩以及坚硬的克拉通块体对应.流体的侵入不仅能够改变断层面上的应力情况,还能降低岩石骨架的岩石力学强度,进而触发地震.1970年云南通海大地震(M7.1)发生在哀牢山一红河断裂带附近的曲江断裂上,其震源处于高速度、低泊松比异常与低速度、高泊松比异常之间的边界区域,被认为是流体挤压后的应变能积累,最终导致脆性破裂,以至于发生地震.根据本次研究获得的多参数结构图像,结合前人的研究成果,我们认为南北地震带地壳强烈形变与流体侵入是造成该区域地震活动性较高及强震频发的两个主要因素.  相似文献   

3.
为了解2013年11月23日吉林前郭MS5.8和2018年5月28日吉林松原MS5.7地震震源区的地壳速度结构及其孕震环境,本研究收集了东北地区CEA固定台站、NECESSArray流动台站和五大连池WAVESArrary流动台站的连续波形数据,使用背景噪声成像方法,获得了前郭地震和松原地震震源区地壳S波速度结构.结果显示,这两个地震主要发生在相对高波速异常区,而震源区下方存在明显低波速异常,且该低速层呈北北东走向,并有向北延展的趋势.结果表明,震源区相对高波速异常区易于积累能量孕育地震,而震源区下方的低波速异常可能代表了流体作用.这种流体作用会降低断层面有效正应力从而触发地震.这种流体作用可能和太平洋板块深俯冲至我国东北地区(包括前郭地震和松原地震震源区)下方的地幔转换内形成"大地幔楔"结构与动力学密切相关.在"大地幔楔"结构中,由于地幔转换带中滞留板块脱水作用和地幔角流作用,容易形成热湿物质上涌,进而导致松辽盆地中北部岩石圈物质拆沉、地幔流体入侵至中下地壳、然后作用于断裂带,从而导致了两个震源区中强度地震的发生.  相似文献   

4.
张渤地震带及邻区近震体波成像及孕震环境分析   总被引:3,自引:0,他引:3       下载免费PDF全文
针对张渤地震带深部孕震构造环境和地幔岩浆对地壳底侵作用的问题,本文充分收集了华北地区区域数字地震台网176个固定台站记录到的观测数据,应用近震体波层析成像方法,获得了华北地区地壳的P波三维速度结构和VP/VS波速比.依据VP/VS波速比能反映岩石物理性质和流变学特征的特性,并结合人工地震测深剖面、大地电磁测深、地球化学等成果,讨论了张渤地震带的孕震环境和动力学机制等问题.研究结果显示:张渤地震带强震位于高低速异常分界线附近或偏向高速体一侧,震源体下方均存在广泛的低速异常分布,据此推测这些低速异常的存在可能与流体有关,同时也揭示出张渤地震带地震主要分布在上地壳,从震源深度分布及地壳物性结构分析上地壳具有发震的构造背景,也有发震的物性基础,是一个易震层和多震层.张渤地震带地壳内部在20 km下方存在偏高的波速比特征,这些可能反映了在该区域内中、上地壳幔源物质的侵入和热状态的岩体在横向上的变化,可能是地幔岩浆长期底侵作用的结果.由于深部幔源物质的侵入,使得地壳深部流体的供给量增加,在地壳发震层下长期存在的流体会影响断裂带的结构,降低断裂带的强度,使区域应力场发生变化从而导致断裂带上应力的集中,进而引发地震发生.  相似文献   

5.
2021年5月21日漾濞MS6.4地震震源区位于川滇块体西边界的维西—乔后—巍山断裂西侧,该地区近些年来发生了多次中强地震,地震活动较为活跃.对漾濞地震序列重定位和漾濞震源区及邻区的地壳精细结构研究,有助于深入理解漾濞地震的孕震环境、发震机理和破裂过程.本文基于2008年 1月1日到2021 年6 月3日区域固定台站接收到的 36938条Pg和 32111 条 Sg波到时数据,采用新发展的三重差地震层析成像算法(tomoTD)开展了漾濞MS 6.4地震震源区三维速度结构成像与地震重定位研究.结果显示:(1)余震活动主要集中在维西—乔后—巍山断裂的西侧,整体呈现沿北北西向的条带状分布,结合已有走滑型震源机制解特征,揭示了北西向隐伏断裂是发震断层,其北西段表现为倾角较陡、结构相对简单的走滑断裂,南东段由两条分支断裂组成.(2)主震的发生及地震序列分布与地壳速度结构不均匀性有着密切的关系.主震及4级以上的地震发生在高速边界上或高低速过渡区域,余震主要发生在低速、高VP/VS 区,主震上方与下方均显示高VP/VS 异常,推测在区域构造应力场的作用下,应力在孕震区的刚性介质中积累,中下地壳流体(或者部分熔融地壳物质)侵入发震断层区,弱化了漾濞 6.4 级地震的主震区.另外,余震东南侧的低VP/VS 区可能代表介质刚性强,可能阻碍了余震向南东方向继续扩展.(3)结合2013 洱源 5.5 级地震研究结果,推测维西—乔后—巍山断裂西侧可能存在着较大的北北西向隐伏断层.  相似文献   

6.
芦山与汶川地震之间存在约40 km的地震空区.震源区和地震空区的深部构造背景的研究对深入了解中强地震的深部孕育环境及地震空区的地震活动性具有重要科学意义.利用本小组布设的15个临时观测地震台以及21个芦山科考台站和21个四川省地震局固定台站记录的远震数据,用H-K叠加方法得到各个台站的地壳厚度和平均泊松比,并构建了接收函数共转换点(CCP)偏移叠加图像以及反演得到台站下方的S波速度模型.我们的结果揭示了震源区和地震空区地壳结构特征差异:(1)汶川震源区的地壳平均泊松比为~0.28;芦山震源区为~0.29;而地震空区处于泊松比变化剧烈的区域;(2)汶川地震与芦山地震的震源区以西下方的Moho面呈现深度上的突变(这与前人的研究成果基本一致),分别从~44 km突变到~59 km,~40 km突变到~50 km,而地震空区地壳平均厚度呈现渐变性变化;(3)地震空区Moho面下凹且具有低速的上地壳.综合一维S波速度结构和H-k以及CCP的初步结果,这可能显示汶川地震的发震断裂在深部方向上向西倾斜并形成切割整个地壳的大型断裂;芦山地震则可能是由于上、下地壳解耦引起的;而地震空区处于两种地震形成机制控制区域的过渡带中.  相似文献   

7.
利用西藏自治区林芝地区的固定地震台站与南迦巴瓦流动测震台站在2017年11月18日至2017年11月24日记录到的430个余震的直达波走时数据反演得到了震源区的三维P波速度、S波速度结构,并利用三维速度结构对余震进行了重定位.成像结果显示,米林地震震源区在0~5km深度内存在低地震波速度异常;在5~15km深度内,存在高地震波速度异常,该高速异常致使震源区西南侧的地震波速度高于东北侧.重定位结果中,余震呈条带状以NW-SE走向展布,震源深度具有西南方向深、东北方向浅的特征.主震位于11km深度处、高地震波速异常体顶部,余震主要分布在高地震波速度与低地震波速度过渡的区域.对成像结果的分析表明,震源区浅部的低速异常具有低泊松比的特性,与富石英的沉积变质杂岩体-东久杂岩单元的岩性特征有关;深部的速度结构特征则可能反映了发震断层上盘地震波速度高,下盘地震波速度低的介质特性.余震重定位结果与成像结果联合表明:此次地震发震断层从11km深度处,东久杂岩体下方的高地震波速度异常顶部开始破裂,继而在5~15km深度内发生后续破裂,后续破裂的发生区域正处于喜马拉雅构造单元与冈底斯构造单元接触的形变区内.此外,根据地震波速度计算的泊松比反映了震源区持续的低泊松比特征,暗示此次地震与流体活动并无直接关系.  相似文献   

8.
张莹莹  高原  石玉涛  Kelly Liu 《地震学报》2015,37(4):541-1093
利用山东、 河北、 辽宁、 北京、 天津等区域地震台网记录的远震资料, 采用接收函数H-k叠加搜索方法, 得到了张家口—渤海地震活动带及其邻区99个台站下方的地壳厚度和平均泊松比估计. 研究结果显示: 研究区地壳厚度在30—46 km之间显著变化, 地壳厚度与地表地形相关性很好, 呈现出从隆起区到盆地区逐渐减薄的变化趋势; 研究区地壳泊松比具有明显的分区特点, 燕山地区泊松比偏低可能与该地区克拉通破坏方式以拆沉为主有关, 郯城—庐江断裂带山东段出现的泊松比高低混杂现象可能与该地区火山活动有关, 张家口—渤海地震活动带西段与隆起的交汇地带, 部分台站下方地壳的泊松比值偏高, 可能是断裂带附近岩石比较破碎与岩石部分熔融共同作用的结果. 此外, 位于渤海湾中间岛链上的北隍城台下方地壳较厚, 泊松比较低, 其原因尚待查明.   相似文献   

9.
芦山与汶川地震震区地壳上地幔结构及深部孕震环境   总被引:1,自引:0,他引:1  
雅安芦山Mw6.6级地震,是汶川7.9级地震之后龙门山地区的又一次强震,给当地的社会经济发展和生态环境带来了巨大的破坏.本文借助全国地震台网连续波形数据,使用背景噪声层析成像方法和远震接收函数分析方法,获得了震区及其周边地区精细的S波速度结构和地壳厚度、泊松比分布情况,进而分析了龙门山断裂带及邻区深部孕震环境.研究结果发现:(1)龙门山断裂带两侧剪切波速度和地壳厚度有非常显著的差异;(2)芦山和汶川地震均位于地壳厚度和波速结构变化剧烈之处,断层的破裂面和余震的分布均处于地震波横向速度梯度和地壳厚度的横向梯度跳变的地区;地震深度处于从均匀波速结构向非均匀波速结构过渡区域;(3)汶川地震及其余震发生在泊松比较低的地区,而芦山地震及其余震发生在高泊松比的地区.我们推测,横向的显著速度变化,以及地壳厚度和地形高度的巨大差异所产生的重力势能等作用可能在一定程度上构造成了龙门山断裂系上大地震的孕震环境,而龙门山断裂带南北段岩石泊松比的差别,可能是导致汶川地震余震空间分布和芦山地震延后5年发生的原因.  相似文献   

10.
为了理解地震的发生和地壳结构不均匀性的关系,利用南加州地震台网的固定台和临时台所记录的2863个兰德斯余震和区域地震,共计107401个P波和19624个S波高质量的到时数据,采用地震层析成像方法得到了兰德斯地震区P波和S波的精细的三维速度结构和泊松比分布.结果显示,地震的发生和分布与地壳结构的横向不均匀性有密切的关系.总体上看,兰德斯地震区余震成丛分布,并被低速块体截断,其中4级以上地震大多分布于P波高、低速异常过渡区域或偏向高速块体一侧,这可能是因为高速区多属地壳脆性介质,易于造成应力集中,导致地震;反之,低速度区则可能代表破碎程度较高、富含流体或温度较高区域,因而更倾向于产生无震变形.基于兰德斯地震区强震震源位置、地震区P波、S波速度异常与泊松比分布推断,兰德斯地震区可能有流体存在.地壳流体易使地壳岩石弱化,从而引发大地震.  相似文献   

11.
本研究利用中国国家地震台网336个固定地震台站记录的远震波形资料,通过P波接收函数H-κ分析估算了中国华南地区的地壳厚度和地壳平均波速比.研究结果显示,地壳厚度约为25~46 km,整体表现西深东浅的变化特点.研究区地壳厚度变化分别与布格重力异常和地形呈现负相关和正相关.扬子块体东部显示较低地壳波速比(<1.7)可能与地表沉积岩的存在相关.四川盆地内部的平均地壳波速比(约1.71~1.8)与全球大陆地盾/克拉通地区相当.但其周围地区地壳平均波速比明显增高(1.81~1.95),推测可能是克拉通形成过程中岩浆的底侵引起下地壳组分以铁镁质麻粒岩为主.华夏块体西部表现为薄的地壳厚度和低的地壳平均波速比(1.65~1.75),暗示该区基性下地壳物质的缺失可能与增厚的地壳发生拆沉有关.华夏块体东南段呈现出较高的地壳平均波速比(约1.77~1.86),地壳组成以中基性铁镁质为主,推断可能与晚中生代铁镁质岩浆底侵作用密切相关.  相似文献   

12.
贾漯昭  王志铄  张亚琳  邢康 《地震》2017,37(1):112-120
使用扣除同一个地震台站记录到的两个距离相近地震的地震波的相同路径, 把波速比指示范围限制在地震集中区的方法计算波速比值。 利用区域台网记录到的地震波形和震相数据, 对2014年8月~2015年8月间, 安徽金寨地区发生的小震震群进行波速比计算, 分析了该方法在该地区震相数据上的误差影响和应用情况。 得到金寨震群震中集中区的波速比值。 结果显示, 震中集中区波速比略低, 结合当地断层分布情况, 这可能体现了该区域的介质硬度稍高, 应力积累高于周边区域。  相似文献   

13.
本文利用区域地震初至波到时数据,通过地震层析成像研究获得了东北日本俯冲带上地幔(深至约150 km)的P波速度(VP)、S波速度(VS)、VP/VS和P波各向异性结构.结果表明,低速及高VP/VS比异常体主要分布在火山下方的下地壳和地幔楔中,其与低频地震的分布吻合,该区域与俯冲板块脱水所释放的流体及其导致的部分熔融密切相关;俯冲的太平洋板块内可能由于脱水脆化导致的双层地震带区域则没有表现出整体的高VP/VS值,其可能与俯冲板块内部含水矿物含量有关;俯冲板块内双重地震带区域及上覆地幔楔薄层主要表现为与海沟平行的方位各向异性和正的径向各向异性,其可能是由于含水矿物的脱水使橄榄石晶格结构发生了从A型到B型的变化所引起的.我们研究表明,结合地震波速度和各向异性结构能够加深对俯冲带内水运移过程的认识.  相似文献   

14.
海原弧形构造区地壳三维精细速度结构成像   总被引:4,自引:0,他引:4       下载免费PDF全文
利用海原弧形构造区及周围区域地震台网1970—2015年期间记录的天然地震到时数据,采用双差地震层析成像方法对构造区地壳三维速度结构与地震震源位置进行联合反演,获得了高分辨率的三维VPVS以及VP/VS模型,分析讨论了速度、波速比分布与强震发生以及断裂等之间的关系.结果显示:研究区域内地震主要沿断裂呈弧状展布,速度在横向分布上具有较大的差异,波速比变化范围为1.60~1.80,平均值约为1.70.大型断裂诸如海原—六盘山断裂带、青铜峡—固原断裂带等位于高速与低速的过渡带,断裂两侧地震波速差异较大.研究区内历史强震多处于高低速过渡区域,海原强震下方下地壳存在低速、高导薄弱层(25~30 km深度),推测原因主要为流体作用所致.依据相对较低的速度与波速比分布推测研究区地壳主要组成成分为酸性的长英质.速度剖面显示地壳可分为上、下两层,上、下地壳厚度变化由西南向东北逐渐减薄,减薄幅度相近;结合前人研究结果推测构造区地壳增厚模式可能主要为上、下地壳共同增厚.  相似文献   

15.
远震接收函数已广泛用于反演台站下方的结构,然而由于地球的非弹性衰减作用,远震数据较难获得高频接收函数,对浅地表结构约束不足.为了克服这一问题,我们使用近震数据的高频接收函数来研究浅表速度结构,并应用于四川理县西山村滑坡体上3个宽频带地震仪记录到的近震事件.本文发展了接收函数V_P-k(V_P/V_S)叠加方法,结合接收函数H-k叠加和波形反演方法获得了台站下方滑坡体的厚度、S波速度和平均V_P/V_S比,并与钻孔得到的滑坡体厚度进行对比.结果表明,滑坡体具有小尺度的横向不均匀性,台站下方滑坡体的平均V_P/V_S比在2.4~3.1之间变化并且在底层存在78~143m·s~(-1)左右的S波低速层.本文观测到的高V_P/V_S比和底层低的S波速度结构,与电磁法获得的滑坡体底层低的电阻率和底部富水特征一致,表明滑坡体h1底界面的抗剪强度相对较弱,是潜在的滑坡危险区域.本文研究结果表明,利用近震接收函数能有效约束浅表的速度结构,进而能为滑坡灾害治理提供一定的地震学参考.  相似文献   

16.
在喜马拉雅造山带的东缘,雅鲁藏布江缝合带在这里发生急剧转折,南迦巴瓦变质体快速隆起,然而关于东构造结的形成机制一直未有定论.利用围绕南迦巴瓦峰的48个宽频带地震台站记录的远震数据提取P波接收函数,采用改进的H-κ叠加方法和共转换点叠加方法综合研究了东构造结的地壳厚度、波速比分布和地壳结构特征.结果表明:研究区平均地壳厚度为64.03 km,大部分台站介于60.48~66.55 km范围;平均波速比为1.728,主要集中范围为1.696~1.742.东构造结地壳厚度横向变化剧烈,构造结西端和北端厚而中间薄,东构造结核部Moho面呈现上隆的构造形态,东西向上隆幅度约为6~7 km,南北向的上隆超过9~10 km.东构造结核部地壳上隆减薄可能由高密度、高波速的岩石圈撕裂残片拆沉到上地幔软流圈后重力失衡所致.平均波速比超过1.8的高值异常展布于东构造结的两侧,推测为环东构造结的壳内部分熔融体.东构造结地壳上隆减薄和壳内部分熔融的存在很可能均与幔源热物质的上涌有关,而软流圈地幔的上涌则可能由印度板片的撕裂引起.  相似文献   

17.
兴蒙造山带诺敏河火山群地壳厚度与波速比研究   总被引:1,自引:1,他引:0       下载免费PDF全文
利用布设于兴蒙造山带诺敏河火山群地区的宽频带流动地震台站资料,基于接收函数方法,获取了该地区的地壳厚度与波速比值.研究结果显示,该地区的地壳厚度介于32~38 km,莫霍面深度在空间上分布特征与五大连池为中心的火山带分布具有较好的一致性:沿着火山带延展方向地壳较薄.该地区的波速比介于1.74~1.84,波速比在空间上与地壳厚度变化具有一致性:高波速比主要集中于靠近五大连池火山带地区,向诺敏河火山和小古里河火山延展.研究认为:诺敏河火山与五大连池火山带可能具有相同的岩浆来源,可能与富钾岩石圈地幔拆沉作用造成的地幔热物质上涌有关.研究区地壳厚度与波速比呈现负相关关系,暗示该地区可能发生过岩浆底侵作用.  相似文献   

18.
朱洪翔  田有  刘财  冯晅 《地球物理学报》2018,61(9):3664-3675
接收函数方法被广泛地应用于地壳上地幔结构的研究中,H-κ叠加方法是其中最常用的方法之一.对于布设在基岩区台站计算的接收函数,H-κ叠加方法可以准确地估计台站下方地壳厚度和平均波速比,但是对于沉积盆地地区计算的接收函数,由于低速沉积层内会产生多次波混响,干扰甚至覆盖接收函数中莫霍面的转换波和多次波震相,从而影响H-κ叠加结果的准确性.为准确估计沉积盆地地区地壳结构,本文提出使用预测反褶积方法去除接收函数中低速沉积层内多次波混响,其中预测步长由接收函数归一化自相关函数获得,物理意义为沉积层内S波双程走时.合成接收函数和实测接收函数试验表明,本文提出的预测反褶积方法可以有效地去除沉积层多次波混响,并结合改进的H-κ叠加方法可以准确地估计沉积层下覆地壳厚度和平均波速比.相比于其他去除接收函数多次波混响的方法,本文提出的预测反褶积方法具有参数设定简单、运算量小、震相幅值较大等特点,适用于大批量数据处理.  相似文献   

19.
基于我们布设的探测深俯冲的中国东北地震台阵NECsaids台阵和固定地震台长时段的观测记录及NECESSArray流动台阵共计152个台站数据提取得到的33752条P波接收函数,采用H-κ叠加分析和共转换点(CCP)叠加成像等方法进行统一分析处理,并汇集他人接收函数研究结果得到中国东北东南部地区迄今为止最高分辨率的地壳厚度和平均波速比分布图像.对中国东北东南部地区不同构造体的地壳特征综合分析研究表明:研究区不同陆块的地壳属性存在明显差别,张广才岭地块中南部的地壳厚度和波速比与华北克拉通东北缘相当,地壳厚度同地表地形之间显示有明显的正相关关系;松辽地块东南缘地壳最薄、波速比最高,地壳厚度同壳内波速比之间显示出明显的负相关关系;兴凯地块西部地区的地壳结构表现为稍厚的地壳厚度和研究区内最低的壳内波速比,其地壳厚度同壳内波速比之间亦显示出明显的负相关关系;佳木斯地块西南缘在具有"正常"的壳内波速比同时地壳最厚.研究区内的郯庐断裂带北延段在切穿其下Moho面的同时表现出南北分段的特征:北段(44.4°N—47°N)两分支之下的Moho面整体下凹,而南段(41.5°N—43.3°N)两分支之下的Moho面则整体上隆.长白山天池火山下方表现为Moho面下凹沉落及高壳内波速比特征,推测其壳内岩浆囊很可能存在于火山口东北隅至少10 km的范围内.  相似文献   

20.
针对构造演化历史甚为复杂的兴蒙造山带的地壳结构探测较为有限的问题,本研究利用2016-2018年布设在兴蒙造山带西南部的第3期NECsaids地震台阵和固定地震台及NECESSArray的流动台数据,采用时间域最大熵谱反褶积方法提取到研究区的24027条高质量P波接收函数,进而采用H-κ和共转换点(CCP)叠加方法以及分层剥离转换震相分析地壳各向异性的方法来辨识研究区的地壳结构特征.分析表明:研究区地壳各向异性整体呈近NW向和近EW向,与区域速度场和板块绝对运动方向一致,表明现今地壳结构主要受控于太平洋板块构造域;大兴安岭-太行山重力梯级带是明显的地壳结构差异过渡区,其西侧地壳厚度明显高于东侧且呈降低趋势,而平均波速比则呈现东西两侧盆地较高的分布特征;局部地区的高波速比和不同的各向异性特征则显示受到新生代火山和断裂活动的改造作用,二连盆地的厚地壳和复杂的各向异性特征则展现出古亚洲洋闭合及蒙古-鄂霍茨克和古太平洋等多期构造作用的可能影响.  相似文献   

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

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