首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 109 毫秒
1.
陈国兴  金丹丹  朱姣  李小军 《岩土力学》2015,36(6):1721-1736
针对基岩明显起伏、土层非均匀分布的典型河口盆地场地,考虑土体非线性特征,采用黏弹性人工边界模拟无限域对地震波动的影响,建立大尺度精细化二维有限元模型,分析了盆地地表地震动幅值、频谱、持时、传递函数特征,探讨了基岩起伏土层的地震动聚集效应及盆地边缘效应。结果表明:(1)盆地近地表土层表现出不同程度的地震动放大效应,且随土层深度增加呈非单调递减特征,基岩突变处地震动聚集效应明显,盆地两侧产生较为显著的边缘效应;场地中、长周期地震动的放大作用显著;(2)多遇地震、偶遇地震和罕遇地震水平时,场地卓越周期依次介于0.35~0.65 s、0.40~0.75 s和0.50~1.05 s之间;给出了盆地地表PGA(地表峰值加速度)、卓越周期均值等值线图及地表加速度反应谱放大因子建议值,地表设计地震动参数amax(地震影响系数)与Tg(特征周期)明显大于现行《建筑抗震设计规范》取值;(3)盆地特殊位置地表地震动持时得到不同幅度增长,且与输入地震动特性相关;(4)该盆地对0.5~2.0 Hz频段基岩地震动的放大效应比较显著,对小于0.2 Hz或大于2.5 Hz的基岩地震动,该盆地地震动放大效应不明显;(5)福州城区及其邻近区域地震动放大效应普遍较大。大尺度二维非线性分析一定程度上能合理反映微地形起伏、土层分布及土体非线性对地震波传播过程的影响。  相似文献   

2.
针对新型核电工程结构AP1000核岛结构设计地基中的5类非坚硬岩场地,即硬岩场地、软岩场地、上限软-中等土场地、软-中等土场地和软土场地,采用一维土层场地模型开展场地土和计算基底条件对设计地震动影响计算分析。分析中,场地模型的计算基底剪切波速分别取为700、1 100、2 438 m/s,计算基底输入地震动分别选择基于核电建设相关技术文件和规范规定的反应谱RG1.60谱、AP1000谱和HAD101/01谱(5个阻尼比)合成的人工地震动时程。计算分析表明:非坚硬岩场地会导致场地地震动峰值加速度及频谱特性显著变化,场地越软影响程度越显著;除软土场地外,场地对地震动峰值加速度和反应谱的影响均为放大作用,软土场地对地震动较低频段反应谱有放大作用,但对峰值加速度和较高频段反应谱具有强烈的减小作用;对于各类场地,计算基底及其剪切波速的变化均会导致地表地震动峰值及频谱特性明显甚至显著变化,其影响程度与计算基底剪切波速成正比;随着场地由硬变软,计算基底剪切波速的变化对场地地震动的影响程度大为减小,至软土场地几乎不产生影响。考虑到场地类型及计算基底选取对场地地震动的显著影响,我国核电厂建设引用AP1000标准设计时应合理分析场地的适宜性。  相似文献   

3.
远场大地震作用下大尺度深软场地的非线性地震效应分析   总被引:1,自引:0,他引:1  
战吉艳  陈国兴  刘建达  李小军 《岩土力学》2013,34(11):3229-3238
基于ABAQUS软件自行研制的并行计算显式算法集群平台,针对苏州城区典型地层剖面,建立了大尺度深软场地的二维精细化非线性有限元分析模型,对人工地震波和大地震远场地震动作用下深软场地的非线性地震效应进行了比较研究。研究结果表明:(1)与人工地震波作用时深软场地的地表峰值加速度放大效应相比,大地震远场地震波作用时的放大效应尤为显著,由于土介质的横向不均匀性及其非线性,使不同地表的峰值加速度放大效应存在明显的变异性。(2)深软场地对周期小于0.3 s的高频地震波均具有显著的滤波效应;大地震远场地震波作用时,深软场地对周期0.85~1.65 s的长周期地震波的放大效应非常显著,但对2.5~7.0 s的长周期地震波呈现出明显的滤波效应。(3)地震动峰值加速度PGA值沿土层深度和横向的分布形态呈现出明显的高低起伏现象,在不同成因的土层更迭面附近及土介质横向不均匀性显著的区域,地震波的局部聚焦放大和过滤减小现象尤为明显,且大地震远场地震动作用时,20 m以浅土层的PGA值呈现出非常显著的放大效应。(4)地震波的频谱特性、土层的横向不均匀性对深软场地地表加速度反应谱? 谱的谱形有显著影响;给出了描述加速度反应谱沿土层深度变化特征的三维谱形曲线,可以直观地展示出深软场地中细长地下结构地震反应可能存在类共振现象的土层深度。  相似文献   

