首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 187 毫秒
1.

我们在安宁河断层带西侧的花岗片麻岩和碎裂岩中发现假玄武玻璃.并开展了假玄武玻璃的摩擦滑动实验,为国际首次.在化学成分和微观结构分析的基础上,本研究使用粉碎后的假玄武玻璃粉末作为模拟断层泥,在有效正应力250 MPa、温度100~500℃、孔隙水压30 MPa和干燥条件下开展了摩擦滑动实验.为获得假玄武玻璃摩擦系数的速度依赖性参数(a-b),实验中的轴向加载速率在1 μm·s-1,0.2 μm·s-1,0.04 μm·s-1之间进行切换.实验结果显示,在100~500℃假玄武玻璃的摩擦系数随着温度增加而增大,含孔隙水压条件下摩擦系数比干燥条件大.在加载速度为1~0.2 μm·s-1条件下,100~300℃温度范围内,假玄武玻璃表现为速度强化型稳定摩擦滑动;在300~400℃的温度范围内,假玄武玻璃的摩擦滑动行为表现出由速度强化向速度弱化转变的趋势,在400~500℃时其表现出显著的速度弱化行为,从而在该温度范围内具备自发的地震成核条件.本实验还发现假玄武玻璃的滑动稳定性与加载速率有关,在相同温度条件下以及100~300℃的温度范围内,假玄武玻璃在实验室标速条件下(1~0.2 μm·s-1)表现出速度强化的摩擦滑动行为;而在低速(0.2~0.04 μm·s-1)条件下,表现出速度弱化.因此在100~500℃的温度条件下,降低剪切滑动速度可以促进假玄武玻璃的不稳定滑动.对比花岗岩与假玄武玻璃的摩擦强度与滑动稳定性表明,在100~500℃(5~30 km深度),假玄武玻璃和花岗岩都具备地震成核的条件.假玄武玻璃的摩擦强度比花岗岩高,由此推测,被假玄武玻璃愈合的同震破裂,原地再次发生同震滑动的可能性小,但如果断裂带内存在假玄武玻璃,不仅能增加断层不稳定滑动和地震的风险,而且会促进同震动态破裂传播.

  相似文献   

2.

为了探讨在俯冲带构造加载环境中叶蛇纹石矿物的摩擦力学性质, 我们使用叶蛇纹石作为实验样品, 开展了高压流体条件下的摩擦滑动实验研究.实验在气体介质的高温高压三轴实验系统中进行, 实验中的有效正应力为30 MPa, 孔隙流体压力为100 MPa, 温度范围为100~500℃.实验过程中为了得到叶蛇纹石摩擦强度系数的速度依赖性, 我们将轴向加载速率在1.0 μm·s-1、0.2 μm·s-1和0.04 μm·s-1之间进行切换.实验结果发现在实验温度范围内叶蛇纹石的摩擦强度系数随着温度的升高表现出系统的降低趋势, 摩擦系数从100℃的0.81降低到500℃的0.41.在高于300℃的温度条件下, 叶蛇纹石的剪切力学曲线表现出显著的位移弱化现象.在300℃的温度条件下, 叶蛇纹石在0.2~1.0 μm·s-1的加载速率范围内表现出速度强化的摩擦行为, 而在0.04 μm·s-1的加载速率下表现为速度弱化.在100~200℃以及400~500℃的温度条件下皆表现为速度强化.我们的实验研究表明在俯冲带地幔楔附近的高压流体条件和低有效正压力环境下, 叶蛇纹石矿物可以促成慢滑移现象, 但是其背后的物理机制可能与已有的认知(不稳定的摩擦滑动由叶蛇纹石脱水过程所控制)不同.这为我们探讨俯冲带慢滑移现象的物理过程和机制提供了新的思路.

  相似文献   

3.

