浅层地震多波联合勘探方法.pdf

《浅层地震多波联合勘探方法.pdf》由会员分享，可在线阅读，更多相关《浅层地震多波联合勘探方法.pdf（25页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201911063256.6 (22)申请日 2019.10.31 (71)申请人 山东电力工程咨询院有限公司 地址 250014 山东省济南市历下区闵子骞 路106号 (72)发明人 孙旭张学强杨庆义付康伟 蔡伟刘浩孙长帅刘灿 李仁杰王冰焦学文任传健 于天文郭翔卢向星 (74)专利代理机构 济南圣达知识产权代理有限 公司 37221 代理人 李圣梅 (51)Int.Cl. G01V 1/30(2006.01) (54)发明名称 浅层地震多波联合勘探方法 (57)摘要 本公开提出。

2、了浅层地震多波联合勘探方法， 建立多波联合勘探观测系统进行数据采集； 基于 所采集的数据， 针对异常体探测， 联合使用联合 地震映像法、 瑞雷波相邻道相位差法提取相速度 剖面两种处理方法进行异常体探测工作， 并对结 果进行比对约束， 判断异常体的空间分布， 排除 假异常； 针对S波速度结构分层， 联合折射波层析 成像、 面波频散曲线反演两种处理方法进行联合 反演， 获取精度更高的S波速度结构。 相比现有的 单方法勘探观测系统， 其能够充分的利用地下的 地震波信息， 达到一次作业， 多波多方法并用的 目的。 权利要求书2页 说明书11页 附图11页 CN 110687602 A 2020.01.。

3、14 CN 110687602 A 1.浅层地震多波联合勘探方法， 其特征是， 包括： 建立多波联合勘探观测系统进行数据采集； 基于所采集的数据， 针对异常体探测， 联合使用联合地震映像法、 瑞雷波相邻道相位差 法提取相速度剖面两种处理方法进行异常体探测工作， 并对结果进行比对约束， 判断异常 体的空间分布， 排除假异常； 针对S波速度结构分层， 联合折射波层析成像、 面波频散曲线反演两种处理方法进行联 合反演， 获取精度更高的S波速度结构。 2.如权利要求1所述的浅层地震多波联合勘探方法， 其特征是， 建立多波联合勘探观测 系统进行数据采集工作， 其建立过程为： 炮点从最左边出发， 逐道间距。

4、移动， 即炮间距等于道间距， 排列不动， 直到炮点达到最 右端； 最左端、 最右端距最近的检波器距离需根据地震映像、 面波、 折射波的野外勘探偏移距 的选择要求来选择； 敲击完反向炮之后， 一个排列结束， 而后移动排列， 排列的移动需要满足排列与排列之 间至少有一道重合。 3.如权利要求1所述的浅层地震多波联合勘探方法， 其特征是， 异常体探测及S波速度 结构分层作为探测任务单独进行， 或同时进行。 4.如权利要求1所述的浅层地震多波联合勘探方法， 其特征是， 针对异常体探测， 具体 为： 基于多波联合勘探观测系统进行数据采集所获得的原始的地震记录， 提取不同偏移距 的地震映像剖面， 而后识别。

5、地震映像剖面上的异常区域； 基于原始的地震记录， 利用相邻道相位差法提取多道频散曲线， 再利用半波长经验公 式得到瑞雷面波相速度剖面， 基于相对速度剖面圈定异常； 基于地震映像剖面、 瑞雷面波相速度剖面结果， 进行比照确认， 从而得到更精确的异常 位置， 排除掉假异常。 5.如权利要求4所述的浅层地震多波联合勘探方法， 其特征是， 计算瑞雷波相速度剖 面， 需要运用多排列正反炮覆盖叠加的地震记录数据进行， 提取数据过程： 在利用多波联合 勘探观测系统采集到的数据中， 确定面波排列的长度与偏移距， 按该偏移距抽出此排列的 正反炮数据， 然后排列移动半个长列长度， 接着抽出正反炮数据， 直至到达测。

6、线末端， 这些 数据组成了后续处理的数据集。 6.如权利要求1所述的浅层地震多波联合勘探方法， 其特征是， 针对S波速度结构分层， 联合折射波层析成像、 面波频散曲线反演两种处理方法进行联合反演， 具体为： 基于多波联合勘探观测系统进行数据采集所获得的原始的地震记录， 抽出相遇或追逐 相遇观测系统， 然后拾取初至， 解释并建立初始模型， 再进行走时层析成像来反演地下纵波 速度， 并根据速度的分布进行分层； 基于原始的地震记录， 确定该测线上面波要做多少个排列， 确定偏移距， 然后抽出确定 好的排列的面波记录； 基于层析成像的结果， 计算出每个排列中心的每层平均纵波速度与层厚， 将其为先验 信息。

