CN109113742B - 一种煤储层现今地应力预测方法 - Google Patents
一种煤储层现今地应力预测方法 Download PDFInfo
- Publication number
- CN109113742B CN109113742B CN201810869778.4A CN201810869778A CN109113742B CN 109113742 B CN109113742 B CN 109113742B CN 201810869778 A CN201810869778 A CN 201810869778A CN 109113742 B CN109113742 B CN 109113742B
- Authority
- CN
- China
- Prior art keywords
- coal
- stress
- static
- rock
- dynamic
- 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
- 239000003245 coal Substances 0.000 title claims abstract description 191
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000011435 rock Substances 0.000 claims abstract description 114
- 230000003068 static effect Effects 0.000 claims abstract description 65
- 238000004364 calculation method Methods 0.000 claims abstract description 37
- 238000012360 testing method Methods 0.000 claims abstract description 23
- 238000011161 development Methods 0.000 claims abstract description 12
- 230000000694 effects Effects 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 11
- 238000002474 experimental method Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000005553 drilling Methods 0.000 claims description 3
- 238000013213 extrapolation Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000002706 hydrostatic effect Effects 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000009933 burial Methods 0.000 claims description 2
- 239000007789 gas Substances 0.000 description 19
- 230000035699 permeability Effects 0.000 description 6
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21C—MINING OR QUARRYING
- E21C39/00—Devices for testing in situ the hardness or other properties of minerals, e.g. for giving information as to the selection of suitable mining tools
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种煤储层现今地应力预测方法，即在煤岩参数实验测试和测井曲线计算基础上，建立静态参数与动态参数之间的量化关系，考虑到煤岩低杨氏模量、高泊松比、低强度、易变形等岩石力学特点，反演获取基于Maxwell煤储层现今地应力模型中的未知参数，建立基于测井曲线的煤储层地应力剖面，实现预测煤储层现今地应力的目的。本发明针对性强，可准确预测煤储层现今地应力，为煤层气有效勘探开发奠定基础，减小风险和成本，可广泛应用于煤层气勘探开发技术领域。
Description
技术领域
本发明涉及煤层气勘探开发技术领域，具体是指一种煤储层现今地应力预测方法。
背景技术
随着油气勘探开发的深入，煤层气、页岩气以及致密油气等非常规油气逐渐显示出巨大潜力。煤层气是一种以吸附状态为主赋存于煤层及其围岩中的非常规天然气，其全球资源量约为256.3×1012m3，中国埋深2000m以内的煤层气地质资源量为36.8×1012m3，其中可采资源量为10.9×1012m3。
现今地应力状态对煤层气勘探开发具有重要影响：(1)煤层渗透率是制约煤层气开发的关键因素之一，其变化对现今地应力非常敏感，明显受地应力控制。低地应力区煤层渗透率高，随深度增加渗透率变化幅度不大，煤层气单井产气量高；而高地应力区煤层渗透率低，随深度增加渗透率急剧减小，单井产气量低。(2)现今地应力影响煤层天然裂缝的有效性，当区域最大主应力方向与煤层中裂缝的优势方向一致时，裂缝受拉张而呈开启状态；而当区域最大主应力方向与煤层中裂缝优势方向垂直时，裂缝受挤压而呈闭合状态。(3)现今地应力控制水力压裂裂缝的形态及其扩展规律，水力压裂时，裂缝总是趋于弱面形成并扩展。在其它条件相同前提下，三轴应力差别不大时，在平面和剖面上都容易形成近似圆形的裂缝；而三轴应力差别越大时，形成的裂缝长轴与短轴差别越大，三轴应力椭球体的形态大致反映出裂缝延伸的形态。煤层水力压裂缝的形态主要受地应力大小、方向以及岩石性质等控制。(4)现今地应力影响煤层气的吸附、解吸、扩散和渗流，从而影响煤层气的产出。(5)现今地应力影响煤层气开发方案设计，是煤层气储层开发井网部署的关键参数。
因此，准确预测煤储层现今地应力状态至关重要。目前，针对煤层现今地应力的预测，现有技术主要包括两种：(1)对常规砂泥岩储层地应力的预测研究已较为深入，对煤层现今地应力的预测仍使用砂泥岩储层地应力预测方法，没有考虑到煤岩具有低杨氏模量、高泊松比、低强度、易变形等力学特点，与常规砂泥岩相比差异明显；(2)在煤储层地应力预测时注意到煤岩易变形等特点，但在其预测方法中现今地应力因素考虑不全面，如：煤层变形受温度影响，温度变化产生附加应力。由此表明，无论使用上述哪种预测方法，其结果与实际必然存在偏差。
发明内容
为解决上述技术问题，本发明提供的技术方案为：一种煤储层现今地应力预测方法，包括以下步骤：
(1)采集钻孔岩芯中煤层上覆不同时代的岩石样品，测试获取古应力值，确定最后期强构造活动的地质时间；
(2)采集钻孔岩芯煤岩样品，测试获取其静态杨氏模量、静态泊松比、静态粘性系数和静态Biot系数；
(3)获取声波时差测井曲线和密度测井曲线，计算煤岩动态杨氏模量、动态泊松比和动态Biot系数；
(4)基于步骤(2)实验和步骤(3)计算结果，构建煤岩动态/静态杨氏模量、动态/静态泊松比、动态/静态Biot系数、静态粘性系数/动态杨氏模量之间的关系，实现基于测井曲线表征煤岩岩石力学参数的目的；
(5)基于煤储层试井分析获取应力参数，结合步骤(3)的密度测井曲线，计算实测深度处的现今地应力；
(6)基于步骤(5)实测储层压力和步骤(3)测井曲线，反演计算伊顿系数，实现基于测井曲线表征煤储层压力的目的；
(7)综合步骤(1)、步骤(4)、步骤(5)和步骤(6)的分析计算结果，反演获取煤储层现今地应力模型中的未知参数，构建基于Maxwell粘弹性体的煤储层现今地应力预测模型，利用测井曲线准确计算煤储层地应力连续剖面，实现预测煤储层现今地应力的目的。
优选地，所述步骤(1)包括以下子步骤：
(1.1)确定煤层发育层位，在钻孔中采集煤层上部不同时代岩石样品，在实验室内加工成50mm×25mm(长度×直径)的标准圆柱塞状；
(1.2)通过岩石声发射法测试Kaiser点，获取不同时期古应力值，确定最后期强构造活动的地质时间，记至现今经历时间为t。
优选地，所述步骤(2)包括以下子步骤：
(2.1)按照国家标准GB/T 23561.9-2009《煤和岩石物理力学性质测定方法》，测试获取煤岩静态杨氏模量(Es)和静态泊松比(μs)；
(2.2)基于Cross-plotting法测试获取煤岩样品的静态Biot系数(αs)；
(2.3)基于煤岩的三角波周期加载，反推确定Maxwell体的静态粘性系数(η)。
优选地，所述步骤(3)包括以下子步骤：
(3.1)获取声波时差测井和密度测井曲线，对其进行校正；
(3.2)计算煤岩的动态杨氏模量(Ed)、动态泊松比(μd)和动态Biot系数(αd)；
其中，煤岩的动态杨氏模量(Ed)计算公式为：
式中：Ed为煤岩动态杨氏模量，ρ为煤岩密度，Δts为横波时差，Δtp为纵波时差，β为单位转换系数；
煤岩动态泊松比(μd)计算公式为：
式中：μd为煤岩动态泊松比，Δts为横波时差，Δtp为纵波时差；
煤岩动态Biot系数(αd)计算公式为：
式中：αd为煤岩动态Biot系数，ρ和ρma分别为煤岩和煤岩骨架密度，νp和νmp分别为煤岩和煤岩骨架纵波速度，νs和νms分别为煤岩和煤岩骨架横波速度。
优选地，所述步骤(4)包括以下子步骤：
(4.1)构建煤岩动态/静态杨氏模量关系，记为：
Es＝f1(Ed) (4)
(4.2)构建煤岩动态/静态泊松比关系，记为：
μs＝f2(μd) (5)
(4.3)构建煤岩动态/静态Biot系数关系，记为：
αs＝f3(αd) (6)
(4.4)构建煤岩静态粘性系数/动态杨氏模量关系，记为：
η＝f4(Ed) (7)
优选地，所述步骤(5)包括以下子步骤：
(5.1)基于煤储层试井分析，获取实测处埋深(z)，储层压力(Po)，破裂压力(Pf)，闭合压力(Pc)和重张压力(Pr)；
(5.2)依据步骤(3)的密度测井曲线，对其积分获得实测深度处的垂向应力，其计算公式为：
式中：Sv为垂向应力，ρ为煤岩密度，h为深度，z为实测处埋深，g为重力加速度；
(5.3)计算实测深度处的水平最大主应力和水平最小主应力，其计算公式为：
式中：Sv为垂向应力，Sh和SH分别为水平最小和最大主应力，T为煤岩抗张强度，Po为储层压力，Pf为破裂压力，Pc为闭合压力，Pr为重张压力。
优选地，所述步骤(6)中的伊顿系数计算公式为：
式中：Sv为垂向应力，Po为储层压力，S为静水压力，Δt为计算点的声波时差，Δtn为计算点对应的正常趋势线上的声波时差，c为伊顿系数。
优选地，所述步骤(7)包括以下子步骤：
(7.1)选取Maxwell体表征煤岩，Maxwell体由一个弹簧和一个阻尼器串联而成，可反映岩石应力应变的时间特性，有效表征煤岩的易变形等特征；
(7.2)综合步骤(1)、步骤(4)、步骤(5)和步骤(6)的分析计算结果，反推获取水平最大主应力和水平最小主应力方向上的应变(εH和εh)和应力附加项(ST和St)，其中，应力附加项主要用于校正由温度等其它因素造成的附加应力，应变和应力附加项均采取平均值处理方式：
式中：n为样品总数，εH和εh为水平最大主应力和水平最小主应力方向上的平均应变，ST和St为水平最大主应力和水平最小主应力方向上的应力附加项平均值；
(7.3)构建带有应力附加项的Maxwell煤储层地应力模型，其计算公式如下：
式中：Sv为垂向应力，Sh和SH分别为水平最小和最大主应力，μs为静态泊松比，Es为静态杨氏模量，αs为静态Biot系数，Po为储层压力，εH和εh为水平最大主应力和水平最小主应力方向上的平均应变，ST和St为水平最大主应力和水平最小主应力方向上的应力附加项平均值，η为Maxwell体的静态粘性系数，t为最后期强构造活动至现今经历时间；
(7.4)公式(12)中垂向应力Sv可通过公式(8)计算，静态参数可通过公式(4)～公式(7)计算，储层压力Po可通过公式(10)计算，平均应变(εH和εh)和应力附加项平均值(ST和St)可由前述步骤反推后取平均值确定，最后期强构造活动至现今经历时间t通过多组声发射实验结果确定，由此，构建基于Maxwell粘弹性体的煤储层现今地应力预测模型，利用测井曲线准确计算煤储层地应力连续剖面，实验预测煤储层现今地应力的目的。
采用以上预测方法后，本发明具有如下优点：
(1)本发明在煤岩参数实验测试和测井曲线计算基础上，建立静态参数与动态参数之间的量化关系，考虑煤岩低杨氏模量、高泊松比、低强度、易变形等岩石力学特点，反演获取基于Maxwell煤储层现今地应力模型中的未知参数，建立基于测井曲线的煤储层地应力剖面；
(2)本发明针对性强，可准确预测煤储层现今地应力，为煤层气有效勘探开发奠定基础，减小风险和成本。
附图说明
图1是一种煤储层现今地应力预测方法的流程示意图。
具体实施方式
下面对本发明做进一步的详细说明。
结合附图1，一种煤储层现今地应力预测方法，包括以下步骤：
(1)采集钻孔岩芯中煤层上覆不同时代的岩石样品，测试获取古应力值，确定最后期强构造活动的地质时间，包括以下子步骤：
(1.1)确定煤层发育层位，在钻孔中采集煤层上部不同时代岩石样品，在实验室内加工成50mm×25mm(长度×直径)的标准圆柱塞状；
(1.2)通过岩石声发射法测试Kaiser点，获取不同时期古应力值，确定最后期强构造活动的地质时间，记至现今经历时间为t；
(2)采集钻孔岩芯煤岩样品，测试获取其静态杨氏模量、静态泊松比、静态粘性系数和静态Biot系数，包括以下子步骤：
(2.1)按照国家标准GB/T 23561.9-2009《煤和岩石物理力学性质测定方法》，测试获取煤岩静态杨氏模量(Es)和静态泊松比(μs)；
(2.2)基于Cross-plotting法测试获取煤岩样品的静态Biot系数(αs)；
(2.3)基于煤岩的三角波周期加载，反推确定Maxwell体的静态粘性系数(η)；
(3)获取声波时差测井曲线和密度测井曲线，计算煤岩动态杨氏模量、动态泊松比和动态Biot系数，包括以下子步骤：
(3.1)获取声波时差测井和密度测井曲线，对其进行校正；
(3.2)计算煤岩的动态杨氏模量(Ed)、动态泊松比(μd)和动态Biot系数(αd)；
其中，煤岩的动态杨氏模量(Ed)计算公式为：
式中：Ed为煤岩动态杨氏模量，ρ为煤岩密度，Δts为横波时差，Δtp为纵波时差，β为单位转换系数；
煤岩动态泊松比(μd)计算公式为：
式中：μd为煤岩动态泊松比，Δts为横波时差，Δtp为纵波时差；
煤岩动态Biot系数(αd)计算公式为：
式中：αd为煤岩动态Biot系数，ρ和ρma分别为煤岩和煤岩骨架密度，νp和νmp分别为煤岩和煤岩骨架纵波速度，νs和νms分别为煤岩和煤岩骨架横波速度；
(4)基于步骤(2)实验和步骤(3)计算结果，构建煤岩动态/静态杨氏模量、动态/静态泊松比、动态/静态Biot系数、静态粘性系数/动态杨氏模量之间的关系，实现基于测井曲线表征煤岩岩石力学参数，包括以下子步骤：
(4.1)构建煤岩动态/静态杨氏模量关系，记为：
Es＝f1(Ed) (4)
(4.2)构建煤岩动态/静态泊松比关系，记为：
μs＝f2(μd) (5)
(4.3)构建煤岩动态/静态Biot系数关系，记为：
αs＝f3(αd) (6)
(4.4)构建煤岩静态粘性系数/动态杨氏模量关系，记为：
η＝f4(Ed) (7)
(5)基于煤储层试井分析获取应力参数，结合步骤(3)的密度测井曲线，计算实测深度处的现今地应力，包括以下子步骤：
(5.1)基于煤储层试井分析，获取实测处埋深(z)，储层压力(Po)，破裂压力(Pf)，闭合压力(Pc)和重张压力(Pr)；
(5.2)依据步骤(3)的密度测井曲线，对其积分获得实测深度处的垂向应力，其计算公式为：
式中：Sv为垂向应力，ρ为煤岩密度，h为深度，z为实测处埋深，g为重力加速度；
(5.3)计算实测深度处的水平最大主应力和水平最小主应力，其计算公式为：
式中：Sv为垂向应力，Sh和SH分别为水平最小和最大主应力，T为煤岩抗张强度，Po为储层压力，Pf为破裂压力，Pc为闭合压力，Pr为重张压力；
(6)基于步骤(5)实测储层压力和步骤(3)测井曲线，反演计算伊顿系数，实现基于测井曲线表征煤储层压力，伊顿系数计算公式为：
式中：Sv为垂向应力，Po为储层压力，S为静水压力，Δt为计算点的声波时差，Δtn为计算点对应的正常趋势线上的声波时差，c为伊顿系数；
(7)综合步骤(1)、步骤(4)、步骤(5)和步骤(6)的分析计算结果，反演获取煤储层现今地应力模型中的未知参数，构建基于Maxwell粘弹性体的煤储层现今地应力预测模型，利用测井曲线准确计算煤储层地应力连续剖面，实现预测煤储层现今地应力的目的，包括以下子步骤：
(7.1)选取Maxwell体表征煤岩，Maxwell体由一个弹簧和一个阻尼器串联而成，可反映岩石应力应变的时间特性，有效表征煤岩的易变形等特征；
(7.2)综合步骤(1)、步骤(4)、步骤(5)和步骤(6)的分析计算结果，反推获取水平最大主应力和水平最小主应力方向上的应变(εH和εh)和应力附加项(ST和St)，其中，应力附加项主要用于校正由温度等其它因素造成的附加应力，应变和应力附加项均采取平均值处理方式：
式中：n为样品总数，εH和εh为水平最大主应力和水平最小主应力方向上的平均应变，ST和St为水平最大主应力和水平最小主应力方向上的应力附加项平均值；
(7.3)构建带有应力附加项的Maxwell煤储层地应力模型，其计算公式如下：
式中：Sv为垂向应力，Sh和SH分别为水平最小和最大主应力，μs为静态泊松比，Es为静态杨氏模量，αs为静态Biot系数，Po为储层压力，εH和εh为水平最大主应力和水平最小主应力方向上的平均应变，ST和St为水平最大主应力和水平最小主应力方向上的应力附加项平均值，η为Maxwell体的静态粘性系数，t为最后期强构造活动至现今经历时间；
(7.4)公式(12)中垂向应力Sv可通过公式(8)计算，静态参数可通过公式(4)～公式(7)计算，储层压力Po可通过公式(10)计算，平均应变(εH和εh)和应力附加项平均值(ST和St)可由前述步骤反推后取平均值确定，最后期强构造活动至现今经历时间t通过多组声发射实验结果确定，由此，构建基于Maxwell粘弹性体的煤储层现今地应力预测模型，利用测井曲线准确计算煤储层地应力连续剖面，实现预测煤储层现今地应力的目的。
Claims (8)
1.一种煤储层现今地应力预测方法，其特征在于，包括以下步骤：
(1)采集钻孔岩芯中煤层上覆不同时代的岩石样品，测试获取古应力值，确定最后期强构造活动的地质时间；
(2)采集钻孔岩芯煤岩样品，测试获取其静态杨氏模量、静态泊松比、静态粘性系数和静态Biot系数；
(3)获取声波时差测井曲线和密度测井曲线，计算煤岩动态杨氏模量、动态泊松比和动态Biot系数；
(4)基于步骤(2)实验和步骤(3)计算结果，构建煤岩动态/静态杨氏模量、动态/静态泊松比、动态/静态Biot系数、静态粘性系数/动态杨氏模量之间的关系，实现基于测井曲线表征煤岩岩石力学参数的目的；
(5)基于煤储层试井分析获取应力参数，结合步骤(3)的密度测井曲线，计算实测深度处的现今地应力；
(6)基于步骤(5)实测储层压力和步骤(3)测井曲线，反演计算伊顿系数，实现基于测井曲线表征煤储层压力的目的；
(7)综合步骤(1)、步骤(4)、步骤(5)和步骤(6)的分析计算结果，反演获取煤储层现今地应力模型中的未知参数，构建基于Maxwell粘弹性体的煤储层现今地应力预测模型，利用测井曲线准确计算煤储层地应力连续剖面，实现预测煤储层现今地应力的目的。
2.根据权利要求1所述的一种煤储层现今地应力预测方法，其特征在于，所述步骤(1)包括以下子步骤：
(1.1)确定煤层发育层位，在钻孔中采集煤层上部不同时代岩石样品，在实验室内加工成长度×直径为50mm×25mm的标准圆柱塞状；
(1.2)通过岩石声发射法测试Kaiser点，获取不同时期古应力值，确定最后期强构造活动的地质时间，记至现今经历时间为t。
3.根据权利要求1所述的一种煤储层现今地应力预测方法，其特征在于，所述步骤(2)包括以下子步骤：
(2.1)按照国家标准GB/T 23561.9-2009《煤和岩石物理力学性质测定方法》，测试获取煤岩静态杨氏模量Es和静态泊松比μs；
(2.2)基于Cross-plotting法测试获取煤岩样品的静态Biot系数αs；
(2.3)基于煤岩的三角波周期加载，反推确定Maxwell体的静态粘性系数η。
4.根据权利要求1所述的一种煤储层现今地应力预测方法，其特征在于，所述步骤(3)包括以下子步骤：
(3.1)获取声波时差测井和密度测井曲线，对其进行校正；
(3.2)计算煤岩的动态杨氏模量Ed、动态泊松比μd和动态Biot系数αd；
其中，煤岩的动态杨氏模量Ed计算公式为：
式中：Ed为煤岩动态杨氏模量，ρ为煤岩密度，Δts为横波时差，Δtp为纵波时差，β为单位转换系数；
煤岩动态泊松比μd计算公式为：
式中：μd为煤岩动态泊松比，Δts为横波时差，Δtp为纵波时差；
煤岩动态Biot系数αd计算公式为：
式中：αd为煤岩动态Biot系数，ρ和ρma分别为煤岩和煤岩骨架密度，vp和vmp分别为煤岩和煤岩骨架纵波速度，vs和vms分别为煤岩和煤岩骨架横波速度。
5.根据权利要求1所述的一种煤储层现今地应力预测方法，其特征在于，所述步骤(4)包括以下子步骤：
(4.1)构建煤岩动态/静态杨氏模量关系，记为：
Es＝f1(Ed) (4)
(4.2)构建煤岩动态/静态泊松比关系，记为：
μs＝f2(μd) (5)
(4.3)构建煤岩动态/静态Biot系数关系，记为：
αs＝f3(αd) (6)
(4.4)构建煤岩静态粘性系数/动态杨氏模量关系，记为：
η＝f4(Ed) (7)。
6.根据权利要求1所述的一种煤储层现今地应力预测方法，其特征在于，所述步骤(5)包括以下子步骤：
(5.1)基于煤储层试井分析，获取实测处埋深z，储层压力Po，破裂压力Pf，闭合压力Pc和重张压力Pr；
(5.2)依据步骤(3)的密度测井曲线，对其积分获得实测深度处的垂向应力，其计算公式为：
式中：Sv为垂向应力，ρ为煤岩密度，h为深度，z为实测处埋深，g为重力加速度；
(5.3)计算实测深度处的水平最大主应力和水平最小主应力，其计算公式为：
式中：Sv为垂向应力，Sh和SH分别为水平最小和最大主应力，T为煤岩抗张强度，Po为储层压力，Pf为破裂压力，Pc为闭合压力，Pr为重张压力。
7.根据权利要求1所述的一种煤储层现今地应力预测方法，其特征在于，所述步骤(6)中的伊顿系数计算公式为：
式中：Sv为垂向应力，Po为储层压力，S为静水压力，Δt为计算点的声波时差，Δtn为计算点对应的正常趋势线上的声波时差，c为伊顿系数。
8.根据权利要求1所述的一种煤储层现今地应力预测方法，其特征在于，所述步骤(7)包括以下子步骤：
(7.1)选取Maxwell体表征煤岩，Maxwell体由一个弹簧和一个阻尼器串联而成，可反映岩石应力应变的时间特性，有效表征煤岩的易变形特征；
(7.2)综合步骤(1)、步骤(4)、步骤(5)和步骤(6)的分析计算结果，反推获取水平最大主应力和水平最小主应力方向上的应变εH、εh和应力附加项ST和St，其中，应力附加项用于校正由温度造成的附加应力，应变和应力附加项均采取平均值处理方式：
式中：n为样品总数，εH和εh分别为水平最大主应力和水平最小主应力方向上的应变，ST和St为应力附加项；
(7.3)构建带有应力附加项的Maxwell煤储层地应力模型，其计算公式如下：
式中：Sv为垂向应力，Sh和SH分别为水平最小和最大主应力，μs为静态泊松比，Es为静态杨氏模量，αs为静态Biot系数，Po为储层压力，εH和εh分别为水平最大主应力和水平最小主应力方向上的应变，ST和St为应力附加项，η为Maxwell体的静态粘性系数，t为最后期强构造活动至现今经历时间；
(7.4)公式(12)中垂向应力Sv可通过公式(8)计算，静态参数可通过公式(4)～公式(7)计算，储层压力Po可通过公式(10)计算，水平最大主应力和水平最小主应力方向上的应变εH和εh、应力附加项平均值ST和St可由前述步骤反推后取平均值确定，最后期强构造活动至现今经历时间t通过多组声发射实验结果确定，由此，构建基于Maxwell粘弹性体的煤储层现今地应力预测模型，利用测井曲线准确计算煤储层地应力连续剖面，实现预测煤储层现今地应力的目的。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810869778.4A CN109113742B (zh) | 2018-08-02 | 2018-08-02 | 一种煤储层现今地应力预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810869778.4A CN109113742B (zh) | 2018-08-02 | 2018-08-02 | 一种煤储层现今地应力预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109113742A CN109113742A (zh) | 2019-01-01 |
| CN109113742B true CN109113742B (zh) | 2019-12-31 |
Family
ID=64851705
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810869778.4A Active CN109113742B (zh) | 2018-08-02 | 2018-08-02 | 一种煤储层现今地应力预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109113742B (zh) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110006568B (zh) * | 2019-04-12 | 2020-01-07 | 中国地质科学院地质力学研究所 | 一种利用岩芯获取三维地应力的方法及获取系统 |
| CN110598244B (zh) * | 2019-07-26 | 2021-03-12 | 中国矿业大学 | 一种瓦斯气团生命周期预测及基于其的瓦斯浓度预测方法 |
| CN110598178B (zh) * | 2019-08-08 | 2023-04-28 | 大庆油田有限责任公司 | 一种油藏现今地应力快速计算方法 |
| CN110727031B (zh) * | 2019-11-18 | 2021-04-13 | 科吉思石油技术咨询(北京)有限公司 | 一种基于三维叠前地震反演结果的地应力获取方法 |
| CN111520125B (zh) * | 2020-05-07 | 2021-11-12 | 中国矿业大学 | 一种褶皱发育区现今地应力方位转向的定量分析方法 |
| CN111537332B (zh) * | 2020-05-21 | 2021-04-02 | 中国矿业大学 | 一种考虑温度效应的深层-超深层岩石力学参数预测方法 |
| CN112412434B (zh) * | 2020-06-02 | 2024-03-26 | 中国海洋石油集团有限公司 | 一种改进的疏松砂岩地应力计算方法 |
| CN112379416B (zh) * | 2020-10-13 | 2024-02-06 | 北京恒泰兴科信息技术有限公司 | 煤岩岩石物理建模预测横波方法、装置及电子设备 |
| CN113361131A (zh) * | 2021-06-25 | 2021-09-07 | 中国矿业大学 | 一种原生结构煤的静态力学参数计算方法 |
| CN113779664A (zh) * | 2021-08-02 | 2021-12-10 | 清能艾科（深圳）能源技术有限公司 | 单井地应力预测方法、装置、电子设备及存储介质 |
| CN114294059B (zh) * | 2021-12-29 | 2024-05-07 | 中煤科工开采研究院有限公司 | 一种坚硬岩层水力压裂层位高精度控制方法 |
| CN116401886B (zh) * | 2023-04-17 | 2023-11-24 | 应急管理部国家自然灾害防治研究院 | 一种适用于川藏地区的快速地应力估算方法 |
| CN117647342B (zh) * | 2024-01-30 | 2024-04-19 | 中国地质大学（北京） | 一种基于声发射小波分析的煤储层地应力测定方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1818337A (zh) * | 2006-03-27 | 2006-08-16 | 天地科技股份有限公司 | 煤矿井下地质力学参数快速测试方法及装备 |
| CN101526009A (zh) * | 2009-04-09 | 2009-09-09 | 西安科技大学 | 围岩失稳声光电集成监测系统及其监测方法 |
| CN103206206A (zh) * | 2012-01-16 | 2013-07-17 | 中国石油化工股份有限公司 | 用于计算超深井地层破裂压力的方法 |
| KR101450504B1 (ko) * | 2014-04-25 | 2014-10-15 | 한국지질자원연구원 | 암석의 역학적 특성을 이용한 채석면의 할석 용이도 평가 방법 |
| AU2014202974A1 (en) * | 2013-06-07 | 2015-01-15 | Sandvik Mining And Construction Oy | Arrangement and method in rock breaking |
| CN106150501A (zh) * | 2016-06-30 | 2016-11-23 | 太原理工大学 | 利用上覆采空区废水原位弱化主控煤岩层的方法 |
| CN107038313A (zh) * | 2017-04-26 | 2017-08-11 | 东北大学 | 基于数值岩芯的分层地应力精细描述方法 |
-
2018
- 2018-08-02 CN CN201810869778.4A patent/CN109113742B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1818337A (zh) * | 2006-03-27 | 2006-08-16 | 天地科技股份有限公司 | 煤矿井下地质力学参数快速测试方法及装备 |
| CN101526009A (zh) * | 2009-04-09 | 2009-09-09 | 西安科技大学 | 围岩失稳声光电集成监测系统及其监测方法 |
| CN103206206A (zh) * | 2012-01-16 | 2013-07-17 | 中国石油化工股份有限公司 | 用于计算超深井地层破裂压力的方法 |
| AU2014202974A1 (en) * | 2013-06-07 | 2015-01-15 | Sandvik Mining And Construction Oy | Arrangement and method in rock breaking |
| KR101450504B1 (ko) * | 2014-04-25 | 2014-10-15 | 한국지질자원연구원 | 암석의 역학적 특성을 이용한 채석면의 할석 용이도 평가 방법 |
| CN106150501A (zh) * | 2016-06-30 | 2016-11-23 | 太原理工大学 | 利用上覆采空区废水原位弱化主控煤岩层的方法 |
| CN107038313A (zh) * | 2017-04-26 | 2017-08-11 | 东北大学 | 基于数值岩芯的分层地应力精细描述方法 |
Non-Patent Citations (3)
| Title |
|---|
| 地应力对煤储层渗透性影响及其控制机理研;刘大锰等;《煤炭科技技术》;20170630;全文 * |
| 沁水盆地郑庄区块地应力发育特征及其地质意义;李叶朋等;《煤炭科学技术》;20171030;全文 * |
| 煤层气储层地应力场宏观分布规律统计分析;陈世达等;<煤炭科学技术>;20180630;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109113742A (zh) | 2019-01-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109113742B (zh) | 一种煤储层现今地应力预测方法 | |
| US10428626B2 (en) | Production estimation in subterranean formations | |
| CN108009705A (zh) | 一种基于支持向量机技术的页岩储层可压性评价方法 | |
| CN105277982B (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
| Yang et al. | Viscoplastic deformation of the Bakken and adjacent formations and its relation to hydraulic fracture growth | |
| CN108304959B (zh) | 提高地层流体压力预测精度的方法 | |
| CN110969699B (zh) | 一种高精度三维地应力建模方法及其应用 | |
| CN113640119B (zh) | 一种确定应力相关岩石动态Biot系数的方法 | |
| CN113126167A (zh) | 一种基于测井资料的井周近井壁地层损伤评价方法 | |
| CN112412434B (zh) | 一种改进的疏松砂岩地应力计算方法 | |
| CN110727031B (zh) | 一种基于三维叠前地震反演结果的地应力获取方法 | |
| CN105929448B (zh) | 一种原状泥页岩地层声波时差响应反演方法 | |
| CN110320562B (zh) | 页岩气储层水平井中声波时差的校正方法 | |
| CN109975189B (zh) | 孔隙型砂岩储层产能预测方法与装置 | |
| CN106019378B (zh) | 一种时移测井曲线的动态重构方法 | |
| CN110671102B (zh) | 一种气井临界出砂压差的确定方法及系统 | |
| CN108195669B (zh) | 油藏围压条件下岩石静态力学参数的矫正及预测方法 | |
| Eide | On shear wave velocity testing in clay | |
| CN113777668A (zh) | 用于砂泥岩互层致密气藏的地应力计算方法及装置 | |
| CN111596377A (zh) | 一种高地应力软岩隧道松动圈的联合测试方法 | |
| CN111257134A (zh) | 剪切应力的连续深度处理方法 | |
| CN110427688A (zh) | 一种基于实测振动的地应力大小预测方法 | |
| Escobar-Segovia et al. | New challenge to optimize wells drilling using the geomechanical model design | |
| RU2797376C1 (ru) | Способ определения трещинного коллектора и способ добычи углеводородов | |
| Su et al. | In-situ Stress Evaluation of Deep Carbonate Reservoir in Block S of Tarim Basin |
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 |