CN105370272A - 页岩气储层的测井评价方法 - Google Patents
页岩气储层的测井评价方法 Download PDFInfo
- Publication number
- CN105370272A CN105370272A CN201510898493.XA CN201510898493A CN105370272A CN 105370272 A CN105370272 A CN 105370272A CN 201510898493 A CN201510898493 A CN 201510898493A CN 105370272 A CN105370272 A CN 105370272A
- Authority
- CN
- China
- Prior art keywords
- resistivity
- logging
- density
- neutron
- well
- 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 abstract description 38
- 238000012937 correction Methods 0.000 claims abstract description 101
- 238000011156 evaluation Methods 0.000 claims abstract description 80
- 239000011435 rock Substances 0.000 claims description 32
- 239000012530 fluid Substances 0.000 claims description 26
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 25
- 239000011707 mineral Substances 0.000 claims description 25
- 230000000704 physical effect Effects 0.000 claims description 23
- 239000000470 constituent Substances 0.000 claims description 22
- 230000002159 abnormal effect Effects 0.000 claims description 14
- 239000004927 clay Substances 0.000 claims description 10
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 8
- 229910052799 carbon Inorganic materials 0.000 claims description 8
- 229910052683 pyrite Inorganic materials 0.000 claims description 7
- 239000011028 pyrite Substances 0.000 claims description 7
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 claims description 7
- 229910021532 Calcite Inorganic materials 0.000 claims description 3
- 238000002637 fluid replacement therapy Methods 0.000 claims description 3
- 239000010453 quartz Substances 0.000 claims description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 7
- 238000009825 accumulation Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 238000001730 gamma-ray spectroscopy Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000012360 testing 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
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
本发明提供一种页岩气储层的测井评价方法，包括：(A)获取页岩气储层的常规测井资料；(B)对获取的常规测井资料进行校正；(C)对校正后的常规测井资料进行一致性处理；(D)通过校正后的常规测井资料以及一致性处理后的常规测井资料进行测井评价，以获得页岩气储层的测井评价参数。在所述方法中，可仅通过对常规测井资料进行校正和一致性处理来对页岩气储层进行测井评价，从而获得准确的页岩气储层的测井评价参数。
Description
技术领域
本发明涉及油气勘探领域，更具体地讲，涉及一种页岩气储层的测井评价方法。
背景技术
页岩气作为一种重要的非常规天然气资源已经逐渐成为油气资源勘探开发的重点。页岩气储层在沉积背景、发育构造部位、分布与聚集规律、储层特征和渗流机理等方面与常规油气存在较大差异，因此，不能使用常规油气储层测井的评价方法来进行页岩气储层的测井评价。
页岩气储层的测井评价的主要目的是获得：矿物组分评价、地球化学参数评价、物性参数评价、含气性评价和可压裂性评价。目前，用于获得矿物组分评价的方法包括：常规测井方法、元素俘获能谱测井(ECS)方法、自然伽玛能谱测井(NGS)方法以及斯伦贝谢公司通过综合应用常规测井方法和ECS方法，来基于随机求解方法获得矿物组分评价的方法等。用于地球化学参数评价的方法主要包括：Schmoker利用密度和自然伽玛测井资料计算有机碳含量(TOC)的方法(Schmoker,1979,1981)、Passey等提出的利用声波曲线和地层电阻率曲线计算TOC的方法(Passeyetal.,1990)、Lewis等提出的利用ECS方法和干酪根转换因子等参数计算TOC的方法(Lewisetal.,2004)、Jacobi等提出的利用密度和核磁共振测井来计算TOC的方法(Jacobietal.,2009)、Pemper等提出的利用脉冲中子和自然伽玛能谱计算TOC的方法(Pemperetal.,2009)以及Khoshnoodkia等提出的利用电阻率、中子、密度和声波测井资料以及基于模糊逻辑技术建立的神经网络来计算TOC的方法(Khoshnoodkiaetal.,2011)。用于获得物性参数评价的方法包括：Herron等提出的利用ECS方法对孔隙度进行计算的方法(Herronetal.,2002)、Abu-Shanab等提出的利用密度与NMR孔隙度相结合的密度-核磁共振(DMR)孔隙度来评价致密储层孔隙度的方法(Abu-Shanabetal.,2005)、毛志强等提出的结合声波时差和核磁共振测井资料计算低孔渗气层真实孔隙度的方法(毛志强等，2010)以及LeCompte利用核磁共振测井(NMR)计算页岩孔隙度的方法(LeCompte,2010)。用于获得含气性评价的方法包括：现场解吸法、等温吸附法以及测井解释法。用于获得可压裂性评价的方法包括：岩石力学实验与阵列声波测井方法以及通过利用阵列声波测井提取地层纵、横波速度来基于密度测井资料计算岩石弹性力学参数和岩石强度参数的方法。
然而，上述页岩气储层的测井评价方法均需要核磁共振、电阻率成像和元素测井等非常规的成像测井资料，仅根据常规测井资料无法通过上述方法获得页岩气储层的准确测井评价。但是，由于我国页岩气储层勘探开发还处于起步阶段，针对页岩气勘探开发及评价的测井资料的积累不足，往往导致无法利用上述方法获得页岩气储层的测井评价。
因此，现有技术中的页岩气储层的测井评价方法无法通过常规测井资料对页岩气储层进行测井评价。
发明内容
本发明的示例性实施例在于提供一种页岩气储层的测井评价方法。所述方法能够克服现有技术中的页岩气储层的测井评价方法无法通过常规测井资料准确地对页岩气储层进行测井评价的缺陷。
根据本发明示例性实施例，提供一种页岩气储层的测井评价方法，包括：(A)获取页岩气储层的常规测井资料；(B)对获取的常规测井资料进行校正；(C)对校正后的常规测井资料进行一致性处理；(D)通过校正后的常规测井资料以及一致性处理后的常规测井资料进行测井评价，以获得页岩气储层的测井评价参数。
可选地，所述常规测井资料包括与岩石相关的测井资料以及与测井相关的测井资料，其中，与岩石相关的测井资料包括密度、声波时差、中子、自然伽马、自然电位、深电阻率、中电阻率和浅电阻率，与测井相关的测井资料包括井径。
可选地，步骤(B)包括：(B1)对密度、声波时差、中子、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径进行深度校正，获得深度校正后的密度、声波时差、中子、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径；(B2)对深度校正后的密度和中子进行进一步校正。
可选地，步骤(B2)包括：(B21)通过深度校正后的中子和密度获得中子-密度交会图，提取中子-密度交会图中预定区域之外的区域的作为异常值区域，并获得与所述异常值区域对应的异常测井段；(B22)获得整个测井段中与所述异常测井段的特性最接近的接近测井段，其中，所述特性包括：由深度校正后的自然伽马表示的岩性、由深度校正后的声波时差或中子表示的物性、以及由深度校正后的深电阻率、中电阻率或浅电阻率表示的流体性质；(B23)使用接近测井段的深度校正后的中子与接近测井段的岩性、物性和流体性质的拟合曲线替换所述异常测井段的深度校正后的中子与异常测井段的岩性、物性和流体性质的拟合曲线，来对所述异常测井段的深度校正后的中子进行进一步校正；(B24)使用接近测井段的深度校正后的密度与接近测井段的岩性、物性和流体性质的拟合曲线替换所述异常测井段的深度校正后的密度与异常测井段的岩性、物性和流体性质的拟合曲线，来对所述异常测井段的深度校正后的密度进行进一步校正。
可选地，步骤(C)包括：(C1)通过使用深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度获得分布稳定的页岩层段；(C2)获得深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度中的每个的去压实趋势；(C3)将深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度中的每个分别减去对应的去压实趋势来获得相应的差值，并通过使用每个差值获得深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度中的每个的直方图；(C4)通过所述分布稳定的页岩层段以及获得的相应直方图来对深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度进行一致性处理。
可选地，所述测井评价参数包括：矿物组分、地球化学参数、物性参数、含气性评价和可压裂性评价，其中，步骤(D)包括：(D1)通过使用一致性处理后的自然伽马以及进一步校正后的中子-密度交会图获得矿物组分中的粘土含量；(D2)通过使用深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度获得矿物组分中的黄铁矿含量和脆性矿物含量、地球化学参数中的有机碳体积组分含量、物性参数中的岩石孔隙度大小以及含气性评价中的含气饱和度大小；(D3)通过将获得的粘土含量、黄铁矿含量、脆性矿物含量、有机碳体积组分含量、岩石孔隙度大小输入到岩石物理建模中，来获得纵波速度、横波速度和密度，通过纵波速度、横波速度和进一步校正后的密度获得可压裂性评价的弹性参数，所述可压裂性评价的弹性参数为杨氏模量和泊松比。
可选地，所述脆性矿物含量包括方解石含量和石英含量，所述岩石物理建模包括骨架混合模型、流体模型、干岩石模型、流体替换技术方程。
在根据本发明示例性实施例的页岩气储层的测井评价方法中，可仅通过对常规测井资料进行校正和一致性处理来对页岩气储层进行测井评价，从而获得准确的页岩气储层的测井评价参数。
附图说明
通过下面结合示例性地示出实施例的附图进行的描述，本发明示例性实施例的上述和其他目的和特点将会变得更加清楚，其中：
图1示出根据本发明示例性实施例的页岩气储层的测井评价方法的流程图；
图2示出根据本发明示例性实施例的页岩气储层的测井评价方法中对常规测井资料进行校正步骤的流程图；
图3示出根据本发明示例性实施例的页岩气储层的测井评价方法中对密度和中子进行进一步校正步骤的流程图；
图4示出根据本发明示例性实施例的对密度和中子进行进一步校正后的中子和密度的曲线图；
图5示出根据本发明示例性实施例的页岩气储层的测井评价方法中对校正后的常规测井资料进行一致性处理步骤的流程图；
图6示出根据本发明示例性实施例的分布稳定的页岩层段的示图；
图7示出根据本发明示例性实施例的对进一步校正后的声波时差进行一致性处理之前和之后的声波时差的分布图；
图8示出根据本发明示例性实施例的对进一步校正后的密度进行一致性处理之前和之后的密度的分布图；
图9示出根据本发明示例性实施例的页岩气储层的测井评价方法中获得测井评价参数步骤的流程图；
图10示出根据本发明示例性实施例的通过迭代获得的孔隙度和饱和度的曲线图。
具体实施方式
以下，将参照附图更充分地描述本发明的示例性实施例，示例性实施例在附图中示出。然而，可以以许多不同的形式实施示例性实施例，并且不应被解释为局限于在此阐述的示例性实施例。相反，提供这些实施例从而本公开将会彻底和完整，并将完全地将示例性实施例的范围传达给本领域的技术人员。
图1示出根据本发明示例性实施例的页岩气储层的测井评价方法的流程图。
在步骤S100，获取页岩气储层的常规测井资料。这里，所述常规测井资料是指已经从现有的测井(包括几年之间甚至几十年之间的测井)中获得的测井资料。
作为示例，所述常规测井资料包括与岩石相关的测井资料以及与测井相关的测井资料，其中，与岩石相关的测井资料包括密度、声波时差、中子、自然伽马、自然电位、深电阻率、中电阻率和浅电阻率，与测井相关的测井资料包括井径。这里，与岩石相关的测井资料中的密度指测井中各个深度的岩石的密度、声波时差指测井中各个深度的岩石的声波、中子指测井中各个深度的岩石的中子信息、自然伽马指测井中各个深度的岩石的自然伽马信息、自然电位指测井中各个深度的岩石的自然电位信息、深电阻率指测井中各个深度的岩石的深电阻率信息、中电阻率指测井中各个深度的岩石的中电阻率信息、浅电阻率指测井中各个深度的岩石的浅电阻率信息。
在步骤S200，对获取的常规测井资料进行校正。这里，为了便于后续处理，以获得精确的测井评价参数，需要对获得的测井资料进行校正处理。
作为示例，可首先对常规测井资料进行深度校正，之后，对常规测井资料中的密度和中子进行进一步校正。
图2示出根据本发明示例性实施例的页岩气储层的测井评价方法中对常规测井资料进行校正步骤的流程图。
在步骤S210，对密度、声波时差、中子、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径进行深度校正，获得深度校正后的密度、声波时差、中子、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径。
这里，由于获得的常规测井资料中的每个均与深度相关，因此，为了有效地使用获得的常规测井资料，需要对常规测井资料进行深度匹配检查，即，深度校正。
在步骤S220，对深度校正后的密度和中子进行进一步校正。这里，为了校正获得的常规测井资料中的异常数据，需要对深度校正后的密度和中子进行进一步校正。
图3示出根据本发明示例性实施例的页岩气储层的测井评价方法中对密度和中子进行进一步校正步骤的流程图。
在步骤S221，通过深度校正后的中子和密度获得中子-密度交会图，提取中子-密度交会图中预定区域之外的区域的作为异常值区域，并获得与所述异常值区域对应的异常测井段。例如，可根据需要预先设定所述预定区域的范围，并将预定区域之外的区域作为异常值区域。
在步骤S222，获得整个测井段中与所述异常测井段的特性最接近的接近测井段，其中，所述特性包括：由深度校正后的自然伽马表示的岩性、由深度校正后的声波时差或中子表示的物性、以及由深度校正后的深电阻率、中电阻率或浅电阻率表示的流体性质。
作为示例，可将岩性与异常测井段的岩性之间的差异在第一预定范围内、物性与异常测井段的物性之间的差异在第二预定范围内且流体性质与异常测井段的流体性质之间的差异在第三预定范围内的测井段确定为接近测井段。这里，可根据实际需要来设定第一预定范围、第二预定范围和第三预定范围。
在步骤S223，使用接近测井段的深度校正后的中子与接近测井段的岩性、物性和流体性质的拟合曲线替换所述异常测井段的深度校正后的中子与异常测井段的岩性、物性和流体性质的拟合曲线，来对所述异常测井段的深度校正后的中子进行进一步校正。
作为示例，可将接近测井段的深度校正后的中子与表示岩性的深度校正后的自然伽马、表示物性的深度校正后的声波时差和表示流体性质的深度校正后的深电阻率的拟合曲线来替换异常测井段的深度校正后的中子与深度校正后的自然伽马、声波时差和深电阻率的拟合曲线。
应该理解，拟合虚线不限于上述拟合曲线，还可以是深度校正后的中子与表示岩性的深度校正后的自然伽马、表示物性的深度校正后的声波时差和表示流体性质的深度校正后的中电阻率的拟合曲线，或者是深度校正后的中子与表示岩性的深度校正后的自然伽马、表示物性的深度校正后的声波时差和表示流体性质的深度校正后的浅电阻率的拟合曲线。
在步骤S224使用接近测井段的深度校正后的密度与接近测井段的岩性、物性和流体性质的拟合曲线替换所述异常测井段的深度校正后的密度与异常测井段的岩性、物性和流体性质的拟合曲线，来对所述异常测井段的深度校正后的密度进行进一步校正。
作为示例，可将接近测井段的深度校正后的密度与表示岩性的深度校正后的自然伽马、表示物性的深度校正后的声波时差和表示流体性质的深度校正后的深电阻率的拟合曲线来替换异常测井段的深度校正后的密度与深度校正后的自然伽马、声波时差和深电阻率的拟合曲线。
应该理解，拟合虚线不限于上述拟合曲线，还可以是深度校正后的密度与表示岩性的深度校正后的自然伽马、表示物性的深度校正后的声波时差和表示流体性质的深度校正后的中电阻率的拟合曲线，或者是深度校正后的密度与表示岩性的深度校正后的自然伽马、表示物性的深度校正后的声波时差和表示流体性质的深度校正后的浅电阻率的拟合曲线。
图4示出根据本发明示例性实施例的对密度和中子进行进一步校正后的中子和密度的曲线图。
如图4所示，图4的左侧为对常规测井资料进行深度校正后的曲线，中间部分和右侧部分分别包括对中子和密度进行进一步校正后获得的曲线。从图中可以看出，对中子和密度进行进一步校正后获得的曲线已经去除了明显异常的部分。
再次参照图1，在步骤S300，对校正后的常规测井资料进行一致性处理。这里，为了便于后续处理，以获得精确的测井评价参数，需要对步骤S200获得的校正后的测井资料进行一致性处理。
图5示出根据本发明示例性实施例的页岩气储层的测井评价方法中对校正后的常规测井资料进行一致性处理步骤的流程图。
在步骤S310，通过使用深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度获得分布稳定的页岩层段。作为示例，可利用多井地层对比方法获得分布稳定的页岩层段。
图6示出根据本发明示例性实施例的分布稳定的页岩层段的示图。
如图6所示，图中虚线框中的部分的常规测井资料(井径差值、深电阻等参数)分布稳定，即，测井资料的曲线比较相似，因此可将虚线框中的部分确定为分布稳定的页岩层段。
再次参照图5，在步骤S320，获得深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度中的每个的去压实趋势。这里，可使用地球物理分析手段中的常规方法获得深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度中的每个的去压实趋势
在步骤S330，将深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度中的每个分别减去对应的去压实趋势来获得相应的差值，并通过使用每个差值获得深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度中的每个的直方图。这里，将去压实趋势应用于对测井资料进行一致性处理的过程中，可以很大提高一致性处理的精度。
在步骤S340，通过所述分布稳定的页岩层段以及获得的相应直方图来对深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度进行一致性处理。
下面，以常规测井资料中的声波时差和密度为例，来示出在步骤S300进行一致性处理的示例。
图7示出根据本发明示例性实施例的对进一步校正后的声波时差进行一致性处理之前和之后的声波时差的分布图。
如图7所示，从图7中的声波-深度的分布图可以看出，与图7中的左上角的一致性处理之前的声波-深度交会图的分布相比，图7中左下角的一致性处理之后的声波-深度交会图的分布变得更加集中。并且，从图7中的声波时差差值直方图可以看出，与图7右上角的一致性处理之前的声波时差差值直方图的分布相比，图7中左下角的一致性处理之后的声波时差差值直方图的分布变得更加集中。
图8示出根据本发明示例性实施例的对进一步校正后的密度进行一致性处理之前和之后的密度的分布图。
如图8所示，从图8中的密度-深度的分布图可以看出，与图8中的左上角的一致性处理之前的密度-深度交会图的分布相比，图8中左下角的一致性处理之后的密度-深度交会图的分布变得更加集中。并且，从图8中的密度差值直方图可以看出，与图8右上角的一致性处理之前的密度差值直方图的分布相比，图8中左下角的一致性处理之后的密度差值直方图的分布变得更加集中。
再次参照图1，在步骤S400，通过校正后的常规测井资料以及一致性处理后的常规测井资料进行测井评价，以获得页岩气储层的测井评价参数。作为示例，所述测井评价参数包括：矿物组分、地球化学参数、物性参数、含气性评价和可压裂性评价。
图9示出根据本发明示例性实施例的页岩气储层的测井评价方法中获得测井评价参数步骤的流程图。
在步骤S410，通过使用一致性处理后的自然伽马以及进一步校正后的中子-密度交会图获得矿物组分中的粘土含量。这里，作为示例，可首先通过使用进一步校正后的中子-密度交会图获得第一粘土含量，然后可通过使用一致性处理后的自然伽马获得与第一粘土含量最接近的第二粘土含量来作为矿物组分中的粘土含量。
在步骤S420，通过使用深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度获得矿物组分中的黄铁矿含量和脆性矿物含量、地球化学参数中的有机碳体积组分含量、物性参数中的岩石孔隙度大小以及含气性评价中的含气饱和度大小。
作为示例，可利用一致性处理后的常规测井资料获得一致性处理后的密度、声波时差、中子、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径中的每个所对应的曲线以及所有曲线对应的方程组。这里，所述方程组中包含测井评价参数所包括的矿物组分、地球化学参数、物性参数、含气性评价和可压裂性评价。为了获得准确的测井评价参数，可对所述方程组进行三轮迭代。具体地讲，由于测井评价参数中的物性参数中的岩石孔隙度以及含气性评价中的含气饱和度互相影响，因此，可首先假设地层中饱含水(即，令饱和度为1)，来求解岩石孔隙度的大小；其次，可将求解出的孔隙度代入方程组中，通过使用求解出的孔隙度大小来求解真实的饱和度大小；最后，可将求解出的孔隙度大小和饱和度大小代入到方程组中，来求解矿物组分中的黄铁矿含量和脆性矿物含量、地球化学参数中的有机碳体积组分含量。这里，作为示例，所述脆性矿物含量包括方解石含量和石英含量。
图10示出根据本发明示例性实施例的通过迭代获得的孔隙度和饱和度的曲线图。
如图10所示，通过三次迭代能够一步步获得更加精确的孔隙度大小和饱和度大小。
在步骤S430，通过将获得的粘土含量、黄铁矿含量、脆性矿物含量、有机碳体积组分含量、岩石孔隙度大小输入到岩石物理建模中，来获得纵波速度、横波速度和密度，通过纵波速度、横波速度和进一步校正后的密度获得可压裂性评价的弹性参数，所述可压裂性评价的弹性参数为杨氏模量和泊松比。作为示例，所述岩石物理建模包括骨架混合模型、流体模型、干岩石模型、流体替换技术方程。
在根据本发明示例性实施例的页岩气储层的测井评价方法中，可仅通过对常规测井资料进行校正和一致性处理来对页岩气储层进行测井评价，从而获得准确的页岩气储层的测井评价参数。
尽管已经参照其示例性实施例具体显示和描述了本发明，但是本领域的技术人员应该理解，在不脱离权利要求所限定的本发明的精神和范围的情况下，可以对其进行形式和细节上的各种改变。
Claims (7)
1.一种页岩气储层的测井评价方法，包括：
(A)获取页岩气储层的常规测井资料；
(B)对获取的常规测井资料进行校正；
(C)对校正后的常规测井资料进行一致性处理；
(D)通过校正后的常规测井资料以及一致性处理后的常规测井资料进行测井评价，以获得页岩气储层的测井评价参数。
2.如权利要求1所述的测井评价方法，其中，所述常规测井资料包括与岩石相关的测井资料以及与测井相关的测井资料，其中，与岩石相关的测井资料包括密度、声波时差、中子、自然伽马、自然电位、深电阻率、中电阻率和浅电阻率，与测井相关的测井资料包括井径。
3.如权利要求2所述的测井评价方法，其中，步骤(B)包括：
(B1)对密度、声波时差、中子、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径进行深度校正，获得深度校正后的密度、声波时差、中子、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径；
(B2)对深度校正后的密度和中子进行进一步校正。
4.如权利要求3所述的测井评价方法，其中，步骤(B2)包括：
(B21)通过深度校正后的中子和密度获得中子-密度交会图，提取中子-密度交会图中预定区域之外的区域的作为异常值区域，并获得与所述异常值区域对应的异常测井段；
(B22)获得整个测井段中与所述异常测井段的特性最接近的接近测井段，其中，所述特性包括：由深度校正后的自然伽马表示的岩性、由深度校正后的声波时差或中子表示的物性、以及由深度校正后的深电阻率、中电阻率或浅电阻率表示的流体性质；
(B23)使用接近测井段的深度校正后的中子与接近测井段的岩性、物性和流体性质的拟合曲线替换所述异常测井段的深度校正后的中子与异常测井段的岩性、物性和流体性质的拟合曲线，来对所述异常测井段的深度校正后的中子进行进一步校正；
(B24)使用接近测井段的深度校正后的密度与接近测井段的岩性、物性和流体性质的拟合曲线替换所述异常测井段的深度校正后的密度与异常测井段的岩性、物性和流体性质的拟合曲线，来对所述异常测井段的深度校正后的密度进行进一步校正。
5.如权利要求4所述的测井评价方法，其中，步骤(C)包括：
(C1)通过使用深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度获得分布稳定的页岩层段；
(C2)获得深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度中的每个的去压实趋势；
(C3)将深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度中的每个分别减去对应的去压实趋势来获得相应的差值，并通过使用每个差值获得深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度中的每个的直方图；
(C4)通过所述分布稳定的页岩层段以及获得的相应直方图来对深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度进行一致性处理。
6.如权利要求5所述的测井评价方法，其中，所述测井评价参数包括：矿物组分、地球化学参数、物性参数、含气性评价和可压裂性评价，
其中，步骤(D)包括：
(D1)通过使用一致性处理后的自然伽马以及进一步校正后的中子-密度交会图获得矿物组分中的粘土含量；
(D2)通过使用深度校正后的声波时差、自然伽马、自然电位、深电阻率、中电阻率、浅电阻率和井径以及进一步校正后的中子和密度获得矿物组分中的黄铁矿含量和脆性矿物含量、地球化学参数中的有机碳体积组分含量、物性参数中的岩石孔隙度大小以及含气性评价中的含气饱和度大小；
(D3)通过将获得的粘土含量、黄铁矿含量、脆性矿物含量、有机碳体积组分含量、岩石孔隙度大小输入到岩石物理建模中，来获得纵波速度、横波速度和密度，通过纵波速度、横波速度和进一步校正后的密度获得可压裂性评价的弹性参数，所述可压裂性评价的弹性参数为杨氏模量和泊松比。
7.如权利要求6所述的测井评价方法，其中，所述脆性矿物含量包括方解石含量和石英含量，所述岩石物理建模包括骨架混合模型、流体模型、干岩石模型、流体替换技术方程。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510898493.XA CN105370272A (zh) | 2015-12-07 | 2015-12-07 | 页岩气储层的测井评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510898493.XA CN105370272A (zh) | 2015-12-07 | 2015-12-07 | 页岩气储层的测井评价方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105370272A true CN105370272A (zh) | 2016-03-02 |
Family
ID=55372830
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510898493.XA Pending CN105370272A (zh) | 2015-12-07 | 2015-12-07 | 页岩气储层的测井评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105370272A (zh) |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105787803A (zh) * | 2016-03-17 | 2016-07-20 | 成都创源油气技术开发有限公司 | 页岩气藏三维地质建模方法 |
| CN105822294A (zh) * | 2016-03-17 | 2016-08-03 | 成都创源油气技术开发有限公司 | 页岩气藏地质评价单井建模方法 |
| CN106560820A (zh) * | 2016-10-09 | 2017-04-12 | 钦州学院 | 一种页岩气储层的测井评价方法 |
| CN106600436A (zh) * | 2016-11-28 | 2017-04-26 | 中国石油集团川庆钻探工程有限公司 | 计算页岩气地层矿物组分含量及孔隙度的方法 |
| CN108303510A (zh) * | 2017-12-25 | 2018-07-20 | 中国石油天然气股份有限公司 | 页岩气储层性能的评价方法、装置及计算机存储介质 |
| CN108301823A (zh) * | 2018-01-19 | 2018-07-20 | 北京捷贝通石油技术股份有限公司 | 一种识别储层油气甜点的方法 |
| CN108362621A (zh) * | 2018-02-06 | 2018-08-03 | 长江大学 | 基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法 |
| CN108647417A (zh) * | 2018-04-28 | 2018-10-12 | 中石化石油工程技术服务有限公司 | 一种确定页岩气储层含气饱和度的简易方法 |
| CN109632597A (zh) * | 2018-11-06 | 2019-04-16 | 广州海洋地质调查局 | 一种水合物与游离气存在状态的显示方法及处理终端 |
| CN109707378A (zh) * | 2019-02-20 | 2019-05-03 | 西北大学 | 一种基于泥浆侵入特性及纵向对比的低阻油层识别方法 |
| CN110703330A (zh) * | 2018-07-10 | 2020-01-17 | 中国石油化工股份有限公司 | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 |
| CN111042811A (zh) * | 2020-01-13 | 2020-04-21 | 中国石油天然气股份有限公司大港油田分公司 | 一种基于敏感参数叠合的页岩油产能评价方法 |
| CN111220522A (zh) * | 2019-10-28 | 2020-06-02 | 中国石油大学(华东) | 一种高泥质细粒沉积物水合物饱和度的岩心刻度测井计算方法 |
| CN113219553A (zh) * | 2021-04-25 | 2021-08-06 | 中国石油大学(北京) | 一种烃源岩总有机碳含量预测方法及装置 |
| CN113236237A (zh) * | 2021-05-11 | 2021-08-10 | 西南石油大学 | 一种基于常规测井曲线的致密砂岩储层有效性评价方法 |
| CN113804853A (zh) * | 2021-09-02 | 2021-12-17 | 中煤矿业集团有限公司 | 煤矿瓦斯含量检测方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103670388A (zh) * | 2013-12-12 | 2014-03-26 | 中国石油天然气股份有限公司 | 一种泥页岩有机碳含量的评价方法 |
| CN104047597A (zh) * | 2014-01-26 | 2014-09-17 | 中国石油大学(华东) | 含油气泥页岩地层测井曲线标准化方法 |
| CN104500049A (zh) * | 2014-10-20 | 2015-04-08 | 成都创源油气技术开发有限公司 | 页岩气地球物理快速评价方法 |
| CA2867580A1 (en) * | 2013-10-15 | 2015-04-15 | Korea Institute Of Geoscience And Mineral Resources (Kigam) | Fracking method at fracking intervals of a horizontal drilling zone in a sweet spot range determined based on measurements of resistivity and density in the horizontal drilling zone |
| CN104533400A (zh) * | 2014-11-12 | 2015-04-22 | 中国海洋石油总公司 | 一种重构测井曲线的方法 |
-
2015
- 2015-12-07 CN CN201510898493.XA patent/CN105370272A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2867580A1 (en) * | 2013-10-15 | 2015-04-15 | Korea Institute Of Geoscience And Mineral Resources (Kigam) | Fracking method at fracking intervals of a horizontal drilling zone in a sweet spot range determined based on measurements of resistivity and density in the horizontal drilling zone |
| CN103670388A (zh) * | 2013-12-12 | 2014-03-26 | 中国石油天然气股份有限公司 | 一种泥页岩有机碳含量的评价方法 |
| CN104047597A (zh) * | 2014-01-26 | 2014-09-17 | 中国石油大学(华东) | 含油气泥页岩地层测井曲线标准化方法 |
| CN104500049A (zh) * | 2014-10-20 | 2015-04-08 | 成都创源油气技术开发有限公司 | 页岩气地球物理快速评价方法 |
| CN104533400A (zh) * | 2014-11-12 | 2015-04-22 | 中国海洋石油总公司 | 一种重构测井曲线的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 周灿灿: "《复杂碎屑岩测井岩石物理与处理评价》", 31 October 2013 * |
| 国家能源局: "《页岩气测井资料处理与解释规范》", 15 October 2014 * |
Cited By (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105822294A (zh) * | 2016-03-17 | 2016-08-03 | 成都创源油气技术开发有限公司 | 页岩气藏地质评价单井建模方法 |
| CN105787803A (zh) * | 2016-03-17 | 2016-07-20 | 成都创源油气技术开发有限公司 | 页岩气藏三维地质建模方法 |
| CN106560820A (zh) * | 2016-10-09 | 2017-04-12 | 钦州学院 | 一种页岩气储层的测井评价方法 |
| CN106600436A (zh) * | 2016-11-28 | 2017-04-26 | 中国石油集团川庆钻探工程有限公司 | 计算页岩气地层矿物组分含量及孔隙度的方法 |
| CN106600436B (zh) * | 2016-11-28 | 2021-04-09 | 中国石油天然气集团有限公司 | 计算页岩气地层矿物组分含量及孔隙度的方法 |
| CN108303510A (zh) * | 2017-12-25 | 2018-07-20 | 中国石油天然气股份有限公司 | 页岩气储层性能的评价方法、装置及计算机存储介质 |
| CN108301823A (zh) * | 2018-01-19 | 2018-07-20 | 北京捷贝通石油技术股份有限公司 | 一种识别储层油气甜点的方法 |
| CN108362621A (zh) * | 2018-02-06 | 2018-08-03 | 长江大学 | 基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法 |
| CN108647417A (zh) * | 2018-04-28 | 2018-10-12 | 中石化石油工程技术服务有限公司 | 一种确定页岩气储层含气饱和度的简易方法 |
| CN110703330A (zh) * | 2018-07-10 | 2020-01-17 | 中国石油化工股份有限公司 | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 |
| CN110703330B (zh) * | 2018-07-10 | 2021-12-24 | 中国石油化工股份有限公司 | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 |
| CN109632597B (zh) * | 2018-11-06 | 2019-08-16 | 广州海洋地质调查局 | 一种水合物与游离气存在状态的显示方法及处理终端 |
| CN109632597A (zh) * | 2018-11-06 | 2019-04-16 | 广州海洋地质调查局 | 一种水合物与游离气存在状态的显示方法及处理终端 |
| CN109707378A (zh) * | 2019-02-20 | 2019-05-03 | 西北大学 | 一种基于泥浆侵入特性及纵向对比的低阻油层识别方法 |
| CN111220522A (zh) * | 2019-10-28 | 2020-06-02 | 中国石油大学(华东) | 一种高泥质细粒沉积物水合物饱和度的岩心刻度测井计算方法 |
| CN111042811A (zh) * | 2020-01-13 | 2020-04-21 | 中国石油天然气股份有限公司大港油田分公司 | 一种基于敏感参数叠合的页岩油产能评价方法 |
| CN113219553A (zh) * | 2021-04-25 | 2021-08-06 | 中国石油大学(北京) | 一种烃源岩总有机碳含量预测方法及装置 |
| CN113236237A (zh) * | 2021-05-11 | 2021-08-10 | 西南石油大学 | 一种基于常规测井曲线的致密砂岩储层有效性评价方法 |
| CN113236237B (zh) * | 2021-05-11 | 2022-05-03 | 西南石油大学 | 一种基于常规测井曲线的致密砂岩储层有效性评价方法 |
| CN113804853A (zh) * | 2021-09-02 | 2021-12-17 | 中煤矿业集团有限公司 | 煤矿瓦斯含量检测方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105370272A (zh) | 页岩气储层的测井评价方法 | |
| CN103454685B (zh) | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 | |
| Yu et al. | Porosity estimation in kerogen-bearing shale gas reservoirs | |
| CN104806232B (zh) | 一种确定碳酸盐岩储层孔隙度下限的方法 | |
| CN104516018B (zh) | 一种地球物理勘探中岩性约束下的孔隙度反演方法 | |
| CN103670388A (zh) | 一种泥页岩有机碳含量的评价方法 | |
| CN105651966A (zh) | 一种页岩油气优质储层评价方法及参数确定方法 | |
| CN106168685A (zh) | 一种页岩气单井地质综合评价方法 | |
| Gherabati et al. | Assessment of hydrocarbon in place and recovery factors in the Eagle Ford Shale play | |
| CN103135135A (zh) | 一种基于疏松砂岩模型进行烃类定量预测的方法和装置 | |
| CN105275456A (zh) | 一种利用测井资料识别优质泥页岩的方法 | |
| CN107133670A (zh) | 一种基于决策树数据挖掘算法的复杂岩性识别方法及系统 | |
| CN104573150A (zh) | 一种针对富含有机质页岩的各向异性岩石物理建模方法 | |
| CN103867198B (zh) | 碳酸盐岩天然气层地层密度判别法 | |
| CN103993871A (zh) | 针对薄互层地层的测井资料标准化处理方法及装置 | |
| CN105527653A (zh) | 一种基于地质信息的虚拟井构建方法 | |
| CN106223942A (zh) | 一种基于测井曲线重构的砾岩油藏泥质含量计算方法 | |
| CN103616724B (zh) | 非常规致密及泥页岩岩心岩性精细准确描述方法 | |
| Clerke et al. | Wireline spectral porosity analysis of the Arab Limestone—From Rosetta Stone to CIPHER | |
| CN104948176A (zh) | 一种基于渗透增大率识别碳酸盐岩储层裂缝的方法 | |
| CN107807410A (zh) | 储层地球物理敏感参数优选方法及系统 | |
| CN106556863A (zh) | 基于深度域叠前角道集的孔隙度预测方法 | |
| CN107688200A (zh) | 针对低幅度构造畸变的速度校正方法 | |
| CN105093304A (zh) | 一种地球物理勘探利用测井曲线自动求取岩性曲线的方法 | |
| CN110118994B (zh) | 一种基于地震反演和机器学习的陆相烃源岩定量预测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180507 Address after: 071000 No. 189, Fan Yang Xi Road, Zhuozhou, Baoding, Hebei Applicant after: BGP INC., CHINA NATIONAL PETROLEUM Corp. Address before: 610213 216, Huayang Road, Huayang Town, Huayang Town, Shuangliu County, Chengdu, Sichuan, Sichuan Science and Technology Department of geophysical exploration company Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| TA01 | Transfer of patent application right | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160302 |
|
| RJ01 | Rejection of invention patent application after publication |