CN109339771B - 一种页岩油气层孔隙压力预测方法及*** - Google Patents
一种页岩油气层孔隙压力预测方法及*** Download PDFInfo
- Publication number
- CN109339771B CN109339771B CN201710648490.XA CN201710648490A CN109339771B CN 109339771 B CN109339771 B CN 109339771B CN 201710648490 A CN201710648490 A CN 201710648490A CN 109339771 B CN109339771 B CN 109339771B
- Authority
- CN
- China
- Prior art keywords
- pore pressure
- elastic
- pressure
- shale
- elastic parameters
- 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 66
- 239000011148 porous material Substances 0.000 title claims abstract description 60
- 229930195733 hydrocarbon Natural products 0.000 title claims abstract description 18
- 150000002430 hydrocarbons Chemical class 0.000 title claims abstract description 18
- 239000004215 Carbon black (E152) Substances 0.000 title claims abstract description 17
- 238000012360 testing method Methods 0.000 claims abstract description 37
- 238000005553 drilling Methods 0.000 claims abstract description 35
- 239000012530 fluid Substances 0.000 claims abstract description 31
- 239000003079 shale oil Substances 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 2
- 238000009826 distribution Methods 0.000 abstract description 7
- 230000008901 benefit Effects 0.000 abstract description 6
- 230000015572 biosynthetic process Effects 0.000 description 30
- 239000008186 active pharmaceutical agent Substances 0.000 description 12
- 239000011435 rock Substances 0.000 description 7
- 239000000571 coke Substances 0.000 description 6
- 206010017076 Fracture Diseases 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 230000002159 abnormal effect Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000005056 compaction Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000011545 laboratory measurement Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000005416 organic matter Substances 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 206010016042 Facial bones fracture Diseases 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 229910052925 anhydrite Inorganic materials 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- OSGAYBCDTDRGGQ-UHFFFAOYSA-L calcium sulfate Chemical compound [Ca+2].[O-]S([O-])(=O)=O OSGAYBCDTDRGGQ-UHFFFAOYSA-L 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000018044 dehydration Effects 0.000 description 1
- 238000006297 dehydration reaction Methods 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- GUJOJGAPFQRJSV-UHFFFAOYSA-N dialuminum;dioxosilane;oxygen(2-);hydrate Chemical compound O.[O-2].[O-2].[O-2].[Al+3].[Al+3].O=[Si]=O.O=[Si]=O.O=[Si]=O.O=[Si]=O GUJOJGAPFQRJSV-UHFFFAOYSA-N 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- HVCNNTAUBZIYCG-UHFFFAOYSA-N ethyl 2-[4-[(6-chloro-1,3-benzothiazol-2-yl)oxy]phenoxy]propanoate Chemical compound C1=CC(OC(C)C(=O)OCC)=CC=C1OC1=NC2=CC=C(Cl)C=C2S1 HVCNNTAUBZIYCG-UHFFFAOYSA-N 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 229910052602 gypsum Inorganic materials 0.000 description 1
- 239000010440 gypsum Substances 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 231100000225 lethality Toxicity 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 229910052901 montmorillonite Inorganic materials 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000002574 poison Substances 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000011800 void material 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
- 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
-
- 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
- E21B47/06—Measuring temperature or pressure
-
- 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
- E21B49/005—Testing the nature of borehole walls or the formation by using drilling mud or cutting data
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
公开了一种页岩油气层孔隙压力预测方法及***。该方法包括：1)基于地震资料，抽取共反射点道集；2)基于所述共反射点道集，进行叠前弹性反演获得钻井液密度测试点处的N种弹性参数；3)对所述钻井液密度测试点处的N种弹性参数值与所述测试点处的孔隙压力值分别进行拟合；4)选取N种弹性参数中拟合优度最高的M个，建立压力预测的M元拟合公式，其中M<N；5)基于所述M元拟合公式，利用所选取的M种弹性参数值计算页岩油气层孔隙压力。相较于传统的API预测方法，根据本发明的页岩油气层孔隙压力预测方法的压力误差值整体较小，且分布相对更加平稳，可靠性强。
Description
技术领域
本发明涉及油气地球物理领域，更具体地，涉及一种页岩油气层孔隙压力预测方法及***。
背景技术
页岩气资源富集首先需要充沛的物质基础，如位于深水陆棚的沉积相带中、页岩发育有规模、富含有机质、脆性矿物含量高、干酪根处于成熟和过成熟阶段等。但对于涪陵地区海相沉积的页岩，在物质基础相对稳定的情况下，页岩储层的保存条件决定了页岩气的产量。这些保存条件包括：良好的顶底板、适中的埋藏深度、远离开启断裂、构造样式平缓、地层压力系数高等。在这些保存条件中，构造运动改造弱、开启断裂不发育是页岩气保存条件的重要因素，而地层压力系数则是页岩气保存条件的综合判断指标，储层高压-超高压往往表明保存条件好，地层能量足。
地层孔隙压力实质是地层孔隙或裂缝中流体(油、气、水)所具有的压力，而异常高压或低压是指孔隙压力高于或者低于静液压力。形成地层压力异常可以有多种因素，比如不平衡压实(快速沉积)、构造挤压(如断层)、水热增压、生烃作用、蒙脱土脱水、浓差与逆浓差、石膏/硬石膏转化、流体密度差异、水势面的不规则性(山脚钻井)、深部气体充填封存箱的分隔和抬升等。而对于焦石坝地区龙马溪-五峰组的页岩地层压力异常，一般认为主要是由于页岩层有机质富集造成巨大的生烃作用，使得优质页岩层段和区域富集了大规模的气体，形成异常高压的机制。实验室的岩心实验测量也证实：对于TOC含量不同的页岩岩样，随着地层压力的增加，其吸附气体的含量也是逐渐增加的。四川盆地东南部及其周缘实际勘探的效果也表明，无论是中国石油区块还是中国石化区块范围，地层流体压力系数越高对应的页岩气体产量往往也是越高，比如齐岳山断裂以西的焦石坝一期区块，其平均压力系数可到1.5，对应的试气产量为11-50万方每天，而齐岳山断裂以东的彭页区块平均压力系数仅为1.0左右，实际试气产量仅为2万方每天。而在焦石坝一期主体与周边多口井的实际测试产量结果与其地层压力系数均有良好的正相关关系，即页岩地层压力系数越高，对应的最终页岩气体产量也越大。可以说，异常高压是海相页岩气富集高产的关键因素。
此外，地层压力的预测可以降低钻井风险、提高效率、降低成本(Gutierrez，2006)，对于页岩甜点品质预测和工程布置施工至关重要。
地层压力预测方法，一般分为钻前地震资料预测、随钻地层压力监测及钻后测井检测三类，在勘探阶段，主要关注的是如何利用地震资料预测异常地层压力的分布特征。国内外关于地层压力预测方法的研究已历经数十年，尽管描述地下压力信息的方法和模型有多种，从各方法的思路原理来看，真正表征地下压力平衡的公式是1926年的Terzaghi公式：
PP＝Po-Pe (1)
其中，Pp代表地层孔隙压力；Po代表上覆地层压力，Pe代表垂直有效应力。上覆压力是指上覆岩石骨架和孔隙空间流体的总体重量引起的压力，与上覆岩层的厚度、骨架密度和空隙流体密度有关。垂直有效应力是指地层骨架或者岩层承受的垂直方向的压实作用产生的应力，也称为骨架应力，不可直接测量。
目前已经有多种模型或方法描述公式3.6.1中的孔隙压力或者有效应力，进而直接或者间接对地层压力进行预测，如等效深度法(1965)、Eaton法(1976)、Bowers法(1995)、Filippone法(1982)、Eberhart-Phillips法(1985)等。这些方法的选择特别是参数的确立都存在极强的区域特征，整体看来，操作相对繁琐。与之相对应，由于上述压力预测模型均为区域性的实验室或者井筒测试参数筛选后所建立，常常与某些敏感参数存在较强的对应关系，比如地层速度、声波时差和密度等参数。因而，在压力资料比较丰富的情况下，可以进行弹性参数筛选，建立适合区域性压力预测的特有模型，比如API预测方法(2015)等。然而这类方法的精度仍然与弹性参数的相关性存在较大关联，且存在极大的区域性。因此，有必要开发一种分布平稳、可靠性强的页岩油气层孔隙压力预测方法及***。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
随着页岩气资源勘探开发的深入，特别是针对我国涪陵页岩气田产能建设区域，页岩油气层的压力预测精度直接决定了页岩气的最终产量。焦石坝南地区龙马溪-五峰组页岩层经历了剧烈的构造变动，地表地形复杂多变，现有的钻探和压力测试资料稀少，难于开展基于压力平衡法的地层压力有效预测，也难于确定有效的压力预测弹性参数和组合形式。为了解决上述问题，本发明基于钻井液密度测试点处计算孔隙压力值与相应纵波阻抗、横波阻抗、纵波速度、横波速度、密度、杨氏模量、泊松比和纵横波速度比等弹性参数的关系，测试不同组合指数拟合下的拟合优度，最终筛选建立了一种页岩油气层孔隙压力地震的多元指数预测方法，对焦石坝南地区页岩层孔隙压力进行预测。
根据本发明的一方面，提出一种页岩油气层孔隙压力预测方法。该方法主要包括：
1)基于地震资料，抽取共反射点道集；
2)基于所述共反射点道集，进行叠前弹性反演获得钻井液密度测试点处的N种弹性参数；
3)对所述钻井液密度测试点处的N种弹性参数值与所述测试点处的孔隙压力值分别进行拟合；
4)选取N种弹性参数中拟合优度最高的M个，建立压力预测的M元拟合公式，其中M<N；
5)基于所述M元拟合公式，利用所选取的M种弹性参数值计算页岩油气层孔隙压力。
优选地，所述N种弹性参数包括纵波阻抗、横波阻抗、纵波速度、横波速度、密度、杨氏模量、泊松比、纵横波速度比。
优选地，在步骤3)中对所述钻井液密度测试点处的N种弹性参数值与所述测试点处的孔隙压力值分别进行指数形式的拟合。
优选地，M元拟合公式为：
其中，Pp为页岩油气层孔隙压力，Xi为对应的弹性参数，为所选择的拟合优度最高的M个弹性参数之一，Ai和Bi为不同弹性参数对应指数拟合公式的系数，Ri为不同弹性参数对应的拟合优度。
优选地，步骤1)包括对地震资料进行保幅处理获得共反射点道集。
根据本发明的另一方面，提出一种页岩油气层孔隙压力预测***。该***包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：
1)基于地震资料，抽取共反射点道集；
2)基于所述共反射点道集，进行叠前弹性反演获得钻井液密度测试点处的N种弹性参数；
3)对所述钻井液密度测试点处的N种弹性参数值与所述测试点处的孔隙压力值分别进行拟合；
4)选取N种弹性参数中拟合优度最高的M个，建立压力预测的M元拟合公式，其中M<N；
5)基于所述M元拟合公式，利用所选取的M种弹性参数值计算页岩油气层孔隙压力。
优选地，所述N种弹性参数包括纵波阻抗、横波阻抗、纵波速度、横波速度、密度、杨氏模量、泊松比、纵横波速度比。
优选地，在步骤3)中对所述钻井液密度测试点处的N种弹性参数值与所述测试点处的孔隙压力值分别进行指数形式的拟合。
优选地，M元拟合公式为：
其中，Pp为页岩油气层孔隙压力，Xi为对应的弹性参数，为所选择的拟合优度最高的M个弹性参数之一，Ai和Bi为不同弹性参数对应指数拟合公式的系数，Ri为不同弹性参数对应的拟合优度。
优选地，步骤1)包括对地震资料进行保幅处理获得共反射点道集。
相较于传统的API预测方法，根据本发明的页岩油气层孔隙压力预测方法的压力误差值整体较小，且分布相对更加平稳，可靠性强。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1为根据本发明的示例性实施方案的页岩油气层孔隙压力预测方法的流程图；
图2a-2h为不同弹性参数与压力的交汇及指数拟合曲线图，其中，图2a中的弹性参数为纵波阻抗Ip，图2b中的弹性参数为横波阻抗Is，图2c中的弹性参数为纵波速度Vp，图2d中的弹性参数为横波速度Vs，图2e中的弹性参数为密度ρ，图2f中的弹性参数为杨氏模量Ymod，图2g中的弹性参数为泊松比σ，图2h中的弹性参数为纵横波速度比Vp/Vs；
图3为实际样点处实际压力与两种方法预测压力对比图，其中，虚线为理论值，
为本发明的预测方法，
为传统的API预测方法；
图4为实际样点处两种方法预测压力经绝对误差对比图，其中，
为本发明的预测方法，
为传统的API预测方法；
图5为焦石坝南优质页岩段依据API预测方法预测的地层压力分布图；
图6为焦石坝南优质页岩段依据本发明方法预测的地层压力分布图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
依据压力平衡原理(公式1)，主流的压力预测方法为首先依据密度信息计算上覆压力，而后选择合适的模型和参数对孔隙压力或有效应力进行模拟，从而表征地层的孔隙压力。地层压力的预测精度取决于资料的准确和详细程度，且其预测手段操作流程复杂，影响因素较多，实际工作中必须分阶段进行。
实际上，地层压力常与某些敏感参数存在较强的对应关系，比如地层速度、声波时差和密度等参数。因而，在压力资料比较丰富的情况下，可以进行弹性参数筛选，建立适合区域性压力预测的特有模型。比如API法：
其中，σ为泊松比，Ip代表纵波阻抗，C为压力系数因子。
然而，不同弹性参数在不同地区对压力信息有不同程度的反映，表征地震速度与孔隙压力关系的Eaton模型(1976)、地震速度与垂直有效应力关系的Bowers模型(1976)、Filippone模型(1982)等都是基于大量数据统计和实验室测量建立的，每一个参数对压力信息的贡献形式都有严格论证。相比之下，基于弹性参数拟合的经验公式往往缺乏实验室论证，应用范围有限。
焦石坝南地区页岩层经历了剧烈的构造变动，地表地形复杂多变，现有的钻探和压力测试资料稀少，难于开展基于压力平衡法的地层压力有效预测，也难于确定有效的压力预测弹性参数和组合形式。为此，本发明基于钻井液密度测试点处计算孔隙压力值与不同弹性参数进行不同类型的拟合，筛选敏感参数建立多元预测模型，从而对页岩层的地层压力进行预测。
以下参照图1详细描述根据本发明示例性实施方案的页岩油气层孔隙压力预测方法。
该方法主要包括：
步骤1：基于地震资料，抽取共反射点道集。
在一个示例中，对地震资料进行保幅处理，抽取高品质共反射点(CRP)道集。
步骤2：基于所述共反射点道集，进行叠前弹性反演获得钻井液密度测试点处的N种弹性参数。
在一个示例中，N种弹性参数包括纵波阻抗、横波阻抗、纵波速度、横波速度、密度、杨氏模量、泊松比、纵横波速度比。本领域技术人员应当理解，N种弹性参数可以是任意适当的弹性参数类型。
步骤3：对所述钻井液密度测试点处的N种弹性参数值与所述测试点处的孔隙压力值分别进行拟合。
页岩钻探一般为平衡钻井，所采用钻井液密度对应重力基本近似于地层压力。某一深度点处钻井液密度产生的压力数据为：
其中，H为上覆岩层的垂直高度，ρ为钻井液密度，g为重力加速度。例如，依据焦石坝南工区内焦页5-7井直井段钻井液密度信息，即可计算出相应的压力数据。
实际工作中，变量间未必都有线性关系，如服药后血药浓度与时间的关系；疾病疗效与疗程长短的关系；毒物剂量与致死率的关系等常呈曲线关系。曲线拟合是指选择适当的曲线类型来拟合观测数据，并用拟合的曲线方程分析两变量间的关系。用连续曲线近似地刻画或比拟离散点组成坐标之间的函数关系，可以通过简单变量直接表征非线性化的资料和问题。
图2a-2h为焦石坝南焦页5-7井垂直井段42个实测钻井液密度点处压力数据和相应位置处的弹性参数交汇图，其中，图2a中的弹性参数为纵波阻抗Ip，图2b中的弹性参数为横波阻抗Is，图2c中的弹性参数为纵波速度Vp，图2d中的弹性参数为横波速度Vs，图2e中的弹性参数为密度ρ，图2f中的弹性参数为杨氏模量Ymod，图2g中的弹性参数为泊松比σ，图2h中的弹性参数为纵横波速度比Vp/Vs，Pp为页岩油气层孔隙压力。可以看出，采用指数形式对不同交汇组合有不同吻合程度的拟合。
拟合的形式则需要根据其与实际非线性点的吻合关系予以评价，一般采用拟合优度来表示。拟合优度是指回归直线对观测值的拟合程度，是表达因变量与所有自变量之间的总体关系。度量拟合优度的统计量是可决系数(亦称确定系数)R2，取值范围是[0，1]。具体R等于回归平方和在总平方和中所占的比率，即回归方程所能解释的因变量变异性的百分比。R2的值越接近1，说明回归直线对观测值的拟合程度越好；反之，R2的值越接近0，说明回归直线对观测值的拟合程度越差。下表为图2a-2h的8个交汇图对应的指数拟合公式和相应拟合优度，其中，Pp为地层压力，也即孔隙压力。
表1不同弹性参数与压力的指数拟合形式和拟合优度
从表1中可以看出，压力-纵波阻抗、压力-纵波速度有最高的拟合优度，高于压力-横波波阻抗和压力-横波速度，因为页岩有机质孔隙内大量生烃，导致含气量升高和地层压力增加，对应的纵波速度往往会降低，而横波为沿岩石骨架传播的剪切波，敏感程度不高。此外，压力-杨氏模量、压力-泊松比均有较强的拟合优度。
步骤4：选取N种弹性参数中拟合优度最高的M个，建立压力预测的M元拟合公式，其中M<N。
选取弹性参数中拟合优度最高的几个，开展压力预测的多元拟合公式构建。例如，由于压力-密度却有最低的拟合优度，纵波阻抗与纵波速度本身又有极强的相关性，故可以只选择拟合优度最高的压力-纵波阻抗组合，结合压力-杨氏模量、压力-泊松比等两外两个较高拟合优度的组合构建三元指数预测模型。
不同弹性参数的指数拟合通用形式为：
其中，Pp为页岩油气层孔隙压力，Xi为弹性参数，Ai和Bi为弹性参数对应指数拟合公式的系数。
进一步引入拟合优度构建页岩油气层孔隙压力地震多元指数预测模型为：
其中，Pp为页岩油气层孔隙压力，Xi为对应的弹性参数，为所选择的拟合优度最高的M个弹性参数之一，Ai和Bi为不同弹性参数对应指数拟合公式的系数，Ri为不同弹性参数对应的拟合优度。
例如，在选择表1中拟合优度最高的纵波阻抗、杨氏模量、泊松比这三个弹性参数的情况下，所建立的三元指数预测公式为：
PP＝180.61·e-2E-7·Ip+30.03·e-2E-11·Ymod+43.27·e-6.059·σ (5)
其中，Ip为纵波阻抗，Ymod为杨氏模量，σ为泊松比。
步骤5：基于所述M元拟合公式，利用所选取的M种弹性参数值计算页岩油气层孔隙压力。
例如，选择表1中拟合优度最高的纵波阻抗、杨氏模量、泊松比这三个弹性参数的情况下，则可以基于公式(5)，利用钻井液密度测试点处的纵波阻抗、杨氏模量、泊松比这三个弹性参数的值来计算在这些测试点处的孔隙压力。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
分别选用传统的API预测方法和本发明提出的地震多元指数压力预测方法(在本示例中为地震三元指数压力预测方法，所选择的是拟合优度最高的纵波阻抗、杨氏模量、泊松比这三个弹性参数)对焦页5-7井垂直井段42个实测钻井液密度点处压力数据进行预测，并与依据钻井液密度计算的压力数据进行对比，如图3所示。API预测方法的参数C依据与实际压力误差最小原则进行选择，具体为C＝2.37E16。
显然，两种预测方法预测的结果大体趋势一致，基本与理论值相同，这是因为二者都是依据纵波阻抗和泊松比建立的公式。但是由于引入了杨氏模量的信息，大致可以看出本发明预测方法的误差更小，尤其是在前18个样点位置，基于API预测方法预测的精度要明显低于本发明方法预测的结果。更加细节的对比可以从图4的绝对误差对比上观测到，尽管传统API预测方法可以调整系数参数以使得预测值的主体分布与理论值相近，但是断续在前12个样点处、第20个样点和第31个样点处存在极大的误差。与之相比，本发明方法预测的压力误差值整体较小，且分布相对更加平稳，可靠性强。
图5和图6为对于焦石坝南龙马溪-五峰组优质页岩层段分别依据API预测方法和本发明提出的地震三元指数预测方法预测的地层压力分布。两种方法所采用的纵波阻抗、杨氏模量和泊松比等弹性参数为依据该区叠前弹性反演而得，叠前反演前利用多次波和线性噪音去除、随机噪音衰减和地表一致性剩余静校正等手段提高道集品质，反演过程中采用严格质控策略特别是加重密度模型约束程度来提高反演精度。
可以看出，两种方法预测的地层异常高压区域均位于平桥断背斜、白马向斜优势部位、乌江背斜带和沙子沱断鼻中北部，以及平桥断背斜以西北的大片区域。然而，两者在细节上有所差别，基于页岩油气层孔隙压力地震三元指数预测方法预测的地层压力变化更加连续，且在上述几个目标地带范围更大，与地质情况更加吻合。特别地，传统API方法预测的结果在焦页8井位置处压力值较低，与实际产量不符，而本发明方法预测的结果显示焦页8井位于平桥断背斜异常高压带的边部，对应高的页岩气产量，与实际生产情况吻合。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (8)
1.一种页岩油气层孔隙压力预测方法，其特征在于，包括：
1)基于地震资料，抽取共反射点道集；
2)基于所述共反射点道集，进行叠前弹性反演获得钻井液密度测试点处的N种弹性参数；
3)对所述钻井液密度测试点处的N种弹性参数值与所述测试点处的孔隙压力值分别进行拟合；
4)选取N种弹性参数中拟合优度最高的M个，建立压力预测的M元拟合公式，其中M<N；
5)基于所述M元拟合公式，利用所选取的M种弹性参数值计算页岩油气层孔隙压力；
M元拟合公式为：
其中，Pp为页岩油气层孔隙压力，Xi为对应的弹性参数，为所选择的拟合优度最高的M个弹性参数之一，Ai和Bi为不同弹性参数对应指数拟合公式的系数，Ri为不同弹性参数对应的拟合优度。
2.根据权利要求1所述的页岩油气层孔隙压力预测方法，其特征在于，所述N种弹性参数包括纵波阻抗、横波阻抗、纵波速度、横波速度、密度、杨氏模量、泊松比、纵横波速度比。
3.根据权利要求1所述的页岩油气层孔隙压力预测方法，其特征在于，在步骤3)中对所述钻井液密度测试点处的N种弹性参数值与所述测试点处的孔隙压力值分别进行指数形式的拟合。
4.根据权利要求1所述的页岩油气层孔隙压力预测方法，其特征在于，步骤1)包括对地震资料进行保幅处理获得共反射点道集。
5.一种页岩油气层孔隙压力预测***，其特征在于，所述***包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：
1)基于地震资料，抽取共反射点道集；
2)基于所述共反射点道集，进行叠前弹性反演获得钻井液密度测试点处的N种弹性参数；
3)对所述钻井液密度测试点处的N种弹性参数值与所述测试点处的孔隙压力值分别进行拟合；
4)选取N种弹性参数中拟合优度最高的M个，建立压力预测的M元拟合公式，其中M<N；
5)基于所述M元拟合公式，利用所选取的M种弹性参数值计算页岩油气层孔隙压力；
M元拟合公式为：
其中，Pp为页岩油气层孔隙压力，Xi为对应的弹性参数，为所选择的拟合优度最高的M个弹性参数之一，Ai和Bi为不同弹性参数对应指数拟合公式的系数，Ri为不同弹性参数对应的拟合优度。
6.根据权利要求5所述的页岩油气层孔隙压力预测***，其特征在于，所述N种弹性参数包括纵波阻抗、横波阻抗、纵波速度、横波速度、密度、杨氏模量、泊松比、纵横波速度比。
7.根据权利要求5所述的页岩油气层孔隙压力预测***，其特征在于，在步骤3)中对所述钻井液密度测试点处的N种弹性参数值与所述测试点处的孔隙压力值分别进行指数形式的拟合。
8.根据权利要求5所述的页岩油气层孔隙压力预测***，其特征在于，步骤1)包括对地震资料进行保幅处理获得共反射点道集。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710648490.XA CN109339771B (zh) | 2017-08-01 | 2017-08-01 | 一种页岩油气层孔隙压力预测方法及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710648490.XA CN109339771B (zh) | 2017-08-01 | 2017-08-01 | 一种页岩油气层孔隙压力预测方法及*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109339771A CN109339771A (zh) | 2019-02-15 |
| CN109339771B true CN109339771B (zh) | 2022-04-12 |
Family
ID=65291230
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710648490.XA Active CN109339771B (zh) | 2017-08-01 | 2017-08-01 | 一种页岩油气层孔隙压力预测方法及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109339771B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113514890B (zh) * | 2021-03-18 | 2022-12-09 | 中国石油大学(华东) | 一种使用地震数据预测地层压力系数的方法、装置及设备 |
| CN113625364B (zh) * | 2021-08-18 | 2022-08-02 | 西南石油大学 | 一种基于双重校正的泥页岩地层孔隙压力的计算方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101163990A (zh) * | 2005-03-08 | 2008-04-16 | 地质力学国际公司 | 用于孔压预测的量化风险评估 |
| CN102640018A (zh) * | 2009-12-18 | 2012-08-15 | 雪佛龙美国公司 | 用于电缆和lwd测井数据的石油物理和地球物理地层评价的工作流程 |
| CN103628866A (zh) * | 2013-11-22 | 2014-03-12 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 获取地层的地层压力系数的方法 |
| CN104267429A (zh) * | 2014-09-30 | 2015-01-07 | 中国石油天然气股份有限公司 | 确定地层压力的方法及装置 |
| CN104698492A (zh) * | 2013-12-09 | 2015-06-10 | 中国石油天然气股份有限公司 | 一种计算异常地层压力的方法 |
| CN105445791A (zh) * | 2015-11-25 | 2016-03-30 | 成都理工大学 | 一种基于多种地震属性的地层孔隙压力预测方法 |
| CN106285642A (zh) * | 2015-05-14 | 2017-01-04 | 中国石油化工股份有限公司 | 一种基于地震资料的地层孔隙压力预测方法 |
| CN106368691A (zh) * | 2015-07-24 | 2017-02-01 | 中国石油化工股份有限公司 | 基于岩石物理地震信息三维异常孔隙压力预测方法 |
| CN106845086A (zh) * | 2016-12-30 | 2017-06-13 | 中国石油天然气集团公司 | 地层压力计算方法及装置 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6751558B2 (en) * | 2001-03-13 | 2004-06-15 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
| GB2404988B (en) * | 2002-04-10 | 2006-04-12 | Schlumberger Technology Corp | Method,apparatus and system for pore pressure prediction in presence of dipping formations |
| US7493241B2 (en) * | 2003-07-23 | 2009-02-17 | Lee Wook B | 3D velocity modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved for curved-ray prestack time migration and for such migration followed by prestack depth migration |
| FR2893421B1 (fr) * | 2005-11-14 | 2007-12-21 | Inst Francais Du Petrole | Methode d'evaluation quantitative des pressions de fluides et de detection des surpressions d'un milieu souterrain. |
| EP2153246B1 (en) * | 2007-05-09 | 2015-09-16 | ExxonMobil Upstream Research Company | Inversion of 4d seismic data |
| CN105242307B (zh) * | 2015-09-22 | 2016-08-31 | 中国石油大学(北京) | 复杂碳酸盐储层地震孔隙度获取方法及装置 |
-
2017
- 2017-08-01 CN CN201710648490.XA patent/CN109339771B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101163990A (zh) * | 2005-03-08 | 2008-04-16 | 地质力学国际公司 | 用于孔压预测的量化风险评估 |
| CN102640018A (zh) * | 2009-12-18 | 2012-08-15 | 雪佛龙美国公司 | 用于电缆和lwd测井数据的石油物理和地球物理地层评价的工作流程 |
| CN103628866A (zh) * | 2013-11-22 | 2014-03-12 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 获取地层的地层压力系数的方法 |
| CN104698492A (zh) * | 2013-12-09 | 2015-06-10 | 中国石油天然气股份有限公司 | 一种计算异常地层压力的方法 |
| CN104267429A (zh) * | 2014-09-30 | 2015-01-07 | 中国石油天然气股份有限公司 | 确定地层压力的方法及装置 |
| CN106285642A (zh) * | 2015-05-14 | 2017-01-04 | 中国石油化工股份有限公司 | 一种基于地震资料的地层孔隙压力预测方法 |
| CN106368691A (zh) * | 2015-07-24 | 2017-02-01 | 中国石油化工股份有限公司 | 基于岩石物理地震信息三维异常孔隙压力预测方法 |
| CN105445791A (zh) * | 2015-11-25 | 2016-03-30 | 成都理工大学 | 一种基于多种地震属性的地层孔隙压力预测方法 |
| CN106845086A (zh) * | 2016-12-30 | 2017-06-13 | 中国石油天然气集团公司 | 地层压力计算方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| Use of reflection tomography to predict pore pressure in overpressured reservoir sands;C.M. Sayers等;《the 2003 SEG Annual Meeting》;20031031;第1-4页 * |
| 地震地层压力预测综述;孙武亮,等;《勘探地球物理进展》;20071231;第30卷(第6期);第428-432页 * |
| 基于三维地震数据的地层压力预测方法研究;谭峰;《中国优秀硕士学位论文全文数据库 基础科学辑》;20170315;第12-32页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109339771A (zh) | 2019-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106951660B (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
| Close et al. | Integrated workflows for shale gas and case study results for the Horn River Basin, British Columbia, Canada | |
| CN110321595B (zh) | 一种测井提取静态品质系数的断层封闭性评价方法 | |
| Liu et al. | Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs | |
| CN105931125B (zh) | 一种致密油分段多簇体积压裂水平井产量预测方法 | |
| Alfataierge et al. | 3D hydraulic-fracture simulation integrated with 4D time-lapse multicomponent seismic and microseismic interpretations, Wattenberg Field, Colorado | |
| CN109339771B (zh) | 一种页岩油气层孔隙压力预测方法及*** | |
| Buijs et al. | An engineered approach to hydraulically fracture the Vaca Muerta shale | |
| Alimahomed* et al. | Impact of pore pressure on modeled hydraulic fracture geometry and well spacing in the East Duvernay Shale Basin, Canada | |
| Sabea et al. | Geological model of the Khabour Reservoir for studying the gas condensate blockage effect on gas production, Akkas Gas Field, Western Iraq | |
| Soltanzadeh* et al. | Application of mechanical and mineralogical rock properties to identify fracture fabrics in the Devonian Duvernay formation in Alberta | |
| CN110320562B (zh) | 页岩气储层水平井中声波时差的校正方法 | |
| Cherian et al. | Understanding development drivers in horizontal wellbores in the Midland Basin | |
| Soltanzadeh et al. | A regional review of geomechanical drilling experience and problems in the Duvernay Formation in Alberta | |
| Abass et al. | Integration of mud logging and wire-line logging to detect overpressure zones: a case study of middle Miocene Kareem Formation in Ashrafi oil field, Gulf of Suez, Egypt | |
| Wüst et al. | Interbedded Carbonate and Calcareous Shales of the Devonian Duvernay Formation of Alberta, Canada: Implications for Completion Due to High Variability of Geomechanical Properties | |
| Peza et al. | 3-D Integrated Workflow for Understanding the Fracture Interference and Its Impact into the Gas Production of the Woodford Shale | |
| Martemyanov et al. | Analytic modelling for wellbore stability analysis | |
| Grazulis | Analysis of stress and geomechanical properties in the Niobrara Formation of Wattenberg Field, Colorado, USA | |
| Zeng et al. | A new method for assessing stage-based hydraulic fracturing quality in Wolfcamp formation of Permian Basin | |
| Hai et al. | The Characteristics of Hydraulic Fracture Growth in Woodford Shale, the Anadarko Basin, Oklahoma | |
| Altaf et al. | 3D Geomechanical Modelling for Unconventional Sweet Spot Identification and Optimized Reservoir Development | |
| CN110297280A (zh) | 一种碳酸盐岩超压空间分布特征的预测方法 | |
| Ezati et al. | A new workflow to quantitative 3D screening of suitable zones for hydraulic fracturing | |
| Chatterjee et al. | Stress modeling for the oil and gas fields of Krishna-Godavari and Cauvery basins, India, using finite element technique |
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 |