7、提供给面波进行频散曲线的反演， 从而得到更精确的S波速度结构。 权利要求书 1/2 页 2 CN 110687602 A 2 7.如权利要求6所述的浅层地震多波联合勘探方法， 其特征是， 抽出数据后， 利用能量 比法自动拾取初至， 拾取不准的需要进行调整， 而后利用t0差数法进行解释， 并利用解释的 结果建立初始模型， 初始模型是利用解释的第一层和基岩面的纵波速度， 从第一层到基岩 面均匀递增建立的。 8.如权利要求6所述的浅层地震多波联合勘探方法， 其特征是， 确定偏移距用以下标 准： 在非常软的地层： Vs100m/s， 偏移距应为15m； 在较软的地层： 100m/sVs300m/s， 。

8、偏 移距应为510m； 在较硬的地层： 300m/sVs500m/s， 偏移距应为1020m； 在非常硬的地 层： 500m/sVs， 偏移距应为2040m。 9.如权利要求6所述的浅层地震多波联合勘探方法， 其特征是， 瑞雷面波波速与S波波 速关系式为： 式中： 为泊松比， 该参数本身代表着地下介质的 “软” 、“硬” 程度， 它本身没有量纲值， 其值范围一般在00.5。 10.如权利要求6所述的浅层地震多波联合勘探方法， 其特征是， 用瑞雷波来反演S波速 度结构， 瑞雷波频散曲线的反演方法为带阻尼的最小二乘反演， 依赖初始模型的设定， 作为 先验信息参与并约束反演的层析成像信息， 参与的方。

9、式是在面波排列中心位置， 按一定半 径搜索每层在该处的平均纵波速度， 然后利用纵波与横波速度转换关系得到该处横波速度 作为初始模型的横波速度， 在中心位置的分层结果作为初始模型的层厚， 从而进行反演。 权利要求书 2/2 页 3 CN 110687602 A 3 浅层地震多波联合勘探方法 技术领域 0001 本公开涉及勘探技术领域， 特别是涉及浅层地震多波联合勘探方法。 背景技术 0002 浅层地震勘探作为一种高精度、 无损高效的地球物理勘探技术越来越多的进入到 工程物探， 城市勘探中。 目前， 常用的浅层地震勘探方法有反射波、 折射波、 面波、 地震映像 等几种方法， 其中反射波法主要用于地。

10、层的划分， 折射波法用于揭示地下速度结构特征， 面 波法用于提交测区S波速度结构资料， 地震映像则常用于异常体、 断裂的圈定与识别。 几种 方法都有其优势， 勘探中也往往根据探测的任务选用合适的方法， 但是单一的方法也存在 不足与局限性。 0003 发明人在研究中发现， 单一的勘探方法往往只能提供一部分勘探所需的资料， 并 且由于浅层地质情况的复杂性， 勘探的结果往往具有多解性， 单方法难以对这一多解性进 行约束以得到精度更高的结果。 具体来说， 地震映像法对异常体的横向位置有较好的探测 效果， 但纵向的分辨率不佳， 且难以进行地层分层工作。 折射波法无法提交S波速度资料。 瑞 雷面波法中频散。

11、曲线的反演结果往往受初始模型的影响较大， 往往需要施加一定的先验信 息才能得到较为准确的结果。 另外， 在实际的勘探工作中， 比如对地下破碎带这样的低速异 常体进行探测上， 由于第四纪覆盖层的地质情况复杂， 噪声的干扰较多， 单一的勘探方法得 到的结果往往是假的低速异常体， 提交的勘探结果也无法得到有效的比对检验。 0004 因此， 研究浅层地震多波联合勘探技术对高精度野外勘探来说具有十分重要的意 义。 发明内容 0005 本说明书实施方式的目的是提供浅层地震多波联合勘探方法， 得到精度更高的勘 探成果。 0006 本说明书实施方式提供浅层地震多波联合勘探方法， 通过以下技术方案实现： 000。

12、7 包括： 0008 建立多波联合勘探观测系统进行数据采集； 0009 基于所采集的数据， 针对异常体探测， 联合使用联合地震映像法、 瑞雷波相邻道相 位差法提取相速度剖面两种处理方法进行异常体探测工作， 并对结果进行比对约束， 判断 异常体的空间分布， 排除假异常； 0010 针对S波速度结构分层， 联合折射波层析成像、 面波频散曲线反演两种处理方法进 行联合反演， 获取精度更高的S波速度结构。 0011 进一步的技术方案， 建立多波联合勘探观测系统进行数据采集工作， 其建立过程 为： 0012 炮点从最左边出发， 逐道间距移动， 即炮间距等于道间距， 排列不动， 直到炮点达 到最右端； 说。

