CN110187379A - 一种基于tsp法隧道超前地质预报效果的测试方法 - Google Patents
一种基于tsp法隧道超前地质预报效果的测试方法 Download PDFInfo
- Publication number
- CN110187379A CN110187379A CN201910425988.9A CN201910425988A CN110187379A CN 110187379 A CN110187379 A CN 110187379A CN 201910425988 A CN201910425988 A CN 201910425988A CN 110187379 A CN110187379 A CN 110187379A
- Authority
- CN
- China
- Prior art keywords
- density
- wave velocity
- tunnel
- model
- fault
- 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.)
- Pending
Links
- 230000000694 effects Effects 0.000 title claims abstract description 13
- 238000000034 method Methods 0.000 title abstract description 35
- 238000010998 test method Methods 0.000 title abstract description 6
- 238000004088 simulation Methods 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 14
- 230000007480 spreading Effects 0.000 claims abstract description 8
- 239000011435 rock Substances 0.000 claims description 98
- 239000000945 filler Substances 0.000 claims description 37
- 238000001514 detection method Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 9
- 238000003325 tomography Methods 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 210000002706 plastid Anatomy 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 8
- 230000005284 excitation Effects 0.000 description 7
- 238000010276 construction Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000009412 basement excavation Methods 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 4
- 239000004575 stone Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000003874 inverse correlation nuclear magnetic resonance spectroscopy Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000001028 reflection method Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
- 230000005641 tunneling Effects 0.000 description 1
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/01—Measuring or predicting earthquakes
-
- 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/282—Application of seismic models, synthetic seismograms
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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
一种基于TSP法隧道超前地质预报效果的测试方法，以有效提高TSP法预报效果，为相关规范的修订完善和日常大量预报工作的开展指明方向。包括如下步骤：①建立三维空间模型，建立带隧道模型的正演数值模拟方程；将一阶应力‑运动方程进行空间离散和时间离散；将空间离散为网格后，将各个变量按照交错网格的方式赋予结点；利用雷克子波设置为震源模型；②对模型内隧道断面轮廓边界线进行处理，对位于自由边界上的网格进行处理，使其与实际情况相吻合；③对三维空间模型边界进行匹配层处理；④建立不同规模、产状的单个或多个断层破碎带、单个或多个不同规模岩溶的数值模拟模型，将模型带入有限差分数值模拟方程中进行数值模拟计算，根据正演数值模拟建立不良地质体模型；⑤将获得的正演三分量数据进行反演计算，根据反演计算结果和所建立的不良地质体模型的比对分析。
Description
技术领域
本发明涉及隧道超前地质预报，特别涉及一种基于TSP法隧道超前地质预报效果的测试方法。
背景技术
随着国民经济的发展，隧道工程20km、30km的长大隧道不断涌现。如此规模的长大隧道，难免在施工过程中遇到复杂的工程地质及水文地质问题。因隧道工程属于地下建筑工程，其工程地质、水文地质条件复杂多变，在隧道修建过程中，会遇到构造、岩溶、人为坑洞、煤系地层、高地应力等等不良地质体，易引起隧道涌水突泥、有毒有害气体溢出、软岩变形及硬岩岩爆等各种工程问题。
勘察设计阶段，受到经济、工期以及现阶段勘察技术水平的限制，想要在勘察设计阶段仅仅依靠地表的地质勘察工作就准确地查明上述不良地质体的形态特性，并根据这些特性准确地预测出隧道开挖过程中可能会出现的不良地质体的位置、形态、规模等特征是十分困难的。施工阶段，为了进一步探明不良地质体的分布规律，降低隧道施工风险和运营风险，综合超前地质预报技术已经被作为一项工序贯穿隧道掘进的全过程，并且已经体现出了一定的社会效益和经济效益。
我国隧道超前地质预报物探工作中，TSP(Tunnel Seismic Predication)法是地震波反射法中应用最为普遍的一种地球物理探测方法。《铁路隧道超前地质预报技术规程》(Q/CR 9217-2015)中，仅对地震波反射法的探测距离、可被探测的断层倾角和走向做出了规定，对断层的规模简单模糊描述为“可被探测的规模”。纵观各种文献，未见对该方法的有效探测距离、断层倾角、走向、预报分辨率、前后多个不良地质体等情况下预报效果的理论研究成果，更多的是基于大量预报案例进行的总结和统计，缺乏TSP法相应的预报效果数值模拟测试方法。
发明内容
本发明所要解决的技术问题是提供一种基于TSP法隧道超前地质预报效果的测试方法，以有效提高TSP法预报效果，为相关规范的修订完善和日常大量预报工作的开展指明方向。
本发明解决上这技术问题所采用的技术方案如下：
本发明一种基于TSP法隧道超前地质预报效果的测试方法，包括如下步骤：
①建立三维空间模型，并对模型进行交错网格有限差分，建立带隧道模型的正演数值模拟方程；将一阶应力-运动方程进行空间离散和时间离散；将空间离散为网格后，将各个变量按照交错网格的方式赋予结点；利用雷克子波设置为震源模型；
②对模型内隧道断面轮廓边界线进行处理，对位于自由边界上的网格进行处理，使其与实际情况相吻合；
③对三维空间模型边界进行匹配层处理；
④建立不同规模、产状的单个或多个断层破碎带、单个或多个不同规模岩溶的数值模拟模型，将模型带入有限差分数值模拟方程中进行数值模拟计算，获得TSP正演三分量数据，根据正演数值模拟建立不良地质体模型；
⑤将获得的正演三分量数据进行反演计算，根据反演计算结果和所建立的不良地质体模型的比对分析，对预报分辨率、断层可被探测的具体规模大小、多个不良地质体情况下的预报效果、有效探测距离等进行测试。
所述步骤③中，在研究区域边界加入吸收层，厚度设置为10或者20；将速度分解为平行和垂直两个分量；将水平分量的速度按照笛卡尔坐标系的三个方向进行分解；将分解后的分量利用一阶应力-速度方程进行有限差分逼近；将分量带入速度方程和应力方程，得到三维各项同性介质中完美匹配层吸收边界条件应力方程组及速度方程组的交错网格有限差分格式。
本发明的有益效果是，提出了一种专门针对地下隧道三维空间TSP法的三维正反演数值模拟方法，形成了典型不良地质体的正反演图谱，发明了提高TSP法预报效果的一整套数值模拟测试方法，为后期相关规范的修订完善和日常大量预报工作的开展指明了方向；结论准确，通过引入合适的边界处理方法，使计算结果符合地震波的传播规律；计算高效，本发明采用交错网格有限差分法作为计算方法，其计算效率高、易于并行，在模型大小确定的情况下，其计算消耗时间是恒定的；包含齐全的地震波场，信息含量丰富，并且计算效率高而成为了研究地震波场的一种行之有效的方法；测试重点包括有效预报距离、随预报距离远近的预报具体分辨率、不同产状下倾角走向可被预报的极限角度(临界角)、多个不良地质体情况下的预报效果计算，可利用本发明研究隧道超前地质预报中地震波在各种复杂地质情况下的传播特征与规律。
附图说明
本说明书包括如下98幅附图：
图1是应力分量及速度分量分布图；
图2是隧道边界网格示意图；
图3～图48是典型不良地质体模型的平面示意图；
图49～图94正演模型对应数据反演成果图；
图95～图98是实施例的TSP法预报成果图件。
具体实施方式
下面结合附图和实施案例对本发明进一步说明。
我国隧道超前地质预报物探工作中，TSP(Tunnel Seismic Predication)法是地震波反射法中应用最为普遍的一种地球物理探测方法。《铁路隧道超前地质预报技术规程》(Q/CR 9217-2015)中，仅对地震波反射法的探测距离、可被探测的断层倾角和走向做出了规定，对断层的规模简单模糊描述为“可被探测的规模”。纵观各种文献，未见对该方法的有效探测距离、断层倾角、走向、预报分辨率、前后多个不良地质体等情况下预报效果的理论研究成果，更多的是基于大量预报案例进行的总结和统计，缺乏TSP法相应的预报效果数值模拟测试方法。
本发明种基于TSP法隧道超前地质预报效果的测试方法，包括如下步骤：
①建立三维空间模型，并对模型进行交错网格有限差分，建立带隧道模型的正演数值模拟方程；将一阶应力-运动方程进行空间离散和时间离散；参照图1，将空间离散为网格后，将各个变量按照交错网格的方式赋予结点；利用雷克子波设置为震源模型。
②对模型内隧道断面轮廓边界线进行处理，对位于自由边界上的网格进行处理，使其与实际情况相吻合。对于边界形状的描述，因为用于建模的网格是离散的，因此对网格要求有限制。图2是表示在空间4阶差分情况下，对隧道断面轮廓线凸起和凹陷所需要的最少网格数目。因为有限差分的精度是与网格细分程度以及空间差分的阶的有限差分算子相关。
下面对图2边界上的离散灰色网格的七种类型分别进行说明，利用其判定规则在模拟建模的时候进行判别：
(1)H-：表示水平分界面(Horizontal)的网格。判定时，在水平方向，左右都必须有H-，OL-，OR-，IL-或者IR-的灰色网格与之相邻。处理时，规则与平面自由分界面的处理完全一致，保证应力σzz为零。
(2)VL-：表示真空在其正左边(Vertical boundary grid-point with vacuum tothe Light)的网格。判定时，在垂直方向，上下都必须有VL-，IL-或者OL-的灰色网格与之相邻。处理时，其方式与水平分界面网格H-类似，应当保证应力为σxx为零。
(3)IL-：表示真空部分的内角在其左上角(Inner corner with vacuum above toThe Left)的网格点。判定时，在垂直方向，其上方必须有OL-或者VL-的灰色网格与之相邻；在水平方向，其左边必须有H-或者OL-灰色网格与之相邻。处理时，由于σxx与σzz都属于自由边界上，因此应保证其都被设置为零。
(4)OL-：表示真空部分的外角位于其左方(Outer corner grid-point withvacuum to the Left)。判定时，水平方向必须有一个IL-或者H-灰色网格位于其右边且垂直方向必须有IL-或者VL-灰色网格在其下方。处理时，因为σxz位于真空内部，因此必须保证其为零。如果该网格类型与IL-网格类型相邻，那么质子速度分量必须被设置为零，以保证稳定性和获得较为精确的解。
(5)VR-：表示真空部分位于其正右方(Vertical boundary grid-point withvacuum to the Right)。判定时，垂直方向，上下都必须有VR-，IR-或者OR-灰色网格与之相邻。处理的方式与VL-类似。
(6)IR-：表示真空部分的内角位于其右上方(Inner corner with vacuum aboveto The Right)，判定时，在垂直方向，其上方必须有OR-或者VR-的灰色网格与之相邻；在水平方向，其左边必须有H-或者OR-灰色网格与之相邻。处理时，由于σxx与σzz都属于自由边界上，因此应保证其都被设置为零。
(7)OR-：表示真空部分的外角位于其右方(Outer corner grid-point withvacuum to the Right)。判定时，水平方向必须有一个IR-或者H-灰色网格位于其右边且垂直方向必须有IR-或者VR-灰色网格在其下方。处理时，因为σxz位于真空内部，因此必须保证其为零。如果该网格类型与IR-网格类型相邻，那么质子速度分量必须被设置为零，以保证稳定性和获得较为精确的解。
将自由分界面上的网格密度的导数设置为零，那么粒子速度分量在每一次迭代更新以后都会自动为零。不过计算会因为被自由边界网格上方的真空网格消耗一部分产生不必要的计算损失。对于应力的更新，对于VR-，OR-这样的点，其速度分量在真空内，但是其由应力更新的粒子速度不为零，有必要设置其分量为零，然后再用速度更新其应力。这样设置以后，自由边界条件得以全部满足。
③对三维空间模型边界进行匹配层处理。
在研究区域边界加入吸收层，厚度设置为10或者20；将速度分解为平行和垂直两个分量；将水平分量的速度按照笛卡尔坐标系的三个方向进行分解；将分解后的分量利用一阶应力-速度方程进行有限差分逼近；将分量带入速度方程和应力方程，得到三维各项同性介质中完美匹配层吸收边界条件应力方程组及速度方程组的交错网格有限差分格式。
④建立不同规模、产状的单个或多个断层破碎带、单个或多个不同规模岩溶的数值模拟模型，将模型带入有限差分数值模拟方程中进行数值模拟计算，获得TSP正演三分量数据，根据正演数值模拟建立典型不良地质体模型；
⑤将获得的正演三分量数据进行反演计算，根据反演计算结果和所建立的典型不良地质体模型的比对分析，对预报分辨率、断层可被探测的具体规模大小、多个不良地质体情况下的预报效果、有效探测距离等进行测试。
所述步骤④中，所建立典型不良地质体模型包括断层模型、溶洞模型中的至少一种。
所述断层模型至少包括如下断层模型中的一种：
参照图3，断层模型1：断层破碎带宽度5m，位于掌子面前方50～55m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图4，断层模型2：断层破碎带宽度5m，位于掌子面前方100～105m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图5，断层模型3：断层破碎带宽度5m，位于掌子面前方150～155m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图6，断层模型4：断层破碎带宽度10m，位于掌子面前方50～60m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图7，断层模型5：断层破碎带宽度10m，位于掌子面前方100～110m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图8，断层模型6：断层破碎带宽度10m，位于掌子面前方150～160m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；；
参照图9，断层模型7：共有两条断层破碎带，宽度均为5m，分别位于掌子面前方50～55m、100～105m处，均与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图10，断层模型8：共有两条断层破碎带，宽度均为5m，分别位于掌子面前方100～105m、150～155m处，均与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图11，断层模型9：共有两条断层破碎带，宽度均为10m，分别位于掌子面前方50～60m、100～110m处，均与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图12，断层模型10：共有两条断层破碎带，宽度均为10m，分别位于掌子面前方100～110m、150～160m处，均与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图13，断层模型11：断层破碎带宽度为10m，位于掌子面前方50～60m处，走向与隧道轴线正交，倾角45°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图14，断层模型12：断层破碎带宽度为10m，位于掌子面前方50～60m处，走向与隧道轴线正交，倾角35°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图15，断层模型13：断层破碎带宽度为10m，位于掌子面前方50～60m处，走向与隧道轴线正交，倾角25°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图16，断层模型14：断层破碎带宽度为10m，位于掌子面前方50～60m处，与隧道轴线斜交55°，倾角90°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图17，断层模型15：断层破碎带宽度为10m，位于掌子面前方50～60m处，与隧道轴线斜交45°，倾角90°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图18，断层模型16：断层破碎带宽度为10m，位于掌子面前方50～60m处，与隧道轴线斜交35°，倾角90°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图9，断层模型17：断层破碎带宽度为10m，位于掌子面前方50～60m处，与隧道轴线斜交45°，倾角35°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
参照图20，断层模型18：断层破碎带宽度为5m，位于掌子面前方50～55m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
参照图21，断层模型19：断层破碎带宽度为5m，位于掌子面前方100～105m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
参照图22，断层模型20：断层破碎带宽度为5m，位于掌子面前方150～155m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
参照图23，断层模型21：断层破碎带宽度为10m，位于掌子面前方50～60m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
参照图24，断层模型22：断层破碎带宽度为10m，位于掌子面前方100～110m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
参照图25，断层模型23：断层破碎带宽度为10m，位于掌子面前方150～160m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
参照图26，断层模型24：共有2条断层破碎带，宽度均为5m，分别位于掌子面前方50～55m、100m～105m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
参照图27，断层模型25：共有2条断层破碎带，宽度均为5m，分别位于掌子面前方100～105m、150m～155m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
参照图28，断层模型26：共有2条断层破碎带，宽度均为10m，分别位于掌子面前方50～60m、100m～110m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
参照图29，断层模型27：共有2条断层破碎带，宽度均为10m，分别位于掌子面前方100～110m、150m～160m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3。
所述溶洞模型至少包括如下溶洞模型中的一种：
参照图30，溶洞模型1：溶洞为球状，半径为10m，位于掌子面前方40～60m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图31，溶洞模型2：溶洞为球状，半径为10m，位于掌子面前方65～85m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图32，溶洞模型3：溶洞为球状，半径为10m，位于掌子面前方90～110m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3。
参照图33，溶洞模型4：溶洞为球状，半径为10m，位于掌子面前方115～135m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图34，溶洞模型5：溶洞为球状，半径为10m，位于掌子面前方140～160m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图35，溶洞模型6：共2个溶洞，呈球状，半径均为7m，位于掌子面前方43～57m、68～82m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图36，：共2个溶洞，呈球状，半径均为7m，位于掌子面前方68～82m、93～107m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图37，：共2个溶洞，呈球状，半径均为7m，位于掌子面前方93～107m、118～132m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图38，：共2个溶洞，呈球状，半径均为7m，位于掌子面前方118～132m、143～157m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图39，：洞为球状，半径为4m，位于掌子面前方46～54m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图40，：溶洞为球状，半径为4m，位于掌子面前方71～79m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图41，溶洞模型12：溶洞为球状，半径为4m，位于掌子面前方96～104m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图42，溶洞模型13：溶洞为球状，半径为4m，位于掌子面前方121～129m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图43，溶洞模型14：溶洞为球状，半径为4m，位于掌子面前方146～154m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图44，溶洞模型15：溶洞呈椭球状，半轴分别为沿隧道轴线方向x＝10m，水平方向、垂直于轴线y＝5m，竖直方向、垂直于轴线z＝5m，位于掌子面前方40～60m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图45，溶洞模型16：溶洞呈椭球状，半轴分别为沿隧道轴线方向x＝10m，水平方向垂直于轴线y＝5m，竖直方向、垂直于轴线z＝5m，位于掌子面前方65～85m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
参照图46，溶洞模型17：溶洞呈球状，半径为10m，位于掌子面前方40～60m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝2768m/s，VS＝1600m/s，密度ρ＝0.05g/cm3；
参照图47，溶洞模型18：溶洞呈球状，半径为10m，位于掌子面前方65～85m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝2768m/s，VS＝1600m/s，密度ρ＝0.05g/cm3；
参照图48，溶洞模型19：共2个溶洞，呈球状，半径均为7m，位于掌子面前方43～57m、68～82m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝2768m/s，VS＝1600m/s，密度ρ＝0.05g/cm3。
断层模型1～断层模型27的正演数据经TSP数据处理软件，反演结果如图49～图75所示，溶洞模型1～溶洞模型19的正演数据经TSP数据处理软件，反演结果如图76～图94所示。
本发明提出一种专门针对地下隧道三维空间TSP法的三维正反演数值模拟方法，形成了典型不良地质体的正反演图谱，发明了提高TSP法预报效果的一整套数值模拟测试方法，为后期相关规范的修订完善和日常大量预报工作的开展指明了方向；结论准确，通过引入合适的边界处理方法，使计算结果符合地震波的传播规律；计算高效，本发明采用交错网格有限差分法作为计算方法，其计算效率高、易于并行，在模型大小确定的情况下，其计算消耗时间是恒定的；包含齐全的地震波场，信息含量丰富，并且计算效率高而成为了研究地震波场的一种行之有效的方法；测试重点包括有效预报距离、随预报距离远近的预报具体分辨率、不同产状下倾角走向可被预报的极限角度(临界角)、多个不良地质体情况下的预报效果计算，可利用本发明研究隧道超前地质预报中地震波在各种复杂地质情况下的传播特征与规律。
实施例1—破碎带预报
1、工程概况
坡录元隧道位于云桂铁路(广西段)YGZQ-4标段内，隧道正洞起讫里程为DK270+335～DK282+240段，全长11905m。本次预报里程段内，下伏基岩为泥岩夹砂岩、砂岩夹泥岩、砂岩、页岩互层。普遍含钙质，局部夹灰岩透镜体。
2、观测***的布置
在隧道DK274+418的左边墙和右边墙位置分别布置一个地震波信息接收孔，孔径为50mm。在DK274+434～DK274+461段的左边墙位置，按约1.4m的间距布置20个激发孔分别激发地震波，激发孔孔深1.5m左右，孔径约45mm，孔向下倾斜约15，每个激发孔装填的药量为66g。激发孔和接收孔基本保持在同一个高度上。
3、数据处理后用于资料解释的成果图件
坡录元隧道TSP反射层位及物理力学参数二维成果图见图95所示。
4、结论
本次预报时掌子面里程为：DK274+463，预报里程范围为DK274+463～DK274+580段(即掌子面前方117m)。
(1)DK274+463～DK274+480段：围岩较破碎，岩质较硬；
(2)DK274+480～DK274+520段：围岩破碎，岩质较硬，该段围岩节理裂隙强烈发育，含水(呈滴水状～雨淋状水)，受节理裂隙密集带影响，该段围岩稳定性较差，提请施工时该段注意加强支护防止坍塌；
(3)DK274+520～DK274+580段：围岩较破碎，岩质较硬。
5、开挖验证情况
隧道开挖在DK274+500附近发生较大规模的塌方，所幸未造成人员伤亡。
实施案2—岩溶预报
1、工程概况
古家山隧道全长2066m，本次预报里程范围内隧道埋深40～109m，由二叠系茅口栖霞组灰岩组成，属坚硬岩类，裂隙较发育，岩体为中厚层状结构，岩体较完整。地下水位于洞身以下，该段隧道主要位于垂直渗流带内，岩溶以垂直向岩溶裂隙、溶槽、溶沟发育为主，少量水平向岩溶裂隙，地下水以季节性股状及淋雨状出露。
2、观测***的布置
在隧道出口端K19+327处的左边墙和右边墙位置分别布置一个地震波信息接收孔用于安装检波器，孔径为50mm。在K19+312.5～K19+279段的左边墙位置，按约1.5m的间距布置24个激发孔分别激发地震波，激发孔孔深1.5m左右，孔径约45mm，孔向下倾斜约15°，每个激发孔装填的药量为100g。激发孔和接收孔基本保持在同一个高度上。
3、数据处理后用于资料解释的成果图件
图96是2号检波器SV波速度扫描图像，图97是2号检波器SV波深度偏移图像，图98是古家山隧道TSP反射层位及物理力学参数二维成果图。
4、结论
本次预报时检波器里程为K19+327,掌子面(即隧道开挖工作面)里程为K19+266，预报发现在检波器前方65m附近存在一个低速异常区域，在深度偏移图中发现横波能量有衰减并且能通过，但是显示能量衰减不正常。另外结合二维成果图中的K19+259处杨氏模量数值突然降低(代表围岩强度降低)，泊松比数值突然升高(代表围岩流塑性增加)综合分析：推测在K19+259～K19+251段围岩极破碎，岩溶强烈发育，含水，该段存在溶洞。提请施工时采用超前水平钻孔进行核实。成果图件见图96～98所示。
5、开挖验证情况
隧道开挖掌子面到达K19+259时，在掌子面上出现一个大型溶洞。内部充填物为块石土。后经揭示发现该溶洞占据了整个掌子面，但是在K19+259处由于掌子面的左半部分***后留有约70cm厚的岩墙，挡住了溶洞内块石土的继续塌落。
Claims (5)
1.TSP超前地质预报效果测试方法，包括如下步骤：
①建立三维空间模型，并对模型进行交错网格有限差分，建立带隧道模型的正演数值模拟方程；将一阶应力-运动方程进行空间离散和时间离散；将空间离散为网格后，将各个变量按照交错网格的方式赋予结点；利用雷克子波设置为震源模型；
②对模型内隧道断面轮廓边界线进行处理，对位于自由边界上的网格进行处理，使其与实际情况相吻合；
③对三维空间模型边界进行匹配层处理；
④建立不同规模、产状的单个或多个断层破碎带、单个或多个不同规模岩溶的数值模拟模型，将模型带入有限差分数值模拟方程中进行数值模拟计算，获得TSP正演三分量数据，根据正演数值模拟建立典型不良地质体模型；
⑤将获得的正演三分量数据进行反演计算，根据反演计算结果和所建立的典型不良地质体模型的比对分析，对预报分辨率、断层可被探测的具体规模大小、多个不良地质体情况下的预报效果、有效探测距离等进行测试。
2.如权利要求1所述的TSP超前地质预报效果测试方法，其特征是：所述步骤③中，在研究区域边界加入吸收层，厚度设置为10或者20；将速度分解为平行和垂直两个分量；将水平分量的速度按照笛卡尔坐标系的三个方向进行分解；将分解后的分量利用一阶应力-速度方程进行有限差分逼近；将分量带入速度方程和应力方程，得到三维各项同性介质中完美匹配层吸收边界条件应力方程组及速度方程组的交错网格有限差分格式。
3.如权利要求1所述的TSP超前地质预报效果测试方法，其特征是：所述步骤④，所建立典型不良地质体模型包括断层模型、溶洞模型中的至少一种。
4.如权利要求3所述的TSP超前地质预报效果测试方法，其特征是：所述断层模型至少包括如下断层模型中的一种：
断层模型1：断层破碎带宽度5m，位于掌子面前方50～55m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型2：断层破碎带宽度5m，位于掌子面前方100～105m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型3：断层破碎带宽度5m，位于掌子面前方150～155m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型4：断层破碎带宽度10m，位于掌子面前方50～60m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型5：断层破碎带宽度10m，位于掌子面前方100～110m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型6：断层破碎带宽度10m，位于掌子面前方150～160m处，与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；；
断层模型7：共有两条断层破碎带，宽度均为5m，分别位于掌子面前方50～55m、100～105m处，均与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型8：共有两条断层破碎带，宽度均为5m，分别位于掌子面前方100～105m、150～155m处，均与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型9：共有两条断层破碎带，宽度均为10m，分别位于掌子面前方50～60m、100～110m处，均与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型10：共有两条断层破碎带，宽度均为10m，分别位于掌子面前方100～110m、150～160m处，均与隧道成正交关系；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型11：断层破碎带宽度为10m，位于掌子面前方50～60m处，走向与隧道轴线正交，倾角45°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型12：断层破碎带宽度为10m，位于掌子面前方50～60m处，走向与隧道轴线正交，倾角35°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型13：断层破碎带宽度为10m，位于掌子面前方50～60m处，走向与隧道轴线正交，倾角25°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型14：断层破碎带宽度为10m，位于掌子面前方50～60m处，与隧道轴线斜交55°，倾角90°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型15：断层破碎带宽度为10m，位于掌子面前方50～60m处，与隧道轴线斜交45°，倾角90°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型16：断层破碎带宽度为10m，位于掌子面前方50～60m处，与隧道轴线斜交35°，倾角90°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型17：断层破碎带宽度为10m，位于掌子面前方50～60m处，与隧道轴线斜交45°，倾角35°；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4000m/s，VS＝2312m/s，密度ρ＝2.3g/cm3；
断层模型18：断层破碎带宽度为5m，位于掌子面前方50～55m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
断层模型19：断层破碎带宽度为5m，位于掌子面前方100～105m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
断层模型20：断层破碎带宽度为5m，位于掌子面前方150～155m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
断层模型21：断层破碎带宽度为10m，位于掌子面前方50～60m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
断层模型22：断层破碎带宽度为10m，位于掌子面前方100～110m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
断层模型23：断层破碎带宽度为10m，位于掌子面前方150～160m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
断层模型24：共有2条断层破碎带，宽度均为5m，分别位于掌子面前方50～55m、100m～105m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
断层模型25：共有2条断层破碎带，宽度均为5m，分别位于掌子面前方100～105m、150m～155m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
断层模型26：共有2条断层破碎带，宽度均为10m，分别位于掌子面前方50～60m、100m～110m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3；
断层模型27：共有2条断层破碎带，宽度均为10m，分别位于掌子面前方100～110m、150m～160m处，与隧道轴线斜交60°，倾角45°，破碎带含水；围岩纵、横波速度分别为VP＝5200m/s，VS＝3006m/s，密度ρ＝2.65g/cm3，破碎带纵、横波速度分别为VP＝4200m/s，VS＝2246m/s，密度ρ＝2.4g/cm3。
5.如权利要求3所述的TSP超前地质预报效果测试方法，其特征是：所述溶洞模型至少包括如下溶洞模型中的一种：
溶洞模型1：溶洞为球状，半径为10m，位于掌子面前方40～60m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型2：溶洞为球状，半径为10m，位于掌子面前方65～85m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型3：溶洞为球状，半径为10m，位于掌子面前方90～110m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3。
溶洞模型4：溶洞为球状，半径为10m，位于掌子面前方115～135m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型5：溶洞为球状，半径为10m，位于掌子面前方140～160m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型6：共2个溶洞，呈球状，半径均为7m，位于掌子面前方43～57m、68～82m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型7：共2个溶洞，呈球状，半径均为7m，位于掌子面前方68～82m、93～107m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型8：共2个溶洞，呈球状，半径均为7m，位于掌子面前方93～107m、118～132m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型9：共2个溶洞，呈球状，半径均为7m，位于掌子面前方118～132m、143～157m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型10：洞为球状，半径为4m，位于掌子面前方46～54m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型11：溶洞为球状，半径为4m，位于掌子面前方71～79m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型12：溶洞为球状，半径为4m，位于掌子面前方96～104m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型13：溶洞为球状，半径为4m，位于掌子面前方121～129m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型14：溶洞为球状，半径为4m，位于掌子面前方146～154m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型15：溶洞呈椭球状，半轴分别为沿隧道轴线方向x＝10m，水平方向、垂直于轴线y＝5m，竖直方向、垂直于轴线z＝5m，位于掌子面前方40～60m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型16：溶洞呈椭球状，半轴分别为沿隧道轴线方向x＝10m，水平方向垂直于轴线y＝5m，竖直方向、垂直于轴线z＝5m，位于掌子面前方65～85m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝3000m/s，VS＝1604m/s，密度ρ＝1g/cm3；
溶洞模型17：溶洞呈球状，半径为10m，位于掌子面前方40～60m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝2768m/s，VS＝1600m/s，密度ρ＝0.05g/cm3；
溶洞模型18：溶洞呈球状，半径为10m，位于掌子面前方65～85m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝2768m/s，VS＝1600m/s，密度ρ＝0.05g/cm3；
溶洞模型19：共2个溶洞，呈球状，半径均为7m，位于掌子面前方43～57m、68～82m处；围岩纵、横波速度分别为VP＝3500m/s，VS＝2023m/s，密度ρ＝2.6g/cm3，溶洞填充物的纵、横波速度分别为VP＝2768m/s，VS＝1600m/s，密度ρ＝0.05g/cm3。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910425988.9A CN110187379A (zh) | 2019-05-22 | 2019-05-22 | 一种基于tsp法隧道超前地质预报效果的测试方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910425988.9A CN110187379A (zh) | 2019-05-22 | 2019-05-22 | 一种基于tsp法隧道超前地质预报效果的测试方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110187379A true CN110187379A (zh) | 2019-08-30 |
Family
ID=67717089
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910425988.9A Pending CN110187379A (zh) | 2019-05-22 | 2019-05-22 | 一种基于tsp法隧道超前地质预报效果的测试方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110187379A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117849312A (zh) * | 2024-01-04 | 2024-04-09 | 深地科学与工程云龙湖实验室 | 一种基于岩石物理实验预测隧道挖掘后地表塌陷的方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102183790A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 基于时空双变网格的弹性波正演模拟技术 |
| CN102253415A (zh) * | 2011-04-19 | 2011-11-23 | 中国石油大学(华东) | 基于裂缝等效介质模型的地震响应模式建立方法 |
| CN102495434A (zh) * | 2011-11-25 | 2012-06-13 | 成都畅达通地下工程科技发展有限公司 | 地下工程超前地质预报的方法 |
| CN104483702A (zh) * | 2014-12-30 | 2015-04-01 | 中国科学院地质与地球物理研究所 | 一种适用于非均匀运动水体的地震正演模拟方法 |
| US20160034612A1 (en) * | 2014-07-30 | 2016-02-04 | Chevron U.S.A. Inc. | Re-ordered Interpolation and Convolution for Faster Staggered-Grid Processing |
| CN107367772A (zh) * | 2017-08-29 | 2017-11-21 | 西南石油大学 | 一种预测隧道前方不良地质的超前预报方法 |
| CN109188512A (zh) * | 2018-09-17 | 2019-01-11 | 中国石油大学(华东) | 基于非规则扇形网格剖分的起伏隧道空间正演模拟***及方法 |
-
2019
- 2019-05-22 CN CN201910425988.9A patent/CN110187379A/zh active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102183790A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 基于时空双变网格的弹性波正演模拟技术 |
| CN102253415A (zh) * | 2011-04-19 | 2011-11-23 | 中国石油大学(华东) | 基于裂缝等效介质模型的地震响应模式建立方法 |
| CN102495434A (zh) * | 2011-11-25 | 2012-06-13 | 成都畅达通地下工程科技发展有限公司 | 地下工程超前地质预报的方法 |
| US20160034612A1 (en) * | 2014-07-30 | 2016-02-04 | Chevron U.S.A. Inc. | Re-ordered Interpolation and Convolution for Faster Staggered-Grid Processing |
| US9928315B2 (en) * | 2014-07-30 | 2018-03-27 | Chevron U.S.A. Inc. | Re-ordered interpolation and convolution for faster staggered-grid processing |
| CN104483702A (zh) * | 2014-12-30 | 2015-04-01 | 中国科学院地质与地球物理研究所 | 一种适用于非均匀运动水体的地震正演模拟方法 |
| CN107367772A (zh) * | 2017-08-29 | 2017-11-21 | 西南石油大学 | 一种预测隧道前方不良地质的超前预报方法 |
| CN109188512A (zh) * | 2018-09-17 | 2019-01-11 | 中国石油大学(华东) | 基于非规则扇形网格剖分的起伏隧道空间正演模拟***及方法 |
Non-Patent Citations (1)
| Title |
|---|
| 朱夏乐: "隧道地震超前地质预报的三维有限差分正演模拟", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117849312A (zh) * | 2024-01-04 | 2024-04-09 | 深地科学与工程云龙湖实验室 | 一种基于岩石物理实验预测隧道挖掘后地表塌陷的方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Peterson et al. | Applications of algebraic reconstruction techniques to crosshole seismic data | |
| Pan et al. | 3D scene and geological modeling using integrated multi-source spatial data: Methodology, challenges, and suggestions | |
| Agliardi et al. | Structurally-controlled instability, damage and slope failure in a porphyry rock mass | |
| Zhang et al. | Determination of statistical discontinuity persistence for a rock mass characterized by non-persistent fractures | |
| CN101930083B (zh) | 井间电磁波层析成像多孔对联合反演方法 | |
| CN112485823B (zh) | 高效综合超前地质预报方法 | |
| Voeckler et al. | Estimating regional-scale fractured bedrock hydraulic conductivity using discrete fracture network (DFN) modeling | |
| CN116703244B (zh) | 采矿沉陷区治理效果与综合评价方法 | |
| Giustiniani et al. | Reflection seismic sections across the Geothermal Province of Tuscany from reprocessing CROP profiles | |
| CN113376695B (zh) | 一种适用于煤层底板复杂陷落柱的全波形反演方法 | |
| CN110705168A (zh) | 构造应力场的模拟方法 | |
| Petrov et al. | Tectonophysics of hydrothermal ore formation: an example of the Antei Mo–U deposit, Transbaikalia | |
| CN110133748B (zh) | 一种集成碱交代型铀矿深部成矿信息的方法 | |
| CN110187379A (zh) | 一种基于tsp法隧道超前地质预报效果的测试方法 | |
| CN111859687A (zh) | 一种刻画含铀砂层地质结构的混合地质建模方法及*** | |
| Guyoton et al. | Interrelation between induced seismic instabilities and complex geological structure | |
| CN112734076A (zh) | 一种硬岩型铀矿资源大比例尺定位预测方法 | |
| Gao et al. | Site-scale bedrock fracture modeling of a spent fuel reprocessing site based on borehole group in Northwest, China | |
| Nurpeisova et al. | FORECAST CHANGES IN THE GEODYNAMIC REGIME OF GEOLOGICAL ENVIRONMENT DURING LARGE-SCALE SUBSOIL DEVELOPMENT | |
| Lønøy et al. | Paleokarst reservoir architecture derived from forward modelling of modern cave systems-a case study | |
| Chen et al. | Deformation patterns and failure mechanisms of soft-hard-interbedded anti-inclined layered rock slope in Wolong open-pit coal mine | |
| Espin et al. | Stability assessment of a bedding rock slope using Q-slope and seismic tomography: A case study in the Ecuadorian Amazon | |
| Grechishnikova | Natural fracture networks in unconventional reservoirs for improved development strategies: Niobrara play case study | |
| Liu et al. | Detailed shallow structure of the seismogenic fault of the 1976 Ms7. 8 Tangshan earthquake, China | |
| Igbokwe et al. | Geometries, Driving Factors, And Connectivity Of Background Fractures At The Latemar Carbonate Platform (N. Italy): Relevance For Subsurface Reservoir Modelling |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190830 |