川西地区台阵环境噪声瑞利波相速度层析成像

2006年中国地震局地质研究所地震动力学国家重点实验室在川西地区(26°N~32°N,100°E~105°E)布设了由297台宽频带数字地震仪组成的流动观测台阵.利用该密集台阵29°N以北156个台站2007年1~12月份的地震环境噪声记录和互相关技术,我们得到了所有台站对的面波经验格林函数和瑞利波相速度频散曲线,并进一步反演得到了观测台阵下方2~35 s周期的瑞利波相速度分布图像．本文结果表明,观测台阵覆盖的川滇地块、松潘-甘孜地块和四川盆地的地壳速度结构存在显著差异,具体表现为：(1)短周期（2~8 s）相速度分布与地表构造特征相吻合,作为川滇地块、松潘-甘孜地块和四川盆地之间的边界断裂,龙门山断裂带和鲜水河断裂带对上述三个地块上地壳的速度结构具有明显的控制作用,四川盆地前陆低速特征表明相应区域存在较厚的（约10 km）沉积盖层;(2)中周期（12~18 s）相速度分布表明,川滇地块和松潘-甘孜地块中上地壳速度结构存在明显的不均匀横向变化,并形成了尺度不同且高、低速相间的分块结构,而四川盆地中地壳整体上已经表现出相对高速;(3)长周期（25~35 s）相速度分布表明,松潘-甘孜地块,特别是川滇地块中下地壳表现为广泛的明显低速异常,意味着它们的中下地壳相对软弱,而四川盆地的中下地壳呈现整体性的相对高速,意味着四川盆地具有相对坚硬的中下地壳,并且以汶川地震的震中为界,龙门山断裂带的地壳结构显示了北段为高速异常,南段为低速异常的分段特征．     

基于地震探测的盐源盆地构造特征及动力学意义

基于覆盖盐源盆地的短周期天然地震台阵和布设的一条人工地震测线所获得的地震数据,从中提取地震能量属性,并通过地震层析成像获得该地区的浅部地震速度结构,继而对短周期地震台阵一个月的噪声数据进行互相关得到经验格林函数,再通过时频分析获得相速度频散曲线,反演获得不同深度的S波速度分布.研究结果显示,盐源盆地地震特征主要分为三层...     

汶川Ms 8.0地震:川西流动地震台阵观测数据的初步分析

2006年10月,中国地震局地质研究所地震动力学国家重点实验室,在中国四川西部(100°~105°E,26;N)布设了由297个宽频带地震台组成的密集流动地震观测台阵。截至2008年6月,川西台阵已记录远震事件(mb>5.0,30;)690个。5.12汶川MS8.0地震的发生为检验川西台阵观测系统的观测能力和工作状况提供了机会。文中对5.12汶川地震及其强余震序列观测记录进行了初步分析。结果表明:1)已发布的5.12汶川地震及其强余震参数有必要进一步修正,它们的震中位置可能存在8~24km误差,主震震源深度为19km的估计可能更接近实际情况;2)2008年5月16日理县(MS5.9)地震的波场分析表明:该地震面波不发育,它不可能是一个浅源地震,地震动速度波场垂直及水平分量峰值的异常增大与断层密切有关,在不考虑地形地貌、土质条件等因素的情况下,震中距200~250km范围内的地震动速度峰值异常增大可达正常衰减的4倍以上;3)利用接收函数方法对四川盆地和松潘-甘孜地块地壳结构的初步分析表明:四川成都北侧地壳厚度达到46km,四川盆地地壳从西向东逐渐增厚,下地壳显示了较为坚硬的特征,松潘-甘孜地     

华北地区地震环境噪声特征研究

利用华北流动地震台阵观测的垂直分向连续波形数据, 通过计算功率谱密度和相应的概率密度函数, 对华北地区地震环境噪声特征进行了分析研究. 结果表明, 东部平原和沉积盆地2 Hz以上的高频环境噪声水平与全球新高噪声模型（NHNM）相近, 周期3——18 s的平均噪声水平低于NHNM和新低噪声模型（NLNM）的平均值; 山区及西部高原的高频噪声水平明显低于NHNM, 周期1——18 s的噪声水平大多明显低于NHNM和NLNM的平均值; 不同区域18 s以上周期的噪声水平差异相对较小. 流动地震台阵部分台站的环境噪声存在明显的昼夜变化, 个别台站噪声水平明显高于周边台站, 表明这些台站受人类活动干扰较大. 不同台站的噪声水平分析表明, 将台站布设在摆坑内, 能在一定程度上降低高频和低频段的噪声水平. 台站环境噪声特征的研究结果可为流动地震台阵观测数据质量的定量评估, 观测期间的台站优化调整等提供重要依据.      