为探索断层岩石摩擦特性对于断层力学性质的影响,我们采集了龙门山汶茂断裂韧性剪切带中的富含层状硅酸盐矿物的糜棱岩样品进行了水热条件下摩擦滑动实验研究.实验在三轴压机之上完成,实验温度为100~600℃,有效正压力100 MPa,孔隙水压分别为30 MPa和130 MPa.为获得摩擦滑动的稳定性参数(a-b),剪切滑移速率在1.22 μm·s-1,0.244 μm·s-1和0.0488 μm·s-1之间切换.实验发现在200~500℃的温度范围内,摩擦系数随着温度的增加而显著增大(约0.56~0.72).在200~300℃范围内,随温度的升高糜棱岩的摩擦滑动表现出由稳定的速度强化向不稳定速度弱化转变的趋势.在有效正压力不变的情况下,孔隙水压的增大会促进糜棱岩的摩擦滑动在500~600℃温度范围内由不稳定的速度弱化向稳定的速度强化的转变.实验给出的断层在原地深度处的脆性和塑性变形机制的转变,有助于理解断层深部的地震成核机制以及成核的温压条件.

  相似文献   

4.
田平  何昌荣 《地球物理学报》2020,63(12):4440-4450

为了研究下地壳主要造岩矿物之一的辉石在摩擦滑动中的力学性质,本文选择辉石中的透辉石作为实验样品,在约200 MPa有效正应力(30 MPa孔隙水压)及100~600℃的热水条件下进行了速度阶跃摩擦滑动实验,在滑动开始后将剪切滑动速率在0.122 μm·s-1和1.22 μm·s-1之间进行切换来确定摩擦滑动对速度变化的响应.通过实验发现:透辉石在102~200℃表现为速度强化,随着温度升高,在约215℃附近从速度强化转变为速度弱化.在310℃以上,速度弱化的程度受温度的影响不大,而且摩擦本构参数a值和b值随温度的变化趋势相同.对变形样品的显微观测发现,与其速度强化行为相对应的微观变形主要以数条平行的R1剪切为主,卸样时造成的裂纹集中在R1剪切内部,形态长而曲折;与其速度弱化行为相对应的断层泥内部,剪切带变窄且数量明显增加,在内部互相连接,在边界附近多发育B剪切,卸样时形成的裂纹大量存在于B剪切带与R1剪切附近,形态短且平直.

  相似文献   

5.
文中采用尺寸为50mm×50mm×150mm的花岗闪长岩、大理岩以及砂岩3种不同的岩石样品,开展了超声尾波和声发射同步观测研究。结果显示:1)尾波变化与声发射演化过程存在良好的对应关系,声发射频度增加时尾波的变化随之改变,尤其是声发射频度增加的早期阶段尾波便出现变化。这预示了分析尾波可获得岩石早期的损伤信息。2)不同变形阶段,尾波变化的物理机制不同。加载初期,尾波变化存在明显的散射体改变特征;随后,在线弹性变形阶段则以波速变化为主;加载后期,散射体的改变增加,且与波速变化共存,散射体的改变效应与岩石微破裂的程度有关。3)随着载荷的增加,波速增加的幅度总体呈逐渐减小的趋势,这与通过直达波获得的认识基本一致。4)岩石变形过程中产生的微破裂会改变尾波变化的物理机制,散射体效应将显著增强。同时,声发射波形将对超声尾波造成干扰,因此在仅利用尾波资料分析岩石损伤时需要关注相关的问题。总之,超声尾波变化和声发射均可反映岩石内部的损伤情况,且不同的变形阶段尾波的变化机制不同,二者的联合观测可起到相互验证的作用,有利于提高观测结果的可靠性。  相似文献   

6.

流体运移时携带热量,会产生温度变化.温度作为地下流体的示踪剂得到了广泛应用,但通常需要事先确定热扩散系数.实际上,基于浅层地壳不同深度的周期性温度-时间序列,利用流体运动对振幅和相位的影响,可以同步获取热扩散系数和地下流体运移信息.本文以新疆喀什地区5个钻孔不同深度的基岩温度数据为基础,获取了热扩散系数和流体运移特征.主要成果有:(1)获得了不同测点的热扩散系数α,以及视热扩散系αAαΦ(即单独通过振幅或相位获得的热扩散系数).其中,ααAαΦ的值分别为1.52~8.91、0.79~1.71和(1.53~33.1)×10-6 m2·s-1.另外,当流体热效应不明显时,通过相位获得的视热扩散系αΦ更接近真实的热扩散系数α.(2)获得了不同测点的流体流向和流速信息.测点的流体流向均向上,靠近天山地区的测点流体流速为(0.10~1.94)×10-7 m·s-1,靠近昆仑山地区的测点流体流速为(8.56~9.71)×10-7 m·s-1,不同地区测点流体流速的差异可能与区域水文地质环境有关.总之,通过多深度的连续基岩温度观测,有望获得浅层地壳的热扩散系数及流体运移特征.

  相似文献   

