zhenbo

ISSN 0253-3782 CN 11-2021/P

2014年36卷5期

显示方式:          |     

线性单事件地震定位方法的改进及应用尝试
金明培, 汪荣江, 戴仕贵
2014, 36(5): 757-769. doi: 10.3969/j.issn.0253-3782.2014.05.001
[摘要](40) [HTML全文] (18)
摘要:
根据大震速报和快速预警实际需要,首先对Inglada线性单事件定位方法进行了适当的改进, 使其在仅有P波到时数据的情况下也能快速定位,且求解过程仅需简单迭代而不用奇异值分解; 其次, 尝试将改进后的方法从单层均匀模型引入到分层均匀模型中的近源台网定位情形,并通过单层均匀和分层均匀两种不同模型的理论实验讨论了该方法的可行性和适用范围; 最后整合了质量高且分布较好的距离2008年汶川MS8.0地震震中最近的强震、微震, 以及川西流动台阵等观测记录资料,对汶川MS8.0主震初始破裂点的时空参数进行了多种模型的定位实验. 结果表明, 改进后的线性单事件定位方法简单、快捷、易用,可广泛应用于近源地震定位,尤其是用于无法得到S波到时的中强以上直至巨大地震的速报、地震现场流动台网的快速定位以及地震的快速预警等.
2013年芦山MS7.0地震前后姑咱台四分量钻孔应变时频特征分析
刘琦, 张晶, 池顺良, 闫伟
2014, 36(5): 770-779. doi: 10.3969/j.issn.0253-3782.2014.05.002
[摘要](43) [HTML全文] (15)
摘要:
2013年4月20日芦山MS7.0地震前, 四川省姑咱台四分量钻孔应变仪再度记录到类似于2008年汶川MS8.0地震前的高频应变畸变信号. 为了对该畸变信号各频率成分的时频演化特征进行细致分析, 讨论应变异常的可信度及其与芦山地震的相关性, 本文利用S变换方法对应变数据进行了处理. 结果表明, 研究时段内时频域中共出现两簇高能量异常, 一簇开始于2012年10月并持续了约4个月, 另一簇则开始于芦山MS7.0地震前数天. 同时段远震的影响及施工干扰并不能构成该异常的主体, 因此该异常可信度较高. 异常中包含了周期由数分钟至数小时的信号成分, 其中较短周期的信号成分震前出现较晚, 震后衰减较快. 震后异常信号的整体能量开始逐步衰减, 异常信号成分减少, 至2013年8月末应变观测记录已基本恢复至正常状态, 整个应变异常显示的时频演化过程与芦山MS7.0地震具有较高相关性.
利用图像信息方法研究芦山MS7.0地震前川滇及附近地区的图像异常
张小涛, 张永仙, 夏彩韵, 吴永加, 余怀忠
2014, 36(5): 780-789. doi: 10.3969/j.issn.0253-3782.2014.05.003
[摘要](42) [HTML全文] (17)
摘要:
图像信息方法是一种基于统计物理学的地震预测方法. 该方法通过严格的统计检验给出地震活动状态显著偏离平均状态的地区(即显著平静或显著活跃的地区), 称为“地震热点”. 以往研究结果表明该方法是一个预测效能较高的中长期地震预测方法. 本文利用中国地震台网中心提供的1970年以来的地震目录, 在剔除余震、 选取完备震级的基础上, 基于前期研究工作结果中预测效能较高的计算参数, 即网格尺度为2°×2°、 预测时间窗为9年, 分析2008年汶川MS8.0地震后到2013年芦山MS7.0地震前川滇地区的图像信息演化过程. 结果表明, 2009—2010年间芦山MS7.0地震震中附近的鲜水河断裂带、 龙门山断裂带和安宁河断裂带交汇区(29°—31°N, 101°—105°E)存在明显的地震热点, 2011—2012年该地震热点缩小至(29°—31°N, 101°—103°E)范围, 2012年以后消失, 继而于2013年发生芦山MS7.0地震. 本文还对比研究了网格尺度为1°×1°时的图像信息演化过程, 结果显示2009—2010年芦山地震附近的新津—成都—德阳断裂南部、 龙泉山西缘断裂南段、 马边—盐津断裂带、 峨边断裂带、 天全—萦经断裂带(29°—30°N, 103°—104°E)也出现过地震热点, 2011年以后该热点消失, 2013年发生芦山MS7.0地震. 本文讨论了计算参数对预测结果的影响, 并根据目前存在的地震热点分析了川滇及附近地区今后强震发生的可能地点.
2012年11月20日宁夏永宁MS4.6 地震震源机制解与发震构造
曾宪伟, 谢祖军, 莘海亮
2014, 36(5): 790-799. doi: 10.3969/j.issn.0253-3782.2014.05.004
[摘要](39) [HTML全文] (17)
摘要:
2012年11月20日在宁夏银川市永宁县与兴庆区交界处发生MS4.6地震,为了更好地了解此次地震的发震构造,首先采用Hypo2000绝对定位方法得到该地震的震中位置及余震分布;然后采用CAP方法反演了此次地震的震源机制解和震源深度. 反演结果表明,永宁MS4.6地震是一个带有少量逆冲分量的右旋走滑地震.该地震矩震级为MW4.3,最佳双力偶解为:节面Ⅰ走向11°,倾角74°,滑动角171°;节面Ⅱ走向103°,倾角81°,滑动角16°.最佳震源深度为8km左右.从该地震震中和震源机制解以及震源深度剖面分布来看,这次地震很可能发生在银川隐伏主断层西侧的次级断层上.
基于W震相技术的全球强震(MW≥6.5)矩心矩张量自动反演系统评估
赵旭, DuputelZacharie, 黄志斌, 陈宏峰, 赵博
2014, 36(5): 800-809. doi: 10.3969/j.issn.0253-3782.2014.05.005
[摘要](49) [HTML全文] (17)
摘要:
使用最近构建的虚拟全球地震台网记录的波形数据,采用OpenMP的并行编程技术优化原有算法, 研制了基于W震相技术的全球强震(MW≥6.5)矩心矩张量自动反演系统.为了评估该系统的准确度和时效性,将离线自动测定的2008年1月—2013年7月全球140次地震(MW6.5—9.0)的矩心矩张量与全球矩心矩张量工作组(GCMT)的结果进行了比较.结果表明:该系统可准确测定全球MW≥6.5地震的矩张量,绝大多数地震矩震级与GCMT给出的矩震级呈现出良好的线性趋势,两者之差ΔMW的标准方差约0.13,ΔMW位于区间(-0.2,0.2)的地震占总数的96%;地震矩6个分量分别与GCMT相应的结果沿对角线近线性分布,多数地震矩心水平位置与GCMT给出的矩心水平位置比较接近,两者间大圆弧距离位于区间(0,50km)的地震次数占总数的84%;在台站覆盖较均匀的条件下,该系统能够实现震后25—40分钟自动准确测定全球MW≥6.5地震的矩心矩张量.
大别-苏鲁及其邻近地区基于背景 噪声的勒夫波群速度成像
叶庆东, 丁志峰, 郑晨, 吕苗苗, 吴萍萍, 陈浩朋
2014, 36(5): 810-825. doi: 10.3969/j.issn.0253-3782.2014.05.006
[摘要](35) [HTML全文] (18)
摘要:
大别-苏鲁造山带是中国大陆东部地区最重要的构造之一. 为了研究该地区的地壳上地幔速度结构,本文收集了国家数字地震台网和中国区域地震台网的山东、河南、安徽、江苏和湖北等省的144个宽频带地震台连续两年(2009年5月—2011年5月)的水平向地震记录(E分量和N分量)数据,首先对台站对之间E-E,E-N,N-N和N-E分量进行互相关,然后分别对这4个互相关分量采用相位权重叠加法进行叠加,最后旋转到横向分量(T-T)获得勒夫波经验格林函数(EGF);用频时分析(FTAN)方法获得4000余条勒夫波群速度频散曲线,并进一步反演得到了周期为6—40s的勒夫波群速度分布图.结果表明,周期为6—10s的勒夫波群速度分布与地表构造特征相吻合.大别造山带、苏鲁造山带、湖北西部隆起均表现为高速;华北盆地发育,表现为大面积的低速;江汉盆地、南襄盆地、合肥盆地等因其规模不同而显示不同程度的低速.在周期为6—30s的勒夫波群速度分布图上,大别和苏鲁地区均显示高速,已有的研究结果中地壳的低速并没有得到反映.其原因一方面可能与勒夫波群速度纵向分辨不高有关,另一方面高压变质岩深度分布可能比已有研究结果给出的要深. 6—30 s的分布图上郯庐断裂带及其邻近地区表现为不同程度的高速,可能与该地区白垩纪以来处于拉张构造体制,地幔物质受到扰动,造成物质上涌有关.
广东及其邻域噪声面波层析成像
沈玉松, 康英
2014, 36(5): 826-836. doi: 10.3969/j.issn.0253-3782.2014.05.007
[摘要](47) [HTML全文] (25)
摘要:
通过收集广东及其邻域104个固定地震台站近10个月的垂直分量连续波形数据资料, 使用地震背景噪声互相关格林函数方法, 获得了大部分台站对的背景噪声互相关曲线. 基于这些对称叠加的互相关曲线, 利用时频分析方法, 进一步提取了该地区周期为5—40 s的基阶瑞雷波群速度频散曲线. 其噪声来源分析结果显示: 广东及其邻域的噪声场来源有很强的方向性, 短周期(5—10 s)噪声主要来自东南方向, 范围基本与海岸线分布一致, 可能是由于近海水陆相互作用产生的;较长周期 (15—30 s) 噪声主要来自三大洋的方位. 以这些提取的噪声面波资料为基础, 采用噪声面波层析成像方法反演得到了该地区周期为5—28 s的瑞雷波群速度层析成像图, 从该图可以看出, 广东及其邻域地下结构的横向变化总体较小, 沉积层厚度较薄, 地壳中可能普遍存在一个低速层;从研究区历史地震的分布及其表层地质构造的发育特征来看, 地震主要分布在高、低速过渡带附近, 表明面波群速度与地震之间具有较强的耦合关系;从群速度的低速异常特征来看, 广东及其邻域普遍分布的温泉和高地热主要受深部构造的控制和影响.
利用接收函数研究渤海湾盆地沉积层结构
武岩, 丁志峰, 朱露培
2014, 36(5): 837-849. doi: 10.3969/j.issn.0253-3782.2014.05.008
[摘要](42) [HTML全文] (20)
摘要:
对布设于渤海湾盆地地区的70个台站2006年9月—2009年9月记录到的895个远震数据进行接收函数计算,发现位于沉积层上的台站记录的接收函数震相较为复杂,由沉积层基底的转换波、多次波及莫霍面转换波、多次波共同构成,且随着沉积层厚度的增加,沉积层基底转换波、多次波与莫霍面转换波、多次波震相相互叠加,使之无法分辨出莫霍面转换波及多次波到时.为此,本文利用相邻算法对沉积层模型进行搜索, 以获取与实际波形误差最小的波形相对应的模型,即台站下方的沉积层结构.结果显示:① 渤海湾盆地由NNE向展布的坳陷带与隆起带组成,冀中坳陷带的沉积层厚度为3—6 km, 坳陷带内沿NE-NNE向分布着带状凹陷带与凸起带相间的次级构造,而且其东侧靠近沧县隆起区的凹陷带沉积层基底埋深最大,体现了渤海湾盆地基底受正断层控制而形成的半地堑(地堑)-半地垒(地垒)构造;② 隆起带的波速比vP/vS大于坳陷带的波速比,推测与隆起带下第三系地层的缺失有关, 坳陷带内高波速比地区与地热田的分布相对应, 说明沉积层内高温环境下波速比较高;③ 隆起带内的S波平均速度较坳陷带内的S波平均速度小,且冀中坳陷带内沉积层较厚地区对应的S波平均速度比其它地区更大, 说明渤海湾地区沉积层厚度与沉积层内的S波平均速度成反比关系,推测与下第三系地层的厚度有关,且下第三系地层越厚的地区,S波平均速度越大.以上结果表明,利用接收函数研究的结果与通过地质、地热等方面的研究结果存在着良好的对应关系.
用接收函数方法研究中国境内地壳结构
张攀, 朱良保, 陈浩朋, 王清东, 杨颖航
2014, 36(5): 850-861. doi: 10.3969/j.issn.0253-3782.2014.05.009
[摘要](42) [HTML全文] (22)
摘要:
利用中国数字地震台网30个台站的高质量宽频带远震数据,采用H-k叠加搜索法对中国境内的地壳结构进行研究,获得了研究区内的地壳厚度和vP/vS分布特征.结果表明, 中国境内的vP/vS值介于1.6—1.9之间,地壳厚度变化剧烈,在29—81 km之间.100°—110°E之间存在一个地壳厚度陡变带, 将中国分为东西两个部分.东部地壳厚度相对均匀,为31—36 km, 西部地区地壳厚度相对较厚且变化较大,中部地区地壳厚度为34—49 km.总的看来,青藏高原地区地壳最厚,可达81 km;天山、准噶尔盆地和内蒙古地区地壳厚度次之;华南地区地壳最薄.另外,中国大陆地壳平均波速比为1.738(σ=0.253),比全球大陆平均波速比1.78(σ=0.269)低.较低的波速比可能暗示中国境内地壳低速层的存在或者铁镁质成分的缺失.
南极长城站地壳和上地幔各向异性分析
刘莎, 田宝峰, 姜旭东, 徐志强, 郑钰, 杨建思
2014, 36(5): 862-871. doi: 10.3969/j.issn.0253-3782.2014.05.010
[摘要](42) [HTML全文] (20)
摘要:
利用我国第24次和第25次南极科学考察队于2008年2月—2010年3月南极长城站记录到的地震事件数据进行剪切波分裂研究. 选取近震事件对Sg波进行剪切波分裂计算,结果表明快波偏振方向有两个,分别为北东向和近南北向; 慢波延迟时间的范围为1.45—5.17 ms/km,平均值为3.54 ms/km.同时选取长城站记录到的远震数据SKS波震相进行剪切波分裂计算,得出上地幔快波偏振方向优势取向为北东向, 慢波延迟时间平均值为1.60 s. 剪切波分裂结果显示长城站地区地壳和上地幔具有明显的各向异性, 并显示长城站地区地壳与上地幔快波偏振方向几乎平行,表明壳幔变形的一致关系.另外,地壳和上地幔各向异性的快波偏振方向不仅与长城站附近的海沟方向平行,同时也与绝对板块的运动方向平行.该结果进一步说明了绝对板块的运动是构成上地幔各向异性的主要原因.
基于随温度变化的热系数模拟板块俯冲动力学过程
张晨, 张双喜
2014, 36(5): 872-882. doi: 10.3969/j.issn.0253-3782.2014.05.011
[摘要](40) [HTML全文] (19)
摘要:
热传导系数和热膨胀系数是影响板块俯冲动力学过程的两个重要参数. 由于地球介质的不均匀性,热系数也会随深度发生变化.然而,这种变化在地球动力学模拟研究中往往被忽略.本文针对随温度变化的热传导系数和热膨胀系数, 模拟板块俯冲的动力学过程,分析热系数、黏度对板块俯冲形态的影响及其对应的地幔对流特征.结果表明,依温度变化的热传导系数和热膨胀系数会影响地幔温度及黏度分布,进而改变板块的俯冲角度;黏度是控制板块俯冲动力学演化过程的重要因素;地幔对流受黏度结构的影响,呈现分层对流及局部多个对流环等多种不同形态的对流场特征.
凸起地形对地震动特性的影响
郝明辉, 张郁山
2014, 36(5): 883-894. doi: 10.3969/j.issn.0253-3782.2014.05.012
[摘要](42) [HTML全文] (19)
摘要:
采用基于ABAQUS平台的显式有限元动力学分析方法,结合人工黏弹性边界理论,研究了局部凸起地形对地震动特性(包括反应谱、峰值加速度、峰值速度和峰值位移等)的影响,分析了台地宽度对地形放大效应的影响.结果表明:凸起地形平台段空间点地震动受地形效应影响较大,在零阻尼条件下,其谱比曲线呈双峰特点,.8—0.9s的中长周期段谱比达到一个较大值1.6,在0.08—0.09s的高频段谱比超过2.0,且最大值出现在平台中点;对于凸起地形斜坡段,在大部分周期点处,顶点的谱比高于其它斜坡点,而且在周期超过0.4s的频段,斜坡段观测点的谱比表现出较明显的规律性,即越靠近顶点的观测点,其谱比值越大;坡底段地表不同观测点的谱比基本在脚点与计算边界点(人工边界点)对应的谱比值之间变化,在不同的频段均表现出较明显的规律性.凸起地形平台段宽度对地震动高频成分的放大效应具有较大影响,但只局限在一定的宽度范围内,随着宽度的增大,其对地形放大效应的影响逐渐减弱. 此外,台地宽度的变化对地震动峰值加速度的放大效应有一定的影响,而对峰值速度、峰值位移的放大效应的影响则不明显.
多震相菲涅耳体射线走时同时反演成像
李兴旺, 白超英, 黄国娇
2014, 36(5): 895-909. doi: 10.3969/j.issn.0253-3782.2014.05.013
[摘要](38) [HTML全文] (8)
摘要:
高频假设下的地震射线理论以及相应的地震成像理论表明,在射线稀疏条件下,不可能得到较高分辨率的构造成像;而有限频射线理论更符合实际地震的传播规律,即地震波的走时不仅与中心射线(传统的几何射线)上的速度分布有关,而且与中心射线附近一定范围(称其为第一菲涅耳体)内的速度异常分布有关.鉴于此,本文提出了计算多震相地震波菲涅耳体有限频射线的方法,并定义了走时敏感核函数,同时给出了利用多震相菲涅耳体有限频射线进行速度模型和反射界面同时反演成像的公式.利用多震相走时资料,使用传统射线层析成像方法与有限频射线层析成像方法进行了速度和界面的同时反演成像.结果表明,当射线密度较小时,无论是对速度模型的重建还是对反射界面几何形状的更新,有限频射线层析成像方法均优于传统射线层析成像方法, 而变频有限频射线层析成像则是实际地震层析成像的首选反演算法.
对朝鲜2006年、2009年和2013年3次地下核试验的相对定位
潘常周, 靳平, 徐雄, 王红春, 肖卫国, 沈旭峰, 张诚鎏, 徐恒垒
2014, 36(5): 910-918. doi: 10.3969/j.issn.0253-3782.2014.05.014
[摘要](41) [HTML全文] (18)
摘要:
采用双差定位方法, 对朝鲜2006年、2009年和2013年的3次核试验进行相对定位.根据19个区域地震台和8个远震台的垂直分量波形资料, 计算了多种震相(P波、S波和LR波等)的互相关时差, 并由测量时差反演两两核爆的相对空间位置, 通过随机生成部分台站组合进行扰动实验以估算定位的不确定度.结果表明, 相对于朝鲜2009年核试验, 2013年核试验位于其南257 m、西385 m, 2006年核试验位于其南503 m、东2589 m.当扰动实验中随机选取的台站数目为20个时, 在95%置信度下, 朝鲜2013年和2006年核试验相对于2009年核试验的地表定位结果的不确定度约为150 m.
内蒙古区域地震台网监测能力研究
刘芳, 蒋长胜, 张帆, 杨彦明, 梁莹, 王磊, 苗春兰
2014, 36(5): 919-929. doi: 10.3969/j.issn.0253-3782.2014.05.015
[摘要](43) [HTML全文] (13)
摘要:
国际上新近发展的“基于概率的完整性震级”(PMC)方法,具有可考察地震定位中由于台站人为选择等造成的台网监测能力下降,以及避免传统基于G-R关系的统计算法因地震数目过少而无法评估等优点.本研究利用PMC方法,计算得到内蒙古区域地震台网39个台站对周边地震事件的检测概率及台网检测概率.单台检测概率结果显示:PMC方法能够客观地反映39个台站对地震事件的检测能力;因台网布局等影响,内蒙古区域地震台网中西部和中东部地区的台站检测能力较强,而靠近蒙古、俄罗斯边境的台站, 阿拉善右旗附近地区的台站,以及邻近吉林、黑龙江等地区的台站检测能力较低.合成检测概率结果显示,由于邻省台站的引入,全区80%的地区基于概率的最小完整性震级MP达到2.2左右,其余地区MP达到3.3左右.为提高地震台网监测能力,建议在监测能力较弱的中蒙交界地区、东北部地区,以及阿拉善左旗以西地区适度加密台站,进一步优化台网布局.
基于Copula理论挖掘地磁Z分量震磁信息
张明东, 张文蕾, 曹井泉, 王建国, 谭青, 张玮
2014, 36(5): 930-943. doi: 10.3969/j.issn.0253-3782.2014.05.016
[摘要](34) [HTML全文] (21)
摘要:
利用Copula理论,选取天津地区5个台站的地磁资料,分析计算了1986—2013年的地磁Z值数据,经过在多种模型中比较Kendall与Spearman系数,得出最优模型,并分析不同周期下Kendall系数的时序曲线. 计算结果表明,在MS≥5.0地震前,Kendall系数的时序曲线处于紊乱期,震后其对应性迅速转好,一致性明显. 这是由于地震应力积累导致地下电性结构发生变化,而震后应力得到释放,磁场回归常态造成的. 尤为明显的是2006年文安MS5.1地震前,青光台(距震中最近)Kendall系数异于其它台站出现大幅下滑; 而2012年唐山MS4.8地震前,宁河台(距震中最近)Kendall系数异于其它台站上升幅度不大.
有限元算法在声波方程数值模拟中的频散分析
印兴耀, 周建科, 吴国忱, 梁锴
2014, 36(5): 944-1898. doi: 10.3969/j.issn.0253-3782.2014.05.017
[摘要](43) [HTML全文] (17)
摘要:
针对有限元算法在地震波数值模拟中的数值频散问题,利用集中质量矩阵双线性插值有限元算法,推导了二维声波方程的频散函数.在此基础上采用定量分析方法,对比分析了网格纵横长度比变化时的入射方向、空间采样间隔、地震波频率以及地层速度对数值频散的影响.数值算例和模型正演结果表明:当采用集中质量矩阵双线性插值有限元算法时,为了有效地压制数值频散,在所使用震源子波的峰值频率对应的波长内,采样点数目应不少于20个;减小网格长度的纵横比可以有效地抑制入射角(波传播方向与z轴的夹角)较小的地震波的数值频散;地震波频率越高,传播速度越慢,频散越严重,尤其是当相速度与其所对应的频率比值小于2倍空间采样间隔时,不仅会出现严重的数值频散,还会出现假频现象.
FFT和DFT在地震数据采集器幅频特性测试中的应用
李彩华, 滕云田, 张旸, 胡星星
2014, 36(5): 956-963. doi: 10.3969/j.issn.0253-3782.2014.05.018
[摘要](40) [HTML全文] (17)
摘要:
在进行地震数据采集器的幅频特性测试中,当标定信号频率稍高时(约为采集器采样率1/5以上),采集器输出的正弦波形发生畸变,无法人工读取正弦波峰值.为了解决该问题,本文提出了基于离散傅里叶变换(DFT)的积分恢复算法,将该方法与快速傅里叶变化(FFT)谱分析方法同时应用于一台高精度地震数据采集器的幅频特性测试数据处理.结果表明这两种算法均能标准化、高精度地完成采集器幅频特性标定数据的计算.其中基于DFT的积分恢复算法的计算程序操作简单,能够高效率地完成数据采集器的幅频特性检定工作;并且该方法还可以应用于地震数据采集器的线性度测试.
几种自由边界实施方法在完全 匹配层条件下的对比研究
黄建平, 杨宇, 李振春, 田坤
2014, 36(5): 964-977. doi: 10.3969/j.issn.0253-3782.2014.05.019
[摘要](33) [HTML全文] (17)
摘要:
传统的完全匹配层技术是一种能够较为有效地消除边界反射的边界条件,但是当表层为泊松比较高的自由表面时,该技术可能会产生不稳定的现象.针对传统的完全匹配层技术固有的不稳定和掠射情况下吸收效果不佳等缺陷,发展了多轴完全匹配层、卷积完全匹配层以及将两者结合的多轴卷积完全匹配层等3种边界条件.本文介绍了水平自由表面的不同处理方法以及传统、多轴、卷积和多轴卷积等4种完全匹配层条件的原理,通过二维半无限空间模型的交错网格有限差分正演模拟对比,分析了几种自由边界实施方法在这几种完全匹配层条件下的稳定性,并通过提取单道波形与解析解进行对比,定性分析了水平自由表面几种不同处理方法的准确性以及各自的适用条件. 结果表明,泊松比和水平自由表面实施方法对波场模拟效果及其稳定性有重要影响.
全球大震和中国及邻区中强震地震活动(2014年5—6月)
陈培善
2014, 36(5): 978-979. doi: 10.3969/j.issn.0253-3782.2014.05.020
[摘要](43) [HTML全文] (13)
摘要:
本系统由北京仁和汇智信息技术有限公司设计开发 技术支持: info@rhhz.net 百度统计