13、明书 1/11 页 4 CN 110687602 A 4 0013 最左端、 最右端距最近的检波器距离需根据地震映像、 面波、 折射波的野外勘探偏 移距的选择要求来选择； 0014 敲击完反向炮之后， 一个排列结束， 而后移动排列， 排列的移动需要满足排列与排 列之间至少有一道重合。 0015 进一步的技术方案， 异常体探测及S波速度结构分层作为探测任务单独进行， 或同 时进行。 0016 进一步的技术方案， 针对异常体探测， 具体为： 0017 基于多波联合勘探观测系统进行数据采集所获得的原始的地震记录， 提取不同偏 移距的地震映像剖面， 而后识别地震映像剖面上的异常区域； 0018 基于原。

14、始的地震记录， 利用相邻道相位差法提取多道频散曲线， 再利用半波长经 验公式得到瑞雷面波相速度剖面， 基于相对速度剖面圈定异常； 0019 基于地震映像剖面、 瑞雷面波相速度剖面结果， 进行比照确认， 从而得到更精确的 异常位置， 排除掉假异常。 0020 进一步的技术方案， 计算瑞雷波相速度剖面， 需要运用多排列正反炮覆盖叠加的 地震记录数据进行， 提取数据过程： 在利用多波联合勘探观测系统采集到的数据中， 确定面 波排列的长度与偏移距， 按该偏移距抽出此排列的正反炮数据， 然后排列移动半个长列长 度， 接着抽出正反炮数据， 直至到达测线末端， 这些数据组成了后续处理的数据集。 0021 进。

15、一步的技术方案， 针对S波速度结构分层， 联合折射波层析成像、 面波频散曲线 反演两种处理方法进行联合反演， 具体为： 0022 基于多波联合勘探观测系统进行数据采集所获得的原始的地震记录， 抽出相遇或 追逐相遇观测系统， 然后拾取初至， 解释并建立初始模型， 再进行走时层析成像来反演地下 纵波速度， 并根据速度的分布进行分层； 0023 基于原始的地震记录， 确定该测线上面波要做多少个排列， 确定偏移距， 然后抽出 确定好的排列的面波记录； 0024 基于层析成像的结果， 计算出每个排列中心的每层平均纵波速度与层厚， 将其为 先验信息提供给面波进行频散曲线的反演， 从而得到更精确的S波速度结。

16、构。 0025 进一步的技术方案， 抽出数据后， 利用能量比法自动拾取初至， 拾取不准的需要进 行调整， 而后利用t0差数法进行解释， 并利用解释的结果建立初始模型， 初始模型是利用解 释的第一层和基岩面的纵波速度， 从第一层到基岩面均匀递增建立的。 0026 进一步的技术方案， 确定偏移距用以下标准： 在非常软的地层： Vs100m/s， 偏移距 应为15m； 在较软的地层： 100m/sVs300m/s， 偏移距应为510m； 在较硬的地层： 300m/s Vs500m/s， 偏移距应为1020m； 在非常硬的地层： 500m/sVs， 偏移距应为2040m。 0027 进一步的技术方案，。

17、 瑞雷面波波速与S波波速关系式为： 0028 0029 式中： 为泊松比， 该参数本身代表着地下介质的 “软” 、“硬” 程度， 它本身没有量纲 值， 其值范围一般在00.5。 0030 进一步的技术方案， 用瑞雷波来反演S波速度结构， 瑞雷波频散曲线的反演方法为 带阻尼的最小二乘反演， 依赖初始模型的设定， 作为先验信息参与并约束反演的层析成像 说明书 2/11 页 5 CN 110687602 A 5 信息， 参与的方式是在面波排列中心位置， 按一定半径搜索每层在该处的平均纵波速度， 然 后利用纵波与横波速度转换关系得到该处横波速度作为初始模型的横波速度， 在中心位置 的分层结果作为初始模。

18、型的层厚， 从而进行反演。 0031 进一步的技术方案， 结合地质资料对处理的结果进行进一步的解释： 在异常体的 探测中， 探测到的高速低速异常体有多种解释结果， 这就需要结合前期地质资料进行约束； 0032 在S波速度分层中， 得到的S波速度需要进一步进行后续的处理， 从而用于后续的 解释指导工作。 0033 与现有技术相比， 本公开的有益效果是： 0034 本公开联合使用多种勘探方法， 把地下多种波作为有效波进行多波多方法并用的 联合勘探， 则会有效的缓解上述难题， 一是可以利用单方法的结果进行对比检验从而更精 确的确定地下异常体， 二是有条件进行多种波的联合反演， 从而得到精度、 分辨率。

19、更高的结 果。 另外， 进行一次作业， 多波多方法并用的勘探方式也有助于更充分的利用地下信息， 提 高勘探效率。 0035 本公开提出的一种适用于多波联合勘探方法， 相比现有的单方法勘探观测系统， 其能够充分的利用地下的地震波信息， 达到一次作业， 多波多方法并用的目的。 0036 本公开提出了多波联合勘探技术， 针对异常体探测， 联合使用联合地震映像法、 瑞 雷波相邻道相位差法提取相速度剖面两种处理方法进行异常体探测工作， 通过对结果进行 比对约束， 能够提高对异常体的定位精度， 消除 “虚假” 异常体； 针对S波速度分层来说， 本发 明联合折射波层析成像、 面波频散曲线反演两种处理方法进行。