4.
研究河谷场地地震效应对场地选址和抗震设计具有重要的指导意义。基于离心机振动台试验分析,研究了梯形河谷场地地震动响应规律。结果表明:基岩河谷场地对地震动有一定的放大效应,放大效应随着地形的变化而变化,但放大效应不显著,场地不同位置对反应谱的影响较小;基岩-土模型基岩面地震动放大倍数明显增大,不同输入地震动情况下放大倍数不同,各个场点对频段为0.5~2.5 s的地震动有明显的放大作用,对地震动放大频域范围明显加大,这与纯基岩场地的有明显的不同,虽然各个场点的反应谱形状有一定的差别,但是反应谱的平台值和特征周期相差不大;由于河谷场地地形效应,河谷场地地表峰值加速度随着地形的变化放大倍数随之变化,阶地级数越高放大倍数越大,谷底放大倍数最小,随着输入地震动强度的增加,阶地级数越高反应谱的平台值越高,特征周期越大。  相似文献   

5.
苏永奇  马巍  吴志坚  马尔曼 《冰川冻土》2016,38(4):1090-1098
根据青藏工程走廊北麓河及楚玛尔河场地地震危险性分析结果,合成年超越概率为1.97%、1.00%、0.21%、0.10%、0.04%、0.02%的人造基岩地震波作为输入地震动,结合场地钻孔剖面及波速资料,和已有的冻土动力学研究成果,建立一维模型,通过等效线性化方法进行场地地震反应分析计算,研究了青藏工程走廊多年冻土场地地震动加速度峰值特征及影响因素.研究结果表明,北麓河场地与楚玛尔河场地的人造基岩地震波峰值及持时均存在显著差异,北麓河场地峰值大、持时短,以近震影响为主,楚玛尔河场地峰值小、持时长,以中远震影响为主;多年冻土区场地,夏季场地地震动加速度峰值显著大于冬季,活动层融化对场地地震动加速度峰值有明显的放大效应;冬季场地冻结后,场地地震动加速度峰值随冻土波速增大而减小,最大减小幅度为6.1%,随动剪切模量比减小、阻尼比增大而减小,最大减小幅度为8.9%.活动层的融化有利于放大场地地震动加速度峰值,重大冻土工程抗震设防应予以重视.  相似文献   

6.
针对城市区域沉积盆地三维地震动模拟问题,发挥边界元法处理无限域波动问题的独特优势和快速多极算法(FMM)处理大规模问题的特性,形成适用于宽频、高自由度模型的快速多极子间接边界元法(FMM-IBEM)。通过精度对比和效率验证发现,FMM-IBEM可精确高效地实现三维复杂场地地震动宽频模拟,同时大幅度降低存储量和计算量。采用该方法对某含起伏基岩面的三维复杂形状沉积盆地进行频域和时域分析,结果表明:入射频率、盆地形状和盆地内部的位置点对沉积盆地地表位移放大作用均有不同程度的影响;总体上低频域的放大作用高于高频域;频谱曲线呈现多个共振频段;沉积盆地中部和边缘处地表位移放大效应差别显著。重要工程选址宜避开沉积盆地中地震能量积聚区域。  相似文献   