昆仑山口西M_S 8.1地震同震影响场的数值模拟

根据野外调查得到的地震破裂以及同震位错分布,利用三维有限元模型研究了昆仑山口西MS8.1地震的同震位移场和同震应力场的分布,讨论了震后地震活动与同震应力场变化的关系。模拟的MS8.1地震同震位移场显示,东昆仑断裂带南北两侧的同震位移衰减程度存在差异,断裂带南部的同震位移衰减大于北部同震位移的衰减;同震应力场研究结果则表明,MS8.1地震的同震最大剪切应力主要集中在东昆仑断裂带的周边地区,最大剪切应力等值线沿断裂带密集分布,形成了高梯度带。MS8.1地震的发生使青藏活动地块区的地壳应力发生了剧烈的变化,可能导致研究区的地壳应力场得到松弛。库伦破裂应力的研究表明,昆仑山地震对其后发生的5次强余震有一定的触发作用,昆仑山口西MS8.1地震可能也触发了2003年的德令哈6.6级地震     

秦岭及周边地区背景噪声Love波层析成像

和其它地震成像方法比较,背景噪声成像方法能够较有效的获得地壳以及上地幔的速度结构.本研究将PhaseWeighted-Stack方法应用到了互相关函数的叠加过程中,有效的增强了Love波信号,压制了获得的格林函数中的噪声.本文将这种技术应用到了北京大学与中国地质科学院合作布设在秦岭及其周边地区的(69台)宽频带流动地震台阵,利用台站的水平分量噪声记录提取的Love波频散曲线来获得台站间Love波相、群速度,进而得到了从6 s到30 s周期的Love波相速度二维分布图像.进一步进行了三维SH波速度反演,获得了研究区地壳和上地幔SH波3D速度结构,为研究该地区构造运动和上地幔动力学提供基础观测资料.一个重要的结果显示秦岭造山带以及鄂尔多斯块体外部西北方向在50 km深度呈现较强的SH波低速异常,这个连贯的低速异常带是否对应于青藏高原上地幔物质在鄂尔多斯块体和四川盆地块体之间的向东扩展流动而形成的温度异常或者部分熔融体还有待于进一步验证.     

一种基于密集线性台阵的虚拟地震测深新方法及其应用

本文提出了一种基于密集线性台阵的虚拟地震测深新方法,该方法在引入密集线性台阵数据的基础上,对之前广泛应用的单台虚拟地震测深方法进行了改进.针对复杂构造地区,新方法不再依赖传统的平面波入射假设,而是通过直接提取测线上的Ss直达波与Ss Pmp反射波震相走时,来更准确地约束反射点处Moho面深度.我们将新方法应用于跨龙门山构造带布设的密集流动地震台阵数据,对扬子克拉通与青藏高原构造边界带下方的地壳厚度进行了研究.结果表明扬子克拉通西部地壳相对较薄,地壳在龙门山构造带内部沿北西向迅速增厚.该区域地壳厚度变化特征可能与自新生代以来印度—欧亚板块碰撞所导致的青藏高原隆升过程有关.与单台虚拟地震测深方法相比,新方法对于复杂构造区域地壳厚度,如跨龙门山构造带,有着更可靠的约束.     

分离三分量远震接收函数的多道最大或然性反褶积方法

迄今为止 ,广泛应用的分离接收函数方法主要基于单个台站的观测数据 ,并依赖等效震源假定的有效性。当地壳内存在明显的间断面时 ,等效震源假定不适于宽频带地震波形数据 ,接收函数的估计将在不同程度上遭到破坏。解决这个问题的出路在于避免使用等效震源假定。依据多道最大或然性反褶积原理 ,给出了利用多道观测数据分离接收函数的方法。我们的方法避免了等效震源假定 ,这对于在复杂构造环境条件下改善接收函数的估计具有实际价值。合成地震图的数值检验结果验证了给出的方法。利用该方法 ,我们得到了跨越大别山造山带 ,从江西大冶附近的大箕铺 (30°2 0′N ,1 1 5°0 3′E)到河南兰考附近的崔林 (34°4 0′N ,1 1 4°4 9′E) ,总长度约 5 0 0km的流动地震台阵剖面各台站的三分量接收函数估计。与基于单台数据的接收函数径向分量比较表明 ,两种方法给出的结果可能存在不可忽略的差别     

