CN111929728A - 一种三维三分量超前精细化地质预报方法 - Google Patents
一种三维三分量超前精细化地质预报方法 Download PDFInfo
- Publication number
- CN111929728A CN111929728A CN202010812964.1A CN202010812964A CN111929728A CN 111929728 A CN111929728 A CN 111929728A CN 202010812964 A CN202010812964 A CN 202010812964A CN 111929728 A CN111929728 A CN 111929728A
- Authority
- CN
- China
- Prior art keywords
- wave
- dimensional
- amplitude
- waves
- rock
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000001914 filtration Methods 0.000 claims abstract description 29
- 230000008569 process Effects 0.000 claims abstract description 18
- 238000010586 diagram Methods 0.000 claims abstract description 16
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 229920000642 polymer Polymers 0.000 claims abstract description 10
- 238000010521 absorption reaction Methods 0.000 claims abstract description 8
- 230000000694 effects Effects 0.000 claims abstract description 8
- 239000000463 material Substances 0.000 claims abstract description 8
- 239000011435 rock Substances 0.000 claims description 44
- 238000005422 blasting Methods 0.000 claims description 38
- 238000012545 processing Methods 0.000 claims description 36
- 238000012937 correction Methods 0.000 claims description 28
- 230000005284 excitation Effects 0.000 claims description 20
- 238000009826 distribution Methods 0.000 claims description 15
- 230000004044 response Effects 0.000 claims description 13
- 230000005540 biological transmission Effects 0.000 claims description 11
- 230000035939 shock Effects 0.000 claims description 11
- 238000009412 basement excavation Methods 0.000 claims description 10
- 238000010606 normalization Methods 0.000 claims description 9
- 238000001228 spectrum Methods 0.000 claims description 8
- 230000003321 amplification Effects 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 6
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 6
- 238000010183 spectrum analysis Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 230000005012 migration Effects 0.000 claims description 4
- 238000013508 migration Methods 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 4
- 238000013528 artificial neural network Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000005553 drilling Methods 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 6
- 239000007822 coupling agent Substances 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 5
- 238000004880 explosion Methods 0.000 abstract description 5
- 238000013277 forecasting method Methods 0.000 abstract description 2
- 230000009471 action Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000011202 physical detection method Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000001028 reflection method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种三维三分量超前精细化地质预报方法，属于隧道工程技术领域，首先建立观测***，然后振幅均衡，接着进行直达波拾取，然后提取纵、横波波速，然后干扰去除、滤波处理，接着制备二维纵、横波速度与反射界面图和三维纵波速度与反射界面图，最后对探测成果进行解释，其中采用高分子凝胶作耦合剂，孔口采用高吸收衰减材料封堵，可以压制干扰波，保证地震波信号，同时采用回线开路触发计时方式，把回线绑扎在震源上，震源***产生振动传播的同时炸断回线触发采集，与***产生过程中的***延迟没有关系，在绑扎回线不脱落震源的条件下基本上不存在触发时间误差，通过二维和三维的精细化地质预报，能达到隧道动态设计施工、节约成本的效果。
Description
技术领域
本发明涉及隧道工程技术领域，特别涉及一种三维三分量超前精细化地质预报方法。
背景技术
超前地质预报(TGP)或隧道超前地质预报，是在隧道开挖时，对掌子面前方及其周边(主要是铁路隧道)的围岩与地层情况做出超前预报。常用的隧道超前地质预报方法主要分为地质方法和物探方法，其中物探方法是根据隧道掌子面前方岩体的物性差异而开展的数据信息化采集、处理和分析的物理探测方法。地震波反射法是常用的超前地质预报物探方法之一，通过激发地震波并接受地震反射波的方法来推演反射界面的位置、地层波速等参数，进而起到探测隧道掌子面前方围岩地质情况的作用，具有分辨率高、探测深度大、对构造反应灵敏的特点。
目前，常用的地下工程超前地质预报的划分还没有统一的标准，多数均按物探的探测距离划分为长期、短期预报两个阶段，有的虽然提出了宏观超前地质预报概念，但对其理论、步骤和方法的研究较少，实际应用的实例较少，在超前地质预报技术综合应用中，仅以两种或两种先进物探手段进行超前探测后就做出超前地质预报的现象较为常见，在众多的隧道施工中，工作人员不断的优化组合，但是现有的超前精细化地质预报还是存在预测精准度一般，对优化设计施工的效果一般。
发明内容
本发明的目的在于提供一种三维三分量超前精细化地质预报方法，具有能达到隧道动态设计施工、节约成本的优点，以解决上述背景技术中提出的问题。
为实现上述目的，本发明提供如下技术方案：
一种三维三分量超前精细化地质预报方法，包括如下步骤：
S1：建立观测***，即选取超前预报***仪器，并在开展隧道超前地质预报工作前，需要在隧道开挖方向洞壁钻孔；
S2：振幅均衡，野外原始记录或水平叠加剖面上，往往会出现浅、中、深层或道间反射能量差异较大，如不做处理就不便于显示或影响叠加效果；
S3：直达波拾取，进行直达波拾取时，对观测***的检测到的***震动波与普通震动波在岩石中的传输速度进行线性动校正，再进行直达波拾取；
S4：提取纵、横波波速，计算机对采集地震记录进行数据长度设置，调零和平均振幅谱处理，使波形信号处于同一量级水平，进行波形对比和频谱分析；利用噪音信号与有效信号频带响应值差异，设定带通滤波的上限和下限值，采用振幅补偿等非线性放大技术对信号进行处理；利用速度分析结果进行地震剖面的时深转换，获得偏移处理后的偏移结果；
S5：干扰去除、滤波处理；
S6：二维纵、横波速度与反射界面图；
S7：三维纵波速度与反射界面图；
S8；探测成果解释，解释主要以纵波波速和纵波反射界面为主。
优选的，所述S1中为采集到较好的TGP原始数据，要求隧道内能产生震动的一切工序均停工，检波器安装时采用高分子凝胶耦合，通过采用高吸声材料封堵接收孔来压制声波噪声；在隧道开挖方向右壁预设均匀分布的18个激发孔位，和1个接收孔位，激发孔间距为1.5m，最小偏移距为20m，接收孔和激发孔高度均约1.2m，孔径均为50mm，孔深均约2.0m。
优选的，所述S2中的振幅均衡过程包括如下步骤：
S201导入三维***震动观测***数据；
S202按照隧道面元位置，对三维***震动观测***数据中的激发孔位与接收孔位信息进行重新排序；
S203计算三维***震动***的每一组炮检对的振幅系数，其中将每个接收孔位的振幅乘以最大振幅平均值与最大振幅绝对值比值的平方根，得到均衡后振幅；
S204对计算出的每一组炮检对的振幅系数进行累加；
S205绘制面元振幅分布图；
S206对绘制的面元振幅分布图进行统计分析。
优选的，所述S4中对每一个接收孔位接收到的***震动反射波进行预处理，预处理包括对反射波数据中的横波与纵波进行时间校正，得到预处理后的反射波模拟数据，以及对应每一个面元的横波速度和纵波速度除以归一化因子进行归一化处理，再将预处理的反射波模拟数据与归一化处理的横波速度和纵波速度一起构成标签数据并输入专用神经网络进行处理，得到二次处理结果，再将二次处理结果乘以归一化因子，得到最终的横波速度和纵波速度。
优选的，所述S5中干扰去除、滤波处理过程包括计算机对经过处理的横波数据和纵波数据进行数据长度设置、调零和平均振幅谱处理，使横波信号与纵波信号处于同一量级水平，然后进行波形对比和频谱分析，利用噪音信号与有效信号频带响应值差异，设定带通滤波的上限和下限值，采用振幅补偿等非线性放大技术对信号进行处理。
优选的，所述S6中，***震动波以大角度入射透过高阻抗层，达到下伏的低阻抗层即使能量比较的弱，在发生广角反射时，在底面也能接受到较强的信号，在数据处理中，为了保留广角反射波，在存在广角反射的区域，不能首先切除初至折射波，而是在校正后的的地震波上进行频率波数域滤波、二维滤波或线性干扰剔除发来压制相关的干扰，提高信噪比，获得较好的叠加剖面。
优选的，所述S7中三维纵波速度与反射界面图在布置时需要去除掌子面或者隧道侧面壁出现塌陷的区域，若掌子面或者隧道侧面壁出现塌陷，则此处的反射界面图通过相邻的图面平滑连接。
优选的，所述S8中的探测成果解释，包括二维解释、三维解释和预测结果综合预报。
优选的，所述S3中，对观测***的检测到的***震动波与普通震动波在岩石中的传输速度进行线性动校正包括：
对所述岩石进行结构判断，判断所述岩石是否为分布均匀的介质；
当所述岩石为分布均匀介质时，根据所述***震动波与普通震动波在岩石中深度以及所述岩石的组成确定动校正量，然后依据动校正量来消除差值进行线性动校正；
当所述岩石为分布不均匀介质时，将所述岩石按照结构分布分成多段均匀介质，分别确定每段岩石介质的动校正量，然后对***震动波与普通震动波在岩石中的传输速度进行分段线性动校正。
优选的，所述设定带通滤波的上限和下限值时，通过下述公式确定带通滤波的上限和下限值：
f＝f0-f0*exp(lg(g+h)*ln10)÷exp(ln2)
上述公式中，F为所述带通滤波的上限，f为所述带通滤波的下限，f0为有效信号的平均频率，g为所述噪音信号的频带响应值，h为所述有效信号频带响应值，ln和lg均为对数函数，e为自然常数。
与现有技术相比，本发明的有益效果是：本发明提出的三维三分量超前精细化地质预报方法，首先建立观测***，然后振幅均衡，接着进行直达波拾取，然后提取纵、横波波速，然后干扰去除、滤波处理，接着制备二维纵、横波速度与反射界面图和三维纵波速度与反射界面图，最后对探测成果进行解释，其中采用高分子凝胶作耦合剂，孔口采用高吸收衰减材料封堵，可以压制干扰波，保证地震波信号，同时采用回线开路触发计时方式，把回线绑扎在震源上，震源***产生振动传播的同时炸断回线触发采集，与***产生过程中的***延迟没有关系，在绑扎回线不脱落震源的条件下基本上不存在触发时间误差，通过二维和三维的精细化地质预报，能达到隧道动态设计施工、节约成本的效果。
附图说明
图1为本发明的超前精细化地质预报流程框图；
图2为本发明的建立观测***的原理示意图；
图3为本发明的振幅均衡流程框图；
图4为本发明的提取纵、横波波速流程框图；
图5为本发明的纵、横波反射界面俯视示意图；
图6为本发明的三维纵波速度与反射界面图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅附图1-图6：一种三维三分量超前精细化地质预报方法，包括如下步骤：
第一步：建立观测***，即先选取合适的超前预报***仪器，其中TGP超前地质预报***采用三分量速度型检波器，三分量速度型检波器具有高灵敏度、高指向性、与隧道围岩传播地震波频率相适应的宽频带特性，最大程度地保证了地震波波形的完整和纵横波信息的丰富明确性，在开展隧道超前地质预报工作前，需要在隧道开挖方向洞壁钻孔，其中包括激发孔和接收孔，激发孔相比接收孔更加靠近掌子面，激发孔的数量要远多于接收孔的数量，这样通过接收孔就能接收到***震动的反射波，其中在隧道开挖方向右壁预设均匀分布的18个激发孔位，和1个接收孔位，激发孔间距为1.5m，最小偏移距为20m，接收孔和激发孔高度均约1.2m，孔径均为50mm，孔深均约2.0m，其中为采集到较好的TGP原始数据，要求隧道内能产生震动的一切工序均停工，检波器安装时采用高分子凝胶耦合，通过采用高吸声材料封堵接收孔来压制声波噪声，即TGP超前地质预报***的检波器采用高分子凝胶作耦合剂，三分量检波器通过高分子凝胶直接与钻孔围岩接触，地震信号通过软线引出与仪器连接，孔口采用高吸收衰减材料封堵，起到压制干扰波，保证地震波信号的作用，同时TGP超前地质预报***采用回线开路触发计时方式，把回线绑扎在震源上，震源***产生振动传播的同时炸断回线触发采集，与***产生过程中的***延迟没有关系，在绑扎回线不脱落震源的条件下基本上不存在触发时间误差。
第二步：振幅均衡，野外原始记录或水平叠加剖面上，往往会出现浅、中、深层或道间反射能量差异较大，如不做处理就不便于显示或影响叠加效果，其中振幅均衡过程包括如下步骤：首先导入三维***震动观测***数据，其中三维地***震动观测***数据是指***震动数据采集中激发孔和接收孔的布设情况，接着按照隧道面元位置，对三维***震动观测***数据中的激发孔位与接收孔位信息进行重新排序，重新排序包括计算三维***震动观测***的每一组炮检对的中点坐标；将相同的中点坐标的炮检对放在一起组成一个面元，按照面元位置，对三维***震动观测***数据的炮点与检波点(激发孔和接收孔)信息进行重新排序，这里炮点与检波点信息可包括三维观测***炮检对的属性分布，例如，覆盖次数、炮检距和方位角等，接着计算三维***震动***的每一组炮检对的振幅系数，其中将每个接收孔位的振幅乘以最大振幅平均值与最大振幅绝对值比值的平方根，得到均衡后振幅，接着对计算出的每一组炮检对的振幅系数进行累加，接着绘制面元振幅分布图，最后对绘制的面元振幅分布图进行统计分析。
第三步：直达波拾取，进行直达波拾取时，对观测***的检测到的***震动波与普通震动波在岩石中的传输速度进行线性动校正，再进行直达波拾取。
第四步：提取纵、横波波速，计算机对采集地震记录进行数据长度设置，调零和平均振幅谱处理，使波形信号处于同一量级水平，进行波形对比和频谱分析；利用噪音信号与有效信号频带响应值差异，设定带通滤波的上限和下限值，采用振幅补偿等非线性放大技术对信号进行处理，利用速度分析结果进行地震剖面的时深转换，获得偏移处理后的偏移结果；其中每一个接收孔位接收到的***震动反射波进行预处理，预处理包括对反射波数据中的横波与纵波进行时间校正，得到预处理后的反射波模拟数据，以及对应每一个面元的横波速度和纵波速度除以归一化因子进行归一化处理，再将预处理的反射波模拟数据与归一化处理的横波速度和纵波速度一起构成标签数据并输入专用神经网络进行处理，得到二次处理结果，再将二次处理结果乘以归一化因子，得到最终的横波速度和纵波速度。
第五步：干扰去除、滤波处理，干扰去除、滤波处理过程包括计算机对经过处理的横波数据和纵波数据进行数据长度设置、调零和平均振幅谱处理，使横波信号与纵波信号处于同一量级水平，然后进行波形对比和频谱分析，利用噪音信号与有效信号频带响应值差异，设定带通滤波的上限和下限值，采用振幅补偿等非线性放大技术对信号进行处理。
第六步：二维纵、横波速度与反射界面图，其中***震动波以大角度入射透过高阻抗层，达到下伏的低阻抗层即使能量比较的弱，在发生广角反射时，在底面也能接受到较强的信号，在数据处理中，为了保留广角反射波，在存在广角反射的区域，不能首先切除初至折射波，而是在校正后的的地震波上进行频率波数域滤波、二维滤波或线性干扰剔除发来压制相关的干扰，提高信噪比，获得较好的叠加剖面。
第七步：三维纵波速度与反射界面图，其中三维纵波速度与反射界面图在布置时需要去除掌子面或者隧道侧面壁出现塌陷的区域，若掌子面或者隧道侧面壁出现塌陷，则此处的反射界面图通过相邻的图面平滑连接。
第七步；探测成果解释，解释主要以纵波波速和纵波反射界面为主，纵波波速和纵波反射界面包括二维解释、三维解释和预测结果综合预报。
本发明提出的三维三分量超前精细化地质预报方法，首先建立观测***，开展隧道超前地质预报工作前，需要在隧道开挖方向洞壁钻孔，其中包括激发孔和接收孔，激发孔相比接收孔更加靠近掌子面，激发孔的数量要远多于接收孔的数量，TGP超前地质预报***的检波器采用高分子凝胶作耦合剂，三分量检波器通过高分子凝胶直接与钻孔围岩接触，地震信号通过软线引出与仪器连接，孔口采用高吸收衰减材料封堵，起到压制干扰波，保证地震波信号的作用，同时TGP超前地质预报***采用回线开路触发计时方式，把回线绑扎在震源上，震源***产生振动传播的同时炸断回线触发采集，与***产生过程中的***延迟没有关系，在绑扎回线不脱落震源的条件下基本上不存在触发时间误差，然后振幅均衡，野外原始记录或水平叠加剖面上，往往会出现浅、中、深层或道间反射能量差异较大，如不做处理就不便于显示或影响叠加效果，其中振幅均衡过程首先导入三维***震动观测***数据，然后接着按照隧道面元位置，对三维***震动观测***数据中的激发孔位与接收孔位信息进行重新排序，接着计算三维***震动***的每一组炮检对的振幅系数，接着对计算出的每一组炮检对的振幅系数进行累加，接着绘制面元振幅分布图，最后对绘制的面元振幅分布图进行统计分析；振幅均衡后进行直达波拾取，进行直达波拾取时，对观测***的检测到的***震动波与普通震动波在岩石中的传输速度进行线性动校正，再进行直达波拾取；然后提取纵、横波波速，计算机对采集地震记录进行数据长度设置，调零和平均振幅谱处理，使波形信号处于同一量级水平，进行波形对比和频谱分析；然后干扰去除、滤波处理，接着制备二维纵、横波速度与反射界面图和三维纵波速度与反射界面图，最后对探测成果进行解释，解释主要以纵波波速和纵波反射界面为主，纵波波速和纵波反射界面包括二维解释、三维解释和预测结果综合预报。
综上所述，本发明提出的三维三分量超前精细化地质预报方法，首先建立观测***，然后振幅均衡，接着进行直达波拾取，然后提取纵、横波波速，然后干扰去除、滤波处理，接着制备二维纵、横波速度与反射界面图和三维纵波速度与反射界面图，最后对探测成果进行解释，其中采用高分子凝胶作耦合剂，孔口采用高吸收衰减材料封堵，可以压制干扰波，保证地震波信号，同时采用回线开路触发计时方式，把回线绑扎在震源上，震源***产生振动传播的同时炸断回线触发采集，与***产生过程中的***延迟没有关系，在绑扎回线不脱落震源的条件下基本上不存在触发时间误差，通过二维和三维的精细化地质预报，能达到隧道动态设计施工、节约成本的效果。
本发明的另一个具体实施例中，所述S3中，对观测***的检测到的***震动波与普通震动波在岩石中的传输速度进行线性动校正包括：
对所述岩石进行结构判断，判断所述岩石是否为分布均匀的介质；
当所述岩石为分布均匀介质时，根据所述***震动波与普通震动波在岩石中深度以及所述岩石的组成确定动校正量，然后依据动校正量来消除差值进行线性动校正；
当所述岩石为分布不均匀介质时，将所述岩石按照结构分布分成多段均匀介质，分别确定每段岩石介质的动校正量，然后对***震动波与普通震动波在岩石中的传输速度进行分段线性动校正。
有益效果：上述技术方案对观测***的检测到的***震动波与普通震动波在岩石中的传输速度进行线性动校正时，根据所述岩石的结构是否为均匀的介质来确定动校正量，而且所述动校正量与所述***震动波与普通震动波在岩石中深度以及所述岩石的组成还有关，通过上述技术方案对观测***的检测到的***震动波与普通震动波在岩石中的传输速度进行线性动校正时根据岩石的构造进行校正可使得计算的动校正量更加准确，进而使得校正效果更佳。
所述设定带通滤波的上限和下限值时，通过下述公式确定带通滤波的上限和下限值：
f＝f0-f0*exp(lg(g+h)*ln10)÷exp(ln2)
上述公式中，F为所述带通滤波的上限，f为所述带通滤波的下限，f0为有效信号的平均频率，g为所述噪音信号的频带响应值，h为所述有效信号频带响应值，ln和lg均为对数函数，e为自然常数。
有益效果：在上述技术方案中，通过设置带通滤波的上限值和下限值使得带通滤波过滤掉高频信号和低频信号，保留下来存在于上限值和下限值之间的中频信号，使得过滤后的信号频率比较缓和，有效降低了高频信号和低频信号对信号整体带来的偏差影响，此外，确定设置带通滤波的上限值和下限值是在有效信号的平均频率的基础上充分考虑到噪音信号的频带响应值和有效信号频带响应值，使得通过带通滤波的波中混入的干扰因素较少，可以有效降低信号信息的偏差率，使得信号能够准确传达所传递的信息。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种三维三分量超前精细化地质预报方法，其特征在于，包括如下步骤：
S1：建立观测***，即选取超前预报***仪器，并在开展隧道超前地质预报工作前，需要在隧道开挖方向洞壁钻孔；
S2：振幅均衡，野外原始记录或水平叠加剖面上，往往会出现浅、中、深层或道间反射能量差异较大，如不做处理就不便于显示或影响叠加效果；
S3：直达波拾取，进行直达波拾取时，对观测***的检测到的***震动波与普通震动波在岩石中的传输速度进行线性动校正，再进行直达波拾取；
S4：提取纵、横波波速，计算机对采集地震记录进行数据长度设置，调零和平均振幅谱处理，使波形信号处于同一量级水平，进行波形对比和频谱分析；利用噪音信号与有效信号频带响应值差异，设定带通滤波的上限和下限值，采用振幅补偿等非线性放大技术对信号进行处理；利用速度分析结果进行地震剖面的时深转换，获得偏移处理后的偏移结果；
S5：干扰去除、滤波处理；
S6：二维纵、横波速度与反射界面图；
S7：三维纵波速度与反射界面图；
S8：探测成果解释，解释主要以纵波波速和纵波反射界面为主。
2.根据权利要求1所述的一种三维三分量超前精细化地质预报方法，其特征在于，所述S1中为采集到较好的TGP原始数据，要求隧道内能产生震动的一切工序均停工，检波器安装时采用高分子凝胶耦合，通过采用高吸声材料封堵接收孔来压制声波噪声；在隧道开挖方向右壁预设均匀分布的18个激发孔位，和1个接收孔位，激发孔间距为1.5m，最小偏移距为20m，接收孔和激发孔高度均约1.2m，孔径均为50mm，孔深均约2.0m。
3.根据权利要求1所述的一种三维三分量超前精细化地质预报方法，其特征在于，所述S2中的振幅均衡过程包括如下步骤：
S201导入三维***震动观测***数据；
S202按照隧道面元位置，对三维***震动观测***数据中的激发孔位与接收孔位信息进行重新排序；
S203计算三维***震动***的每一组炮检对的振幅系数，其中将每个接收孔位的振幅乘以最大振幅平均值与最大振幅绝对值比值的平方根，得到均衡后振幅；
S204对计算出的每一组炮检对的振幅系数进行累加；
S205绘制面元振幅分布图；
S206对绘制的面元振幅分布图进行统计分析。
4.根据权利要求1所述的一种三维三分量超前精细化地质预报方法，其特征在于，所述S4中对每一个接收孔位接收到的***震动反射波进行预处理，预处理包括对反射波数据中的横波与纵波进行时间校正，得到预处理后的反射波模拟数据，以及对应每一个面元的横波速度和纵波速度除以归一化因子进行归一化处理，再将预处理的反射波模拟数据与归一化处理的横波速度和纵波速度一起构成标签数据并输入专用神经网络进行处理，得到二次处理结果，再将二次处理结果乘以归一化因子，得到最终的横波速度和纵波速度。
5.根据权利要求1所述的一种三维三分量超前精细化地质预报方法，其特征在于，所述S5中干扰去除、滤波处理过程包括计算机对经过处理的横波数据和纵波数据进行数据长度设置、调零和平均振幅谱处理，使横波信号与纵波信号处于同一量级水平，然后进行波形对比和频谱分析，利用噪音信号与有效信号频带响应值差异，设定带通滤波的上限和下限值，采用振幅补偿等非线性放大技术对信号进行处理。
6.根据权利要求1所述的一种三维三分量超前精细化地质预报方法，其特征在于，所述S6中，***震动波以大角度入射透过高阻抗层，达到下伏的低阻抗层即使能量比较的弱，在发生广角反射时，在底面也能接受到较强的信号，在数据处理中，为了保留广角反射波，在存在广角反射的区域，不能首先切除初至折射波，而是在校正后的的地震波上进行频率波数域滤波、二维滤波或线性干扰剔除发来压制相关的干扰，提高信噪比，获得较好的叠加剖面。
7.根据权利要求1所述的一种三维三分量超前精细化地质预报方法，其特征在于，所述S7中三维纵波速度与反射界面图在布置时需要去除掌子面或者隧道侧面壁出现塌陷的区域，若掌子面或者隧道侧面壁出现塌陷，则此处的反射界面图通过相邻的图面平滑连接。
8.根据权利要求1所述的一种三维三分量超前精细化地质预报方法，其特征在于，所述S8中的探测成果解释，包括二维解释、三维解释和预测结果综合预报。
9.根据权利要求1所述的一种三维三分量超前精细化地质预报方法，其特征在于，所述S3中，对观测***的检测到的***震动波与普通震动波在岩石中的传输速度进行线性动校正包括：
对所述岩石进行结构判断，判断所述岩石是否为分布均匀的介质；
当所述岩石为分布均匀介质时，根据所述***震动波与普通震动波在岩石中深度以及所述岩石的组成确定动校正量，然后依据动校正量来消除差值进行线性动校正；
当所述岩石为分布不均匀介质时，将所述岩石按照结构分布分成多段均匀介质，分别确定每段岩石介质的动校正量，然后对***震动波与普通震动波在岩石中的传输速度进行分段线性动校正。
10.根据权利要求5所述的一种三维三分量超前精细化地质预报方法，其特征在于，所述设定带通滤波的上限和下限值时，通过下述公式确定带通滤波的上限和下限值：
f＝f0-f0*exp(lg(g+h)*ln10)÷exp(ln2)
上述公式中，F为所述带通滤波的上限，f为所述带通滤波的下限，f0为有效信号的平均频率，g为所述噪音信号的频带响应值，h为所述有效信号频带响应值，ln和lg均为对数函数，e为自然常数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010812964.1A CN111929728A (zh) | 2020-08-13 | 2020-08-13 | 一种三维三分量超前精细化地质预报方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010812964.1A CN111929728A (zh) | 2020-08-13 | 2020-08-13 | 一种三维三分量超前精细化地质预报方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111929728A true CN111929728A (zh) | 2020-11-13 |
Family
ID=73311686
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010812964.1A Withdrawn CN111929728A (zh) | 2020-08-13 | 2020-08-13 | 一种三维三分量超前精细化地质预报方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111929728A (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112305591A (zh) * | 2020-10-10 | 2021-02-02 | 中国地质大学（北京） | 隧道超前地质预报方法、计算机可读存储介质 |
| CN113126161A (zh) * | 2021-05-27 | 2021-07-16 | 深圳市岩土综合勘察设计有限公司 | 一种基于冲击波预测溶洞埋深和尺寸的方法及*** |
| CN113484910A (zh) * | 2021-06-28 | 2021-10-08 | 西安建筑科技大学 | 基于地震干涉法的隧道超前地质预报方法及*** |
| CN115932945A (zh) * | 2022-10-14 | 2023-04-07 | 扎赉诺尔煤业有限责任公司 | 掘进巷道剩余煤厚地震多波多分量的探测方法 |
| US11668848B2 (en) | 2021-06-24 | 2023-06-06 | Saudi Arabian Oil Company | Method and system for seismic imaging using S-wave velocity models and machine learning |
| CN117706614A (zh) * | 2023-12-19 | 2024-03-15 | 四川中铁二院环保科技有限公司 | 一种隧道超前地质预报***及其超前地质预报方法 |
-
2020
- 2020-08-13 CN CN202010812964.1A patent/CN111929728A/zh not_active Withdrawn
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112305591A (zh) * | 2020-10-10 | 2021-02-02 | 中国地质大学（北京） | 隧道超前地质预报方法、计算机可读存储介质 |
| CN112305591B (zh) * | 2020-10-10 | 2022-04-29 | 中国地质大学（北京） | 隧道超前地质预报方法、计算机可读存储介质 |
| CN113126161A (zh) * | 2021-05-27 | 2021-07-16 | 深圳市岩土综合勘察设计有限公司 | 一种基于冲击波预测溶洞埋深和尺寸的方法及*** |
| US11668848B2 (en) | 2021-06-24 | 2023-06-06 | Saudi Arabian Oil Company | Method and system for seismic imaging using S-wave velocity models and machine learning |
| CN113484910A (zh) * | 2021-06-28 | 2021-10-08 | 西安建筑科技大学 | 基于地震干涉法的隧道超前地质预报方法及*** |
| CN113484910B (zh) * | 2021-06-28 | 2024-03-01 | 西安建筑科技大学 | 基于地震干涉法的隧道超前地质预报方法及*** |
| CN115932945A (zh) * | 2022-10-14 | 2023-04-07 | 扎赉诺尔煤业有限责任公司 | 掘进巷道剩余煤厚地震多波多分量的探测方法 |
| CN115932945B (zh) * | 2022-10-14 | 2024-04-02 | 扎赉诺尔煤业有限责任公司 | 掘进巷道剩余煤厚地震多波多分量的探测方法 |
| CN117706614A (zh) * | 2023-12-19 | 2024-03-15 | 四川中铁二院环保科技有限公司 | 一种隧道超前地质预报***及其超前地质预报方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111929728A (zh) | 一种三维三分量超前精细化地质预报方法 | |
| CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
| CN108983284B (zh) | 一种适用于海上斜缆数据的f-p域鬼波压制方法 | |
| WO2017024523A1 (zh) | 一种射线弹性参数的反演方法 | |
| WO2009081210A1 (en) | Method of and apparatus for exploring a region below a surface of the earth | |
| CN102176054B (zh) | 近地表综合信息处理解释方法 | |
| CN109188520B (zh) | 薄储层厚度预测方法及装置 | |
| CN107544087B (zh) | 一种测量近地表地层品质因子的方法及装置 | |
| CN116520419B (zh) | 一种热流体裂缝通道识别方法 | |
| CN106526678A (zh) | 一种反射声波测井的波场分离方法及装置 | |
| CN104330826A (zh) | 一种去除复杂地表条件下多种噪音的方法 | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| CN115373023A (zh) | 一种基于地震反射和车辆噪声的联合探测方法 | |
| CN105298482B (zh) | 一种井旁构造反射波与井壁直达波的分离方法 | |
| CN109188528A (zh) | 井间弹性波层析成像***和方法 | |
| Nardone et al. | Shallow velocity structure and site effects at Mt. Vesuvius, Italy, from HVSR and array measurements of ambient vibrations | |
| Satoh et al. | Nonlinear behaviour of scoria soil sediments evaluated from borehole records in eastern Shizuoka prefecture, Japan | |
| CN113514889B (zh) | 一种提升海洋深反射地震数据中低频信号能量的处理方法 | |
| CN111691876B (zh) | 一种利用声波测井对邻井成像的方法、装置及存储介质 | |
| CN111060969B (zh) | 一种井控q补偿方法 | |
| CN103940568B (zh) | 一种检测坚硬顶板卸压***效果的方法 | |
| CN113296149A (zh) | 一种基于微动探测技术的滑坡体地层分布情况探测方法 | |
| Yang et al. | Revealing the shallow soil structure of the Yigong Lake in the Tibetan Plateau using a portable distributed acoustic sensing interrogator | |
| Glassmoyer et al. | Source parameters and effects of bandwidth and local geology on high-frequency ground motions observed for aftershocks of the northeastern Ohio earthquake of 31 January 1986 | |
| Zheng et al. | Seismic characterization of reservoirs with variable fracture spacing by double focusing Gaussian beams |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20201113 |