CN114458305A - 孔隙结构系数确定方法及装置 - Google Patents
孔隙结构系数确定方法及装置 Download PDFInfo
- Publication number
- CN114458305A CN114458305A CN202011208849.XA CN202011208849A CN114458305A CN 114458305 A CN114458305 A CN 114458305A CN 202011208849 A CN202011208849 A CN 202011208849A CN 114458305 A CN114458305 A CN 114458305A
- Authority
- CN
- China
- Prior art keywords
- pore
- pore structure
- target area
- resistivity
- determining
- 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
- 239000011148 porous material Substances 0.000 title claims abstract description 291
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000012530 fluid Substances 0.000 claims abstract description 48
- 238000004364 calculation method Methods 0.000 claims abstract description 37
- 239000011435 rock Substances 0.000 claims abstract description 23
- 238000005755 formation reaction Methods 0.000 claims description 41
- 238000010839 reverse transcription Methods 0.000 claims description 41
- 230000015572 biosynthetic process Effects 0.000 claims description 40
- 238000004590 computer program Methods 0.000 claims description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 15
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 12
- 239000011707 mineral Substances 0.000 claims description 12
- 239000008398 formation water Substances 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 8
- 230000003628 erosive effect Effects 0.000 claims description 4
- 239000000706 filtrate Substances 0.000 claims description 4
- 239000011800 void material Substances 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 17
- MQTOSJVFKKJCRP-BICOPXKESA-N azithromycin Chemical compound O([C@@H]1[C@@H](C)C(=O)O[C@@H]([C@@]([C@H](O)[C@@H](C)N(C)C[C@H](C)C[C@@](C)(O)[C@H](O[C@H]2[C@@H]([C@H](C[C@@H](C)O2)N(C)C)O)[C@H]1C)(C)O)CC)[C@H]1C[C@@](C)(OC)[C@@H](O)[C@H](C)O1 MQTOSJVFKKJCRP-BICOPXKESA-N 0.000 description 13
- 229960004099 azithromycin Drugs 0.000 description 13
- 238000002474 experimental method Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 230000005611 electricity Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 235000020681 well water Nutrition 0.000 description 1
- 239000002349 well water 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种孔隙结构系数确定方法及装置，该方法包括：根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型；确定孔隙流体和孔隙结构导致的同类测井信息的差异；根据电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解；根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解；根据解析解和理论解，确定孔喉比的最优解；根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。将孔隙流体和孔隙结构分析与阿奇公式相结合，求取目标区域的孔隙结构系数，相较于理论值，提高了孔隙结构系数的准确性。
Description
技术领域
本发明涉及致密油气藏勘探开发技术领域，尤其涉及一种孔隙结构系数确定方法及装置。
背景技术
含水饱和度指地层中水占岩心孔隙空间的百分比，是评价地层含气性的重要标志。常用经典阿奇公式计算，该公式中的孔隙结构系数是关系到计算准确性的重要参数，在一般地层计算中均选用区域性常量。在复杂碳酸盐岩地层中，如四川盆地高磨地区，多样的地层孔隙结构(孔隙大小组合、连通性)导致地层强烈的非均质性，引起孔隙结构系数发生明显变化，依旧采用经典阿奇公式选取区域性常量计算，就会出现计算值与实验值偏差较大的问题，使得确定的孔隙结构系数准确性低。
发明内容
本发明实施例提供一种孔隙结构系数确定方法，用以提高孔隙结构系数的准确性，该方法包括：
获取目标区域的测井数据和岩电实验数据；
根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型；
根据目标区域的测井数据，确定孔隙流体和孔隙结构导致的同类测井信息的差异；
根据所述电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解；
根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解；
根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解；
根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。
本发明实施例还提供一种孔隙结构系数确定装置，用以提高孔隙结构系数的准确性，该装置包括：
数据获取模块，用于获取目标区域的测井数据和岩电实验数据；
电阻率计算模型确定模块，用于根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型；
同类测井信息差异模块，用于根据目标区域的测井数据，确定孔隙流体和孔隙结构导致的同类测井信息的差异；
解析解确定模块，用于根据所述电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解；
理论解确定模块，用于根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解；
孔喉比最优解确定模块，用于根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解；
孔隙结构系数确定模块，用于根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述孔隙结构系数确定方法。
本发明实施例也提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述孔隙结构系数确定方法的计算机程序。
本发明实施例中，通过获取目标区域的测井数据和岩电实验数据；根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型；根据目标区域的测井数据，确定孔隙流体和孔隙结构导致的同类测井信息的差异；根据电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解；根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解；根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解；根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。通过结合岩电实验和测井多种地球物理手段，将孔隙流体和孔隙结构分析与阿奇公式相结合，求取目标区域的孔隙结构系数，相较于利用阿奇公式和岩电实验数据所求得的理论值，更加贴合实际地层中多样的孔隙结构，从而提高孔隙结构系数的准确性。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中孔隙结构系数确定方法的示意图。
图2为本发明具体实施例中步骤106的具体实现方法示意图。
图3为本发明具体实施例中步骤107的具体实现方法示意图。
图4为本发明具体实施例中单井岩心孔喉比与孔喉加权均值、饱和度中值半径之间的相关关系示意图。
图5为本发明具体实施例中根据单井岩心数据所建立的孔喉比与孔喉加权均值的回归关系拟合示意图。
图6为本发明具体实施例中对δ函数模拟权系数的标定过程示意图。
图7为本发明具体实施例中测井计算孔隙结构系数m与含水饱和度反演的理论解m*之间相关性的标定示意图。
图8为本发明具体实施例中以两种孔隙结构系数所构建的单井含水饱和度计算模型的对比分析图。
图9为本发明实施例中孔隙结构系数确定装置的示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例提供了一种孔隙结构系数确定方法，用以提高孔隙结构系数的准确性，如图1所示，该方法包括：
步骤101：获取目标区域的测井数据和岩电实验数据；
步骤102：根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型；
步骤103：根据目标区域的测井数据，确定孔隙流体和孔隙结构导致的同类测井信息的差异；
步骤104：根据电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解；
步骤105：根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解；
步骤106：根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解；
步骤107：根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。
由图1所示流程可以得知，本发明实施例中，通过获取目标区域的测井数据和岩电实验数据；根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型；根据目标区域的测井数据，确定孔隙流体和孔隙结构导致的同类测井信息的差异；根据电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解；根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解；根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解；根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。通过结合岩电实验和测井多种地球物理手段，将孔隙流体和孔隙结构分析与阿奇公式相结合，求取目标区域的孔隙结构系数，相较于利用阿奇公式和岩电实验数据所求得的理论值，更加贴合实际地层中多样的孔隙结构，从而提高孔隙结构系数的准确性。
具体实施时，首先目标区域的测井数据和岩电实验数据。测井数据是指测井所得的地层参数值，例如泥质含量、孔隙流体和孔隙结构等。岩电实验数据是指实验室中对岩心进行岩电实验所得到的数据，岩电实验作为岩石物理研究的一个重要手段，主要通过测量岩心的孔隙度、电阻率和饱和度等参数来求取阿奇公式中的4个关键参数，进而准确地计算地层含油气饱和度。
获取目标区域的测井数据后，根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型。具体实施例中所建立的电阻率计算模型为：
其中，RT代表深地层电阻率，深地层的深度范围在距离井筒1m以上；RXO代表浅地层电阻率，浅地层的深度范围在距离井筒1m以内；RTmai代表深地层中第i种的矿物电阻率；Vmai代表地层中第i种矿物的体积含量；RTsh代表深地层中泥质的电阻率；Vsh代表地层中泥质的体积含量；RTf代表地层孔隙中流体的电阻率；φf代表地层的孔隙度；
代表孔隙结构对测井值的贡献因子，目前测井信息度量孔隙的范围一般认为从um到mm；σ、υ分别代表孔隙直径与长度；RXOmai代表浅地层中第i种的矿物电阻率；RXOsh代表浅地层中泥质的电阻率；RXOf代表浅地层孔隙中流体的电阻率。
同时，根据目标区域的测井数据，确定孔隙流体和孔隙结构导致的同类测井信息的差异。根据测井解释评价基本理论，孔隙流体和孔隙结构会导致同类测井信息之间产生差异，即：
其中，⊕代表测井曲线与孔隙结构存在一种相互作用，ξ代表同类测井信息之间的差异，如声波中子孔隙度的差异φac/φcnl，地层电阻率差异RT/RXO，纵横波速度比Vp/Vs，斯通利波能量衰减量α等。
且经典阿奇公式为：
其中，SW代表含水饱和度；n代表胶结指数；φ代表地层孔隙度；m代表孔隙结构系数；Rt代表地层电阻率；a、b均为区域常量；Rw代表地层水电阻率。
根据上述电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解，得到目标区域的孔隙结构系数的解析解为：
其中，m代表目标区域的孔隙结构系数的解析解；a、b均为区域常量；Rw代表地层水电阻率；f(Vsh)代表泥质含量；RT代表深地层电阻率，深地层的深度范围在距离井筒1m以上；RXO代表浅地层电阻率，浅地层的深度范围在距离井筒1m以内；Rmf代表泥浆滤液电阻率；n0代表胶结指数，为区域常量；φ代表总孔隙度；δ代表孔隙流体和孔隙结构导致的同类测井信息的差异的总和；λ代表大小孔径比；r代表孔径；rc代表大小孔径分界，也称为孔喉比。
其中，按照如下公式，计算孔隙流体和孔隙结构导致的同类测井信息的差异的总和δ的值：
其中，w1+w2+w3＝1，w1，w2，w3代表不同测井信息的权系数；α代表斯通利波能量衰减量；β代表溶蚀孔洞面孔率；φfv代表裂缝孔隙度；φac代表声波孔隙度；φcnl代表中子孔隙度。
根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解，得到目标区域的孔隙结构系数的理论解为：
其中，m*代表目标区域的孔隙结构系数的理论解；a、b均为区域常量；Rw代表地层水电阻率；C(Sw)代表岩心含水饱和度；RT*代表岩心取样深度段的地层电阻率；φ*代表岩心取样深度段的孔隙度。
根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解，如图2所示，包括：
步骤201：联立目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，得到大小孔径比λ的理论值；
步骤202：根据大小孔径比λ的理论值，确定孔喉比的最优解。
具体实施时：
已知
可求取孔喉比rc的最优解。
确定孔喉比的最优解后，根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数，如图3所示，包括：
步骤301：根据孔喉比的最优解和岩电实验数据中的岩心孔喉加权均值，确定孔喉比与孔喉加权均值的关联关系；
步骤302：根据孔喉比与孔喉加权均值的关联关系和目标区域的测井数据中的孔喉加权均值测井数据，确定孔喉比的测井数值；
步骤303：根据孔喉比与孔喉加权均值的关联关系，对孔隙流体和孔隙结构导致的同类测井信息的差异的总和δ的权系数w1，w2，w3进行标定，确定δ的值；
步骤304：根据孔喉比的测井数值、δ的值和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。
步骤301具体实施时，在进行多次试验后，发现孔喉比rc与岩心孔喉加权均值ravg有很好的相关关系。对单井岩心数据为例，如图4所示，孔喉比rc与岩心孔喉加权均值ravg的相关系数0.962，要明显高于孔喉比rc与饱和度中值半径的相关系数0.798。
孔喉比rc与孔喉加权均值ravg有正相关关系，可以表达为：rc＝ρ×ravg τ，如图5所示，根据单井岩心数据所建立的孔喉比与孔喉加权均值的关联关系为：rc＝1.169×ravg 0.9573。
根据上述关联关系、目标区域的测井数据中的孔喉加权均值测井数据，就能够确定孔喉比的测井数值。
步骤303具体实施时，如图6所示，对权系数w1，w2，w3进行多次标定，第一组权系数W1＝0.6，W2＝0.2，W3＝0.3得到的相关系数为0.92，第二组权系数W1＝0.2，W2＝0.6，W3＝0.3得到的相关系数为0.932，第三组权系数W1＝0.2，W2＝0.2，W3＝0.6得到的相关系数为0.95，第四组权系数W1＝0，W2＝0，W3＝1得到的相关系数为0.962。选取能够使得岩心孔喉加权均值ravg与最优孔喉比(也可称为最优大小孔径分界)rc的相关系数最高的权系数取值，得到W1＝0，W2＝0，W3＝1，从而能够确定δ的值。
将孔喉比的测井数值、δ的值代入目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。如图7所示，为本发明具体实施例中，根据本发明提供的方法确定的测井计算孔隙结构系数m，以及根据岩心含水饱和度反演得到的孔隙结构系数的理论解m*之间的标定结果。
具体实施时，得到目标区域的孔隙结构系数后，还可以利用目标区域的孔隙结构系数和经典阿奇公式，求得含水饱和度SW。如图8所示，为利用本发明实施例确定的孔隙结构系数求得的含水饱和度，与根据阿奇公式求得的含水饱和度之间的比较，其中，经典阿奇公式中的岩电参数取值为：a＝1，b＝1，m＝2.28，n＝1.79；Rw≈0.021Ω·m(130°)。
从图8的含水饱和度SW标定的结果可知，本发明具体实例中所建立的含水饱和度模型与岩心含水饱和度匹配更好，相关系数更高，证明本发明实施例所确定的孔隙结构系数的准确性更高。
基于同一发明构思，本发明实施例还提供一种孔隙结构系数确定装置，由于孔隙结构系数确定装置所解决问题的原理与孔隙结构系数确定方法相似，因此孔隙结构系数确定装置的实施可以参见孔隙结构系数确定方法的实施，重复之处不再赘述，具体结构如图9所示：
数据获取模块901，用于获取目标区域的测井数据和岩电实验数据；
电阻率计算模型确定模块902，用于根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型；
同类测井信息差异模块903，用于根据目标区域的测井数据，确定孔隙流体和孔隙结构导致的同类测井信息的差异；
解析解确定模块904，用于根据电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解；
理论解确定模块905，用于根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解；
孔喉比最优解确定模块906，用于根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解；
孔隙结构系数确定模块907，用于根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。
具体实施例中，电阻率计算模型为：
其中，RT代表深地层电阻率，深地层的深度范围在距离井筒1m以上；RXO代表浅地层电阻率，浅地层的深度范围在距离井筒1m以内；RTmai代表深地层中第i种的矿物电阻率；Vmai代表地层中第i种矿物的体积含量；RTsh代表深地层中泥质的电阻率；Vsh代表地层中泥质的体积含量；RTf代表地层孔隙中流体的电阻率；φf代表地层的孔隙度；
代表孔隙结构对测井值的贡献因子，目前测井信息度量孔隙的范围一般认为从um到mm；σ、υ分别代表孔隙直径与长度；RXOmai代表浅地层中第i种的矿物电阻率；RXOsh代表浅地层中泥质的电阻率；RXOf代表浅地层孔隙中流体的电阻率。
具体实施时，目标区域的孔隙结构系数的解析解为：
其中，m代表目标区域的孔隙结构系数的解析解；a、b均为区域常量；Rw代表地层水电阻率；f(Vsh)代表泥质含量；RT代表深地层电阻率，深地层的深度范围在距离井筒1m以上；RXO代表浅地层电阻率，浅地层的深度范围在距离井筒1m以内；Rmf代表泥浆滤液电阻率；n0代表胶结指数，为区域常量；φ代表总孔隙度；δ代表孔隙流体和孔隙结构导致的同类测井信息的差异的总和；λ代表大小孔径比；r代表孔径；rc代表大小孔径分界。
其中，按照如下公式，计算孔隙流体和孔隙结构导致的同类测井信息的差异的总和δ的值：
其中，w1+w2+w3＝1，w1，w2，w3代表不同测井信息的权系数；α代表斯通利波能量衰减量；β代表溶蚀孔洞面孔率；φfv代表裂缝孔隙度；φac代表声波孔隙度；φcnl代表中子孔隙度。
具体实施例中，目标区域的孔隙结构系数的理论解为：
其中，m*代表目标区域的孔隙结构系数的理论解；a、b均为区域常量；Rw代表地层水电阻率；C(Sw)代表岩心含水饱和度；RT*代表岩心取样深度段的地层电阻率；φ*代表岩心取样深度段的孔隙度。
具体实施时，孔喉比最优解确定模块906具体用于：
联立目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，得到大小孔径比的理论值；
根据大小孔径比的理论值，确定孔喉比的最优解。
具体实施例中，孔隙结构系数确定模块907具体用于：
根据孔喉比的最优解和岩电实验数据中的岩心孔喉加权均值，确定孔喉比与孔喉加权均值的关联关系；
根据孔喉比与孔喉加权均值的关联关系和目标区域的测井数据中的孔喉加权均值测井数据，确定孔喉比的测井数值；
根据孔喉比与孔喉加权均值的关联关系，对孔隙流体和孔隙结构导致的同类测井信息的差异的总和δ的权系数w1，w2，w3进行标定，确定δ的值；
根据孔喉比的测井数值、δ的值和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述孔隙结构系数确定方法。
本发明实施例还提供一种计算机可读存储介质，存储有执行上述孔隙结构系数确定方法的计算机程序。
综上所述，本发明实施例提供的孔隙结构系数确定方法及装置具有如下优点：
通过获取目标区域的测井数据和岩电实验数据；根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型；根据目标区域的测井数据，确定孔隙流体和孔隙结构导致的同类测井信息的差异；根据电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解；根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解；根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解；根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。通过结合岩电实验和测井多种地球物理手段，将孔隙流体和孔隙结构分析与阿奇公式相结合，求取目标区域的孔隙结构系数，相较于利用阿奇公式和岩电实验数据所求得的理论值，更加贴合实际地层中多样的孔隙结构，从而提高孔隙结构系数的准确性。
本领域内的技术人员应明白，本发明的实施例可提供为方法、装置、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、装置和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (16)
1.一种孔隙结构系数确定方法，其特征在于，包括：
获取目标区域的测井数据和岩电实验数据；
根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型；
根据目标区域的测井数据，确定孔隙流体和孔隙结构导致的同类测井信息的差异；
根据所述电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解；
根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解；
根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解；
根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。
2.如权利要求1所述的方法，其特征在于，所述电阻率计算模型为：
其中，RT代表深地层电阻率，深地层的深度范围在距离井筒1m以上；RXO代表浅地层电阻率，浅地层的深度范围在距离井筒1m以内；RTmai代表深地层中第i种的矿物电阻率；Vmai代表地层中第i种矿物的体积含量；RTsh代表深地层中泥质的电阻率；Vsh代表地层中泥质的体积含量；RTf代表地层孔隙中流体的电阻率；φf代表地层的孔隙度；
代表孔隙结构对测井值的贡献因子，目前测井信息度量孔隙的范围一般认为从um到mm；σ、υ分别代表孔隙直径与长度；RXOmai代表浅地层中第i种的矿物电阻率；RXOsh代表浅地层中泥质的电阻率；RXOf代表浅地层孔隙中流体的电阻率。
3.如权利要求1所述的方法，其特征在于，所述目标区域的孔隙结构系数的解析解为：
其中，m代表目标区域的孔隙结构系数的解析解；a、b均为区域常量；Rw代表地层水电阻率；f(Vsh)代表泥质含量；RT代表深地层电阻率，深地层的深度范围在距离井筒1m以上；RXO代表浅地层电阻率，浅地层的深度范围在距离井筒1m以内；Rmf代表泥浆滤液电阻率；n0代表胶结指数，为区域常量；φ代表总孔隙度；δ代表孔隙流体和孔隙结构导致的同类测井信息的差异的总和；λ代表大小孔径比；r代表孔径；rc代表大小孔径分界。
4.如权利要求3所述的方法，其特征在于，按照如下公式，计算孔隙流体和孔隙结构导致的同类测井信息的差异的总和δ的值：
其中，w1+w2+w3＝1，w1，w2，w3代表不同测井信息的权系数；α代表斯通利波能量衰减量；β代表溶蚀孔洞面孔率；φfv代表裂缝孔隙度；φac代表声波孔隙度；φcnl代表中子孔隙度。
5.如权利要求4所述的方法，其特征在于，根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数，包括：
根据孔喉比的最优解和岩电实验数据中的岩心孔喉加权均值，确定孔喉比与孔喉加权均值的关联关系；
根据孔喉比与孔喉加权均值的关联关系和目标区域的测井数据中的孔喉加权均值测井数据，确定孔喉比的测井数值；
根据孔喉比与孔喉加权均值的关联关系，对孔隙流体和孔隙结构导致的同类测井信息的差异的总和δ的权系数w1，w2，w3进行标定，确定δ的值；
根据孔喉比的测井数值、δ的值和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。
6.如权利要求1所述的方法，其特征在于，所述目标区域的孔隙结构系数的理论解为：
其中，m*代表目标区域的孔隙结构系数的理论解；a、b均为区域常量；Rw代表地层水电阻率；C(Sw)代表岩心含水饱和度；RT*代表岩心取样深度段的地层电阻率；φ*代表岩心取样深度段的孔隙度。
7.如权利要求1所述的方法，其特征在于，根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解，包括：
联立目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，得到大小孔径比的理论值；
根据大小孔径比的理论值，确定孔喉比的最优解。
8.一种孔隙结构系数确定装置，其特征在于，包括：
数据获取模块，用于获取目标区域的测井数据和岩电实验数据；
电阻率计算模型确定模块，用于根据目标区域的测井数据，建立岩石骨架、泥质含量、孔隙流体和孔隙结构作为影响因子的电阻率计算模型；
同类测井信息差异模块，用于根据目标区域的测井数据，确定孔隙流体和孔隙结构导致的同类测井信息的差异；
解析解确定模块，用于根据所述电阻率计算模型、同类测井信息的差异和经典阿奇公式，确定目标区域的孔隙结构系数的解析解；
理论解确定模块，用于根据目标区域的岩电实验数据，确定目标区域的孔隙结构系数的理论解；
孔喉比最优解确定模块，用于根据目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，确定孔喉比的最优解；
孔隙结构系数确定模块，用于根据孔喉比的最优解和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。
9.如权利要求8所述的装置，其特征在于，所述电阻率计算模型为：
其中，RT代表深地层电阻率，深地层的深度范围在距离井筒1m以上；RXO代表浅地层电阻率，浅地层的深度范围在距离井筒1m以内；RTmai代表深地层中第i种的矿物电阻率；Vmai代表地层中第i种矿物的体积含量；RTsh代表深地层中泥质的电阻率；Vsh代表地层中泥质的体积含量；RTf代表地层孔隙中流体的电阻率；φf代表地层的孔隙度；
代表孔隙结构对测井值的贡献因子，目前测井信息度量孔隙的范围一般认为从um到mm；σ、υ分别代表孔隙直径与长度；RXOmai代表浅地层中第i种的矿物电阻率；RXOsh代表浅地层中泥质的电阻率；RXOf代表浅地层孔隙中流体的电阻率。
10.如权利要求8所述的装置，其特征在于，所述目标区域的孔隙结构系数的解析解为：
其中，m代表目标区域的孔隙结构系数的解析解；a、b均为区域常量；Rw代表地层水电阻率；f(Vsh)代表泥质含量；RT代表深地层电阻率，深地层的深度范围在距离井筒1m以上；RXO代表浅地层电阻率，浅地层的深度范围在距离井筒1m以内；Rmf代表泥浆滤液电阻率；n0代表胶结指数，为区域常量；φ代表总孔隙度；δ代表孔隙流体和孔隙结构导致的同类测井信息的差异的总和；λ代表大小孔径比；r代表孔径；rc代表大小孔径分界。
11.如权利要求10所述的装置，其特征在于，按照如下公式，计算孔隙流体和孔隙结构导致的同类测井信息的差异的总和δ的值：
其中，w1+w2+w3＝1，w1，w2，w3代表不同测井信息的权系数；α代表斯通利波能量衰减量；β代表溶蚀孔洞面孔率；φfv代表裂缝孔隙度；φac代表声波孔隙度；φcnl代表中子孔隙度。
12.如权利要求11所述的装置，其特征在于，所述孔隙结构系数确定模块具体用于：
根据孔喉比的最优解和岩电实验数据中的岩心孔喉加权均值，确定孔喉比与孔喉加权均值的关联关系；
根据孔喉比与孔喉加权均值的关联关系和目标区域的测井数据中的孔喉加权均值测井数据，确定孔喉比的测井数值；
根据孔喉比与孔喉加权均值的关联关系，对孔隙流体和孔隙结构导致的同类测井信息的差异的总和δ的权系数w1，w2，w3进行标定，确定δ的值；
根据孔喉比的测井数值、δ的值和目标区域的孔隙结构系数的解析解，确定目标区域的孔隙结构系数。
13.如权利要求8所述的装置，其特征在于，所述目标区域的孔隙结构系数的理论解为：
其中，m*代表目标区域的孔隙结构系数的理论解；a、b均为区域常量；Rw代表地层水电阻率；C(Sw)代表岩心含水饱和度；RT*代表岩心取样深度段的地层电阻率；φ*代表岩心取样深度段的孔隙度。
14.如权利要求8所述的装置，其特征在于，所述孔喉比最优解确定模块具体用于：
联立目标区域的孔隙结构系数的解析解和目标区域的孔隙结构系数的理论解，得到大小孔径比的理论值；
根据大小孔径比的理论值，确定孔喉比的最优解。
15.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至7任一所述方法。
16.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至7任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011208849.XA CN114458305A (zh) | 2020-11-03 | 2020-11-03 | 孔隙结构系数确定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011208849.XA CN114458305A (zh) | 2020-11-03 | 2020-11-03 | 孔隙结构系数确定方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114458305A true CN114458305A (zh) | 2022-05-10 |
Family
ID=81404209
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011208849.XA Pending CN114458305A (zh) | 2020-11-03 | 2020-11-03 | 孔隙结构系数确定方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114458305A (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090292472A1 (en) * | 2004-12-17 | 2009-11-26 | Schlumberger Technology Corporation | Method for the Characterization of Geological Formations |
| CN101930082A (zh) * | 2009-06-24 | 2010-12-29 | 中国石油集团川庆钻探工程有限公司 | 采用电阻率资料进行储层流体类型判别的方法 |
| CN103278436A (zh) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
| US20180274364A1 (en) * | 2015-12-30 | 2018-09-27 | Petrochina Company Limited | Method and Device for Determining Karst Development Degree of Reservoir, Computer Readable Storage Medium and Device |
| CN110046414A (zh) * | 2019-04-03 | 2019-07-23 | 长江大学 | 用于识别优势渗流段的储层参数模型构建方法及设备 |
| CN111562629A (zh) * | 2019-02-13 | 2020-08-21 | 中国石油天然气股份有限公司 | 一种基于等效孔隙截面指数的饱和度确定方法及装置 |
-
2020
- 2020-11-03 CN CN202011208849.XA patent/CN114458305A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090292472A1 (en) * | 2004-12-17 | 2009-11-26 | Schlumberger Technology Corporation | Method for the Characterization of Geological Formations |
| CN101930082A (zh) * | 2009-06-24 | 2010-12-29 | 中国石油集团川庆钻探工程有限公司 | 采用电阻率资料进行储层流体类型判别的方法 |
| CN103278436A (zh) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
| US20180274364A1 (en) * | 2015-12-30 | 2018-09-27 | Petrochina Company Limited | Method and Device for Determining Karst Development Degree of Reservoir, Computer Readable Storage Medium and Device |
| CN111562629A (zh) * | 2019-02-13 | 2020-08-21 | 中国石油天然气股份有限公司 | 一种基于等效孔隙截面指数的饱和度确定方法及装置 |
| CN110046414A (zh) * | 2019-04-03 | 2019-07-23 | 长江大学 | 用于识别优势渗流段的储层参数模型构建方法及设备 |
Non-Patent Citations (5)
| Title |
|---|
| TIANZHI TANG ET AL.: "Electrical characteristics of rocks in fractured and caved reservoirs", JOURNAL OF GEOPHYSICS AND ENGINEERING, vol. 14, no. 6, 31 December 2017 (2017-12-31), pages 1437 - 1444, XP020322157, DOI: 10.1088/1742-2140/aa8528 * |
| 李潮流;李长喜;: "特低渗透率砂岩储集层电学性质研究", 测井技术, no. 03, 20 June 2010 (2010-06-20), pages 233 - 237 * |
| 王海涛;张晋言;: "基于孔喉腔模型研究孔隙结构对于多孔介质孔隙度指数的影响", 地球物理学进展, no. 02, 15 April 2019 (2019-04-15), pages 218 - 228 * |
| 陈杰: "基于电阻率测井资料研究致密砂岩孔隙结构特征", 西南石油学院硕士学位论文集, 31 December 2006 (2006-12-31) * |
| 陈杰;刘向君;成竹;杨旭明;刘洪;周改英;: "利用电阻率测井资料研究砂岩孔隙结构特征", 西南石油学院学报, no. 06, 30 December 2005 (2005-12-30), pages 5 - 7 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111425193B (zh) | 一种基于聚类分析测井岩石物理相划分的储层可压性评价方法 | |
| CN108897066B (zh) | 碳酸盐岩裂缝密度定量预测方法及装置 | |
| CN112526107B (zh) | 裂缝型致密砂岩储层甜点识别与定量表征的方法 | |
| CN107808032B (zh) | 页岩气储层含水饱和度计算方法及系统 | |
| US20130282286A1 (en) | System and method for calibrating permeability for use in reservoir modeling | |
| CN106355571B (zh) | 一种白云岩储层质量的确定方法及装置 | |
| CN112149282A (zh) | 一种井中天然气水合物饱和度岩石物理计算方法及系统 | |
| US11435494B1 (en) | Method for predicting oil accumulation depth limit of deep and ultra-deep marine carbonate reservoirs | |
| CN112922591B (zh) | 泥页岩储层岩相“甜点”预测方法及系统 | |
| CN106323835B (zh) | 确定非均质碳酸盐岩储层胶结指数的方法 | |
| CN110344822B (zh) | 含钙质致密储层含水饱和度的确定方法 | |
| CN115390155A (zh) | 针对非均质储层的测井解释方法、装置、电子设备及介质 | |
| CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
| CN105401937B (zh) | 一种基于孔隙结构的饱和度指数预测方法 | |
| CN114458305A (zh) | 孔隙结构系数确定方法及装置 | |
| RU2479714C1 (ru) | Способ получения трехмерного распределения проницаемости пласта | |
| CN112505154B (zh) | 泥页岩储层矿物成分含量解析与岩相识别表征方法 | |
| CN112946739B (zh) | 裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法 | |
| CN116930023A (zh) | 致密砂岩相控分类孔隙度测井精细解释方法和装置 | |
| CN115144929A (zh) | 油气藏储量确定方法、装置、存储介质及电子设备 | |
| CN114565116A (zh) | 一种低渗储层物性参数的预测方法和装置 | |
| CN112230278A (zh) | 渗流场特征参数确定方法及装置 | |
| CN113622903B (zh) | 储层改造方法及装置 | |
| CN112946754B (zh) | 储层孔隙度预测方法及装置 | |
| CN111751900B (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 |