7.

野外地质调查结果显示,断层带常富集碳质.断层带中碳的分布结构是影响断层带电导率特征的一种重要参数.本文在室温、室内湿度和2 MPa正应力条件下,对不同石墨含量(3,5,6和7 wt%)的石英-石墨混合断层泥模拟样品开展了滑动速率介于500 μm·s-1~1 m·s-1的摩擦实验及相应的电导率测量,以期研究断层运动对碳分布结构的影响以及断层带电性特征对碳含量及分布的响应情况.结果显示,摩擦滑动能够显著地改变样品的电性特征(电导率大小及其各向异性).在平行滑动面方向(径向),样品电导率随着滑动位移的增加快速增加,在滑动约数十厘米之后,其电导率基本达到稳定状态;在垂直滑动面方向(轴向),样品电导率基本不随摩擦滑动速率和滑动距离而变化.SEM显微结构观测显示,摩擦滑动所引起的电导率各向异性直接反映了石墨分布结构的变化.该研究结果深化了对地震断裂带浅部电性特征的认识,为野外断层带大地电磁测深资料的解释提供了约束,同时对于了解含碳断层的力学性质和弱矿物相在剪切变形中的分布特征及其演化过程等方面也具有重要意义.

  相似文献   

8.

为探究脆塑性转化带断层的力学性质和滑动稳定性,本文采用干燥的Carrara大理岩预切断层(saw-cut)样品,在气体介质三轴高温岩石力学实验仪上开展了摩擦实验研究,实验温度70~400℃,围压30~100 MPa,位移速率在0.08 μm·s-1,0.4 μm·s-1,2 μm·s-1之间切换.实验力学数据揭示,不同围压下Carrara大理岩断层摩擦系数随温度变化规律不同:低围压(30 MPa)下,摩擦系数随温度升高先增大后减小,中高围压(≥70 MPa)下摩擦系数则表现为随温度先减小后增大.断层摩擦滑动行为在100~300℃的范围内表现出由稳定的速度强化转化为不稳定的速度弱化,且在400℃左右重新转变为稳定的速度强化.实验后断层滑动面形貌和微观结构分析表明,稳定滑动断层面为高反射镜面,擦痕清晰;黏滑断层面为有光泽的凹凸不平的表面;最高围压下蠕滑的断层面粗糙无光泽,擦痕不可辨别.本文认为受温度激活的塑性变形过程逐步主导了岩石变形,对断层激活发生不稳定滑动至关重要,而高围压则会抑制断层的不稳定滑动.本研究结果不仅为识别野外断层黏滑和蠕滑提供了实验证据,而且为探究强震孕育和发生机理等提供重要的参数.

  相似文献   

9.
云龙  郭彦双  马瑾 《地震地质》2011,33(2):356-368
在实验室利用96通道应变记录采集系统和分布式多通道瞬态信号采集系统,观测了预切5.拐折断层的标本在变形失稳过程中应变场和声发射事件的时空演化.实验在双轴伺服加载系统上进行.在Y方向按位移控制方式加载,位移速率先后取0.5μm/s、1μm/s、0.5μm/s和0.1μm/s.观测得到:1)标本沿断层发生周期性的黏滑失稳,...  相似文献   

10.