20、联合反演， 具体是将层析成 像的结果作为先验信息提供给面波频散曲线反演， 有效的提高了反演的精度。 附图说明 0037 构成本公开的一部分的说明书附图用来提供对本公开的进一步理解， 本公开的示 意性实施例及其说明用于解释本公开， 并不构成对本公开的不当限定。 0038 图1是本公开实施例子的多波联合勘探观测系统图； 0039 图2是本公开实施例子的探测异常体与S波速度结构分层的多波联合勘探工作流 程图； 0040 图3(a)-图3(c)是本公开实施例子的验证低速异常体的探测效果， 进行三层层状 介质含低速体模型试算中， 模型示例图以及加入异常体位置后的三层层状介质含低速体相 速度剖面与映像结果。

21、图； 0041 图4(a)-图4(c)是本公开实施例子的验证高速异常体的探测效果， 进行三层层状 介质含高速体模型试算中， 模型示例图以及加入异常体位置后的三层层状介质含高速体相 速度剖面与映像结果图； 0042 图5(a)-图5(b)是本公开实施例子的验证S波速度分层的探测效果， 进行倾斜地层 模型试算中， 模型与层析成像图； 0043 图6(a)-图6(d)是本公开实施例子的验证S波速度分层的探测效果， 进行倾斜地层 模型试算中， 倾斜地层模型多个排列反演结果图； 0044 图7(a)-图7(b)是本公开实施例子的验证S波速度分层的探测效果， 进行三层层状 介质带倾斜地层模型试算中， 模型。

22、与层析成像图； 说明书 3/11 页 6 CN 110687602 A 6 0045 图8(a)-图8(d)是本公开实施例子的验证S波速度分层的探测效果， 进行三层层状 介质带倾斜地层模型试算中， 三层层状介质带倾斜地层模型多个排列反演结果图。 具体实施方式 0046 应该指出， 以下详细说明都是例示性的， 旨在对本公开提供进一步的说明。 除非另 有指明， 本文使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常 理解的相同含义。 0047 需要注意的是， 这里所使用的术语仅是为了描述具体实施方式， 而非意图限制根 据本公开的示例性实施方式。 如在这里所使用的， 除非上下文另外明。

23、确指出， 否则单数形式 也意图包括复数形式， 此外， 还应当理解的是， 当在本说明书中使用术语 “包含” 和/或 “包 括” 时， 其指明存在特征、 步骤、 操作、 器件、 组件和/或它们的组合。 0048 名词解释： 道间距： 相邻地震道中心点之间的距离； 炮间距： 相邻激发点之间的距 离； 偏移距： 炮点到最近接收道中心点的距离。 0049 实施例子一 0050 该实施例公开了浅层地震多波联合勘探方法， 浅层地震多波联合勘探技术， 首先 需要根据勘探区域地质情况， 勘探任务， 确定合适的偏移距， 建立多波联合勘探观测系统进 行数据采集工作。 本发明中， 多波联合观测系统示意图如图1所示。 。

24、0051 在具体实施例子中， 多波联合观测系统在组成上包括检波器、 地震记录仪、 放炮装 置、 线缆以及配套的数据采集处理软件等。 采集的是检波器接收的地震波数据。 传输至地震 记录仪上， 再传输到笔记本等电子计算机上， 在计算机上进行数据处理。 0052 检波器沿测线排列， 震源逐道间距激发。 需要注意的是， 炮点敲击的偏移距要满足 折射波相遇观测系统对偏移距的要求。 如图， 炮点左边出发， 逐道间距移动， 即炮间距等于 道间距， 排列不动， 直到炮点达到最右端。 当炮点敲击完最后一道检波器后， 可以选择不按 道间距逐点移动， 而是根据面波、 折射波的野外勘探要求选择合适的偏移距敲击反向炮。。

25、 敲 击完反向炮之后， 一个排列结束， 而后移动排列， 排列的移动需要满足排列与排列之间至少 有一道重合， 这样有利于后续的折射波层析成像。 0053 多波联合勘探技术， 在利用多波联合勘探观测系统采集到数据后， 需要根据勘探 任务采取相应的联合勘探技术进行进一步的处理。 具体的主要包含以下步骤： 0054 步骤一， 根据勘探区域地质情况， 勘探任务， 确定合适的偏移距， 建立多波联合勘 探观测系统进行数据采集工作。 0055 步骤二： 根据勘探任务， 确定需要使用的多波联合勘探技术方法， 然后利用多波联 合勘探技术处理步骤一获得的地震资料； 0056 地震资料具体的数据格式是地震数据采集中常。

