CN110579797A - 一种泥页岩储层含气量地球物理定量预测方法 - Google Patents
一种泥页岩储层含气量地球物理定量预测方法 Download PDFInfo
- Publication number
- CN110579797A CN110579797A CN201910823702.2A CN201910823702A CN110579797A CN 110579797 A CN110579797 A CN 110579797A CN 201910823702 A CN201910823702 A CN 201910823702A CN 110579797 A CN110579797 A CN 110579797A
- Authority
- CN
- China
- Prior art keywords
- shale reservoir
- gas content
- data volume
- wave impedance
- longitudinal wave
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 9
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 9
- 239000011435 rock Substances 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 6
- 239000007789 gas Substances 0.000 description 68
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 238000012360 testing method Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000003345 natural gas Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000000611 regression analysis Methods 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000010249 in-situ analysis Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Mining & Mineral Resources (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种泥页岩储层含气量地球物理定量预测方法，包括以下步骤：1)获取目标工区泥页岩储层的地质参数，所述地质参数包括孔隙度和总有机碳含量；2)获取待测泥页岩储层的岩石物理参数；所述岩石物理参数包括纵波阻抗和密度；3)根据获取目标工区泥页岩储层的地质参数和岩石物理参数，建立泥页岩储层含气量预测模型；4)基于地震叠前同时反演，通过专业软件获取纵波阻抗数据体和密度数据体；5)将得到的纵波阻抗Ip数据体和密度ρ数据体代入预测模型中，即获得待测泥页岩储层的总含气量预测值。本发明预测结果符合地质事实，具有含气量定量预测结果准确、可靠的特点，适用性较为广泛，可以直接服务于页岩气勘探开发决策。
Description
技术领域
本发明涉及天然气勘探技术，尤其涉及一种泥页岩储层含气量地球物理定量预测方法。
背景技术
页岩气是近年来引起广泛关注度的一种非常规天然气资源，主要以游离态和吸附态赋存于暗色有机质泥页岩中。游离态页岩气主要分布于泥页岩储层的孔隙或裂隙之中，而吸附态页岩气则主要附着于有机质。泥页岩储层的含气量是页岩气选区评价和页岩气产能评估的重要参数，对页岩气勘探开发具有重要现实意义。
关于泥页岩储层含气量预测的方法，目前主要有原位解析测试法、等温吸附试验法、测井解释法、地震预测法等。其中，原位解析测试法虽然能够反映样品的实际含气量，但受取心方式和提钻时间影响，损失气估算误差较大。等温吸附试验法获得的是页岩的最大吸附气量，未考虑游离气量，与实际含气量有较大差别。测井解释法是通过测井数据结合岩心实验建立含气量的测井解释模型，可以预测页岩的含气量，但对于含气量较低的井误差较大，使得页岩含气量总体远小于实际含气量。而且，因客观实验测试数据分布范围有限，以上三种方法预测的含气量无法覆盖全区。地震预测法可以解决区域覆盖性问题，但常用于含气量预测的叠后波阻抗反演技术或地震波频谱衰减技术在泥页岩含气量预测方面精度往往较低，这是因为叠后波阻抗和地震频谱客观上会受到储层岩性、物性、厚度、地震波品质等因素综合影响。专利ZL201410359764.X公布的一种利用地震数据进行泥页岩含气量预测方法，即采用的是叠后反演手段。而且，现有公布的泥页岩含气量地震预测方法中普遍没有将游离气和吸附气区分开来，仅见正在实审中的申请号为201610796174.2的专利公布了一种针对页岩气游离气含量预测的方法。因此，现有方法获得的页岩气含量普遍面临着计算值偏小或误差较大或区域覆盖面不足等问题，为客观指导页岩气勘探开发实践带来了一定程度的不确定性。
发明内容
本发明要解决的技术问题在于针对现有技术中的缺陷，提供一种泥页岩储层含气量地球物理定量预测方法。
本发明解决其技术问题所采用的技术方案是：一种泥页岩储层含气量地球物理定量预测方法，包括以下步骤：
1)获取目标工区泥页岩储层的地质参数，所述地质参数包括用于确定游离气含量的孔隙度参数和用于确定吸附气含量的总有机碳含量参数；
2)获取待测泥页岩储层的岩石物理参数；所述岩石物理参数包括纵波阻抗和密度；
3)根据获取目标工区泥页岩储层的地质参数和岩石物理参数，建立泥页岩储层含气量预测模型，具体如下：
Qg＝0.8032×Φ+0.4563×TOC-1.7840；
Φ＝19.3350-0.001415×Ip；
TOC＝46.0233-16.7344×ρ；
其中，Qg为总含气量，Φ为孔隙度，TOC为总有机碳含量，Ip为纵波阻抗，ρ为密度；
4)基于地震叠前同时反演，通过专业软件获取纵波阻抗数据体和密度数据体；
5)将步骤4)得到的纵波阻抗Ip数据体和密度ρ数据体代入步骤3)所述的预测模型中，即获得待测工区泥页岩储层的总含气量预测值。
所述步骤4)中，地震叠前同时反演，通过专业软件获取纵波阻抗Ip数据体和密度ρ数据体是基于地质、测井和地震资料共同约束，在不同角道集叠加处理的基础上，选择合适的子波、低频模型、反演参数与流程，通过InterWell或Jason软件获取纵波阻抗数据体和密度数据体。
本发明产生的有益效果是：
1、针对现有技术中页岩气含量计算值偏小或误差较大或区域覆盖面不足，难以客观指导页岩气勘探开发，本发明分别从游离气、吸附气两个相态计算泥页岩储层中天然气的总量，符合泥页岩储层含气性特征，计算结果更加贴近地质事实。
2、由于原始地震资料的旅行时和振幅等信息保留更加丰富，叠前同时反演在含气性预测过程中可获得分辨率较高的岩石物理参数数据体，因而含气量预测结果精度也相应更高。
3、本发明在综合考虑游离气和吸附气的条件下，建立了总含气量与地质参数以及岩石物理参数之间的内在关联，并充分利用地震数据区域覆盖面广的特点，联合地震叠前同时反演，获取泥页岩储层总含气量，具有原始数据客观性强、参数关系内在逻辑强、地震预测方法先进、人为影响范围有限等特点，能够实现泥页岩储层总含气量的高精度预测，具有广泛的适用性。该预测结果可以直接服务于页岩气勘探开发决策，也可以为后续的页岩气生产提供高效勘探开发指导。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1为本发明的泥页岩储层含气量地球物理定量预测方法的流程图；
图2为实施例中泥页岩储层纵波阻抗数据体A井和B井的联井剖面示意图；
图3为实施例中泥页岩储层密度数据体A井和B井的联井剖面示意图；
图4为实施例中泥页岩储层的总含气量示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
如图1所示，一种泥页岩储层含气量地球物理定量预测方法，包括以下步骤：
1)获取待测泥页岩储层的地质参数；确定游离气含量的主控地质参数为孔隙度，吸附气含量的主控地质参数为总有机碳含量，然后利用回归分析建立总含气量(Qg)与孔隙度(Φ)和总有机碳含量(TOC)之间的计算关系：Qg＝0.8032×Φ+0.4563×TOC-1.7840。
2)获取待测泥页岩储层的岩石物理参数；确定与孔隙度、总有机碳含量相关的岩石物理参数分别为纵波阻抗、密度，然后利用回归分析建立孔隙度(Φ)与纵波阻抗(Ip)、总有机碳含量(TOC)与密度(ρ)之间的关系组：
Φ＝19.3350-0.001415×Ip；
TOC＝46.0233-16.7344×ρ。
3)根据步骤1)和2)，即得到纵波阻抗(Ip)和密度(ρ)表征总含气量(Qg)的预测模型：
Qg＝-7.636×ρ-0.0011365×Ip+34.7463。
4)基于地震叠前同时反演，通过专业软件InterWell或Jason获取纵波阻抗(Ip)数据体和密度(ρ)数据体。
地震叠前同时反演是本领域现有技术，可以基于地质、测井和地震资料共同约束，在不同角道集叠加处理的基础上，选择合适的子波、低频模型、反演参数与流程，通过InterWell或Jason等专业软件获取所需的敏感岩石物理参数数据体。
结果如图2和图3所示。
图2和图3显示A井和B井的纵波阻抗(Ip)数据体剖面、密度数据体(ρ)剖面，从图中可以看出，页岩层纵波阻抗(Ip)数据体和密度(ρ)数据体横向连续、稳定，纵向分辨率高。
5)将步骤4)得到的纵波阻抗(Ip)数据体和密度(ρ)数据体代入步骤3)所述的预测模型Qg＝-7.636×ρ-0.0011365×Ip+34.7463中，即获得待测工区泥页岩储层的总含气量，结果如图4所示。
从图4中可以看出，页岩气工区目的泥页岩储层含气量范围为1.91m3/t～4.75m3/t，A、B、C、D井控制的主体构造部位含气量普遍高于3.25m3/t，并存在明显内部分异特征。试气之后运用一点法获得的无阻流量结果显示，A、B、C、D四口井的无阻流量均大于17×104m3/d，普遍高于外围其它井，含气量预测结果与实际试气情况吻合度较高。
应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (2)
1.一种泥页岩储层含气量地球物理定量预测方法，其特征在于，包括以下步骤：
1)获取目标工区泥页岩储层的地质参数，所述地质参数包括用于确定游离气含量的孔隙度参数和用于确定吸附气含量的总有机碳含量参数；
2)获取待测泥页岩储层的岩石物理参数；所述岩石物理参数包括纵波阻抗和密度；
3)根据获取目标工区泥页岩储层的地质参数和岩石物理参数，建立泥页岩储层含气量预测模型，具体如下：
Qg＝0.8032×Φ+0.4563×TOC-1.7840；
Φ＝19.3350-0.001415×Ip；
TOC＝46.0233-16.7344×ρ；
其中，Qg为总含气量，Φ为孔隙度，TOC为总有机碳含量，Ip为纵波阻抗，ρ为密度；
4)基于地震叠前同时反演，通过专业软件获取待测工区泥页岩储层的纵波阻抗数据体和密度数据体；
5)将步骤4)得到的纵波阻抗Ip数据体和密度ρ数据体代入步骤3)所述的预测模型中，即获得待测工区泥页岩储层的总含气量预测值。
2.根据权利要求1所述的泥页岩储层含气量地球物理定量预测方法，其特征在于，所述步骤4)中，地震叠前同时反演，通过专业软件获取纵波阻抗数据体和密度数据体是基于地质、测井和地震资料共同约束，在不同角道集叠加处理的基础上，选择合适的子波、低频模型、反演参数与流程，通过InterWell或Jason软件获取纵波阻抗数据体和密度数据体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910823702.2A CN110579797A (zh) | 2019-09-02 | 2019-09-02 | 一种泥页岩储层含气量地球物理定量预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910823702.2A CN110579797A (zh) | 2019-09-02 | 2019-09-02 | 一种泥页岩储层含气量地球物理定量预测方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110579797A true CN110579797A (zh) | 2019-12-17 |
Family
ID=68812421
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910823702.2A Pending CN110579797A (zh) | 2019-09-02 | 2019-09-02 | 一种泥页岩储层含气量地球物理定量预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110579797A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113534263A (zh) * | 2021-07-13 | 2021-10-22 | 广州海洋地质调查局 | 一种不依赖测井资料的含油气饱和度预测方法 |
| CN113960659A (zh) * | 2021-10-14 | 2022-01-21 | 中国矿业大学 | 一种地震岩石物理驱动的煤层气储层含气量预测方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110068788A1 (en) * | 2009-09-23 | 2011-03-24 | Chanh Cao Minh | Estimating porosity and fluid volume |
| WO2013023011A2 (en) * | 2011-08-10 | 2013-02-14 | Schlumberger Canada Limited | Logging in gas shale and other unconventional reservoirs |
| CN104573344A (zh) * | 2014-12-25 | 2015-04-29 | 中国海洋石油总公司 | 一种通过测井数据获取页岩储层含气量的方法 |
| CN105927218A (zh) * | 2016-05-20 | 2016-09-07 | 中国石油大学(北京) | 一种陆相泥页岩储层含气量预测方法及装置 |
| CN106251232A (zh) * | 2016-08-12 | 2016-12-21 | 中国石油天然气股份有限公司 | 确定页岩含气量的方法和装置 |
| WO2018017930A1 (en) * | 2016-07-22 | 2018-01-25 | Gas Sensing Technology Corp. | In situ evaluation of gases and liquids low permeability reservoirs |
| CN108107182A (zh) * | 2016-11-24 | 2018-06-01 | 中国地质大学（北京） | 一种页岩气中游离气和吸附气含量判定方法 |
-
2019
- 2019-09-02 CN CN201910823702.2A patent/CN110579797A/zh active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110068788A1 (en) * | 2009-09-23 | 2011-03-24 | Chanh Cao Minh | Estimating porosity and fluid volume |
| WO2013023011A2 (en) * | 2011-08-10 | 2013-02-14 | Schlumberger Canada Limited | Logging in gas shale and other unconventional reservoirs |
| CN104573344A (zh) * | 2014-12-25 | 2015-04-29 | 中国海洋石油总公司 | 一种通过测井数据获取页岩储层含气量的方法 |
| CN105927218A (zh) * | 2016-05-20 | 2016-09-07 | 中国石油大学(北京) | 一种陆相泥页岩储层含气量预测方法及装置 |
| WO2018017930A1 (en) * | 2016-07-22 | 2018-01-25 | Gas Sensing Technology Corp. | In situ evaluation of gases and liquids low permeability reservoirs |
| CN106251232A (zh) * | 2016-08-12 | 2016-12-21 | 中国石油天然气股份有限公司 | 确定页岩含气量的方法和装置 |
| CN108107182A (zh) * | 2016-11-24 | 2018-06-01 | 中国地质大学（北京） | 一种页岩气中游离气和吸附气含量判定方法 |
Non-Patent Citations (1)
| Title |
|---|
| 罗明生 等: ""页岩气‘甜点’地震预测技术在威远地区的应用", 《2015年全国天然气学术年会》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113534263A (zh) * | 2021-07-13 | 2021-10-22 | 广州海洋地质调查局 | 一种不依赖测井资料的含油气饱和度预测方法 |
| CN113960659A (zh) * | 2021-10-14 | 2022-01-21 | 中国矿业大学 | 一种地震岩石物理驱动的煤层气储层含气量预测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
| WO2019062655A1 (zh) | 薄夹层的确定方法和装置 | |
| CN110318745B (zh) | 一种沉积微相约束下粒径岩性测井评价方法 | |
| CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
| Hunt et al. | Quantitative estimate of fracture density variations in the Nordegg with azimuthal AVO and curvature: A case study | |
| CN113837501A (zh) | 基于图神经网络的页岩总有机碳预测方法、系统和设备 | |
| CN116168224A (zh) | 基于成像砾石含量的机器学习岩相自动识别方法 | |
| CN110579797A (zh) | 一种泥页岩储层含气量地球物理定量预测方法 | |
| CN107817260A (zh) | 一种泥页岩高频层序识别方法 | |
| CN110529106B (zh) | 一种利用测井资料确定煤层显微组分含量的方法 | |
| CN113216945B (zh) | 一种致密砂岩储层渗透性定量评价方法 | |
| CN112946782B (zh) | 一种致密油气储渗体地震精细刻画方法 | |
| CN112185469B (zh) | 一种预测海域天然气水合物有利聚集区的方法 | |
| CN106568918B (zh) | 页岩有机碳含量toc预测方法 | |
| AU2017279838B1 (en) | Method for classifying deep rock geofacies based on data mining | |
| CN112505754A (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
| CN110118994B (zh) | 一种基于地震反演和机器学习的陆相烃源岩定量预测方法 | |
| CN113484907B (zh) | 一种预测不同类型储层平面上分布的方法 | |
| CN113806674B (zh) | 古河道纵向尺度的量化方法、装置、电子设备及存储介质 | |
| CN110454155B (zh) | 一种应用磁化率测井方法确定第四纪地层年代的方法 | |
| CN112394392B (zh) | 对烃源岩分布情况进行评价的方法和装置 | |
| CN112180443B (zh) | 页岩气二维地震甜点区优选方法及装置 | |
| CN113126155A (zh) | 一种针对分布于煤岩间受强反射影响的砂岩储层预测方法 | |
| CN114609666B (zh) | 一种页岩薄储层的预测方法、装置、设备及存储介质 | |
| Amanipoor | Providing a subsurface reservoir quality maps in oil fields by geostatistical methods |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191217 |
|
| RJ01 | Rejection of invention patent application after publication |