CN109113742A - 一种煤储层现今地应力预测方法 - Google Patents
一种煤储层现今地应力预测方法 Download PDFInfo
- Publication number
- CN109113742A CN109113742A CN201810869778.4A CN201810869778A CN109113742A CN 109113742 A CN109113742 A CN 109113742A CN 201810869778 A CN201810869778 A CN 201810869778A CN 109113742 A CN109113742 A CN 109113742A
- Authority
- CN
- China
- Prior art keywords
- coal
- stress
- coal petrography
- coal seam
- static
- 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
- 239000003245 coal Substances 0.000 title claims abstract description 182
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000003068 static effect Effects 0.000 claims abstract description 60
- 239000011435 rock Substances 0.000 claims abstract description 22
- 238000012360 testing method Methods 0.000 claims abstract description 19
- 238000002474 experimental method Methods 0.000 claims abstract description 5
- 238000005259 measurement Methods 0.000 claims description 18
- 238000004458 analytical method Methods 0.000 claims description 12
- 230000000694 effects Effects 0.000 claims description 8
- 238000012512 characterization method Methods 0.000 claims description 6
- 238000005553 drilling Methods 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000002706 hydrostatic effect Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 7
- 230000035699 permeability Effects 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001179 sorption measurement 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分别为煤岩和煤岩骨架密度，νp和νmp分别为煤岩和煤岩骨架纵波速度，νs和νms分别为煤岩和煤岩骨架横波速度。
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 true CN109113742A (zh) | 2019-01-01 |
| CN109113742B 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) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110006568A (zh) * | 2019-04-12 | 2019-07-12 | 中国地质科学院地质力学研究所 | 一种利用岩芯获取三维地应力的方法及其获取系统 |
| CN110598178A (zh) * | 2019-08-08 | 2019-12-20 | 大庆油田有限责任公司 | 一种油藏现今地应力快速计算方法 |
| CN110598244A (zh) * | 2019-07-26 | 2019-12-20 | 中国矿业大学 | 一种瓦斯气团生命周期预测及基于其的瓦斯浓度预测方法 |
| CN110727031A (zh) * | 2019-11-18 | 2020-01-24 | 科吉思石油技术咨询(北京)有限公司 | 一种基于三维叠前地震反演结果的地应力获取方法 |
| CN111520125A (zh) * | 2020-05-07 | 2020-08-11 | 中国矿业大学 | 一种褶皱发育区现今地应力方位转向的定量分析方法 |
| CN112379416A (zh) * | 2020-10-13 | 2021-02-19 | 北京恒泰兴科信息技术有限公司 | 煤岩岩石物理建模预测横波方法、装置及电子设备 |
| CN112412434A (zh) * | 2020-06-02 | 2021-02-26 | 中国海洋石油集团有限公司 | 一种改进的疏松砂岩地应力计算方法 |
| CN111537332B (zh) * | 2020-05-21 | 2021-04-02 | 中国矿业大学 | 一种考虑温度效应的深层-超深层岩石力学参数预测方法 |
| CN113361131A (zh) * | 2021-06-25 | 2021-09-07 | 中国矿业大学 | 一种原生结构煤的静态力学参数计算方法 |
| CN113779664A (zh) * | 2021-08-02 | 2021-12-10 | 清能艾科（深圳）能源技术有限公司 | 单井地应力预测方法、装置、电子设备及存储介质 |
| CN114294059A (zh) * | 2021-12-29 | 2022-04-08 | 中煤科工开采研究院有限公司 | 一种坚硬岩层水力压裂层位高精度控制方法 |
| CN116401886A (zh) * | 2023-04-17 | 2023-07-07 | 应急管理部国家自然灾害防治研究院 | 一种适用于川藏地区的快速地应力估算方法 |
| CN117647342A (zh) * | 2024-01-30 | 2024-03-05 | 中国地质大学（北京） | 一种基于声发射小波分析的煤储层地应力测定方法 |
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 |
|---|
| 刘大锰等: "地应力对煤储层渗透性影响及其控制机理研", 《煤炭科技技术》 * |
| 李叶朋等: "沁水盆地郑庄区块地应力发育特征及其地质意义", 《煤炭科学技术》 * |
| 陈世达等: "煤层气储层地应力场宏观分布规律统计分析", <煤炭科学技术> * |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110006568A (zh) * | 2019-04-12 | 2019-07-12 | 中国地质科学院地质力学研究所 | 一种利用岩芯获取三维地应力的方法及其获取系统 |
| CN110598244A (zh) * | 2019-07-26 | 2019-12-20 | 中国矿业大学 | 一种瓦斯气团生命周期预测及基于其的瓦斯浓度预测方法 |
| CN110598244B (zh) * | 2019-07-26 | 2021-03-12 | 中国矿业大学 | 一种瓦斯气团生命周期预测及基于其的瓦斯浓度预测方法 |
| CN110598178B (zh) * | 2019-08-08 | 2023-04-28 | 大庆油田有限责任公司 | 一种油藏现今地应力快速计算方法 |
| CN110598178A (zh) * | 2019-08-08 | 2019-12-20 | 大庆油田有限责任公司 | 一种油藏现今地应力快速计算方法 |
| CN110727031A (zh) * | 2019-11-18 | 2020-01-24 | 科吉思石油技术咨询(北京)有限公司 | 一种基于三维叠前地震反演结果的地应力获取方法 |
| CN111520125A (zh) * | 2020-05-07 | 2020-08-11 | 中国矿业大学 | 一种褶皱发育区现今地应力方位转向的定量分析方法 |
| CN111520125B (zh) * | 2020-05-07 | 2021-11-12 | 中国矿业大学 | 一种褶皱发育区现今地应力方位转向的定量分析方法 |
| CN111537332B (zh) * | 2020-05-21 | 2021-04-02 | 中国矿业大学 | 一种考虑温度效应的深层-超深层岩石力学参数预测方法 |
| CN112412434A (zh) * | 2020-06-02 | 2021-02-26 | 中国海洋石油集团有限公司 | 一种改进的疏松砂岩地应力计算方法 |
| CN112412434B (zh) * | 2020-06-02 | 2024-03-26 | 中国海洋石油集团有限公司 | 一种改进的疏松砂岩地应力计算方法 |
| CN112379416A (zh) * | 2020-10-13 | 2021-02-19 | 北京恒泰兴科信息技术有限公司 | 煤岩岩石物理建模预测横波方法、装置及电子设备 |
| CN112379416B (zh) * | 2020-10-13 | 2024-02-06 | 北京恒泰兴科信息技术有限公司 | 煤岩岩石物理建模预测横波方法、装置及电子设备 |
| CN113361131A (zh) * | 2021-06-25 | 2021-09-07 | 中国矿业大学 | 一种原生结构煤的静态力学参数计算方法 |
| CN113779664A (zh) * | 2021-08-02 | 2021-12-10 | 清能艾科（深圳）能源技术有限公司 | 单井地应力预测方法、装置、电子设备及存储介质 |
| CN114294059A (zh) * | 2021-12-29 | 2022-04-08 | 中煤科工开采研究院有限公司 | 一种坚硬岩层水力压裂层位高精度控制方法 |
| CN114294059B (zh) * | 2021-12-29 | 2024-05-07 | 中煤科工开采研究院有限公司 | 一种坚硬岩层水力压裂层位高精度控制方法 |
| CN116401886A (zh) * | 2023-04-17 | 2023-07-07 | 应急管理部国家自然灾害防治研究院 | 一种适用于川藏地区的快速地应力估算方法 |
| CN116401886B (zh) * | 2023-04-17 | 2023-11-24 | 应急管理部国家自然灾害防治研究院 | 一种适用于川藏地区的快速地应力估算方法 |
| CN117647342A (zh) * | 2024-01-30 | 2024-03-05 | 中国地质大学（北京） | 一种基于声发射小波分析的煤储层地应力测定方法 |
| CN117647342B (zh) * | 2024-01-30 | 2024-04-19 | 中国地质大学（北京） | 一种基于声发射小波分析的煤储层地应力测定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109113742B (zh) | 2019-12-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109113742A (zh) | 一种煤储层现今地应力预测方法 | |
| Presti et al. | Damping ratio of soils from laboratory and in situ tests | |
| CN105277982B (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
| CN105182424B (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
| CN104698492B (zh) | 一种计算异常地层压力的方法 | |
| CN101231346A (zh) | 应用地震波速度预测煤、岩体物理力学参数方法 | |
| CN103867197A (zh) | 复杂岩性天然气层声波时差判别法 | |
| CN110058323A (zh) | 一种致密砂岩地层脆性指数计算方法 | |
| CN101738637B (zh) | 一种基于速度随频率变化信息的油气检测方法 | |
| CN107728205B (zh) | 一种地层压力预测方法 | |
| CN109509111A (zh) | 探井地层压力的预测方法及系统 | |
| Song et al. | Hydraulic fracturing simulation of complex fractures growth in naturally fractured shale gas reservoir | |
| CN105317435A (zh) | 一种水平井裂缝识别方法 | |
| CN110727031B (zh) | 一种基于三维叠前地震反演结果的地应力获取方法 | |
| CN101315428A (zh) | 基于横波速度的高孔隙碎屑岩沉积地层的高分辨率勘探方法 | |
| CN108304959A (zh) | 提高地层流体压力预测精度的方法 | |
| Wang et al. | Effects of in-situ stress and joint on permeability of the coal bed in Linfen block, southeastern Ordos Basin, China | |
| RU2620822C1 (ru) | Способ выбора бурового раствора для строительства наклонных и горизонтальных скважин, пробуренных в неустойчивых глинистых отложениях | |
| CN110173259A (zh) | 一种斯通利波能量补偿校正的方法 | |
| CN103698494B (zh) | 一种确定岩性圈闭含油气饱和度的方法及装置 | |
| CN113126167A (zh) | 一种基于测井资料的井周近井壁地层损伤评价方法 | |
| Gonzalez et al. | 3D modeling of multistage hydraulic fractures and two-way-coupling geomechanics/fluid-flow simulation of a horizontal well in the nikanassin tight gas formation, Western Canada sedimentary basin | |
| Buntoro et al. | Validation of shale brittleness index calculation from wireline log of well BETRO-001 by using XRD test results and uniaxial test as parameters for determining potential of shale hydrocarbon-brown shale of Pematang Group Formation, Central Sumatra Basin, Bengkalis Trough | |
| Li et al. | In situ stress and mechanical properties of unconventional gas reservoirs | |
| Altowairqi et al. | Measuring ultrasonic characterisation to determine the impact of TOC and the stress field on shale gas anisotropy |
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 |