利用sPn震相测定芦山MS7.0级地震余震的震源深度

利用南北地震带南段密集流动地震台阵的观测数据,采用波形互相关方法拾取Pn波走时,应用滑动时窗相关法识别sPn震相,通过sPn与Pn震相之间的走时差测定了芦山地震序列中28个M_L4.0级以上余震的震源深度.结果表明,震源深度集中在10~20 km范围内,垂直余震带的北西-南东向震源深度剖面揭示,余震分布表现出西深东浅的特点,倾角大约为39°.这些余震在空间上具有较好的线性分布特征,推测可能发生在与主震有关的破裂面上或邻近位置,由此推测主震的破裂面倾角大约为39°.根据余震的空间分布特征,认为芦山地震的发震断层并非双石-大川断裂,可能是其东侧的一条隐伏断层.     

中国东北地区噪声层析成像

中国东北地区是中国唯一的深震孕育区,获取该区的壳幔结构,对于研究板块俯冲、深震以及火山活动等有重要的科学意义.本文利用该区159个固定台站2011年1月至2012年6月和27个流动台站2011年1月至2011年6月间的垂向波形连续记录,计算了台站间的预估格林函数,并采用基于连续小波变换的时频分析方法,测量了双台路径上瑞雷波的群速度和相速度频散曲线.通过质量控制和筛选,最终得到了2204条路径上周期5~40 s的群/相速度频散曲线.检测板测试表明,反演结果可以达到2°×2°的分辨.利用Ditmar & Yanovskaya反演方法,我们得到了研究区（105°E-135°E,39°N-52°N）周期8~30 s的瑞雷波的群速度和相速度分布图.不同周期的群/相速度分布图,反映了不同深度S波速度的横向变化情况.研究结果显示：中国东北地区的地壳上地幔S波速度结构存在横向非均匀性.短周期（如8 s）的群/相速度分布与地表地质构造具有明显的相关性,具体来说,山区显示为高速,沉积盆地显示为低速;随着周期的增大（如15 s,22 s）,地形的控制作用相对减弱;较长周期（如30 s）的群/相速度分布与地壳厚度密切相关.     

利用背景噪声互相关研究汶川地震震源区地震波速度变化

利用2007年3月至2009年3月四川数字地震台网的宽频带连续波形资料,通过计算地震背景噪声互相关提取台站对间的经验格林函数,在0.1～0.5 Hz频带下测量每天经验格林函数与参考经验格林函数的走时偏移,进而得到各台站对在该时段内的相对地震波速度变化.结果表明,2008年5月12日汶川M_s8.0级地震造成了震源区地震波速度的急剧降低,最大降幅达0.4%;大致以安县为界,余震带西南部地区在汶川主震后波速降即达到最大值,而东北部地区的最大波速降一般出现在主震后的1~4个月,相对地震波速度变化的这种分段特性与地震序列的时空分布特征有较好的对应关系;在震源区外围的四川盆地也观测到了震后波速降低,而川西高原内部则没有出现显著的波速变化.进一步的分析和计算结果表明主震的静态应力变化和强地面运动引起的地表破坏都不能很好地解释震后波速的急剧降低,地震导致的断层区内部结构破坏和周边介质应力状态改变可能是波速变化的主要原因.     

TWO-DIMENSIONAL WHOLE CYCLE SIMULATION OF SPONTANE-OUS RUPTURE OF THE 2008 WENCHUAN EARTHQUAKE USING THE CONTINUOUS-DISCRETE ELEMENT METHOD

