地震折射波法在郑州市西区浅层勘探中的应用.pdfVIP

第 2 8卷第 1 期 2 0 0 6年 3月 地震地 质 S EI S MOLOGY AND GEOLOGY V0 l _2 8No1 Ma r 20 06 地震折射波法在郑州市西 区浅层勘探 中的应用 吴 怡 嘉世旭 段永红 徐朝繁 杨卓欣 1 ) 郑州市地 震局 工程勘察事务 所 ， 郑州4 5 0 0 0 7 2 ) 中国地震局地球物理勘探中心， 郑州4 5 0 0 0 2 摘要 在城市活断层勘探中， 对于面波干扰大、 地震反射波法难以开展工作的区域， 可尝试利用 地震折射波法进行探测 ， 并对折射波法探测的原始记录采用时间项、 差异时距曲线和有限差分成像等 方法进行综合计算 、 分析， 以探索折射波法在城市活断层勘探中的应用成效。文中针对郑州市须水断 层西段浅层地震折射波法勘探记录， 利用时间项、 差异时距 曲线和有限差分成像等计算方法 ， 获取剖 面速度结构与界面构造； 综合震相特征、 计算结果等资料确定主要地层的界面深度和构造特征， 3种 方法都取 得了相近 的结论 。后又通过在测线上 4个钻 孔资 料 的验 证 ， 认为 3种方 法的计算 结果 与钻 孔资料相吻合， 说明折射波法勘探在城市活断层探测 中的应用是可行的。 关键 词 城市活断层 地震折 射波法 郑州市西 区 中圈分类号：P 6 3 1 4 2 6 文献标识码： A 文章编号 ： 0 2 5 3 4 9 6 7 ( 2 0 0 6 ) O 1 0 0 8 40 9 O 引言 郑州市东侧为广袤的华北平原 ， 西侧为 山地隆起 ， 位于河淮盆地西北缘 。区内大片区域 由 第四系覆盖 ， 基岩仅 出露在西北部和西南部隆起 区。E W 向构造体 系是 区域 内形成 最早 的构 造 ， 主要 由一些总体走向 E W 的压性、 压扭性高角度断裂构造 ， 与之伴生的轴向近 E W 的褶皱构 造组成。在第四纪地貌上， 区内黄河呈 E W 向于郑州北侧通过 ， 主要断裂乔楼 一白沙断层经须 水 、 郑州市 区 E W 向延伸( 图 1 ) 。 为完成郑州市的活断层勘探 ， 进一步探测这些覆盖于第 四系下的隐伏断层的规模 、 产状、 形 成时间、 是否活动等， 在郑州市东、 西区布置了若干条浅层地震反射勘探测线进行探测工作。在 郑州市东区 ， 主要地震波组清晰， 获得了较丰富的地震波组信息， 取得了较好 的地震反射结果。 但在郑州市西 区， 虽然仪器和探测方法一致 ， 也试验 了不 同的观测 系统 ， 由于地表速度、 第 四系 基底埋深 、 潜水面埋深等构造环境的影响 ， 造成 了强烈 的面波 ， 始终没有接收到清晰 的、 可连续 追踪的反射界面， 反射震相的到时区域被面波干扰压制 。为此 ， 提出利用地震折射 波法进行试 验探测 ， 以探索在 目标层位埋深较浅 、 地下潜水 面较深 的黄土覆盖区利用折射波法进行 活断层 勘探的可行性与可靠性。 我们对郑州市须水断层西段进行浅层地震折射 波法勘探 ， 利用多种计算方法 ， 进一步 确定 其不同地质层位、 断层性质及其构造特征， 结果表明折射波法在浅层勘探中有其独特的优越性。 收稿日期) 2 0 0 41 1 2 5收稿 ， 2 0 0 51 0 0 8改回。 河南省地质局 ， 1 9 8 0 ， 中华人民共和国区域地质调查报告( 郑州幅) 。 维普资讯 1 期 吴怡等 ： 地震折射波法在郑州市西区浅层勘探中的应用 8 5 团 圆 地震测线及编号 钻孔编号及位置 圈 困 断层 新生界等深线 图 1 郑州市西区断层分布及地震折射测线位置 F ig 1 Dis t r ib u t io n o f f a u l t s a n d l o ca t io n o f s e is mic r e f r a ct io n p r o f il e s in W e s t e rn Z h e n g z h o u 1 观测系统和震相分析 初至波走时资料在折射剖面的资料解释中与地壳中康氐界面、 莫霍面等界面的反射波组一 样是主动源地震探测中的一种重要数据。初至波震相具有震相明确、 识别简单可靠 的特点 。可 以把初至波看成从上地壳 回折到地面的回折波 ， 这是以假定上地壳地震波速度随深度增加而增 加为前提条件的。事实上这 种假定与实际情况非常接近 ， 特别是当地表为沉积地层或风化层 时 ， 由于地层的压实程度随深度而增加 ， 风化作用随深度 的增加而减弱或者裂 隙随深度增大而 减小 ， 地震波速随深度的增加而相应地加大 。初至波时距 曲线对随深度变化 的速度梯度非常敏 感 ， 能得到准确的浅部速度结构 ， 利用初至波进行走时反演是一种较为准确的方法 ， 特别是对反 射波施工难度大、 面波干涉严重的地区更是一种好的替代方法 。 地震折射测线在郑州市西区共布置了 3条 ， 分别为 X XQ 1 、 X S Q1和 X G Q1 ， 长度均为 2 k m。 使用德国 D MT公 司生产的 S U N MI T数字地震仪和国产 3 8 H Z检波器进行接收 ， 采用了相遇和追 逐相结合 的观测系统 ， 炮 间距 1 0 01 3 0 m， 接收点间距 为 5 m， 6 0道或 8 O道接收 ， 最大接收长度 维普资讯 8 6 地震地质 2 8 卷 为 6 0 0 m。这里以穿越须水断层西段的 X X Q1 测线为主进行重点讨论。 X XQ 1 测线位于须水西， 白寨以北 ， s N向垂直穿越乔楼 一白沙断裂须水段 ， 测线长 2 k m， 共 进行 3 3次激发 ， 最大接收距约 6 0 0 m。通过对地震数据 的回放整理 ， 发现在炮检距约 7 01 0 0 m 以远出现 1 组视速度约 1 7 k m s ， 可稳定追踪至接收段末端的强震相波组， 为叙述方便， 这里以 P q命名该波组 。为更好地 分析震相 ， 提取到时 ， 这里采用 1 02 0 H z ( 低频端 ) 和 1 2 01 5 0 H z ( 高频端) 带通滤波 ， 以 1 7 k m s 为折合速度绘制记 录截面图( 图 2 ) 。 善 屋 留 罢 厦 翟 桩 号 图2 郑州市西区地震折射 X XQ 1 测线折合记录截面图 F ig 2 R e co r d i n g s e ct io n o f X X Q I p r o f i l e in w e s t e r n r e g i o n o f Z h e n g z h o u 桩号数据表示炮点在测线上 的距离 ( m) 图 2显示 ， 区域内第四系沉积厚约 4 01 0 0 m 的底部 反射震相应 出现在近炮点 折合到时 1 0 01 8 0 m s 的区间， 但该区间被面波和低频 波强烈干扰 ， 通过数字滤波 、 去噪及面波切 除等处 理方法也没能提取出有效反射震相。折射波速度大 ， 作为初至波可清晰、 可靠地追踪。 X X Q 1测线虽然仅长 2 k m， 但折射震相也显示出 3段不同的到时特征 ： 测线南段 ， 约 2 2 0 0 m 桩号以南各炮 ， 在炮检距约 7 O 4 0 0 m的区间段能对初至波组 P q进行平稳可靠 的追踪 ， 该波组 震相折合到时约 1 0 01 2 0 ms ， 视速度随炮检距 的加大略有增加 ， 约为 1 81 9 k m s ， 显示 了较 为平稳的界 面构造 ( 图 2 a ) ； 测线 中偏北约 2 4 0 02 9 0 0 m 区段 P q波到 时明显滞 后 ， 约 1 2 0 维普资讯 1期 吴 怡等： 地震折射波法在郑州市西区浅层勘探中的应用 8 7 1 4 0 m s ， 视 速度偏低 ， 走 时 曲线 呈现 出下 凹特征 ( 图 2 b ) ； 测线 北段 ， P q波组 到时又提前 至约 1 0 01 2 0 m s ， 近炮点走时滞后 ， 随炮检距增加走 时超前 ， 显示 出由南向北界面上倾 的特征。另 外 ， 近炮点初至波组随炮检距视速度 的增加而迅速增加 ( 0 51 0 k m s ) ， 显示为强速度梯度介 质结构特征 。 X X Q1 测线 的 P q波组到时主要显示了 2个特征 ： 一是为可长距离追踪 、 视速度约 1 7 k m s 的稳定震相波组 ； 二是该波组到时在 2 4 0 0 2 9 0 0 m区段相对于两侧到时明显 滞后 ， 显示 了与 之相关的界面构造以及须水断层的可能存在。X S Q 1 和 X G Q 1 测线 P q 震相特征与之类似， 追踪 相对平稳 。 通过资料分析 ， 区内 3条测线都显示了可长距离追踪 、 视速度约 1 7 k m s 的 P q强震相为本 次折射勘探的主要波组和 目的层 ， 该波组到时的超前、 滞后或局部突跳 ， 较好地体现了与之相关 的 P q界面的构造形态及可能存在的断层 。下面是根据 P q到时、 利用有限差分速度反演成像 、 时间项反演方法获取的剖面速度结构与界面构造 。 2 初至波有限差分层析成像 初至波有 限差分层析成像方法将初始速度模型网格化 ， 通过射线追踪得到研究区所有射线 路径经过的网格单元到达地表的理论走时， 与实际观测走时比较， 用拉普拉斯插值方法得到研 究 区域所有射线网格单元的慢度修正值。这样得到的速度模型作为下次迭代的初始模型 ， 多次 迭代后 ， 当走时残差达到一稳定值时， 停止迭代 ， 得到的模型即为反演结果模型。原理在有关文 献 中已有详细介绍( H o l e ， 1 9 9 2 ； 段永红等， 2 0 0 2 ) ， 这里主要介绍计算结果 。 图 3是 X X Q1 测线初至波有限差分成像结果 ， 包括 3方面内容 ： 图 3 a是有 限差分反演速度 结构 ； 图 3 b是射线轨迹密度分布， 可用来判定速度结构不同区域的可靠性 ； 图 3 c 是有限差分算 法理论时距曲线与实际时距曲线的吻合情况。 图 3 c 中， 粗线条是实际走时 ， 细线条是理论走时。X X Q 1测线的 3 3次炮点激发 ， 其中有 l 0 炮点为相对长距离( 约6 0 0 m ) 观测接收， 其余为2 0 0 4 0 0 m观测接收， 主要观测段覆盖次数达 6 次 。由图 3 c 可 以看 出， 理论时距曲线与实际时距曲线能很好地吻合。 从图 3 b可见， 测线两端即 3 4 0 0 m的桩号几乎没有射线分布 ， 深度 6 5 m 以上 的射线分布均匀且密集 ， 6 5 m以下 的不同区域射线密度差异较大。射线密度稀疏 区域 的结果 只作为参考或不作解释。可解释区间主要为 1 4 0 0 3 3 5 0 m桩号间和 6 5 m深度以上范 围， 在 图 3 a中用虚线表示。 图 3 a速度结构显示 ： 在约 6 0 m深度处出现了上下速度 明显变化 的界面特征 ， 该深度以上的 介质岩性是速度等值线密集的强梯度结构 ， 以下是速度 等值线迅速变缓的弱梯度结构 ， 显示了 良好 的界面构造 ， 该界面附近的速度约 1 7 k m s ， 与震相分析中讨论的区内可长距离追踪 、 视速 度约 1 7 k m s的 P q折射震相对应 ， 说明 P q震相 的确较好 地反映了与之相关的界 面的构造特 征 ， 为了叙述方便 ， 这里命名为 P q界面； P q界面总的趋势是从小桩号到大桩号( 由南向北) 缓慢 加深 ， 深约 5 5 6 5 m， 并在 1 8 0 01 9 0 0 m桩号 、 2 9 0 03 0 5 0 m桩号段 明显凸起 ， 凸起处相对高 差分别约为 8 m和 1 2 m， 埋深分别约为 4 7 m和 5 2 m； 特别要指出的是在 2 5 0 0 2 8 0 0 m桩号段 P q 界面迅速下陷， 构造形态南缓北陡 ， 最深处在 1 7 0 0 m桩号， 深约 7 5 m， 与两侧形成约 l 52 0 m的 断差 。前面震相分析也显示了该局部区域走时的明显滞后 ， 可能与 X X Q 1测线穿越垂直的须水 维普资讯 8 8 地震地质 2 8卷 0 2 O 目 一 4 0 点 聪一 B 0 一 l 0 0 1 4 0 0 1 6 0 0 1 8 0 0 2 0 0 0 2 2 0 0 2 4 0 0 2 6 0 0 2 8 0 0 3 0 0 0 3 2 0 0 3 4 0 0 测线长度 _ 0 1 0 2 0 3 0 - 4 0 5 0 6 0 7 0 - 8 0 - 9 0 0 6 0 8 L2 1 4 1 6 1 8 速度 k _ S - 1 2 2 测线长度 - 测线条数 根 测线长度 _ 图 3 X X Q !测线有限差分初 至波反 演结果 F i g 3 I n v e r s io n r e s u l t s o f X X Q 1 p r o fil e u s i n g fin i t e d iff e r e n ce m e t h o d a反演速度结构 ； b射线密度分布 ；c 走 时曲线拟合 断层相关。图 3 a显示更深处断层仍有延伸 ， 但射线稀疏 ， 只作参考使用。 P q界面以上速度等值线密度稳定 ， 震相分析也显示仅为炮检距 5 0 m以内的接收 ，未能形成 目、 馘翳 目， | 匿窨 维普资讯 1 期 吴怡 等 ： 地震折射波法在郑州市西 区浅层勘探 中的应用 8 9 有效的超浅层折射追踪 ， 难 以进一步细划分 ， 总体显示为强速度梯度结构 。中间层速度等值线 横 向变化较小， 下陷断层构造主要体现在 P q界面。 X S Q1 测线在 1 7 2 0 m桩号为相对低点 ， 深约 8 2 m， 由该低点 向南北两侧缓慢抬升 ， 界面形态 南陡北缓， 在 1 4 0 0 m， 2 9 0 0 m桩号深约 6 5 m， 该相对低点下凹构造是否为上街断层 ， 须进一步勘 探和钻井证实。 X G Q 1测线在约 8 5 01 O 0 0 m桩号 P q界面有一 明显凹陷， 速度等值线 与西侧形成 约 1 5 2 0 m的高差 ， 产状西陡东缓 ， X G Q 1测线西侧 8 5 0 1 O 0 0 m桩号段的明显凹陷暗示了古荥隐伏断 层存在的可能性 ， 须进一步勘探 、 钻井确定 。 3 时间项方法 为了使结果更加可靠 ， 对同一种资料应用不同的方法进行计算 分析 ， 下面是利用时间项法 反演界面构造 。1 9 5 7年 S ch e id e g g e r 等人提出了时间项法 ， 该方法可以利用首波走 时资料反演 炮点、 观测点的时间项和首波沿界面的滑行速度 。它不要求爆破点与观测点成规则分布 ， 并可 利用多次爆破的资料 ， 因而， 特别适合于大量爆破资料 的分析处理。该方法要求爆破点与观测 点至少有一个重合 ， 这时即能保证问题的解答是惟一 的( S ch e id e g g e r e t a 1 ， 1 9 5 7 ； 陈颐等 ， 1 9 7 4 ； 段永红等， 2 0 0 2 ) 。 时间项反演结果包括 2方面内容 ， 一是时间项反演走时拟合图， 显示了时间项法计算理论 走时与实测走时的吻合情况， 二是时间项法反演的深度界面构造。 图4 a ， b是 X X Q 1测线时间项法反演在 2 4 0 0 3 O 0 0 m桩号段凹陷构造的走时拟合图。图 中空圈表示任意给定的初始模型走时曲线 ， 实圈表示实测走时， 叉圈表示理论走时。由图可见 ， 时间项法反演理论走时很好地拟合了实测走时， 大部分几乎与实测走时重合 ， 个别点拟合误差 约 5 ms 。图4 a的 2 4 0 0 m桩号炮点两侧 P q走时南超北滞差异 ， 相同炮检距走时差达 2 0 m s ， 显示 了该炮点向北 P q 界面下倾造成走时滞后 ； 图4 b的 3 O 0 0 m桩号炮点 P q走时平稳 ， 整个接收段 折合走时几乎不变甚至略有滞后 ， 显示了该炮点向南 P q界面下倾 ， 走时滞后。 图 4 c 是 X X Q 1测线时间项法反演 P q震相波组的界 面构造埋深及产状 。与有 限差分结果 相比， 取得了较好的一致性 ， P q界面埋深约以5 0 6 0 m为轴 ， 北端略深， 突出显示 了 1 8 0 0 m桩 号、 3 O 0 0 m桩号附近的凸起及 2 7 0 0 m桩号 的凹陷 ， 2 2 0 0 m桩号的局部凸起在有限差分结果中 仅有轻微体现。特别要强调 的是 2 7 0 0 m桩号的凹陷构造 ， 该凹陷大约位于 2 3 5 0 3 O 0 0 m桩 号之间 ， 最深约 7 0 m， 与两侧形成约 2 0 m的断差 ， 产状两侧基本对称 ， 有限差分结果为南缓北陡， 但规模及埋深基本相同 ， 进一步证实了须水断层存在的可能性 。 X S Q1 测线结果显示 ， P q界面构造简单 ， 测线中段 由南 向北为单一上倾 ， 界 面埋深约 7 3 4 5 m。 在约 1 6 0 0 m， 2 9 5 0 m 桩号界 面出现 了相对低点 ( 约 7 3 m) 及 高点 ( 约 4 5 m) 埋深 ， 北端 2 9 5 0 m桩号的凸起产状平缓 ， 与两侧 2 6 0 0 m， 3 3 0 0 m桩号处相对高差约 1 5 m， 南端 1 6 0 0 m桩号 的凹陷南侧界面迅速抬升 ， 与 1 4 5 0 m桩号处形成近 2 0 m的断差。X S Q 1 测线 P q界面 2 3 0 0 m桩 号以西埋深 5 5 4 0 m， 总趋势 由西向东缓慢抬升 ， 中间有 局部小 的起伏。相对凹陷在约 9 0 0 1 O 0 0 m桩号处 ， 该凹陷东侧较为平缓 ， 西侧界面抬升 ， 在约 8 5 0 m桩号埋深 4 3 m， 形成约 1 2 m的 断差。测线东侧 2 3 0 0 m桩号 以东界 面迅 速抬升 ， 2 4 0 0 m桩号埋深 2 5 m， 形成 约 1 5 m 的断差。 X S Q 1 ， X G Q 1 测线时间项法结果与有限差分相 比， 构造特征相似 ， 埋深略浅。 维普资讯 地震地质 2 8 卷 昌 窿 翟 目 嫩 测线长度 k 0 初始走时 实测走时 +计算走时 图 4 郑州市西 区须水测线 ( X X Q 1 ) 时间项反演结果 F ig 4 I n v e r s i o n r e s u l t s o f X X Q 1 p r o fil e u s in g t ime t e r m s m e t h o d in w e s t e r n Z h e n g z h o u a桩号 2 4 0 0 m炮点测线时间项反演走时拟合 ；b桩号 3 0 0 0 m炮点测线时间项反演走时拟合 ； c 测线时 间项反 演结果 4 综合解释 两种方法都取得了相似的界面深度和构造特征 ， 由于时间项法是利用待测定界面以上至地 表层的平均速度计算该界面 的埋深及起伏变化 ， 所 以这种方法得到的结 果对界面埋深较为敏 感 ， 但过分依赖上层平均速度以及平均速度的横 向变化会影响解 的稳定性。有限差分成像主要 对速度敏感 ， 它利用所有炮点的初至波到时联合反演 ， 得到的是纵横 向非均匀构造条件下的二 维速度构造剖面 ， 该方法对界面不如时间项 敏感 ， 它通过速度等值线 的梯度变化来确定界面。 所 以综合应用两种方法 ， 利用有限差分反演以确定较为准确的速度结构 ， 计算出相对准确 的平 均速度 ， 使时间项法在一定的约束速度条件下进行界面计算 ， 用所得 到的界面反过来进一步确 定有限差分成像界面构造 ， 可使浅层地震折射勘探结果更加可靠 。 以有限差分速度结构为主 ， 其中界面埋深构造参考 了时间项法的结果 ， 并利用该区域钻井 勘探及地质资料 ， 绘制了 X X Q 1折射测线综合解释图( 图 5 ) 。 根据区域地质调查报告得到 ， 测线所在的小区域 内新生代沉积厚度 由西南深 2 0 0 m 向东北 逐渐加深为 1 4 0 0 m， 钻井勘探显示区域内第 四系沉积( Q) 厚 4 48 2 m， 其下为新近纪地层 ( N) 。 这次地震折射勘探 ， 折射波组 P q是区域内惟一可长距离稳定追踪 的震相波组 ， 结果显示区域内 介质速度约 1 7 k m s的构造界面在深约 4 0 8 0 m的区间变化 ， 上下层介质速度梯度明显不同 ， 维普资讯 1 期 吴怡等 ： 地震折射波法在郑州市西区浅层勘探中的应用 9 1 其构造埋深与第 四系沉积厚度相近 ， 很可能反映了区域第 四系沉 积厚度 ， P q波组是第 四系底 部 、 新近系顶部的折射首波。 桩号 三 三 推 测 断 层 回 速 度 k n s ， 圈 钻 进 勘 探 图 5 郑州市西区 X X Q 1折射测线综合解释 图 F i g 5 S y n t h e s i s e x p l a n a t i o n r e s u l t s o f X X Q 1 p r o f i l e o f w e s t e r n Z h e n g z h o u 桩号数据表示炮点在测线上的距离( n 1 ) 如 图5 ， X X Q 1 测线在 2 7 0 0 m桩号显示了较为明显 的断层迹象 ， P q界面在约 2 5 5 0 2 7 5 0 m 桩号处迅速下陷 ， 埋深约 8 0 m， 产状南缓北陡， 与两侧形成约 2 5 m一3 0 m的断差。为了确定该构 造是否存在 ， 在相应区域布置了 4个钻孑 L 进一步详细勘探 ， 这里仅 以位于 凹陷处的X X Z K一 3 钻 孔 ( 2 6 0 0 m桩号 ) 与 X X Z K一 4钻孔 ( 2 8 5 0 m桩号) 进行分析 ， 可见深约 1 5 m以上为相 同的粉土、 粉质黏土 ； 约 2 0 m 以下就显示 了明显的不同， X X Z K一 3在 2 0 m， 3 2 m， 4 8 m， 6 2 m， 7 3 m取样 ， 岩性 分别是粉质黏土 、 粉砂 、 粉砂 、 细砂岩、 黏土( 坚硬一硬 塑) ， 而X X Z K一 4 在相 同深度位置分别是 黏土、 细砂岩 、 砂砾岩 、 粉质粉土( 坚硬一硬塑) 及细砂 ( 密实 ) 。北侧 凸起 ( X X Z K一4 ) 地层介质 岩性与南侧凹陷区( X X Z K一3 ) 明显不同， 北侧凸起区主要以细砂岩、 砂砾岩等强度较大的岩性 构成， 南侧凹陷区则主要以粉砂 、 黏土等相对松散的岩性构成 ， 较好地体现了断层两侧介质岩性 的差异及构造产状特征 ， 进一步证实了 E W 走向的须水断层的存在。 X S Q 1测线南侧约 1 7 7 0 m桩号处 ， P q界面迅速加深至 7 6 m， 明显下 凹、 南陡北缓 ， 根据上伏 速度等值线得到 ： “ 疑似断层” 向 S 倾斜 ， 相对断差约 1 5 m。X G Q 1测线西侧约 9 0 0 m桩号， P q界 面明显下凹， 西陡东缓 ， 界面埋深约 5 6 m， 相对断差约 1 6 m， “ 疑似断层” 向 w 倾斜。这 2条测线 所显示的 P q界面畸变凹陷规模不大， 是否与上街 、 荥古断层有关需进一步工作 ， 但至少在测线 的这 2个地方出现 了 P q界面的局部凹陷。 通过对郑州西区浅层地震折射波法勘探资料的处理解释 ， 显示了地震折射波法在浅层勘探 中的有效性。折射波法野外施工要求 比反射法相对宽松 ， 炮点密度 、 覆盖次数都相对较低， 勘探 成本较低。另外该方法对介质速度的特殊敏感是其它方法不 可替代 的。利用反射法对构造界 面 的敏感和折射法对速度的敏感 ， 联合使用这 2种技术将会使城市活断层地震勘探更加可靠与 有效 。但对那些地震面波干扰大， 反射法难 以开展的区域 ， 地震折射波法不失为一种好 的选择。 维普资讯 9 2 地震地质 2 8 卷 参 考 文 献 陈颐， 陈运泰1 9 7 4 晋中南地区地壳结构的初步研究 J 地球物理学报，1 7 ( 3 ) ： 1 8 6 1 9 9 C HEN Yo n g，C HEN Y u nt a i 1 9 7 4 A p r e l imin a r y r e s e a r ch o n cr u s t a l v e l o c it y s t r u ct u r e in ce n t r a l s o u t h e r n r e g io n o f S h a n x i J C h i n e s e J G e o p h y s ，1 7 ( 3 ) ： 1 8 6 1 9 9 ( in C h in e s e ) 段永红， 张先康 ， 方盛明 2 0 0 2 华北地区上地壳结构的三维有限差分层析成像 J 地球物理学报， 4 5 ( 3 ) ： 3 6 2 - 3 6 9 D UAN Y o n g - h o n g ，Z HANG Xia n - k a n g，F ANG S h e n g - min g 2 0 0 2 T h r e e - d ime n s io n a l f in it e d if f e r e n ce t o mo g r a p h y o f v e l o cit y s t ruct u r e o f t h e u p p e r crus t i n N o A h C h i n a J C h i n e s e J G e o p h y s ， 4 5 ( 3 ) ： 3 6 2 3 6 9( i n C h i n e s e ) H o l e J A 1 9 9 2 N o n l i n e a r h i 【g h - r e s o l u t i o n t h r e e - d i me n s i o n a l s e is mi c t r a v e l t im e t o mo g r a p h y J J G e o p h y s R e s ， 1 3 5： 6 5 53- 65 6 2 S ch e id e g g e r A E， W il l mo r e P L 1 9 5 7T h e u s e o f a l e a s t s q u a r e me t h o d f o r t h e in t e r p r e t a t io n o f d a t a f r o m s e is mic s u r - v e y s J G e o p h y s i cs ， 2 2 ( 1 ) ： 9 2 2 T HE AP P LI CATI ON OF RE FRACTI ON M ETHoD TO S HALLOW EXPLoRATI ON OF W ES TERN ZHENGZHOU AREA W U Yi J I A S h i x u DUAN Yo n g h o n g XU Ch a o f a n YANG Z h u o x in ) 1 ) D e p a r t m e n t o f E n g i n e e r i n g S u r v e y ， E a r t h q u a k e A d min i s t r a t i o n of Z h e n g z h o u C it y ， Z h e n g z h o u 4 5 0 0 0 7 ，C h i n a 2 ) G e o p h y s i ca l E x p l o r a t i o n C e n t e r ， C h i n a E a r t h q u a k e A d mi n i s t r a t io n ， Z h e n g z h o u 4 5 0 0 0 2 ，C h in a Abs t r a ct S e is mic r e f r a ct io n me t ho d is t e n t a t i v e l y u s e d in t h e e x p l o r a t io n o f ur b a n a ct iv e f a u l t s wh e r e s t r o n g in t e rfe r e n ce o f s u rfa ce wa v e s e x is t s a n d i t is d iffi cu l t t o u s e s e is mic r e f l e ct io n s u r v e yOr ig in a l s e is mic r e f r a ct i o n d a t a a r e ca l cu l a t e d a nd in v e r s e d b y us i ng t ime t e r ms，C urve o f d if f e r e n t ia l t ime d i s t a n ce a n d fini t e dif