共和盆地位于青藏高原东北缘,是我国重要的干热岩地热资源赋存区之一.成因机制研究是干热岩地热资源研究中最为基础与核心的工作之一,也是干热岩地热资源潜力精细评价和合理开发利用的重要依据.本文基于98块采自共和—贵德盆地岩石样品的放射性生热率数据,分析了盆地主要岩性岩石的放射性生热率分布特征,讨论了共和盆地干热岩地热资源的成因机制,并初步建立了干热岩地热资源的成因模式.研究表明,共和盆地恰卜恰地热区沉积岩(以泥岩和粉砂质泥岩为主)放射性生热率为1.21~2.02 μW·m-3,平均值为1.67±0.29 μW·m-3;以花岗岩和花岗岩闪长岩等为主的基底岩石的放射性生热率介于1.17~5.81 μW·m-3之间,平均值为3.20±1.07 μW·m-3.贵德盆地扎仓寺地热区沉积岩(以砂岩和泥质砂岩为主)放射性生热率为1.83~2.40 μW·m-3,平均值为2.13±0.23 μW·m-3;基底花岗质岩石的放射性生热率介于0.92~6.49 μW·m-3之间,平均值为2.81±1.40 μW·m-3.测试数据显示共和—贵德盆地基底花岗岩不存在高放射性生热率异常.但是,新生代以来印度—欧亚板块持续性陆-陆碰撞作用造成的壳内放射性生热元素富集层增厚,导致了花岗岩放射性生热率的热贡献量同步增大(30.3~40.5 mW·m-2),因此,花岗岩放射性生热可为共和盆地干热岩地热资源提供稳定的壳内热源基础.基于放射性生热率数据和热流配分研究,结合研究区已有地质-地球物理研究资料,本文认为壳内部分熔融层作为附加热源为共和盆地干热岩地热资源提供了重要的附加热流贡献.在此基础上,本文初步构建了共和盆地干热岩地热资源成因模式:加厚地壳花岗岩放射性生热与壳内部分熔融层供热.

  相似文献   

11.
本文针对深水环境下中深层偏移成像质量差的问题,考虑海水速度变化对中深层偏移成像质量的影响,从在大水深中加入深海声道模型入手,分析在偏移成像当中海水速度的不同选取对水平层状介质、倾斜层介质以及较复杂介质模型偏移成像质量的影响,通过在同一模型上改变海水速度进行成像,分析中深层成像效果可以得出:在大水深反射资料数据处理当中应该考虑真实的海水速度进行成像处理,否则会由于海水速度的选取不当而造成偏移成像层位的空间位置和中深层层位几何形态的变化.  相似文献   

12.
三门湾沿海声层析潮流观测实验   总被引:3,自引:0,他引:3       下载免费PDF全文
2009年9月6日至9日在三门湾进行了沿海声层析(Coastal Acoustic Tomography,CAT)潮流观测实验.实验由7台沿海声层析仪组网进行,并分别由渔船定点抛锚于7个站位.实验期间,还进行了定点ADCP(Acoustic Doppler Current Profiler)观测.通过建立逆模式对声传播时间差进行解析,引入权重因子,用L-curve法确定阻尼因子的最佳值,继而根据阻尼最小二乘法得到流速的最佳解.根据逆模式得到的流速分布可知该海区的潮流以半日潮(M2)为主,M2潮流椭圆呈东南-西北走向,潮流基本都是顺着水道,即涨潮为西北流向,退潮为东南流向.西北向与东南向最大流速分别为1.03 m·s-1和1.09 m·s-1.实验期间该区域的余流是从湾外流入湾内,平均流速约为0.05 m·s-1.CAT与定点ADCP流速的东分量和北分量的均方差均小于0.18 m·s-1.这样大面积的潮流和余流水平分布的同步观测,用传统观测手段很难实现.通过以上结果可以得出,沿海声层析技术可以作为一种新的测流方法对强潮海区进行大面积潮流观测,可在我国沿海的海洋环境监测等方面发挥重要作用.  相似文献   

13.
对大地震前的扰动现象的研究有助于认识地震孕育的动力学过程,震前重力扰动已成为关注的热点之一. 对大地震前的扰动现象的研究有助于认识地震孕育的动力学过程,震前重力扰动已成为关注的热点之一. 为检验日本Mw9.0级地震是否存在震前扰动现象,本研究利用全球超导重力仪记录到的地震前后7天内20组秒采样数据进行分析. 经潮汐、大气改正等处理去除仪器的漂移及残余潮汐效应,得到非潮汐重力变化曲线.结果表明大部分振幅大于30×10-8 m·s-2的曲线反映了全球Mw≥6级地震引起的高频波动信号,其中11组数据在3月9日Mw7.3级前震之前出现了扰动现象.震前扰动可分解为三个频段,其中,低于0.1 Hz和高于0.18 Hz的分量分别反映了地震波动信号及非构造信息,中间频段(0.118~0.18 Hz)信号能够较大程度地压制地震波动信号、并同时保留异常扰动信息.它的振幅在3月7日10时之前基本保持约1×10-8 m·s-2,之后开始逐渐增大,到3月9日7.3级前震前后达到最大,此后振荡衰减,振幅保持约(5~10)×10-8 m·s-2,直至主震发生.中间频段信号的变化特征与主震前的应力迁移过程以及实验记录到的地震成核过程有许多相似之处;不过,震前重力异常是否与主震前的应力加速积累有关,仍待进一步研究.  相似文献   