7.
采用有限元离散模型、黏弹性人工边界和中心差分积分格式相结合的方法,开展三维凹陷地形场地地震反应计算,分析了P波和SV波入射下均匀弹性半空间三维棱锥形凹陷、半球凹陷和棱柱形凹陷地形场地对场地地震动的影响特征与差异。研究表明:(1)不同形状的凹陷地形对场地地震动的影响差异显著,即棱柱凹陷地形>半球凹陷地形>棱锥凹陷地形,但不同形状凹陷地形影响之间的差异更多地表现在中频范围内;(2)凹陷地形对地震波具有复杂的散射作用,无论是垂直入射的P波还是垂直入射的SV波(仅竖向或水平向地震动输入),在凹陷区及附近均能产生两个方向很显著的地震动,且在SV波垂直入射下的凹陷散射更为强烈,其产生的竖向地震动幅值甚至超过水平向的;(3)凹陷地形可能会产生地震动的地形边缘效应,棱柱形凹陷边缘处地震动被显著地放大,而半球凹陷地形、棱锥形凹陷地形中这一效应相对较弱。研究揭示了凹陷边缘陡峭程度对场地地震动有显著的影响,而且凹陷地形边缘效应可能存在并与凹陷边缘陡峭程度密切相关。这一研究结果对实际工程问题分析给出了重要提示:建立凹陷地形的简化分析模型时需要特别关注凹陷边缘部位的合理处理。研究结果可应用于山区工程...  相似文献   

8.
多地貌单元复合场地非线性地震效应特征二维分析   总被引:2,自引:0,他引:2  
金丹丹  陈国兴  董菲蕃 《岩土力学》2014,35(6):1818-1824
针对由河漫滩、阶地、波状起伏的丘陵地带和残丘等组成的典型复合地形场地,基于ABAQUS软件的显式有限元并行计算集群平台,考虑土体的非线性特征,采用黏弹性人工边界模拟切除的无限域对地震波动的影响,建立了大尺度复合地形真实场地的二维有限元模型,并将基岩地震动视为SH波,分析了不同地貌单元复合场地地震效应特征。结果表明:地形差异引起地表各测点峰值加速度(PGA)及峰值位移(PGD)不同;在相同地震波入射下,河漫滩区地表PGA大于Ⅰ级阶地区域,但其PGD明显小于Ⅰ级阶地区域,且河漫滩区的峰值位移差动较小;场地地表水平加速度反应谱呈双峰及多峰现象,加速度在0.4~1.2 Hz频段放大较多,且加速度反应谱特性与输入地震动相关;特定频段(0.5~1.75 Hz)地震动放大效应与聚集效应比较明显,频率小于0.2 Hz及大于2 Hz时,场地地震动放大特征不显著。低丘区的敏感频率段宽于河漫滩区;观测点的位置影响相应的水平加速度持时大小,局部地形变化将引起地震动持时差异,且持时和输入地震动特性相关。  相似文献   

9.
高山峡谷复杂地形是我国西部地区常见的地形和场地条件,大量的工程(桥梁、大坝等)修建在这类场地上。实际震害调查结果表明不规则地形对地震动具有明显的放大作用,对边坡的稳定性和建筑物的安全性构成不利的影响。因此研究高山峡谷复杂地形的地震动放大效应具有重要的工程价值。本文针对浅切割的高山峡谷复杂地形(山体顶和峡谷底的高差在100~500 m范围内),基于边界积分方程法获得场地任意点的地震动,详细探讨了入射波类型(P波和SV波)、入射波频率、入射角、山体和峡谷的几何参数对浅切割的高山峡谷复杂地形的地震动放大效应的影响,发现山体的存在明显改变场地地震动的空间分布;高山峡谷复杂地形对地表地震动的放大作用与入射波频率密切相关;不同的地震波类型对应不同的地震动空间分布模式;峡谷深度对地震动放大效应的影响取决于入射波的类型;不对称的几何形状导致地震动的空间分布也呈现明显的不对称,并且SV波入射时的不对称性强于P波;当地震波斜入射时,峡谷背向震源侧的地震动远大于面向震源侧的地震动。本文的研究方法可以获得考虑高山峡谷复杂地形地震动放大效应的场地任意位置的地震动,为边坡的稳定性分析提供更符合实际的地震输入。  相似文献   