26、用的数据格式， 比如segy格式。 数 据是以二进制的形式存储， 具体数据内容包括地震道数、 采样点数、 采样频率等元数据以及 检波器接收到的每一地震道数据。 0057 步骤三： 针对异常体探测， 联合地震映像时间剖面， 瑞雷面波相速度剖面处理技 术， 判断异常体的空间分布， 排除假异常； 0058 步骤四： 针对S波速度结构分层， 联合折射波层析成像， 瑞雷面波频散曲线反演进 行联合反演， 获取精度更高的S波速度结构； 说明书 4/11 页 7 CN 110687602 A 7 0059 在具体实施例子中， 异常体探测及S波速度结构两者可以同时进行， 也可以根据地 质条件与勘探任务择一进行。。

27、 0060 步骤五： 结合勘探现场的前期地质先验资料、 勘探成果， 对步骤三、 四获得的结果 进行进一步解释。 针对异常体探测， 结合测区资料， 推断出异常体可能所属地质体， 从而指 导后续的勘探； 针对S波速度分层， 获取更精确的S波速度信息后， 可以进一步更精确的推断 出地基承载力等一些岩土参数， 指导后续的工作。 0061 所述步骤一中， 通过勘探区域地质情况， 勘探任务， 确定合适的偏移距， 建立多波 联合勘探观测系统进行数据采集工作， 其建立过程为： 0062 多波联合观测系统示意图如图1所示。 0063 检波器沿测线排列， 震源逐道间距激发。 需要注意的是， 炮点敲击的偏移距要满足。

28、 折射波相遇观测系统对偏移距的要求。 0064 炮点从最左边出发， 逐道间距移动， 即炮间距等于道间距， 排列不动， 直到炮点达 到最右端。 当炮点敲击完最后一道检波器后， 可以选择不按道间距逐点移动， 而是根据面 波、 折射波的野外勘探要求选择合适的偏移距敲击反向炮。 最左端、 最右端距最近的检波器 距离需根据地震映像、 面波、 折射波的野外勘探偏移距的选择要求来选择。 敲击完反向炮之 后， 一个排列结束， 而后移动排列， 排列的移动需要满足排列与排列之间至少有一道重合， 一般移动半个排列的长度。 0065 所述步骤二中， 针对不同的勘探任务， 利用相应的多波联合勘探技术处理获取的 地震资料。

29、。 多波联合勘探技术主要针对的是异常体的探测与S波速度结构分层这两个勘探 任务， 可以单独进行， 也可以同时完成这两个任务。 0066 所述步骤三中， 针对异常体探测， 联合地震映像时间剖面， 瑞雷面波相速度剖面处 理技术， 判断异常体的空间分布， 其过程为： 0067 (3-1)预处理观测系统采集的原始的地震记录， 按需要重建观测系统， 方便后续处 理。 0068 在该实施例中， 因为要统一坐标系， 方便后续的处理与绘图， 需要重建观测系统。 后续处理是针对重建后的数据进行处理， 按照异常体探测或者S波速度分层的工作流程进 行后续数据处理， 之前的数据不参与后续的处理。 0069 (3-2)。

30、利用(3-1)中整理的地震记录， 提取不同偏移距的地震映像剖面， 而后识别 地震映像剖面上的异常区域。 0070 (3-3)利用(3-1)中整理的地震记录， 计算瑞雷面波相速度剖面， 计算的方法为利 用相邻道相位差法提取多道频散曲线， 然后再利用半波长经验公式得到相速度剖面， 最后 通过分析相对速度剖面圈定异常。 0071 (3-4)联合步骤(3-2)和(3-3)的结果， 即地震映像剖面、 瑞雷波相速度剖面结果， 进行比照确认， 从而既可以得到更精确的异常位置， 也可以排除掉假异常。 0072 所述步骤四中， 针对S波速度结构分层， 联合折射波层析成像， 瑞雷面波频散曲线 反演进行联合反演， 。

31、其过程为： 0073 (4-1)整理原始的地震记录， 按需要重建观测系统， 方便后续处理。 0074 (4-2)利用步骤(4-1)中整理的地震记录， 抽出相遇或追逐相遇观测系统， 然后拾 取初至， 解释并建立初始模型， 再进行走时层析成像来反演地下纵波速度， 并根据速度的分 说明书 5/11 页 8 CN 110687602 A 8 布进行分层。 0075 在该步骤中， 检波器接收到的地震波可以展示成地震波随着时间在不停的振动， 初至就是这些振动在时间轴上第一次起跳的位置。 0076 反演出来的地下速度是一个随着深度变化的函数， 往往变化较快的地方存在地层 性质的变化， 也即存在不同地层的分界。