14.
机载SAR对海探测时, 探测范围小和时空匹配难等局限使其无法借助风条纹和辅助资料反演海面风矢量.本文在仿真研究CMOD5.N地球物理模型参数的函数关系, 实例分析机载SAR探测图像中距离向均值曲线变化规律的基础上, 发现相同风向、风速条件下, CMOD5.N模型构建的标准曲线和探测图像的距离向均值曲线遵循统一的归一化雷达截面随入射角变化规律, 且两者具有良好的相关性.据此, 本文提出将距离向均值曲线与标准曲线逐条匹配, 采用相关系数判定两者的相关程度, 选择使得相关系数绝对值最大的标准曲线作为最优匹配曲线, 进而直接确定风向和风速的海面风矢量反演方法.机载SAR飞行探测实验结果表明, 海面风矢量反演结果与浮标观测结果的均方根误差为风向11.3°, 风速0.9 m·s-1, 高于反演精度指标要求, 原因在于该方法既避免了机载SAR探测图像中斑点噪声的影响, 又不会产生局部最优解, 提高了海面风矢量反演精度.  相似文献   

15.
针对波速分层的区域岩体,在异向波速模型的基础上,对垂向上的应力波按岩体波速值大小作分段区别,推导震源应力波走时关系式,建立分层速度定位目标函数,基于此提出一种由参数准备、层速度反演、微震定位三个模块组成的分层速度定位模型SV,并采用遗传算法进行优化求解.然后,对分层速度定位模型在已构建微震监测系统的白鹤滩水电站左岸岩质边坡进行验证.微震事件重定位结果表明,分层速度定位模型定位微震事件的最大、最小和平均偏离层内错动带程度指标较单一速度模型分别降低了57.17%、36.51%和57.35%,证明了定位模型在波速分层的区域岩体微震定位应用中比单一速度定位模型更加合理可靠.  相似文献   

16.

通过布置于龙门山断裂带中段、龙门山山前玉皇观区域的地震观测台站阵列接收地震数据,研究该区域的地震动放大效应和地下地质结构.观测阵列共10台宽频带地震仪,分布在玉皇观河口冲积扇区域.分别采用参考场址谱比法(RSSR)和HV谱比法(HVSR)计算64个高信噪比近震数据的振幅谱比函数,结果显示在玉皇观区域具有较明显的地震动放大效应,并且局部场址效应显著.以S06场址为例,建立近地表地震地质模型,通过SH波放大效应正演模拟研究该场址的地震动放大模式.RSSR与HVSR的结果表明,两者所计算的场址放大效应主频一致,但是HVSR的放大峰值却比RSSR的放大峰值大一倍左右,表明HVSR的结果可能包含了波场在近地表低速层之下传播路径的改造作用.另外,采用27个远震P波的接收函数计算了该区域地壳上地幔S波速度结构.接收函数研究结果显示玉皇观地区的莫霍面深度为44 km,沉积盖层、结晶地壳和上地幔的S波速度分别为2.5 km·s-1、3.5 km·s-1和4.5 km·s-1.观测阵列台站之间的接收函数反演结果一致性较好,说明本研究区域范围内地形地貌等近地表结构因素的相对变化对接收函数的影响不大.

  相似文献   

17.
黏滑实验的震级评估和应力降分析   总被引:2,自引:0,他引:2       下载免费PDF全文
本文通过三种结构模型的黏滑地震模拟实验,利用高频速度连续观测系统获得了地震失稳过程的速度特征,讨论了最大位移量的选取方法,估算了实验室黏滑型地震的矩震级,探讨了黏滑类型、应力降大小与震级的关系.结果表明,黏滑型地震的应力降过程可能包含一次到多次高频振荡,对应若干次黏滑子事件.高频振荡的摆动幅度很大,包含有多种频率成分,峰值速度0.003~0.008 m·s-1.初步估计黏滑型地震的震级范围为-4.4~-3级,断层构造面的差异对各种黏滑模型的地震震级分布有明显影响.总体来看应力降与地震震级没有明显相关性,决定地震震级的主要因素应当是震源尺度.  相似文献   