The May 12, 2008 M_S7.9 Wenchuan earthquake is ranked as one of the most devastating natural disasters ever occurred in modern Chinese history. The Longmenshan Fault(LMSF) zone is the seismogenic source structure, which consists of three sub-parallel faults, i.e., the Guanxian-Jiangyou Fault(GJF) in the frontal, the Yingxiu-Beichuan Fault(YBF) in the central fault and the Wenchuan-Maowen Fault(WMF) in the back of the LMSF. In this study, geological survey and seismic profiles are used to constrain the faults geometry and medium parameters. Three visco-elastic finite element models of the LMSF with different main faults are established. From the phase of interseismic stress accumulation to coseismic stress release and postseismic adjustment, the Wenchuan earthquake is simulated using Continuous-Discrete Element Method(CDEM). Modeling results show that before the 2008 Wenchuan earthquake, the GJF becomes unstable due to the interaction between its unique fault geometry and the tectonic stress loading. In the fault geometry model, the GJF is the most gently dipped fault among the three faults, which in return makes it having the smallest normal stress and the greatest shear stress. The continuous shear stress loading finally meets the fault failure criteria and the Wenchuan earthquake starts to initiate on the GJF at the depth of 15~20km. The earthquake rupture then propagated to the YBF. At the same time, due to the GJF and YBF rupture, the interseismic stress accumulation has been greatly reduced, causing the WMF failed to rupture. Although the stress accumulation in the WMF has been reduced significantly after the earthquake, yet it has not been released completely, which means that the WMF likely has with high seismic risk after the 2008 Wenchuan earthquake. We also find that the stress perturbation caused by gently dipping segment of the fault can promote the passive rupture in the steeply dipping segment, making the upper limit of dip angles larger than traditional assumption.     

Internal structure of Longmenshan fault zone at Hongkou outcrop,Sichuan, China,that caused the 2008 Wenchuan earthquake

This paper reports the internal structures of the Beichuan fault zone of Longmenshan fault system that caused the 2008 Wenchuan earthquake, at an outcrop in Hongkou, Sichuan province, China. Present work is a part of comprehensive project of Institute of Geology, China Earthquake Administration, trying to understand deformation processes in Longmenshan fault zones and eventually to reproduce Wenchuan earthquake by modeling based on measured mechanical and transport properties. Outcrop studies could be integrated with those performed on samples recovered from fault zone drilling, during the Wenchuan Earthquake Fault Scientific Drilling (WFSD) Project, to understand along-fault and depth variation of fault zone properties. The hanging wall side of the fault zone consists of weakly-foliated, clayey fault gouge of about 1 m in width and of several fault breccia zones of 30–40 m in total width. We could not find any pseudotachylite at this outcrop. Displacement during the Wenchuan earthquake is highly localized within the fault gouge layer along narrower slipping-zones of about 10 to 20 mm in width. This is an important constraint for analyzing thermal pressurization, an important dynamic weakening mechanism of faults. Overlapping patterns of striations on slickenside surface suggest that seismic slip at a given time occurred in even narrower zone of a few to several millimeters, so that localization of deformation must have occurred within a slipping zone during coseismic fault motion. Fault breccia zones are bounded by thin black gouge layers containing amorphous carbon. Fault gouge contains illite and chlorite minerals, but not smectite. Clayey fault gouge next to coseismic slipping zone also contains amorphous carbon and small amounts of graphite. The structural observations and mineralogical data obtained from outcrop exposures of the fault zone of the Wenchuan earthquake can be compared with those obtained from the WFSD-1 and WFSD-2 boreholes, which have been drilled very close to the Hongkou outcrop. The presence of carbon and graphite, observed next to the slipping-zone, may affect the mechanical properties of the fault and also provide useful information about coseismic chemical changes.     

DISCUSSING THE INFLUENCE OF MAJOR EARTHQUAKE SEQUENCE IN THE LONGMENSHAN AREA ON THE LUSHAN EARTHQUAKE FROM STRESS EVOLUTION

Calculating the coseismic static Coulomb stress change induced by an earthquake and interseismic stress change permits to explain the distribution of aftershocks, the earthquake sequence and other seismic observations. Four earthquakes greater than M7 have occurred in the Longmenshan area before the 2013 Lushan earthquake since 1900. This paper analyzes the influence of these four events on the Lushan earthquake, the stress evolution after the Lushan earthquake accompany with strong earthquake sequence on Longmenshan Fault, and the stress state of the gap between the Lushan and Wenchuan earthquakes. To address these issues would help future seismic risk assessment in the region. We construct a three dimensional finite element model based on the geological structure, the deep inversion results of density and velocity, and the GPS and the stress observation data. The simulation results show that the annual variation rate of Coulomb stress is higher on the Xianshuihe fault and southern segment of the Longmenshan fault, which is consistant with the regional seismicity. The coseismic Coulomb stresses induced by Kangding, Songpan, and Wenchuan earthquakes at the Lushan earthquake epicenter is greater than 0, implying that the three earthquakes may promote the occurrence of the Lushan earthquake, especially the Wenchuan earthquake. The viscous relaxation is remarkable which cannot be ignored in the analysis of stress evolution. From the stress evolution of this area, we can find that the gap between the Wenchuan and Lushan earthquakes is still at a relatively high stress level after the Lushan earthquake.     

