CN107037502A - 一种基于flac 3d数值模拟的地应力场反演方法 - Google Patents
一种基于flac 3d数值模拟的地应力场反演方法 Download PDFInfo
- Publication number
- CN107037502A CN107037502A CN201710245023.2A CN201710245023A CN107037502A CN 107037502 A CN107037502 A CN 107037502A CN 201710245023 A CN201710245023 A CN 201710245023A CN 107037502 A CN107037502 A CN 107037502A
- Authority
- CN
- China
- Prior art keywords
- stress
- model
- flac
- measuring point
- map
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000004088 simulation Methods 0.000 title claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000001514 detection method Methods 0.000 claims abstract description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 238000005096 rolling process Methods 0.000 claims description 2
- 239000003245 coal Substances 0.000 description 12
- 239000011435 rock Substances 0.000 description 7
- 238000012360 testing method Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种基于FLAC 3D数值模拟的地应力场反演方法,包括:依据实际工程地质条件,选取测点进行地应力测试;根据区域构造发育情况,建立合适的数值模型;根据测点地应力测试结果,定义模型边界条件;对数值模型进行运算处理,输出模型应力云图;根据地应力测点位置,结合模型应力云图,对区域地应力进行分析。本发明依据测点地应力实测结果,利用FLAC 3D数值模拟软件建立地应力场反演模型,之后对模型加载,满足其生成的应力场和测点应力状态基本一致,从而反演整个工作面地应力场的分布特征,更能直观合理的反映实际地应力场,且操作简单,易于实际应用。
Description
技术领域
本发明涉及煤炭地质工程领域的地应力场反演方法,具体涉及一种基于FLAC 3D数值模拟的地应力场反演方法。
背景技术
地应力是指天然岩土体未受任何扰动的应力,目前对地应力的实测方法主要有数理统计法、灰色理论法、人工智能法等,这些方法大多基于实测地应力值的基础上进行的计算方法,不能直观反映实际的地应力场分布特征。本发明基于这个原因,依据测点实际地应力值,利用FLAC 3D数值模拟建立合理的地应力场反演模型,更能直观反映实际地应力场分布特征。
发明内容
1.本发明的目的
鉴于此,本发明针对煤炭地质工程领域地应力场反演方法的不足,提供一种更能直观反映实际地应力场分布特征的地应力场反演方法。
2.本发明的技术方案
为实现上述目的,本发明一种基于FLAC 3D数值模拟的地应力场反演方法,该方法包括:步骤A,依据实际工程地质条件,选取测点进行地应力测试;步骤B,根据区域构造发育情况,建立合适的数值模型;步骤C,根据测点地应力测试结果,定义模型边界条件;步骤D,对数值模型进行运算处理,输出模型应力云图;步骤E,根据地应力测点位置,结合模型应力云图,对区域地应力场进行分析。
本发明依据测点地应力实测结果,利用FLAC 3D数值模拟软件建立地应力场反演模型,之后对模型加载,满足其生成的应力场和测点应力状态基本一致,反演整个工作面地应力场的分布特征,更能直观合理的反映实际地应力场,且操作简单,易于实际应用。
附图说明
此处用来说明的附图是为了对本发明的进一步解释和说明,为本申请的一部分,但并不能限定本发明。
图1为本发明基于FLAC 3D数值模拟的地应力场反演方法流程图。
图2为本发明实例某矿11113工作面地应力场反演地质模型。
图3A、3B及图3C分别为本发明实例模型初始应力反演云图、剖面Z方向应力反演云图及剖面X方向应力反演云图。
具体实施方案
下面结合附图及本发明所应用的实例,对本发明进行进一步说明。
图1为本发明基于FLAC 3D数值模拟的地应力场反演方法流程图。如图1所示,该方法包括:
步骤A,依据实际工程地质条件,选取测点进行地应力测试;
步骤B,根据区域构造发育情况,建立合适的数值模型;
步骤C,根据测点地应力测试结果,定义模型边界条件;
步骤D,对数值模型进行运算处理,输出模型应力云图;
步骤E,根据地应力测点位置,结合模型应力云图,对区域地应力场进行分析。
以下结合具体实例,对上述每一步骤进行解释说明。
步骤A,依据实际工程地质条件,选取测点进行地应力测试。
利用声发射(AE)法测试地应力,对定向试件进行单轴压缩,同时接收其声发射信号。根据凯塞效应,当作用力达到某一临界值时,声发射活动会突然增多,此临界应力值即为岩芯试件先前所受的应力,本发明实例中的试验采取的岩样是从钻孔中采取具有明显层面的岩芯样,因此可以知道采样处岩层层面走向,具体操作步骤是:
①现场取定向岩芯;
②按一定工艺要求在室内加工声发射试样;
③在选定参数下对试样加载并记录AE率,分别绘出AE事件率—时间曲线和应力—时间曲线;
④以K点和响应曲线斜率陡增点的选值决定各试样的正应力,而后由各方向的正应力计算主应力大小和方向。
⑤数据整理
按照以上所示方法整理试验数据来判断取样方向上的正应力值,测试结果为12.9MPa,然后根据下文的公式来计算各测点的主应力值及方向。
因为取样时,各岩样的夹角为45°,可采用以下公式计算最大、最小主应力的大小和方向:
式中:分别是平行、垂直和45°方向的正应力;σ1为平面最大主应力;σ3为平面最小主应力。并且规定应力以压为正;φ角以由主应力σ1逆时针转到方向为正。各测点的应力值见表1。
表1各测点的应力大小和方向
步骤B,根据区域构造发育情况,建立合适的数值模型。本发明实例中,根据测点地应力实测结果及岩层物理参数,运用FLAC 3D数值模拟建立地应力场反演地质模型。具体:
计算模型是某矿11113工作面的地质剖面图和地质平面图经过适当简化确定而成。11113工作面为东翼A组煤采区首采面,位于3煤层第一阶段,回采上限位于3煤层80m防水煤柱和风氧化带露头之间。该面设计走向长度870.0m,倾向长度108m,据资料分析,本区段煤岩层为一单斜构造,地层倾向北东,地层倾角变化较大在∠6~26°平均16°。在下顺槽拨门位置发育DF9-1正断层,该断层落差10m,断层倾向N,倾角为70°。拟模型介质类型共分9种,其中包括砂岩组、泥岩组、煤层、断层带等。
根据某矿区11113工作面地质资料,结合断层带附近工程地质条件特征,选取-650m水平为背景,断层带厚度据揭露显示取值3m,以3煤层为特征煤层,嵌入地质模型,煤厚5m。
设计模型几何尺寸为长300m、宽200m、高120m,共分为7200000个单元格,煤岩层近似水平,初始地应力12.46Mpa,侧压系数取0.25,主要以构造应力为主。本发明实例所建的地应力场反演地质模型,如图2所示。
步骤C,根据测点地应力测试结果,定义模型边界条件。本发明实例中,在吸取前人的经验条件下,初步将边界条件定义为混合边界条件,具体如下:
(1)位移边界
①将模型底部定义为固定支承,限制节点在X和Z两个方向的运动(X、Z均为水平方向);②将模型左侧定义为滚支承,限制节点在X方向的运动。
(2)应力边界
①在模型顶部施加向下的垂向应力;②在模型X、Y方向施加边界应力。
以3煤层底板等高线图为底图,而待拟合的点位于等高线-490与-400水平上,将地应力拟合点位置投放其中,模拟仅以-490水平地应力测试结果为例。
步骤D,对数值模型进行运算处理,输出模型应力云图。本发明实例中,通过不断调整边界荷载的大小值,并计算到最终结果能满足在两个测点的模拟值与实测值最为接近。在此应力作用下模型计算至基本平衡,应力不再作大的调整时结束运算。此时通过输入FLAC命令语句输出模型初始体力分布图,作切片得出最大、小主应力图、X、Z方向应力云图,来综合分析模型受力特征。输出结果见图3A-3C。
步骤E,根据地应力测点位置,结合模型应力云图,对区域地应力场进行分析。
由图3A、图3B、图3C可知,通过模型模拟、输出图件分析可知:Z方向与X方向的应力分布较为均一,在断层带附近有所变化,但整体上呈现等幅变化趋势;最大与最小主应力变化从图件显示为正常地应力分布;主应力迹线图及其细化图上体现出了模型受力作用下内部主应力的变化大小与方向性,图中所示“L”型显示了所在单元的主应力的大小和方向,其中较长的线段表示最大主应力,较小线段的表示最小主应力,并且其值与长度成正比,直观反应了整体或局部的应力状态。
本发明依据测点地应力实测结果,利用FLAC 3D数值模拟软件建立地应力场反演模型,之后对模型加载,满足其生成的应力场和测点应力状态基本一致,反演整个工作面地应力场的分布特征,更能直观合理的反映实际地应力场,且操作简单,易于实际应用。
以上所述具体实例,对本发明的目的,过程和有益效果进行详细说明,并不用于限定本发明的限定范围,凡在本发明的精神原则之内,所做的任何修改、等同替换等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种研究一种基于FLAC 3D数值模拟的地应力场反演方法,其特征在于,该方法包括:
步骤A,依据实际工程地质条件,选取测点进行地应力测试;
步骤B,根据区域构造发育情况,建立合适的数值模型;
步骤C,根据测点地应力测试结果,定义模型边界条件;
步骤D,对数值模型进行运算处理,输出模型应力云图;
步骤E,根据地应力测点位置,结合模型应力云图,对区域地应力进行分析。
2.如权利要求1所述的基于FLAC 3D数值模拟的地应力场反演方法,其特征在于,步骤B中,根据区域构造发育情况,建立合适的数值模型。
3.如权利要求1所述的基于FLAC 3D数值模拟的地应力场反演方法,其特征在于,步骤C中,根据测点地应力测试结果,定义模型边界条件。具体如下:
(1)位移边界
①将模型底部定义为固定支承,限制节点在X和Z两个方向的运动(X、Z均为水平方向);②将模型左侧定义为滚支承,限制节点在X方向的运动。
(2)应力边界
①在模型顶部施加向下的垂向应力;②在模型X、Y方向施加边界应力。
4.如权利要求1所述的基于FLAC 3D数值模拟的地应力场反演方法,其特征在于,步骤D中,对数值模型进行运算处理,输出模型应力云图。通过不断调整边界荷载的大小值,并计算到最终结果能满足在两个测点的模拟值与实测值最为接近。在此应力作用下模型计算至基本平衡,应力不再作大的调整时结束运算。此时通过输入FLAC命令语句输出模型初始体力分布图,作切片得出最大、小主应力图、X、Z方向应力云图,来综合分析模型受力特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710245023.2A CN107037502A (zh) | 2017-04-14 | 2017-04-14 | 一种基于flac 3d数值模拟的地应力场反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710245023.2A CN107037502A (zh) | 2017-04-14 | 2017-04-14 | 一种基于flac 3d数值模拟的地应力场反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107037502A true CN107037502A (zh) | 2017-08-11 |
Family
ID=59535131
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710245023.2A Pending CN107037502A (zh) | 2017-04-14 | 2017-04-14 | 一种基于flac 3d数值模拟的地应力场反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107037502A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108645552A (zh) * | 2018-05-28 | 2018-10-12 | 中国地质大学(北京) | 一种基于双应力***的地应力测量方法 |
CN110162882A (zh) * | 2019-05-23 | 2019-08-23 | 山东东山新驿煤矿有限公司 | 一种煤矿巷道穿越断层区围岩精准控制方法 |
CN111881612A (zh) * | 2020-08-05 | 2020-11-03 | 武汉市政工程设计研究院有限责任公司 | 一种正应力和剪应力不同权重二维应力场反演方法及*** |
-
2017
- 2017-04-14 CN CN201710245023.2A patent/CN107037502A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108645552A (zh) * | 2018-05-28 | 2018-10-12 | 中国地质大学(北京) | 一种基于双应力***的地应力测量方法 |
CN110162882A (zh) * | 2019-05-23 | 2019-08-23 | 山东东山新驿煤矿有限公司 | 一种煤矿巷道穿越断层区围岩精准控制方法 |
CN111881612A (zh) * | 2020-08-05 | 2020-11-03 | 武汉市政工程设计研究院有限责任公司 | 一种正应力和剪应力不同权重二维应力场反演方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106485015B (zh) | 一种矿井断层影响范围的确定方法 | |
CN104866653B (zh) | 一种获取地下三维密度结构的方法 | |
CN113901681B (zh) | 一种全寿命周期页岩气储层双甜点三维可压性评估方法 | |
CN106484958B (zh) | 一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的*** | |
CN107066771B (zh) | 一种平推式滑坡运动距离计算方法及应用 | |
CN108693572B (zh) | 一种基于三维建模的地应力场反演方法 | |
Ning et al. | Study of the in situ stress field in a deep valley and its influence on rock slope stability in Southwest China | |
CN104850897A (zh) | 一种基于地震信息的煤与瓦斯突出的预测方法 | |
CN106291755B (zh) | 一种隆升区低级序断层发育规律定量预测方法 | |
CN106326528A (zh) | 一种露天矿端帮地下开采诱发地表裂缝分布规律预测方法 | |
CN110162882A (zh) | 一种煤矿巷道穿越断层区围岩精准控制方法 | |
CN107037502A (zh) | 一种基于flac 3d数值模拟的地应力场反演方法 | |
Ceccato et al. | Constraints upon fault zone properties by combined structural analysis of virtual outcrop models and discrete fracture network modelling | |
CN107366534A (zh) | 粗化渗透率的确定方法和装置 | |
CN116703244B (zh) | 采矿沉陷区治理效果与综合评价方法 | |
CN107728205B (zh) | 一种地层压力预测方法 | |
Liu et al. | Quantitative multiparameter prediction of fractured tight sandstone reservoirs: a case study of the Yanchang Formation of the Ordos Basin, Central China | |
Han et al. | Recognizing fracture distribution within the coalbed methane reservoir and its implication for hydraulic fracturing: A method combining field observation, well logging, and micro-seismic detection | |
Marchegiani et al. | Scaling properties of the dimensional and spatial characteristics of fault and fracture systems in the Majella Mountain, central Italy | |
CN108665146A (zh) | 大型露天矿山边坡稳定性评价精度确定方法 | |
CN110705168A (zh) | 构造应力场的模拟方法 | |
CN109064560A (zh) | 一种地铁隧道施工过程中地质条件模型化方法 | |
Fang et al. | Prediction method and distribution characteristics of in situ stress based on borehole deformation—a case study of coal measure stratum in shizhuang block, qinshui basin | |
CN110133748A (zh) | 一种集成碱交代型铀矿深部成矿信息的方法 | |
CN106154326A (zh) | 一种纵弯褶皱裂缝密度评价的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170811 |