10.
钟紫蓝  史跃波  李锦强  赵密  杜修力 《岩土力学》2022,(7):2015-2024+2033
研究等效线性化方法中土体动力参数不确定性对场地地震反应的影响,提出了可考虑土体动力特征参数间相关性的土体动剪切模量和动阻尼曲线随机样本生成方法,基于Matlab开发一维等效线性化场地地震反应分析计算程序,用于开展场地地震反应随机动力分析。以Ⅱ类成层场地为例建立一维自由场分析模型,采用不同地震设防水准下露头基岩场地加速度反应谱为目标反应谱,合成人工地震记录折半后作为场地下卧基岩处的输入。计算结果表明,考虑土体动力特征参数不确定性对场地地震反应具有较大的影响,且影响程度与地震动强度、频谱成分及场地基本周期均密切相关;场地最大峰值加速度和地表加速度反应谱的波动范围随地震动强度等级的增加而变大,且最大峰值应变和最大峰值加速度的波动范围可达10%和14%;目标反应谱平台段及场地基本周期处对应的场地地表加速度反应谱的波动范围超过20%。  相似文献   

11.
陈国兴  李磊  丁杰发  赵凯 《岩土力学》2020,41(9):3056-3065
地震基岩深度和土体动力本构模型的选取对核岛场地地震效应评价结果的合理性具有重要影响。以拟建某沿海核电厂深度470 m沉积土夹火山岩层场地的3个钻孔剖面为研究对象,采用一维等效线性波传分析(ELA)法、基于Matasovic本构模型和Davidenkov-Chen-Zhao(DCZ)本构模型的一维非线性分析(NLA)法,选取不同剪切波速的5个岩土层作为地震基岩,研究了输入地震动特性、地震基岩深度和土体动力本构模型的选取对巨厚沉积土夹火山岩层场地非线性地震反应特性的影响。结果表明:(1)以浅层硬岩夹层或深部土层作为地震基岩,NLA法计算的5%阻尼比的地表谱加速度SA的短周期部分较之ELA法的计算值大,但两者计算的地表SA谱的长周期部分几乎一致;(2)基于Matasovic模型和DCZ模型的NLA法计算的地表SA谱谱形和峰值加速度随深度的变化趋势基本一致;(3)从NLA法计算的地表峰值加速度和累积绝对速度而言,以剪切波速约2 500 m/s的浅层硬岩夹层作为地震基岩是适宜的。  相似文献   

12.
考虑核电厂地基-基础的动力相互作用,应用显式动力有限差分法分析了地震作用下极软岩、较软岩、坚硬岩上核电厂建筑结构基础的地震响应特征,比较了岩石坚硬程度对基础加速度反应谱的影响。研究表明:随着岩石坚硬程度的提高,核电厂建筑物结构基础的地震响应有增加的趋势;在周期轴上,基础处的加速度反应谱曲线会随着岩石坚硬程度的提高逐渐向短周期(高频段)方向移动。在高频段,建造于较坚硬岩石上基础结构的加速度反应谱值偏大;在中等频段,建造于较软岩石上基础结构的加速度反应谱值偏大;在低频段,岩石坚硬程度对加速度反应谱的影响不显著。  相似文献   

13.
王志亮  阳栋 《岩土力学》2013,34(Z2):400-406
核电作为清洁型能源具有广阔的发展前景,研究核电站的地震响应具有重要意义。建立了AP1000核岛的有限元模型,进行模态分析获取核岛的振动形态。用无限元模拟远场地基,有限元模拟近场地基,并考虑基础埋置及土与结构相互作用,开展了不同地基条件下的核岛地震响应分析。结果表明:随着地基刚度的增加,核岛节点的响应加速度幅值增大,而加速度反应谱峰值对应的周期减小,核岛节点应力增大,但其相对位移减小;考虑土与结构相互作用之后,软弱地基的响应加速度峰值大幅减小,地基加速度反应谱峰值对应的周期增加,土与结构相互作用对此类地基上核岛的地震输入有较大影响;对于坚硬岩石地基,土与结构相互作用对地震输入的影响可以忽略。  相似文献   