汶川地震断裂带东北端浅部结构的人工地震探测

结合汶川地震断裂带动态监测,利用快速响应探测系统,开展了断层带浅部结构人工地震探测.针对地震断裂带动态监测条件下的复杂波场和低信噪比的情况,在f-k波场分离的基础上,分别利用了折射波共中心点成像、面波速度反演、反射波叠加成像方法,进行了浅层断层和构造成像处理,并对处理结果进行了综合解释,给出了断裂带浅部断层分布和速度特征.为汶川地震龙门山断裂带东北端动态监测提供了基础结构信息,所发展的断裂带快速响应探测技术对于地震应急动态监测具有重要意义.     

Stress changes on major faults caused by <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript>7.9 Wenchuan earthquake,May 12, 2008

On May 12, 2008, a magnitude 7.9 earthquake ruptured the Longmenshan fault system in Sichuan Province, China, collapsing buildings and killing tens of thousands people. As predicted, aftershocks may last for at least one year, and moreover, large aftershocks are likely to occur. Therefore, it is critical to outline the areas with potential aftershocks before reconstruction and re-settling people as to avoid future disasters. It is demonstrated that the redistribution of stress induced by an earthquake should trigger successive seismic activity. Based on static stress triggering theory, we calculated the coseismic stress changes on major faults induced by the Wenchuan earthquake, with elastic dislocation theory and the multilayered crustal model. We also discuss the stress distribution and its significance for future seismic activity under the impact of the Wenchuan earthquake. It is shown that coulomb failure stress (CFS) increases obviously on the Daofu-Kangding segment of the Xianshuihe Fault, the Maqu and Nanping segment of the Eastern Kunlun Fault, the Qingchuan Fault, southern segment of the Minjiang Fault, Pengxian-Guanxian Fault, Jiangyou-Guangyuan Fault, and Jiangyou-Guanxian Fault. The increased stress raises the probability of earthquake occurrence on these faults. Since these areas are highly populated, earthquake monitoring and early disaster alarm system are needed. CFS increases with a magnitude of 0.03–0.06 MPa on the Qingchuan Fault, which is close to the northern end of the rapture of Wenchuan earthquake. The occurrence of some strong aftershocks, including three events with magnitude higher than 5.0, indicates that the seismic activities have been triggered by the main shock. Aftershocks seem to migrate northwards. Since the CFS change on the Lueyang-Mianxian Fault located on the NEE of the Qingchuan Fault is rather small (±0.01 MPa), the migration of aftershocks might be terminated in the area near Hanzhong City. The CFS change on the western Qinling Fault is around 10 Pa, and the impact of static triggering can be neglected. The increment of CFS on the Pengxian-Guanxian Fault and Beichuan-Yingxiu Fault southwest to the main rupture is 0.005–0.015 MPa, which would facilitate earthquake triggering in these areas. Very few aftershocks in these areas indicate that the accumulated stress has not been released sufficiently. High seismic risk is predicated in these areas due to co-seismic CFS loading. The Wenchuan earthquake released the accumulated CFS on the Fubianhe Fault, the Huya Fault, the Ha’nan-Qingshanwan Fault, and the Diebu-Bailongjiang Fault. The decrement of CFS changes on the Longquanshan Fault east to Chengdu City is about 0.002 MPa. The seismic activity will be depressed by decrement of CFS on these faults. Supported by Knowledge Innovation Program of Chinese Academy of Sciences (Grant No. KZCX-SW-153), National Natural Science Foundation of China (Grant Nos. 40574011 and 40474028)     

利用主动震源检测汶川地震余震引起的浅层波速变化

为了监测汶川地震后断裂带附近的波速变化,2008年6月,我们利用快速响应探测系统,在龙门山断裂带北端进行了为期3周的连续观测实验.实验使用电动落锤作为主动震源,GPS连续同步授时的高精度数采作为接收系统.利用记录到的高时间精度观测数据,结合互相关时延检测技术,计算了P波和面波走时的相对变化,来表征地下波速的变化.我们观测到两次面波波速的异常变化,变化幅度最高达到2%,远远大于气压变化所能引起的波速变化.结合该时段内地震活动记录,我们认为面波的波速变化是由在附近发生的两次地震事件的同震效应引起的,该观测结果同理论计算结果相符.     

