CN111460601A - 一种基于岩石物理建模的正交各向异性地层地应力预测方法 - Google Patents
一种基于岩石物理建模的正交各向异性地层地应力预测方法 Download PDFInfo
- Publication number
- CN111460601A CN111460601A CN201911239424.2A CN201911239424A CN111460601A CN 111460601 A CN111460601 A CN 111460601A CN 201911239424 A CN201911239424 A CN 201911239424A CN 111460601 A CN111460601 A CN 111460601A
- Authority
- CN
- China
- Prior art keywords
- gpa
- shale
- rock
- orthotropic
- model
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 34
- 239000011148 porous material Substances 0.000 claims abstract description 29
- 230000003068 static effect Effects 0.000 claims abstract description 18
- 238000006243 chemical reaction Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 59
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 55
- 239000011707 mineral Substances 0.000 claims description 55
- 239000000203 mixture Substances 0.000 claims description 46
- 239000012530 fluid Substances 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 16
- 239000000463 material Substances 0.000 claims description 16
- 238000012360 testing method Methods 0.000 claims description 16
- 239000005416 organic matter Substances 0.000 claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 11
- 229920006395 saturated elastomer Polymers 0.000 claims description 9
- 238000002474 experimental method Methods 0.000 claims description 7
- 239000002023 wood Substances 0.000 claims description 6
- 230000006835 compression Effects 0.000 claims description 5
- 238000007906 compression Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 4
- IAZDPXIOMUYVGZ-UHFFFAOYSA-N Dimethylsulphoxide Chemical compound CS(C)=O IAZDPXIOMUYVGZ-UHFFFAOYSA-N 0.000 claims description 3
- 230000001364 causal effect Effects 0.000 claims description 3
- 238000005056 compaction Methods 0.000 claims description 3
- 239000008398 formation water Substances 0.000 claims description 3
- 230000002706 hydrostatic effect Effects 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000005336 cracking Methods 0.000 abstract description 2
- 239000004927 clay Substances 0.000 description 8
- 239000002245 particle Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 229910021532 Calcite Inorganic materials 0.000 description 3
- 229910000514 dolomite Inorganic materials 0.000 description 3
- 239000010459 dolomite Substances 0.000 description 3
- 239000010433 feldspar Substances 0.000 description 3
- 239000011028 pyrite Substances 0.000 description 3
- 229910052683 pyrite Inorganic materials 0.000 description 3
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 3
- 239000010453 quartz Substances 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000013351 cheese Nutrition 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- 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
-
- 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/42—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators in one well and receivers elsewhere or vice versa
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/10—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices using induction coils
-
- 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
- E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
- E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Food Science & Technology (AREA)
- Geophysics (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Medicinal Chemistry (AREA)
- Mining & Mineral Resources (AREA)
- Electromagnetism (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- Acoustics & Sound (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于岩石物理建模的正交各向异性地层地应力预测方法，包括以下步骤：获取页岩物性参数，建立正交各向异性岩石物理模型，计算刚度系数、动静态岩石弹性参数、Biot系数，基于Eaton法利用纵波时差来预测地层孔隙压力，最后计算正交各向异性页岩地层最大、最小水平主应力。本发明构建了基于岩石物理建模的地应力评价体系，提高了静态岩石弹性参数的转换精度和地应力的预测精度，解决了正交各向异性页岩地层的地应力评价难题。利用该发明计算的正交各向异性地层地应力可以为后续的正交各向异性页岩地层井壁稳定性分析、压裂起裂分析等提供基础数据。
Description
技术领域
本发明涉及一种基于岩石物理建模的正交各向异性地层地应力预测方法，属于岩石物理技术领域。
背景技术
地应力是指存在于地壳中的应力，其来源主要包括上覆岩层重力、构造应力、水压附加应力、地温附加应力、化学附加应力等，通常采用垂向地应力、最大水平地应力、最小水平地应力三个参数来表征。地应力是影响钻井、水力压裂的关键基础参数之一，地应力在石油工程中有广泛的应用。就钻井工程而言，地应力是井壁稳定性分析的重要参数之一；就水力压力而言，地应力是影响裂缝起裂压力和扩展压力的关键性参数。
目前，石油工程领域研究地应力的方法很多，可以分为室内测试、矿场测试、测井解释以及其它方法等四个大类：①室内测试主要包括差应变测试、Kaiser声发射测试、波速各向异性测试、古地磁定向测试等；②矿场测试主要包括地破实验、水压致裂法、应力解除法、应力恢复法等；③测井解释主要包括声波测井解释、井壁崩落法、井壁诱导缝解释等；④而其它方法包括地质资料分析、地震资料预测、地应力数值模拟等。其中，矿场测试和室内测试的精度最高、测井解释次之、其它方法最低。但是，矿场测试和室内测试成本高、耗时长，而且只能获得取芯深度点的地应力大小和方向，所测得的数据十分有限；而测井解释具有纵向分辨率高、测井资料相对连续的特点，可以解释得到沿井深相对连续分布的地应力剖面，再结合多种方法和数据进行测井解释剖面的刻度，在一定程度上提升了地应力计算精度，这使得该方法被现场广泛采用。通过测井资料计算地层地应力时，垂向地应力通过中子密度测井资料进行积分求取，因此，通过测井资料解释求取地应力的难点在于求取最大和最小水平地应力。目前，最大和最小水平地应力的常规测井解释模型主要包括单轴应变模型、黄氏模型、组合弹簧模型和倾斜地层模型，但是，这些常规模型都是以各向同性介质为基础，但实际页岩地层是各向异性介质，因此，常规测井解释无法准确求取地层地应力。
对于各向异性地层中地应力的测井解释，Thiercelin和Plumb在各向同性地应力模型的基础上，区分了垂直和水平方向上力学性质的差别，建立了横观各向同性地层地应力模型；Higgins等采用横向各向同性地层地应力模型，对Baxter页岩进行地应力评价，进一步证实了各向异性的显著影响；Sayers、Gholami及邓金根等相继采用横向各向同性地层地应力模型开展了页岩地应力评价。但是，深层页岩地质力学环境复杂，深层页岩层理和裂隙十分发育，水平向层理表现出垂直横观各向同性特性、而垂直裂缝又表现出水平横观各向同性特性，综合表现出显著的正交各向异性特征，而现有的地应力评价模型主要是基于具有轴对称的横观各向同性本构，与正交各向异性地层特性不能完全匹配，其原因可以归结为以下三个方面：①现有技术在缺少阵列声波测井资料的情况下，很难准确确定横观各向同性地层岩石的5个独立刚度系数，对于更加复杂的正交各向异性问题，尚缺少有效表征正交各向异性地层岩石的12个独立刚度系数，这使得正交各向异性地层地应力的测井评价根本无法实现，常规测井资料更是无法实现地应力的准确评价；②杨氏模量和泊松比是影响各向异性地层地应力大小的重要参数，但是，利用测井、地震测量的波速计算的杨氏模量及泊松比为动态弹性参数，计算地应力需要静态弹性参数，页岩在垂直于层理面方向和平行与层理面方向的岩石弹性特征差异较大，而当前的岩石弹性参数动静态转换关系没有考虑这一特性，有可能导致较低的静态弹性参数转换精度；③Biot系数是评价地应力的关键参数之一，现有的地应力计算模型中，无法考虑Biot系数的各向异性特征，其主要原因是无法利用测井、地震等地球物理方法直接测量计算Biot系数所需的相应参数。因此，现有的横观各向同性地层地应力模型与正交各向异性地层特性不能完全匹配，无法准确计算正交各向异性地层的地应力。
发明内容
针对上述问题，本发明主要是克服现有技术中的不足之处，提出一种基于岩石物理建模的正交各向异性地层地应力预测方法。
本发明解决上述技术问题所提供的技术方案是：一种基于岩石物理建模的正交各向异性地层地应力预测方法，包括以下步骤：
步骤S1、通过测井、录井以及室内测试得到页岩物性参数，所述页岩物性参数包括矿物组成及含量、地层流体组成及含量、孔隙度、饱和度；
步骤S2、建立正交各向异性岩石物理模型；
步骤S3、根据正交各向异性岩石物理模型计算正交各向异性的刚度系数；
步骤S4、根据刚度系数与岩石弹性参数的关系式，计算动态岩石弹性参数；
式中：E1为沿最大水平主应力方向的杨氏模量，GPa；E2为沿最小水平主应力方向的杨氏模量，GPa；E3为沿垂向地应力方向的杨氏模量，GPa；ν12为垂直于最大水平主应力方向，沿最小水平主应力方向的泊松比，无因次；ν13为垂直于最大水平主应力方向，沿垂向地应力方向的泊松比，无因次；ν21为垂直于最小水平主应力方向，沿最大水平主应力方向的泊松比，无因次；ν23为垂直于最小水平主应力方向，沿垂向地应力方向的泊松比，无因次；ν31为垂直于垂向地应力方向，沿最大水平主应力方向的泊松比，无因次；ν32为垂直于垂向地应力方向，沿最小水平主应力方向的泊松比，无因次；
步骤S5、利用动态岩石弹性参数，并结合岩石弹性参数的动静态转换关系来计算静态岩石弹性参数；
步骤S6、根据下式计算正交各向异性岩石物理模型中各向异性的Biot系数；
式中：
为干正交各向异性页岩的刚度系数，GPa；
为VTI背景介质的刚度系数，GPa；α11、α22、α33为Biot系数，无因次；
步骤S7、基于Eaton法利用纵波时差来预测地层孔隙压力；
Pp＝σv-(σv-pw)(ACn/AC)c
式中：σv为垂向地应力，MPa；pw为地层水静液柱压力，MPa；ACn为该深度点正常压实趋势线纵波时差，μs/ft；AC为实际的纵波时差，μs/ft；Pp为地层孔隙压力，MPa；
步骤S8、利用以下关系式，计算正交各向异性页岩地层最大、最小水平主应力；
式中：σH为最大水平主应力，MPa；σh为最小水平主应力，MPa；εH为最大水平主应力方向所对应的最大水平主应变，无因次；εh为最小水平主应力方向所对应的最小水平主应变，无因次。
进一步的技术方案是，所述正交各向异性岩石物理模型的具体建立过程为：
步骤S21、利用Hashin-Shtrikman界限计算基质矿物混合物的等效弹性模量；
步骤S22、利用各向同性SCA模型和各向同性DEM模型，将无机质孔隙加入到基质矿物混合物中得到干基质，建立干基质页岩岩石物理模型并计算干基质的等效弹性模量；
步骤S23、利用各向异性SCA模型和各向异性DEM模型计算有机质矿物混合物的等效弹性刚度张量；
步骤S24、利用CL系数描述有机质矿物的成层分布，结合Bond变换正态分布和CL系数形成得到取向分布的有机质矿物混合物等效弹性刚度张量；
步骤S25、利用各向异性DEM模型，将有机质孔隙加入到取向分布的有机质矿物混合物中得到干有机质，建立干有机质页岩岩石物理模型并计算干有机质的等效弹性刚度张量；
步骤S26、利用各向异性DEM模型，将干基质添加到干有机质中，得到包含空孔隙的干页岩，建立干VTI页岩石物理模型并计算干页岩的刚度张量；
步骤S27、以干VTI页岩基质为背景，利用Schoenberg线性滑动模型的各向异性形式，将垂向分布的裂缝添加到其中，得到干ORT页岩岩石物理模型；
步骤S28、根据含气饱和度、含水饱和度，利用Wood公式计算孔隙内气-水混合物的体积模量，并转换成刚度张量；并结合干页岩呈ORT性质的情况，利用Brown-Korringa模型，将混合流体添加到干页岩中得到饱和流体页岩，建立饱和流体页岩岩石物理模型。
进一步的技术方案是，所述步骤S21的具体过程是：利用Hashin-Shtrikman界限计算基质矿物混合物体积模量和剪切模量的上下限，并取上下限平均作为基质矿物混合物的体积模量和剪切模量。
进一步的技术方案是，所述Hashin-Shtrikman界限计算基质矿物混合物的公式如下：
式中：KHS+为基质矿物混合物的等效体积模量上限，GPa；KHS-为基质矿物混合物的等效体积模量下限，GPa；μHS+为基质矿物混合物的等效剪切积模量上限，GPa；μHS-为基质矿物混合物的等效剪切模量下限，GPa。
进一步的技术方案是，所述步骤S22中各向同性SCA模型的计算公式如下：
式中：vi为第i种材料的体积分数，无因次；P*i和Q*i为第i种材料的几何因子，无因次；
为等效的体积模量，GPa；
为等效的剪切模量，GPa；
各向同性DEM模型的计算公式如下：
式中：K1、μ1分别为背景介质的体积模量和剪切模量，GPa；K2、μ2分布为包裹体的体积模量和剪切模量，GPa；v为包裹体的含量，无因次。
进一步的技术方案是，所述步骤S23中各向异性SCA模型的计算公式如下：
式中：
为SCA模型等效的刚度张量，GPa；
为第n相物质的Eshelby刚度张量，无因次；
为四阶单位刚度张量，无因；
为第n相物质的刚度张量，GPa；vn为第n相物质的体积分数，无因次；
各向异性DEM模型的计算公式如下：
式中：
为背景介质刚度张量，GPa；
为包裹体刚度张量，GPa；
为包裹体的Eshelby刚度张量，无因次；
为四阶单位刚度张量，无因次；v被添加项的体积，小数。
进一步的技术方案是，所述步骤S27中各向异性Schoenberg线性滑动模型计算公式如下：
式中：C为添加裂缝后的干ORT页岩刚度矩阵，GPa；c11b、c12b、c13b、c33b、c44b、c66b为不含垂向裂缝的干VTI页岩刚度系数，GPa；ΔN、ΔV、ΔH描述裂缝特征的弱度，无因次。
进一步的技术方案是，所述步骤S28中Wood公式如下所示：
式中：Kg为气体的体积模量，GPa；Kw为水的体积模量，GPa；Kf为混合流体的等效体积模量，GPa；sg为含气饱和度，无因次；μf为混合流体的等效体积模量，GPa；
Brown-Korringa模型的计算公式如下：
式中：
为饱和岩石的柔度张量，GPa-1；
为干岩石骨架的柔度张量，GPa-1；
为岩石基质的柔度张量，GPa-1；βfl、βgr分布为流体和岩石基质的可压缩系数，GPa-1；φ为孔隙度，无因次。
进一步的技术方案是，所述步骤S5中岩石弹性参数的动静态转换关系通过测量的波速计算岩心的动态岩石弹性参数并和三轴压缩实验得到的静态岩石弹性参数建立得到。
本发明具有以下有益效果：本发明构建了基于岩石物理建模的地应力评价体系，提高了静态岩石弹性参数的转换精度和地应力的预测精度，解决了正交各向异性页岩地层的地应力评价难题。利用该发明计算的正交各向异性地层地应力可以为后续的正交各向异性页岩地层井壁稳定性分析、压裂起裂分析等提供基础数据。
附图说明
图1是本发明的流程框图；
图2是本发明的建模流程图；
图3是x1井矿物及流体的测井解释图；
图4是岩石弹性参数的动、静态转换图；
图5是x1井岩石物理建模预测结果以及同测井测量结果的对比图；
图6是岩石静态弹性参数和各向异性Biot系数图；
图7是孔隙压力和各类模型计算的地应力图。
具体实施方式
下面结合实施例和附图对本发明做更进一步的说明。
如图1所示，本发明的一种基于岩石物理建模的正交各向异性地层地应力预测方法，具体实施过程如下：
步骤S1、结合测井、录井以及室内测试，得到建模所需的页岩物性参数，其x1井矿物及流体的测井解释图如图3所示，所述页岩物性参数包括矿物组成及含量、地层流体组成及含量、孔隙度、饱和度；
利用测井、录井以及室内测试等手段，获取页岩的矿物组分及相对含量，矿物组分包括基质矿物和有机质矿物两大类，基质矿物主要包括石英、长石、方解石、白云石、黄铁矿等，有机质矿物主要黏土和干酪根；
再利用室内测试或测井，获取页岩孔隙度，以及孔隙中所含的页岩气及水对应的包含度；页岩孔隙中的气体分为吸附气和游离气，基于Xu-White模型的原理，将孔隙分为无机质孔隙和有机质孔隙两类，其中无机质孔隙包含了粒内孔、粒间孔以及微裂缝；利用游离气和吸附气相对含量计算无机质孔隙度和有机质孔隙度；
吸附气的含量可由兰格缪尔等温吸附方程进行计算，具体公式如下：
式中，ga为吸附气含量，m3/t；Vlc为在油藏温度下经过TOC校正后的兰格缪尔体积，m3/t；Plt为油藏温度下的兰格缪尔压力，Mpa；p为储层压力，Mpa。
游离气的含量可由以下方程进行计算：
式中，gf为游离气含量，m3/t；ψ为常数，无因次；ρ为体积密度，g/cm3；Bg为气体体积压缩系数，无因次；φe为有效孔隙度，无因次；Sw为有效孔隙度，无因次。
有机质孔隙度和无机质孔隙度的计算方式如下：
式中，φo为有机质孔隙度，无因次；φg为无机质孔隙度，无因次；φ为总孔隙度，无因次；ga为吸附气含量，m3/t；gf为游离气含量，m3/t；
步骤S2、建立正交各向异性岩石物理模型，其具体过程如图2所示；
步骤S21、利用Hashin-Shtrikman界限计算基质矿物混合物的等效弹性模量；
利用Hashin-Shtrikman界限计算基质矿物混合物体积模量和剪切模量的上下限，并取上下限平均作为基质矿物混合物的体积模量和剪切模量；基质矿物是指除了黏土和干酪根以外的矿物，主要包括：石英、长石、方解石、白云石、黄铁矿等；
所述Hashin-Shtrikman界限计算基质矿物混合物的公式如下：
式中：KHS+为基质矿物混合物的等效体积模量上限，GPa；KHS-为基质矿物混合物的等效体积模量下限，GPa；μHS+为基质矿物混合物的等效剪切积模量上限，GPa；μHS-为基质矿物混合物的等效剪切模量下限，GPa；
本实施例中页岩岩石矿物及流体弹性模量如表1所示：
表1
| 物性参数 | K(GPa) | μ(GPa) | 密度(g/cm<sup>3</sup>) |
| 石英 | 37 | 4 | 2.65 |
| 方解石 | 76.8 | 32 | 2.71 |
| 长石 | 75.6 | 25.6 | 2.63 |
| 白云石 | 94.9 | 45 | 2.87 |
| 黄铁矿 | 147.4 | 132.5 | 4.93 |
| 黏土 | 25 | 9 | 2.55 |
| 干酪根 | 2.9 | 2.7 | 1.3 |
| 水 | 2.2 | 0 | 1.04 |
| 气 | 0.01 | 0 | 0.1 |
步骤S22、利用各向同性SCA模型和各向同性DEM模型，将无机质孔隙加入到基质矿物混合物中得到干基质，建立干基质页岩岩石物理模型并计算干基质的等效弹性模量；
具体是给定页岩基质岩石的临界孔隙度，利用各向同性SCA模型计算干岩石在临界孔隙度时候的体积模量和剪切模量，然后再利用各向同性DEM模型将孔隙度逐步调整至真实无机质孔隙度φg，以计算真实无机质孔隙度时的干页岩基质的体积模量和剪切模量；
各向同性SCA模型的计算公式如下：
式中：vi为第i种材料的体积分数，无因次；P*i为第i种材料的几何因子1，无因次；Q*i为第i种材料的几何因子2，无因次；
为等效的体积模量，GPa；
为等效的剪切模量，GPa；
各向同性DEM模型的计算公式如下：
式中：K1、μ1分别为背景介质的体积模量和剪切模量，GPa；K2、μ2分布为包裹体的体积模量和剪切模量，GPa；v为包裹体的含量，无因次；初始条件满足K*(0)＝K1，μ*(0)＝μ1；
步骤S23、利用各向异性SCA模型和各向异性DEM模型计算有机质矿物混合物的等效弹性刚度张量；
具体的是首先利用各向异性SCA模型计算干酪根和黏土含量各占50％情况下的弹性特征，然后用各向异性DEM模型将两者的体积含量调整至对应的百分比，这既保证了干酪根和黏土的相互连通，又避免了由于干酪根和黏土添加顺序的不同而导致混合物弹性模量的不对称；
各向异性SCA模型的计算公式如下：
式中：
为SCA模型等效的刚度张量，GPa；
为第n相物质的Eshelby刚度张量，无因次；
为四阶单位刚度张量，无因；
为第n相物质的刚度张量，GPa；vn为第n相物质的体积分数，无因次；
各向异性DEM模型的计算公式如下：
式中：
为背景介质刚度张量，GPa；
为包裹体刚度张量，GPa；
为包裹体的Eshelby刚度张量，无因次；
为四阶单位刚度张量，无因次；v被添加项的体积，小数；
步骤S24、利用CL系数描述有机质矿物的成层分布，结合Bond变换正态分布和CL系数形成得到取向分布的有机质矿物混合物等效弹性刚度张量；
步骤S241、利用扫描电镜资料，确定有机质颗粒随对称轴的偏转角度，给出表征有机质颗粒偏转角度的排列分布系数CL；
步骤S242、对于偏转到单个角度下的有机质颗粒，利用Bond变换将完全定向排列的有机质颗粒旋转到对应的偏转角度，得到该角度下对应的有机质颗粒的刚度张量；
步骤S243、将分布函数上所有排列角度的有机质颗粒的刚度矩阵利用V-R-H模型进行叠加，得到不同角度叠加的最终的有机质等效弹性刚度张量；
步骤S25、利用各向异性DEM模型，将有机质孔隙加入到取向分布的有机质矿物混合物中得到干有机质，建立干有机质页岩岩石物理模型并计算干有机质的等效弹性刚度张量；
步骤S26、利用各向异性DEM模型，将干基质添加到干有机质中，得到包含空孔隙的干页岩，建立干VTI页岩石物理模型并计算干页岩的刚度张量；
步骤S27、以干VTI页岩基质为背景，利用Schoenberg线性滑动模型的各向异性形式，将垂向分布的裂缝添加到其中，得到干ORT页岩岩石物理模型；
各向异性Schoenberg线性滑动模型计算公式如下：
式中：C为添加裂缝后的干ORT页岩刚度矩阵，GPa；c11b、c12b、c13b、c33b、c44b、c66b为不含垂向裂缝的干VTI页岩刚度系数，GPa；ΔN、ΔV、ΔH描述裂缝特征的弱度，无因次；
步骤S28、根据含气饱和度、含水饱和度，利用Wood公式计算孔隙内气-水混合物的体积模量，并转换成刚度张量；并结合干页岩呈VTI性质的情况，利用Brown-Korringa模型，将混合流体添加到干页岩中得到饱和流体页岩，建立饱和流体页岩的岩石物理模型；
Wood公式如下所示：
式中：Kg为气体的体积模量，GPa；Kw为水的体积模量，GPa；Kf为混合流体的等效体积模量，GPa；sg为含气饱和度，无因次；μf为混合流体的等效体积模量，GPa；
Brown-Korringa模型的计算公式如下：
式中：
为饱和岩石的柔度张量，GPa-1；
为干岩石骨架的柔度张量，GPa-1；
为岩石基质的柔度张量，GPa-1；βfl、βgr分布为流体和岩石基质的可压缩系数，GPa-1；φ为孔隙度，无因次；
步骤S3、根据正交各向异性岩石物理模型计算正交各向异性的刚度系数，其结果如图5(e)所示，刚度系数包括：c11、c12、c13、c22、c23、c33、c44、c55、c66；
步骤S4、根据刚度系数与岩石弹性参数的关系式，计算动态岩石弹性参数，其弹性参数包括3个杨氏模量、6个泊松比以及3个剪切模量，其中杨氏模量和泊松比与；
式中：E为杨氏模量，GPa；ν为泊松比，无量纲；μ为剪切模量，GPa；
步骤S5、利用动态岩石弹性参数，并结合岩石弹性参数的动静态转换关系来计算静态岩石弹性参数；
具体的是在利用波速测量实验测量页岩岩心波速的同时，对岩心开展了三轴岩石力学压缩实验，其结果如图4所示，用测量的波速计算岩心的动态岩石弹性参数并和三轴压缩实验得到的静态岩石弹性参数建立相应的动静态转换关系；
步骤S6、根据下式计算正交各向异性岩石物理模型中各向异性的Biot系数；
式中：
为干正交各向异性页岩的刚度系数，GPa；
为VTI背景介质的刚度系数，GPa；α11、α22、α33为Biot系数，无因次；
其岩石静态弹性参数和各向异性Biot系数如图6所示；
步骤S7、基于Eaton法利用纵波时差来预测地层孔隙压力；
Pp＝σv-(σv-pw)(ACn/AC)c
式中：σv为垂向地应力，MPa；pw为地层水静液柱压力，MPa；ACn为该深度点正常压实趋势线纵波时差，μs/ft；AC为实际的纵波时差，μs/ft；Pp为地层孔隙压力，MPa；
步骤S8、利用以下关系式，计算正交各向异性页岩地层最大、最小水平主应力，其结果如图7所示；
式中：σH为最大水平主应力，MPa；σh为最小水平主应力，MPa；εH为最大水平主应力方向所对应的最大水平主应变，无因次；εh为最小水平主应力方向所对应的最小水平主应变，无因次。
本实施例得到的孔隙压力和各类模型计算的地应力如图7所示；
图7(b)为不同地应力模型计算的最大水平主应力，图中A点为通过压裂施工曲线反算的最大水平主应力，其值为88.26MPa。B点为3161.8米处井下岩心声发射实验(基于Kaiser效应)获取的最大水平主应力，其值为95.21MPa。在B点对应的深度点，ORT地应力模型预测结果为98.98MPa，VTI地应力模型的预测结果为100.15MPa，ISO地应力模型的预测结果为84.79MPa，三种模型预测的地应力与实测结果的误差分别为3.96％、5.19％、10.94％，其中ORT模型的精度最高，ISO模型的精度最差。
图7(c)为不同地应力模型计算的最小水平主应力，图中A点为通过压裂施工曲线反算的最大水平主应力，其值为64.70MPa。B点为3161.8米处井下岩心声发射实验(基于Kaiser效应)获取的最小水平主应力，其值为74.13MPa。在B点对应的深度点，ORT地应力模型预测结果为75.55MPa，VTI地应力模型的预测结果为78.40MPa，ISO地应力模型的预测结果为68.09MPa，三种模型预测的地应力与实测结果的误差分别为1.92％、5.76％、8.15％，其中ORT模型的精度最高，ISO模型的精度最差。
以上所述，并非对本发明作任何形式上的限制，虽然本发明已通过上述实施例揭示，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容做出些变动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。
Claims (9)
1.一种基于岩石物理建模的正交各向异性地层地应力预测方法，其特征在于，包括以下步骤：
步骤S1、通过测井、录井以及室内测试得到页岩物性参数，所述页岩物性参数包括矿物组成及含量、地层流体组成及含量、孔隙度、饱和度；
步骤S2、建立正交各向异性岩石物理模型；
步骤S3、根据正交各向异性岩石物理模型计算正交各向异性的刚度系数；
步骤S4、根据刚度系数与岩石弹性参数的关系式，计算动态岩石弹性参数；
式中：E1为沿最大水平主应力方向的杨氏模量，GPa；E2为沿最小水平主应力方向的杨氏模量，GPa；E3为沿垂向地应力方向的杨氏模量，GPa；ν12为垂直于最大水平主应力方向，沿最小水平主应力方向的泊松比，无因次；ν13为垂直于最大水平主应力方向，沿垂向地应力方向的泊松比，无因次；ν21为垂直于最小水平主应力方向，沿最大水平主应力方向的泊松比，无因次；ν23为垂直于最小水平主应力方向，沿垂向地应力方向的泊松比，无因次；ν31为垂直于垂向地应力方向，沿最大水平主应力方向的泊松比，无因次；ν32为垂直于垂向地应力方向，沿最小水平主应力方向的泊松比，无因次；
步骤S5、利用动态岩石弹性参数，并结合岩石弹性参数的动静态转换关系来计算静态岩石弹性参数；
步骤S6、根据下式计算正交各向异性岩石物理模型中各向异性的Biot系数；
式中：
为干正交各向异性页岩的刚度系数，GPa；
为VTI背景介质的刚度系数，GPa；α11、α22、α33为Biot系数，无因次；
步骤S7、基于Eaton法利用纵波时差来预测地层孔隙压力；
Pp＝σv-(σv-pw)(ACn/AC)c
式中：σv为垂向地应力，MPa；pw为地层水静液柱压力，MPa；ACn为该深度点正常压实趋势线纵波时差，μs/ft；AC为实际的纵波时差，μs/ft；Pp为地层孔隙压力，MPa；
步骤S8、利用以下关系式，计算正交各向异性页岩地层最大、最小水平主应力；
式中：σH为最大水平主应力，MPa；σh为最小水平主应力，MPa；εH为最大水平主应力方向所对应的最大水平主应变，无因次；εh为最小水平主应力方向所对应的最小水平主应变，无因次。
2.根据权利要求1所述的一种基于岩石物理建模的正交各向异性地层地应力预测方法，其特征在于，所述正交各向异性岩石物理模型的具体建立过程为：
步骤S21、利用Hashin-Shtrikman界限计算基质矿物混合物的等效弹性模量；
步骤S22、利用各向同性SCA模型和各向同性DEM模型，将无机质孔隙加入到基质矿物混合物中得到干基质，建立干基质页岩岩石物理模型并计算干基质的等效弹性模量；
步骤S23、利用各向异性SCA模型和各向异性DEM模型计算有机质矿物混合物的等效弹性刚度张量；
步骤S24、利用CL系数描述有机质矿物的成层分布，结合Bond变换正态分布和CL系数形成得到取向分布的有机质矿物混合物等效弹性刚度张量；
步骤S25、利用各向异性DEM模型，将有机质孔隙加入到取向分布的有机质矿物混合物中得到干有机质，建立干有机质页岩岩石物理模型并计算干有机质的等效弹性刚度张量；
步骤S26、利用各向异性DEM模型，将干基质添加到干有机质中，得到包含空孔隙的干页岩，建立干VTI页岩石物理模型并计算干页岩的弹性张量；
步骤S27、以干VTI页岩基质为背景，利用Schoenberg线性滑动模型的各向异性形式，将垂向分布的裂缝添加到其中，得到干ORT页岩岩石物理模型；
步骤S28、根据含气饱和度、含水饱和度，利用Wood公式计算孔隙内气-水混合物的体积模量，并转换成刚度张量；并结合干页岩呈ORT性质的情况，利用Brown-Korringa模型，将混合流体添加到干页岩中得到饱和流体页岩，建立饱和流体的正交各向异性岩石物理模型。
3.根据权利要求2所述的一种基于岩石物理建模的正交各向异性地层地应力预测方法，其特征在于，所述步骤S21的具体过程是：利用Hashin-Shtrikman界限计算基质矿物混合物体积模量和剪切模量的上下限，并取上下限平均作为基质矿物混合物的体积模量和剪切模量。
4.根据权利要求2所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述Hashin-Shtrikman界限计算基质矿物混合物的公式如下：
式中：KHS+为基质矿物混合物的等效体积模量上限，GPa；KHS-为基质矿物混合物的等效体积模量下限，GPa；μHS+为基质矿物混合物的等效剪切积模量上限，GPa；μHS-为基质矿物混合物的等效剪切模量下限，GPa。
5.根据权利要求2所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述步骤S22中各向同性SCA模型的计算公式如下：
式中：vi为第i种材料的体积分数，无因次；P*i为第i种材料的几何因子1，无因次；Q*i为第i种材料的几何因子2；
为等效的体积模量，GPa；
为等效的剪切模量，GPa；
各向同性DEM模型的计算公式如下：
式中：K1、μ1分别为背景介质的体积模量和剪切模量，GPa；K2、μ2分别为包裹体的体积模量和剪切模量，GPa；v为包裹体的体积分数，无因次。
6.根据权利要求2所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述步骤S23中各向异性SCA模型的计算公式如下：
式中：
为SCA模型等效的刚度张量，GPa；
为第n相物质的Eshelby刚度张量，无因次；
为四阶单位刚度张量，无因；
为第n相物质的刚度张量，GPa；vn为第n相物质的体积分数，无因次；
各向异性DEM模型的计算公式如下：
式中：
为背景介质刚度张量，GPa；
为包裹体刚度张量，GPa；
为包裹体的Eshelby刚度张量，无因次；
为四阶单位刚度张量，无因次；v被添加项的体积，小数。
7.根据权利要求2所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述步骤S27中各向异性Schoenberg线性滑动模型计算公式如下：
式中：C为添加裂缝后的干ORT页岩刚度矩阵，GPa；c11b、c12b、c13b、c33b、c44b、c66b为不含垂向裂缝的干VTI页岩刚度系数，GPa；ΔN、ΔV、ΔH描述裂缝特征的弱度，无因次。
8.根据权利要求2所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述步骤S28中Wood公式如下所示：
式中：Kg为气体的体积模量，GPa；Kw为水的体积模量，GPa；Kf为混合流体的等效体积模量，GPa；sg为含气饱和度，无因次；μf为混合流体的等效体积模量，GPa；
Brown-Korringa模型的计算公式如下：
式中：
为饱和岩石的柔度张量，GPa-1；
为干岩石骨架的柔度张量，GPa-1；
为岩石基质的柔度张量，GPa-1；βfl、βgr分别为流体和岩石基质的可压缩系数，GPa-1；φ为孔隙度，无因次。
9.根据权利要求1所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述步骤S5中岩石弹性参数的动静态转换关系通过测量的波速计算岩心的动态岩石弹性参数并和三轴压缩实验得到的静态岩石弹性参数建立得到。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911239424.2A CN111460601B (zh) | 2019-12-06 | 2019-12-06 | 基于岩石物理建模的正交各向异性地层地应力预测方法 |
| US17/106,270 US11106843B2 (en) | 2019-12-06 | 2020-11-30 | Systems and methods for determining in-situ stresses based on orthotropic rock physics model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911239424.2A CN111460601B (zh) | 2019-12-06 | 2019-12-06 | 基于岩石物理建模的正交各向异性地层地应力预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111460601A true CN111460601A (zh) | 2020-07-28 |
| CN111460601B CN111460601B (zh) | 2020-11-17 |
Family
ID=71683148
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911239424.2A Active CN111460601B (zh) | 2019-12-06 | 2019-12-06 | 基于岩石物理建模的正交各向异性地层地应力预测方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US11106843B2 (zh) |
| CN (1) | CN111460601B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112525803A (zh) * | 2021-02-10 | 2021-03-19 | 西南石油大学 | 基于黏土壳状模型的砂岩孔隙度有效应力系数计算方法 |
| CN113189648A (zh) * | 2021-04-30 | 2021-07-30 | 西南石油大学 | 一种正交各向异性页岩脆性指数预测方法 |
| CN113189647A (zh) * | 2021-04-30 | 2021-07-30 | 西南石油大学 | 一种横观各向同性页岩地层脆性指数预测方法 |
| CN113255174A (zh) * | 2021-07-15 | 2021-08-13 | 西南石油大学 | 考虑岩石动态强度和混合破碎模式的钻齿力学计算方法 |
| CN113449427A (zh) * | 2021-07-01 | 2021-09-28 | 长江大学 | 一种确定页岩压力分布的方法 |
| CN113835119A (zh) * | 2021-09-18 | 2021-12-24 | 中国石油大学(华东) | 基于方位扩展弹性阻抗的页岩储层水平地应力预测方法 |
| CN115435936A (zh) * | 2022-09-20 | 2022-12-06 | 成都理工大学 | 一种深层致密储层地应力大小的计算方法 |
| CN116227158A (zh) * | 2023-01-13 | 2023-06-06 | 福建永福电力设计股份有限公司 | 一种海上风电吸力桩与导管架分界面等效刚度计算方法 |
| CN116401886A (zh) * | 2023-04-17 | 2023-07-07 | 应急管理部国家自然灾害防治研究院 | 一种适用于川藏地区的快速地应力估算方法 |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113570270B (zh) * | 2021-08-03 | 2023-05-09 | 东莞理工学院 | 一种定量建立流域演化弹性系数的方法 |
| CN113960659B (zh) * | 2021-10-14 | 2023-02-03 | 中国矿业大学 | 一种地震岩石物理驱动的煤层气储层含气量预测方法 |
| CN114088516B (zh) * | 2021-11-23 | 2023-07-18 | 贵州大学 | 一种基于岩心损伤各向异性的单孔声发射地应力测量方法 |
| CN114280272A (zh) * | 2021-12-13 | 2022-04-05 | 华能煤炭技术研究有限公司 | 煤岩岩石物理参数分析方法、***及电子设备 |
| CN115618512B (zh) * | 2022-10-12 | 2024-04-26 | 北京建筑大学 | 岩石节理各向异性特征量化方法、装置、设备及存储介质 |
| CN116025345B (zh) * | 2022-12-16 | 2024-04-26 | 中国石油天然气集团有限公司 | 一种横观各向同性地层静态有效应力系数的确定方法 |
| CN117647342B (zh) * | 2024-01-30 | 2024-04-19 | 中国地质大学（北京） | 一种基于声发射小波分析的煤储层地应力测定方法 |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100121623A1 (en) * | 2008-11-12 | 2010-05-13 | Terra Nova Sciences Llc | Methods and systems for constructing and using a subterranean geomechanics model spanning local to zonal scale in complex geological environments |
| US20130199787A1 (en) * | 2010-10-27 | 2013-08-08 | Bruce A. Dale | Method and System for Fracture Stimulation |
| CN105317430A (zh) * | 2014-07-29 | 2016-02-10 | 中国石油化工股份有限公司 | 测定各向异性地层水平井井壁围岩周向应力的方法 |
| CN105370268A (zh) * | 2015-10-23 | 2016-03-02 | 中国石油天然气集团公司 | 一种优化水平井分段压裂参数的方法和装置 |
| CN106295119A (zh) * | 2016-07-26 | 2017-01-04 | 中国石油大学(华东) | 一种页岩气地层地应力计算方法 |
| CN106597531A (zh) * | 2015-10-16 | 2017-04-26 | 中国石油化工股份有限公司 | 含垂直裂缝的页岩的波场传播特征的正演模拟方法 |
| CN106772673A (zh) * | 2016-11-29 | 2017-05-31 | 西南石油大学 | 一种页岩储层地应力预测建模方法 |
| CN108071389A (zh) * | 2017-12-10 | 2018-05-25 | 西南石油大学 | 一种页岩强度准则及页岩气井井眼坍塌压力预测方法 |
| CN108459346A (zh) * | 2017-02-20 | 2018-08-28 | 中国石油化工股份有限公司 | 一种页岩水平层理缝密度地震预测方法 |
| CN108595769A (zh) * | 2018-03-28 | 2018-09-28 | 北京航空航天大学 | 一种基于优化算法的减振器刚度模拟方法 |
| CN108952700A (zh) * | 2018-08-21 | 2018-12-07 | 西南石油大学 | 一种各向异性地层井壁破裂压力确定方法 |
| CN108957544A (zh) * | 2018-08-16 | 2018-12-07 | 中国石油大学(北京) | 近地表各向异性参数的测量方法、装置、地震计及介质 |
| CN109871507A (zh) * | 2017-12-05 | 2019-06-11 | 中国矿业大学(北京) | 正交各向异性煤层裂隙绝对渗透率计算方法 |
-
2019
- 2019-12-06 CN CN201911239424.2A patent/CN111460601B/zh active Active
-
2020
- 2020-11-30 US US17/106,270 patent/US11106843B2/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100121623A1 (en) * | 2008-11-12 | 2010-05-13 | Terra Nova Sciences Llc | Methods and systems for constructing and using a subterranean geomechanics model spanning local to zonal scale in complex geological environments |
| US20130199787A1 (en) * | 2010-10-27 | 2013-08-08 | Bruce A. Dale | Method and System for Fracture Stimulation |
| CN105317430A (zh) * | 2014-07-29 | 2016-02-10 | 中国石油化工股份有限公司 | 测定各向异性地层水平井井壁围岩周向应力的方法 |
| CN106597531A (zh) * | 2015-10-16 | 2017-04-26 | 中国石油化工股份有限公司 | 含垂直裂缝的页岩的波场传播特征的正演模拟方法 |
| CN105370268A (zh) * | 2015-10-23 | 2016-03-02 | 中国石油天然气集团公司 | 一种优化水平井分段压裂参数的方法和装置 |
| CN106295119A (zh) * | 2016-07-26 | 2017-01-04 | 中国石油大学(华东) | 一种页岩气地层地应力计算方法 |
| CN106772673A (zh) * | 2016-11-29 | 2017-05-31 | 西南石油大学 | 一种页岩储层地应力预测建模方法 |
| CN108459346A (zh) * | 2017-02-20 | 2018-08-28 | 中国石油化工股份有限公司 | 一种页岩水平层理缝密度地震预测方法 |
| CN109871507A (zh) * | 2017-12-05 | 2019-06-11 | 中国矿业大学(北京) | 正交各向异性煤层裂隙绝对渗透率计算方法 |
| CN108071389A (zh) * | 2017-12-10 | 2018-05-25 | 西南石油大学 | 一种页岩强度准则及页岩气井井眼坍塌压力预测方法 |
| CN108595769A (zh) * | 2018-03-28 | 2018-09-28 | 北京航空航天大学 | 一种基于优化算法的减振器刚度模拟方法 |
| CN108957544A (zh) * | 2018-08-16 | 2018-12-07 | 中国石油大学(北京) | 近地表各向异性参数的测量方法、装置、地震计及介质 |
| CN108952700A (zh) * | 2018-08-21 | 2018-12-07 | 西南石油大学 | 一种各向异性地层井壁破裂压力确定方法 |
Non-Patent Citations (2)
| Title |
|---|
| QING-DONG ZENG等: ""Numerical study of hydraulic fracture propagation accounting for rock anisotropy"", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 * |
| 桂俊川等: ""正交各向异性岩石弹性参数的空间展布"", 《西南石油大学学报（自然科学版）》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112525803A (zh) * | 2021-02-10 | 2021-03-19 | 西南石油大学 | 基于黏土壳状模型的砂岩孔隙度有效应力系数计算方法 |
| CN113189648B (zh) * | 2021-04-30 | 2022-03-11 | 西南石油大学 | 一种正交各向异性页岩脆性指数预测方法 |
| CN113189648A (zh) * | 2021-04-30 | 2021-07-30 | 西南石油大学 | 一种正交各向异性页岩脆性指数预测方法 |
| CN113189647A (zh) * | 2021-04-30 | 2021-07-30 | 西南石油大学 | 一种横观各向同性页岩地层脆性指数预测方法 |
| CN113189647B (zh) * | 2021-04-30 | 2022-03-11 | 西南石油大学 | 一种横观各向同性页岩地层脆性指数预测方法 |
| CN113449427A (zh) * | 2021-07-01 | 2021-09-28 | 长江大学 | 一种确定页岩压力分布的方法 |
| CN113449427B (zh) * | 2021-07-01 | 2022-08-16 | 长江大学 | 一种确定页岩压力分布的方法 |
| CN113255174A (zh) * | 2021-07-15 | 2021-08-13 | 西南石油大学 | 考虑岩石动态强度和混合破碎模式的钻齿力学计算方法 |
| CN113835119A (zh) * | 2021-09-18 | 2021-12-24 | 中国石油大学(华东) | 基于方位扩展弹性阻抗的页岩储层水平地应力预测方法 |
| CN115435936A (zh) * | 2022-09-20 | 2022-12-06 | 成都理工大学 | 一种深层致密储层地应力大小的计算方法 |
| CN115435936B (zh) * | 2022-09-20 | 2023-10-10 | 成都理工大学 | 一种深层致密储层地应力大小的计算方法 |
| CN116227158A (zh) * | 2023-01-13 | 2023-06-06 | 福建永福电力设计股份有限公司 | 一种海上风电吸力桩与导管架分界面等效刚度计算方法 |
| CN116227158B (zh) * | 2023-01-13 | 2024-01-09 | 福建永福电力设计股份有限公司 | 一种海上风电吸力桩与导管架分界面等效刚度计算方法 |
| CN116401886A (zh) * | 2023-04-17 | 2023-07-07 | 应急管理部国家自然灾害防治研究院 | 一种适用于川藏地区的快速地应力估算方法 |
| CN116401886B (zh) * | 2023-04-17 | 2023-11-24 | 应急管理部国家自然灾害防治研究院 | 一种适用于川藏地区的快速地应力估算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111460601B (zh) | 2020-11-17 |
| US20210173976A1 (en) | 2021-06-10 |
| US11106843B2 (en) | 2021-08-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111460601B (zh) | 基于岩石物理建模的正交各向异性地层地应力预测方法 | |
| CN111460602B (zh) | 基于岩石物理建模的横观各向同性地层地应力预测方法 | |
| Chen et al. | Experimental study of permeability change of organic-rich gas shales under high effective stress | |
| CN108071389B (zh) | 一种页岩气井井眼坍塌压力预测方法 | |
| CN103867198B (zh) | 碳酸盐岩天然气层地层密度判别法 | |
| CN104007482A (zh) | 一种基于各向异性有效场的泥页岩岩石物理模型方法 | |
| CN205538580U (zh) | 裂隙介质***渗透张量室内测定装置 | |
| Shafer | Recent advances in core analysis | |
| CN110501758B (zh) | 一种砂砾岩储层纵向连续脆性指数预测方法 | |
| CN109458176A (zh) | 碳酸盐岩储层压力的预测方法及其应用 | |
| CN105527652A (zh) | 一种岩石脆性的测井方法和装置 | |
| CN108304959B (zh) | 提高地层流体压力预测精度的方法 | |
| Pathi | Factors affecting the permeability of gas shales | |
| Crawford et al. | Determining static elastic anisotropy in shales from sidewall cores: impact on stress prediction and hydraulic fracture modeling | |
| CN100455768C (zh) | 确定扩径层段地层密度的方法 | |
| CN112412434B (zh) | 一种改进的疏松砂岩地应力计算方法 | |
| Karev et al. | Experimental studies of the deformation, destruction and filtration in rocks: a review | |
| CN110909486B (zh) | 一种正交各向异性页岩岩石物理模型的建立方法 | |
| Du | Anisotropic rock poroelasticity evolution in ultra‐low permeability sandstones under pore pressure, confining pressure, and temperature: experiments with Biot's coefficient | |
| CN108875115B (zh) | 一种确定岩石强度的方法 | |
| Al-Malikee et al. | Indirect prediction of rock elasticity and compressibility strength using well log data at selected sites within Rumaila Oilfield, Southern Iraq | |
| CN104712299B (zh) | 适合气井控水增气压裂的设计方法 | |
| Cai et al. | Comprehensive evaluation of rock mechanical properties and in-situ stress in tight sandstone oil reservoirs | |
| CN108195669B (zh) | 油藏围压条件下岩石静态力学参数的矫正及预测方法 | |
| Davies et al. | Determination of geomechanical properties of a typical Niger Delta reservoir rock using geophysical well logs |
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 | ||
| OL01 | Intention to license declared | ||
| OL01 | Intention to license declared |