32、面。 后续会根据这一特点判断进行分层。 0077 (4-3)利用步骤(4-1)中整理的地震记录， 确定该测线上面波要做多少个排列， 确 定偏移距， 然后抽出确定好的排列的面波记录。 0078 (4-4)从步骤(4-3)中层析成像的结果中， 计算出每个排列中心的每层平均纵波速 度与层厚， 将其为先验信息提供给面波进行频散曲线的反演， 从而得到更精确的S波速度结 构。 0079 步骤(3-3)中， 计算瑞雷波相速度剖面， 需要运用多排列正反炮覆盖叠加的地震记 录数据进行， 提取数据过程： 在利用多波联合勘探观测系统采集到的数据中， 确定面波排列 的长度与偏移距， 按该偏移距抽出此排列的正反炮数据，。

33、 然后排列移动半个长列长度， 接着 抽出正反炮数据， 直至到达测线末端， 这些数据组成了后续处理的数据集。 0080 步骤(4-2)中， 抽出数据后， 利用能量比法自动拾取初至， 拾取不准的需要进行手 动调整， 而后利用t0差数法进行解释， 并利用解释的结果建立初始模型， 初始模型是利用解 释的第一层和基岩面的纵波速度， 从第一层到基岩面均匀递增建立的。 0081 步骤(4-3)中， 确定偏移距可以用以下标准： 在非常软的地层(Vs100m/s)， 偏移距 应为15m； 在较软的地层(100m/sVs300m/s)， 偏移距应为510m； 在较硬的地层(300m/s Vs500m/s)， 偏移。

34、距应为1020m； 在非常硬的地层(500m/sVs)， 偏移距应为2040m。 0082 步骤(4-4)中， 瑞雷面波波速与S波波速关系式为： 0083 0084 式中： 为泊松比， 该参数本身代表着地下介质的 “软” 、“硬” 程度， 它本身没有量纲 值， 其值范围一般在00.5。 泊松比越大， 表示地下介质越松软； 泊松比越小， 则表示越刚 硬。 尚未完全固结的第四系地层的泊松比一般在0.350.48左右。 可以看到瑞雷波波速与S 波波速有关， 并且在浅层， 瑞雷波速度与S波速度差距不多， 因此可以用瑞雷波来反演S波速 度结构。 瑞雷波频散曲线的反演方法为带阻尼的最小二乘反演， 比较依赖。

35、初始模型的设定， 作为先验信息参与并约束反演的层析成像信息， 参与的方式是在面波排列中心位置， 按一 定半径搜索每层在该处的平均纵波速度， 然后利用纵波与横波速度转换关系得到该处横波 速度作为初始模型的横波速度， 在中心位置的分层结果作为初始模型的层厚， 从而进行反 演。 0085 所述步骤五中， 需要结合地质资料对处理的结果进行进一步的解释。 在异常体的 探测中， 探测到的高速低速异常体可能有多种解释结果， 这就需要结合前期地质资料进行 约束， 比如在采空区， 探测的低速异常体往往是地下的空洞。 在S波速度分层中， 得到的S波 速度需要进一步进行后续的处理， 比如处理得到地基承载力等岩土参数。

36、从而用于后续的解 释指导工作。 0086 本公开一种新的适用于多波联合勘探的野外观测方法， 能够达到一次作业， 多波 多方法并用的目的。 另外， 针对异常体探测、 S波速度分层这两个勘探任务发明了新的联合 说明书 6/11 页 9 CN 110687602 A 9 勘探技术。 具体而言， 就异常体探测来说， 本发明联合使用联合地震映像法、 瑞雷波相邻道 相位差法提取相速度剖面两种处理方法进行异常体探测工作， 并对结果进行比对约束； 就S 波速度分层来说， 本发明联合折射波层析成像、 面波频散曲线反演两种处理方法进行联合 反演， 从而得到精度更高的结果。 本发明可操作性强， 应用快速， 野外勘探。

37、人员可快速上手。 0087 更为详细的例子： 0088 多波联合勘探技术， 在利用多波联合勘探观测系统采集到数据后， 需要根据勘探 任务采取相应的联合勘探技术进行进一步的处理。 勘探任务是异常体探测的工作流程如图 2。 主要包含以下步骤： 0089 (1)整理原始地震记录， 根据需要重建观测系统， 特别要注意坐标、 偏移距的设置， 设置错误往往会影响后续的工作， 比如地震映像剖面的提取； 0090 (2)利用(1)中整理好的地震记录， 提取不同偏移距的地震映像剖面， 而后识别地 震映像剖面上的异常区域； 0091 后续结合实际的地质资料， 解释该异常可能对应的地质体， 比如高速异常体对应 的土。