基于库仑应力变化分析巴颜喀拉地块东端的强震相互关系

本研究基于分层黏弹介质模型,考虑同震位错效应和震后黏滞松弛效应,分析巴颜喀拉地块东端1976年松潘地震序列、2008年汶川8.0级地震、2013年芦山7.0级地震和2017年九寨沟7.0级地震等多次大地震的可能存在的触发关系,计算大地震引起的周边各主要断裂的库仑应力变化.结果显示,1976年松潘地震序列各次地震间关系密切,存在明显的相继触发作用;综合考虑同震和震后效应,汶川8.0级地震对同属于龙门山断裂带的芦山7.0级地震有触发作用,且震后效应影响不可忽略;1976年地震序列,特别是1976年8月16日7.2级地震促进了2017年8月8日九寨沟7.0级地震的发生;汶川地震对九寨沟地震的影响研究中,采用不同的汶川地震同震位错模型,计算结果有差异.综合考虑多次大地震对周边断裂带的影响,龙门山断裂带南段、鲜水河断裂带中南段、平武—青川断裂北段、灌县—安县断裂北段、文县断裂的累积库仑应力增加显著,巴颜喀拉地块东端的东昆仑断裂带东段、迭部—白龙江断裂带西段以及金沙江断裂带库仑应力亦有所增加.综合考虑各重要断裂带已有的大地震危险性分析结果和库仑应力变化计算结果,龙门山断裂带南段、鲜水河断裂带中南段、东昆仑断裂带玛沁—玛曲段和金沙江断裂带的发震紧迫性有所增强,需引起关注.     

The seismicity and tectonic stress field characteristics of the Longmenshan fault zone before the Wenchuan MS8.0 earthquake

The seismicity of Longmenshan fault zone and its vicinities before the 12 May 2008 Wenchuan M_S8.0 earthquake is studied. Based on the digital seismic waveform data observed from regional seismic networks and mobile stations, the focal mechanism solutions are determined. Our analysis results show that the seismicities of Longmenshan fault zone before the 12 May 2008 Wenchuan earthquake were in stable state. No obvious phenomena of seismic activity intensifying appeared. According to focal mechanism solutions of some small earthquakes before the 12 May 2008 Wenchuan earthquake, the direction of principal compressive stress P-axis is WNW-ESE. The two hypocenter fault planes are NE-striking and NW-striking. The plane of NE direction is among N50°?70°E, the dip angles of fault planes are 60°?70° and it is very steep. The faultings of most earthquakes are dominantly characterized by dip-slip reverse and small part of faultings present strike-slip. The azimuths of principal compressive stress, the strikes of source fault planes and the dislocation types calculated from some small earthquakes before the 12 May 2008 Wenchuan earthquake are in accordance with that of the main shock. The average stress field of micro-rupture along the Longmenshan fault zone before the great earthquake is also consistent with that calculated from main shock. Zipingpu dam is located in the east side 20 km from the initial rupture area of the 12 May 2008 Wenchuan earthquake. The activity increment of small earthquakes in the Zipingpu dam is in the period of water discharging. The source parameter results of the small earthquakes which occurred near the initial rupture area of the 12 May 2008 Wenchuan earthquake indicate that the focal depths are 5 to 14 km and the source parameters are identical with that of earthquake.     

Active source monitoring at the Wenchuan fault zone: coseismic velocity change associated with aftershock event and its implication

With the improvement of seismic observation system, more and more observations indicate that earthquakes may cause seismic velocity change. However, the amplitude and spatial distribution of the velocity variation remains a controversial issue. Recent active source monitoring carried out adjacent to Wenchuan Fault Scientific Drilling (WFSD) revealed unambiguous coseismic velocity change associated with a local M_s5.5 earthquake. Here, we carry out forward modeling using two-dimensional spectral element method to further investigate the amplitude and spatial distribution of observed velocity change. The model is well constrained by results from seismic reflection and WFSD coring. Our model strongly suggests that the observed coseismic velocity change is localized within the fault zone with width of ~120 m rather than dynamic strong ground shaking. And a velocity decrease of ~2.0 % within the fault zone is required to fit the observed travel time delay distribution, which coincides with rock mechanical experiment and theoretical modeling.