地震前后重力场变化的典型震例分析

第22卷第4期东北地震研究Vol.22No.4 2006年12月SEISMOLOGICAL RESEARCH OF NORTHEAST CHINA Dec.2006地震前后重力场变化的典型震例分析李清林!李天生!秦建增!谢汝一(中国地震局地球物理勘探中心9河南郑州450002)摘要!系统总结了中国地震局地球物理勘探中心利用流动重力资料进行地震预报的实践经验9研究了8次地震的孕震环境\重力场变化\成因分析和预报实况等0结果表明9一次较强的地震前9处于不同地段的重力测点和测段都要出现上升-下降-上升或下降-上升-下降的较明显的中\短期异常变化9重力场的变化图象较清晰地反映了地震前后重力场的演化过程0关键词!重力场9重力场变化9震例中图分类号!P315.7P312.7 文献标识码!A文章编号!1004-2865"2006#04-0014-110引言中国地震局地球物理勘探中心是我国地震系统较早开展流动重力测量和地震实验预报的单位之一9自上世纪70年代开始从事该项工作以来9积累了较丰富的观测资料和多次成功的震例0先后在华北地区\南北地震带\西北地区\西南地区\连云港-郑州-西宁\京津唐张渤地区\西藏自治区\内蒙古呼和浩特-包头地区\四川盐源和松潘地区\云南龙陵地区\山西地区和冀豫鲁地区开展了流动重力\流动地磁监测工作0本文较系统地整理了30多年来8次地震前后重力场变化的典型震例0震级选取以中国地震局分析预报中心地震月报目录为准01龙陵!S7.4主震后几次强余震1.1地震活动概况1976年5月29日在云南省龙陵地区相继发生了7.3级\7.4级强烈地震07月21日\8月2日\8月4日9在老震区又相继发生了一次6.6级\两次4~5级强余震01.2地质构造概况龙陵-潞西的东南隅为局部重力负异常区9该异常区是由怒江\龙陵\瑞丽及畹町三条大断裂所围限的花岗岩体展布的区域0龙陵测点位于龙陵-瑞丽大断裂带上9潞西测点位于该断裂带的东南侧0收稿日期!2006.0l作者简介!李清林(1951-)9男9高级工程师.主要研究方向为地震监测预报\综合地球物理勘探.李清林等:地震前后重力场变化的典型震例分析!期图2盐源地震前后重力段差值变化曲线fig.2The curve of the gravity difference vaiue before and after Yanyuan earthguake1.3重力场变化7.3级\7.4级主震发生后9地球物理勘探中心立即按中国地震局的安排9奔赴震区进行重磁复测和人工地震测深工作O 剖面跨越震中区9穿过了局部布格重力异常区和异常梯级带0 6.6级强余震震中距潞西4Okm 9距龙陵15km 0为了监测强余震前后的重力场变化9在该剖面上的大龙-龙陵\龙陵-J 863\J 863-潞西\潞西-遮放四个测段上进行了短周期重力观测0图1为J 863-潞西测段强余震前后的重力段差值随时间的变化曲线0可以看出9在1976年7月21日6.6级和8月2日5.O 级强余震之前短短几天时间内9该测段重力段差值出现了明显的震前上升-下降的变化过程9变化量达到1OOx1Onm ~S -20图1龙陵强余震前后重力段差值变化曲线fig.1The curve of the gravity difference vaiue before and after Longiing Strong afterShock2盐源!S 6.7\!S 6.4地震2.1地震活动概况1976年11月7日和12月13日9相继在四川省盐源县发生了6.7级和6.4级两次较强地震0震前布设的流动重力复测剖面距震中较近9最近的盐源测点距震中盐塘区仅45km 9西昌测点距震中132km(程式等91986)O2.2地质与地球物理概况该区的安宁河深大断裂与南北向密集的布格重力异常梯级带相对应9且多处出现拐弯\扭曲\挤压和局部的正负相间的重力异常O 老凉山地块内部整体性较好9地震活动相对平静9故将昭觉测点看做正常场O2.3重力场变化1976年5月用两台ZSM-!型和一台ZSM-"型重力仪进行复测O 1976年5~7月的观测结果与1972年建点基值相比发现:西昌-龙恩日达测段段差值上升193x1Onm ~S -29龙恩日达-昭觉测段段差值上升55x1Onm ~S -29其它测段未出现显著变化9见图2O 从1976年6月起对两个复测剖面(共12个测段)9!"!!东北地震研究!!卷进行了加密观测0取得了较完整的盐源地震前后重力场变化资料0可以看出9沿重力梯度带南北向剖面的变化不大9而横切重力梯度带的东西向剖面出现较大异常9且在1976年8月~9月变化最为剧烈02.4预报实况根据上述重力场变化情况9结合该区流动地磁和定点地磁观测资料9对盐源地震做出了较好的短期预报9受到四川省地震局和地方地震局的表彰和奖励03北京高丽营!S 4.0地震3.1地震活动概况1978年10月3日9在北京市高丽营地区发生了!"4.0地震03.2地质构造和地球物理概况高丽营 小汤山测点位于太行山山前断裂带中的高丽营超壳断裂与南口-孙河第四纪活动性断裂交汇部位附近0顺义和高丽营测点同处于局部重力异常区内9由地壳深部活动引起的重力变化是同步的9因此9反映在两点间的重力段差值变化亦小0而高丽营和小汤山测点则分别位于局部重力异常区的边缘及重力异常梯级带上9由地壳深部活动引起的重力变化会有明显差异9因此反映在两点间的重力段差值变化也就比较明显0图3高丽营!S 4.0地震前后重力段差值变化曲线Fig.3The curve Of the gravity difference vaiue befOre and after GaOiiyeng !S 4.0earthguake3.3重力场变化1977年就开始在北京地区进行流动重力复测0距高丽营4.0级地震震中最近的两个测段为顺义-高丽营和高丽营-小汤山0使用的仪器为加拿大CG-2型重力仪 李天生91989)0图3为两侧段重力段差值随时间变化曲线0由图3可以看出9自1978年4月开始到1979年1月9在这两个测段上采用了观测周期短达一个月左右的加密观测0 4.0级地震前9高丽营-小汤山的段差值从1978年4月开始就出现了30~50x10nm .S -2的异常变化9变化量均超过3倍中误差9并持续了半年之久0相邻的顺义-高丽营测段没有出现相应的明显变化04大同-阳高!L 6.1地震4.1地震概况!"李清林等:地震前后重力场变化的典型震例分析!期1989年10月19日1时9在山西省大同-阳高间发生了6.1级地震9在此之前的10月18日22时和之后的19日2时9还分别发生了5.7级和5.6级两次地震9为震群型的地震活动04.2地壳深部构造特征京西北地区为一地壳局部隆起区9在地壳上部和下部均存在着低速夹层9这些是地壳深部物质变异的深部构造条件 李清林等91991 0在大同-阳高之间发现地壳至地幔存在深大断裂带9其重要证据为pm 波及其一系列反射波在相应部位均发生了明显扭曲9莫霍界面落差约4km 0壳内存在低速层9其埋深在10~20km 之间0这些资料一方面说明9大同-阳高地震区下部存在着发震的地壳深部构造背景 另一方面9该区低速层和京西北地区 北京西山\延庆\宣化一带 的低速层埋深一致 另外大同-阳高地震的震源深度为10~15km 0这些资料可能说明9该地震的震源区与京西北地区的地壳深部有良好的连通9震源区地下应力场的变化以及地震活动必将对京西北地区的地壳构造运动\地球物理场的变化产生影响04.3重力场变化由图4营城子测点为代表的京西北地区的重力点值变化曲线可见9自1987年5月开始9该点出现了超过2倍中误差的异常量9并且异常幅度逐渐增大9至1988年4月异常值达70x10nm S -297月23日在其西北200km 处发生阳原4.8级地震0震后9该点点值出现了仅25x10nm S -2的下降 在重力值尚未恢复到正常范围时91989年该点重力值再次增大9且超过了阳原地震前的异常量级9达95x10nmS -20之后9发生了大同-阳高6.1级地震0图4大同-阳高!"6.1地震前后重力点值变化曲线Fig.4The curve of the gravity point vaiue before and after Datong-yanggao !L 6.1earthguake高丽营\阳坊\东园\怀来\张家口\柴沟堡等测点的重力趋势性异常变化形态基本一致9只是在起伏形态和量级上有所差别0这说明京西北地区处在同一个构造应力场作用下9只是因为各测点所处的具体构造部位不同9其重力变化量级就有相应的差别05太原!L 4.5和灵石!L 4.#地震5.1地震概况1993年6月26日在山西太原附近发生了!"4.$级地震91994年4月26日在山西灵石附近发生了!L 4.8级地震0这两次地震均发生在中国地震局1993年初确定的地震重点监视区和地震危险区内9并位于地球物理勘探中心的重力监测网内9为研究这两次地震前后的重力场变化提供了可靠的资料 李清林等91997 05.2地壳构造与地热场特征5.2.1地壳构造特征晋中南地区北部9主要由两条"#向的交城断裂和太谷断裂控制9在太原盆地内9发育一些次一级的"#和"$向断裂0在太原盆地的南部9有汾阳-洪洞断裂控制0晋中南地区的地壳深部构造特征则表现为9太原盆地东侧为地壳界面上隆区9西侧为地壳界面!"!!东北地震研究!!卷下坳区9太原\清徐\介休一带的地壳厚度约为40km O 这种地壳界面大起伏的变化特征9反映了山西断陷带具有大陆裂谷性质05.2.2地热特征山西断陷带各断陷盆地内部9具有较高的大地热流值和地温梯度值0高热流值的断陷盆地差异运动明显9深部活动性强9地震活动水平高;反之9低热流值的断陷盆地周围9尤其是隆起区9深部活动性不强9地热活动和地震活动水平较低05.3重力场变化以沁水测点为起算零点9推算各测点的重力值0重力点值平均中误差在 8X10nm .S -2以内0自建点至1994年8月9该测网共进行10次测量0图5太原!"4.#!灵石!"4.$地震前后部分重力点值变化曲线Fig.5The curve of Some gravity point vaiue before and after Taiyuan !L 4.3and LingShi !L 4.4earthguake 由图5可以看出9榆次\太原\晋祠\文水\杏花村\大麦郊几个测点的重力值从1992年初就开始缓慢下降0整个测区北部的重力值呈整体下降9趋势性十分明显0下降幅度较大的测点集中分布在测区北环的西北侧9即沿交城断裂附近9下降幅度达-35X10nm .S -20在重力场下降过程中9于1993年6月26日发生了太原4.3级地震0地震发生时9复测人员正在测区工作9震前就已完成了测网北环的测量9震后立即重新测量该测网的北环9结果表明重力场仍继续下降0到1993年10月测量时9重力点值基本恢复到建点时的水平0之后重力场再次下降9在下降变化的过程中于1994年4月26日又发生了灵石4.4级地震9到1994年8月重力场再次恢复到建点时的水平0图6为祁县-文水\榆次-太原测段段差值变化曲线0该变化趋势和点值的变化趋势基本相同9也是在重力段差值下降变化的过程中发生地震0上述两测段则距震中最近9且又是垂直于构造线走向9因此对地震的反映较明显0图6太原!"4.#!灵石!"4.$地震前后部分测段重力差值变化曲线Fig.6The curve of the gravity difference vaiue in Some meaSuring Sectionbefore and after Taiyuan !"4.#and LingShi !"4.$earthguake !"!"李清林等1地震前后重力场变化的典型震例分析4期从该区重力场的变化可以看出以下特点1(1地震孕育过程中a重力点值呈下降变化a距震中越近a其变化幅度越大O震后几天a仍维持下降趋势S尔后a重力场逐渐回升O重力场下降-回升的变化过程a反映了该区应力场积累-释放-恢复的完整过程O2)该区南部重力点值的变化和北部的重力点值变化不一致a北部呈现出明显的趋势性下降a 而南部则不然O这说明距震中较远的山西南部和距震中较近的北部有着不同的变化特征O3)异常过程和异常量级与一个4~5级地震的水平基本相当O6张北!L6.2级地震6.1地震概况1998年1月10日11时50分a在河北省张家口西北的张北-尚义间发生M L6.2地震(4l.l N, ll4.3E)O该地震波及范围较广a西至包头呼和浩特a东到秦皇岛天津a南到石家庄太原都有不同程度的震感O其中张家口地区受灾面积约860km2O张北县8个乡76个行政村房屋倒塌及受损约80%a2万余户无家可归O尚义县5个乡房屋倒塌及受损约1万余间O地震造成约50余人死亡a2000多人受伤a其中重伤250余人O自张北6.2级主震发生后至3月底a在震中区发生数百次余震a其中2.0~2.9级117次a 3.0~3.9级31次4.0~4.9级19次O余震集中在主震区东侧约20km的范围内O主震震源深度10km(李清林等a1999)O6.2地壳深浅构造特征6.2.1震中区构造特征张北6.2级地震位于张家口-宣化盆地的北台地上O震中区为晚第三纪玄武岩大面积出露区a在其附近还可以见三大片第四纪玄武岩出露和8处第四纪火山口遗迹O尚义-汉塔坝断裂在震中附近通过O该断裂为燕山山前断裂带的一部分a是宣化盆地延(庆)-怀(来)盆地与燕山山脉的边界断裂a为晚第三纪以来活动性较明显的活动性断裂O6.2.2晋北构造特征从大的区域构造背景上看a震中区位于阴山燕山太行山三山交界处a是阴山构造带燕山构造带和山西断陷带的结合部O从应力场的观点看a该区的应力场较易集中a且受阴山燕山构造带和山西断陷带构造运动的影响强烈O6.2.3地壳深部构造特征在张北震区下方a深10~15km处有一低速体a张北6.2级地震震源即位于低速体的上界面处O不难推测a低速度薄层的软弱可塑性为该区的水平运动提供了产生水平自由滑动面的有利条件S而其切穿莫霍面的深大断裂则为该区的垂直运动和物质上涌创造了条件O受到来自太平洋板块俯冲带的影响和山西断陷带燕山阴山构造带的横向推挤作用a导致该区应力集中a更为该区的水平运动创造了良好的条件O6.3重力场变化张北6.2级地震震中距晋北测网最近的测点(大同)直线距离约150km a距最远的测点(太原)约380km O1997年1月建点至1998年1月10日张北地震前a晋北测网进行的三期重力测量分别为1997年1月(16个点) 1997年9月(全部测点) 1997年11月(全部测点)O6.3.1重力点值变化图7给出了晋北地区部分测点点值变化曲线O测量中a1997年9月相对于1997年1月a该区重力场无明显变化a点值段差值的变化均在误差范围以内O而1997年11月相对于1997年9月a重力场出现了明显的大面积的上升变化a变化量超过40x10im~S-2的测点共20个a占晋北地区总测点数的!!东北地震研究!!卷62.5%9且均分布在晋北测区的中北部0变化最大的测点为朔州\杜庄\砂河镇\神池\龙宫\原平9其变化量分别为92\73\73\70\70\72x10nm .S -20这6个测点均位于恒山和六棱山山前断裂带及其附近0图7张北!"6.2地震前后部分重力点值变化曲线Fig.7The curve of the gravity point vaiue in Some meaSuring Section before and after Zhangbei !L 6.2earthguake这次地震前后的重力场变化形态和1976年唐山7.8级地震\1981年道孚6.9级地震前后重力场变化形态基本相似0这些地震前9该区重力场持续上升变化9震后重力场逐渐恢复到基值水平06.3.2重力段差值变化张北地区的重力段差值也有不同程度的变化9图8中给出了部分测段的重力段差值变化0段差值最大变化量达74x10nm .S -20其中一些变化较明显的测段跨越了六棱山山前断裂带(大同-杜庄段43x10nm .S -209恒山山前断裂带(中庄铺-大营35x10nm .S -209系舟山山前断裂带(阳曲-忻州74x10nm .S -20及其附近06.3.3重力场平面变化震前的1997年9月相对于1997年1月的平面和立体变化图像呈现出等值线稀疏9相对较平坦9重力场变化梯度较小的特征0而1997年11月相对于1997年9月\1998年2月相对于1997年9月9其平面和立体变化图像呈现出等值线密集9奇峰突起9重力场变化梯度大小分异的特征0变化梯度较大的有两个条带9它们是朔州-大同-杜庄条带和娄烦-忻州-原平条带07山西临猗!L 5.0级地震7.1地震活动概况1998年7月11日19时04分9在山西临猗-永济间发生了!L 5.0地震(34 57'"9ll0 38'#09这是!"李清林等 地震前后重力场变化的典型震例分析!期1991年大同5.8级地震之后山西断陷盆地南部所发生的唯一的一次达到!L 5.0的地震9是山西地震带经l989~l993年的活动高潮期9l994~l997年的平静期之后的一次中强地震活动事件0该次地震为主-余震型9震源深度22km 0余震衰减较快0震源机制解的结果表明9其主压应力方向与山西乃至华北的主压应力方向基本一致9震源错动为走滑类型(李清林等92000D07.2地壳深浅构造特征山西中南部有三个主要断陷盆地 太原盆地 临汾盆地和运城盆地9其走向都是"#或""#向9由一系列"#或""#向的深浅断层所圈限0临猗地震发生在运城盆地9其发震构造为双泉-临猗断裂9是控制运城盆地北侧沉降和变形的主要断裂带0从地壳深部构造看9临猗地震震中位于太行山地壳厚度变异带牛状凸起处9地壳厚度为36.0km 07.3重力场变化7.3.l 点值变化自l997年9月开始9该区l9个测点的重力点值出现了明显的 范围较大的变化9变化量在-40~-ll5>l0nm ~S -2之间9占全部测点 5l 个D 总数的37%9且多集中在山西南部地区0距震中最近的测点是临猗9仅20km 左右0点值变化详见图907.3.2段差值变化有6个测段的段差值 绛县-东镇 东镇-闻喜 东镇-侯马 介休-灵石 介休-平遥 介休-孝义D 出现了明显的变化9最大变化量为-92>l0nm ~S -20段差值变化详见图l007.3.3趋势面变化该区重力点值平面变化等值线图 趋势面D 呈现出西南 东北负 正的反向变化9且中间的临汾地区变化梯度呈现出明显增大的态势0图8张北!L 6.2地震前后部分测段重力段差值变化曲线Fig.8The curve of the gravity difference vaiue in Some meaSuring Section before and after Zhangbei !L 6.2earthguake!"!!!!东北地震研究!!卷图9临猗!5.0级地震前后重力点值变化曲线LFig.9The curve of the gravity point vaiue before and after Lingi M L5.0earthguake图10临猗!5.0级地震前后重力段差值变化曲线LFig.10The curve of the gravity difference vaiue before and after Lingi!L5.0earthguake7.3.4震后该区重力逐渐恢复!到1998年底!基本恢复到基值水平"7.3.5预报实况地球物理勘探中心根据上述重力场变化情况9在1998年度地震趋势研究报告中做出了明确的预报意见9震后对应较好9受到中国地震局的通报表扬08河南原阳!L4.5级地震8.1地震概况2002年11月7日9在河南省原阳县发生了一次!L4.5地震0该次地震的微观震中位于原阳县官厂乡东偏北34 57'"9113 55'# 9宏观震中位于原阳县官厂乡全屋村\新厂村\刘固村一带9震中烈度V度强9震源深度18km0全屋村\新厂村部分房屋窗户玻璃有震碎现象9刘固村小学教室有掉砖和裂缝现象李清林等92003 08.2地质构造概况从区域地质构造分析9在震中区的西部是太行山山前断裂带9北部是汤阴地堑和与之相关的北北东向的汤西\汤东断裂9南部为黄河及其近东西向的黄河断裂9东部是北西向的新乡-商丘断裂0震中位于开封-原阳凹陷的中心9震中区新生代沉积厚度约为2900m08.3重力场变化8.3.1点值变化新乡\南新庄和开封测点分别位于新乡-商丘断裂的两侧9它们表现出不同的活动方式9新乡\南新庄测点自1996年12月开始下降9而开封测点则在基值附近震荡9到1999年10月新乡测点点值的下降接近最低点9而开封测点刚开始下降9与此同时南新庄测点则继续保持下降的趋势9这可能反!"李清林等地震前后重力场变化的典型震例分析!期映了新乡-商丘断裂活动的信息9开封\东明测点的点值经过2年多的下降7于2002年l0月7重力发生反向上升的转折变化7在上升的过程中7于2002年ll月7日发生了原阳!4.5级地震0图ll是豫鲁L测区中比较典型的四个测点的点值变化曲线0图ll原阳!4.5级地震前后重力点值变化曲线LFig.ll The curve of the gravity point vaiue before and after Yuanyang!L4.5earthguake图l2原阳!4.5地震前后重力段差值变化曲线LFig.l2The curve of the gravity difference vaiue before and after Yuanyang!L4.5earthguake8.3.2段差值变化!新乡-南新庄\南新庄-郑州测段分别位于新乡-商丘断裂的两侧7新乡-南新庄的段差值多年来一直呈趋势上升7南新庄-郑州的段差值多年来一直呈趋势下降9于2002年l0月复测时7其趋势异常没有出现十分明显的转折变化9之后7发生了原阳地震0四间房-滑县\长垣-慈周寨两测段则自l999即出现了明显的下降变化7这种变化一直持续到2002年6月7之后7出现反向变化7期间发生了原阳地震0图l2是该区部分重力测段段差值变化曲线09重力变化机理的讨论9.l与形变有关地震是震源区介质发生的快速破裂错动或是原断层的失稳0在地震的孕育过程中7由于应力\应变的积累7孕震区内介质的物理\力学性质会发生不同程度的变化0因此7地壳形变会导致地面的空间变化7孕震区内的介质密度也会发生不同程度的变化9在大地震的形成过程中7软流层还可能存在着物质的迁移等0所有这些都将导致重力场的变化7它们可以通过地表的重力测量反映出来09.2与断层蠕动有关地震之前7由于地壳内部应力的长期作用7导致活动性断层产生蠕动变化7这种蠕动变化以地面形变和物质质量变化两种形态表现出来0我国大部分地区现今构造活动的水平分量大大超过了垂直分量7一种解释是深部(可能为孕震层)的应变活动引起了地表垂直形变的巨大变化(李清林!"!!东北地震研究22卷等1997)O蠕动变化同时也是地壳内部物质的流动它控制着地块的水平与垂直运动流动运动导致应力的累积式增长最终引起破裂和地震的发生李清林等1999)0从这8次地震震例看有变化的测点和测段多位于该区主要断层及其深断裂带附近O这些部位均为构造脆弱部位也是温差变化大的部位水平和垂直地温梯度大温差热应力高差异活动强烈是推动断层蠕动的主要力源之一0这种现象表明无论是深部的地壳运动或是近地表的水平和垂直运动都将在深大断裂带附近形成应力的再积累从而导致断裂带及其附近的重力场有较明显的变化09.3与质量迁移有关陈运泰等1980)从理论上计算了质量迁移的重力效应认为某些大地震可能与地壳和上地幔内的质量迁移有关所观测到的重力变化大部分是质量迁移引起的0上述重力场变化大部分应该是由该区地壳深部质量迁移和断层蠕动同时引起的重力场变化0参考文献1程式等.1976年盐源地震前后流动重力的异常变化J.地壳形变与地震VOI.6NO.3p229-2341986.2李天生.北京高丽营4.0级地震前后重力异常变化J.华北地震科学VOI.7NO.1p73-781989.3李清林等.京西北地区重力场变化与大同-阳高地震.华北地震科学VOI.9NO.1p22-271991.4李清林等.太原和灵石地震前后的重力场变化及其成因初探J.地壳形变与地震VOI.17NO.4p52-561997.5李清林等.张北6.2级地震及晋北地区的重力场变化J.地震学刊NO.4p7-151999.6李清林等.临猗M L5.0级地震前后的重力场变化及预报过程J.华北地震科学VOI.18NO.4p13-192000.7李清林等.河南原阳2002年11月7日M L4.5地震前的重力场变化.科技学术论文集地震研究与工程抗震M.原子能出版社p111~1152003.8陈运泰等.1975年海城地震与1976年唐山地震前后的重力场变化J.地震学报19802(1)p21~30.ANALYSIS ON THE GRAVITY FIELD CHANGE OF THE TYPICAL CASES BEFORE AND AFTER EARTHOUAKELI Oing-Iin LI Tian-sheng OIN Jian-zeng XIE Ru-yi(Research Center of Exploration Geophysics,China Eagthguake Administration,henan Zhengzhou450002,China)Abstract The practicaI experience on earthguake prediction by mobiIe gravity measurement data of Re-search Center of ExpIoration Geophysics China Earthguake Administration were summed up in detaiI and the seismogenic circument gravity fieId change cause of formation prediction resuIt and so on about8 earthguakes were studied.The resuIt showed that the different gravity observationaI point appeared medi-um-short term anomoIies obviousIy which is rising-faIIing-rising or faIIing-rising-faIIing before a stroger event.The evoIution process before and after a event was refIected cIearIy by the gravity change image.Key words gravity fieId gravity fieId change earthguake case。