38、层中的孤石， 低速异常体对应的破碎带等。 0092 (3)利用(1)中整理好的地震记录， 提取不同偏移距的地震映像剖面， 而后识别地 震映像剖面上的异常区域； 0093 (4)利用(1)中整理好的地震记录， 计算瑞雷面波相速度剖面， 计算的方法为利用 相邻道相位差法提取多道频散曲线， 然后再利用半波长经验公式得到相速度剖面， 最后通 过分析相对速度剖面圈定异常。 0094 (5)比对(3)与(4)的结果， 更精确地确定异常位置。 0095 三层层状介质含低速体模型试算。 0096 模型大小为100m40m， 在埋深6m的地方存在一个大小为8m5m的低速异常体， 模 拟岩溶， 破碎带等地质上常见。

39、的低速体， 其上顶埋深为3.5m， 下定埋深为8.5m， 层状介质第 一层埋深为5m， 第二层埋深为14m。 模型的弹性参数如图3(a)所示。 0097 观测系统为多波联合勘探观测系统， 排列的长度为40m， 偏移距为10m， 震源从0m位 置开始激发， 逐道间距移动， 排列的移动长度为半个排列长。 0098 经上述步骤处理后， 得到的结果如图3(b)， 图3(c)所示。 从图3(c)地震映像看到， 横坐标为4555m这一段范围存在同相轴的错乱， 错乱的同相轴向下偏移。 从图3(b)相速度 剖面图中可以看到， 在横坐标为50m， 埋深为5m处存在一个低速的异常体， 该异常体的x方向 空间位置与。

40、地震映像结果对应的很好。 图中红色虚线的位置即是模型中设定的低速异常体 位置， 可以看到相速度剖面上反映的低速异常较好的吻合了模型中的低速体， 虽然其在纵 向上的位置有所偏移， 但总的来说效果令人满意， 特别是当相速度剖面结果与地震映像的 结果进行比对， 两者的结果都有一致的反映时， 更验证了对异常体探测的结果是可靠的。 两 者比对也能提供更准确的异常体空间位置。 0099 三层层状介质含高速体模型试算。 0100 模型大小为100m40m， 在埋深6m的地方存在一个大小为6m5m的椭圆状高速异 常体， 模拟孤石等地质上常见的高速体， 其上顶埋深为3.5m， 下定埋深为8.5m。 模型的弹性 。

41、参数如图4(a)所示。 0101 观测系统为多波联合勘探观测系统， 排列的长度为40m， 偏移距为10m， 震源从0m位 说明书 7/11 页 10 CN 110687602 A 10 置开始激发， 逐道间距移动， 排列的移动长度为半个排列长。 0102 经上述步骤处理后， 得到的结果如图4(b)， 图4(c)所示。 从图4(c)地震映像看到， 横坐标为4555m这一段范围存在同相轴的错乱， 错乱的同相轴向向偏移。 从图4(b)相速度 剖面图中可以看到， 在横坐标为50m， 埋深为5m处存在一个低速的异常体， 该异常体的x方向 空间位置与地震映像结果对应的很好。 从相速度剖面图中可以看到， 在。

42、横坐标为50m， 埋深 为5m处存在一个低速的异常体， 该异常体的x方向空间位置与地震映像结果对应的很好。 但 是该相速度剖面结果却显示下层的相速度比上层速度高， 这与实际的模型相悖， 推测是由 于该高速异常体孤石的存在导致相邻道相位差法提取频散曲线时， 异常体下方出现了虚假 的低速。 图中红色虚线的位置即是模型中设定的低速异常体位置， 可以看到不管是横向还 是纵向， 相速度剖面上反映的高速异常很好的吻合了模型中的高速体， 另外， 将相速度剖面 结果与地震映像的结果进行比对， 两者的结果都有一致的反映， 两者比对也能提供更准确 的异常体空间位置， 同时这也更验证了联合勘探对高速异常体探测的结果。

43、是可靠的。 0103 勘探任务是S波速度结构分层的工作流程也如图2。 主要包含以下步骤： 0104 (1)整理原始地震记录， 根据需要重建观测系统； 0105 (2)利用(1)中整理的地震记录， 抽出相遇或追逐相遇观测系统， 组织成相应的折 射波数据； 0106 抽出相遇观测系统， 先抽出炮点第一道检波器位置处激发的地震记录， 再抽出炮 点在距最后一道检波器位置处激发的地震记录， 这两个地震记录组合成的数据即为相遇观 测系统采集的数据， 其代表的观测系统即为相遇观测系统。 追逐相遇观测系统也是一样， 只 是在上述相遇观测系统上， 规定第一道检波器到最后一道检波器方向为正向， 分别抽出炮 点在第。

44、一道检波器反向一定长度位置处激发的地震记录以及炮点在最后一道检波器正向 相同长度位置处激发的地震记录， 四个地震记录共同构成追逐相遇观测系统。 0107 (3)利用能量比法自动拾取初至， 并根据经验手动调整自动拾取有问题的初至； 0108 (4)利用t0差数法进行解释工作， 得到覆盖层与基岩面的深度与纵波速度， 再根据 这些信息建立初始模型， 具体的建立方法是速度沿深度递增模型， 遵循公式： 0109 0110 (5)利用SIRT， 即联合迭代重建法进行走时层析成像反演地下纵波速度， 得到速 度-深度剖面； 0111 (6)利用(1)中的地震记录， 确定好偏移距， 该测线要做的面波排列数， 抽。

