CN113820750B - 基于弹塑性力学定量预测泥岩构造裂缝的方法 - Google Patents
基于弹塑性力学定量预测泥岩构造裂缝的方法 Download PDFInfo
- Publication number
- CN113820750B CN113820750B CN202110925418.3A CN202110925418A CN113820750B CN 113820750 B CN113820750 B CN 113820750B CN 202110925418 A CN202110925418 A CN 202110925418A CN 113820750 B CN113820750 B CN 113820750B
- Authority
- CN
- China
- Prior art keywords
- stress
- mudstone
- fracture
- rock
- crack
- 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 48
- 239000011435 rock Substances 0.000 claims abstract description 107
- 238000010276 construction Methods 0.000 claims abstract description 23
- 238000002474 experimental method Methods 0.000 claims abstract description 18
- 238000011084 recovery Methods 0.000 claims abstract description 15
- 239000002131 composite material Substances 0.000 claims abstract description 14
- 238000011161 development Methods 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims description 43
- 238000009826 distribution Methods 0.000 claims description 22
- 238000005553 drilling Methods 0.000 claims description 20
- 238000011160 research Methods 0.000 claims description 18
- 238000005336 cracking Methods 0.000 claims description 17
- 238000001125 extrusion Methods 0.000 claims description 16
- 230000006378 damage Effects 0.000 claims description 14
- 230000003068 static effect Effects 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000004088 simulation Methods 0.000 claims description 10
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 9
- 239000011707 mineral Substances 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000010008 shearing Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000008901 benefit Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 6
- 238000012669 compression test Methods 0.000 claims description 5
- 238000012512 characterization method Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 4
- 239000004927 clay Substances 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000013517 stratification Methods 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 238000004134 energy conservation Methods 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 238000010297 mechanical methods and process Methods 0.000 claims description 2
- 238000007493 shaping process Methods 0.000 claims description 2
- 239000004575 stone Substances 0.000 claims description 2
- 210000005056 cell body Anatomy 0.000 claims 2
- 210000004027 cell Anatomy 0.000 claims 1
- 239000003208 petroleum Substances 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 6
- 239000000463 material Substances 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 238000005452 bending Methods 0.000 description 2
- 238000009933 burial Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 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 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 241001415846 Procellariidae Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000013028 emission testing Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910052602 gypsum Inorganic materials 0.000 description 1
- 239000010440 gypsum Substances 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Abstract
本申请提供一种基于弹塑性力学定量预测泥岩构造裂缝的方法，涉及石油勘探开发技术领域。其中，该定量预测泥岩构造裂缝的方法包括：单井厚层泥岩段构造裂缝的分析和统计；进行岩石力学实验及力学参数测井解释；计算脆性指数及恢复古力学参数；泥岩屈服‑破坏复合强度准则建立；泥岩构造裂缝力学参数模型推导；地质模型构建及古应力状态恢复；三维古应力场模拟及裂缝定量预测。本申请技术方案解决了脆性指数约束下的泥岩屈服‑破裂复合强度准则建立的难题，适合于任何含泥岩隔层型碎屑岩储层的裂缝定量预测工作。
Description
技术领域
本申请涉及石油勘探开发技术领域，尤其涉及一种基于弹塑性力学定量预测泥岩构造裂缝的方法。
背景技术
随着油气勘探开发的不断进行与开采技术的持续进步，复杂断块油气藏、岩性油气藏、非常规致密砂岩油气藏以及裂缝性油气藏受到人们越来越多的关注；由于我国含油气盆地以陆相沉积地层为主，频繁的湖平面变化，造成了垂向上砂岩、泥岩的频繁交互变化，在多期构造运动作用下，不同厚度的泥岩段内也发育了一系列不同力学性质的裂缝，统计发现主要以低角度剪切缝为主，且裂缝内极少充填，两壁常出现光滑的摩擦镜面。
目前，针对裂缝的大量研究主要集中储层段内，逐渐形成了较为成熟的地质、地球物理和力学识别和预测方法，但对泥岩和页岩裂缝的研究仍停留在单井描述和平面预测阶段，相对于砂岩储层，泥岩由于掺杂了灰质、硅质、石膏等矿物成分，在受力变形过程中表现出不同程度的复杂塑性力学行为，给井距大、取芯少、泥岩厚地区的裂缝预测工作带来了很大难题。
发明内容
本申请实施例的目的是提供一种基于弹塑性力学定量预测泥岩构造裂缝的方法，以实现对泥岩和页岩裂缝的预测研究。
为解决上述技术问题，本申请实施例提供如下技术方案：
本申请提供一种基于弹塑性力学定量预测泥岩构造裂缝的方法，该基于弹塑性力学定量预测泥岩构造裂缝的方法包括：单井厚层泥岩段构造裂缝的分析和统计；进行岩石力学实验及力学参数测井解释；计算脆性指数及恢复古力学参数；泥岩屈服-破坏复合强度准则建立；泥岩构造裂缝力学参数模型推导；地质模型构建及古应力状态恢复；三维古应力场模拟及裂缝定量预测。
在本申请的一些变更实施方式中，所述单井厚层泥岩段构造裂缝的分析和统计，具体为：
收集现场成像测井资料和钻井取芯资料，统计泥岩段地层产状和构造裂缝参数；
绘制单井裂缝发育柱状图，明确厚层泥岩段裂缝分布的差异性，针对不同裂缝发育程度的深度点采样进行镜下薄片观察，识别矿物类型及含量，分析其对裂缝密度的影响；
其中，所述地层产状包括层理倾角和裂缝产状，所述裂缝参数包括裂缝线密度参数、裂缝体密度参数、裂缝开度参数及裂缝长度参数。
在本申请的一些变更实施方式中，所述进行岩石力学实验，具体为：
在泥岩厚层段，钻取柱塞岩样，且定义岩心钻取方向与轴向垂直，进行巴西劈裂试验、三轴压缩试验和直剪试验，测试并计算泥岩的Ⅰ型和Ⅱ型断裂韧性，同时记录各类试验的应力-应变曲线；
其中，所述三轴压缩试验用于获取岩石力学参数，所述岩石力学参数包括密度、弹性模量、泊松比、抗压强度、屈服强度、内摩擦系数及内摩擦角参数，所述巴西劈裂试验用于获得岩石的抗压强度，所述直剪试验用于获取岩石的抗剪强度、内摩擦系数以及内摩擦角；
在具有层理面的泥岩段，钻取柱塞岩样，且定义岩心钻取方向与轴向平行，进行巴西劈裂试验和直剪试验，同时记录各类试验的应力-应变曲线；
其中，所述巴西劈裂试验用于获得所述层理面的抗拉强度，所述直剪试验用于获得所述层理面的抗剪强度、内摩擦系数和内摩擦角。
在本申请的一些变更实施方式中，所述力学参数测井解释，具体为：
通过全波列测井和偶极横波测井计算得到横波时差曲线，并与常规测井资料进行相关性分析，拟合泥岩段的横波时差预测公式(1)，并通过公式(2)、(3)、(4)、(5)、(6)计算单井动态岩石力学参数；
Δts＝0.948Δtp+57.42 (1)
式中：Δts为横波时差(μs/m)，Δtp为纵波时差(μs/m)；
Sc＝Ed0.008Vsh+0.0045(1-Vsh) (4)
式中：Ed为动态杨氏模量(GPa)，ρb为体积密度(g/cm3)，μd为动态泊松比，Sc为单轴抗压强度(MPa)，C为岩石内聚力(MPa)，为层理面的内摩擦角(°)，Vsh为泥质的体积含量；
基于所述单井动态岩石力学参数结果进行断裂韧性计算，得出针对泥岩的Ⅰ型和Ⅱ型断裂韧性计算模型(7)(8)；
KⅡC＝0.0466P+0.1674St-0.1851 (8)
式中：P为围压(MPa)，St为互层岩石的抗拉强度(MPa)。
在本申请的一些变更实施方式中，所述计算脆性指数及恢复古力学参数，具体为：
通过公式(9)(10)对所述单井动态岩石力学模型进行连续矫正，制作单井静态力学参数图；
ES＝2.902Ed-14.92 (9)
μS＝4.59μd-1.07 (10)
式中：ES为静态杨氏模量(GPa)，μS为静态泊松比；
结合密度测井资料、阵列声波或偶极声波所测量的纵、横波数据和杨氏模量、泊松比参数，优选泊-杨法根据公式(11)计算泥岩段的脆性指数；
式中：BIe为弹性参数脆性指数，Emax和Emin分别为某一层段杨氏模量的最大和最小值，μmax和μmin分别为某一层段泊松比的最大和最小值，E为杨氏模量GPa；
收集地震叠前反演得到的研究区的P波阻抗、S波阻抗和密度数据体，并以单井测井解释结果为约束，计算现今应力条件下的弹性模量、泊松比，并运算得到泥岩段的脆性指数体；
收集所述研究区的深度域地震解释剖面，所述剖面垂直于主体构造走向，且所述剖面包含目的层和断裂解释结果，采用构造平衡剖面法恢复断距和二维褶皱形态，确定构造强烈活动期或关键造缝期；
建立所述研究区的泥岩层段构造地质模型，在划分断层相对级别的前提下，以一级、二级断层为控制，采用断层相关褶皱理论进行关键造缝期的古构造形态恢复；
通过三轴力学实验拟合建立不同围压条件下的力学参数变化曲线，将围压等效为埋藏深度下的上覆地层压力，根据所述古构造形态恢复的结果得到关键造缝期的泥岩层段古埋深，结合力学参数变化曲线计算该时期的古力学参数和脆性指数三维体，并按照脆性指数的分布趋势和区间特征将其进行等级划分。
在本申请的一些变更实施方式中，所述泥岩屈服-破坏复合强度准则建立，具体为：
判断三向应力状态下层理面的滑移和破坏情况；
判断三向压应力状态下层理面的滑移和破坏情况；
判断基岩的屈服和破裂状态，在三向挤压应力状态下，采用弹塑性屈服准则和塑形硬化准则判断岩石的屈服临界状态和破裂状态；
判断基岩中的裂缝倾角。
在本申请的一些变更实施方式中，所述判断三向应力状态下层理面的滑移和破坏情况，具体为：
当最小主应力σ3<0或σ3为张应力时且与地层的夹角为45°-90°时，采用公式(12)判断层理面是否发生滑移和破坏情况；
式中：Ft为岩石发生张破裂的应力条件，FS为岩石发生剪破裂的应力条件，σt为层理面的抗张强度(MPa)，μ为岩石的泊松比，P为围压(MPa)；
所述判断三向压应力状态下层理面的滑移和破坏情况，具体为：
当水平挤压力σ1与地层的夹角为30°-60°时，采用公式(13)判断层理面是否发生滑移和破坏；
式中：cj为层理面的内聚力(MPa)，为层理面的内摩擦角(°)，θh为水平挤压力σ1与地层的夹角(°)；
所述采用弹塑性屈服准则和塑形硬化准则判断岩石的屈服临界状态和破裂状态，具体为：
采用公式(14)判断岩石的屈服临界状态和破裂状态；
式中：I1为应力第一不变量，J2为应力第二不变量，a和K为与岩石内摩擦角和内聚力相关的实验常数，为岩石的内摩擦角(°)。
在本申请的一些变更实施方式中，所述泥岩构造裂缝力学参数模型推导，具体为：
选取岩石内表征单元体进行分析并简化，得到公式(15)；
式中：w总为当前应力状态下所述单元体受力吸收的总应变能密度(J/m3)，wd为岩石屈服后发生塑性变形所消耗的塑性应变能密度(J/m3)，w为存储的可释放弹性应变能密度(J/m3)，we为想产生裂缝必须克服的弹性应变能密度(J/m3)，V表征所述单元体的体积(m3)，Sf为新增裂缝的表面积(m2)，J为产生单位面积裂缝所需要的能量(J/m2)；
根据能量守恒原理和最大应变能密度理论，得到公式(16)；
式中：Dvf为所述单元体内缝隙的体密度(m2/m3)，E为性模量(GPa)，μ为泊松比；
推导建立三向挤压应力状态下泥岩和层理面裂缝密度、开度的力学模型，采用拉-剪复合破裂、莫尔－库仑破裂准则和德鲁克-普拉格准则进行岩石破裂判别，若满足破裂条件，通过公式(17)进行裂缝参数计算；
式中：J0为零围压时裂缝表面能(J/m2)，ε0为最大弹性张应变，σ3为最小张应变，K0为比例系数，w为应变能密度(J/m3)，wf为裂缝应变能密度(J/m3)，we为新增单位裂缝表面积必须克服的应变能(J/m3)，σp为岩石破裂应力(MPa)，Dlf为所述单元体内缝隙的线密度(l/m)，b为裂缝开度m，E为杨氏模量GPa，为岩石的内摩擦角(°)，θ为破裂角(°)，C0为岩石的内聚力MPa，L1、L3为最大主应力、最小主应力方向的单元网络长度m，b为裂缝开度m；
其中，所述最大弹性张应变为对应于0.85σc时的应变，所述最小张应变为对应0.45σc时的张应变。
在本申请的一些变更实施方式中，所述地质模型构建及古应力状态恢复，具体为：
收集地震解释、地质分层和测井解释结果，构建所述研究区的断层模型以及储层地质模型，以进行面模型重构，并通过面模型重构进行三维模型构建，完成三维地质力学模型的建立；
对所述研究区的钻井岩心进行声发射试验取样，每组样品至少包括4个试样；
将加工好的所述试样进行重复加载实验，加载速率为0.1MPa/s，测出所述试样在受载过程中的声发射信号随载荷的变化曲线，在第二次加载声发射载荷变化曲线上，找出声发射抹录不尽点，参照所述抹录不尽点的载荷值，在第一次加载声发射载荷曲线中，确定出凯塞尔点，取所述凯塞尔点和抹录不尽点载荷的平均值，即为岩心在地下所受的最大正应力，所述凯塞尔点的数量代表构造运动的期次；
根据构造形迹力学方法，进行裂缝的分期配套，进行裂缝走向统计分析，找出共轭裂缝的优势组系，以其共轭夹角平分线作为关键造缝期古应力场最大主应力方向；
其中，在所述声发射试验取样中，当试样的数量为4个时，其中一个所述试样取自品行井眼轴线的方向，其余所述试样取自水平面内互成45°角的三个方向。
在本申请的一些变更实施方式中，所述三维古应力场模拟及裂缝定量预测，具体为：
基于所述三维地质力学模型和关键造缝期古应力测试结果，设置模型力学边界，设定加载速率为0.1MPa/s，进行弹塑性有限元三维应力场模拟，并反复尝试，直到古应力场的最大主应力分量与次级挤压断层分布趋势、最小主应力分量与次级伸展断层分布一致、剪应力分量与走滑断层分布趋势一致为止；
基于得到的破裂准则、屈服准则、裂缝密度、开度和产状力学模型通过软件获取每个节点的最大主应力、最下主应力、中间主应力、剪应力、最大主应变、最下主应变、中间主应变、屈服强度参数，并经过计算得到每个节点裂缝的体密度、线密度和开度值；
以基岩力学层为单位进行单井预测裂缝线密度、开度与实测值对比，当平均吻合率低于90％时，修正裂缝力学参数、调整边界条件重新模拟三维应力场，直至平均吻合率达到90％以上为止，完成厚层泥岩段构造裂缝的三维定量预测。
相较于现有技术，本申请提供的基于弹塑性力学定量预测泥岩构造裂缝的方法，该方法解决了脆性指数约束下的泥岩屈服-破裂复合强度准则建立的难题，适合于任何含泥岩隔层型碎屑岩储层的裂缝定量预测工作；将任何软、硬岩石中的裂缝产生都归结为可释放弹性应变能的释放，推导建立了泥岩裂缝参数的力学模型，依托有限元模拟平台，可有效预测碎屑岩储层中构造裂缝的空间分布特征，包括裂缝密度、开度和产状，为泥岩裂缝性油气藏的勘探开发提供了地质依据，为避免钻井过程中井漏和井壁失稳事故的发生提供了保障，减少了勘探开发的风险和成本。
附图说明
通过参考附图阅读下文的详细描述，本申请示例性实施方式的上述以及其他目的、特征和优点将变得易于理解。在附图中，以示例性而非限制性的方式示出了本申请的若干实施方式，相同或对应的标号表示相同或对应的部分，其中：
图1示意性地示出了本发明实施例提供的种基于弹塑性力学定量预测泥岩构造裂缝的方法的流程示意图；
图2示意性地示出了本发明实施例提供的种基于弹塑性力学定量预测泥岩构造裂缝的方法的泥岩p-q塑性屈服曲线中参数对初始屈服面的影响示意图；
图3示意性地示出了本发明实施例提供的种基于弹塑性力学定量预测泥岩构造裂缝的方法的泥岩内变量等效塑性剪应变和硬化函数的关系示意图；
图4示意性地示出了本发明实施例提供的种基于弹塑性力学定量预测泥岩构造裂缝的方法的单位体积泥岩样品中能量的分配关系示意图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施方式。虽然附图中显示了本公开的示例性实施方式，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。
需要注意的是，除非另有说明，本申请使用的技术术语或者科学术语应当为本申请所属领域技术人员所理解的通常意义。
实施例一
参考附图1-附图4所示，本发明的实施例一提出一种基于弹塑性力学定量预测泥岩构造裂缝的方法，该基于弹塑性力学定量预测泥岩构造裂缝的方法包括：
步骤1：单井厚层泥岩段构造裂缝的分析和统计；
步骤2：进行岩石力学实验及力学参数测井解释；
步骤3：计算脆性指数及恢复古力学参数；
步骤4：泥岩屈服-破坏复合强度准则建立；
步骤5：泥岩构造裂缝力学参数模型推导；
步骤6：地质模型构建及古应力状态恢复；
步骤7：三维古应力场模拟及裂缝定量预测。
具体的，本方案针对塑性变形复杂的厚层泥岩段，采用岩心观察和测井解释进行单井构造裂缝的统计分析，采用岩石力学实验和弹塑性力学理论推导建立泥岩屈服破裂准则以及裂缝参数表征模型，采用地震解释和地层对比构建三维地质模型，在关键造缝期应力状态恢复的前提下，定量预测三维古应力场分布及裂缝空间分布。
进一步的，步骤1：单井厚层泥岩段构造裂缝的分析和统计，具体包括：
步骤1.1：由于成岩缝多数被矿物充填，收集现场成像测井(FMI)资料和钻井取芯资料，统计泥岩段地层产状和构造裂缝参数，包括层理倾角和裂缝产状、线密度、体密度、开度、长度参数；
步骤1.2：基于观察、统计结果，绘制单井裂缝发育柱状图，明确厚层泥岩段裂缝分布的差异性，针对不同裂缝发育程度的深度点采样进行镜下薄片观察，识别矿物类型及含量，分析其对裂缝密度的影响。
进一步的，步骤2：进行岩石力学实验及力学参数测井解释，具体包括：
步骤2.1：针对泥岩厚层段，钻取柱塞样，定义岩芯钻取方向与轴向垂直，参照岩石力学试验规范，岩样加工成尺寸设定为5cm×2.5cm×2.5cm的圆柱形岩样，并开展巴西劈裂试验、单轴/三轴压缩实验及直剪试验；其中，三轴压缩实验所用的岩样采用相近深度、相近岩性的三块为一组，并设置不同的围压值(0MPa、5MPa、10MPa……60MPa)，获取岩石力学参数，包括密度、弹性模量、泊松比、抗压强度、屈服强度、内摩擦系数和内摩擦角；通过巴西劈裂试验测试获得岩石的抗拉强度，同时在泥岩内设置初始裂缝，长度为2a，裂缝与加载方向夹角为θ，测试并计算出泥岩的I型和II型断裂韧性，记录各类岩石力学实验的应力-应变曲线；
步骤2.2：针对具有层理面的泥岩段，定义岩芯钻取方向与岩芯轴向平行，开展巴西劈裂试验和直剪试验；其中，通过巴西劈裂试验测试获得层理面的抗拉强度，通过变围压下的直剪试验获得层理面的抗剪强度、内摩擦系数和内摩擦角，同时记录各类岩石力学实验应力-应变曲线；
步骤2.3：针对常规测井资料缺少横波时差曲线的特点，用全波列测井和偶极横波测井计算得到横波时差曲线，并与常规测井资料进行相关性分析，拟合泥岩段的横波时差预测公式，即公式(1)：
Δts＝0.948Δtp+57.42 (1)
式中：Δts为横波时差(μs/m)，Δtp为纵波时差(μs/m)；
进而，可计算单井动态岩石力学参数，具体为：通过公式(2)计算动态弹性模量、通过公式(3)计算动态泊松比、通过公式(4)计算单轴抗压强度、通过公式(5)计算岩石内聚力以及通过公式(6)计算层理面的内摩擦角，公式(2)-(6)分别为：
Sc＝Ed0.008Vsh+0.0045(1-Vsh) (4)
式中：Ed为动态杨氏模量(GPa)，ρb为体积密度(g/cm3)，μd为动态泊松比，Sc为单轴抗压强度(MPa)，C为岩石内聚力(MPa)，为层理面的内摩擦角(°)，Vsh为泥质的体积含量(小数)；
步骤2.4：基于测井计算力学参数结果进行断裂韧性计算，优选出针对泥岩的I型和II型断裂韧性计算模型，分别为下述公式(7)和(8)；
KⅡC＝0.0466P+0.1674St-0.1851 (8)
式中：P为围压(MPa)，St为互层岩石的抗拉强度(MPa)。
进一步的，步骤3：计算脆性指数及恢复古力学参数，具体包括：
步骤3.1：由测井曲线计算出的动态力学参数和实验得到的静态参数会存在一定误差，静态力学参数更符合实际工程需要，因此，需要优选拟合公式分析动静转化关系，泥岩静态弹性模量采用公式(9)，泥岩静态泊松比采用公式(10)，用这两个公式对测井计算的力学参数进行连续校正，并制作单井静态力学参数图，公式(9)和(10)分别为：
ES＝2.902Ed-14.92 (9)
μS＝4.59μd-1.07 (10)
式中：ES为静态杨氏模量(GPa)，μS为静态泊松比；
步骤3.2：结合密度测井资料、阵列声波或偶极声波所测量的纵、横波数据和杨氏模量、泊松比参数，优选泊-杨法根据公式(11)计算泥岩段的脆性指数，公式(11)为：
式中：BIe为弹性参数脆性指数，Emax和Emin分别为某一层段杨氏模量的最大和最小值，μmax和μmin分别为某一层段泊松比的最大和最小值，E为杨氏模量GPa。
步骤3.3：收集地震叠前反演得到的研究区P波阻抗、S波阻抗和密度数据体，在此基础上以单井测井解释结果为约束，计算现今应力条件下的弹性模量、泊松比，并运算得到泥岩段的脆性指数体；
步骤3.4：收集研究区深度域地震解释剖面，剖面尽量垂直于主体构造走向，剖面包含目的层和断裂解释结果，采用构造平衡剖面法恢复断距和二维褶皱形态，确定构造强烈活动期或关键造缝期；
步骤3.5：建立研究区泥岩层段构造地质模型，在划分断层相对级别的前提下，以一级、二级断层为控制，采用断层相关褶皱理论进行关键造缝期古构造形态恢复，具体地，采用弯曲去褶皱法、恢复到基准面法恢复伸展构造区古构造，采用“拼板恢复”法恢复走滑构造区古构造，采用斜剪切法、弯曲滑动法和断层平行流法恢复挤压区古构造；
步骤3.6：通过三轴力学实验拟合建立不同围压条件下的力学参数变化曲线，将围压等效为埋藏深度下的上覆地层压力，根据古构造恢复结果可以得到关键造缝期的泥岩层段古埋深，结合力学参数变化曲线计算出此时期的古力学参数和脆性指数三维体，进一步按照脆性指数的分布趋势和区间特征，将其划分为三个等级：强(脆性)、中等(脆-塑性)和弱(塑性)。
进一步的，步骤4：泥岩屈服-破坏复合强度准则建立，具体包括：
根据大量岩心观察结果，当泥岩段富含碳酸盐岩、硅质矿物时，即脆性指数较大时，常显示出明显的层理构造，将泥岩段的变形和破坏方式分为层理面滑移型和基质破坏型两种基本类型；
步骤4.1：首先判断三向应力状态下层理面的滑移、破坏情况，当最小主应力σ3<0或张应力时与地层的夹角为45°-90°时，泥岩位于层面拉张破裂或剪切破裂可能范围内，但油气田地下的实际地层倾角极少能超过75°，根据平行四边形法则，泥岩层理面上既有张应力也有剪应力，更倾向于产生拉-剪复合破裂，采用公式(12)判断层理面是否发生滑移和破坏情况；
式中：Ft为岩石发生张破裂的应力条件，FS为岩石发生剪破裂的应力条件，σt为层理面的抗张强度(MPa)，μ为岩石的泊松比，P为围压(MPa)，即σ3；
步骤4.2：判断三向压应力状态下层理面的滑移、破坏情况，当水平挤压力σ1与地层的夹角为30°-60°时，泥岩段位于层面剪切破裂可能范围内，此时采用莫尔－库仑剪切准则判别层理面是否发生滑移破坏，采用公式(13)判断层理面是否发生滑移和破坏；
式中：cj为层理面的内聚力(MPa)，为层理面的内摩擦角(°)，θh为水平挤压力σ1与地层的夹角(°)；
步骤4.3：判断基岩，即泥岩的屈服、破裂状态，三向挤压应力状态下，采用弹塑性屈服准则(德鲁克-普拉格准则)和塑性硬化准则判别岩石的屈服临界状态和破裂状态，德鲁克-普拉格屈服准则认为当材料的最大形状变形能(又称畸变能)达到某一常数时，材料就屈服并进入塑性状态，采用公式(14)判断岩石的屈服临界状态和破裂状态：
式中：I1为应力第一不变量，J2为应力第二不变量，a和K为与岩石内摩擦角和内聚力相关的实验常数， 为岩石的内摩擦角(°)；
步骤4.4：根据步骤4.3，当泥岩发生屈服后，仍然有可能发生破裂，依据脆性指数等级分两种情况进行判断，当泥岩脆性指数位于弱等级时，符合理想弹塑性本构材料，进入初始无限塑性状态即预示岩石很难发生规模性破坏，主要以发生持续塑性变形为主，而当泥岩脆性指数位于强或中等等级时，符合硬化弹塑性本构材料，破坏面是硬化面和软化面的分界面，或者说破坏是塑性变形过程发展的最终结果，这样，可以定义岩石峰值应力p-q塑性屈服曲线，参考附图2所示，根据公式(18)确定出适合泥岩的屈服-破坏准则；
式中：αp为硬化函数，即塑性变形破坏过程中屈服面将随着硬化函数演变，C0和A分别是岩石内聚力和内摩擦系数，通过以上力学实验获得，其中， 为岩石的内摩擦角(°)，这里为保证参数A无量纲，设置P0＝1MPa，又/>σ1p为不同围压σ3下的泥岩峰值应力与前面砂岩的σc相当，σij为加载过程中的应力偏量(MPa)，δij为克罗内克符号；
步骤4.5：对于泥岩屈服-破坏准则中的硬化函数αp，其公式(19)如下：
式中：为硬化函数初始值，可由初始屈服点予以确定，参考附图2所示，当αp＝1表示岩石达到了宏观破坏时的硬化状态，故其取值范围为/>b为硬化参数，表示岩石的塑性硬化率，可基于硬化函数αp和塑性内变量γp的关系来确定，参考附图3所示，b值越小，αp变化越快，且逐渐趋近于1，/>为塑性应变增量偏张量，/>为塑性总应变量，εp为塑性体积塑性应变量，γp为等效剪切塑性应变量，这些应变和应力参数可在Ansys弹塑性应力场模拟过程中直接提取；
步骤4.6：判断基质(泥岩)中裂缝的倾角，在三维空间中，裂缝的倾角和走向要采取投影计算的方法来确定，Ansys坐标系的X轴与大地坐标的X轴(东)重合，Z轴与Y轴负方向(南)重合，Y和Z轴重合，因此若裂缝面法线方向向量在整体坐标系下的方向余玄已经确定为将/>投影到XOZ平台，其投影线与Z轴负方向的夹角为αz＝arctan(-l/n)，走向-角α可确定：若0≤αz<90°，α＝90°-αz，若-90°<αz<0，α＝(-90°-αz)+360°，在Ansys三维坐标系中，裂缝倾角是裂缝表面与XY平面的夹角，也就是平面lx+my+nz＝0与平面y＝0之间的夹角αdip(0°≤αdip≤90°)，其公式(20)如下：/>
进一步的，步骤5：泥岩构造裂缝力学参数模型推导，具体包括：
步骤5.1：大量岩石力学实验结果表明：可释放弹性应变能是促使岩石破坏的根本原因，泥岩内主应力达到单轴抗压强度σc的0.85倍时，大量微裂缝开始集结并产生宏观裂缝，因此可以把该时刻的应变能称作要产生裂缝而必须克服的最大弹性应变能，于是新增裂缝表面积需要的能量就等于当前单元体内应变能与要产生裂缝而必须克服的最大弹性应变能之差，选取岩石内表征单元体进行分析(REV)，在适当简化的前提下，参考附图4所示，得到公式(15)；
式中：w总为当前应力状态下所述单元体受力吸收的总应变能密度(J/m3)，wd为岩石屈服后发生塑性变形所消耗的塑性应变能密度(J/m3)，w为存储的可释放弹性应变能密度(J/m3)，we为想产生裂缝必须克服的弹性应变能密度(J/m3)，V表征所述单元体的体积(m3)，Sf为新增裂缝的表面积(m2)，J为产生单位面积裂缝所需要的能量(J/m2)；
步骤5.2：由能量守恒原理和最大应变能密度理论，泥岩新增裂缝总表面积等于裂缝体密度，于是有公式(16)；
式中：Dvf为所述单元体内缝隙的体密度(m2/m3)，即裂缝体的总表面积与单元体体积之比，E为杨氏模量GPa，μ为泊松比；
步骤5.3：推导建立三向挤压应力状态下泥岩和层理面裂缝密度、开度的力学模型，采用拉-剪复合破裂、莫尔－库仑破裂准则和德鲁克-普拉格准则进行岩石破裂判别，若满足破裂条件，用公式(17)进行裂缝参数计算；
/>
式中：J0为零围压时裂缝表面能(J/m2)，ε0为最大弹性张应变，对应于0.85σc时的应变，σ3为最小张应变，对应0.45σc时的张应变，K0为比例系数，w为应变能密度(J/m3)，wf为裂缝应变能密度(J/m3)，we为新增单位裂缝表面积必须克服的应变能(J/m3)，σp为岩石破裂应力(MPa)，Dlf为所述单元体内缝隙的线密度(l/m)，b为裂缝开度m，为内摩擦角(°)，θ为破裂角(°)，其中，破裂压力可由公式(17)或公式(22)换算得到，对应0.85σc这一关键点。
进一步的，步骤6：地质模型构建及古应力状态恢复，具体包括：
步骤6.1：以Petrel地质建模软件为平台，收集地震解释、地质分层和测井解释结果，构建研究区断层模型以及储层地质模型，将层面和断层数据导入AutoCad软件进行面模型重构，并将面模型再次导出到Ansys软件中进行三维模型构建，划分网格后按照“形心”搜索法将井-震反演获得的古力学参数离散赋值到三维模型的每个网格内，完成三维地质力学模型的建立；
步骤6.2：对研究区钻井岩心进行声发射试验取样，钻井位置最好平均分布在研究区平面上或呈四边形，每组样品至少包括4个试样，其中1个取自垂直方向(平行井眼轴线)，另3个取自水平面内互成45°角的3个方向；
步骤6.3：声发射试验所用设备由伺服岩石刚性试验机和声发射测试系统组成，将加工好的试样在室内进行重复加载实验，加载速率为0.1MPa/s，测出岩样受载过程中的声发射信号随载荷的变化曲线，在第二次加载声发射载荷变化曲线上，找出声发射抹录不净点，参照抹录不尽点载荷值，在第一次加载声发射载荷曲线中，确定出凯塞尔点，取凯塞尔点和抹录不尽点载荷的平均值，即为岩心在地下所受的最大正应力，凯塞尔点的数量即代表了构造运动的期次；
步骤6.4：根据构造形迹力学方法，进行裂缝的分期配套，采用玫瑰花图进行裂缝走向统计分析，找出共轭裂缝的优势组系，以其共轭夹角平分线作为关键造缝期古应力场最大主应力方向。
进一步的，步骤7：三维古应力场模拟及裂缝定量预测，具体包括：
步骤7.1：基于Ansys中建立的三维地质力学模型和关键造缝期古应力测试结果，设置模型力学边界，设定加载速率为0.1MPa/s，进行弹塑性有限元三维应力场模拟，并反复尝试，直到古应力场的最大主应力分量与次级挤压断层分布趋势、最小主应力分量与次级伸展断层分布一致、剪应力分量与走滑断层分布趋势一致为止：
步骤7.2：将上述破裂准则、屈服准则、裂缝密度、开度和产状力学模型编写成APDL语言程序，植入到Ansys软件平台中，在古应力场模拟的基础上，提取每个节点的最大主应力、最下主应力、中间主应力、剪应力、最大主应变、最下主应变、中间主应变、屈服强度参数，经过计算后可得到每个节点裂缝的体密度、线密度和开度值；
步骤7.3：以岩石力学层为单位进行单井预测裂缝线密度、开度与实测值对比，如果平均吻合率低于90％，则修正裂缝力学参数、调整边界条件重新模拟三维应力场，直到平均吻合率达到90％以上为止，最终完成厚层泥岩段构造裂缝的三维定量预测。
根据上述所列，本发明实施例提出一种基于弹塑性力学定量预测泥岩构造裂缝的方法，能够解决脆性指数约束下的泥岩屈服-破裂复合强度准则建立的难题，适合于任何含泥岩隔层型碎屑岩储层的裂缝定量预测工作；将任何软、硬岩石中的裂缝产生都归结为可释放弹性应变能的释放，推导建立了泥岩裂缝参数的力学模型，依托有限元模拟平台，可有效预测碎屑岩储层中构造裂缝的空间分布特征，包括裂缝密度、开度和产状，为泥岩裂缝性油气藏的勘探开发提供了地质依据，为避免钻井过程中井漏和井壁失稳事故的发生提供了保障，减少了勘探开发的风险和成本。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。
Claims (8)
1.一种基于弹塑性力学定量预测泥岩构造裂缝的方法，其特征在于，包括：
单井厚层泥岩段构造裂缝的分析和统计；
进行岩石力学实验及力学参数测井解释；
计算脆性指数及恢复古力学参数；
泥岩屈服-破坏复合强度准则建立；
泥岩构造裂缝力学参数模型推导；
地质模型构建及古应力状态恢复；
三维古应力场模拟及裂缝定量预测；
其中，所述泥岩屈服-破坏复合强度准则建立，具体为：
判断三向应力状态下层理面的滑移和破坏情况；
判断三向压应力状态下层理面的滑移和破坏情况；
判断基岩的屈服和破裂状态，在三向挤压应力状态下，采用弹塑性屈服准则和塑形硬化准则判断岩石的屈服临界状态和破裂状态；
判断基岩中的裂缝倾角；
其中，所述判断三向应力状态下层理面的滑移和破坏情况，具体为：
当最小主应力σ3<0或σ3为张应力时且与地层的夹角为45°-90°时，采用公式(12)判断层理面是否发生滑移和破坏情况；
式中：Ft为岩石发生张破裂的应力条件，FS为岩石发生剪破裂的应力条件，σt为层理面的抗张强度MPa，μ为岩石的泊松比，P为围压MPa，为岩石的内摩擦角°，C为岩石内聚力MPa；
所述判断三向压应力状态下层理面的滑移和破坏情况，具体为：
当水平挤压力σ1与地层的夹角为30°-60°时，采用公式(13)判断层理面是否发生滑移和破坏；
式中：cj为层理面的内聚力MPa，为层理面的内摩擦角°，θh为水平挤压力σ1与地层的夹角°；
所述采用弹塑性屈服准则和塑形硬化准则判断岩石的屈服临界状态和破裂状态，具体为：
采用公式(14)判断岩石的屈服临界状态和破裂状态；
式中：f为岩石屈服临街状态和破裂状态的应力条件，I1为应力第一不变量，J2为应力第二不变量，a和K为与岩石内摩擦角和内聚力相关的实验常数，为岩石的内摩擦角°，σ2为中间主应力MPa，σx为x方向正应力MPa，σy为y方向正应力MPa，σz为z方向正应力MPa，τxy为平行于xy面方向的剪应力MPa，τyz为平行于yz面方向的剪应力MPa，τxz为平行于xz面方向的剪应力MPa。
2.根据权利要求1所述的基于弹塑性力学定量预测泥岩构造裂缝的方法，其特征在于，
所述单井厚层泥岩段构造裂缝的分析和统计，具体为：
收集现场成像测井资料和钻井取芯资料，统计泥岩段地层产状和构造裂缝参数；
绘制单井裂缝发育柱状图，明确厚层泥岩段裂缝分布的差异性，针对不同裂缝发育程度的深度点采样进行镜下薄片观察，识别矿物类型及含量，分析其对裂缝密度的影响；
其中，所述地层产状包括层理倾角和裂缝产状，所述裂缝参数包括裂缝线密度参数、裂缝体密度参数、裂缝开度参数及裂缝长度参数。
3.根据权利要求1所述的基于弹塑性力学定量预测泥岩构造裂缝的方法，其特征在于，
所述进行岩石力学实验，具体为：
在泥岩厚层段，钻取柱塞岩样，且定义岩心钻取方向与轴向垂直，进行巴西劈裂试验、三轴压缩试验和直剪试验，测试并计算泥岩的Ⅰ型和Ⅱ型断裂韧性，同时记录各类试验的应力-应变曲线；
其中，所述三轴压缩试验用于获取岩石力学参数，所述岩石力学参数包括密度、弹性模量、泊松比、抗压强度、屈服强度、内摩擦系数及内摩擦角参数，所述巴西劈裂试验用于获得岩石的抗压强度，所述直剪试验用于获取岩石的抗剪强度、内摩擦系数以及内摩擦角；
在具有层理面的泥岩段，钻取柱塞岩样，且定义岩心钻取方向与轴向平行，进行巴西劈裂试验和直剪试验，同时记录各类试验的应力-应变曲线；
其中，所述巴西劈裂试验用于获得所述层理面的抗拉强度，所述直剪试验用于获得所述层理面的抗剪强度、内摩擦系数和内摩擦角。
4.根据权利要求3所述的基于弹塑性力学定量预测泥岩构造裂缝的方法，其特征在于，
所述力学参数测井解释，具体为：
通过全波列测井和偶极横波测井计算得到横波时差曲线，并与常规测井资料进行相关性分析，拟合泥岩段的横波时差预测公式(1)，并通过公式(2)、(3)、(4)、(5)、(6)计算单井动态岩石力学参数；
Δts＝0.948Δtp+57.42 (1)
式中：Δts为横波时差μs/m，Δtp为纵波时差μs/m；
Sc＝Ed0.008Vsh+0.0045(1-Vsh) (4)
式中：Ed为动态杨氏模量GPa，ρb为体积密度g/cm3，μd为动态泊松比，Sc为单轴抗压强度MPa，C为岩石内聚力MPa，为层理面的内摩擦角°，Vsh为泥质的体积含量；
基于所述单井动态岩石力学参数结果进行断裂韧性计算，得出针对泥岩的Ⅰ型和Ⅱ型断裂韧性计算模型(7)(8)；
KⅡC＝0.0466P+0.1674St-0.1851 (8)
式中：P为围压MPa，St为互层岩石的抗拉强度MPa。
5.根据权利要求4所述的基于弹塑性力学定量预测泥岩构造裂缝的方法，其特征在于，
所述计算脆性指数及恢复古力学参数，具体为：
通过公式(9)(10)对单井动态岩石力学模型进行连续矫正，制作单井静态力学参数图；
ES＝2.902Ed-14.92 (9)
μS＝4.59μd-1.07 (10)
式中：ES为静态杨氏模量GPa，μS为静态泊松比；
结合密度测井资料、阵列声波或偶极声波所测量的纵、横波数据和杨氏模量、泊松比参数，选择泊-杨法根据公式(11)计算泥岩段的脆性指数；
式中：BIe为弹性参数脆性指数，Emax和Emin分别为某一层段杨氏模量的最大和最小值，μmax和μmin分别为某一层段泊松比的最大和最小值，E为杨氏模量GPa；
收集地震叠前反演得到的研究区的P波阻抗、S波阻抗和密度数据体，并以单井测井解释结果为约束，计算现今应力条件下的杨氏模量、泊松比，并运算得到泥岩段的脆性指数体；
收集所述研究区的深度域地震解释剖面，所述剖面垂直于主体构造走向，且所述剖面包含目的层和断裂解释结果，采用构造平衡剖面法恢复断距和二维褶皱形态，确定构造强烈活动期或关键造缝期；
建立所述研究区的泥岩层段构造地质模型，在划分断层相对级别的前提下，以一级、二级断层为控制，采用断层相关褶皱理论进行关键造缝期的古构造形态恢复；
通过三轴力学实验拟合建立不同围压条件下的力学参数变化曲线，将围压等效为埋藏深度下的上覆地层压力，根据所述古构造形态恢复的结果得到关键造缝期的泥岩层段古埋深，结合力学参数变化曲线计算该时期的古力学参数和脆性指数三维体，并按照脆性指数的分布趋势和区间特征将其进行等级划分。
6.根据权利要求5所述的基于弹塑性力学定量预测泥岩构造裂缝的方法，其特征在于，
所述泥岩构造裂缝力学参数模型推导，具体为：
选取岩石内表征单元体进行分析并简化，得到公式(15)；
式中：w总为当前应力状态下所述单元体受力吸收的总应变能密度J/m3，wd为岩石屈服后发生塑性变形所消耗的塑性应变能密度J/m3，w为存储的可释放弹性应变能密度J/m3，we为想产生裂缝必须克服的弹性应变能密度J/m3，V表征所述单元体的体积m3，Sf为新增裂缝的表面积m2，J为产生单位面积裂缝所需要的能量J/m2；
根据能量守恒原理和最大应变能密度理论，得到公式(16)；
式中：Dvf为所述单元体内缝隙的体密度m2/m3，μ为泊松比；
推导建立三向挤压应力状态下泥岩和层理面裂缝密度、开度的力学模型，采用拉-剪复合破裂、莫尔－库仑破裂准则和德鲁克-普拉格准则进行岩石破裂判别，若满足破裂条件，通过公式(17)进行裂缝参数计算；
式中：J0为零围压时裂缝表面能J/m2，ε0为最大弹性张应变，σ3为最小张应变，K0为比例系数，w为应变能密度J/m3，wf为裂缝应变能密度J/m3，we为新增单位裂缝表面积必须克服的应变能J/m3，σp为岩石破裂应力MPa，Dlf为所述单元体内缝隙的线密度J/m，为岩石的内摩擦角°，θ为破裂角°，C0为岩石的内聚力MPa，L1、L3为最大主应力、最小主应力方向的单元网络长度m，b为裂缝开度m；
其中，所述最大弹性张应变为对应于0.85σc时的应变，所述最小张应变为对应0.45σc时的张应变，σc为压缩强度极限MPa。
7.根据权利要求6所述的基于弹塑性力学定量预测泥岩构造裂缝的方法，其特征在于，
所述地质模型构建及古应力状态恢复，具体为：
收集地震解释、地质分层和测井解释结果，构建所述研究区的断层模型以及储层地质模型，以进行面模型重构，并通过面模型重构进行三维模型构建，完成三维地质力学模型的建立；
对所述研究区的钻井岩心进行声发射试验取样，每组样品至少包括4个试样；
将加工好的所述试样进行重复加载实验，加载速率为0.1MPa/s，测出所述试样在受载过程中的声发射信号随载荷的变化曲线，在第二次加载声发射载荷变化曲线上，找出声发射抹录不尽点，参照所述抹录不尽点的载荷值，在第一次加载声发射载荷曲线中，确定出凯塞尔点，取所述凯塞尔点和抹录不尽点载荷的平均值，即为岩心在地下所受的最大正应力，所述凯塞尔点的数量代表构造运动的期次；
根据构造形迹力学方法，进行裂缝的分期配套，进行裂缝走向统计分析，找出共轭裂缝的优势组系，以其共轭夹角平分线作为关键造缝期古应力场最大主应力方向；
其中，在所述声发射试验取样中，当试样的数量为4个时，其中一个所述试样取自品行井眼轴线的方向，其余所述试样取自水平面内互成45°角的三个方向。
8.根据权利要求7所述的基于弹塑性力学定量预测泥岩构造裂缝的方法，其特征在于，
所述三维古应力场模拟及裂缝定量预测，具体为：
基于所述三维地质力学模型和关键造缝期古应力测试结果，设置模型力学边界，设定加载速率为0.1MPa/s，进行弹塑性有限元三维应力场模拟，并反复尝试，直到古应力场的最大主应力分量与次级挤压断层分布趋势、最小主应力分量与次级伸展断层分布一致、剪应力分量与走滑断层分布趋势一致为止；
基于得到的破裂准则、屈服准则、裂缝密度、开度和产状力学模型通过软件获取每个节点的最大主应力、最下主应力、中间主应力、剪应力、最大主应变、最下主应变、中间主应变、屈服强度参数，并经过计算得到每个节点裂缝的体密度、线密度和开度值；
以基岩力学层为单位进行单井预测裂缝线密度、开度与实测值对比，当平均吻合率低于90％时，修正裂缝力学参数、调整边界条件重新模拟三维应力场，直至平均吻合率达到90％以上为止，完成厚层泥岩段构造裂缝的三维定量预测。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110925418.3A CN113820750B (zh) | 2021-08-12 | 2021-08-12 | 基于弹塑性力学定量预测泥岩构造裂缝的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110925418.3A CN113820750B (zh) | 2021-08-12 | 2021-08-12 | 基于弹塑性力学定量预测泥岩构造裂缝的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113820750A CN113820750A (zh) | 2021-12-21 |
| CN113820750B true CN113820750B (zh) | 2024-04-16 |
Family
ID=78913167
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110925418.3A Active CN113820750B (zh) | 2021-08-12 | 2021-08-12 | 基于弹塑性力学定量预测泥岩构造裂缝的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113820750B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115267905B (zh) * | 2022-07-28 | 2023-12-19 | 中国石油大学(华东) | 一种复杂构造区钻井工程中裂缝性井漏预测方法 |
| CN117235839B (zh) * | 2022-09-16 | 2024-04-19 | 中国葛洲坝集团三峡建设工程有限公司 | 基于自平衡应力场分析的井筒水泥环安全载荷计算方法 |
| CN115856996B (zh) * | 2022-11-08 | 2023-06-30 | 应急管理部国家自然灾害防治研究院 | 一种地震数值预测预报方法和系统 |
| CN115618658B (zh) * | 2022-12-16 | 2023-03-14 | 新疆石油管理局有限公司 | 基于弹塑性本构的地层破裂压力动态预测方法及装置 |
| CN116359013A (zh) * | 2023-03-31 | 2023-06-30 | 长江水利委员会长江科学院 | 基于多数据源分析的隧洞围岩力学参数取值方法 |
| CN116822330B (zh) * | 2023-05-19 | 2024-02-20 | 四川大学 | 支撑剂弹塑性嵌入过程分析方法、装置、设备及存储介质 |
| CN116337628B (zh) * | 2023-05-31 | 2024-03-22 | 中国地质大学（北京） | 三维应力下活动断裂粘滑发震模拟方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104200039A (zh) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | 一种构造裂缝产状定量预测方法 |
| CN105259035A (zh) * | 2015-10-26 | 2016-01-20 | 中国石油大学(华东) | 一种岩石材料时效-弹塑性力学本构模型的构建方法 |
| CN106019405A (zh) * | 2016-05-13 | 2016-10-12 | 赵向原 | 储层裂缝建模方法及系统 |
| CN107489417A (zh) * | 2017-07-25 | 2017-12-19 | 中国海洋石油总公司 | 基于xrf录井的岩矿组成反演法砂泥岩岩性识别方法 |
| CN112363226A (zh) * | 2020-11-06 | 2021-02-12 | 大庆油田有限责任公司 | 一种非常规油气有利区地球物理预测方法 |
| CN113534291A (zh) * | 2021-07-20 | 2021-10-22 | 中国石油大学(华东) | 岩石力学层约束下的低渗透储层不同尺度裂缝定量预测方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7649804B2 (en) * | 2007-10-25 | 2010-01-19 | Schlumberger Technology Corporation | In-situ determination of yield stress state of earth formations |
| US8498853B2 (en) * | 2009-07-20 | 2013-07-30 | Exxonmobil Upstream Research Company | Petrophysical method for predicting plastic mechanical properties in rock formations |
| US20150205006A1 (en) * | 2010-03-25 | 2015-07-23 | Schlumberger Technology Corporation | Downhole modeling using inverted pressure and regional stress |
| US10619469B2 (en) * | 2016-06-23 | 2020-04-14 | Saudi Arabian Oil Company | Hydraulic fracturing in kerogen-rich unconventional formations |
| US10816440B2 (en) * | 2017-02-20 | 2020-10-27 | Conocophillips Company | Rock mechanical properties from drill cuttings |
-
2021
- 2021-08-12 CN CN202110925418.3A patent/CN113820750B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104200039A (zh) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | 一种构造裂缝产状定量预测方法 |
| CN105259035A (zh) * | 2015-10-26 | 2016-01-20 | 中国石油大学(华东) | 一种岩石材料时效-弹塑性力学本构模型的构建方法 |
| CN106019405A (zh) * | 2016-05-13 | 2016-10-12 | 赵向原 | 储层裂缝建模方法及系统 |
| CN107489417A (zh) * | 2017-07-25 | 2017-12-19 | 中国海洋石油总公司 | 基于xrf录井的岩矿组成反演法砂泥岩岩性识别方法 |
| CN112363226A (zh) * | 2020-11-06 | 2021-02-12 | 大庆油田有限责任公司 | 一种非常规油气有利区地球物理预测方法 |
| CN113534291A (zh) * | 2021-07-20 | 2021-10-22 | 中国石油大学(华东) | 岩石力学层约束下的低渗透储层不同尺度裂缝定量预测方法 |
Non-Patent Citations (8)
| Title |
|---|
| A Mud Weight Window Investigation Based on the Mogi-Coulomb Failure Criterion and Elasto-Plastic Constitutive Mode;Mirzaghorbanali, A等;Petroleum Science and Technology;第29卷(第02期);第121-131页 * |
| Elasto-plastic fracturing model for wellbore stability using non-penetrating fluids;Bernt S. Aadnøy等;Journal of Petroleum Science and Engineering;第45卷(第3-4期);第179–192页 * |
| 基于试验的花岗岩渐进破坏本构模型研究;刘泉声;胡云华;刘滨;;岩土力学;第30卷(第02期);第289-296页 * |
| 岩石力学层与构造裂缝发育关系研究;赵乐强等;山东科技大学学报（自然科学版）;第37卷(第01期);第36-45页 * |
| 新场气田须五段岩石力学参数特征及测井解释;徐浩;谢润成;杨松;姚勇;张冲;周文;单钰铭;;科学技术与工程(第08期);第19-21页 * |
| 砂泥岩间互地层破裂准则选取及在裂缝穿透性评价中的应用;卢虎胜;张俊峰;李世银;邹军;刘爱武;;中国石油大学学报(自然科学版)(第03期);第14-19页 * |
| 赵乐强等.岩石力学层与构造裂缝发育关系研究.山东科技大学学报（自然科学版）.2018,第37卷(第01期),第36-45页. * |
| 非常规油气储层脆性指数测井评价方法研究进展;赖锦;王贵文;范卓颖;陈晶;王抒忱;周正龙;范旭强;;石油科学通报(第03期);第331-332页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113820750A (zh) | 2021-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113820750B (zh) | 基于弹塑性力学定量预测泥岩构造裂缝的方法 | |
| CN113534291B (zh) | 岩石力学层约束下的低渗透储层不同尺度裂缝定量预测方法 | |
| CN104200039A (zh) | 一种构造裂缝产状定量预测方法 | |
| CN115267905B (zh) | 一种复杂构造区钻井工程中裂缝性井漏预测方法 | |
| CN105653811B (zh) | 高压变电站深回填区旋挖钻进中的入岩深度确定方法 | |
| CN106372297A (zh) | 砂土岩溶地层中盾构与溶洞安全垂直距离的确定方法 | |
| CN114755310A (zh) | 一种裂缝性储层岩石力学层演变规律预测方法 | |
| Sobhaniaragh et al. | Computational modelling of multi-stage hydraulic fractures under stress shadowing and intersecting with pre-existing natural fractures | |
| CN104153768A (zh) | 一种评价花岗岩储层储集性能的方法 | |
| CN106199712B (zh) | 一种确定压裂套管变形区域的方法及装置 | |
| CN115903078A (zh) | 水力压裂层位及其确定方法、水力压裂位置的确定方法 | |
| CN117555029A (zh) | 一种油气储层地质力学层单井定量划分方法 | |
| Zhang | In situ stress inversion and distribution characteristics of tunnel based on numerical simulation and neural network technology | |
| Xia et al. | Study on shear strength characteristics of columnar jointed basalt based on in-situ direct shear test at Baihetan Hydropower Station | |
| Feng et al. | Development characteristics and quantitative prediction of multiperiod fractures in superdeep thrust-fold belt | |
| CN116305462A (zh) | 一种考虑爆破与富水扰动影响的岩石Hoek-Brown准则参数确定方法 | |
| Pierce et al. | Synthetic rock mass applications in mass mining | |
| Shahid et al. | Determining the rock brittle index (BI) using multivariate regression (a case study) | |
| Sugawara et al. | Rock Stress' 03: Proceedings of the Third International Symposium on Rock Stress, Kumamoto, Japan, 4-6 November 2003 | |
| Gao et al. | An Evaluation Method of Rock Brittleness Based on the Prepeak Crack Initiation and Postpeak Stress Drop Characteristics | |
| Xu et al. | Geomechanical characteristics and stimulation of Dibei deep tight sandstone reservoirs in the Kuqa DeprEssion of Tarim Basin | |
| Lian et al. | Simulation of ground stress field and fracture anticipation with effect of pore pressure | |
| Li et al. | Sensitivity analysis in the estimation of mechanical parameters of engineering rock mass based on the Hoek–Brown criterion | |
| CN112418597B (zh) | 一种适用于页岩气老区开发调整井的储层可压性评价方法 | |
| Lin et al. | Prediction of natural fracture distribution characteristics in tight oil and gas reservoirs based on paleotectonic stress field |
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 |