CN112709566B - 一种碳酸盐岩气藏产层的预测方法 - Google Patents
一种碳酸盐岩气藏产层的预测方法 Download PDFInfo
- Publication number
- CN112709566B CN112709566B CN201911025173.8A CN201911025173A CN112709566B CN 112709566 B CN112709566 B CN 112709566B CN 201911025173 A CN201911025173 A CN 201911025173A CN 112709566 B CN112709566 B CN 112709566B
- Authority
- CN
- China
- Prior art keywords
- well
- predicted
- time difference
- fitting
- lateral resistivity
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 89
- 238000000034 method Methods 0.000 title claims abstract description 50
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 44
- 238000005070 sampling Methods 0.000 claims abstract description 66
- 238000012216 screening Methods 0.000 claims abstract description 4
- 230000014509 gene expression Effects 0.000 claims description 33
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000005856 abnormality Effects 0.000 claims description 6
- 238000011161 development Methods 0.000 abstract description 7
- 239000007789 gas Substances 0.000 description 21
- 230000004044 response Effects 0.000 description 8
- 230000000694 effects Effects 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 150000002430 hydrocarbons Chemical class 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- -1 carbonate hydrocarbon Chemical class 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004062 sedimentation 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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 And Detection Of Objects (AREA)
Abstract
本发明涉及一种碳酸盐岩气藏产层的预测方法，属于测井技术领域。该预测方法包括：建立目标单井所在区域的井段中补偿中子和拟合深侧向电阻率之间关系、补偿中子和拟合声波时差之间关系的表现形式；确定待预测井段各深度采样点的补偿中子和拟合深侧向电阻率的关系、补偿中子和拟合声波时差的关系；计算各深度采样点的拟合深侧向电阻率值和拟合声波时差值；比较待预测井段储层段中各深度采样点的拟合深侧向电阻率值与深侧向电阻率值大小、拟合声波时差值和声波时差值大小，根据比较结果对产层进行筛选。本发明的预测方法，能够相对准确、可靠的利用裸眼井常规测井资料预测出储层段能否产出，为碳酸盐岩气藏制定合理高效的开发技术对策提供依据。
Description
技术领域
本发明涉及一种碳酸盐岩气藏产层的预测方法，属于测井技术领域。
背景技术
测井是在勘探和开发油气、煤及金属矿体的过程中，利用各种仪器测量井下岩层的物理参数及井的技术状况，分析所记录的资料，进行地质和工程方面的研究。
在石油钻井时，在钻到设计井深深度后都必须进行完井电测，以获得各种石油地质及工程技术资料，作为完井和开发油田的原始资料，这种测井习惯上称为裸眼常规测井，裸眼井常规测井测量的主要目的是获取各类油气藏的岩性、物性、含油气性等关键技术参数。
在油气井下完套管后，从其投产后到报废整个生产过程中，利用各种测井仪器进行井下测试以获取相应地下信息的测井称为生产测井，生产测井的主要目的是测量测井油气井生产过程中流动剖面，掌握各油气层射孔层段产出或吸入流体的性质和流量，以便对油气井产状和油层开采特征做出评价。
套管井中的产出剖面测井、饱和度测井等生产测井是现代测井技术中识别产层的重要技术手段，利用生产井中温度、压力、涡轮转速等参数的变化能够相对准确的解释出产层、非产层和产层的产量，但因受生产投资成本的限制，套管井中的产出剖面测井、饱和度测井等生产测井只对一部分井进行过测井，其余的井的产出状况只能靠预测，因此仅仅基于产剖测井资解释的产层段为油气藏剩余油气的挖潜和开发方案的制定等方面的研究提供的基础资料是不全面的，而利用裸眼井常规测井资料预测碳酸盐岩气藏储层能否产出是支撑上述这些研究的全面开展是非常重要的技术手段。
对于非均质性较弱的油气藏，裸眼井常规测井解释获得的储层的孔隙度、渗透率、含油气饱和度等参数值的大小和生产测井解释储层段的能否产出往往是对应的，具体表现为在同一个射孔层段内，高孔隙度、高渗透率的油气储层段是产出层段的特征；然而碳酸盐岩气藏储层段由于受沉积作用、成岩作用、构造作用的共同作用，储层的孔隙类型、溶蚀程度、裂缝发育程度复杂多样，储层的非均质性较强，导致产层与非产层的孔隙的大小对应性差、储层段能否产出层间差异较大，利用裸眼井常规测井资料解释的孔隙度等物性参数表征产层与非产层非常困难。
现有技术利用裸眼井常规测井资料预测碳酸盐岩气藏产层的方法只处于利用简单的常规测井响应特征分析阶段，因碳酸盐岩气藏的产层往往表现为声波时差增大、电阻率减小、补偿中子减小等特征，利用这些特征的存在与否可预示着储层段能否产出，但这些参数增大或减小等特征无固定、统一的标准，使得这些分析结果似是而非，利用裸眼井常规测井资料实现预测碳酸盐岩气藏产层的方法的研究目前尚属空白对碳酸盐岩气藏合理有效开发带来障碍。
发明内容
本发明的目的是提供一种碳酸盐岩气藏产层的预测方法，能够相对准确、可靠地预测出储层段中的产层。
为了实现以上目的，本发明所采用的技术方案是：
一种碳酸盐岩气藏产层的预测方法，包括以下步骤：
1)对目标单井所在区域内若干同时具有常规测井资料和生产测井资料的单井进行分析，建立目标单井所在区域的井段中补偿中子和拟合深侧向电阻率之间关系的表现形式、补偿中子和拟合声波时差之间关系的表现形式；
2)利用目标单井待预测井段各深度采样点的补偿中子、深侧向电阻率和声波时差，采用数据拟合的方法确定待预测井段各深度采样点的两表现形式中待定参数，进而得到目标单井待预测井段各深度采样点的补偿中子和拟合深侧向电阻率的关系、补偿中子和拟合声波时差的关系；
3)获取目标单井待预测井段各深度采样点的补偿中子值、深侧向电阻率值和声波时差值分别代入对应深度采样点的补偿中子和拟合深侧向电阻率的关系、补偿中子和拟合声波时差的关系，得到各深度采样点的拟合深侧向电阻率值和拟合声波时差值；
4)比较待预测井段的储层段中各深度采样点的拟合深侧向电阻率值与深侧向电阻率值之间大小、拟合声波时差值和声波时差值之间的大小，筛选出待预测井段中拟合深侧向电阻率值大于深侧向电阻率值的井段，若该筛选出的井段中存在拟合声波时差值小于声波时差值的异常，则认为该筛选出的井段能够产出，为预测的产层；若不存在所述异常，则筛选出的井段为非产层。
本发明的碳酸盐岩气藏产层的预测方法，通过分析目标单井待预测井段中储层段的拟合深侧向电阻率与深侧向电阻率之间、拟合声波时差与声波时差之间的差异特征来预测储层段中的产层，该碳酸盐岩气藏产层的预测方法能够相对准确、可靠的利用裸眼井常规测井资料预测出储层段能否产出，为碳酸盐岩气藏制定合理高效的开发技术对策提供了依据。
上述碳酸盐气藏产层的预测方法，还可以将步骤2)替换为：获取目标单井待预测井段各深度采样点的补偿中子值、深侧向电阻率值和声波时差值，利用获取到的数据采用数据拟合的方法确定待预测井段各深度采样点的两表现形式中待定参数的值，进而得到目标单井待预测井段各深度采样点的补偿中子和拟合深侧向电阻率的数学关系、补偿中子和拟合声波时差的数学关系；将步骤3)替换为：将目标单井待预测井段各深度采样点的补偿中子值代入对应深度采样点的补偿中子和拟合深侧向电阻率的数学关系、补偿中子和拟合声波时差的数学关系中，得到各深度采样点的拟合深侧向电阻率值和拟合声波时差值。
步骤4)中，待预测井段的储层段为目标单井裸眼井常规测井资料解释的在目标单井待预测井段中的储层段。
步骤1)中，所述生产测井资料为套管井产出剖面测井资料。
步骤1)中，先确定目标单井所在区域的储层段中非产层补偿中子和深侧向电阻率的关系、补偿中子和声波时差的关系，再建立目标单井所在区域井段中补偿中子和拟合深侧向电阻率之间关系、补偿中子和拟合声波时差之间关系的表现形式；储层段中非产层补偿中子和深侧向电阻率的关系、补偿中子和声波时差的关系均通过交会图确定；储层段中非产层补偿中子和深侧向电阻率的关系通过储层段中非产层补偿中子和深侧向电阻率的交会图确定；储层段中非产层补偿中子和声波时差的关系通过储层段中非产层补偿中子和声波时差的交会图确定。
经过对目标单井所在区域的储层段中产层与非产层在补偿中子、深侧向电阻率、声波时差测井响应特征上的差异进行分析，将目标单井所在区域的储层段中非产层补偿中子和深侧向电阻率之间关系的表现形式作为目标单井所在区域的井段补偿中子和拟合深侧向电阻率之间关系的表现形式，将目标单井所在区域的储层段中非产层补偿中子和声波时差之间关系的表现形式作为目标单井所在区域的井段补偿中子和拟合声波时差之间关系的表现形式；上述目标单井所在区域的储层段是对同时具有常规测井资料和生产测井资料的单井进行分析确定的。
补偿中子和拟合深侧向电阻率之间关系的表现形式、补偿中子和拟合声波时差之间关系的表现形式均为含有待定参数的关系表达式。
步骤2)中的两表现形式是指步骤1)中确定的补偿中子和拟合深侧向电阻率之间关系的表现形式、补偿中子和拟合声波时差之间关系的表现形式。
补偿中子和拟合深侧向电阻率之间关系的表现形式中的待定参数、补偿中子和拟合声波时差之间关系的表现形式中的待定参数均可通过数据拟合的方法确定，如采用最小二乘法进行数据拟合确定。
目标单井待预测井段补偿中子和拟合深侧向电阻率之间关系的表现形式为：lgRdc＝a·ΦN+b；式中，Rdc为拟合深侧向电阻率，ΦN为补偿中子，a、b为待定参数。
待定参数a的计算方法为：
式中，ΦNi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的补偿中子值，Rdi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的深侧向电阻率值，m为固定的迭代窗长，n为固定的迭代步长。
待定参数b的计算方法为：
式中，ΦNi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的补偿中子值，Rdi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的深侧向电阻率值，m为固定的迭代窗长，n为固定的迭代步长。
目标单井待预测井段补偿中子和拟合声波时差之间关系的表现形式为：Δtc＝c·ΦN+d；式中，Δtc为拟合声波时差，ΦN为补偿中子，c、d为待定参数。
待定参数c的计算方法为：
式中，ΦNi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的补偿中子值，Δti为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的声波时差值，m为固定的迭代窗长，n为固定的迭代步长。
待定参数d的计算方法为：
式中，ΦNi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的补偿中子值，Δti为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的声波时差值，m为固定的迭代窗长，n为固定的迭代步长。
固定迭代窗长m是指待预测井段中的一个固定长度为m的滑动窗口，滑动窗口是滑窗迭代算法中的一个术语。
固定迭代步长n是指每次固定窗长m向下移动长度为n的距离。
附图说明
图1为实施例中碳酸盐岩气藏产层的预测方法的流程图；
图2为实施例中A气田M4井常规测井组合成果、生产测井解释成果和产层段预测效果图。
具体实施方式
以下结合中国A气田的碳酸盐岩气藏对本发明的技术方案作进一步的说明。
实施例1
本实施例的碳酸盐岩气藏产层的预测方法，流程图见图1，包括以下步骤：
1)取目标单井所在区域内若干同时具有裸眼井常规测井资料和生产测井资料的单井；
利用所取单井的裸眼井常规测井资料和生产测井资料，分析各所取单井裸眼井常规测井解释的储层段中产层与非产层在补偿中子、深侧向电阻率和声波时差测井响应特征上的差异，确定目标单井所在区域储层段的非产层中补偿中子和深侧向电阻率的关系类型、补偿中子和声波时差的关系类型，依据所确定的关系类型对应得到目标单井所在区域储层段中非产层补偿中子和深侧向电阻率之间关系的表现形式、补偿中子和声波时差之间的关系的表现形式，作为目标单井待预测井段中补偿中子和拟合深侧向电阻率之间关系的表现形式、补偿中子和拟合声波时差之间关系的表现形式：
在生产测井解释的碳酸盐岩气藏的产层中，往往具有饱含气、有裂缝或溶蚀等沟通作用、孔隙连通性好等特征，产层往往是沿着优势渗流通道产出的储层，产层相对于非产层，补偿中子因储层含气的挖掘效应导致其测井响应特征值更小，声波时差因储层含气的导致其测井响应特征值更大，裂缝或溶蚀的沟通作用因钻井液泥浆侵入的影响导致产层的深侧向电阻率的测井响应特征值更低，同时当储层中裂缝为低角度缝或斜交缝时声波时差表现为增大的特征，通过对A气田碳酸盐岩气藏的12口共同具有裸眼井常规测井资料、解释成果和套管井产出剖面测井解释成果的井资料分析表明，每口井在非产层中补偿中子和深侧向电阻率、补偿中子和声波时差的交会图中，非产层段补偿中子和深侧向电阻率的交会关系都表现为较好的对数关系，补偿中子和声波时差的交会关系都表现为较好的正线性关系，而在产层段中因补偿中子和深侧向电阻率的减小、声波时差的增大，产层段的相关特征和非产层的差异明显，因此可建立A气田碳酸盐岩气藏单井储层段中非产层补偿中子和深侧向电阻率之间的交会关系表达式、补偿中子和声波时差之间的交会关系表达式，补偿中子和深侧向电阻率的交会关系表达式为：
lg Rd＝a·ΦN+b (1)
补偿中子和声波时差的交会关系表达式为：
Δt＝c·ΦN+d (2)
式(1)、(2)中，Rd为深侧向电阻率，单位为Ω·m，Δt为声波时差，单位为μs/ft，ΦN为补偿中子，单位为P.U，a、b、c、d均为待定参数，式(1)和(2)中的关系在A气田碳酸盐岩气藏中每口井中都是适用的。
在A气田碳酸盐岩气藏的储层中，每口井在不同深度处的的非产层段中式(1)和式(2)中的待定参数a、b、c、d是变化的，但补偿中子和深侧向电阻率之间关系的表现形式、补偿中子和声波时差之间关系的表现形式是不变的，无论在非产层和产层段都可以利用式(1)和式(2)的关系确定出拟合深侧向电阻率和拟合声波时差，只是在非产层段深侧向电阻率和拟合深侧向电阻率、声波时差和拟合声波时差是相同的，因此在目标单井待预测井段中，补偿中子和拟合深侧向电阻率之间关系的表现形式为：
lg Rdc＝a·ΦN+b (3)
补偿中子和拟合声波时差之间关系的表现形式分别为：
Δtc＝c·ΦN+d (4)
式(3)(4)中，Rdc为拟合深侧向电阻率，单位为Ω·m，Δtc为拟合声波时差，单位为μs/ft，式(3)和(4)中的关系在A气田碳酸盐岩气藏中每口井中同样都是适用的。
2)利用目标单井待预测井段中各深度采样点对应的补偿中子、深侧向电阻率和声波时差，采用最小二乘法进行拟合得到待预测井段各深度采样点的补偿中子和拟合深侧向电阻率之间关系的表现形式、补偿中子和拟合声波时差之间关系的表现形式中的待定参数，进而得到待预测井段各深度采样点的补偿中子和拟合深侧向电阻率之间的关系、补偿中子和拟合声波时差之间的关系：
本实施例选取中国A气田M4井来进行产层预测处理方法的说明，M4井为中国A气田的一口碳酸盐岩气藏生产井，M4井的3850～3912m段常规测井组合成果、生产测井解释成果和产层段预测效果图如图2所示，图2中的第三道中的虚线、第四道中的虚线、第五道中的虚线分别为M4井的补偿中子ΦN、深侧向电阻率Rd和声波时差Δt测井曲线，第七道为裸眼井常规测井解释的储层段，第八道为生产测井解释的产层段，因M4井具有裸眼井常规测井资料、解释成果和套管井产出剖面测井解释成果，可以用于预测产层的验证。因此，为了更好地证明本发明的预测方法的准确性，实施例中选取M4井3850～3912米井段为待预测井段。
利用M4井的各深度采样点的补偿中子ΦN、深侧向电阻率Rd和声波时差Δt，采用最小二乘法确定式(3)和式(4)的中的待定参数；在待预测井段中，选取长度m为固定的迭代窗长(单位为米)，长度n为固定的迭代步长(单位为米)，待定参数a、b、c、d的表达式分别为：
式(5)～(8)中，ΦNi为第i个深度采样点的补偿中子值，单位为P.U，Rdi为第i个深度采样点的深侧向电阻率值，单位为Ω·m，Δti为第i个深度采样点的声波时差值，单位为μs/ft。
由(3)、(4)、(5)、(6)、(7)、(8)即可确定待预测井段各深度采样点的拟合深侧向电阻率和补偿中子之间的关系、拟合声波时差和补偿中子之间的关系。
3)根据M4井的测井资料获取该井待预测井段各深度采样点的的补偿中子ΦN(单位为P.U)、深侧向电阻率Rd(单位为Ω·m)和声波时差Δt(单位为μs/ft)测井响应特征值，将获取得到的数据分别代入待预测井段各深度采样点的补偿中子和拟合深侧向电阻率之间的关系、补偿中子和拟合声波时差之间的关系，得到各深度采样点拟合深侧向电阻率值和拟合声波时差值；利用待预测井段的各深度采样点拟合深侧向电阻率值和拟合声波时差值可以得到待预测井段的拟合深侧向电阻率曲线和拟合声波时差曲线，图2中第四道的实线和第五道的实线分别为确定的拟合深侧向电阻率曲线和拟合声波时差曲线。
4)比较单井待预测井段中储层段的拟合深侧向电阻率与深侧向电阻率之间、拟合声波时差和声波时差之间的大小，筛选出储层段中拟合深侧向电阻率大于深侧向电阻率的井段，若该筛选出的井段中存在拟合声波时差小于声波时差的异常，则认为该筛选出的井段能够产出，为预测的产层；目标单井待预测井段的储层段为目标单井裸眼井常规测井资料解释的在目标单井待预测井段中的储层段：
本实施例在M4井待预测井段3850～3912米通过分析各储层段中拟合深侧向电阻率与深侧向电阻率之间、拟合声波时差和声波时差之间的大小来预测产层，由于产层段的深侧向电阻率因储层连通性增强、钻井液的侵入等表现出减小的特征，拟合深侧向电阻率大于深侧向电阻率曲线，同时产层段的声波时差由于裂缝、含气等因素导致的增大特征，局部产层段表现拟合声波时差小于声波时差曲线，因此，可综合依据拟合深侧向电阻率与深侧向电阻率之间、拟合声波时差和声波时差之间的差异预测储层段能否产出，图2中第九道即为利用本发明的方法预测的M4井产层段，图2中第八道为M4井的生产测井产出剖面测井解释出的产层段，通过对比分析表明，本发明方法利用裸眼井常规测井预测的碳酸盐岩储层中的产层段与产出剖面解释出的产层段能够很好的对应，证明了该碳酸盐岩气藏产层的预测方法是可行、有效的，能够为碳酸盐岩气藏制定合理高效的开发技术对策提供技术支撑。
作为本发明的实施例2，目标单井所在区域内同时具有裸眼井常规测井资料和生产测井资料的单井的数量较少时，在目标单井所在区域内取的单井的数量也可以为一口。
作为本发明的实施例3，还可以通过其他数据拟合的方法确定各深度采样点的拟合深侧向电阻率和补偿中子之间关系、拟合声波时差和补偿中子之间关系的待定参数a、b、c、d的值。
作为本发明的实施例4，还可以将实施例1中由步骤3)可确定出单井待预测井段中的拟合深侧向电阻率曲线和拟合声波时差曲线，可通过分析单井待预测井段中储层段的拟合深侧向电阻率与深侧向电阻率之间、拟合声波时差和声波时差之间的差异特征来预测储层段能否产出，筛选出储层段中拟合深侧向电阻率大于深侧向电阻率的井段，若该筛选出的井段中存在拟合声波时差小于声波时差的异常，则认为该筛选出的井段能够产出，为预测的产层；否则，认为该筛选出的井段不能够产出，为非产层。
作为本发明的实施例5，还以在实施例1的步骤3)和步骤4)中，直接根据M4井的测井资料获取该井待预测井段各深度采样点的补偿中子ΦN、深侧向电阻率Rd和声波时差Δt测井响应特征值(即补偿中子值、深侧向电阻率值和声波时差值)，采用最小二乘法确定式(3)和式(4)的中的待定参数；将待预测井段各深度采样点的补偿中子ΦN(单位为P.U)、深侧向电阻率Rd(单位为Ω·m)和声波时差Δt(单位为μs/ft)测井响应特征值带入式(5)、(6)、(7)、(8)中确定出各深度采样点处待定参数a、b、c、d的值，然后得到的待预测井段各深度采样点的拟合深侧向电阻率和补偿中子之间的数学关系、拟合声波时差和补偿中子之间数学关系，将目标单井待预测井段各深度采样点的补偿中子值代入对应深度采样点的补偿中子和拟合深侧向电阻率的数学关系、补偿中子和拟合声波时差的数学关系中，得到各深度采样点的拟合深侧向电阻率值和拟合声波时差值。
作为本发明的实施例6，还可以将实施例1中确定目标单井所在区域储层段的非产层中补偿中子和深侧向电阻率的关系类型、补偿中子和声波时差的关系类型的方法替换为其他常用的方法。
Claims (3)
1.一种碳酸盐岩气藏产层的预测方法，其特征在于：包括以下步骤：
1)对目标单井所在区域内若干同时具有常规测井资料和生产测井资料的单井进行分析，建立目标单井所在区域的井段中补偿中子和拟合深侧向电阻率之间关系的表现形式、补偿中子和拟合声波时差之间关系的表现形式；
2)利用目标单井待预测井段各深度采样点的补偿中子、深侧向电阻率和声波时差，采用数据拟合的方法确定待预测井段各深度采样点的两表现形式中待定参数，进而得到目标单井待预测井段各深度采样点的补偿中子和拟合深侧向电阻率的关系、补偿中子和拟合声波时差的关系；
3)获取目标单井待预测井段各深度采样点的补偿中子值、深侧向电阻率值和声波时差值分别代入对应深度采样点的补偿中子和拟合深侧向电阻率的关系、补偿中子和拟合声波时差的关系，得到各深度采样点的拟合深侧向电阻率值和拟合声波时差值；
4)比较待预测井段的储层段中各深度采样点的拟合深侧向电阻率值与深侧向电阻率值之间大小、拟合声波时差值和声波时差值之间的大小，筛选出待预测井段中拟合深侧向电阻率值大于深侧向电阻率值的井段，若该筛选出的井段中存在拟合声波时差值小于声波时差值的异常，则认为该筛选出的井段能够产出，为预测的产层；若不存在所述异常，则筛选出的井段为非产层；
所述目标单井待预测井段补偿中子和拟合深侧向电阻率之间关系的表现形式为：lgRdc＝a·ΦN+b；式中，Rdc为拟合深侧向电阻率，ΦN为补偿中子，a、b为待定参数；
待定参数a的计算方法为：
式中，ΦNi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的补偿中子值，Rdi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的深侧向电阻率值，m为固定的迭代窗长，n为固定的迭代步长；
待定参数b的计算方法为：
式中，ΦNi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的补偿中子值，Rdi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的深侧向电阻率值，m为固定的迭代窗长，n为固定的迭代步长；
所述目标单井待预测井段补偿中子和拟合声波时差之间关系的表现形式为：Δtc＝c·ΦN+d；式中，Δtc为拟合声波时差，ΦN为补偿中子，c、d为待定参数；
待定参数c的计算方法为：
式中，ΦNi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的补偿中子值，Δti为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的声波时差值，m为固定的迭代窗长，n为固定的迭代步长；
待定参数d的计算方法为：
式中，ΦNi为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的补偿中子值，Δti为碳酸盐岩气藏目标单井待预测井段中第i个深度采样点对应的声波时差值，m为固定的迭代窗长，n为固定的迭代步长。
2.根据权利要求1所述的碳酸盐岩气藏产层的预测方法，其特征在于：
将步骤2)替换为：获取目标单井待预测井段各深度采样点的补偿中子值、深侧向电阻率值和声波时差值，利用获取到的数据采用数据拟合的方法确定待预测井段各深度采样点的两表现形式中待定参数的值，进而得到目标单井待预测井段各深度采样点的补偿中子和拟合深侧向电阻率的数学关系、补偿中子和拟合声波时差的数学关系；
将步骤3)替换为：将目标单井待预测井段各深度采样点的补偿中子值代入对应深度采样点的补偿中子和拟合深侧向电阻率的数学关系、补偿中子和拟合声波时差的数学关系中，得到各深度采样点的拟合深侧向电阻率值和拟合声波时差值。
3.根据权利要求1或2所述的碳酸盐岩气藏产层的预测方法，其特征在于：
步骤4)中，待预测井段的储层段为目标单井裸眼井常规测井资料解释的在目标单井待预测井段中的储层段。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911025173.8A CN112709566B (zh) | 2019-10-25 | 2019-10-25 | 一种碳酸盐岩气藏产层的预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911025173.8A CN112709566B (zh) | 2019-10-25 | 2019-10-25 | 一种碳酸盐岩气藏产层的预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112709566A CN112709566A (zh) | 2021-04-27 |
| CN112709566B true CN112709566B (zh) | 2024-03-19 |
Family
ID=75540836
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911025173.8A Active CN112709566B (zh) | 2019-10-25 | 2019-10-25 | 一种碳酸盐岩气藏产层的预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112709566B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104747163A (zh) * | 2013-12-31 | 2015-07-01 | 中国石油天然气股份有限公司 | 一种在致密砂岩中识别储层裂缝的方法及装置 |
| CN105221133A (zh) * | 2015-08-20 | 2016-01-06 | 中国石油化工股份有限公司 | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 |
| WO2016161914A1 (zh) * | 2015-04-07 | 2016-10-13 | 四川行之智汇知识产权运营有限公司 | 一种利用地质和测井信息预测储层成岩相的方法 |
| CN107795320A (zh) * | 2016-09-05 | 2018-03-13 | 中国石油化工股份有限公司 | 一种水平井碳酸盐岩储层参数的计算方法 |
| CN109425912A (zh) * | 2017-09-05 | 2019-03-05 | 中国石油化工股份有限公司 | 一种碳酸盐岩油藏储层有效裂缝的识别方法 |
| CN109931054A (zh) * | 2018-12-27 | 2019-06-25 | 西南石油大学 | 致密砂岩储层压力的预测方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020177954A1 (en) * | 1994-03-17 | 2002-11-28 | Vail William Banning | Processing formation resistivity measurements obtained from within a cased well used to quantitatively determine the amount of oil and gas present |
-
2019
- 2019-10-25 CN CN201911025173.8A patent/CN112709566B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104747163A (zh) * | 2013-12-31 | 2015-07-01 | 中国石油天然气股份有限公司 | 一种在致密砂岩中识别储层裂缝的方法及装置 |
| WO2016161914A1 (zh) * | 2015-04-07 | 2016-10-13 | 四川行之智汇知识产权运营有限公司 | 一种利用地质和测井信息预测储层成岩相的方法 |
| CN105221133A (zh) * | 2015-08-20 | 2016-01-06 | 中国石油化工股份有限公司 | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 |
| CN107795320A (zh) * | 2016-09-05 | 2018-03-13 | 中国石油化工股份有限公司 | 一种水平井碳酸盐岩储层参数的计算方法 |
| CN109425912A (zh) * | 2017-09-05 | 2019-03-05 | 中国石油化工股份有限公司 | 一种碳酸盐岩油藏储层有效裂缝的识别方法 |
| CN109931054A (zh) * | 2018-12-27 | 2019-06-25 | 西南石油大学 | 致密砂岩储层压力的预测方法 |
Non-Patent Citations (3)
| Title |
|---|
| 五百梯气田石炭系储层裂缝测井识别研究;周平;曹子颜;崔恒远;王亦然;;河北地质大学学报(04);全文 * |
| 低电阻率油气层测井识别技术;赵佐安, 何绪全, 唐雪萍;天然气工业(04);全文 * |
| 塔河地区碳酸盐岩裂缝型储层的测井评价技术;陈冬;魏修成;;石油物探(02);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112709566A (zh) | 2021-04-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108713089B (zh) | 基于钻孔流体和钻探录井估计地层性质 | |
| CN106761677B (zh) | 页岩气水平井单井产能的测井预测方法 | |
| CN104514552A (zh) | 一种煤层气储层识别与丰度预测的方法 | |
| US9366776B2 (en) | Integrated formation modeling systems and methods | |
| CN109543915B (zh) | 基于测井信息对全井段烃源岩总有机碳含量的识别方法 | |
| CN113419284B (zh) | 一种基于聚类分析测井岩石物理相双甜点识别方法 | |
| CN103821495B (zh) | 测井方法 | |
| CN112709566B (zh) | 一种碳酸盐岩气藏产层的预测方法 | |
| CN112560246A (zh) | 一种目标井散点地层压力系数的预测方法 | |
| CHAFEET et al. | STUDY OF PETROPHYSICAL PROPERTIES OF MISHRIF AND YAMAMA FORMATIONS AT SELECTED FIELDS, SOUTH IRAQ. | |
| CN110017136B (zh) | 一种基于视水层电阻率的水淹层识别与产水率预测方法 | |
| CN111894571B (zh) | 一种基于岩性扫描测井资料的流体性质识别方法 | |
| Chong et al. | Application of machine learning to characterize gas hydrate reservoirs in Mackenzie Delta (Canada) and on the Alaska north slope (USA) | |
| CN113006779B (zh) | 一种砾岩地层岩石可钻性级值的评价方法和校正方法 | |
| CN108732621B (zh) | 一种基于ffc-电阻率的随钻精细时深预测方法 | |
| CN108756867B (zh) | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 | |
| Leary et al. | Crosswell seismic applications to highly heterogeneous tight gas reservoirs | |
| Slim et al. | Quantification of Downhole Production Profiles Using Production, Acoustic, and Temperature Logging: A Case Study from a High-GOR Producer to Design an Optimized Production Strategy | |
| Aissou et al. | Variations in the present-day stress field and its implications for hydrocarbon development in the Hassi Terfa field, Hassi Messaoud, Algeria | |
| Nwosu et al. | Computer processed interpretation of geophysical logs of an oil field Niger Delta sedimentary basin, onshore, Nigeria | |
| CN114895364B (zh) | 基于温压耦合电阻率约束的深部地温场预测方法及装置 | |
| CN114895365B (zh) | 高温高压岩心电阻率标定的深部地温场预测方法及装置 | |
| Habeeb | Introduction and Investigation into Oil Well Logging Operations | |
| Rezaeeparto et al. | Prediction of Total Organic Carbon (TOC) Utilizing ΔlogR and Artificial Neural Network (ANN) Methods and Geochemical Facies Determination of Kazhdumi Formation in One of the Fields-Southwest of Iran | |
| Elkhateeb | Integrated Formation Evaluation Study in Cliff Head Field, Northern Perth Basin, Western Australia |
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 |