CN113050169B - 一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法 - Google Patents
一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法 Download PDFInfo
- Publication number
- CN113050169B CN113050169B CN202110291115.0A CN202110291115A CN113050169B CN 113050169 B CN113050169 B CN 113050169B CN 202110291115 A CN202110291115 A CN 202110291115A CN 113050169 B CN113050169 B CN 113050169B
- Authority
- CN
- China
- Prior art keywords
- rock mass
- hole
- wave velocity
- longitudinal wave
- anisotropy
- 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.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims abstract description 97
- 238000004458 analytical method Methods 0.000 title claims abstract description 19
- 238000012614 Monte-Carlo sampling Methods 0.000 title claims abstract description 14
- 238000012360 testing method Methods 0.000 claims abstract description 84
- 238000000034 method Methods 0.000 claims abstract description 35
- 238000005259 measurement Methods 0.000 claims abstract description 13
- 238000005070 sampling Methods 0.000 claims abstract description 8
- 230000008859 change Effects 0.000 claims abstract description 7
- 238000010998 test method Methods 0.000 claims abstract description 5
- 239000000523 sample Substances 0.000 claims description 37
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 230000014509 gene expression Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 238000011065 in-situ storage Methods 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241000197727 Euscorpius alpha Species 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000006467 substitution reaction Methods 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/449—Statistical methods not provided for in G01N29/4409, e.g. averaging, smoothing and interpolation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
-
- 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/303—Analysis for determining velocity profiles or travel times
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Probability & Statistics with Applications (AREA)
- Signal Processing (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法,包括:在岩体的垂直层理面打设三个测试孔,测试孔呈等边三角形分布;采用单孔测试方法、跨孔组合方法,获得各测试孔的轴向纵波波速、径向纵波波速;利用蒙特卡罗随机抽样方法,计算轴向纵波波速、径向纵波波速的标准差和平均值,并随机生成N个轴向纵波波速、N个径向纵波波速的变量样本;根据岩体各向异性系数概念,确定N个岩体各向异性系数;统计N个岩体各向异性系数大于某一确定各向异性系数的超越概率;改变某一确定各向异性系数,绘制超越概率随某一确定各向异性系数的变化曲线,确定岩体各向异性系数。本发明的方法比传统方法操作更加简便,测量结果更加真实可靠。
Description
技术领域
本发明涉及岩体工程技术领域,更具体的涉及一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法。
背景技术
高质量的地质勘察是工程安全顺利建设的重要保证。而在漫长的地质历史过程中,岩体经过了构造变动、风化作用以及卸荷作用等各种内外力地质作用的破坏和改造,形成了复杂的物理力学特性。实际地层中岩体的弹性模量、泊松比、抗拉强度以及抗压强度在空间上分布往往是不均匀的,正确认识岩体的不均匀各向异性对工程的安全高效建设具有重要意义。
已有的岩体各向异性测量方法一般采用室内试验或者现场测量装置获得。采用室内试验的方式,需要钻取并制备大量岩心试样,耗费时间长,花费大,不能充分反映原位工程岩体赋存的地应力环境。通过现场装置测试的方式能够在原位状态下测量岩体的各向异性,但测试装置制作成本大,操作复杂,并不适用于地质条件复杂的工程。然而,岩体弹性波波速探测法是一种常用的无损检测岩体质量的方法,通过声波发射***将电信号转换为声波信号,声波信号在岩体内部传播,并由声波接受***将声波信号转换为电信号,随后利用数据处理***对电信号进行处理得到可用的声波数据。根据显示的声波数据可以得到超声波在岩体内的传播时间以及传播速度,进而反映岩体的松弛特征以及完整性。
工程岩体各个方向的物理力学性质不同,不确定因素较多,具有一定的随机性,而现有的岩体各向异性的测量方法,往往是给出一个或者多个确定值,这仅能表示测点位置岩体的各向异性,并不能代表整个区域岩体的各向异性情况,具有一定的局限性。
发明内容
本发明实施例提供了一种基于蒙特卡罗分析,利用单孔、跨孔声波组合检测方式来快速确定岩体各向异性的方法,旨在解决现有岩体各向异性确定方法中成本高、操作复杂以及检测结果不准确的问题。
本发明实施例提供一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法,包括:
在岩体的垂直层理面打设三个测试孔,测试孔呈等边三角形分布;其中,所述三个测试孔包括:测试孔A、测试孔B、测试孔C;
采用单孔测试方法,获得各测试孔的轴向纵波波速;及采用跨孔组合方法,获得各测试孔的径向纵波波速;
利用蒙特卡罗随机抽样方法,计算轴向纵波波速的标准差和平均值,并随机生成N个轴向纵波波速的变量样本;及利用蒙特卡罗随机抽样方法,计算径向纵波波速的标准差和平均值,并随机生成N个径向纵波波速的变量样本;
根据岩体各向异性系数概念,结合N个轴向纵波波速的变量样本、N个径向纵波波速的变量样本,确定N个岩体各向异性系数;
统计N个岩体各向异性系数大于某一确定各向异性系数的超越概率;改变某一确定各向异性系数,绘制超越概率随某一确定各向异性系数的变化曲线,确定岩体各向异性系数。
进一步地,
所述测试孔的半径大于超声波仪器测量探头的最大半径1.5~2.0倍;
所述等边三角形的边长为1.0~1.5m。
进一步地,采用单孔测试方法,获得测试孔A的轴向纵波波速,具体包括:
向测试孔A中注水,直至测试孔A孔口有水流出;
将超声波发射探头置于测试孔A孔底、第一超声波接收探头置于测试孔A孔底上方位置处,操作声波仪进行测试;
对声波仪显示器中的数据进行采集并保存,且对每一对测定读数三次,最大读数之差不大于3%,记录此时波速为测试孔A的轴向纵波波速CAZ。
进一步地,采用跨孔组合方法,获得测试孔A的径向纵波波速,具体包括:
将超声波发射探头置于测试孔A孔底、第二超声波接收探头置于测试孔B孔底、第三超声波接收探头置于测试孔C孔底,操作声波仪进行测试;
对声波仪显示器中的数据进行采集并保存,且对每一对测定读数三次,最大读数之差不大于3%,记录此时波速CA→B;
将超声波发射探头置于测试孔B孔底、第二超声波接收探头置于测试孔A孔底、第三超声波接收探头置于测试孔C孔底,操作声波仪进行测试;
对声波仪显示器中的数据进行采集并保存,且对每一对测定读数三次,最大读数之差不大于3%,记录此时波速CB→A;
对波速CA→B和波速CB→A求平均值,确定径向纵波波速:
其中,β为岩体各向同性面与水平面的夹角。
进一步地,所述根据岩体各向异性系数概念,结合N个轴向纵波波速的变量样本、N个径向纵波波速的变量样本,确定N个岩体各向异性系数;具体包括:
岩体各向异性系数ηd由轴线方向上动弹性模量Edz比垂直轴线方向上的动弹性模量Edx来确定,如下式:
岩体的动弹性模量Ed,如下式:
式中,ρ为岩体密度;ν为泊松比;Vp为岩体的纵波波速;
根据上述两式,确定岩体的各向异性系数表达式:
将N个轴向纵波波速的变量样本、N个径向纵波波速的变量样本代入上式,获得N个岩体各向异性系数ηd;Ciz为随机生成的轴向纵波波速;Cix为随机生成的径向纵波波速。
进一步地,所述统计N个岩体各向异性系数大于某一确定各向异性系数的超越概率,具体包括:
统计大于某一确定各向异性系数η的样本数n,则大于某一确定岩体各向异性系数η的超越概率为:
增加样本总数N的值,直到超越概率P不再随着N的增加而变化,此时超越概率P即为大于某一确定岩体各向异性系数η的超越概率。
本发明实施例提供一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法,与现有技术相比,其有益效果如下:
本发明涉及一种原位岩体各向异性的测量方法,更具体的是基于蒙特卡罗分析方法,利用单孔、跨孔声波组合检测的方式快速确定岩体各向异性系数。蒙特卡罗分析法是一种采用随机抽样统计来估算结果的计算方法,通过蒙特卡罗分析法对数据进行处理,能够计算大于某一确定各向异性值的概率,得到工程岩体各向异性值的取值范围,更贴近工程现场实际情况,增加计算结果的可靠度。本发明的方法适用于工程建设中岩体各向异性的合理、高效的测量,比传统方法操作更加简便,测量结果更加真实可靠。
附图说明
图1为本发明实施例提供的测试孔分布平面图;
图2为本发明实施例提供的测试孔分布剖面图;
图3为本发明实施例提供的基于蒙特卡罗抽样的岩体各向异性系数概率分析方法流程图;
图4为本发明实施例提供的超越概率P随各向异性值η的变化规律。
附图标记:
1-第一超声波信号接收探头;2-超声波信号发射探头;3-第二超声波信号接收探头;4-第三超声波信号接收探头;5-测试孔A;6-测试孔B;7-测试孔C;8-层理岩体。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法,该方法包括:
步骤1:在岩体的垂直层理面打设三个测试孔,测试孔呈等边三角形分布;其中,所述三个测试孔包括:测试孔A、测试孔B、测试孔C。
步骤2:采用单孔测试方法,获得各测试孔的轴向纵波波速;及采用跨孔组合方法,获得各测试孔的径向纵波波速。
步骤3:利用蒙特卡罗随机抽样方法,计算轴向纵波波速的标准差和平均值,并随机生成N个轴向纵波波速的变量样本;及利用蒙特卡罗随机抽样方法,计算径向纵波波速的标准差和平均值,并随机生成N个径向纵波波速的变量样本。
步骤4:根据岩体各向异性系数概念,结合N个轴向纵波波速的变量样本、N个径向纵波波速的变量样本,确定N个岩体各向异性系数。
步骤5:统计N个岩体各向异性系数大于某一确定各向异性系数的超越概率;改变某一确定各向异性系数,绘制超越概率随某一确定各向异性系数的变化曲线,确定岩体各向异性系数。
本发明的主要原理如下:
(1)岩体的各向异性系数可由轴线方向上动弹性模量比径向的动弹性模量来确定,如式(1)
而岩体的动弹性模量Ed可由式(2)表示:
式中,ρ为岩体密度;ν为泊松比;Cp为岩体的纵波波速。
(2)假设η为某一确定岩体各向异性参数,其中Ciz,Cix为轴线方向以及径向的纵波波速,由于岩体本身物理力学性质的复杂性,这里将Ciz,Cix视为服从正态分布的随机变量。令α=η-ηd,则当α<0时,表明由随机变量Ciz,Cix求得的岩体各向异性参数ηd要大于确定岩体各向异性参数η。假设随机变量的样本总数为N,其中有n个样本的计算结果小于0,即α<0,则为保证概率P有足够的精确度,增加样本总数N的数量,直到概率P不再随N的增加而敏感,此时P即为超过确定各向异性值η的概率。同理,改变η的值,即可得到超过不同各向异性值的概率。
对于上述步骤1~5的具体说明如下:
(1)在待测工程岩体的临空壁面垂直层理面钻设三个测试孔,分别标注为A测孔、B测孔、C测孔。三个测孔的孔口中心点连线构成一等边三角形,边长一般为1.0~1.5m,测试孔半径一般要大于超声波仪器的测量探头的最大半径1.5~2.0倍,参见图1和图2。
(2)利用单孔测试的方法获得各个测孔轴线方向的纵波波速,分别记为CAZ、CBZ、CCZ,并计算各个测孔轴线方向的纵波波速的平均值Mz和标准差σz。利用空间穿透测试的方法获得测孔A到测孔B、测孔B到测孔C、测孔C到测孔A间岩体径向的纵波波速,分别记为CAX、CBX、CCX,并计算各水平方向纵波波速的平均值Mx和标准差σx。
(3)假设各随机变量(Ciz,Cix)服从正态分布,根据各个测孔轴线方向的纵波波速的平均值Mz和标准差σz,随机生成N个轴线方向的纵波波速的变量样本,同理,根据各径向纵波波速的平均值Mx和标准差σx,随机生成N个径向的纵波波速的变量样本。这里按照正态分布的方式生成随机变量样本,并不局限于正态分布,也可以是其他分布函数,如韦伯分布。
(4)将N个变量样本代入到式(2)中,即可得到N个各向异性值,再将计算得到的各向异性值带入式(3)中,统计α<0的样本数量n,则可得超过某一确定各向异性值的概率增加样本总数N的数量,直到概率P不再随N的增加而变化,此时P即为超过确定各向异性值η的概率,记为超越概率。
式中:η为某一确定各向异性值;ηd为计算得到的各向异性值;Ciz为随机生成的轴向纵波波速;Cix为随机生成的径向纵波波速。
(5)改变η的值,得到超过不同各向异性值的超越概率P,以不同各向异性值为横轴,以超越概率P为纵轴,绘制超越概率随各向异性值变化的曲线,参见图4,根据变化曲线即可得到小于某一概率的岩体各向异性值。
实施例分析:
下面结合具体实施例,对本发明进行进一步说明,具体包括如下步骤:
(1)通过量角器或罗盘确定层理岩体各向同性面与水平面的夹角β,在层理岩体8的临空壁面上标记出A测孔5、B测孔6、C测孔7三个测孔的中心点位置,测孔间垂直间距为1.0~1.5m,三个测孔中心点连线构成一等边三角形,这样可以保证在进行跨孔测试时,超声波信号能够被信号接收***清晰识别。根据标记中心点位置,利用钻机垂直层理面钻设A测孔5、B测孔6、C测孔7,钻设深度设置为5m。
(2)向A测孔5注水,直至孔口有水流出,若孔壁有漏水现象,则须在测试过程中不断注水,保证后续测试中超声波接收探头3和超声波接收探头4始终被水浸没。将超声波发射探头2缓缓置于测试孔孔底,超声波接收探头1置于距孔底上方一定位置处,操作声波仪进行测试,当声波仪显示器中出现稳定、清晰的波形后,对数据进行采集并保存,每一对测定应读数三次,最大读数之差不宜大于3%,记录此时纵波波速为CAZ。同理,分别测量B测孔6、C测孔7中的轴向纵波波速,并记为CBZ、CCZ。
(3)将超声波发射探头2置于A测孔5孔底,超声波接收探头3和超声波接收探头4分别置于B测孔6以及C测孔7孔底,操作声波仪进行测试,当声波仪显示器中出现稳定、清晰的波形后,对数据进行采集并保存,每一对测定应读数三次,最大读数之差不宜大于3%,记录此时纵波波速为CA→B。同理,改变发射探头2的位置于B测孔、C测孔,测量CB→A、CB→C、CC→A、CC→B。考虑到现场地质条件复杂,为提高测量的准确度,这里取CA→B、CB→A的平均值,则径向纵波波速可表示为取CB→C、CC→B的平均值,则径向纵波波速可表示为取CC→A、CA→C的平均值,则径向纵波波速可表示为
(4)由现场测试得到的各径向纵波波速以及轴向纵波波速,分别计算其平均值M和标准差σ,具体见表1。根据计算得到的标准差σ与平均值M,按照正态分布随机生成20组样本数据,将样本数据代入式(2)计算ηd。
表1现场实测纵波波速
(5)假设η=1.4,将ηd、η代入式(3),统计α<0的样本数量。由表2可知α<0的样本数量为11,样本总数为20,则P=11/20=0.55。由于样本总量较小,所以此时的P并不为一个恒定值,增加样本总数,并计算P值,直到P值趋近于稳定,这里计算了样本总数为20、100、500、1000、5000、10000时的P值,具体见表3。由表3可知,当随机样本总数足够大时,P值趋近于0.376,故各向异性值大于1.4的概率为0.376。改变η值,可以得到超过不同各向异性值η的概率,以η为横坐标,超过η值的超越概率P为纵坐标,绘制超越概率P随各向异性值η值的变化曲线,如附图3所示。由附图3可知,随着η值的增加,超越概率P逐渐减小,当η为1.28时,超过η值的概率为1,当η大于1.52时,超过η值的概率为0.01,综上,可以认为该工程岩体的各向异性值在1.28~1.52之间。
表2随机生成的样本数据
表3不同样本量对应的超越概率P
以上公开的仅为本发明的几个具体实施例,本领域的技术人员可以对本发明实施例进行各种改动和变型而不脱离本发明的精神和范围,但是,本发明实施例并非局限于此,任何本领域的技术人员能思之的变化都应落入本发明的保护范围内。
Claims (6)
1.一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法,其特征在于,包括:
在岩体的垂直层理面打设三个测试孔,测试孔呈等边三角形分布;其中,所述三个测试孔包括:测试孔A、测试孔B、测试孔C;
采用单孔测试方法,获得各测试孔的轴向纵波波速;及采用跨孔组合方法,获得各测试孔的径向纵波波速;
利用蒙特卡罗随机抽样方法,计算轴向纵波波速的标准差和平均值,并随机生成N个轴向纵波波速的变量样本;及利用蒙特卡罗随机抽样方法,计算径向纵波波速的标准差和平均值,并随机生成N个径向纵波波速的变量样本;
根据岩体各向异性系数概念,结合N个轴向纵波波速的变量样本、N个径向纵波波速的变量样本,确定N个岩体各向异性系数;
统计N个岩体各向异性系数大于某一确定各向异性系数的超越概率;改变某一确定各向异性系数,绘制超越概率随某一确定各向异性系数的变化曲线,确定岩体各向异性系数。
2.如权利要求1所述的基于蒙特卡罗抽样的岩体各向异性系数概率分析方法,其特征在于,
所述测试孔的半径大于超声波仪器测量探头的最大半径1.5~2.0倍;
所述等边三角形的边长为1.0~1.5m。
3.如权利要求1所述的基于蒙特卡罗抽样的岩体各向异性系数概率分析方法,其特征在于,采用单孔测试方法,获得测试孔A的轴向纵波波速,具体包括:
向测试孔A中注水,直至测试孔A孔口有水流出;
将超声波发射探头置于测试孔A孔底、第一超声波接收探头置于测试孔A孔底上方位置处,操作声波仪进行测试;
对声波仪显示器中的数据进行采集并保存,且对每一对测定读数三次,最大读数之差不大于3%,记录此时波速为测试孔A的轴向纵波波速CAZ。
4.如权利要求1所述的基于蒙特卡罗抽样的岩体各向异性系数概率分析方法,其特征在于,采用跨孔组合方法,获得测试孔A的径向纵波波速,具体包括:
将超声波发射探头置于测试孔A孔底、第二超声波接收探头置于测试孔B孔底、第三超声波接收探头置于测试孔C孔底,操作声波仪进行测试;
对声波仪显示器中的数据进行采集并保存,且对每一对测定读数三次,最大读数之差不大于3%,记录此时波速CA→B;
将超声波发射探头置于测试孔B孔底、第二超声波接收探头置于测试孔A孔底、第三超声波接收探头置于测试孔C孔底,操作声波仪进行测试;
对声波仪显示器中的数据进行采集并保存,且对每一对测定读数三次,最大读数之差不大于3%,记录此时波速CB→A;
对波速CA→B和波速CB→A求平均值,确定径向纵波波速:
其中,β为岩体各向同性面与水平面的夹角。
5.如权利要求1所述的基于蒙特卡罗抽样的岩体各向异性系数概率分析方法,其特征在于,所述根据岩体各向异性系数概念,结合N个轴向纵波波速的变量样本、N个径向纵波波速的变量样本,确定N个岩体各向异性系数;具体包括:
岩体各向异性系数ηd由轴线方向上动弹性模量Edz比垂直轴线方向上的动弹性模量Edx来确定,如下式:
岩体的动弹性模量Ed,如下式:
式中,ρ为岩体密度;ν为泊松比;Vp为岩体的纵波波速;
根据上述两式,确定岩体的各向异性系数表达式:
将N个轴向纵波波速的变量样本、N个径向纵波波速的变量样本代入上式,获得N个岩体各向异性系数ηd;Ciz为随机生成的轴向纵波波速;Cix为随机生成的径向纵波波速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110291115.0A CN113050169B (zh) | 2021-03-18 | 2021-03-18 | 一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110291115.0A CN113050169B (zh) | 2021-03-18 | 2021-03-18 | 一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113050169A CN113050169A (zh) | 2021-06-29 |
CN113050169B true CN113050169B (zh) | 2021-11-05 |
Family
ID=76513286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110291115.0A Active CN113050169B (zh) | 2021-03-18 | 2021-03-18 | 一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113050169B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104374827A (zh) * | 2013-08-12 | 2015-02-25 | 中国科学院武汉岩土力学研究所 | 横观各向同性岩体原位动弹性模量的各向异性系数的测量方法 |
CN106872579A (zh) * | 2017-02-13 | 2017-06-20 | 长江勘测规划设计研究有限责任公司 | 正态分布拟合岩体波速划分岩体质量分级的方法 |
CN109655903A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | 页岩层横波速度预测方法及*** |
CN111764369A (zh) * | 2020-06-18 | 2020-10-13 | 河海大学 | 一种融合试验与监测数据挖掘的高坝岩体卸荷松弛评价方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090070042A1 (en) * | 2007-09-11 | 2009-03-12 | Richard Birchwood | Joint inversion of borehole acoustic radial profiles for in situ stresses as well as third-order nonlinear dynamic moduli, linear dynamic elastic moduli, and static elastic moduli in an isotropically stressed reference state |
US7649804B2 (en) * | 2007-10-25 | 2010-01-19 | Schlumberger Technology Corporation | In-situ determination of yield stress state of earth formations |
WO2014084751A1 (en) * | 2012-11-30 | 2014-06-05 | Schlumberger Holdings Limited | A method for processing acoustic waveforms |
-
2021
- 2021-03-18 CN CN202110291115.0A patent/CN113050169B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104374827A (zh) * | 2013-08-12 | 2015-02-25 | 中国科学院武汉岩土力学研究所 | 横观各向同性岩体原位动弹性模量的各向异性系数的测量方法 |
CN106872579A (zh) * | 2017-02-13 | 2017-06-20 | 长江勘测规划设计研究有限责任公司 | 正态分布拟合岩体波速划分岩体质量分级的方法 |
CN109655903A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | 页岩层横波速度预测方法及*** |
CN111764369A (zh) * | 2020-06-18 | 2020-10-13 | 河海大学 | 一种融合试验与监测数据挖掘的高坝岩体卸荷松弛评价方法 |
Non-Patent Citations (1)
Title |
---|
不同损伤程度下岩石力学参数变化的声波测试;严鹏;《岩土力学》;20151231;第36卷(第12期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113050169A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111208198A (zh) | 一种岩体实时波速测定及质量评价的方法 | |
Herron et al. | Ultrasonic velocities and crystalline anisotropy in the ice core from Dye 3, Greenland | |
CN110082429B (zh) | 一种结合机器学习的隧道衬砌无损检测的辅助判定方法 | |
WO2006016520A1 (ja) | 音響トモグラフィを用いた樹木内部の非破壊検査方法及び装置 | |
CN111487315A (zh) | 一种隧道衬砌厚度和脱空的声频无损检测方法 | |
CN109612611B (zh) | 一种超声检测结构永存应力的方法 | |
JP2004150946A (ja) | ボール打撃によるコンクリート剛性の非破壊測定装置および方法 | |
CN104374827A (zh) | 横观各向同性岩体原位动弹性模量的各向异性系数的测量方法 | |
CN106525978A (zh) | 一种利用剪切模量变化计算软土结构性扰动度的方法 | |
CN105606463A (zh) | 一种基于中智函数的岩体结构面抗剪强度综合评价方法 | |
CN114880755B (zh) | 铁路隧道的围岩分级方法、装置、设备及可读存储介质 | |
CN109459497B (zh) | 一种含平行裂缝横向各向同性岩石的声波速度计算方法 | |
CN113050169B (zh) | 一种基于蒙特卡罗抽样的岩体各向异性系数概率分析方法 | |
Song et al. | Prediction and analysis of geomechanical properties of jimusaer shale using a machine learning approach | |
CN117270053B (zh) | 一种基于偶极声波测井的地层横波慢度和衰减计算方法 | |
CN113158315A (zh) | 一种基于静力触探数据的岩土体参数三维非平稳条件随机场建模方法 | |
CN109709610B (zh) | 一种岩石裂缝探测方法及*** | |
Govindaraj | A new approach for estimation of properties of metamorphic rocks | |
CN110231401B (zh) | 一种测点间信息关联方法及土木结构无损检测方法 | |
Yu et al. | Acoustic emission enabled particle size estimation via low stress-varied axial interface shearing | |
CN113622905B (zh) | 一种基于多因素综合分析的页岩储层脆性评价方法 | |
Niu et al. | Laboratory small-strain stiffness measurement using distributed acoustic sensing | |
CN112611805B (zh) | 一种基于衰减系数的评价围岩松动圈范围的方法 | |
CN110806444B (zh) | 基于浅地层剖面仪与svm的海底底质识别与分类方法 | |
CN109100801B (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 |