14.
沙牌坝址基岩场地地震动输入参数研究   总被引:2,自引:0,他引:2  
钟菊芳  温世亿  胡晓 《岩土力学》2011,32(2):387-392
重大水利水电工程地震动输入参数必须根据专门的地震危险性分析结果来确定。目前由地震危险性分析得到的一致概率反应谱具有包络的意义,不能反映实际地震的频谱特性,输入“一致概率反应谱”可能导致地震作用偏大;拟合设计反应谱人工生成地震动加速度时程的频率非平稳性也没有得到很好解决。为了解决这些问题,得到与坝址地震危险性一致、具体地震的输入参数,结合沙牌大坝提出了一套适用于重大水利水电工程基岩场地地震动输入参数确定方法:通过以有效峰值加速度为参数的概率地震危险性计算分析,确定坝址不同超越概率下的有效峰值加速度及对坝址贡献最大的潜在震源区;在最大贡献潜在震源内利用震级空间联合分布概率最大法确定坝址设定地震,依据加速度反应谱衰减关系确定与坝址设定地震对应的设计反应谱;根据设定地震结果和时变功率谱模型参数衰减关系确定时变功率谱,将时变功率谱和最小相位谱按三角级数叠加法进行强度和频率非平稳地震加速度时程合成。在对沙牌坝址区域的地震活动性及地震构造环境分析评价的基础上,采用上述方法,得到了坝址基岩场地不同超越概率下的有效峰值加速度、设计反应谱、强度和频率非平稳地震加速度时程等地震动输入参数。  相似文献   

15.
通过对取自海南昌江核电厂一期工程主厂区的黑云母花岗岩,在RMT-150B多功能全自动刚性岩石伺服试验机进行单轴不同频率的循环加卸载试验,研究花岗岩动应变-应力滞回曲线、动弹性模量和阻尼比同循环频率之间的规律,然后对破坏断裂断口进行细观电镜扫描试验,分析了不同频率下花岗岩断口细观形貌与结构特征。结果表明,花岗岩在循环荷载作用下的加卸应力-应变曲线并不重合,而是形成一个封闭的滞回环;当循环频率在0.01~1.0 Hz范围内变化时,随着频率增大,花岗岩滞回环面积、动弹模和阻尼比都增大;通过细观数字图像试验分析发现,循环荷载作用下,随着频率的增大,岩石断口细观裂纹总的数目、面积、分形维数均变大。文中给出不同频率循环荷载作用下岩石破裂的细观损伤力学机制,研究成果对海南昌江核电厂的地震反应分析和场地安全性评价有着重要参考价值。  相似文献   

16.
This paper highlights the seismic microzonation carried out for a nuclear power plant site. Nuclear power plants are considered to be one of the most important and critical structures designed to withstand all natural disasters. Seismic microzonation is a process of demarcating a region into individual areas having different levels of various seismic hazards. This will help in identifying regions having high seismic hazard which is vital for engineering design and land-use planning. The main objective of this paper is to carry out the seismic microzonation of a nuclear power plant site situated in the east coast of South India, based on the spatial distribution of the hazard index value. The hazard index represents the consolidated effect of all major earthquake hazards and hazard influencing parameters. The present work will provide new directions for assessing the seismic hazards of new power plant sites in the country. Major seismic hazards considered for the evaluation of the hazard index are (1) intensity of ground shaking at bedrock, (2) site amplification, (3) liquefaction potential and (4) the predominant frequency of the earthquake motion at the surface. The intensity of ground shaking in terms of peak horizontal acceleration (PHA) was estimated for the study area using both deterministic and probabilistic approaches with logic tree methodology. The site characterization of the study area has been carried out using the multichannel analysis of surface waves test and available borehole data. One-dimensional ground response analysis was carried out at major locations within the study area for evaluating PHA and spectral accelerations at the ground surface. Based on the standard penetration test data, deterministic as well as probabilistic liquefaction hazard analysis has been carried out for the entire study area. Finally, all the major earthquake hazards estimated above, and other significant parameters representing local geology were integrated using the analytic hierarchy process and hazard index map for the study area was prepared. Maps showing the spatial variation of seismic hazards (intensity of ground shaking, liquefaction potential and predominant frequency) and hazard index are presented in this work.  相似文献   