45、出相应的 面波排列， 并以各个面波排列的中心为中心位置， 按一定的搜索半径， 搜索该位置处的各层 厚度、 每层的平均纵波速度作为后续处理的先验信息； 0112 (7)将(6)中得到的每层层厚， 纵波速度作为先验信息参与进面波频散曲线的反演 中， 具体参与的方法是将层厚作为初始模型， 将纵波速度作为每次正演时候的准确值。 反演 的方法是带阻尼的最小二乘反演； 0113 (8)结合实际地质资料对结果进行解释。 0114 倾斜地层模型试算。 0115 模型大小为120m40m， 一共两层介质， 模拟一个倾斜的地层， 地层左侧界面深度 为5.5m， 右侧界面深度为18.5m。 模型的弹性参数如图5(a。

46、)所示。 说明书 8/11 页 11 CN 110687602 A 11 0116 观测系统为排列长度120m， 炮点从0m开始激发， 间隔20m移动激发， 地面全排列接 收， 然后再抽出相应的地震记录， 面波排列长度为40m， 偏移距为10m。 0117 经上述步骤处理后， 得到的层析成像结果如图5(b)所示。 设置搜索半径为10m， 得 到的先验信息如表1所示。 0118 表1倾斜地层先验信息表 0119 0120 将以上信息作为先验信息进行频散曲线反演得到的结果如图6(a)-图6(d)。 图6 (a)-图6(d)中排列一、 二、 三、 四的横坐标中心位置分别为30m， 50m， 70m，。

47、 90m。 相比于未加先 验信息的反演结果， 施加先验信息的S波速度反演结果更符合真实的模型， 随着排列向右推 进， 虽然幅度和实际有差距， 第二层地层的厚度也逐渐增长， 因为这是一个倾斜地层模型， 这也与模型中逐渐变厚的第一层介质对应的很好。 具体的误差分析如表2。 0121 表2倾斜地层反演误差表 0122 0123 通过误差表可以看到施加先验信息的反演结果精度提升了很多。 相比未加先验信 息的结果， 施加先验信息的反演结果误差总体少了7.2。 这也说明， 施加先验信息的S波速 度结构反演结果往往优于未施加先验信息的结果。 0124 三层介质倾斜地层模型试算。 0125 模型大小为120m。

48、40m， 第一层层状介质埋深为4.8m， 第二层模拟一个倾斜的地 层， 地层左侧界面深度为9.5m， 右侧界面深度为21.4m。 模型的弹性参数如图7(a)所示。 0126 观测系统为排列长度120m， 炮点从0m开始激发， 间隔20m移动激发， 地面全排列接 收， 然后再抽出相应的地震记录， 面波排列长度为40m， 偏移距为10m。 说明书 9/11 页 12 CN 110687602 A 12 0127 经上述步骤处理后， 得到的层析成像结果如图7(b)所示。 设置搜索半径为10m， 得 到的先验信息如表3所示。 0128 表3三层介质倾斜地层先验信息表。 0129 0130 将以上信息作。

49、为先验信息进行频散曲线反演得到的结果如图8(a)-图8(d)。 图8 (a)-图8(d)中排列一、 二、 三、 四的横坐标中心位置分别为30m， 50m， 70m， 90m。 相比于未提供 先验信息的频散曲线反演结果。 并且施加先验信息的反演结果， 其S波速度与反演得到的层 厚度往往更加贴近真实的模型， 这也表明施加先验信息的反演结果更加精确。 相比于前面 较为简单的模型， 该模型的未加先验信息的频散曲线反演很容易收敛到与真实模型差异较 大的极小点上， 导致结果与真实模型相差较大， 施加先验信息后， 反演更容易收敛到与真实 模型差异相近的极小值上， 这也是反演经过精度提高的原因。 对于前面提到。

50、的随着排列的 前进， 瑞雷面波速度与深度曲线显示的深度速度越来越小， 施加先验信息的反演结果也显 示了第二层介质的厚度在逐渐的增厚， 这也与其有所对应， 但是增厚的幅度小于真实模型 中倾斜的地层， 其原因可能是反演的过程往往受多个因素影响， 其结果也是多个影响因素 的共同体现， 真实模型的层厚变化幅度不一定会一致的体现在反演的结果上。 具体的误差 分析如表4。 0131 表4三层介质带倾斜地层反演误差表 0132 说明书 10/11 页 13 CN 110687602 A 13 0133 可以看到下相比于未加先验信息的结果， 施加先验信息的结果的误差下降了16 左右， 引起误差的大头在层厚上，。