CN108693572B - 一种基于三维建模的地应力场反演方法 - Google Patents
一种基于三维建模的地应力场反演方法 Download PDFInfo
- Publication number
- CN108693572B CN108693572B CN201810200389.2A CN201810200389A CN108693572B CN 108693572 B CN108693572 B CN 108693572B CN 201810200389 A CN201810200389 A CN 201810200389A CN 108693572 B CN108693572 B CN 108693572B
- Authority
- CN
- China
- Prior art keywords
- stress
- stratum
- model
- geological 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000005070 sampling Methods 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 12
- 238000011160 research Methods 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims description 19
- 239000011435 rock Substances 0.000 claims description 18
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 238000005553 drilling Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 239000011229 interlayer Substances 0.000 claims description 3
- 239000010410 layer Substances 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 8
- 230000006378 damage Effects 0.000 description 5
- 238000007789 sealing Methods 0.000 description 4
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 description 2
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000013201 Stress fracture Diseases 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- 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/14—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 using acoustic emission techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种基于三维建模的地应力场反演方法,该方法只需要在工程现场进行取样工作即可,有效减轻了现场工作量,可以实现多个测点取样,实现地应力多测点精细测量;可建立工程大尺度含单一或复杂构造的三维数值模型,通过并行计算对反演研究区域的地应力分布特征进行反演;可得到直观、可视化的整体地应力分布情况,有效提取并查看模型不同水平、不同倾斜面、不同构造处的垂直应力、最大水平主应力、最小水平主应力值及方向等信息。
Description
技术领域
本发明涉及地质勘测领域,特别是涉及一种基于三维建模的地应力场反演方法。
背景技术
地应力是存在于地壳中的未受工程扰动的天然应力,含构造的地应力场往往更加复杂多变。长期以来的实践表明,地应力是引起隧道、煤矿巷道围岩、支护结构变形和破坏以及矿井动力现象的根本作用力。当井巷开挖后,初始应力场和变形场受到开挖扰动的影响而重新分配,造成围岩损伤,严重时导致失稳、垮塌和破坏,而当井巷处于地质构造的影响范围内时,构造对应力场的影响对巷道失稳破坏则尤为明显。目前地应力现场测试方法有水压致裂法、应力解除法、钻孔崩落法等,实验室测试有声发射法、岩石微裂隙闭合变形临界应力法等,但现场岩体应力测试是一项费时费力、耗费财力的复杂工作,一定程度上限制了现场岩体应力测试的推广应用,同时实测值只能代表测点附近一定范围内的应力情况,有限的测点并不能反映区域整体的应力场。在此情况下寻找能直观有效的反应地下工程整体地应力特征的方法非常有必要,通过三维数值模拟软件进行含构造地应力反演是其中一种有效的途径,但目前现场实测地应力普遍测点较少,同时利用三维数值模拟软件对地应力进行反演时,所建立的三维模型为了平衡计算速度与计算精度,往往尺寸较小,难以达到工程大尺度,实际计算效果与精度会受到局限,同时,方便、快捷的建立含构造的三维数值模型以及地应力场进行有效的反演也是一个难题。
发明内容
本发明主为解决现有问题的不足之处而提供一种基于三维建模的地应力场反演方法。
为解决上述技术问题,本发明采用的一个技术方案是:提供一种基于三维建模的地应力场反演方法,该方法的步骤包括:选定地质模型的范围,在所述地质模型区域边界以及内部的各岩层布置多个测点,对每一测点分别进行采样,并标记采样样品的原始方位,将每一采样样品置于空间直角坐标系中,沿坐标系内指定方向在采样样品上钻取;应用伺服压力试验机和声发射设备,对每一钻取的试样进行物理力学性质的测试与声发射采集,计算各个测点对应的垂直应力、最大水平主应力和最小水平主应力的大小和方向;其中,物理力学性质的测试至少得到试样的抗压强度、弹性模量和泊松比参数;根据地质模型的地质报告,建立地质模型工程尺度的断层、褶曲以及二者组合构造的三维地质模型,对构造的三维地址模型划分网格并确定边界条件、载荷条件;对三维地质模型进行反演计算,得到整个地质模型区域的应力场分布,将地质模型的多个测点与反演模拟结果对应点的计算应力大小与方向进行比较,利用边界荷载反分析法验证地应力反演效果,得到应力反演结果。
其中,在沿坐标系内指定方向在采样样品上钻取的步骤中,钻取试样的方向为沿X、Y、Z、X45°Y、Y45°Z、Z45°X空间六个方向,每一方向钻取的试件数量至少为3个。
其中,采样样本数量至少为8个,体积不小于300mm3,试件体积大于等于Φ25mm×50mm,对每一精磨加工以符合岩石力学试验标准。
其中,根据地质模型的地质报告,建立地质模型工程尺度的断层、褶曲以及二者组合构造的三维地质模型的步骤包括:确定整个地质模型研究范围以及范围内各地层岩石物理力学参数,构造类型、规模、数量以及分布情况;划分地质模型的地层模型截面;其中,对于断层模型,依据各地层厚度建立地层截面;对于褶曲模型,根据褶曲的幅高、区内长度、两翼倾角等参数,选取褶曲各地层边界位置的关键点,用曲线依次连接各地层选定的边界点,形成褶曲地层边界;对于断层和褶曲的组合构造,同褶曲地层模型的建立;添加断层元素,建立地质模型截面;对于断层模型,依据断层的倾角、断距、张开度、落差参数,在地层截面的基础上添加断层,并对地层进行相应的错动;对于褶曲模型,无需再添加断层;对于断层和褶曲组合构造,在褶曲地层截面的基础上同样添加对应的断层,并进行褶曲地层的错动,完成组合构造模型的建立;将建好的各模型截面的图形导入MIDAS软件,生成地层、褶曲、断层相应的封闭面,对封闭面进行实体生成,完成整个对应的含构造的三维地质模型的建立。
其中,在对构造的三维地质模型划分网格并确定边界条件、载荷条件的步骤中,包括步骤:在三维地质模型的构造区域以及层间接触面精细划分网格,其他位置一般划分网格;将三维地质模型底面垂直位移固定,上表面自由,四个侧面法向位移固定作为初始边界条件,以三维地质模型上覆岩层的重力作为初始载荷条件,对应模拟自重应力过程;根据在三维地质模型边界测点的实测应力值大小作为水平应力的载荷条件,三维地质模型底面垂直位移固定,相邻两个侧面法向位移固定,另两个侧面施加对应的地应力实测值,对应模拟构造应力过程。
其中,利用并行计算机对MIDAS模拟软件中的三维地质模型进行反演计算,得到整个地质模型的应力场分布,将每一测点与反演模拟结果对应点的计算应力大小与方向进行比较,利用边界荷载反分析法验证地应力反演效果,调整边界载荷以使各测点反演值不断逼近实测值,当各测点反演值与实测值二者之间的误差达到10%以内时,反演结果即为最终的结果。
其中,根据计算结果查看不同位置垂直应力、最大水平主应力、最小水平主应力等值线分布图、计算结果曲线图数据,并对整体应力场分布、构造等因素对应力场分布的影响进行分析。
区别于现有技术,本发明的基于三维建模的地应力场反演方法只需要在工程现场进行取样工作即可,有效减轻了现场工作量,可以实现多个测点取样,实现地应力多测点精细测量;可建立工程大尺度含单一或复杂构造的三维数值模型,通过并行计算对反演研究区域的地应力分布特征进行反演;可得到直观、可视化的整体地应力分布情况,有效提取并查看模型不同水平、不同倾斜面、不同构造处的垂直应力、最大水平主应力、最小水平主应力值及方向等信息。
附图说明
图1是本发明提供的一种基于三维建模的地应力场反演方法的流程示意图。
图2是本发明提供的一种基于三维建模的地应力场反演方法的逻辑示意图。
图3是本发明提供的一种基于三维建模的地应力场反演方法的地质模型的工程尺度的断层建模的建模示意图。
图4是本发明提供的一种基于三维建模的地应力场反演方法的地质模型的工程尺度的褶曲建模的建模示意图。
图5是本发明提供的一种基于三维建模的地应力场反演方法的地质模型的工程尺度的复杂构造建模的建模示意图。
图6是本发明提供的一种基于三维建模的地应力场反演方法的三维地质模型整体垂直主应力、最大水平主应力、最小水平主应力等值线图。
图7是本发明提供的一种基于三维建模的地应力场反演方法的断层倾斜方向的垂直主应力、最大水平主应力、最小水平主应力等值线图。
图8是本发明提供的一种基于三维建模的地应力场反演方法的褶曲剖面垂直主应力、最大水平主应力、最小水平主应力等值线图。
具体实施方式
下面结合具体实施方式对本发明的技术方案作进一步更详细的描述。显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例,都应属于本发明保护的范围。
参阅图1和图2,图1是本发明提供的一种基于三维建模的地应力场反演方法的流程示意图;图2是本发明提供的一种基于三维建模的地应力场反演方法的逻辑示意图。该方法的步骤包括:
S110:选定地质模型的范围,在地质模型区域边界以及内部的各岩层布置多个测点,对每一测点分别进行采样,并标记采样样品的原始方位,将每一采样样品置于空间直角坐标系中,沿坐标系内指定方向在采样样品上钻取。
采样样本数量至少为8个,钻取试样的方向为沿X、Y、Z、X45°Y、Y45°Z、Z45°X空间六个方向,每一方向钻取的试件数量至少为3个。具体的,首先确定垂直方向的方位,再在上面寻找一个近水平的面,确定正北方向N,以N作为X方向,建立X、Y、Z空间坐标系,并分别在XOY平面、XOZ平面和YOZ平面上确定与两轴夹角为45°的方向X45°Y、Y45°Z、Z45°X。然后按照X、Y、Z、X45°Y、Y45°Z、Z45°X空间六个方向钻取试件,每个方向至少钻取3个,试件尺寸不小于Φ25mm×50mm,试件精磨加工,标准要求符合岩石力学试验标准。
S120:应用伺服压力试验机和声发射设备,对每一钻取的试样进行物理力学性质的测试与声发射采集,计算各个测点对应的垂直应力、最大水平主应力和最小水平主应力的大小和方向;其中,物理力学性质的测试至少得到试样的抗压强度、弹性模量和泊松比参数。
应用伺服压力试验机和声发射设备,根据声发射能量在时间上的分布特征(第一次陡增点)确定每个试件凯塞效应特征点对应的应力值,对相同方向上的K值取平均值,得到、、、、、六个方向的值,根据相应的弹性力学计算公式,求得主应力量值、、及其方向,并最终确定各测点垂直应力、最大水平主应力和最小水平主应力的大小及方向。
S130:根据地质模型的地质报告,建立地质模型工程尺度的断层、褶曲以及二者组合构造的三维地质模型,对构造的三维地址模型划分网格并确定边界条件、载荷条件。
工程尺度的断层建模、褶曲建模及复杂构造建模示意图如图3-5所示根据地质模型的地质报告,建立地质模型工程尺度的断层、褶曲以及二者组合构造的三维地质模型的步骤包括:
确定整个地质模型研究范围以及范围内各地层岩石物理力学参数,构造类型、规模、数量以及分布情况;
划分地质模型的地层模型截面;其中,对于断层模型,依据各地层厚度建立地层截面;对于褶曲模型,根据褶曲的幅高、区内长度、两翼倾角等参数,选取褶曲各地层边界位置的关键点,用曲线依次连接各地层选定的边界点,形成褶曲地层边界;对于断层和褶曲的组合构造,同褶曲地层模型的建立。
添加断层元素,建立地质模型截面;对于断层模型,依据断层的倾角、断距、张开度、落差参数,在地层截面的基础上添加断层,并对地层进行相应的错动;对于褶曲模型,无需再添加断层;对于断层和褶曲组合构造,在褶曲地层截面的基础上同样添加对应的断层,并进行褶曲地层的错动,完成组合构造模型的建立;
将建好的各模型截面的图形导入MIDAS软件,生成地层、褶曲、断层相应的封闭面,对封闭面进行实体生成,完成整个对应的含构造的三维地质模型的建立。
在对构造的三维地质模型划分网格并确定边界条件、载荷条件的步骤中,包括步骤:
在三维地质模型的构造区域以及层间接触面精细划分网格,其他位置一般划分网格;
将三维地质模型底面垂直位移固定,上表面自由,四个侧面法向位移固定作为初始边界条件,以三维地质模型上覆岩层的重力作为初始载荷条件,对应模拟自重应力过程;
根据在三维地质模型边界测点的实测应力值大小作为水平应力的载荷条件,三维地质模型底面垂直位移固定,相邻两个侧面法向位移固定,另两个侧面施加对应的地应力实测值,对应模拟构造应力过程。
S140:对三维地质模型进行反演计算,得到整个地质模型区域的应力场分布,将地质模型的多个测点与反演模拟结果对应点的计算应力大小与方向进行比较,利用边界荷载反分析法验证地应力反演效果,得到应力反演结果。
利用并行计算机对MIDAS模拟软件中的三维地质模型进行反演计算,得到整个地质模型的应力场分布,将每一测点与反演模拟结果对应点的计算应力大小与方向进行比较,利用边界荷载反分析法验证地应力反演效果,调整边界载荷以使各测点反演值不断逼近实测值,当各测点反演值与实测值二者之间的误差达到10%以内时,反演结果即为最终的结果。
根据计算结果查看不同位置垂直应力、最大水平主应力、最小水平主应力等值线分布图、计算结果曲线图数据,并对整体应力场分布、构造等因素对应力场分布的影响进行分析。具体的,三维地质模型整体、断层倾斜方向及褶曲剖面垂直主应力、最大水平主应力、最小水平主应力等值线图如图6-8所示。
区别于现有技术,本发明的基于三维建模的地应力场反演方法只需要在工程现场进行取样工作即可,有效减轻了现场工作量,可以实现多个测点取样,实现地应力多测点精细测量;可建立工程大尺度含单一或复杂构造的三维数值模型,通过并行计算对反演研究区域的地应力分布特征进行反演;可得到直观、可视化的整体地应力分布情况,有效提取并查看模型不同水平、不同倾斜面、不同构造处的垂直应力、最大水平主应力、最小水平主应力值及方向等信息。
以上仅为本发明的实施方式,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (6)
1.一种基于三维建模的地应力场反演方法,其特征在于,包括:
选定地质模型的范围,在所述地质模型区域边界以及内部的各岩层布置多个测点,对每一测点分别进行采样,并标记采样样品的原始方位,将每一采样样品置于空间直角坐标系中,沿坐标系内指定方向在采样样品上钻取;
应用伺服压力试验机和声发射设备,对每一钻取的试样进行物理力学性质的测试与声发射采集,计算各个测点对应的垂直应力、最大水平主应力和最小水平主应力的大小和方向;其中,物理力学性质的测试至少得到试样的抗压强度、弹性模量和泊松比参数;
确定整个地质模型研究范围以及范围内各地层岩石物理力学参数,构造类型、规模、数量以及分布情况;
划分地质模型的地层模型截面;其中,对于断层模型,依据各地层厚度建立地层截面;对于褶曲模型,根据褶曲的幅高、区内长度、两翼倾角参数,选取褶曲各地层边界位置的关键点,用曲线依次连接各地层选定的边界点,形成褶曲地层边界;对于断层和褶曲的组合构造,同褶曲地层模型的建立;
添加断层元素,建立地质模型截面;对于断层模型,依据断层的倾角、断距、张开度、落差参数,在地层截面的基础上添加断层,并对地层进行相应的错动;对于褶曲模型,无需再添加断层;对于断层和褶曲组合构造,在褶曲地层截面的基础上同样添加对应的断层,并进行褶曲地层的错动,完成组合构造模型的建立;
将建好的各模型截面的图形导入MIDAS软件,生成地层、褶曲、断层相应的封闭面,对封闭面进行实体生成,完成整个对应的含构造的三维地质模型的建立,对构造的三维地址模型划分网格并确定边界条件、载荷条件;
对三维地质模型进行反演计算,得到整个地质模型区域的应力场分布,将地质模型的多个测点与反演模拟结果对应点的计算应力大小与方向进行比较,利用边界荷载反分析法验证地应力反演效果,得到应力反演结果。
2.根据权利要求1所述的基于三维建模的地应力场反演方法,其特征在于,在沿坐标系内指定方向在采样样品上钻取的步骤中,钻取试样的方向为沿X、Y、Z、X45°Y、Y45°Z、Z45°X空间六个方向,每一方向钻取的试件数量至少为3个。
3.根据权利要求1所述的基于三维建模的地应力场反演方法,其特征在于,采样样本数量至少为8个,体积不小于300mm3,试件体积大于等于Φ25mm×50mm,对每一精磨加工以符合岩石力学试验标准。
4.根据权利要求1所述的基于三维建模的地应力场反演方法,其特征在于,在对构造的三维地质模型划分网格并确定边界条件、载荷条件的步骤中,包括步骤:
在三维地质模型的构造区域以及层间接触面精细划分网格,其他位置一般划分网格;
将三维地质模型底面垂直位移固定,上表面自由,四个侧面法向位移固定作为初始边界条件,以三维地质模型上覆岩层的重力作为初始载荷条件,对应模拟自重应力过程;
根据在三维地质模型边界测点的实测应力值大小作为水平应力的载荷条件,三维地质模型底面垂直位移固定,相邻两个侧面法向位移固定,另两个侧面施加对应的地应力实测值,对应模拟构造应力过程。
5.根据权利要求1所述的基于三维建模的地应力场反演方法,其特征在于,利用并行计算机对MIDAS模拟软件中的三维地质模型进行反演计算,得到整个地质模型的应力场分布,将每一测点与反演模拟结果对应点的计算应力大小与方向进行比较,利用边界荷载反分析法验证地应力反演效果,调整边界载荷以使各测点反演值不断逼近实测值,当各测点反演值与实测值二者之间的误差达到10%以内时,反演结果即为最终的结果。
6.根据权利要求5所述的基于三维建模的地应力场反演方法,其特征在于,根据计算结果查看不同位置垂直应力、最大水平主应力、最小水平主应力等值线分布图,计算结果曲线图数据,并对整体应力场分布、构造因素对应力场分布的影响进行分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810200389.2A CN108693572B (zh) | 2018-03-12 | 2018-03-12 | 一种基于三维建模的地应力场反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810200389.2A CN108693572B (zh) | 2018-03-12 | 2018-03-12 | 一种基于三维建模的地应力场反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108693572A CN108693572A (zh) | 2018-10-23 |
CN108693572B true CN108693572B (zh) | 2020-07-03 |
Family
ID=63844497
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810200389.2A Active CN108693572B (zh) | 2018-03-12 | 2018-03-12 | 一种基于三维建模的地应力场反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108693572B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110162882B (zh) * | 2019-05-23 | 2023-04-07 | 山东东山新驿煤矿有限公司 | 一种煤矿巷道穿越断层区围岩精准控制方法 |
CN110823704B (zh) * | 2019-10-15 | 2020-11-06 | 山东大学 | 一种tbm搭载的岩渣在线抗压强度测试系统及方法 |
CN111538108A (zh) * | 2020-06-12 | 2020-08-14 | 鞍钢集团矿业有限公司 | 基于原生地应力场的矿岩节理构造模拟方法 |
CN112415582B (zh) * | 2020-10-15 | 2022-07-01 | 山东大学 | 一种结构和波速随机布设的三维速度地质建模方法 |
CN112800649B (zh) * | 2021-01-15 | 2022-08-23 | 中国石油大学(北京) | 一种基于有限元数值模拟软件的三维断层模型建立方法 |
CN113935215B (zh) * | 2021-10-13 | 2022-08-23 | 西南石油大学 | 断层发育油气藏地应力反演方法 |
CN114564775B (zh) * | 2022-02-14 | 2022-10-11 | 中国地质大学(北京) | 适用于采空区上方修建隧道的双尺度耦合数值模拟方法 |
CN114861519B (zh) * | 2022-03-07 | 2023-06-30 | 成都理工大学 | 复杂地质条件下初始地应力场加速优化反演方法 |
CN116341293B (zh) * | 2023-05-30 | 2023-09-01 | 煤炭科学研究总院有限公司 | 三维应力场反演方法及装置 |
CN116341292B (zh) * | 2023-05-30 | 2023-08-08 | 煤炭科学研究总院有限公司 | 应力场的获取方法、装置、设备及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101162177A (zh) * | 2007-11-14 | 2008-04-16 | 南京银茂铅锌矿业有限公司 | 一种测量地应力的方法 |
CN103091727A (zh) * | 2013-01-11 | 2013-05-08 | 上海交通大学 | 基于实测结果及线性规划法的地应力场ddm反演方法 |
CN105866835A (zh) * | 2016-03-28 | 2016-08-17 | 中国石油大学(华东) | 一种基于地应力分布的断层三维封闭性定量评价方法 |
CN106018107A (zh) * | 2016-05-20 | 2016-10-12 | 重庆大学 | 声发射法三维地应力测试方法 |
CN106875471A (zh) * | 2017-01-13 | 2017-06-20 | 山东科技大学 | 煤系含或隔水层三维可视化建模方法 |
-
2018
- 2018-03-12 CN CN201810200389.2A patent/CN108693572B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101162177A (zh) * | 2007-11-14 | 2008-04-16 | 南京银茂铅锌矿业有限公司 | 一种测量地应力的方法 |
CN103091727A (zh) * | 2013-01-11 | 2013-05-08 | 上海交通大学 | 基于实测结果及线性规划法的地应力场ddm反演方法 |
CN105866835A (zh) * | 2016-03-28 | 2016-08-17 | 中国石油大学(华东) | 一种基于地应力分布的断层三维封闭性定量评价方法 |
CN106018107A (zh) * | 2016-05-20 | 2016-10-12 | 重庆大学 | 声发射法三维地应力测试方法 |
CN106875471A (zh) * | 2017-01-13 | 2017-06-20 | 山东科技大学 | 煤系含或隔水层三维可视化建模方法 |
Non-Patent Citations (1)
Title |
---|
深埋大断层构造区三维地应力场反演分析;黄耀光等;《煤矿开采》;20140630;第19卷(第3期);第23-28、8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108693572A (zh) | 2018-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108693572B (zh) | 一种基于三维建模的地应力场反演方法 | |
Yazdani et al. | Displacement-based numerical back analysis for estimation of rock mass parameters in Siah Bisheh powerhouse cavern using continuum and discontinuum approach | |
Guo et al. | Simulation of three-dimensional tectonic stress fields and quantitative prediction of tectonic fracture within the Damintun Depression, Liaohe Basin, northeast China | |
Kulatilake et al. | Evaluation of rock slope stability for Yujian River dam site by kinematic and block theory analyses | |
Poulsen et al. | Convergence of synthetic rock mass modelling and the Hoek–Brown strength criterion | |
Glover et al. | Fluid flow in synthetic rough fractures and application to the Hachimantai geothermal hot dry rock test site | |
Hosseini et al. | Passive seismic velocity tomography on longwall mining panel based on simultaneous iterative reconstructive technique (SIRT) | |
Zhang et al. | Cooperative monitoring and numerical investigation on the stability of the south slope of the Fushun west open-pit mine | |
Nikadat et al. | Numerical modelling of stress analysis around rectangular tunnels with large discontinuities (fault) by a hybridized indirect BEM | |
Tal et al. | An improved numerical manifold method for simulation of sequential excavation in fractured rocks | |
CN108918682B (zh) | 深切河谷斜坡岩体现今天然地应力室内测试分析方法 | |
CN114662360A (zh) | 一种利用多软件联合构建复杂地质三维地应力场反演方法 | |
Tan et al. | Influence of complicated faults on the differentiation and accumulation of in-situ stress in deep rock mass | |
Galindo et al. | An accessible calculation method of the bearing capacity of shallow foundations on anisotropic rock masses | |
Gao et al. | Can we estimate far-field stress using the mean of local stresses? An examination based on numerical simulations | |
CN116542050A (zh) | 一种节理岩体开挖稳定性的多尺度模拟方法 | |
CN114491774B (zh) | 一种深部背斜构造及地层结构三维数值模型构建方法 | |
CN111365072B (zh) | 一种裂隙岩体巨型洞室的导洞勘察结构及方法 | |
Feng et al. | Efficient modeling of the gravity anomaly caused by a sedimentary basin with lateral variable density contrast and its application in basement relief estimation | |
Karagöz et al. | Investigation of wedge sliding potential using different methods in two open pit albite mines excavated in gneiss in western Türkiye | |
Swyer et al. | Using boundary element modeling of fault slip to predict patterns of stress perturbation and related fractures in geothermal reservoirs and explore parameter uncertainty | |
Vitali | Tunnel behavior under complex anisotropic conditions | |
Konietzky | Numerical stress field modelling for underground structures | |
Jing et al. | Three-Dimensional Mechanical Information Model and Its Application to Large Underground Water-Sealed Storage Caverns for Oil or Gas | |
Qinjie et al. | In-situ Stress Measurement and its Engineering Application in Deep Coal Mines: A Case Study in the Xinji Coalfield of China. |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20181023 Assignee: Shanxi Chenshi Electromechanical Technology Co.,Ltd. Assignor: Taiyuan University of Technology Contract record no.: X2022140000012 Denomination of invention: An inversion method of geostress field based on 3d modeling Granted publication date: 20200703 License type: Common License Record date: 20221214 |