17.
Nuclear power plants are designed to prevent the hazardous effects of the earthquakes and any external events to keep the safety of the plant. Ninety-one shallow seismic refraction profiles were performed to determine shear wave velocity of the engineering layers at the site of El Dabaa area that is situated to the northern coastline of Egypt for seismic hazard microzonation evaluation according to hazard index values. A microzonation is a procedure of delineating an area into individual zones having different ranks of numerous seismic hazards. This will aid in classifying areas of high seismic risk which is vigorous for industrial design of nuclear structures. The site response analysis requires the characterization of subsurface materials considering local subsurface profiles of the site. Site classification of the area under investigation was undertaken using P- and S-waves and available borehole data. The studied nuclear power plant site has been characterized as per NEHRP site classification using an average velocity of transverse wave (V s 30 ) of depth 30 m which acquired from seismic survey. This site was categorized into two site classes: the major one is “site class B,” and the minor one is “site class A.” The attenuation coefficient, the damping ratio and the liquefaction potential are geotechnical parameters which were derived from P- and S-waves, and have their major effects on the seismic hazard contribution. 1D ground response analysis was carried out in the places of seismic profiles inside the site for estimating the amount of ground quaking using peak ground acceleration (PGA), site amplification, predominant frequency and spectral accelerations on the surface of ground by the DEEPSOIL software package. Seven factors (criteria) deliberated to assess the earthquake hazard index map are: (1) the peak ground acceleration at the bedrock, (2) the amplification of the site, (3) the liquefaction potential, (4) the main frequency of the earthquake signal, (5) the average V s of the first 30 m from the ground surface, (6) the depth to the groundwater and (7) the depth to the bedrock. These features were exemplified in normalized maps after uniting them to 0–1 scores according to some criteria by the minimum and maximum values as linear scaling points. Multi-criteria evaluation is an application of multi-criteria decision analysis theory that used for developing a seismic hazard index map for a nuclear power plant site at El Dabaa area in ArcGIS 10.1 software. Two models of decision making were used in this work for seismic hazard microzonation. The analytic hierarchy process model was applied to conduct the relative weights of the criteria by pairwise comparison using Expert Choice Software. An earthquake hazard index map was combined using Weighted Linear Combination model of the raster weighted overlay tool of ArcGIS 10.1. The results indicated that most of the study site of the nuclear power plant is a region of low to moderate hazard; its values are ranging between 0.2 and 0.4.  相似文献   

18.
火山岩的密度和地震波传播速度一般高于相同深度的沉积围岩,与围岩具有较明显的密度和速度差,由于火山岩与围岩界面形态及产状复杂,导致地震波成像困难.笔者结合澧水火山岩盆地地震勘探试验工作,讨论了适合于火山岩盖层下寻找深部金属矿的地震勘探技术.采用小道间距、长排列、高覆盖次数、火山岩顶界面(或火山岩内部)激发的采集方式,有助于获取火山岩内部及其下伏地层界面产生的反射波.  相似文献   

19.
勘探实践表明,混合沉积逐渐成为渤海海域中深部勘探的理想领域.通过大量钻井资料统计,揭示了渤海海域混合沉积的成因条件,认为渤海海域混积岩的形成受气候与水体环境、物源、古地貌条件以及沉积样式的控制.模型正演表明不同厚度以及不同地层岩性组合的混积岩地震响应特征不同.实例表明,在水体气候环境、物源类型、古地貌特征等成因环境约束的条件下,通过单井分析、地震相分析、微古地貌恢复、正演分析、属性分析、叠前反演等多种技术综合的以“成因模式找规律,地震技术描分布”的混积型储层综合预测技术组合在渤海海域混合沉积储层预测中是可行的,其中不同厚度及岩性组合混积型储层采取不同的储层预测关键技术对渤海海域混积岩勘探有较好的指导意义.   相似文献   

20.
以29口钻孔的重矿物百分含量数据来分析样本,通过Q型聚类分析、因子分析方法恢复不同物源体系的发育范围及其母岩类型,并以不同物源体系的边界范围为限定条件,结合ZTR指数及其等值线图,对珠江口盆地惠州凹陷M层物源特征进行了系统的研究。研究结果表明:研究区发育3大物源沉积区,分别为东北-北部物源沉积区、西北物源沉积区和东南物源沉积区,且各区重矿物组合特征明显,主次物源分明,整体表现为母岩类型复杂的特征。东北-北部沉积物主要母岩为酸性岩浆岩,次要母岩为变质岩和沉积岩,西北沉积物主要母岩为中性及基性岩浆岩,次要母岩为酸性岩浆岩和变质岩,东南沉积物主要母岩为中性及基性岩浆岩,次要母岩为变质岩。东北-北部物源区和西北物源区为远源沉积的产物,而东南物源区为近源沉积的产物。东北-北部物源区沉积物来自北方和东北方向,而西北物源区和东南物源区沉积物分别来自西北方向和东南方向。  相似文献   

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

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