CN107180135A - 一种洞穴型岩溶岩体数值模型构建方法 - Google Patents
一种洞穴型岩溶岩体数值模型构建方法 Download PDFInfo
- Publication number
- CN107180135A CN107180135A CN201710378892.2A CN201710378892A CN107180135A CN 107180135 A CN107180135 A CN 107180135A CN 201710378892 A CN201710378892 A CN 201710378892A CN 107180135 A CN107180135 A CN 107180135A
- Authority
- CN
- China
- Prior art keywords
- solution cavity
- karst
- karstification
- msub
- mrow
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
- G01V9/007—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00 by detecting gases or particles representative of underground layers at or near the surface
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种洞穴型岩溶岩体有限元数值模型构建方法，包括：根据工程勘察与物探确定：体岩溶率、溶洞位置范围、溶洞大小分布规律及其概率密度函数；基于蒙特卡罗法确定溶洞位置、确定溶洞大小生成溶洞，通过体岩溶率检验模型是否符合要求，最后将岩石物理力学参数赋予介质区域，将已知溶洞内包含物的物理力学参数赋予溶洞区域，得到待分析的洞穴型岩溶岩体有限元数值模型。本发明提供的洞穴型岩溶岩体模型可以灵活地描述实际工程中溶洞的分布特征，完整、准确地描述有溶洞分布的岩溶岩体，具有较强适应性，对于评价岩溶区地基、考察岩溶区渗漏规律、研究岩溶区隧道建设等方面具有十分重要的意义。
Description
技术领域
本发明涉及地下工程和岩土工程技术领域，具体地说，涉及一种洞穴型岩溶岩体有限元数值模型的构建方法。
背景技术
可溶岩除了分布有软弱结构面以外，还分布有影响岩体稳定的岩溶结构体，从而形成了岩溶岩体。岩溶结构体具有不规则的特性，因此可以认为岩溶岩体是一种极不均质的岩体。这些岩溶结构体的分布，成为了影响岩体工程力学性质和稳定性的主要因素。
岩溶分为洞穴型、裂隙型和管道型三个形态，其中洞穴型是指发育的干溶洞或填充型溶洞，裂隙型是指由各种构造裂隙经溶蚀形成的岩溶裂隙，管道型是指岩溶裂隙经进一步溶蚀扩大呈汇流的的管道特征。
岩体内部由于岩溶发育而存在缺陷，这些大小不等，位置随机、非随机的岩溶缺陷的存在，改变了岩体的力学特性，降低了岩体的承载能力。现有技术中，许多研究只是针对单个溶洞，只能用来分析溶洞对周围岩体的影响。然而岩溶的形态大小不同，分布复杂，当前建立的模型并不能准确地描述岩溶发育地区的岩体。
基于上述情况，需要一种能够准确地描述岩溶岩体的有限元数值模型的构建方法。
发明内容
为了解决上述问题，同时由于岩溶分为洞穴型、裂隙型、管道型三个形态，本发明提供了一种能够准确地描述洞穴型岩溶岩体的有限元数值模型的构建方法。
本发明通过以下技术方案实现：
洞穴型岩溶岩体数值模型构建方法，其特征在于包括以下步骤：
确定构建岩体数值模型的工程范围，并对工程范围进行勘察、钻探和物探；
(1)获得体岩溶率或线岩溶率，通过工程钻探，获得溶洞、岩石、溶隙长度或钻孔穿过可溶岩的长度参数，并计算得到体岩溶率或线岩溶率；
(2)获得溶洞位置范围，通过工程勘察和物探，确定溶洞的分布范围或可溶岩的分布范围；
(3)获得溶洞直径分布规律，通过工程勘察获得溶洞直径数据或岩溶发育情况，确定该地区溶洞直径的分布规律，获得其分布特征的概率密度函数；
(4)建立数值模型，基于蒙特卡罗法确定溶洞位置、溶洞直径并生成溶洞，通过体岩溶率检验模型是否符合要求；
(5)将岩石物理力学参数赋予介质区域，将已知溶洞内包含物的物理力学参数赋予溶洞区域，得到所述工程范围的岩溶岩体有限元数值模型。
所述获得体岩溶率是通过工程地质钻探获得相关溶洞体积和岩石总体积，并通过体岩溶率＝溶洞体积/岩石总体积×100％计算得到；所述线岩溶率是通过工程地质钻探获得相关溶洞长度或溶隙长度、钻孔穿过可溶岩的长度并通过线岩溶率＝钻孔所遇溶洞长度或溶隙长度/钻孔穿过可溶岩的长度×100％计算得到。
所述溶洞分布范围或可溶岩的分布范围以溶洞位置概率模型数学表达式表示为：
其中，ci(i＝1,2,3)是溶洞中心的空间坐标，ai、bi(i＝1,2,3)表示岩体空间边界，R表示[0,1]上的随机数。
所述溶洞直径分布规律可以通过统计分析现场溶洞直径数据获得，或者通过岩溶发育程度选择分布规律。
通过岩溶发育程度选择溶洞直径分布规律的概率密度函数是：
其中，α是形状参数，β是尺度参数，d是溶洞直径，D是最大溶洞直径，直径单位为米。
当岩溶弱发育时，α＝1.5,β＝1.3,0＜d≤4；岩溶中等发育时，α＝5.3,β＝1.2,0＜d≤10；岩溶强烈发育时，α＝3.7,β＝2.4,0＜d≤15。
所述建立数值模型是采用有限元软件内置语言进行编程，生成符合体岩溶率、溶洞位置和溶洞大小要求的溶洞区域，包括以下步骤：
第一步：参数的生成，确定体岩溶率或线岩溶率，设置模型生成范围区间为X米×Y米×Z米，将岩体区域Ω沿边划分为l×h×k个小区域Δ，对小区域进行编号并将其中心坐标存入数组中其编号对应的位置；将溶洞中心空间坐标ci(i＝1,2,3)的确定转化为小区域编号的确定，在溶洞位置概率模型得到的条件下由随机数发生器确定小区域编号，以溶洞直径分布规律的概率密度函数得到的分布规律为准则生成溶洞直径的随机数列，并处理为溶洞半径r的随机数列；
第二步：确定溶洞生成的区间，设置一个过滤数组用来记录小区域Δ是否可用，过滤数组初始元素均为1，在产生溶洞之后将溶洞占用的小区域对应的元素值修改为0；在第n个溶洞产生之后将元素值为1的小区域收集作为新的抽样区间Φn，以保证有效溶洞的生成效率；
将溶洞半径随机数列导入ANSYS中，如有限元软件ANSYS15.0中，在确定第n个溶洞的位置时已知其半径值rn，遍历抽样区间Φ，将边界rn范围内的小区域对应的过滤数组元素值修改为0；为了避免溶洞交错和后续网格划分的健康，设置网格尺寸elemsize，遍历生成区间Φ将距离已经生成的第i个溶洞ri+rn+elemsize(0＜i＜n)范围内的小区域对应的过滤数组元素值修改为0；最后在借助过滤数组形成生成区间Φn-1并确定第n个溶洞的位置；
第三步：生成溶洞并计算体岩溶率，在确定的位置生成溶洞，计算此时已生成的溶洞的体积之和与模型体积的比值，如未达到体岩溶率的要求，则开始下一个溶洞的生成，如果已达到体岩溶率要求，则程序结束并返回体岩溶率值；或，在生成区间为0时，程序结束并返回体岩溶率值。
本发明的有益效果：相较于现有的溶洞有限元数值模型建立方法，本发明提供的洞穴型岩溶岩体模型可以灵活地描述实际工程中溶洞的分布特征，更加完整、准确地描述有溶洞分布的岩溶岩体，具有较强适应性，对于评价岩溶区地基、考察岩溶区渗漏规律、研究岩溶区隧道建设等方面具有十分重要的意义。
附图说明
图1是本发明的流程示意图。
图2是步骤S104的具体流程示意图。
图3是岩溶发育程度与溶洞直径分布规律的对应关系示意图，图中A表示岩溶弱发育的溶洞直径分布规律，B表示岩溶中等发育的溶洞直径分布规律，C表示岩溶强烈发育的溶洞直径分布规律，纵坐标表示概率密度，横坐标表示溶洞直径，单位：米。
图4是步骤S105的程序框图。
图5是选择域缩减示意图。
图6是本发明的一个实施例的数值模型图。
具体实施方式
下面通过实施例对本发明进行具体的描述，本实施例只用于对发明进行进一步的说明，但不能理解为对本发明保护范围的限制，本领域的技术人员根据上述发明的内容作出的一些非本质的改进和调整也属于本发明保护的范围。
本发明洞穴型岩溶岩体数值模型构建方法包括：
确定构建岩体数值模型的工程范围，并对工程范围进行勘察、钻探和物探；
(1)获得体岩溶率或线岩溶率，通过工程钻探，获得溶洞、岩石、溶隙长度或钻孔穿过可溶岩的长度参数，并计算得到体岩溶率或线岩溶率；
(2)获得溶洞位置范围，通过工程勘察和物探，确定溶洞的分布范围或可溶岩的分布范围；
(3)获得溶洞直径分布规律，通过工程勘察获得溶洞大小数据或岩溶发育情况，确定该地区溶洞直径的分布规律，获得其分布特征的概率密度函数；
(4)建立数值模型，基于蒙特卡罗法确定溶洞位置、溶洞直径并生成溶洞，通过体岩溶率检验模型是否符合要求；
(5)将岩石物理力学参数赋予介质区域，将已知溶洞内包含物的物理力学参数赋予溶洞区域，得到所述工程范围的岩溶岩体有限元数值模型。
结合图1至图5，对本发明方法说明如下：
图1是本发明的流程示意图。
步骤S101，对实际工程进行勘察、钻探和物探等，并对得到的资料进行统计分析。
步骤S102，确定该岩溶区域的体岩溶率。
体岩溶率是评价该区域岩溶发育程度的核心指标。勘查方法一般使用钻探法。体岩溶率＝溶洞体积/岩石总体积×100％，在体岩溶不容易得到时，可以用线岩溶率代替体岩溶率。线岩溶率＝钻孔所遇溶洞、溶隙长度/钻孔穿过可溶岩的长度×100％。
步骤S103，通过勘察和探测资料，确定溶洞在地层中的分布范围。
溶洞发育在可溶性岩层，溶洞的分布范围基本上就是可溶性岩层的分布范围。
由于地质复杂性，溶洞基本呈随机均匀分布。确定溶洞位置的ci(i＝1,2,3)是溶洞中心的空间坐标，溶洞位置概率模型数学表达式为：
其中，ci(i＝1,2,3)是溶洞中心的空间坐标，ai、bi(i＝1,2,3)表示岩体空间边界，R表示[0,1]上的随机数。
步骤S104，确定溶洞直径的统计规律，具体细节见图2。由于是洞穴型岩溶，溶洞可近似成球体。
图2是步骤S104的具体流程示意图。通过对现场溶洞直径数据统计分析，得到分布规律。在没有现场溶洞直径的统计数据时采取步骤S204，由该地区岩溶发育程度选定岩溶大小分布规律。
岩溶发育程度评定表见表1。岩溶发育程度与溶洞直径分布规律对应关系见图3，溶洞直径概率密度函数表达式为：α是形状参数，β是尺度参数，d是溶洞直径，D是最大溶洞直径，直径单位为米。岩溶发育程度不同，表达式参数不同：岩溶弱发育时，α＝1.5,β＝1.3,0＜d≤4；岩溶中等发育时，α＝5.3,β＝1.2,0＜d≤10；岩溶强烈发育时，α＝3.7,β＝2.4,0＜d≤15。
表1 岩溶发育程度评定表
步骤S105，使用有限元软件内置语言进行编程，生成符合体岩溶率、溶洞位置和溶洞大小要求的溶洞区域，具体编程细节见图4，对图4解读如下。
第一步：参数的生成。体岩溶率已在步骤S102中确定。设置模型生成区间为X米×Y米×Z米，将岩体区域Ω沿边划分为l×h×k个小区域Δ，对小区域进行编号并将其中心坐标存入数组中其编号对应的位置。这样将溶洞中心空间坐标ci(i＝1,2,3)的确定转化为小区域编号的确定，在步骤S103得到的条件下由随机数发生器确定小区域编号；以步骤S104得到的分布规律为准则生成溶洞直径的随机数列，并处理为溶洞半径r的随机数列。
第二步：确定溶洞生成的区间。由于溶洞的产生会占用空间，使溶洞生成的区间发生变化(见图5)，为提高溶洞产生效率，所以设置一个过滤数组用来记录小区域Δ是否可用。过滤数组初始元素均为1，在产生溶洞之后将溶洞占用的小区域对应的元素值修改为0。这样在第n个溶洞产生之后将元素值为1的小区域收集作为新的抽样区间Φn，以保证有效溶洞的生成效率。
将溶洞半径随机数列导入ANSYS中，因此在确定第n个溶洞的位置时已知其半径值rn，为防止溶洞超出模型边界，遍历抽样区间Φ，将边界rn范围内的小区域对应的过滤数组元素值修改为0。为了避免溶洞交错和后续网格划分的健康，设置网格尺寸elemsize，遍历生成区间Φ将距离已经生成的第i个溶洞ri+rn+elemsize范围内的小区域对应的过滤数组元素值修改为0。最后在借助过滤数组形成生成区间Φn-1并确定第n个溶洞的位置。
第三步：生成溶洞并计算体岩溶率。在确定的位置生成溶洞，计算此时已生成的溶洞的体积之和与模型体积的比值。如未达到体岩溶率的要求，则开始下一个溶洞的生成；如果已达到体岩溶率要求，则程序结束并返回体岩溶率值。(在生成区间为0时，程序同样结束)
步骤S106，在步骤S105得到的模型中，将岩石物理力学参数赋予介质区域，将已知溶洞内包含物的物理力学参数赋予溶洞区域，得到洞穴型岩溶岩体有限元数值模型。
图6是本发明一个实施例的有限元数值模型图。
在本实施例中，首先进行步骤S101得到相关资料。进行步骤S102，根据钻孔所遇溶洞、溶隙长度和钻孔穿过可溶岩的长度计算取该研究地层区域体岩溶率为3％。
进行步骤S103，根据工程勘察确定溶洞大多分布在地表以下5～26m深度范围内。由于地质复杂性，溶洞基本呈随机均匀分布。确定溶洞位置的ci(i＝1,2,3)是溶洞中心的空间坐标，溶洞位置概率模型数学表达式为：
其中，ai、bi(i＝1,2,3)表示岩体空间边界；R表示[0,1]上的随机数。
进行步骤S104，根据勘察资料及其岩溶率的取值，确定该研究地层区域属于岩溶弱发育，因此选定其概率密度函数数学表达式为：
其中，α＝1.5，β＝1.3，D＝4。
需要说明的是，本发明中溶洞大小的分布规律服从伽马分布，但本发明不限于此，根据现场数据分析得到的其他分布也是合理的，例如正态分布。
进行步骤S105。首先，由步骤S103中公式(3)的岩体空间边界确定岩溶岩体的生成范围，将岩体区域沿边划分为若干个小区域Δ，对小区域进行编号并将其中心坐标存入数组中其编号对应的位置。通过随机数发生器选取小区域编号，并读取其中心坐标作为即将生成的溶洞中心空间坐标ci(i＝1,2,3)。然后，以步骤S104得到的分布规律为准则生成溶洞直径的随机数列，并处理为溶洞半径r的随机数列。最后，将数据导入ANSYS15.0，运行程序进行溶洞的生成。
进行步骤S106，将岩石物理力学参数赋予介质区域，将溶洞内包含物的物理力学参数赋予溶洞区域，最终得到该研究地层区域的岩溶岩体数值模型。
如图6所示，图中长方体为可溶岩层，图中各球体是随机生成的溶洞区域。从最终岩溶岩体数值模型中可以看出，溶洞的参数由概率密度函数生成与实际情况表现出很好的统计相似性，随机大小的溶洞随机分布在岩体空间中充分体现了岩体中溶洞的空间分布特征，能完整、准确的描述岩溶岩体。在该岩溶岩体数值模型的基础上，可以对该研究岩溶地层的地基稳定性、岩体质量和地下水渗流规律等进行研究。
Claims (6)
1.一种洞穴型岩溶岩体数值模型构建方法，其特征在于包括以下步骤：
确定构建岩体数值模型的工程范围，并对工程范围进行勘察、钻探和物探；
(1)获得体岩溶率或线岩溶率，通过工程钻探，获得溶洞、岩石、溶隙长度或钻孔穿过可溶岩的长度参数，并计算得到体岩溶率或线岩溶率；
(2)获得溶洞位置范围，通过工程勘察和物探，确定溶洞的分布范围或可溶岩的分布范围；
(3)获得溶洞直径分布规律，通过工程勘察获得溶洞大小数据或岩溶发育情况，确定该地区溶洞直径的分布规律，获得其分布特征的概率密度函数；
(4)建立数值模型，基于蒙特卡罗法确定溶洞位置、溶洞直径并生成溶洞，通过体岩溶率检验模型是否符合要求；
(5)将岩石物理力学参数赋予介质区域，将已知溶洞内包含物的物理力学参数赋予溶洞区域，得到所述工程范围的岩溶岩体有限元数值模型。
2.根据权利要求1所述洞穴型岩溶岩体数值模型构建方法，其特征在于：所述获得体岩溶率是通过工程地质钻探获得相关溶洞体积和岩石总体积，并通过体岩溶率＝溶洞体积/岩石总体积×100％计算得到；所述线岩溶率是通过工程地质钻探获得相关溶洞长度或溶隙长度、钻孔穿过可溶岩的长度并通过线岩溶率＝钻孔所遇溶洞长度或溶隙长度/钻孔穿过可溶岩的长度×100％计算得到。
3.根据权利要求2所述洞穴型岩溶岩体数值模型构建方法，其特征在于：所述溶洞分布范围或可溶岩的分布范围以溶洞位置概率模型数学表达式表示为：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>f</mi>
<mo>(</mo>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>R</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
其中，ci(i＝1,2,3)是溶洞中心的空间坐标，ai、bi(i＝1,2,3)表示岩体空间边界，R表示[0,1]上的随机数。
4.根据权利要求3所述的洞穴型岩溶岩体数值模型构建方法，其特征在于：所述溶洞直径分布规律通过统计分析现场溶洞直径数据获得，或者通过岩溶发育程度选择分布规律；所述通过岩溶发育程度选择溶洞直径分布规律的概率密度函数是：
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msup>
<mi>&beta;</mi>
<mi>&alpha;</mi>
</msup>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msup>
<mi>d</mi>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>&beta;</mi>
<mi>d</mi>
</mrow>
</msup>
<mo>,</mo>
<mn>0</mn>
<mo>&le;</mo>
<mi>d</mi>
<mo>&le;</mo>
<mi>D</mi>
</mrow>
其中，α是形状参数，β是尺度参数，d是溶洞直径，D是最大溶洞直径，直径单位为米。
5.根据权利要求4所述的洞穴型岩溶岩体数值模型构建方法，其特征在于：当岩溶弱发育时，α＝1.5,β＝1.3,0＜d≤4；岩溶中等发育时，α＝5.3,β＝1.2,0＜d≤10；岩溶强烈发育时，α＝3.7,β＝2.4,0＜d≤15。
6.根据权利要求1至5任一项所述的洞穴型岩溶岩体数值模型构建方法，其特征在于：所述建立数值模型是采用有限元软件内置语言进行编程，生成符合体岩溶率、溶洞位置和溶洞大小要求的溶洞区域，包括以下步骤：
第一步：参数的生成，确定体岩溶率或线岩溶率，设置模型生成范围区间为X米×Y米×Z米，将岩体区域Ω沿边划分为l×h×k个小区域Δ，对小区域进行编号并将其中心坐标存入数组中其编号对应的位置；将溶洞中心空间坐标ci(i＝1,2,3)的确定转化为小区域编号的确定，在溶洞位置概率模型得到的条件下由随机数发生器确定小区域编号，以溶洞直径分布规律的概率密度函数得到的分布规律为准则生成溶洞直径的随机数列，并处理为溶洞半径r的随机数列；
第二步：确定溶洞生成的区间，设置一个过滤数组用来记录小区域Δ是否可用，过滤数组初始元素均为1，在产生溶洞之后将溶洞占用的小区域对应的元素值修改为0；在第n个溶洞产生之后将元素值为1的小区域收集作为新的抽样区间Φn，以保证有效溶洞的生成效率；
将溶洞半径随机数列导入ANSYS中，在确定第n个溶洞的位置时已知其半径值rn，遍历抽样区间Φ，将边界rn范围内的小区域对应的过滤数组元素值修改为0；为了避免溶洞交错和后续网格划分的健康，设置网格尺寸elemsize，遍历生成区间Φ将距离已经生成的第i个溶洞ri+rn+elemsize(0＜i＜n)范围内的小区域对应的过滤数组元素值修改为0；最后在借助过滤数组形成生成区间Φn-1并确定第n个溶洞的位置；
第三步：生成溶洞并计算体岩溶率，在确定的位置生成溶洞，计算此时已生成的溶洞的体积之和与模型体积的比值，如未达到体岩溶率的要求，则开始下一个溶洞的生成，如果已达到体岩溶率要求，则程序结束并返回体岩溶率值；或，在生成区间为0时，程序结束并返回体岩溶率值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710378892.2A CN107180135B (zh) | 2017-05-25 | 2017-05-25 | 一种洞穴型岩溶岩体数值模型构建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710378892.2A CN107180135B (zh) | 2017-05-25 | 2017-05-25 | 一种洞穴型岩溶岩体数值模型构建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107180135A true CN107180135A (zh) | 2017-09-19 |
| CN107180135B CN107180135B (zh) | 2020-08-07 |
Family
ID=59832703
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710378892.2A Expired - Fee Related CN107180135B (zh) | 2017-05-25 | 2017-05-25 | 一种洞穴型岩溶岩体数值模型构建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107180135B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107843712A (zh) * | 2017-10-27 | 2018-03-27 | 中国石油大学(华东) | 一种岩溶水垂向掘洞机制的模拟装置 |
| CN108090313A (zh) * | 2018-02-05 | 2018-05-29 | 东北大学 | 一种复杂岩石裂隙模型建模识别方法 |
| CN112626950A (zh) * | 2020-11-24 | 2021-04-09 | 山东省交通规划设计院有限公司 | 一种结合土拱效应对岩溶路基的处置方法 |
| CN112883464A (zh) * | 2021-01-26 | 2021-06-01 | 南通大学 | 一种软岩隧道开挖引起围岩大变形的不确定性预测方法 |
| CN115014887A (zh) * | 2022-04-14 | 2022-09-06 | 贵阳市城市轨道交通集团有限公司 | 一种浅表地层分布式岩溶模型的制备方法和装置 |
| CN117973082A (zh) * | 2024-03-28 | 2024-05-03 | 山东省地矿工程勘察院(山东省地质矿产勘查开发局八〇一水文地质工程地质大队) | 一种岩溶发育区模拟预测方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103577642A (zh) * | 2013-11-06 | 2014-02-12 | 中国建筑股份有限公司 | 一种隐伏溶洞与新建山岭隧道间安全距离确定方法 |
| US8682629B2 (en) * | 2010-05-25 | 2014-03-25 | Schlumberger Technology Corporation | Multi-phasic dynamic karst reservoir numerical simulator |
-
2017
- 2017-05-25 CN CN201710378892.2A patent/CN107180135B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8682629B2 (en) * | 2010-05-25 | 2014-03-25 | Schlumberger Technology Corporation | Multi-phasic dynamic karst reservoir numerical simulator |
| CN103577642A (zh) * | 2013-11-06 | 2014-02-12 | 中国建筑股份有限公司 | 一种隐伏溶洞与新建山岭隧道间安全距离确定方法 |
Non-Patent Citations (2)
| Title |
|---|
| LAURENT EISENLOHR .ETC: "Numerical versus statistical modelling of natural response of a", 《JOURNL OF HYDROLOGY》 * |
| 赵明阶 等: "岩溶尺寸对隧道围岩稳定性影响的模型试验研究", 《岩石力学与工程学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107843712A (zh) * | 2017-10-27 | 2018-03-27 | 中国石油大学(华东) | 一种岩溶水垂向掘洞机制的模拟装置 |
| CN108090313A (zh) * | 2018-02-05 | 2018-05-29 | 东北大学 | 一种复杂岩石裂隙模型建模识别方法 |
| CN108090313B (zh) * | 2018-02-05 | 2021-03-19 | 东北大学 | 一种复杂岩石裂隙模型建模识别方法 |
| CN112626950A (zh) * | 2020-11-24 | 2021-04-09 | 山东省交通规划设计院有限公司 | 一种结合土拱效应对岩溶路基的处置方法 |
| CN112883464A (zh) * | 2021-01-26 | 2021-06-01 | 南通大学 | 一种软岩隧道开挖引起围岩大变形的不确定性预测方法 |
| CN115014887A (zh) * | 2022-04-14 | 2022-09-06 | 贵阳市城市轨道交通集团有限公司 | 一种浅表地层分布式岩溶模型的制备方法和装置 |
| CN117973082A (zh) * | 2024-03-28 | 2024-05-03 | 山东省地矿工程勘察院(山东省地质矿产勘查开发局八〇一水文地质工程地质大队) | 一种岩溶发育区模拟预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107180135B (zh) | 2020-08-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107180135A (zh) | 一种洞穴型岩溶岩体数值模型构建方法 | |
| CN101894189B (zh) | 煤层底板突水评价的新方法 | |
| Jeannin et al. | Speleological investigations | |
| CN104727813B (zh) | 一种源储共生型地层的孔隙度测量方法 | |
| CN104749617A (zh) | 一种多尺度裂缝储层正演模型建立方法 | |
| CN104865614A (zh) | 基于可变骨架参数的复杂储层流体识别方法 | |
| CN103075150A (zh) | 一种原孔位多次应力解除法的地应力测试方法 | |
| CN103065051A (zh) | 一种对岩体自动进行分级分区的方法 | |
| CN104297810A (zh) | 一种地层元素测井中净非弹性散射伽马能谱的获取方法 | |
| CN105134185A (zh) | 储层流体性质识别方法 | |
| CN109710891A (zh) | 一种基于流动单元分类计算砂岩型铀矿砂体渗透率的方法 | |
| Chesnaux et al. | Regional fracture network permeability using outcrop scale measurements | |
| CN105931125A (zh) | 一种致密油分段多簇体积压裂水平井产量预测方法 | |
| Lv et al. | Elaborate simulation and predication of the tunnel drainage effect on karst groundwater field and discharge based on Visual MODFLOW | |
| CN105182421B (zh) | 一种定量评价地层脆性特征的方法 | |
| CN201137491Y (zh) | 脉冲中子双谱流体饱和度测井仪 | |
| CN107515430A (zh) | 一种地震法探测盐湖卤水的方法 | |
| CN112302623B (zh) | 一种基于测井资料的花岗岩潜山风化壳结构划分方法 | |
| CN103698494B (zh) | 一种确定岩性圈闭含油气饱和度的方法及装置 | |
| Kurtuluş | Istanbul geotechnical downhole arrays | |
| CN113960288B (zh) | 一种源-汇系统定量评价方法 | |
| Lou et al. | A new method for detecting karst and groundwater by 3D seismic wave: Case study of the karst tunnel in Zhangjihuai Railway, China | |
| CN114135277B (zh) | 基于地化特征随钻感知的隧道超前地质预报方法及系统 | |
| CN112069444B (zh) | 一种利用测井资料计算储集层试井渗透率方法及计算机 | |
| CN111397663B (zh) | 依据地下水温度来判别岩溶发育程度的方法 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200807 Termination date: 20210525 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |