CN112049623B - 一种中基性火山岩束缚水出水的判定方法和装置 - Google Patents
一种中基性火山岩束缚水出水的判定方法和装置 Download PDFInfo
- Publication number
- CN112049623B CN112049623B CN201910488394.2A CN201910488394A CN112049623B CN 112049623 B CN112049623 B CN 112049623B CN 201910488394 A CN201910488394 A CN 201910488394A CN 112049623 B CN112049623 B CN 112049623B
- Authority
- CN
- China
- Prior art keywords
- nuclear magnetic
- porosity
- volcanic
- movable
- target
- 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 200
- 238000000034 method Methods 0.000 title claims abstract description 48
- 230000005291 magnetic effect Effects 0.000 claims abstract description 184
- 239000011435 rock Substances 0.000 claims abstract description 90
- 238000005259 measurement Methods 0.000 claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims description 44
- 239000012530 fluid Substances 0.000 claims description 41
- 238000010586 diagram Methods 0.000 claims description 31
- 238000005481 NMR spectroscopy Methods 0.000 claims description 30
- 238000001228 spectrum Methods 0.000 claims description 24
- 238000009826 distribution Methods 0.000 claims description 21
- 230000005311 nuclear magnetism Effects 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000005315 distribution function Methods 0.000 claims description 8
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 7
- 239000008398 formation water Substances 0.000 claims description 7
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims description 6
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 claims description 6
- 229910052710 silicon Inorganic materials 0.000 claims description 6
- 239000010703 silicon Substances 0.000 claims description 6
- 238000005094 computer simulation Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 8
- 239000003921 oil Substances 0.000 description 12
- 238000004590 computer program Methods 0.000 description 10
- 229910052500 inorganic mineral Inorganic materials 0.000 description 8
- 239000011707 mineral Substances 0.000 description 8
- 238000012360 testing method Methods 0.000 description 8
- 230000005298 paramagnetic effect Effects 0.000 description 7
- 239000011148 porous material Substances 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000010936 titanium Substances 0.000 description 1
- 229910052719 titanium Inorganic materials 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
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请公开了一种中基性火山岩束缚水出水的判定方法和装置。该方法包括：获取目标岩心的测量结果，其中，目标岩心为待检测中基性火山岩的岩心，测量结果为对目标岩心进行核磁孔隙度和T2截止值测量所得到的结果数据；依据测量结果确定目标岩心的核磁可动束缚水饱和度和核磁孔隙度；依据预设关系式对目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定待检测中基性火山岩是否出水，其中，预设关系式用于表示核磁可动束缚水饱和度和核磁孔隙度之间的线性关系。通过本申请，解决了相关技术中无法对火山岩等复杂地层进行精准勘探的问题。
Description
技术领域
本申请涉及石油勘探领域，具体而言，涉及一种中基性火山岩束缚水出水的判定方法和装置。
背景技术
在油气勘探中，流体类型判别是测井解释中的一个重点参数，地层是否出水与油藏认识及后期开发息息相关。
随着油气勘探的不断深入，火山岩等复杂地层逐渐成为勘探重点。但是，与碎屑岩油气藏相比，火山岩储集层一般埋藏深度大，且岩石类型变化多样，非均质性强，火山岩储层的孔隙结构更为复杂。而上述因素均对电阻率的变化程度有不可忽略的影响，因此，以电阻率法为依托的流体性质判别方法并不适用火山岩等复杂地形的勘探。
即，目前国内外的常规流体类型判别方法对于火山岩储层判别误差较大，并不适用。
针对相关技术中无法对火山岩等复杂地层进行精准勘探的问题，目前尚未提出有效的解决方案。
发明内容
本申请提供一种中基性火山岩束缚水出水的判定方法和装置，以解决相关技术中无法对火山岩等复杂地层进行精准勘探的问题。
根据本申请的一个方面，提供了一种中基性火山岩束缚水出水的判定方法。该方法包括：获取目标岩心的测量结果，其中，所述目标岩心为待检测中基性火山岩的岩心，所述测量结果为对所述目标岩心进行核磁孔隙度和T2截止值测量所得到的结果数据；依据所述测量结果确定所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度；依据预设关系式对所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定所述待检测中基性火山岩是否出水，其中，所述预设关系式用于表示所述核磁可动束缚水饱和度和所述核磁孔隙度之间的线性关系。
可选的，依据预设关系式对所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定所述待检测中基性火山岩是否出水包括：依据所述预设关系式建立火山岩地层出水界限模型，其中，所述火山岩地层出水界限模型为以所述核磁可动束缚水饱和度和核磁孔隙度为坐标轴的二维坐标系模型；将所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度带入至所述火山岩地层出水界限模型中，以确定所述待检测中基性火山岩是否出水；其中，在所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度位于所述火山岩地层出水界限模型中预设关系式第一侧的情况下，则确定所述目标岩心出水；在所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度位于所述火山岩束缚水出水界限模型中预设关系式第二侧的情况下，则确定所述目标岩心未出水。
可选的，所述预设关系式为：φt＝h*Sfswi+k，其中，h、k为预设系数，φt为核磁孔隙度，Sfswi为核磁可动束缚水饱和度。
可选的，依据所述测量结果确定所述目标岩心的核磁可动束缚水饱和度为：依据所述目标岩心的核磁孔隙度和T2截止值测量所得到的结果数据，确定所述目标岩心的核磁共振T2谱分布图；依据所述目标岩心的核磁共振T2谱分布图，确定所述目标岩心的可动流体饱和度和几何平均值；确定所述核磁可动束缚水饱和度的计算模型，并将所述目标岩心的可动流体饱和度和几何平均值带入所述计算模型中，得到所述目标岩心的核磁可动束缚水饱和度。
可选的，所述核磁可动束缚水饱和度的计算模型为：
Sfswi＝a*loge(T2lm)+b*Sff+c，
其中，a、b、d为预设系数，Sfswi为核磁可动束缚水饱和度，T2lm为几何平均值，Sff为可动流体饱和度。
可选的，依据所述目标岩心的核磁共振T2谱分布图，确定所述目标岩心的可动流体饱和度包括：确定所述可动流体饱和度的计算公式；结合所述目标岩心的核磁共振T2谱分布图，通过所述可动流体饱和度的计算公式，确定所述目标岩心的可动流体饱和度，其中，所述可动流体饱和度的计算公式为：
Sff为可动流体饱和度，T2max为最大横向弛豫时间，ATS为最小横向弛豫时间，S(T)为与横向弛豫时间相关的孔隙度分布函数，T2c为T2截止值，Wsi为ECS测井得到地层硅元素的百分含量，Wfe为ECS测井得到地层铁元素的百分含量，Wti为ECS测井得到地层钛元素的百分含量，d、e、f、g为预设参数。
可选的，依据所述测量结果确定所述目标岩心的核磁孔隙度为：依据所述目标岩心的核磁孔隙度和T2截止值测量所得到的结果数据，确定所述目标岩心的核磁共振T2谱分布图；确定所述核磁孔隙度的计算模型；结合所述目标岩心的核磁共振T2谱分布图，通过所述核磁孔隙度的计算模型，确定所述目标岩心的核磁孔隙度，其中，所述核磁孔隙度的计算公式为：φt为核磁孔隙度，T2max为最大横向弛豫时间，ATS为最小横向弛豫时间，S(T)为与横向弛豫时间相关的孔隙度分布函数。
根据本申请的另一方面，提供了一种中基性火山岩束缚水出水的判定装置。该装置包括：获取单元，用于获取目标岩心的测量结果，其中，所述目标岩心为待检测中基性火山岩的岩心，所述测量结果为对所述目标岩心进行核磁孔隙度和T2截止值测量所得到的结果数据；确定单元，用于依据所述测量结果确定所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度；判断单元，用于依据预设关系式对所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定所述待检测中基性火山岩是否出水，其中，所述预设关系式用于表示所述核磁可动束缚水饱和度和所述核磁孔隙度之间的线性关系。
通过本申请，采用对目标岩心进行核磁孔隙度和T2截止值测量以得到测量结果，此时基于获取到的测量结果进一步确定目标岩心的核磁可动束缚水饱和度和核磁孔隙度，此时，可以预设关系式对上述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定待检测中基性火山岩是否出水。解决了相关技术中相关技术中无法对火山岩等复杂地层进行精准勘探的问题。进而达到了对岩石类型多样性高、非均质性强、孔隙结构复杂性高的复杂地层(如：火山岩)进行精准勘探的效果。
附图说明
构成本申请的一部分的附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
图1是根据本申请实施例提供的中基性火山岩束缚水出水的判定方法的流程图；
图2是根据本申请实施例提供的一种可选的火山岩地层出水界限模型的示意图；
图3是根据本申请实施例提供的一种可选的核磁共振T2谱分布图；
图4是根据本申请实施例提供的核磁可动束缚饱和度的比对示意图；
图5是根据本申请实施例提供的T2截止量的对比示意图；
图6是根据本申请实施例提供的一种可选的地层测量数据的示意图；以及
图7是根据本申请实施例提供的中基性火山岩束缚水出水的判定装置的示意图。
具体实施方式
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本申请的实施例。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
根据本申请的实施例，提供了一种中基性火山岩束缚水出水的判定方法。
图1是根据本申请实施例的中基性火山岩束缚水出水的判定方法的流程图。如图1所示，该方法包括以下步骤：
步骤S102，获取目标岩心的测量结果，其中，所述目标岩心为待检测中基性火山岩的岩心，所述测量结果为对所述目标岩心进行核磁孔隙度和T2截止值测量所得到的结果数据。
步骤S104，依据所述测量结果确定所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度。
步骤S106，依据预设关系式对所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定所述待检测中基性火山岩是否出水，其中，所述预设关系式用于表示所述核磁可动束缚水饱和度和所述核磁孔隙度之间的线性关系。
需要说明的是：上述预设关系式为申请人基于大量的中基性火山岩岩心样品的岩石物理实验数据、地层元素测井数据、以及核磁测井数据建立出的数学模型所得出的。也即，申请人通过上述大量数据(研究区试油资料)分析出：可以通过该中基性火山岩的核磁束缚水饱和度和核磁孔隙度之间的相关关系，判断中基性火山岩是否出水时。同时，申请人通过上述大量数据(研究区试油资料)中分析得出，上述中基性火山岩的核磁束缚水饱和度和核磁孔隙度之间的具体相关关系式，即预设关系式。
需要说明的是：上述研究区试油资料中至少包含：具有代表性的火山岩岩心样品24块的测量结果，其中，该测量结果是按照《岩心分析方法(SY/T5336-2006)》和《岩样核磁共振参数实验室测量规范(SY/T6490-2014)》标准流程，对上述具有代表性的火山岩岩心样品24块进行核磁共振实验所得。
还需要说明的是：上述研究区试油资料至少包含：气层、油层、油水同层、水层。
也即，本申请实施例提供的中基性火山岩束缚水出水的判定方法，通过对目标岩心进行核磁孔隙度和T2截止值测量以得到测量结果，此时基于获取到的测量结果进一步确定目标岩心的核磁可动束缚水饱和度和核磁孔隙度，此时，可以预设关系式对上述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定待检测中基性火山岩是否出水。解决了相关技术中相关技术中无法对火山岩等复杂地层进行精准勘探的问题。进而达到了对岩石类型多样性高、非均质性强、孔隙结构复杂性高的复杂地层(如：火山岩)进行精准勘探的效果。
还需要说明的是：步骤S102还可通过如下步骤得以实现：连接测量装置，并从测量装置中读取目标岩心的测量结果，其中，上述测量装置用于对岩心等物质进行分析测量。
可选的，测量装置可以预先将测量结果存储于存储器中，此时，步骤102还可通过如下步骤得以实现：从指定存储器中读取目标岩心的测量结果，其中，指定存储器中存储有测量装置对目标岩心进行测量处理后所得的测量结果。
可选的，在本申请实施例提供的中基性火山岩束缚水出水的判定方法中，依据预设关系式对所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定所述待检测中基性火山岩是否出水包括：依据所述预设关系式建立火山岩地层出水界限模型，其中，所述火山岩地层出水界限模型为以所述核磁可动束缚水饱和度和核磁孔隙度为坐标轴的二维坐标系模型；将所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度带入至所述火山岩地层出水界限模型中，以确定所述待检测中基性火山岩是否出水；其中，在所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度位于所述火山岩地层出水界限模型中预设关系式第一侧的情况下，则确定所述目标岩心出水；在所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度位于所述火山岩束缚水出水界限模型中预设关系式第二侧的情况下，则确定所述目标岩心未出水。
在一个可选的示例中，所述预设关系式为：φt＝h*Sfswi+k，其中，h、k为预设系数，φt为核磁孔隙度，Sfswi为核磁可动束缚水饱和度。
举例说明：如图2所示，建立以核磁孔隙度为纵坐标，以核磁可动束缚水饱和度为横坐标的火山岩地层出水界限模型，其中，h选用0.95，k选用13.8。此时，岩层数据位于上述火山岩地层出水界限模型中的直线关系式之上时，该岩层会出水(包含水层和/或水油同层)。
需要说明的是：上述核磁可动束缚水饱和度的单位选用％。
可选的，在本申请实施例提供的中基性火山岩束缚水出水的判定方法中，所述目标岩心的核磁可动束缚水饱和度可以通过如下方式确定：步骤A1，依据所述目标岩心的核磁孔隙度和T2截止值测量所得到的结果数据，确定所述目标岩心的核磁共振T2谱分布图；步骤A2，依据所述目标岩心的核磁共振T2谱分布图，确定所述目标岩心的可动流体饱和度和几何平均值；步骤A3，确定所述核磁可动束缚水饱和度的计算模型，并将所述目标岩心的可动流体饱和度和几何平均值带入所述计算模型中，得到所述目标岩心的核磁可动束缚水饱和度。
在一个可选的示例中，所述核磁可动束缚水饱和度的计算模型为：
Sfswi＝a*loge(T2lm)+b*Sff+c，其中，a、b、d为预设系数，Sfswi为核磁可动束缚水饱和度，T2lm为几何平均值，Sff为可动流体饱和度。
需要说明的是：如图3所示，核磁可动束缚水饱和度即为离心核磁共振T2谱分布图中大于T2截止值的孔隙占总孔隙体积。经研究发现：核磁可动束缚水饱和度与核磁共振T2谱分布图中的可动流体饱和度和几何平均值具有很好相关性，并建立核磁可动束缚水饱和度的计算模型，其中，该核磁可动束缚水饱和度的计算模型如下：
Sfswi＝2.337*loge(T2lm)+0.207*Sff+7。
进一步的，图4为由上述模型计算所得的核磁可动束缚水饱和度，与实验分析所得的核磁可动束缚水饱和度的对比图，如图4所示，模型计算所得的核磁可动束缚水饱和度与实验分析所得的核磁可动束缚水饱和度一致性好，计算精度较高。即采用本申请实施例所提供的核磁可动束缚水饱和度的计算模型，可以计算出来高精准度的核磁可动束缚水饱和度，进而达到精准判断岩层是否出水。
进一步的，上述目标岩心的可动流体饱和度可以通过如下步骤确定：步骤a1，确定所述可动流体饱和度的计算公式；步骤a2，结合所述目标岩心的核磁共振T2谱分布图，通过所述可动流体饱和度的计算公式，确定所述目标岩心的可动流体饱和度，其中，所述可动流体饱和度的计算公式为：
Sff为可动流体饱和度，T2max为最大横向弛豫时间，ATS为最小横向弛豫时间，S(T)为与横向弛豫时间相关的孔隙度分布函数，T2c为T2截止值。
在一个可选的示例中，T2c＝d*Wsi+e*Wfe+f*Wti+g，将该公式代入上述可动流体饱和度的计算公式中，可得：其中，Wsi为ECS测井得到地层硅元素的百分含量，Wfe为ECS测井得到地层铁元素的百分含量，Wti为ECS测井得到地层钛元素的百分含量。
需要说明的是：在常规测井过程中，岩石骨架及顺磁性矿物往往会对核磁共振测井产生干扰影响。为了避免上述岩石骨架及顺磁性矿物对核磁共振测井产生的干扰影响，本申请实施例建立T2截止值的计算模型，以便削弱强顺磁性矿物和部分弱顺磁性矿物对T2截止值的干扰，其中，上述T2截止值的计算模型选用强顺磁性矿物和部分弱顺磁性矿物为基础参数，且弱磁性矿物包含钛，非磁性元素为硅，强顺磁性矿物为铁。
具体的，所建立的T2截止值的计算模型如下所示：
T2c＝-14.759*Wsi+11.862*Wfe-944.668*Wti+11.0268。
进一步的，图5为由上述模型计算所得的T2截止值与实验分析所得的T2截止值的对比图，如图5所示，模型计算所得的T2截止值与实验分析所得的T2截止值一致性好，计算精度较高。即采用本申请实施例所提供的核磁可动束缚水饱和度的计算模型，可以计算出来高精准度的T2截止值，进而达到精准判断岩层是否出水。
进一步的，将上述T2截止值的计算模型代入上述可动流体饱和度的计算公式中，得到如下公式：
可选的，在本申请实施例提供的中基性火山岩束缚水出水的判定方法中，所述目标岩心的核磁孔隙度可以通过如下方式确定：步骤B1，依据所述目标岩心的核磁孔隙度和T2截止值测量所得到的结果数据，确定所述目标岩心的核磁共振T2谱分布图；步骤B2，确定所述核磁孔隙度的计算模型；步骤B3，结合所述目标岩心的核磁共振T2谱分布图，通过所述核磁孔隙度的计算模型，确定所述目标岩心的核磁孔隙度，其中，所述核磁孔隙度的计算公式为：φt为核磁孔隙度，T2max为最大横向弛豫时间，ATS为最小横向弛豫时间，S(T)为与横向弛豫时间相关的孔隙度分布函数。
最后，采用本申请实施例所提供的中基性火山岩束缚水出水的判定方法，可得出如图6所示的火山岩地层出水判断效果初始图，其中，图中第一道为深度道，第二道为核磁T2谱，第三道为ECS测井得到地层铁元素的百分含量，第四道为ECS测井得到地层硅元素的百分含量，第五道为ECS测井得到地层钛元素的百分含量，第六道为核磁测井数据计算T2截止值，第七道为核磁测井数据计算T2几何平均值，第八道为核磁测井数据计算核磁孔隙度，第九道为核磁测井数据计算核磁可动束缚水饱和度，第十道为试油道。
其中，在图6所表示的地层数据信息中，深度为3398-3414m的试油段的核磁可动束缚水饱和度平均值为22.8，核磁孔隙度平均值为7.8，其中，该数据组点落至图2中的直线关系式之上，因此，深度为3398-3414m的试油段出水，且试油结论为油水同层。
需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
本申请实施例还提供了一种中基性火山岩束缚水出水的判定装置，需要说明的是，本申请实施例的中基性火山岩束缚水出水的判定装置可以用于执行本申请实施例所提供的用于中基性火山岩束缚水出水的判定方法。以下对本申请实施例提供的中基性火山岩束缚水出水的判定装置进行介绍。
图7是根据本申请实施例的中基性火山岩束缚水出水的判定装置的示意图。如图7所示，该装置包括：获取单元71、确定单元73和判断单元75。
获取单元71，用于获取目标岩心的测量结果，其中，所述目标岩心为待检测中基性火山岩的岩心，所述测量结果为对所述目标岩心进行核磁孔隙度和T2截止值测量所得到的结果数据；
确定单元73，用于依据所述测量结果确定所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度；
判断单元75，用于依据预设关系式对所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定所述待检测中基性火山岩是否出水，其中，所述预设关系式用于表示所述核磁可动束缚水饱和度和所述核磁孔隙度之间的线性关系。
可选地，在本申请实施例提供的中基性火山岩束缚水出水的判定装置中，所述判断单元75包括：建立模块，用于依据所述预设关系式建立火山岩地层出水界限模型，其中，所述火山岩地层出水界限模型为以所述核磁可动束缚水饱和度和核磁孔隙度为坐标轴的二维坐标系模型；第一确定模块，用于将所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度带入至所述火山岩地层出水界限模型中，以确定所述待检测中基性火山岩是否出水；其中，在所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度位于所述火山岩地层出水界限模型中预设关系式第一侧的情况下，则确定所述目标岩心出水；在所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度位于所述火山岩束缚水出水界限模型中预设关系式第二侧的情况下，则确定所述目标岩心未出水。
可选地，在本申请实施例提供的中基性火山岩束缚水出水的判定装置中，所述预设关系式为：φt＝h*Sfswi+k，其中，h、k为预设系数，φt为核磁孔隙度，Sfswi为核磁可动束缚水饱和度。
可选地，在本申请实施例提供的中基性火山岩束缚水出水的判定装置中，所述确定单元73为：第二确定模块，用于依据所述目标岩心的核磁孔隙度和T2截止值测量所得到的结果数据，确定所述目标岩心的核磁共振T2谱分布图；第三确定模块，用于依据所述目标岩心的核磁共振T2谱分布图，确定所述目标岩心的可动流体饱和度和几何平均值；获取模块，用于确定所述核磁可动束缚水饱和度的计算模型，并将所述目标岩心的可动流体饱和度和几何平均值带入所述计算模型中，得到所述目标岩心的核磁可动束缚水饱和度。
可选地，在本申请实施例提供的中基性火山岩束缚水出水的判定装置中，所述核磁可动束缚水饱和度的计算模型为：Sfswi＝a*loge(T2lm)+b*Sff+c，其中，a、b、d为预设系数，Sfswi为核磁可动束缚水饱和度，T2lm为几何平均值，Sff为可动流体饱和度。
可选地，在本申请实施例提供的中基性火山岩束缚水出水的判定装置中，确定模块包括：第一确定子模块，用于确定所述可动流体饱和度的计算公式；第二确定子模块，用于结合所述目标岩心的核磁共振T2谱分布图，通过所述可动流体饱和度的计算公式，确定所述目标岩心的可动流体饱和度，其中，所述可动流体饱和度的计算公式为：Sff为可动流体饱和度，T2max为最大横向弛豫时间，ATS为最小横向弛豫时间，S(T)为与横向弛豫时间相关的孔隙度分布函数，T2c为T2截止值，Wsi为ECS测井得到地层硅元素的百分含量，Wfe为ECS测井得到地层铁元素的百分含量，Wti为ECS测井得到地层钛元素的百分含量，d、e、f、g为预设参数。
可选地，在本申请实施例提供的中基性火山岩束缚水出水的判定装置中，所述确定单元73为：第四确定模块，用于依据所述目标岩心的核磁孔隙度和T2截止值测量所得到的结果数据，确定所述目标岩心的核磁共振T2谱分布图；第五确定模块，用于确定所述核磁孔隙度的计算模型；第六确定模块，用于结合所述目标岩心的核磁共振T2谱分布图，通过所述核磁孔隙度的计算模型，确定所述目标岩心的核磁孔隙度，其中，所述核磁孔隙度的计算公式为：φt为核磁孔隙度，T2max为最大横向弛豫时间，ATS为最小横向弛豫时间；S(T)为与横向弛豫时间相关的孔隙度分布函数。
本申请实施例提供的中基性火山岩束缚水出水的判定装置，通过获取单元71获取目标岩心的测量结果，其中，所述目标岩心为待检测中基性火山岩的岩心，所述测量结果为对所述目标岩心进行核磁孔隙度和T2截止值测量所得到的结果数据；确定单元73依据所述测量结果确定所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度；判断单元75依据预设关系式对所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定所述待检测中基性火山岩是否出水，其中，所述预设关系式用于表示所述核磁可动束缚水饱和度和所述核磁孔隙度之间的线性关系，解决了相关技术中相关技术中无法对火山岩等复杂地层进行精准勘探的问题。进而达到了对岩石类型多样性高、非均质性强、孔隙结构复杂性高的复杂地层(如：火山岩)进行精准勘探的效果。
所述中基性火山岩束缚水出水的判定装置包括处理器和存储器，上述获取单元71、确定单元73、判断单元75等均作为程序单元存储在存储器中，由处理器执行存储在存储器中的上述程序单元来实现相应的功能。
处理器中包含内核，由内核去存储器中调取相应的程序单元。内核可以设置一个或以上，通过调整内核参数来达到对岩石类型多样性高、非均质性强、孔隙结构复杂性高的复杂地层(如：火山岩)进行精准勘探的效果。
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)，存储器包括至少一个存储芯片。
本发明实施例提供了一种存储介质，其上存储有程序，该程序被处理器执行时实现所述中基性火山岩束缚水出水的判定方法。
本发明实施例提供了一种处理器，所述处理器用于运行程序，其中，所述程序运行时执行所述中基性火山岩束缚水出水的判定方法。
本发明实施例提供了一种设备，设备包括处理器、存储器及存储在存储器上并可在处理器上运行的程序，处理器执行程序时实现所述中基性火山岩束缚水出水的判定方法中所包含的步骤。本文中的设备可以是服务器、PC、PAD、手机等。
本申请还提供了一种计算机程序产品，当在数据处理设备上执行时，适于执行初始化有所述中基性火山岩束缚水出水的判定方法中所包含的步骤的程序。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (9)
1.一种中基性火山岩束缚水出水的判定方法，其特征在于，包括：
获取目标岩心的测量结果，其中，所述目标岩心为待检测中基性火山岩的岩心，所述测量结果为对所述目标岩心进行核磁孔隙度和T2截止值测量所得到的结果数据；
依据所述测量结果确定所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度；
依据预设关系式对所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定所述待检测中基性火山岩是否出水，其中，所述预设关系式用于表示所述核磁可动束缚水饱和度和所述核磁孔隙度之间的线性关系；
其中，所述预设关系式为：
φt＝h*Sfswi+k，其中，h、k为预设系数，φt为核磁孔隙度，Sfswi为核磁可动束缚水饱和度。
2.根据权利要求1所述的方法，其特征在于，依据预设关系式对所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定所述待检测中基性火山岩是否出水包括：
依据所述预设关系式建立火山岩地层出水界限模型，其中，所述火山岩地层出水界限模型为以所述核磁可动束缚水饱和度和核磁孔隙度为坐标轴的二维坐标系模型；
将所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度带入至所述火山岩地层出水界限模型中，以确定所述待检测中基性火山岩是否出水；
其中，在所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度位于所述火山岩地层出水界限模型中预设关系式第一侧的情况下，则确定所述目标岩心出水；在所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度位于所述火山岩地层出水界限模型中预设关系式第二侧的情况下，则确定所述目标岩心未出水；
其中，所述预设关系式第一侧为所述火山岩地层出水界限模型中预设关系式的上侧，所述预设关系式第二侧为所述火山岩地层出水界限模型中预设关系式的下侧。
3.根据权利要求1-2中任意一项所述的方法，其特征在于，依据所述测量结果确定所述目标岩心的核磁可动束缚水饱和度为：
依据所述目标岩心的核磁孔隙度和T2截止值测量所得到的结果数据，确定所述目标岩心的核磁共振T2谱分布图；
依据所述目标岩心的核磁共振T2谱分布图，确定所述目标岩心的可动流体饱和度和几何平均值；
确定所述核磁可动束缚水饱和度的计算模型，并将所述目标岩心的可动流体饱和度和几何平均值带入所述计算模型中，得到所述目标岩心的核磁可动束缚水饱和度。
4.根据权利要求3所述的方法，其特征在于，所述核磁可动束缚水饱和度的计算模型为：
Sfswi＝a*loge(T2lm)+b*Sff+c，其中，a、b、d为预设系数，Sfswi为核磁可动束缚水饱和度，T2lm为几何平均值，Sff为可动流体饱和度。
5.根据权利要求3所述的方法，其特征在于，依据所述目标岩心的核磁共振T2谱分布图，确定所述目标岩心的可动流体饱和度包括：
确定所述可动流体饱和度的计算公式；
结合所述目标岩心的核磁共振T2谱分布图，通过所述可动流体饱和度的计算公式，确定所述目标岩心的可动流体饱和度，其中，所述可动流体饱和度的计算公式为：
Sff为可动流体饱和度，T2max为最大横向弛豫时间，ATS为最小横向弛豫时间，S(T)为与横向弛豫时间相关的孔隙度分布函数，T2c为T2截止值，Wsi为ECS测井得到地层硅元素的百分含量，Wfe为ECS测井得到地层铁元素的百分含量，Wti为ECS测井得到地层钛元素的百分含量，d、e、f、g为预设参数。
6.根据权利要求1-2中任意一项所述的方法，其特征在于，依据所述测量结果确定所述目标岩心的核磁孔隙度为：
依据所述目标岩心的核磁孔隙度和T2截止值测量所得到的结果数据，确定所述目标岩心的核磁共振T2谱分布图；
确定所述核磁孔隙度的计算模型；
结合所述目标岩心的核磁共振T2谱分布图，通过所述核磁孔隙度的计算模型，确定所述目标岩心的核磁孔隙度，其中，所述核磁孔隙度的计算公式为：
φt为核磁孔隙度，T2max为最大横向弛豫时间，ATS为最小横向弛豫时间，S(T)为与横向弛豫时间相关的孔隙度分布函数。
7.一种中基性火山岩束缚水出水的判定装置，其特征在于，包括：
获取单元，用于获取目标岩心的测量结果，其中，所述目标岩心为待检测中基性火山岩的岩心，所述测量结果为对所述目标岩心进行核磁孔隙度和T2截止值测量所得到的结果数据；
确定单元，用于依据所述测量结果确定所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度；
判断单元，用于依据预设关系式对所述目标岩心的核磁可动束缚水饱和度和核磁孔隙度进行判断处理，以确定所述待检测中基性火山岩是否出水，其中，所述预设关系式用于表示所述核磁可动束缚水饱和度和所述核磁孔隙度之间的线性关系；
其中，所述预设关系式为：
φt＝h*Sfswi+k，其中，h、k为预设系数，φt为核磁孔隙度，Sfswi为核磁可动束缚水饱和度。
8.一种存储介质，其特征在于，所述存储介质包括存储的程序，其中，所述程序执行权利要求1至5中任意一项所述的中基性火山岩束缚水出水的判定方法。
9.一种处理器，其特征在于，所述处理器用于运行程序，其中，所述程序运行时执行权利要求1至5中任意一项所述的中基性火山岩束缚水出水的判定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910488394.2A CN112049623B (zh) | 2019-06-05 | 2019-06-05 | 一种中基性火山岩束缚水出水的判定方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910488394.2A CN112049623B (zh) | 2019-06-05 | 2019-06-05 | 一种中基性火山岩束缚水出水的判定方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112049623A CN112049623A (zh) | 2020-12-08 |
| CN112049623B true CN112049623B (zh) | 2024-05-28 |
Family
ID=73609493
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910488394.2A Active CN112049623B (zh) | 2019-06-05 | 2019-06-05 | 一种中基性火山岩束缚水出水的判定方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112049623B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113153284B (zh) * | 2021-04-30 | 2023-06-30 | 中国石油天然气股份有限公司 | 束缚水饱和参数的确定方法、装置、设备和存储介质 |
| CN113236366B (zh) * | 2021-06-25 | 2022-06-07 | 中国矿业大学 | 一种水害风险低场核磁共振立体监测方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104990854A (zh) * | 2015-07-06 | 2015-10-21 | 中国石油天然气股份有限公司 | 确定束缚水饱和度的方法及装置 |
| CN105114064A (zh) * | 2015-08-04 | 2015-12-02 | 中国石油天然气股份有限公司 | 确定致密砂岩储层饱和度的方法 |
| CN106443800A (zh) * | 2016-07-18 | 2017-02-22 | 中国石油天然气集团公司 | 一种中基性火山岩的核磁孔隙度校正方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106093350B (zh) * | 2016-08-10 | 2018-05-04 | 中国石油天然气股份有限公司 | 确定非均质碳酸盐岩储层饱和度指数的方法 |
-
2019
- 2019-06-05 CN CN201910488394.2A patent/CN112049623B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104990854A (zh) * | 2015-07-06 | 2015-10-21 | 中国石油天然气股份有限公司 | 确定束缚水饱和度的方法及装置 |
| CN105114064A (zh) * | 2015-08-04 | 2015-12-02 | 中国石油天然气股份有限公司 | 确定致密砂岩储层饱和度的方法 |
| CN106443800A (zh) * | 2016-07-18 | 2017-02-22 | 中国石油天然气集团公司 | 一种中基性火山岩的核磁孔隙度校正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112049623A (zh) | 2020-12-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9244188B2 (en) | System and method for estimating a nuclear magnetic resonance relaxation time cutoff | |
| CN110687612B (zh) | 吸附油和游离油含量连续表征的页岩油分析方法及装置 | |
| EP2859385A1 (en) | Methods of investigating formation samples using nmr data | |
| CN112049623B (zh) | 一种中基性火山岩束缚水出水的判定方法和装置 | |
| CN109630090B (zh) | 储盖组合测井评价方法和装置 | |
| CN107688037A (zh) | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 | |
| CN111208565B (zh) | 基于kt模型的孔缝参数反演方法、装置及存储介质 | |
| Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
| CN110782187B (zh) | 基于含水率的枯竭油气藏型储气库动态评价方法及系统 | |
| CN110470581B (zh) | 确定储层应力敏感程度的方法、装置及存储介质 | |
| CN111208566B (zh) | 基于sca模型的孔缝参数反演方法、装置及存储介质 | |
| Haghighi et al. | A permeability predictive model based on hydraulic flow unit for one of Iranian carbonate tight gas reservoir | |
| CN111241652B (zh) | 一种确定地层原油粘度的方法及装置 | |
| CN111812716A (zh) | 页岩气储层总有机碳含量叠前定量预测方法、装置及设备 | |
| CN111077574A (zh) | 一种确定地层弹性参数的方法、装置及系统 | |
| CN111350497A (zh) | 一种基于含水率的低饱和度油层动态评价方法 | |
| CN113552036A (zh) | 中低孔渗储层油水相对渗透率的确定方法及装置 | |
| Goda et al. | Wettability quantification–prediction of wettability for australian formations | |
| Ghamdi et al. | Emerging Techniques in Measuring Capillary Pressure and Permeability Using NMR and AI | |
| CN111535796A (zh) | 一种测井参数的获取方法和装置 | |
| CN111999331A (zh) | 基于核磁测试的毛管压力确定方法及装置 | |
| CN112987122A (zh) | 一种泥质含量计算方法、装置、电子设备及存储介质 | |
| Wang et al. | Determination of Total organic carbon (TOC) content in shale formations by analyzing well logging measurements | |
| CN111856571B (zh) | 测井频段速度向地震频段速度转换方法、装置及设备 | |
| CN114486620B (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 |