18.
利用日本GPS网探测2011年Tohoku海啸引发的电离层扰动   总被引:2,自引:1,他引:1       下载免费PDF全文

海平面的海啸波会产生大气重力波进而引发电离层扰动.本文利用日本GPS总电子含量数据来探测2011年3月11日Tohoku海啸引发的电离层扰动.观测结果表明,在日本上空的电离层中存在两种重力波信号,分别由海平面的海啸波以及地震破裂过程产生.地震产生的电离层重力波分布在震中周围(包括海洋上空以及远离海洋的区域),而海啸引发的电离层重力波主要分布在海洋上空.地震产生的电离层重力波具有不同的水平速度,包括约210 m·s-1以及170 m·s-1,其频率为1.5 mHz;而海啸引发的电离层重力波水平速度快于前者,约为280 m·s-1,其频率为1.0 mHz.此外,海啸引发电离层重力波与海平面上的海啸波有相似的水平速度、方向、运行时间、波形以及频率等传播特征.本文的研究将电离层中的海啸信号与地震信号区分开来,进一步确认电离层对海啸波的敏感性.

  相似文献   

19.
基于镭同位素分布的黄海和东海垂直混合速率计算   总被引:4,自引:0,他引:4       下载免费PDF全文
在黄海和东海采样测定了水体中的镭同位素分布,用平流扩散模型描述镭同位素分布,最小二乘方法计算了垂直涡动扩散系数和上升流或下降流流速.结果给出北黄海中部、南黄海中部、浙江沿岸和台湾北部海域存在上升流,流速分别为0.46×10-3cm·s-1、0.17×10-3~1.39×10-3cm·s-1、2.02×10-3~3.04×10-3 cm·s-1和1.06×10-3~2.51×10-3 cm·s-1.北黄海中部和东海东北部存在下降流.流速分别为-2.30×10-3 cm·s-1和-0.61×10-3~-2.10×10-3 cm·s-1.计算同时给出的垂直涡动扩散系数为5.84~48.2 cm2·s-1,平均值为22.3 cm2·s-1.北黄海和浙江沿岸上升流流速与文献的结果一致;北黄海中部存在下降流与文献的结论一致.本研究结果与文献结果一致是对所建立的方法的肯定,也是对文献研究结果的支持.  相似文献   

20.
2015年9月17日6时54分32秒(北京时间)智利中部伊拉佩尔附近(震中31.57°S,71.67°W)发生了一次M_w8.3大地震,在此次地震震中以南约500 km处的马乌莱地区曾于2010年2月27日14时34分11秒发生过一次M_w8.8强震(震中36.12°S,72.90°W),两次地震余震分布区之间有约75 km的地震空区.本文利用远场体波与面波波形,基于有限断层模型,反演了这两次地震的震源破裂过程.结果显示这两次地震均为逆冲型大地震,2015年伊拉佩尔M_w8.3地震的平均滑动角度为107°,平均滑动量为2.43 m,平均破裂速度为1.82 km·s~(-1),标量地震矩为3.28×10~(21)Nm,95%的标量地震矩在104 s内得到了释放.最大滑动量约8 m,位于沿走向75 km,深度8 km处.2010年马乌莱M_w8.8地震的平均滑动角度为109°,平均滑动量为4.95 m,平均破裂速度1.90 km·s~(-1),标量地震矩为1.86×10~(22)Nm,95%的标量地震矩在121 s内得到了释放.最大滑动量约12.5 m,位于沿走向100 km,深度21 km处.2015年伊拉佩尔M_w8.3地震浅部更大的滑动量应该是其引起了较大海啸的一个原因.基于破裂滑动分布,我们计算了这两次地震引起的周边俯冲带上静态库仑应力变化,结果显示两次地震均显著增加了周边俯冲带上的库仑应力,2010年马乌莱地震使得2015.年伊拉佩尔地震震源区附近的库仑应力增加了(0.01~0.15)×10~5Pa,从应力积累的角度看,2010年马乌莱地震有利于2015年伊拉佩尔地震的发生,对后者的发生起到了促进作用.  相似文献   

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

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