CN110344822B - 含钙质致密储层含水饱和度的确定方法 - Google Patents
含钙质致密储层含水饱和度的确定方法 Download PDFInfo
- Publication number
- CN110344822B CN110344822B CN201810287787.2A CN201810287787A CN110344822B CN 110344822 B CN110344822 B CN 110344822B CN 201810287787 A CN201810287787 A CN 201810287787A CN 110344822 B CN110344822 B CN 110344822B
- Authority
- CN
- China
- Prior art keywords
- data
- calcareous
- resistivity
- saturation
- experimental
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 108
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000004364 calculation method Methods 0.000 claims abstract description 96
- 239000011435 rock Substances 0.000 claims abstract description 88
- 238000004458 analytical method Methods 0.000 claims abstract description 32
- 238000012937 correction Methods 0.000 claims abstract description 27
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 claims description 75
- 229910052791 calcium Inorganic materials 0.000 claims description 75
- 239000011575 calcium Substances 0.000 claims description 75
- 239000008398 formation water Substances 0.000 claims description 58
- 238000010606 normalization Methods 0.000 claims description 21
- 239000011148 porous material Substances 0.000 claims description 20
- 238000001514 detection method Methods 0.000 claims description 19
- 238000010839 reverse transcription Methods 0.000 claims description 18
- 239000012530 fluid Substances 0.000 claims description 17
- 230000005611 electricity Effects 0.000 claims description 14
- 239000000523 sample Substances 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 229920006395 saturated elastomer Polymers 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 7
- 239000000463 material Substances 0.000 abstract description 5
- 238000013461 design Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000018109 developmental process Effects 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000007596 consolidation process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 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
-
- 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
技术领域
本发明涉及油气勘探开发领域,特别涉及一种含钙质致密储层含水饱和度的确定方法。
背景技术
油气储层的含水饱和度是油气勘探开发的重要参数,通过确定含水饱和度,有利于对油气的勘探开发。然而,含钙质致密储层包括大量的钙质,其具有电阻率高、补偿密度值高、自然伽马值低、声波时差值低等特点。并且,由于钙质的存在,使得含钙质致密储层的孔隙空间和孔喉半径减小,不易确定其含水饱和度,不利于对含钙质致密储层的勘探开发,因此,有必要提供一种含钙质致密储层含水饱和度的确定方法。
相关技术中,通过以下方法对含钙质致密储层的含水饱和度进行确定:选取岩心试样,通过试验获取岩心试样的孔隙度、地层水电阻率(RW)、用地层水100%饱和后的电阻率值R0。根据上述数据建立胶结指数m和饱和度指数n为定值的阿尔奇含水饱和度计算模型,然后利用阿尔奇含水饱和度计算模型对含钙质致密储层含水饱和度进行确定。
发明人发现相关技术中至少存在以下问题:
相关技术中,并未考虑到钙质体积百分比对致密储层含水饱和度的影响,通过相关技术提供的方法不能精确地确定含钙质致密储层的含水饱和度。
发明内容
本发明实施例提供了一种含钙质致密储层含水饱和度的确定方法,可解决上述技术问题。具体技术方案如下:
本发明实施例提供了一种含钙质致密储层含水饱和度的确定方法,所述方法包括:
获取目标含钙质致密储层的基础声波时差数据、基础深探测电阻率数据、基础钙质体积百分比数据、基础孔隙度数据、岩电分析数据、含钙质岩石物理体积模型;
根据所述基础声波时差数据、所述基础深探测电阻率数据、所述基础钙质体积百分比数据、所述基础孔隙度数据,确定钙质含量计算模型;
根据所述含钙质岩石物理体积模型,确定含钙质岩石物理体积计算模型;
根据所述钙质含量计算模型、所述含钙质岩石物理体积计算模型,确定孔隙度计算模型;
根据所述基础深探测电阻率数据、所述基础钙质体积百分比数据,确定电阻率钙质校正模型;
根据所述岩电分析数据,建立胶结指数、饱和度指数为变动常数的阿尔奇含水饱和度计算模型;
根据实际深探测电阻率数据、实际声波时差数据、所述钙质含量计算模型、所述含钙质岩石物理体积计算模型、所述阿尔奇含水饱和度计算模型、所述孔隙度计算模型、所述电阻率钙质校正模型,确定所述目标含钙质致密储层的含水饱和度。
在一种可能的设计中,所述根据所述基础声波时差数据、所述基础深探测电阻率数据、所述基础钙质体积百分比数据、所述基础孔隙度数据,确定钙质含量计算模型,包括:
将所述基础深探测电阻率数据的对数值和所述基础声波时差数据进行数据归一化处理;
将数据归一化后的所述基础深探测电阻率数据的对数值与数据归一化后的所述基础声波时差数据之比、所述基础钙质体积百分比数据与所述基础孔隙度数据之比进行交会图分析,确定钙质含量计算模型为:
在一种可能的设计中,所述根据所述含钙质岩石物理体积模型,确定含钙质岩石物理体积计算模型,包括:
获取孔隙中流体的声波时差数值、钙质的声波时差数值、岩石骨架声波时差数值;
根据所述含钙质岩石物理体积模型、所述孔隙中流体的声波时差数值、所述钙质的声波时差数值、所述岩石骨架声波时差数值,确定所述含钙质岩石物理体积计算模型为:
其中,φ为目标含钙质致密储层的孔隙度;Vca为钙质体积百分比;Vm为岩石骨架体积百分比;AC为声波时差;DTf为孔隙中流体的声波时差;DTca为钙质的声波时差;DTm为岩石骨架声波时差。
在一种可能的设计中,所述孔隙度计算模型为:
其中,φ为孔隙度;DTf为孔隙中流体的声波时差;DTca为钙质的声波时差;DTm为岩石骨架声波时差;A和B分别为第一系数和第二系数;为数据归一化后的深探测电阻率的对数值与数据归一化后的声波时差之比;AC为声波时差。
在一种可能的设计中,所述根据所述基础深探测电阻率数据、所述基础钙质体积百分比数据,确定电阻率钙质校正模型,包括:
将所述基础深探测电阻率数据的对数值的增大倍数与所述基础钙质体积百分比数据进行交会图分析,确定所述电阻率钙质校正模型;
所述电阻率钙质校正模型为:
在一种可能的设计中,所述根据所述岩电分析数据,建立胶结指数、饱和度指数为变动常数的阿尔奇含水饱和度计算模型,包括:
所述岩电分析数据包括:多个岩心试样的实验孔隙度数据、实验地层水电阻率数据、用地层水100%饱和后的电阻率数据、实验深探测电阻率数据、实验含水饱和度数据;
根据所述实验孔隙度数据、所述实验地层水电阻率数据、所述用地层水100%饱和后的电阻率数据,确定实验胶结指数数据和第一岩性系数;
将所述实验孔隙度数据与所述实验胶结指数数据进行交会图分析,确定变动常数胶结指数为:
其中,m为胶结指数;E、G、H分别为第五系数、第六系数、第七系数;φ为孔隙度;φ1为孔隙度分段分界值;
根据所述实验含水饱和度数据、所述实验深探测电阻率数据、所述用地层水100%饱和后的电阻率数据,确定实验饱和度指数数据和第二岩性系数;
将所述实验孔隙度数据与所述实验饱和度指数数据进行交会图分析,确定变动常数饱和度指数:
n=I×φ+J (7)
其中,n为饱和度指数;I、J分别为第八系数和第九系数;
将变动常数胶结指数m和饱和度指数n、以及第一岩性系数a和第二岩性系数b代入公式(8),得到所述阿尔奇含水饱和度计算模型:
其中,Sw为含水饱和度;a、b分别为第一岩性系数和第二岩性系数;RW、RTJ分别为地层水电阻率、校正后的深探测电阻率;φ代表孔隙度。
在一种可能的设计中,所述根据所述实验孔隙度数据、所述实验地层水电阻率数据、所述用地层水100%饱和后的电阻率数据,确定实验胶结指数数据和第一岩性系数,包括:
将所述用地层100%饱和后的电阻率数据与所述实验地层水电阻率数据之比的对数值,与所述实验孔隙度数据的对数值进行交会,确定所述实验胶结指数数据和第一岩性系数:
其中,R0为用地层水100%饱和后的电阻率;RW实验地层水电阻率;φ为孔隙度;a为第一岩心系数;m为胶结指数。
在一种可能的设计中,所述根据所述实验含水饱和度数据、所述实验深探测电阻率数据、所述用地层水100%饱和后的电阻率数据,确定实验饱和度指数数据和第二岩性系数,包括:
将所述实验深探测电阻率数据与所述用地层水100%饱和后的电阻率数据之比的对数值,与所述实验含水饱和度数据的对数值进行交会,确定所述实验饱和度指数数据和第二岩性系数:
其中,RT为实验深探测电阻率与用地层水100%饱和后的电阻率数据;SW为含水饱和度;b为第二岩性系数;n为饱和指数。
本发明实施例提供的技术方案带来的有益效果是:
本发明实施例提供的含钙质致密储层含水饱和度的确定方法,通过获取目标含钙质致密储层的基础声波时差数据、基础深探测电阻率数据、基础钙质体积百分比数据、基础孔隙度数据、岩电分析数据、含钙质岩石物理体积模型,并根据这些数据和模型确定钙质含量计算模型、含钙质岩石物理体积计算模型、孔隙度计算模型、电阻率钙质校正模型,建立胶结指数、饱和度指数为变动常数的阿尔奇含水饱和度计算模型。最后,根据实际深探测电阻率数据、实际声波时差数据、阿尔奇含水饱和度计算模型、孔隙度计算模型、电阻率钙质校正模型,对含钙质致密储层含水饱和度进行确定。通过上述方法,不仅考虑了钙质对含钙质致密储层含水饱和度的影响,还校正了深探测电阻率,减少了钙质对深探测电阻率的影响。进一步地,通过建立胶结指数、饱和度指数为变动常数的阿尔奇含水饱和度计算模型,还考虑到岩电分析数据对胶结指数和饱和度指数的影响,使得采用该方法可精确地确定含钙质致密储层的含水饱和度,利于对含钙质致密储层的勘探开发。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的含钙质致密储层含水饱和度的确定方法的流程图;
图2是实施例1提供的目标含钙质致密储层的测井曲线图;
图3是实施例1提供的数据归一化后的深探测电阻率数据的对数值与声波时差数据之比,与钙质体积百分比数据与孔隙度数据之比的交会图;
图4是实施例1提供的含钙质岩石物理体积模型图;
图5是实施例1提供的深探测电阻率数据的对数值的增大倍数与钙质体积百分比的交会图;
图6是实施例1提供的胶结指数与孔隙度数据的交会图;
图7是实施例1提供的饱和度指数与孔隙度数据的交会图;
图8是实施例1提供的目标含钙质致密储层饱和度确定后的测井处理成果图。
具体实施方式
除非另有定义,本发明实施例所用的所有技术术语均具有与本领域技术人员通常理解的相同的含义。在对本发明实施方式作进一步地详细描述之前,对理解本发明实施例一些术语给出定义。
本发明中所涉及的“含水饱和度”指的是:水占空隙的体积与岩石孔隙体积之比。
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本发明实施例提供了一种含钙质致密储层含水饱和度的确定方法,如附图1所示,该方法包括:
步骤101、获取目标含钙质致密储层的基础声波时差数据、基础深探测电阻率数据、基础钙质体积百分比数据、基础孔隙度数据、岩电分析数据、含钙质岩石物理体积模型。
步骤102、根据基础声波时差数据、基础深探测电阻率数据、基础钙质体积百分比数据、基础孔隙度数据,确定钙质含量计算模型。
步骤103、根据含钙质岩石物理体积模型,确定含钙质岩石物理体积计算模型。
步骤104、根据钙质含量计算模型、含钙质岩石物理体积计算模型,确定孔隙度计算模型。
步骤105、根据基础深探测电阻率数据、基础钙质体积百分比数据,确定电阻率钙质校正模型。
步骤106、根据岩电分析数据,建立胶结指数、饱和度指数为变动常数的阿尔奇含水饱和度计算模型。
步骤107、根据实际深探测电阻率数据、实际声波时差数据、钙质含量计算模型、含钙质岩石物理体积计算模型、阿尔奇含水饱和度计算模型、孔隙度计算模型、电阻率钙质校正模型,确定目标含钙质致密储层的含水饱和度。
本发明实施例提供的含钙质致密储层含水饱和度的确定方法,通过获取目标含钙质致密储层的基础声波时差数据、基础深探测电阻率数据、基础钙质体积百分比数据、基础孔隙度数据、岩电分析数据、含钙质岩石物理体积模型,并根据这些数据和模型确定钙质含量计算模型、含钙质岩石物理体积计算模型、孔隙度计算模型、电阻率钙质校正模型,建立胶结指数、饱和度指数为变动常数的阿尔奇含水饱和度计算模型。最后,根据实际深探测电阻率数据、实际声波时差数据、阿尔奇含水饱和度计算模型、孔隙度计算模型、电阻率钙质校正模型,确定目标含钙质致密储层的含水饱和度。通过上述方法,不仅考虑了钙质对含钙质致密储层含水饱和度的影响,还校正了深探测电阻率,减少了钙质对深探测电阻率的影响。进一步地,通过建立胶结指数、饱和度指数为变动常数的阿尔奇含水饱和度计算模型,还考虑到岩电分析数据对胶结指数和饱和度指数的影响,使得采用该方法可精确地确定含钙质致密储层的含水饱和度,利于对含钙质致密储层的勘探开发。
以下对本发明实施例提供的含钙质致密储层含水饱和度的确定方法进行详细阐述:
步骤101中所涉及的“基础声波时差数据”和“基础深探测电阻率数据”均由已经得到的历史测井曲线中提取得到。
基础声波时差数据指的是:通过测井仪器探测声波在地层中通过的时间,可以从含钙质部分的测井曲线中提取得到。
基础深探测电阻率数据指的是:通过测井仪器探测地层的电阻率值,可以从含钙质部分的测井曲线中提取得到。
“基础钙质体积百分比数据”和“基础孔隙度数据”均包括两部分:一部分为通过岩心试验获取,另一部分通过测井仪器根据基础声波时差数据和基础深探测电阻率数据计算而间接得到。
步骤102、根据基础声波时差数据、基础深探测电阻率数据、基础钙质体积百分比数据、基础孔隙度数据,确定钙质含量计算模型,包括:
步骤1021、将基础深探测电阻率数据的对数值和基础声波时差数据进行数据归一化处理。
可以理解的是,深探测电阻率数据、声波时差数据、钙质体积百分比数据、孔隙度数据是一一对应的。
由于基础深探测电阻率数据的单位与基础声波时差数据的单位不同,通过数据归一化处理,可将基础深探测电阻率数据、基础声波时差数据的单位统一,为后期建立钙质含量计算模型奠定了基础。
其中,数据归一化处理的方法为本领域所熟知的方法,其包括:把基础深探测电阻率数据的对数值和基础声波时差数据变成0~1之间的小数。
步骤1022、将数据归一化后的基础深探测电阻率数据的对数值与数据归一化后的基础声波时差数据之比、基础钙质体积百分比数据与基础孔隙度数据之比进行交会图分析,确定钙质含量计算模型为:
具体地,以为横坐标,为纵坐标,建立坐标图,根据基础深探测电阻率数据的对数值数据归一化后与数据归一化后的基础声波时差数据之比、基础钙质体积百分比数据与基础孔隙度数据之比,在该坐标图中做出散点,并拟合,确定第一系数A和第二系数B的数值,进而确定钙质含量计算模型。
步骤103中,根据含钙质岩石物理体积模型,确定含钙质岩石物理体积计算模型,包括:
获取孔隙中流体的声波时差数值、钙质的声波时差数值、岩石骨架声波时差数值;
根据含钙质岩石物理体积模型(如附图2所示)、孔隙中流体的声波时差数值、钙质的声波时差数值、岩石骨架声波时差数值,确定含钙质岩石物理体积计算模型为:
其中,φ为目标含钙质致密储层的孔隙度;Vca为钙质体积百分比;Vm为岩石骨架体积百分比;AC为声波时差;DTf为孔隙中流体的声波时差;DTca为钙质的声波时差;DTm为岩石骨架声波时差。
需要说明的是,孔隙中流体的声波时差DTf、钙质的声波时差DTca、岩石骨架声波时差DTm,均为定值,可通过经验获取其具体数值。
步骤104中,根据钙质含量计算模型、含钙质岩石物理体积计算模型,确定孔隙度计算模型。具体地,可将公式(1)、公式(2)结合,推导出孔隙度计算模型为:
其中,φ为孔隙度;DTf为孔隙中流体的声波时差;DTca为钙质的声波时差;DTm为岩石骨架声波时差;A和B分别为第一系数和第二系数;为深探测电阻率数据的对数值数据归一化后与数据归一化后的声波时差之比;AC为声波时差。
步骤105中,根据基础深探测电阻率数据、基础钙质体积百分比数据,确定电阻率钙质校正模型,包括:
将基础深探测电阻率数据的对数值的增大倍数与基础钙质体积百分比数据进行交会图分析,确定电阻率钙质校正模型;
电阻率钙质校正模型为:
具体地,以钙质体积百分比为横坐标,以深探测电阻率的对数值的增大倍数为纵坐标,建立散点图,然后将散点图拟合,得到第三系数C和第四系数D的具体数值,进而确定电阻率钙质校正模型。
通过确定电阻率钙质校正模型,可减小钙质对深探测电阻率的影响,为后期高精度地确定阿尔奇含水饱和度计算模型奠定了基础。
步骤106中,根据岩电分析数据,建立胶结指数、饱和度指数为变动常数的阿尔奇含水饱和度计算模型,包括:
岩电分析数据包括:多个岩心试样的实验孔隙度数据、实验地层水电阻率数据、用地层水100%饱和后的电阻率数据、实验深探测电阻率数据、实验含水饱和度数据。
岩电分析数据根据不同岩心试样分析的钙质体积百分比、孔隙度、100%饱和水电阻率值,以及同一个样品的不同含水比例电阻率值及对应的含水饱和度、水的电阻率值。
步骤1061、根据实验孔隙度数据、实验地层水电阻率数据、用地层水100%饱和后的电阻率数据,确定实验胶结指数数据和第一岩性系数;
其中,根据实验孔隙度数据、实验地层水电阻率数据、用地层水100%饱和后的电阻率数据,确定实验胶结指数数据和第一岩性系数,包括:
将用地层100%饱和后的电阻率数据与实验地层水电阻率数据之比的对数值,与实验孔隙度数据的对数值进行交会,确定实验胶结指数数据和第一岩性系数:
其中,R0为用地层水100%饱和后的电阻率数据;RW实验地层水电阻率数据;φ为孔隙度;a为第一岩心系数;m为胶结指数。
步骤1062、将实验孔隙度数据与实验胶结指数进行交会图分析,确定变动常数胶结指数为:
其中,m为胶结指数;E、G、H分别为第五系数、第六系数、第七系数;φ为孔隙度;φ1为孔隙度分段分界值。
具体地,以孔隙度φ为横坐标,以胶结指数m为纵坐标,建立坐标图。根据实验孔隙度数据、实验胶结指数数据在该坐标图中做出散点,根据散点拟合得到公式(6)。根据公式(6)得到第五系数E、第六系数G、第七系数H的具体数值。
步骤1063、根据实验含水饱和度数据、实验深探测电阻率数据、用地层水100%饱和后的电阻率数据,确定实验饱和度指数数据和第二岩性系数;
根据实验含水饱和度数据、实验深探测电阻率数据、用地层水100%饱和后的电阻率数据,确定实验饱和度指数数据和第二岩性系数,包括:
将实验深探测电阻率数据与用地层水100%饱和后的电阻率数据之比的对数值,与实验含水饱和度数据的对数值进行交会,确定实验饱和度指数数据和第二岩性系数:
其中,RT为实验深探测电阻率数据与用地层水、R0为用地层水100%饱和后的电阻率数据;SW为实验含水饱和度数据;b为第二岩性系数;n为饱和指数。
步骤1064、将实验孔隙度数据与实验饱和度指数进行交会图分析,确定变动常数饱和度指数:
n=I×φ+J (7)
其中,n为饱和度指数;I、J分别为第八系数和第九系数。
具体地,以孔隙度φ为横坐标,以饱和度指数n为纵坐标,建立坐标图。根据实验孔隙度数据、饱和度指数数据在该坐标图中做出散点,根据散点拟合得到公式(7)。根据公式(7)得到第八系数I、第九系数J的具体数值。
将变动常数胶结指数m和饱和度指数n、以及第一岩性系数a和第二岩性系数b代入公式(8),得到阿尔奇含水饱和度计算模型为:
其中,Sw为含水饱和度;a、b分别为第一岩性系数和第二岩性系数;RW、RTJ分别为地层水电阻率、校正后的深探测电阻率;φ代表孔隙度。
由于不同的孔隙度对应不同的胶结指数和饱和度指数,通过将胶结指数和饱和度与孔隙度建立关系,并应用于阿尔奇含水饱和度计算模型中,使该阿尔奇含水饱和度计算模型能够精确地确定含钙致密储层含水饱和度。
步骤107中,根据实际深探测电阻率数据、实际声波时差数据、钙质含量计算模型、含钙质岩石物理体积计算模型、阿尔奇含水饱和度计算模型、孔隙度计算模型、电阻率钙质校正模型,确定目标含钙质致密储层的含水饱和度。
具体地,通过井下测井方法获取实际深探测电阻率数据、实际声波时差数据,并将相应的数据代入钙质含量计算模型、孔隙度计算模型以及含钙质岩石物理体积计算模型中,计算得到模拟孔隙度和模拟钙质体积百分比。根据模拟孔隙度,确定胶结指数m和饱和度指数n。
将实际声波时差数据、模拟钙质体积百分比代入电阻率钙质校正模型中,得到校正后的深探测电阻率。
最后,将胶结指数m、饱和度指数n、校正后的深探测电阻率代入阿尔奇含水饱和度计算模型中,得到含水饱和度。
需要说明的是,实际探测电阻率数据和实际声波时差数据均通过测井仪器探测得到的测井曲线中获取。并且,该实际探测电阻率数据和实际声波时差数据为含钙质部分、非含钙质部分等所有的探测电阻率数据和声波时差数据。
以下将通过具体实施例进一步地描述本发明。
实施例1
本实施例以四川盆地川中凉上段地区某一含钙致密储层为目标含钙致密储层,通过本发明提供的方法对该目标含钙致密储层的含水饱和度进行确定,具体包括以下步骤:
步骤201、根据附图3提供的测井曲线图,提取含钙质储层的RT、AC曲线所对应的数据,作为基础深探测电阻率数据和基础声波时差数据。并且,获取一部分基础钙质体积百分比数据和基础孔隙度数据。通过岩心试样实验获取另一部分基础钙质体积百分比数据和基础孔隙度数据。通过对多个岩心试样进行岩电分析,获取多个岩心试样的实验孔隙度数据、实验地层水电阻率数据、用地层水100%饱和后的电阻率数据、实验深探测电阻率数据、实验含水饱和度数据。获取含钙质岩石物理体积模型,如附图2所示。
步骤202、将基础深探测电阻率数据的对数值与基础声波时差数据进行数据归一化处理。
将深探测电阻率数据的对数值数据归一化后与数据归一化后的声波时差数据之比为横坐标、钙质体积百分比数据与孔隙度数据之比为纵坐标,建立坐标图,如附图4所示,根据基础深探测电阻率数据的对数值数据归一化后与基础声波时差数据之比、基础钙质体积百分比数据与基础孔隙度数据之比,在该坐标图中做出散点,并拟合,确定第一系数A和第二系数B的数值,进而确定钙质含量计算模型为:
步骤203、根据含钙质岩石物理体积模型,参见附图2所示,确定含钙质岩石物理体积计算模型。
获取孔隙中流体的声波时差数值(取189us/ft)、钙质的声波时差数值(取47.5us/ft)、岩石骨架声波时差数值(可取51~58μs/ft,在此取55μs/ft)。
根据含钙质岩石物理体积模型、孔隙中流体的声波时差数值、钙质的声波时差数值、岩石骨架声波时差数值,确定含钙质岩石物理体积计算模型为:
其中,φ为目标含钙质致密储层的孔隙度;Vca为钙质体积百分比;Vm为岩石骨架体积百分比;AC为声波时差值;DTf为孔隙中流体的声波时差值,取189us/ft;DTca为钙质的声波时差值,取47.5us/ft;DTm为岩石骨架声波时差值,可以取51~58μs/ft,在此取55μs/ft。
步骤204、将公式(1)、公式(2)结合,推导出孔隙度计算模型为:
其中,φ为孔隙度;DTf为孔隙中流体的声波时差;DTca为钙质的声波时差;DTm为岩石骨架声波时差;A和B分别为第一系数和第二系数;为数据归一化后的深探测电阻率数据的对数值与数据归一化后的声波时差数据之比;AC为声波时差。
步骤205、以基础钙质体积百分比数据为横坐标,以基础深探测电阻率数据的对数值的增大倍数为纵坐标,建立散点图,如附图5所示,然后将散点拟合,以确定第三系数C、第四系数D的具体数值,进而确定电阻率钙质校正模型。
电阻率钙质校正模型为:
步骤206、获取岩电分析数据,其包括:多个岩心试样的实验孔隙度数据、实验地层水电阻率数据、用地层水100%饱和后的电阻率数据、实验深探测电阻率数据、实验含水饱和度数据。
根据实验孔隙度数据、实验地层水电阻率数据、用地层水100%饱和后的电阻率数据,确定实验胶结指数数据和第一岩性系数;
其中,根据实验孔隙度数据、实验地层水电阻率数据、用地层水100%饱和后的电阻率数据,确定实验胶结指数数据和第一岩性系数,包括:
将用地层100%饱和后的电阻率数据与实验地层水电阻率数据之比的对数值,与实验孔隙度数据的对数值进行交会,确定实验胶结指数数据和第一岩性系数:
其中,R0为用地层水100%饱和后的电阻率数据;RW实验地层水电阻率数据;φ为孔隙度;a为第一岩心系数,这里取值为1;m为胶结指数。
以孔隙度φ为横坐标,以胶结指数m为纵坐标,建立坐标图,附图6所示。根据实验孔隙度数据、实验胶结指数数据在该坐标图中做出散点,并将散点拟合得到第五系数E、第六系数G、第七系数H的具体数值,进而确定公式(6)。
根据实验含水饱和度数据、实验深探测电阻率数据、用地层水100%饱和后的电阻率数据,确定实验饱和度指数数据和第二岩性系数,包括:
将实验深探测电阻率数据与用地层水100%饱和后的电阻率数据之比的对数值,与实验含水饱和度数据的对数值进行交会,确定实验饱和度指数数据和第二岩性系数:
其中,RT为实验深探测电阻率数据与用地层水、R0为用地层水100%饱和后的电阻率数据;SW为实验含水饱和度数据;b为第二岩性系数,这里取值为1;n为饱和指数。
以孔隙度φ为横坐标,以饱和度指数n为纵坐标,建立坐标图,如附图7所示。根据实验孔隙度数据、实验饱和度指数数据在该坐标图中做出散点,根据散点拟合得到第八系数I、第九系数J的具体数值,进而确定公式(7)。
n=-0.26×φ+4.56 (7)
将变动常数胶结指数m和饱和度指数n、以及第一岩性系数a和第二岩性系数b代入公式(8),得到阿尔奇含水饱和度计算模型为:
其中,Sw为含水饱和度;a、b分别为第一岩性系数和第二岩性系数;RW、RTJ分别为地层水电阻率、校正后的深探测电阻率;φ代表孔隙度。
步骤207、测井曲线实际深探测电阻率数据、实际声波时差数据,并将相应的数据代入钙质含量计算模型、孔隙度计算模型以及含钙质岩石物理体积计算模型中,计算得到模拟孔隙度和模拟钙质体积百分比。根据模拟孔隙度,确定对应的胶结指数m和饱和度指数n。
将实际声波时差数据、模拟钙质体积百分比代入电阻率钙质校正模型中,得到校正后的深探测电阻率。
最后,将胶结指数m、饱和度指数n、校正后的深探测电阻率代入阿尔奇含水饱和度计算模型中,得到含水饱和度。
通过上述方法得到测井处理成果图如附图8所示,由附图8可以确定:目标含钙致密储层的2546~2550m、2556~2559m层段为干层。
通过上述方法获取的含水饱和度,与通过岩心实验获取的含水饱和度的吻合度较好,而采用相关技术中未校正的胶结指数和饱和指数得到的阿尔奇模型计算得到的含水饱和度较低。说明本发明实施例提供的含钙致密储层含水饱和度的确定方法比较精确,能够应用于对含钙致密储层含水饱和对的确定。
以上所述仅为本发明的较佳实施例,并不用以限制本发明的保护范围,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种含钙质致密储层含水饱和度的确定方法,其特征在于,所述方法包括:
获取目标含钙质致密储层的基础声波时差数据、基础深探测电阻率数据、基础钙质体积百分比数据、基础孔隙度数据、岩电分析数据、含钙质岩石物理体积模型;
将所述基础深探测电阻率数据的对数值和所述基础声波时差数据进行数据归一化处理;
将数据归一化后的所述基础深探测电阻率数据的对数值与数据归一化后的所述基础声波时差数据之比、所述基础钙质体积百分比数据与所述基础孔隙度数据之比进行交会图分析,确定钙质含量计算模型为:
根据所述含钙质岩石物理体积模型,确定含钙质岩石物理体积计算模型;
根据所述钙质含量计算模型、所述含钙质岩石物理体积计算模型,确定孔隙度计算模型;
根据所述基础深探测电阻率数据、所述基础钙质体积百分比数据,确定电阻率钙质校正模型;
根据所述岩电分析数据,建立胶结指数、饱和度指数为变动常数的阿尔奇含水饱和度计算模型;
根据实际深探测电阻率数据、实际声波时差数据、所述钙质含量计算模型、所述含钙质岩石物理体积计算模型、所述阿尔奇含水饱和度计算模型、所述孔隙度计算模型、所述电阻率钙质校正模型,确定所述目标含钙质致密储层的含水饱和度。
5.根据权利要求3所述的含钙质致密储层含水饱和度的确定方法,其特征在于,所述根据所述岩电分析数据,建立胶结指数、饱和度指数为变动常数的阿尔奇含水饱和度计算模型,包括:
所述岩电分析数据包括:多个岩心试样的实验孔隙度数据、实验地层水电阻率数据、用地层水100%饱和后的电阻率数据、实验深探测电阻率数据、实验含水饱和度数据;
根据所述实验孔隙度数据、所述实验地层水电阻率数据、所述用地层水100%饱和后的电阻率数据,确定实验胶结指数数据和第一岩性系数;
将所述实验孔隙度数据与所述实验胶结指数数据进行交会图分析,确定变动常数胶结指数为:
其中,m为胶结指数;E、G、H分别为第五系数、第六系数、第七系数;φ为孔隙度;φ1为孔隙度分段分界值;
根据所述实验含水饱和度数据、所述实验深探测电阻率数据、所述用地层水100%饱和后的电阻率数据,确定实验饱和度指数数据和第二岩性系数;
将所述实验孔隙度数据与所述实验饱和度指数数据进行交会图分析,确定变动常数饱和度指数:
n=I×φ+J (7)
其中,n为饱和度指数;I、J分别为第八系数和第九系数;
将变动常数胶结指数m和饱和度指数n、以及第一岩性系数a和第二岩性系数b代入公式(8),得到所述阿尔奇含水饱和度计算模型:
其中,Sw为含水饱和度;a、b分别为第一岩性系数和第二岩性系数;RW、RTJ分别为地层水电阻率、校正后的深探测电阻率;φ代表孔隙度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810287787.2A CN110344822B (zh) | 2018-04-03 | 2018-04-03 | 含钙质致密储层含水饱和度的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810287787.2A CN110344822B (zh) | 2018-04-03 | 2018-04-03 | 含钙质致密储层含水饱和度的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110344822A CN110344822A (zh) | 2019-10-18 |
CN110344822B true CN110344822B (zh) | 2022-06-03 |
Family
ID=68172623
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810287787.2A Active CN110344822B (zh) | 2018-04-03 | 2018-04-03 | 含钙质致密储层含水饱和度的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110344822B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111008482B (zh) * | 2019-12-18 | 2023-04-07 | 长江大学 | 一种变质岩储层含水饱和度计算方法及装置 |
CN112800646B (zh) * | 2021-01-08 | 2024-01-23 | 中海石油深海开发有限公司 | 一种修正阿尔奇模型中岩电参数的方法 |
CN115542419A (zh) * | 2022-08-29 | 2022-12-30 | 广州海洋地质调查局 | 海洋天然气水合物钻后综合评价方法、系统和介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1040245A (zh) * | 1988-07-25 | 1990-03-07 | 石油工业部石油勘探开发科学研究院石油地质研究所 | 一种测井信息的解释处理方法 |
FR2661750A1 (fr) * | 1990-05-03 | 1991-11-08 | Petroleo Brasileiro Sa | Procede d'evaluation d'un parametre de diagraphie des sondages. |
GB9924829D0 (en) * | 1998-11-10 | 1999-12-22 | Schlumberger Ltd | Method for interpreting carbonate reservoirs |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6571619B2 (en) * | 2001-10-11 | 2003-06-03 | Schlumberger Technology Corporation | Real time petrophysical evaluation system |
CN103698811B (zh) * | 2013-12-24 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种碳酸盐岩岩石结构组分测井定量识别方法及其用途 |
CN103867197B (zh) * | 2014-04-04 | 2016-07-20 | 中国石油集团川庆钻探工程有限公司 | 复杂岩性天然气层声波时差判别法 |
CN106321092B (zh) * | 2016-09-27 | 2019-04-09 | 中国石油天然气股份有限公司 | 致密储层分类模型建立方法及致密储层分类方法 |
CN107578343B (zh) * | 2017-07-26 | 2020-07-10 | 中国石油天然气股份有限公司 | 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置 |
CN107780923B (zh) * | 2017-11-01 | 2021-04-20 | 中石化石油工程技术服务有限公司 | 一种基于泥质校正的含水饱和度模型的建立、仿真方法 |
-
2018
- 2018-04-03 CN CN201810287787.2A patent/CN110344822B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1040245A (zh) * | 1988-07-25 | 1990-03-07 | 石油工业部石油勘探开发科学研究院石油地质研究所 | 一种测井信息的解释处理方法 |
FR2661750A1 (fr) * | 1990-05-03 | 1991-11-08 | Petroleo Brasileiro Sa | Procede d'evaluation d'un parametre de diagraphie des sondages. |
GB9924829D0 (en) * | 1998-11-10 | 1999-12-22 | Schlumberger Ltd | Method for interpreting carbonate reservoirs |
Non-Patent Citations (6)
Title |
---|
"A油田水淹层测井综合评价方法研究";卜诚;《中国优秀硕士论文全文库工程科技I辑》;20160415(第4期);全文 * |
"华北大王庄复杂砂岩储层钙质校正方法研究";诸葛月英;《测井技术》;20060531(第5期);全文 * |
"基于印度尼西亚公式的饱和度改进模型";柴婧;《石油地质与工程》;20100228;第24卷(第2期);全文 * |
"大庆致密油储层电阻率模型研究";付健;《中国优秀硕士论文全文库工程科技I辑》;20170215(第2期);全文 * |
"罗家地区泥页岩矿物组分含量解释模型构建方法";耿生臣;《油气地质与采收率》;20130131;第20卷(第1期);全文 * |
"致密砂岩导电规律与导电模型实验研究";刘玥;《中国优秀硕士论文全文库工程科技I辑》;20180215(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110344822A (zh) | 2019-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104278991B (zh) | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 | |
CN103670388B (zh) | 一种泥页岩有机碳含量的评价方法 | |
WO2017024530A1 (zh) | 一种计算烃源岩中有机碳含量的方法 | |
CN102640018B (zh) | 用于电缆和lwd测井数据的石油物理和地球物理地层评价的工作流程 | |
CN107917865A (zh) | 一种致密砂岩储层多参数渗透率预测方法 | |
Sun et al. | The theory and application of DEM-Gassmann rock physics model for complex carbonate reservoirs | |
Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
CN109653725A (zh) | 基于沉积微相和岩石相的混积储层水淹程度测井解释方法 | |
CN110344822B (zh) | 含钙质致密储层含水饱和度的确定方法 | |
CN107829731B (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
CN111220522A (zh) | 一种高泥质细粒沉积物水合物饱和度的岩心刻度测井计算方法 | |
CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
CN113216945A (zh) | 一种致密砂岩储层渗透性定量评价方法 | |
CN108412488B (zh) | 快速确定页岩气储层有机孔隙度的测井方法 | |
CN113720745B (zh) | 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法 | |
CN110489769B (zh) | 一种储油层泥质含量计算方法及电子设备 | |
CN116930023A (zh) | 致密砂岩相控分类孔隙度测井精细解释方法和装置 | |
CN112346147A (zh) | 一种基于中子密度孔隙度差的储层评价方法 | |
Li et al. | A novel approach to the quantitative evaluation of the mineral composition, porosity, and kerogen content of shale using conventional logs: A case study of the Damintun Sag in the Bohai Bay Basin, China | |
CN114086938A (zh) | 一种非均质性砂岩储层的含气饱和度预测方法 | |
Amanipoor | Providing a subsurface reservoir quality maps in oil fields by geostatistical methods | |
RU2771802C1 (ru) | Способ дифференциации пустотности неоднородных карбонатных пластов | |
CN110297280A (zh) | 一种碳酸盐岩超压空间分布特征的预测方法 | |
CN117784278B (zh) | 一种致密砂岩气甜点的预测方法和预测系统 | |
CN114458305A (zh) | 孔隙结构系数确定方法及装置 |
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 |