CN111520125A - 一种褶皱发育区现今地应力方位转向的定量分析方法 - Google Patents
一种褶皱发育区现今地应力方位转向的定量分析方法 Download PDFInfo
- Publication number
- CN111520125A CN111520125A CN202010375516.XA CN202010375516A CN111520125A CN 111520125 A CN111520125 A CN 111520125A CN 202010375516 A CN202010375516 A CN 202010375516A CN 111520125 A CN111520125 A CN 111520125A
- Authority
- CN
- China
- Prior art keywords
- fold
- stress
- ground stress
- distance
- drill hole
- 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
- 238000011161 development Methods 0.000 title claims abstract description 20
- 238000004445 quantitative analysis Methods 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000037303 wrinkles Effects 0.000 title claims abstract description 10
- 206010040954 Skin wrinkling Diseases 0.000 title description 6
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 17
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 13
- 238000005553 drilling Methods 0.000 claims abstract description 11
- 238000011160 research Methods 0.000 claims abstract description 10
- 238000003384 imaging method Methods 0.000 claims abstract description 7
- 238000012360 testing method Methods 0.000 claims abstract description 7
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 239000011435 rock Substances 0.000 claims description 17
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 3
- 230000006698 induction Effects 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 10
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000004898 kneading Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种褶皱发育区现今地应力方位转向的定量分析方法,适用于地质学与油气勘探开发领域。首先在发育褶皱的研究区设置垂直钻孔,收集钻孔资料,计算确定钻孔所在处的褶皱中和面深度,获取钻孔地应力测试信息以及成像测井钻井诱导缝与井壁崩落的解译,获取褶皱中和面深度处现今地应力在x轴向和y轴向的分量大小,记为区域现今地应力;根据地震数据体分析获取钻孔不同深度处的地层曲率,计算得到钻孔内距褶皱中和面距离为h处的应力分量大小,记为局部现今地应力;对区域现今地应力和局部现今地应力进行应力矢量合成,矢量合成的方向即为钻孔内距褶皱中和面距离为h处的现今地应力方位,由此实现定量分析现今地应力方向转向。其步骤简单,使用效果好。
Description
技术领域
本发明涉及一种定量分析方法,尤其涉及一种地质学与油气勘探开发领域使用的褶皱发育区现今地应力方位转向的定量分析方法。
背景技术
现今地应力是指目前存在于地壳中的内应力,是一种随着时间在不断变化的“场”。在构造地质学与油气勘探开发领域,开展现今地应力分析有助于深入了解地壳应力状态,指导油气勘探开发以及干热岩资源利用。
近些年,我国深层-超深层碎屑岩天然气勘探在塔里木盆地库车坳陷克拉苏构造带取得重大突破。克拉苏构造带油气资源丰富,是我国“西气东输”工程的重要气源地,目前已发现包括博孜9、克深8在内的多个埋深超过7000m的超深油气藏。
克拉苏构造带内储层地质力学评价发现,随着埋藏深度增加,现今地应力方位发生转向,但目前尚没有形成一种定量化分析褶皱发育区现今地应力方位转向的有效方法,因此,定量确定褶皱发育区现今地应力方位转向,可保障该区深层-超深层工程的安全实施,提升油气资源开发效率。
发明内容
针对上述技术的不足之处,提供一种步骤简单,检测效果好的褶皱发育区现今地应力方位转向的定量分析方法。
为实现上述技术目的,本发明的褶皱发育区现今地应力方位转向的定量分析方法,其步骤如下:
步骤1,在发育褶皱的研究区设置多个垂直钻孔,收集研究区任一钻孔A的常规测井资料,利用声波时差测井曲线、密度测井曲线计算钻孔A不同深度处岩石的杨氏模量和泊松比,根据杨氏模量大小粗化层系,以杨氏模量(Emax-Emin)/10为间隔粗化层系,即Emin~Emin+(Emax-Emin)/10为一个层系,Emin+(Emax-Emin)/10~Emin+(Emax-Emin)/5为另一个层系,以此类推,其中,Emax与Emin分别为杨氏模量最大值和最小值;
步骤2,基于薄板模型,利用步骤1获取的杨氏模量、泊松比和粗化层系厚度,计算确定钻孔A所在处的褶皱中和面深度,获取钻孔A地应力测试信息以及成像测井钻井诱导缝与井壁崩落的解译;
步骤3,根据钻孔A地应力测试、成像测井钻井诱导缝与井壁崩落的解译,获取褶皱中和面深度处现今地应力在x轴向和y轴向的分量大小,记为区域现今地应力Sr;
步骤4,收集研究区包含钻孔A的地震数据体,根据地震数据体分析获取钻孔A不同深度处的地层曲率1/rx和1/ry,1/rx为平行于x轴向的地层曲率,1/ry为平行于y轴向的地层曲率;
步骤5,利用公式:计算得到钻孔内距褶皱中和面距离为h处的应力分量大小σx和σy,记为局部现今地应力Sl,式中,h为钻孔内计算应力的深度点距褶皱中和面距离,σx为钻孔内距褶皱中和面距离为h处的x轴向应力,σy为钻孔内距褶皱中和面距离为h处的y轴向应力,Eh为钻孔内距褶皱中和面距离为h处的岩石杨氏模量,vh为钻孔内距褶皱中和面距离为h处的泊松比;
步骤6,利用公式:S=Sr+Sl,对区域现今地应力Sr和局部现今地应力Sl进行应力矢量合成,矢量合成的方向即为钻孔内距褶皱中和面距离为h处的现今地应力方位,由此,通过矢量合成距褶皱中和面不同距离h处的区域现今地应力和局部现今地应力,即可获取距褶皱中和面不同距离h处的现今地应力方位,实现定量分析现今地应力方向转向。
利用公式:
式中:htop为中和面到卷入褶皱中地层最顶面的距离;H为卷入褶皱中地层的总厚度;n为层系粗化后总层数;hi为第i层的地层厚度;Ei为第i层地层的杨氏模量,υi为第i层地层的泊松比。
有益效果:
本专利针对现今地应力方位随埋深发生转向的实际问题,在确定褶皱中和面位置基础上,获取区域应力,基于弹性薄板弯曲理论获取局部应力,矢量合成后确定卷入褶皱内任意深度的现今地应力方位,通过不同深度现今地应力方向的确定,进而定量分析现今地应力方位转向。本发明在现今地应力方位转向分析方面具较高实用价值,可操作性强,油气勘探开发实践性高。
附图说明
图1为本发明的褶皱发育区现今地应力方位转向的定量分析方法流程图;
具体实施方式
下面结合附图和实施例对本发明做进一步的说明。
如图1所示,本发明的褶皱发育区现今地应力方位转向的定量分析方法,其步骤如下:
步骤1,在发育褶皱的研究区设置多个垂直钻孔,收集研究区任一钻孔A的常规测井资料,利用声波时差测井曲线、密度测井曲线计算钻孔A不同深度处岩石的杨氏模量和泊松比,根据杨氏模量大小粗化层系,以杨氏模量(Emax-Emin)/10为间隔粗化层系,即Emin~Emin+(Emax-Emin)/10为一个层系,Emin+(Emax-Emin)/10~Emin+(Emax-Emin)/5为另一个层系,以此类推,其中,Emax与Emin分别为杨氏模量最大值和最小值;具体的,利用公式:计算岩石杨氏模量E,利用公式:计算岩石泊松比v,式中:ρ为岩石密度,Δts为声波横波时差,Δtp为声波纵波时差,α为单位转换系数;
步骤2,基于薄板模型,利用步骤1获取的杨氏模量、泊松比和粗化层系厚度,计算确定钻孔A所在处的褶皱中和面深度,获取钻孔A地应力测试信息以及成像测井钻井诱导缝与井壁崩落的解译;利用公式:计算褶皱中和面的位置,
式中:htop为中和面到卷入褶皱中地层最顶面的距离;H为卷入褶皱中地层的总厚度;n为层系粗化后总层数;hi为第i层的地层厚度;Ei为第i层地层的杨氏模量,υi为第i层地层的泊松比;
步骤3,根据钻孔A地应力测试、成像测井钻井诱导缝与井壁崩落的解译,获取褶皱中和面深度处现今地应力在x轴向和y轴向的分量大小,记为区域现今地应力Sr;
步骤4,收集研究区包含钻孔A的地震数据体,根据地震数据体分析获取钻孔A不同深度处的地层曲率1/rx和1/ry,1/rx为平行于x轴向的地层曲率,1/ry为平行于y轴向的地层曲率;
步骤5,利用公式:计算得到钻孔内距褶皱中和面距离为h处的应力分量大小σx和σy,记为局部现今地应力Sl,式中,h为钻孔内计算应力的深度点距褶皱中和面距离,σx为钻孔内距褶皱中和面距离为h处的x轴向应力,σy为钻孔内距褶皱中和面距离为h处的y轴向应力,Eh为钻孔内距褶皱中和面距离为h处的岩石杨氏模量,vh为钻孔内距褶皱中和面距离为h处的泊松比;
步骤6,利用公式:S=Sr+Sl,对区域现今地应力Sr和局部现今地应力Sl进行应力矢量合成,矢量合成的方向即为钻孔内距褶皱中和面距离为h处的现今地应力方位,由此,通过矢量合成距褶皱中和面不同距离h处的区域现今地应力和局部现今地应力,即可获取距褶皱中和面不同距离h处的现今地应力方位,实现定量分析现今地应力方向转向。h是钻孔内的某一深度点;h值不一样,步骤5算出的σx,和σy也就不一样,矢量合成时候的方位也就不一样,所以在不同的h就有不同的方向,也就实现了现今地应力方向转向的定量分析。
Claims (3)
1.一种褶皱发育区现今地应力方位转向的定量分析方法,其特征在于步骤如下:
步骤1,在发育褶皱的研究区设置多个垂直钻孔,收集研究区任一钻孔A的常规测井资料,利用声波时差测井曲线、密度测井曲线计算钻孔A不同深度处岩石的杨氏模量和泊松比,根据杨氏模量大小粗化层系,以杨氏模量(Emax-Emin)/10为间隔粗化层系,即Emin~Emin+(Emax-Emin)/10为一个层系,Emin+(Emax-Emin)/10~Emin+(Emax-Emin)/5为另一个层系,以此类推,其中,Emax与Emin分别为岩石杨氏模量的最大值和最小值;
步骤2,基于薄板模型,利用步骤1获取的杨氏模量、泊松比和粗化层系厚度,计算确定钻孔A所在处的褶皱中和面深度,获取钻孔A地应力测试信息以及成像测井钻井诱导缝与井壁崩落的解译;
步骤3,根据钻孔A地应力测试、成像测井钻井诱导缝与井壁崩落的解译,获取褶皱中和面深度处现今地应力在x轴向和y轴向的分量大小,记为区域现今地应力Sr;
步骤4,收集研究区包含钻孔A的地震数据体,根据地震数据体分析获取钻孔A不同深度处的地层曲率1/rx和1/ry,1/rx为平行于x轴向的地层曲率,1/ry为平行于y轴向的地层曲率;
步骤5,利用公式:计算得到钻孔内距褶皱中和面距离为h处的应力分量大小σx和σy,记为局部现今地应力Sl,式中,h为钻孔内计算应力的深度点距褶皱中和面距离,σx为钻孔内距褶皱中和面距离为h处的x轴向应力,σy为钻孔内距褶皱中和面距离为h处的y轴向应力,Eh为钻孔内距褶皱中和面距离为h处的岩石杨氏模量,vh为钻孔内距褶皱中和面距离为h处的泊松比;
步骤6,利用公式:S=Sr+Sl,对区域现今地应力Sr和局部现今地应力Sl进行应力矢量合成,矢量合成的方向即为钻孔内距褶皱中和面距离为h处的现今地应力方位,由此,通过矢量合成距褶皱中和面不同距离h处的区域现今地应力和局部现今地应力,即可获取距褶皱中和面不同距离h处的现今地应力方位,实现定量分析现今地应力方向转向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010375516.XA CN111520125B (zh) | 2020-05-07 | 2020-05-07 | 一种褶皱发育区现今地应力方位转向的定量分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010375516.XA CN111520125B (zh) | 2020-05-07 | 2020-05-07 | 一种褶皱发育区现今地应力方位转向的定量分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111520125A true CN111520125A (zh) | 2020-08-11 |
CN111520125B CN111520125B (zh) | 2021-11-12 |
Family
ID=71913123
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010375516.XA Active CN111520125B (zh) | 2020-05-07 | 2020-05-07 | 一种褶皱发育区现今地应力方位转向的定量分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111520125B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114442191A (zh) * | 2022-01-04 | 2022-05-06 | 成都理工大学 | 基于岩芯和井眼结构恢复确定井下现今地应力方向的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1966934A (zh) * | 2005-11-16 | 2007-05-23 | 中国石油大学(北京) | 一种随钻预测钻头底下地层坍塌压力和破裂压力的方法 |
CN105927216A (zh) * | 2016-04-25 | 2016-09-07 | 中国矿业大学 | 一种褶皱发育区储层裂缝预测方法及装置 |
CN108952700A (zh) * | 2018-08-21 | 2018-12-07 | 西南石油大学 | 一种各向异性地层井壁破裂压力确定方法 |
CN109113742A (zh) * | 2018-08-02 | 2019-01-01 | 中国矿业大学 | 一种煤储层现今地应力预测方法 |
CN109556782A (zh) * | 2019-01-29 | 2019-04-02 | 中国地质科学院地质力学研究所 | 一种基于钻孔套芯孔径变形的快捷地应力测试方法 |
CN110969699A (zh) * | 2018-09-29 | 2020-04-07 | 北京阳光杰科科技股份有限公司 | 一种高精度三维地应力建模方法及其应用 |
-
2020
- 2020-05-07 CN CN202010375516.XA patent/CN111520125B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1966934A (zh) * | 2005-11-16 | 2007-05-23 | 中国石油大学(北京) | 一种随钻预测钻头底下地层坍塌压力和破裂压力的方法 |
CN105927216A (zh) * | 2016-04-25 | 2016-09-07 | 中国矿业大学 | 一种褶皱发育区储层裂缝预测方法及装置 |
CN109113742A (zh) * | 2018-08-02 | 2019-01-01 | 中国矿业大学 | 一种煤储层现今地应力预测方法 |
CN108952700A (zh) * | 2018-08-21 | 2018-12-07 | 西南石油大学 | 一种各向异性地层井壁破裂压力确定方法 |
CN110969699A (zh) * | 2018-09-29 | 2020-04-07 | 北京阳光杰科科技股份有限公司 | 一种高精度三维地应力建模方法及其应用 |
CN109556782A (zh) * | 2019-01-29 | 2019-04-02 | 中国地质科学院地质力学研究所 | 一种基于钻孔套芯孔径变形的快捷地应力测试方法 |
Non-Patent Citations (1)
Title |
---|
鞠玮等: "滇东雨旺区块多煤层条件下现今地应力特征及其对煤层气开发的影响", 《煤炭学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114442191A (zh) * | 2022-01-04 | 2022-05-06 | 成都理工大学 | 基于岩芯和井眼结构恢复确定井下现今地应力方向的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111520125B (zh) | 2021-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105986816B (zh) | 一种用于识别页岩地层甜点的方法 | |
US20170145793A1 (en) | Method For Modeling Stimulated Reservoir Properties Resulting From Hydraulic Fracturing In Naturally Fractured Reservoirs | |
US20170051598A1 (en) | System For Hydraulic Fracturing Design And Optimization In Naturally Fractured Reservoirs | |
US10810331B2 (en) | System for predicting induced seismicity potential resulting from injection of fluids in naturally fractured reservoirs | |
CN108804728B (zh) | 水平井地层储层分级分析方法及计算机可读存储介质 | |
CN106154334A (zh) | 基于网格搜索的井下微地震事件实时反演定位方法 | |
CN105899749A (zh) | 偏斜井眼和水平井眼的地层和结构解释 | |
CN106875471A (zh) | 煤系含或隔水层三维可视化建模方法 | |
CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
CN111239815A (zh) | 基于三维地震属性的砂岩型铀储层成矿沉积要素提取方法 | |
CN105257252A (zh) | 利用测井资料优选页岩气水平井分簇射孔井段的方法 | |
CN103867198A (zh) | 碳酸盐岩天然气层地层密度判别法 | |
CN111520125B (zh) | 一种褶皱发育区现今地应力方位转向的定量分析方法 | |
CN113779811B (zh) | 一种构造变形对现今地应力方位扰动的定量分析方法 | |
CN108459358A (zh) | 沉积盆地基底顶面形态恢复和深部矿体定位预测的新方法 | |
CN104459790A (zh) | 含油气性盆地有效储层的分析方法和装置 | |
Lorenz et al. | Subsurface fracture spacing: Comparison of inferences from slant/horizontal and vertical cores | |
CN114856532A (zh) | 油气储层改造方法、装置和设备 | |
CN109212611A (zh) | 一种基于井约束匹配追踪的砂砾岩有效储层预测方法 | |
CN115563822A (zh) | 一种地应力各向异性多因素综合评价方法 | |
CN113514879B (zh) | “耳朵”层识别方法 | |
CN110794038B (zh) | 基于声波波速各向异性的岩心室内定向方法 | |
CN114740538A (zh) | 基于多元地球物理的矽卡岩型富铁矿深部勘查方法及*** | |
CN109374757A (zh) | 应用声波幅度定量化处理评价火成岩储层有效性的方法 | |
CN110596760A (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 |