CN113586043A - 束缚水饱和参数的确定方法、装置和计算机设备 - Google Patents
束缚水饱和参数的确定方法、装置和计算机设备 Download PDFInfo
- Publication number
- CN113586043A CN113586043A CN202110812655.9A CN202110812655A CN113586043A CN 113586043 A CN113586043 A CN 113586043A CN 202110812655 A CN202110812655 A CN 202110812655A CN 113586043 A CN113586043 A CN 113586043A
- Authority
- CN
- China
- Prior art keywords
- target reservoir
- bound water
- sampling depth
- volume coefficient
- parameter
- 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.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 430
- 238000000034 method Methods 0.000 title claims abstract description 54
- 239000004927 clay Substances 0.000 claims abstract description 271
- 238000005070 sampling Methods 0.000 claims abstract description 191
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 35
- 239000011148 porous material Substances 0.000 claims description 96
- 230000001419 dependent effect Effects 0.000 claims description 16
- 238000012886 linear function Methods 0.000 claims description 12
- 230000015654 memory Effects 0.000 claims description 12
- 238000011156 evaluation Methods 0.000 abstract description 2
- 238000004590 computer program Methods 0.000 description 15
- 239000011435 rock Substances 0.000 description 14
- 238000012545 processing Methods 0.000 description 11
- 230000002093 peripheral effect Effects 0.000 description 10
- 230000001133 acceleration Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 8
- 238000002474 experimental method Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005553 drilling Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000005251 gamma ray Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 125000005587 carbonate group Chemical group 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000002699 waste material 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
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 (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供了一种束缚水饱和参数的确定方法、装置和计算机设备,属于油田储层评价技术领域。方法包括:确定待研究的目标储层和标准井的目标井段;基于目标井段在每个采样深度的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据,确定目标储层的黏土水体积系数和目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量;对于非标准井的每个采样深度,基于第二元素俘获测井数据和第二常规测井数据、目标储层的黏土水体积系数和目标储层在采样深度的微孔隙束缚水体积系数和黏土含量,确定目标储层在采样深度的束缚水饱和参数。该方法提高了确定的束缚水饱和参数的准确度。
Description
技术领域
本申请涉及油田储层评价技术领域,特别涉及一种束缚水饱和参数的确定方法、装置和计算机设备。
背景技术
在储层开发过程中,储层的束缚水饱和度是评价储层的核心内容,关系到储层参数计算、储量计算、油藏描述和开发方案制定。因此,确定储层的束缚水饱和度,对于高效的开发储层十分重要。
相关技术中,确定束缚水饱和度的过程一般为:对待研究储层的标准井进行钻井取芯工作,获得储层的多个代表性的岩样,对岩样进行核磁共振实验,确定岩样的束缚水饱和度;基于确定的岩样的束缚水饱和度,建立与未进行核磁共振实验的非标准井的束缚水饱和度的对应关系,基于该对应关系确定目标储层在非标准井处的束缚水饱和度。由于该方法中的多个代表性的岩样需要事先进行钻井取芯工作来获取,而钻井取芯工作费时费力,从而降低了确定束缚水饱和度的效率。
发明内容
本申请实施例提供了一种束缚水饱和参数的确定方法、装置和计算机设备,能够提高确定束缚水饱和参数的效率。所述技术方案如下:
一方面,提供了一种束缚水饱和参数的确定方法,所述方法包括:
确定待研究的目标储层和标准井的目标井段,所述标准井为所述目标储层上的井,所述目标井段为与所述目标储层的深度对应的井段;
获取所述目标井段在每个采样深度的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据;
基于所述目标井段在每个采样深度的核磁共振测井数据、所述第一元素俘获测井数据和所述第一常规测井数据,确定所述目标储层的黏土水体积系数和所述目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量;
获取所述目标储层上的非标准井在每个采样深度的第二元素俘获测井数据和第二常规测井数据;
对于所述非标准井的每个采样深度,基于所述第二元素俘获测井数据和所述第二常规测井数据、所述目标储层的黏土水体积系数和所述目标储层在所述采样深度的微孔隙束缚水体积系数和黏土含量,确定所述目标储层在所述采样深度的束缚水饱和参数。
在一种可能的实现方式中,所述核磁共振测井数据包括总束缚水饱和参数、泥质束缚水饱和参数、泥质束缚水孔隙参数和第一总孔隙参数,所述第一元素俘获测井数据包括第一干黏土含量,所述第一常规测井数据包括泥质含量;
所述基于所述目标井段在每个采样深度的核磁共振测井数据、所述第一元素俘获测井数据和所述第一常规测井数据,确定所述目标储层的黏土水体积系数和所述目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量,包括:
对于所述目标储层的每个采样深度,确定所述干黏土含量与所述泥质束缚水孔隙参数之和,得到所述目标储层在所述采样深度的黏土含量;
基于所述目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数,确定所述目标储层的黏土水体积系数;
基于所述目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定所述目标储层的最小束缚水体积系数;
对于所述目标储层的每个采样深度,基于所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,确定所述目标储层在所述采样深度的分散粉砂含量;
基于所述目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数,确定所述目标储层的微孔隙束缚水参数;
对于所述目标储层的每个采样深度,基于所述微孔隙束缚水参数、所述黏土水体积系数、所述最小束缚水体积系数,以及所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,确定所述目标储层在所述采样深度的微孔隙束缚水体积系数。
在一种可能的实现方式中,,所述基于所述目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数,确定所述目标储层的黏土水体积系数,包括:
对所述目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数进行拟合,得到所述目标储层的黏土水体积系数。
在一种可能的实现方式中,,所述基于所述目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定所述目标储层的最小束缚水体积系数,包括:
基于所述目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定第一交会图,所述第一交会图的横坐标为黏土含量,纵坐标为总束缚水饱和参数;
确定所述第一交会图上的第一目标点和第二目标点,所述第一目标点的总束缚水饱和参数为最小值,所述第二目标点的总束缚水饱和参数和黏土含量均为最大值;
基于所述第一目标点和所述第二目标点,组成线性函数;
确定所述线性函数与所述纵坐标的交点处的值为所述最小束缚水体积系数。
在一种可能的实现方式中,,所述对于所述目标储层的每个采样深度,基于所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,确定所述目标储层在所述采样深度的分散粉砂含量,包括:
基于所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,通过第一关系数据,得到所述目标储层在所述采样深度的分散粉砂含量,所述第一关系数据的自变量为所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,因变量为所述目标储层在所述采样深度的分散粉砂含量。
在一种可能的实现方式中,,所述基于所述目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数,确定所述目标储层的微孔隙束缚水参数,包括:
对所述目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数进行拟合,得到所述微孔隙束缚水参数。
在一种可能的实现方式中,,所述对于所述目标储层的每个采样深度,基于所述微孔隙束缚水参数、所述黏土水体积系数、所述最小束缚水体积系数,以及所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,确定所述目标储层在所述采样深度的微孔隙束缚水体积系数,包括:
基于所述微孔隙束缚水参数、所述黏土水体积系数、以及所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,通过第二关系数据,得到所述目标储层在所述采样深度的第一微孔隙束缚水体积系数,所述第二关系数据的参数为所述微孔隙束缚水参数,自变量为所述黏土水体积系数和所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,因变量为所述目标储层在所述采样深度的第一微孔隙束缚水体积系数;
若所述第一微孔隙束缚水体积系数在预设系数范围内,则所述第一微孔隙束缚水体积系数为所述微孔隙束缚水体积系数,所述预设系数范围的最小值为所述最小束缚水体积系数;
若所述第一微孔隙束缚水体积系数不在预设系数范围内,则基于所述预设系数范围,对所述第一微孔隙束缚水体积系数进行校正,得到所述微孔隙束缚水体积系数。
在一种可能的实现方式中,所述第二元素俘获测井数据包括第二干黏土含量,所述第二常规测井数据包括第二总孔隙参数;
所述对于所述非标准井的每个采样深度,基于所述第二元素俘获测井数据和所述第二常规测井数据、所述目标储层的黏土水体积系数和所述目标储层在所述采样深度的微孔隙束缚水体积系数和黏土含量,确定所述目标储层在所述采样深度的束缚水饱和参数,包括:
基于所述黏土水体积系数和所述微孔隙束缚水体积系数,生成包含所述黏土水体积系数和所述微孔隙束缚水体积系数、且以所述第二干黏土含量、所述第二总孔隙参数和所述黏土含量为自变量、所述束缚水饱和参数为因变量的第三关系数据;
基于所述干黏土含量、所述第二总孔隙参数和所述黏土含量,通过所述第三关系数据,得到所述束缚水饱和参数。
另一方面,提供了一种束缚水饱和参数的确定装置,所述装置包括:
第一确定模块,用于确定待研究的目标储层和标准井的目标井段,所述标准井为所述目标储层上的井,所述目标井段为与所述目标储层的深度对应的井段;
第一获取模块,用于获取所述目标井段在每个采样深度的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据;
第二确定模块,用于基于所述目标井段在每个采样深度的核磁共振测井数据、所述第一元素俘获测井数据和所述第一常规测井数据,确定所述目标储层的黏土水体积系数和所述目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量;
第二获取模块,用于获取所述目标储层上的非标准井在每个采样深度的第二元素俘获测井数据和第二常规测井数据;
第三确定模块,用于对于所述非标准井的每个采样深度,基于所述第二元素俘获测井数据和所述第二常规测井数据、所述目标储层的黏土水体积系数和所述目标储层在所述采样深度的微孔隙束缚水体积系数和黏土含量,确定所述目标储层在所述采样深度的束缚水饱和参数。
在一种可能的实现方式中,所述核磁共振测井数据包括总束缚水饱和参数、泥质束缚水饱和参数、泥质束缚水孔隙参数和第一总孔隙参数,所述第一元素俘获测井数据包括第一干黏土含量,所述第一常规测井数据包括泥质含量;
所述第二确定模块,包括:
第一确定单元,用于对于所述目标储层的每个采样深度,确定所述干黏土含量与所述泥质束缚水孔隙参数之和,得到所述目标储层在所述采样深度的黏土含量;
第二确定单元,用于基于所述目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数,确定所述目标储层的黏土水体积系数;
第三确定单元,用于基于所述目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定所述目标储层的最小束缚水体积系数;
第四确定单元,用于对于所述目标储层的每个采样深度,基于所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,确定所述目标储层在所述采样深度的分散粉砂含量;
第五确定单元,用于基于所述目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数,确定所述目标储层的微孔隙束缚水参数;
第六确定单元,用于对于所述目标储层的每个采样深度,基于所述微孔隙束缚水参数、所述黏土水体积系数、所述最小束缚水体积系数,以及所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,确定所述目标储层在所述采样深度的微孔隙束缚水体积系数。
在一种可能的实现方式中,所述第二确定单元,用于:
对所述目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数进行拟合,得到所述目标储层的黏土水体积系数。
在一种可能的实现方式中,所述第三确定单元,用于:
基于所述目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定第一交会图,所述第一交会图的横坐标为黏土含量,纵坐标为总束缚水饱和参数;
确定所述第一交会图上的第一目标点和第二目标点,所述第一目标点的总束缚水饱和参数为最小值,所述第二目标点的总束缚水饱和参数和黏土含量均为最大值;
基于所述第一目标点和所述第二目标点,组成线性函数;
确定所述线性函数与所述纵坐标的交点处的值为所述最小束缚水体积系数。
在一种可能的实现方式中,所述第四确定单元,用于:
基于所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,通过第一关系数据,得到所述目标储层在所述采样深度的分散粉砂含量,所述第一关系数据的自变量为所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,因变量为所述目标储层在所述采样深度的分散粉砂含量。
在一种可能的实现方式中,所述第五确定单元,用于:
对所述目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数进行拟合,得到所述微孔隙束缚水参数。
在一种可能的实现方式中,所述第六确定单元,用于:
基于所述微孔隙束缚水参数、所述黏土水体积系数、以及所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,通过第二关系数据,得到所述目标储层在所述采样深度的第一微孔隙束缚水体积系数,所述第二关系数据的参数为所述微孔隙束缚水参数,自变量为所述黏土水体积系数和所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,因变量为所述目标储层在所述采样深度的第一微孔隙束缚水体积系数;
若所述第一微孔隙束缚水体积系数在预设系数范围内,则所述第一微孔隙束缚水体积系数为所述微孔隙束缚水体积系数,所述预设系数范围的最小值为所述最小束缚水体积系数;
若所述第一微孔隙束缚水体积系数不在预设系数范围内,则基于所述预设系数范围,对所述第一微孔隙束缚水体积系数进行校正,得到所述微孔隙束缚水体积系数。
在一种可能的实现方式中,所述第二元素俘获测井数据包括第二干黏土含量,所述第二常规测井数据包括第二总孔隙参数;所述第三确定模块,用于:
基于所述黏土水体积系数和所述微孔隙束缚水体积系数,生成包含所述黏土水体积系数和所述微孔隙束缚水体积系数、且以所述第二干黏土含量、所述第二总孔隙参数和所述黏土含量为自变量、所述束缚水饱和参数为因变量的第三关系数据;基于所述干黏土含量、所述第二总孔隙参数和所述黏土含量,通过所述第三关系数据,得到所述束缚水饱和参数。
另一方面,提供了一种计算机设备,所述计算机设备包括一个或多个处理器和一个或多个存储器,所述一个或多个存储器中存储有至少一条指令,所述至少一条指令由所述一个或多个处理器加载并执行以实现上述任一实现方式所述的束缚水饱和参数的确定方法所执行的操作。
另一方面,提供了一种计算机可读存储介质,所述计算机可读存储介质中存储至少一条指令,所述至少一条指令由处理器加载并执行,以实现上述任一实现方式所述的束缚水饱和参数的确定方法所执行的操作。
另一方面,提供了一种计算机程序产品或计算机程序,所述计算机程序产品或所述计算机程序包括计算机程序代码,所述计算机程序代码存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取所述计算机程序代码,处理器执行所述计算机程序代码,使得所述计算机设备执行上述的束缚水饱和参数的确定方法所执行的操作。
本申请实施例提供的技术方案的有益效果至少包括:
本申请实施例提供了一种束缚水饱和参数的确定方法,该方法基于目标储层上标准井的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据来确定目标储层的黏土水体积系数和目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量,并基于目标储层的黏土水体积系数和目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量来确定具有第二元素俘获测井数据和第二常规测井数据的非标准井的目标储层的束缚水饱和参数,这样,避免了确定束缚水饱和参数时进行钻井取芯工作,省事省力,从而提高了确定束缚水饱和参数的效率。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种束缚水饱和参数的确定方法的流程图;
图2是本申请实施例提供的一种黏土水体积系数拟合图;
图3是本申请实施例提供的一种第一交会图的示意图;
图4是本申请实施例提供的一种微孔隙束缚水参数的拟合图;
图5是本申请实施例提供的一种束缚水饱和参数的确定结果对比图;
图6是本申请实施例提供的一种束缚水饱和参数的确定结果对比图;
图7是本申请实施例提供的一种束缚水饱和参数的确定装置的框图;
图8是本申请实施例提供的一种计算机设备的框图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,下面将结合附图对本申请实施方式作进一步地详细描述。
本申请的说明书和权利要求书及所述附图中的术语“第一”、“第二”、“第三”和“第四”等是用于区别不同对象,而不是用于描述特定顺序。此外,术语“包括”和“具有”以及它们的任意变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、***、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
本申请实施例提供一种束缚水饱和参数的确定方法,参见图1,方法包括:
步骤101:计算机设备确定待研究的目标储层和标准井的目标井段。
其中,标准井为目标储层上的井,目标井段为与目标储层的深度对应的井段。待研究的目标储层可以为低孔隙度低渗透率的砂岩储层。标准井为目标储层上进行了核磁共振测井、元素俘获测井和常规测井的井,主要是在标准井的目标井段上进行核磁共振测井、元素俘获测井和常规测井,目标井段的测井资料丰富。
需要说明的是,目标井段的核磁共振测井为适应该地区模式下的核磁共振测井,其测井模式是能反映储层流体特性的测量模式。常规测井为适应该地区模式下的常规测井,包含自然伽马测井、电阻率测井、补偿中子测井、岩性密度测井、数字声波测井等。元素俘获测井为适应该地区模式下的元素俘获测井,获取的元素俘获测井资料能够反映岩石矿物含量。常规测井、地层元素测井、核磁共振测井按照《石油电缆测井作业技术规范SY/T5600-2010》标准规定的流程进行。
步骤102:计算机设备获取目标井段在每个采样深度的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据。
其中,核磁共振测井数据包括总束缚水饱和参数、泥质束缚水饱和参数、泥质束缚水孔隙参数和第一总孔隙参数;总束缚水饱和参数为总束缚水饱和度,泥质束缚水饱和参数为泥质束缚水饱和度,泥质束缚水孔隙参数为泥质束缚水孔隙度,第一总孔隙参数为第一总孔隙度。第一元素俘获测井数据包括第一干黏土含量,第一常规测井数据包括泥质含量。
计算机设备分别对核磁共振测井资料、元素俘获测井资料和常规测井资料进行处理,得到核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据。计算机设备对第一元素俘获测井资料和第一常规测井资料的处理按照《裸眼井单井测井数据处理流程SY/T5360-2004》、《探井测井数据处理与资料解释技术规范SY/T 6451-2010》和《储层参数的测井计算方法SY/T 5940-2010》标准规定的流程进行,对核磁共振测井资料的处理按照《核磁共振测井资料处理及解释规范SY/T 6617-2016》标准规定的流程进行。
步骤103:计算机设备基于目标井段在每个采样深度的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据,确定目标储层的黏土水体积系数和目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量。
该步骤包括以下步骤(1)-(6)实现:
(1)计算机设备对于目标储层的每个采样深度,确定干黏土含量与泥质束缚水孔隙参数之和,得到目标储层在该采样深度的黏土含量。
计算机设备通过以下公式一,确定干黏土含量与泥质束缚水孔隙参数之和,得到目标储层的黏土含量。
公式一:Vcl=Vcld+φcwb
其中,Vcld为干黏土含量,单位为v/v;Vcl为黏土含量,单位为v/v;φcwb为泥质束缚水孔隙度,单位为v/v。
(2)计算机设备基于目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数,确定目标储层的黏土水体积系数。
在该步骤中,计算机设备对目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数进行拟合,得到目标储层的黏土水体积系数。计算机设备通过公式二对黏土含量和泥质束缚水孔隙参数进行拟合。
公式二:φcwb=tVcl
其中,Vcl为黏土含量,单位为v/v;φcwb为泥质束缚水孔隙度,单位为v/v;t为黏土水体积系数。
参见图2,计算机设备对在每个采样深度的黏土含量和泥质束缚水孔隙参数采用二乘法进行拟合,得到黏土水体积系数为0.322486,其相关性为0.898217,说明黏土含量和泥质束缚水孔隙参数的相关性较好,得到的黏土水体积系数准确度较高。
(3)计算机设备基于目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定目标储层的最小束缚水体积系数。
该步骤可以通过以下步骤A1-A4实现:
A1:计算机设备基于目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定第一交会图,第一交会图的横坐标为黏土含量,纵坐标为总束缚水饱和参数。
参见图3,图3为某个第一交会图的示意图,该图的横坐标为黏土含量,纵坐标为总束缚水饱和参数,将多个黏土含量和多个总束缚水饱和参数的值分别绘制到该图上,得到包括多个黏土含量和多个总束缚水饱和参数的第一交会图。
A2:计算机设备确定第一交会图上的第一目标点和第二目标点,第一目标点的总束缚水饱和参数为最小值,第二目标点的总束缚水饱和参数和黏土含量均为最大值。
其中,第一目标点为总束缚水饱和参数最小的点,即该点为最小总束缚水饱和度位置;第二目标点为黏土含量为1、总束缚水饱和度为1的泥岩点位置。继续参见图3,总束缚水饱和度的最小值为0.2,因此图中的第一目标点为总束缚水饱和度为0.2的点。总束缚水饱和度的最大值为1,黏土含量的最大值为1,即图中的第二目标点为坐标为(1,1)的点。
A3:计算机设备基于第一目标点和第二目标点,组成线性函数。
在该步骤中,计算机设备连接第一目标点和第二目标点,组成线性函数,如图3所示。
A4:计算机设备确定线性函数与纵坐标的交点处的值为最小束缚水体积系数。
继续参见图3,确定线性函数与纵坐标的交点,将该点的总束缚水饱和参数的值作为最小束缚水体积系数,如图3中的最小束缚水体积系数gmin。
(4)计算机设备对于目标储层的每个采样深度,基于黏土水体积系数、目标储层在采样深度的泥质含量和干黏土含量,确定目标储层在采样深度的分散粉砂含量。
需要说明的是,分散粉砂含量为微粒级分散粉砂的含量。在该步骤中,计算机设备基于黏土水体积系数、目标储层在采样深度的泥质含量和干黏土含量,通过第一关系数据,得到目标储层在该采样深度的分散粉砂含量。第一关系数据的自变量为黏土水体积系数、目标储层在该采样深度的泥质含量和干黏土含量,因变量为目标储层在该采样深度的分散粉砂含量。
其中,Vps为分散粉砂含量,单位为v/v;Vsh泥质含量,单位为v/v;Vcld为干黏土含量,单位为v/v,t为黏土水体积系数。
(5)计算机设备基于目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数,确定目标储层的微孔隙束缚水参数。
在该步骤中,计算机设备对目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数进行拟合,得到微孔隙束缚水参数。
其中,微孔隙束缚水参数包括第一参数和第二参数。计算机设备对总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数采用最小二乘法、通过以下公式三进行拟合,得到第一参数和第二参数。
公式三:
其中,Swi为总束缚水饱和度,单位为v/v;Swc为泥质束缚水饱和度,单位为v/v;Vcl为黏土含量,单位为v/v;Vcld为干黏土含量,单位为v/v;Vps为分散粉砂含量,单位为v/v;φ为第一总孔隙度;A为第一参数,B为第二参数。
以图4为例,计算机设备对总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数进行拟合,将作为自变量,即图中的横坐标,将作为因变量,即图中的纵坐标,通过二乘法拟合,得到第一参数A为0.494125,第二参数B为0.395571,拟合的相关系数为0.871886,说明总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数的相关性较高,即得到的第一参数和第二参数的准确度高。
(6)计算机设备对于目标储层的每个采样深度,基于微孔隙束缚水参数、黏土水体积系数、以及目标储层在采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,确定目标储层在该采样深度的微孔隙束缚水体积系数。
该步骤可以通过以下步骤(1)-(2)实现:
(1)计算机设备基于微孔隙束缚水参数、黏土水体积系数、以及目标储层在采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,通过第二关系数据,得到目标储层在该采样深度的第一微孔隙束缚水体积系数,第二关系数据的参数为微孔隙束缚水参数,自变量为黏土水体积系数和目标储层在采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,因变量为目标储层在采样深度的第一微孔隙束缚水体积系数。
其中,g为第一微孔隙束缚水体积系数;Vcld为干黏土含量,单位为v/v;Vps为分散粉砂含量,单位为v/v;φ为第一总孔隙度;A为第一参数,B为第二参数;gmin为最小束缚水体积系数,[gmin,1]为预设系数范围。
(2)计算机设备若第一微孔隙束缚水体积系数在预设系数范围内,则第一微孔隙束缚水体积系数为微孔隙束缚水体积系数。若第一微孔隙束缚水体积系数不在预设系数范围内,则基于预设系数范围,对第一微孔隙束缚水体积系数进行校正,得到微孔隙束缚水体积系数。
其中,预设系数范围的最小值为最小束缚水体积系数,最大值为常数1。校正的方法可以为将最小束缚水体积系数作为微孔隙束缚水体积系数。
步骤104:计算机设备获取目标储层上的非标准井在每个采样深度的第二元素俘获测井数据和第二常规测井数据。
其中,第二元素俘获测井数据包括第二干黏土含量,第二常规测井数据包括第二总孔隙参数,第二总孔隙参数为第二总孔隙度。计算机设备从非标准井的元素俘获测井资料中获取第二元素俘测井数据,以及从非标准井的常规测井资料中获取第二常规测井数据。
需要说明的是,非标准井的每个深度与目标井段的每个深度相同。非标准井的常规测井为适应该地区模式下的常规测井,包含自然伽马测井、电阻率测井、补偿中子测井、岩性密度测井、数字声波测井等。元素俘获测井为适应该地区模式下的元素俘获测井,获取的元素俘获测井资料能够反映岩石矿物含量。其测井资料的处理技术与步骤101-102相同,在此不再赘述。
步骤105:计算机设备对于非标准井的每个采样深度,基于第二元素俘获测井数据和第二常规测井数据、目标储层的黏土水体积系数和目标储层在该采样深度的微孔隙束缚水体积系数和黏土含量,确定目标储层在该采样深度的束缚水饱和参数。
其中,束缚水饱和参数为束缚水饱和度,包括目标储层在非标准井处和标准井处的束缚水饱和度。标准井处的束缚水饱和度可以通过核磁共振测井资料直接得到,也可以通过本申请提供的方法对核磁共振测井得到的束缚水饱和度进行验证。
该步骤可以通过以下步骤(1)-(2)实现:
(1)计算机设备基于黏土水体积系数和微孔隙束缚水体积系数,生成包含黏土水体积系数和微孔隙束缚水体积系数、且以第二干黏土含量和第二总孔隙参数为自变量、束缚水饱和参数为因变量的第三关系数据。
其中,Swi为束缚水饱和度,单位为v/v;g为微孔隙束缚水体积系数;Vcld为第二干黏土含量,单位为v/v;φ为第二总孔隙度;t为黏土水体积系数。Vcl为第二黏土含量,单位为v/v;
(2)计算机设备基于干黏土含量、第二总孔隙参数和黏土含量,通过第三关系数据,得到束缚水饱和参数。
其中,计算机设备将干黏土含量、第二总孔隙参数和黏土含量代入第三关系数据,得到束缚水饱和度。
在本申请实施例中,由于标准井和非标准井均为目标储层上的井,标准井处和非标准井处的储层具有相同或相近的沉积、成岩和构造等地质条件,因此,可以基于标准井的核磁共振测井技术测得的束缚水饱和度、第一元素俘获测井数据和第一常规测井数据,建立与非标准井的第二元素俘获测井数据和第二常规测井数据之间的对应关系,该对应关系为目标储层的泥质分布形态与束缚水饱和度之间的映射关系,符合对目标储层的岩石物理认识,进而基于该对应关系即可确定非标准井处的束缚水饱和度,该方法准确、可靠,从而避免了对非标准井进行核磁共振测井带来的人力、物力和财力的浪费,且避免了对标准井进行钻井取芯的工作,提高了确定束缚水饱和度的效率。
参见图5,图5为通过本申请实施例提供的方法确定的整个目标井段的束缚水饱和参数的结果对比图。图5中第一道为岩性曲线道,其中GR为自然伽玛曲线,单位为API,SP为自然电位曲线,单位为mv。第二道为深度道,深度单位为m。第三道为电阻率道,其中RLLD为深侧向电阻率曲线,单位为Ω.m,RLLS为浅侧向电阻率曲线,单位为Ω.m。第四道为孔隙度道,其中AC为声波曲线,单位为μs/ft,CNL为补偿中子曲线,单位为%,DEN为补偿密度曲线,单位为g/cm3。第五道为标准T2(弛豫时间)谱道,即T2谱曲线。第六道为地层元素道,其中WCLA为黏土重量百分比曲线,单位w/w,WQFM为石英骨架重量百分比曲线,单位w/w,WCAR为碳酸盐岩矿物重量百分比曲线,单位w/w。第七道为地层剖面道,其中VCL为根据地层元素测井获取的黏土含量曲线,单位v/v,VSI为分散粉砂含量曲线,单位:v/v,VSAND为岩石骨架体积含量曲线,单位v/v;PHIE为总孔隙度曲线,单位v/v;第八道为束缚水饱和度比较道,其中new method_SWI为通过本申请实施例提供的方法确定的束缚水饱和度曲线,单位:v/v,CMR_SWI为岩样核磁共振实验得到的束缚水饱和度曲线,通过对比可以看出,通过本申请实施例的方法确定的束缚水饱和度与岩样核磁共振实验确定的束缚水饱和度吻合度高,说明通过本申请实施例的方法确定的束缚水饱和度的准确度高。
参见图6,图6为将通过本申请实施例的方法确定的束缚水饱和度与岩样核磁共振实验确定的束缚水饱和度进行对比的拟合图;横坐标为岩样核磁共振实验确定的束缚水饱和度,纵坐标为同一深度本申请实施例的方法确定的束缚水饱和度,拟合的曲线方程的相关系数R=0.81,说明本申请实施例的方法确定的束缚水饱和度的平均误差小,与岩样核磁共振实验确定的束缚水饱和度高度吻合,进一步说明了本申请实施例的方法确定的束缚水饱和度的准确度高,具有良好的适应性。
本申请实施例提供了一种束缚水饱和参数的确定方法,该方法基于目标储层上标准井的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据来确定目标储层的黏土水体积系数和目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量,并基于目标储层的黏土水体积系数和目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量来确定具有第二元素俘获测井数据和第二常规测井数据的非标准井的目标储层的束缚水饱和参数,这样,避免了确定束缚水饱和参数时进行钻井取芯工作,省事省力,从而提高了确定束缚水饱和参数的效率。
本申请实施例还提供了一种束缚水饱和参数的确定装置,参见图7,装置包括:
第一确定模块701,用于确定待研究的目标储层和标准井的目标井段,标准井为目标储层上的井,目标井段为与目标储层的深度对应的井段;
第一获取模块702,用于获取目标井段在每个采样深度的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据;
第二确定模块703,用于基于目标井段在每个采样深度的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据,确定目标储层的黏土水体积系数和目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量;
第二获取模块704,用于获取目标储层上的非标准井在每个采样深度的第二元素俘获测井数据和第二常规测井数据;
第三确定模块705,用于对于非标准井的每个采样深度,基于第二元素俘获测井数据和第二常规测井数据、目标储层的黏土水体积系数和目标储层在采样深度的微孔隙束缚水体积系数和黏土含量,确定目标储层在采样深度的束缚水饱和参数。
在一种可能的实现方式中,核磁共振测井数据包括总束缚水饱和参数、泥质束缚水饱和参数、泥质束缚水孔隙参数和第一总孔隙参数,第一元素俘获测井数据包括第一干黏土含量,第一常规测井数据包括泥质含量;
第二确定模块703,包括:
第一确定单元,用于对于目标储层的每个采样深度,确定干黏土含量与泥质束缚水孔隙参数之和,得到目标储层在采样深度的黏土含量;
第二确定单元,用于基于目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数,确定目标储层的黏土水体积系数;
第三确定单元,用于基于目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定目标储层的最小束缚水体积系数;
第四确定单元,用于对于目标储层的每个采样深度,基于黏土水体积系数、目标储层在采样深度的泥质含量和干黏土含量,确定目标储层在采样深度的分散粉砂含量;
第五确定单元,用于基于目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数,确定目标储层的微孔隙束缚水参数;
第六确定单元,用于对于目标储层的每个采样深度,基于微孔隙束缚水参数、黏土水体积系数、最小束缚水体积系数,以及目标储层在采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,确定目标储层在采样深度的微孔隙束缚水体积系数。
在一种可能的实现方式中,第二确定单元,用于:
对目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数进行拟合,得到目标储层的黏土水体积系数。
在一种可能的实现方式中,第三确定单元,用于:
基于目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定第一交会图,第一交会图的横坐标为黏土含量,纵坐标为总束缚水饱和参数;
确定第一交会图上的第一目标点和第二目标点,第一目标点的总束缚水饱和参数为最小值,第二目标点的总束缚水饱和参数和黏土含量均为最大值;
基于第一目标点和第二目标点,组成线性函数;
确定线性函数与纵坐标的交点处的值为最小束缚水体积系数。
在一种可能的实现方式中,第四确定单元,用于:
基于黏土水体积系数、目标储层在采样深度的泥质含量和干黏土含量,通过第一关系数据,得到目标储层在采样深度的分散粉砂含量,第一关系数据的自变量为黏土水体积系数、目标储层在采样深度的泥质含量和干黏土含量,因变量为目标储层在采样深度的分散粉砂含量。
在一种可能的实现方式中,第五确定单元,用于:
对目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数进行拟合,得到微孔隙束缚水参数。
在一种可能的实现方式中,第六确定单元,用于:
基于微孔隙束缚水参数、黏土水体积系数、以及目标储层在采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,通过第二关系数据,得到目标储层在采样深度的第一微孔隙束缚水体积系数,第二关系数据的参数为微孔隙束缚水参数,自变量为黏土水体积系数和目标储层在采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,因变量为目标储层在采样深度的第一微孔隙束缚水体积系数;
若第一微孔隙束缚水体积系数在预设系数范围内,则第一微孔隙束缚水体积系数为微孔隙束缚水体积系数,预设系数范围的最小值为最小束缚水体积系数;
若第一微孔隙束缚水体积系数不在预设系数范围内,则基于预设系数范围,对第一微孔隙束缚水体积系数进行校正,得到微孔隙束缚水体积系数。
在一种可能的实现方式中,第二元素俘获测井数据包括第二干黏土含量,第二常规测井数据包括第二总孔隙参数;第三确定模块705,用于:
基于黏土水体积系数和微孔隙束缚水体积系数,生成包含黏土水体积系数和微孔隙束缚水体积系数、且以第二干黏土含量、第二总孔隙参数和黏土含量为自变量、束缚水饱和参数为因变量的第三关系数据;基于干黏土含量、第二总孔隙参数和黏土含量,通过第三关系数据,得到束缚水饱和参数。
图8示出了本申请一个示例性实施例提供的计算机设备800的结构框图。该计算机设备800可以是便携式移动计算机设备,比如:智能手机、平板电脑、MP3播放器(MovingPicture Experts Group Audio Layer III,动态影像专家压缩标准音频层面3)、MP4(Moving Picture Experts Group Audio Layer IV,动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。计算机设备800还可能被称为用户设备、便携式计算机设备、膝上型计算机设备、台式计算机设备等其他名称。
通常,计算机设备800包括有:处理器801和存储器802。
处理器801可以包括一个或多个处理核心,比如4核心处理器、8核心处理器等。处理器801可以采用DSP(Digital Signal Processing,数字信号处理)、FPGA(Field-Programmable Gate Array,现场可编程门阵列)、PLA(Programmable Logic Array,可编程逻辑阵列)中的至少一种硬件形式来实现。处理器801也可以包括主处理器和协处理器,主处理器是用于对在唤醒状态下的数据进行处理的处理器,也称CPU(Central ProcessingUnit,中央处理器);协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中,处理器801可以集成有GPU(Graphics Processing Unit,图像处理器),GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中,处理器801还可以包括AI(Artificial Intelligence,人工智能)处理器,该AI处理器用于处理有关机器学习的计算操作。
存储器802可以包括一个或多个计算机可读存储介质,该计算机可读存储介质可以是非暂态的。存储器802还可包括高速随机存取存储器,以及非易失性存储器,比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中,存储器802中的非暂态的计算机可读存储介质用于存储至少一个指令,该至少一个指令用于被处理器801所执行以实现本申请中方法实施例提供的束缚水饱和参数的确定方法。
在一些实施例中,计算机设备800还可选包括有:***设备接口803和至少一个***设备。处理器801、存储器802和***设备接口803之间可以通过总线或信号线相连。各个***设备可以通过总线、信号线或电路板与***设备接口803相连。具体地,***设备包括:射频电路804、显示屏805、摄像头组件806、音频电路807、定位组件808和电源809中的至少一种。
***设备接口803可被用于将I/O(Input/Output,输入/输出)相关的至少一个***设备连接到处理器801和存储器802。在一些实施例中,处理器801、存储器802和***设备接口803被集成在同一芯片或电路板上;在一些其他实施例中,处理器801、存储器802和***设备接口803中的任意一个或两个可以在单独的芯片或电路板上实现,本实施例对此不加以限定。
射频电路804用于接收和发射RF(Radio Frequency,射频)信号,也称电磁信号。射频电路804通过电磁信号与通信网络以及其他通信设备进行通信。射频电路804将电信号转换为电磁信号进行发送,或者,将接收到的电磁信号转换为电信号。可选地,射频电路804包括:天线***、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路804可以通过至少一种无线通信协议来与其它计算机设备进行通信。该无线通信协议包括但不限于:万维网、城域网、内联网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity,无线保真)网络。在一些实施例中,射频电路804还可以包括NFC(Near Field Communication,近距离无线通信)有关的电路,本申请对此不加以限定。
显示屏805用于显示UI(User Interface,用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏805是触摸显示屏时,显示屏805还具有采集在显示屏805的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器801进行处理。此时,显示屏805还可以用于提供虚拟按钮和/或虚拟键盘,也称软按钮和/或软键盘。在一些实施例中,显示屏805可以为一个,设置在计算机设备800的前面板;在另一些实施例中,显示屏805可以为至少两个,分别设置在计算机设备800的不同表面或呈折叠设计;在另一些实施例中,显示屏805可以是柔性显示屏,设置在计算机设备800的弯曲表面上或折叠面上。甚至,显示屏805还可以设置成非矩形的不规则图形,也即异形屏。显示屏805可以采用LCD(Liquid Crystal Display,液晶显示屏)、OLED(Organic Light-EmittingDiode,有机发光二极管)等材质制备。
摄像头组件806用于采集图像或视频。可选地,摄像头组件806包括前置摄像头和后置摄像头。通常,前置摄像头设置在计算机设备的前面板,后置摄像头设置在计算机设备的背面。在一些实施例中,后置摄像头为至少两个,分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种,以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality,虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中,摄像头组件806还可以包括闪光灯。闪光灯可以是单色温闪光灯,也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合,可以用于不同色温下的光线补偿。
音频电路807可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波,并将声波转换为电信号输入至处理器801进行处理,或者输入至射频电路804以实现语音通信。出于立体声采集或降噪的目的,麦克风可以为多个,分别设置在计算机设备800的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器801或射频电路804的电信号转换为声波。扬声器可以是传统的薄膜扬声器,也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时,不仅可以将电信号转换为人类可听见的声波,也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中,音频电路807还可以包括耳机插孔。
定位组件808用于定位计算机设备800的当前地理位置,以实现导航或LBS(Location Based Service,基于位置的服务)。定位组件808可以是基于美国的GPS(GlobalPositioning System,全球定位***)、中国的北斗***或俄罗斯的伽利略***的定位组件。
电源809用于为计算机设备800中的各个组件进行供电。电源809可以是交流电、直流电、一次性电池或可充电电池。当电源809包括可充电电池时,该可充电电池可以是有线充电电池或无线充电电池。有线充电电池是通过有线线路充电的电池,无线充电电池是通过无线线圈充电的电池。该可充电电池还可以用于支持快充技术。
在一些实施例中,计算机设备800还包括有一个或多个传感器810。该一个或多个传感器810包括但不限于:加速度传感器811、陀螺仪传感器812、压力传感器813、指纹传感器814、光学传感器815以及接近传感器816。
加速度传感器811可以检测以计算机设备800建立的坐标系的三个坐标轴上的加速度大小。比如,加速度传感器811可以用于检测重力加速度在三个坐标轴上的分量。处理器801可以根据加速度传感器811采集的重力加速度信号,控制显示屏805以横向视图或纵向视图进行用户界面的显示。加速度传感器811还可以用于游戏或者用户的运动数据的采集。
陀螺仪传感器812可以检测计算机设备800的机体方向及转动角度,陀螺仪传感器812可以与加速度传感器811协同采集用户对计算机设备800的3D动作。处理器801根据陀螺仪传感器812采集的数据,可以实现如下功能:动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。
压力传感器813可以设置在计算机设备800的侧边框和/或显示屏805的下层。当压力传感器813设置在计算机设备800的侧边框时,可以检测用户对计算机设备800的握持信号,由处理器801根据压力传感器813采集的握持信号进行左右手识别或快捷操作。当压力传感器813设置在显示屏805的下层时,由处理器801根据用户对显示屏805的压力操作,实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。
指纹传感器814用于采集用户的指纹,由处理器801根据指纹传感器814采集到的指纹识别用户的身份,或者,由指纹传感器814根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时,由处理器801授权该用户执行相关的敏感操作,该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器814可以被设置在计算机设备800的正面、背面或侧面。当计算机设备800上设置有物理按键或厂商Logo时,指纹传感器814可以与物理按键或厂商Logo集成在一起。
光学传感器815用于采集环境光强度。在一个实施例中,处理器801可以根据光学传感器815采集的环境光强度,控制显示屏805的显示亮度。具体地,当环境光强度较高时,调高显示屏805的显示亮度;当环境光强度较低时,调低显示屏805的显示亮度。在另一个实施例中,处理器801还可以根据光学传感器815采集的环境光强度,动态调整摄像头组件806的拍摄参数。
接近传感器816,也称距离传感器,通常设置在计算机设备800的前面板。接近传感器816用于采集用户与计算机设备800的正面之间的距离。在一个实施例中,当接近传感器816检测到用户与计算机设备800的正面之间的距离逐渐变小时,由处理器801控制显示屏805从亮屏状态切换为息屏状态;当接近传感器816检测到用户与计算机设备800的正面之间的距离逐渐变大时,由处理器801控制显示屏805从息屏状态切换为亮屏状态。
本领域技术人员可以理解,图8中示出的结构并不构成对计算机设备800的限定,可以包括比图示更多或更少的组件,或者组合某些组件,或者采用不同的组件布置。
本申请实施例还提供了一种计算机可读存储介质,计算机可读存储介质中存储至少一条指令,至少一条指令由处理器加载并执行,以实现上述任一实现方式的束缚水饱和参数的确定方法所执行的操作。
本申请实施例还提供了一种计算机程序产品或计算机程序,计算机程序产品或计算机程序包括计算机程序代码,计算机程序代码存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取计算机程序代码,处理器执行计算机程序代码,使得计算机设备执行上述的束缚水饱和参数的确定方法所执行的操作。
在一些实施例中,本申请实施例所涉及的计算机程序可被部署在一个计算机设备上执行,或者在位于一个地点的多个计算机设备上执行,又或者,在分布在多个地点且通过通信网络互连的多个计算机设备上执行,分布在多个地点且通过通信网络互连的多个计算机设备可以组成区块链***。
本申请实施例提供了一种束缚水饱和参数的确定方法,该方法基于目标储层上标准井的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据来确定目标储层的黏土水体积系数和目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量,并基于目标储层的黏土水体积系数和目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量来确定具有第二元素俘获测井数据和第二常规测井数据的非标准井的目标储层的束缚水饱和参数,这样,避免了确定束缚水饱和参数时进行钻井取芯工作,省事省力,从而提高了确定束缚水饱和参数的效率。
以上仅为本申请的可选实施例,并不用以限制本申请,凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (10)
1.一种束缚水饱和参数的确定方法,其特征在于,所述方法包括:
确定待研究的目标储层和标准井的目标井段,所述标准井为所述目标储层上的井,所述目标井段为与所述目标储层的深度对应的井段;
获取所述目标井段在每个采样深度的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据;
基于所述目标井段在每个采样深度的核磁共振测井数据、所述第一元素俘获测井数据和所述第一常规测井数据,确定所述目标储层的黏土水体积系数和所述目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量;
获取所述目标储层上的非标准井在每个采样深度的第二元素俘获测井数据和第二常规测井数据;
对于所述非标准井的每个采样深度,基于所述第二元素俘获测井数据和所述第二常规测井数据、所述目标储层的黏土水体积系数和所述目标储层在所述采样深度的微孔隙束缚水体积系数和黏土含量,确定所述目标储层在所述采样深度的束缚水饱和参数。
2.根据权利要求1所述的束缚水饱和参数的确定方法,其特征在于,所述核磁共振测井数据包括总束缚水饱和参数、泥质束缚水饱和参数、泥质束缚水孔隙参数和第一总孔隙参数,所述第一元素俘获测井数据包括第一干黏土含量,所述第一常规测井数据包括泥质含量;
所述基于所述目标井段在每个采样深度的核磁共振测井数据、所述第一元素俘获测井数据和所述第一常规测井数据,确定所述目标储层的黏土水体积系数和所述目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量,包括:
对于所述目标储层的每个采样深度,确定所述干黏土含量与所述泥质束缚水孔隙参数之和,得到所述目标储层在所述采样深度的黏土含量;
基于所述目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数,确定所述目标储层的黏土水体积系数;
基于所述目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定所述目标储层的最小束缚水体积系数;
对于所述目标储层的每个采样深度,基于所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,确定所述目标储层在所述采样深度的分散粉砂含量;
基于所述目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数,确定所述目标储层的微孔隙束缚水参数;
对于所述目标储层的每个采样深度,基于所述微孔隙束缚水参数、所述黏土水体积系数、所述最小束缚水体积系数,以及所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,确定所述目标储层在所述采样深度的微孔隙束缚水体积系数。
3.根据权利要求2所述的束缚水饱和参数的确定方法,其特征在于,所述基于所述目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数,确定所述目标储层的黏土水体积系数,包括:
对所述目标储层在每个采样深度的黏土含量和泥质束缚水孔隙参数进行拟合,得到所述目标储层的黏土水体积系数。
4.根据权利要求2所述的束缚水饱和参数的确定方法,其特征在于,所述基于所述目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定所述目标储层的最小束缚水体积系数,包括:
基于所述目标储层在每个采样深度的总束缚水饱和参数和黏土含量,确定第一交会图,所述第一交会图的横坐标为黏土含量,纵坐标为总束缚水饱和参数;
确定所述第一交会图上的第一目标点和第二目标点,所述第一目标点的总束缚水饱和参数为最小值,所述第二目标点的总束缚水饱和参数和黏土含量均为最大值;
基于所述第一目标点和所述第二目标点,组成线性函数;
确定所述线性函数与所述纵坐标的交点处的值为所述最小束缚水体积系数。
5.根据权利要求2所述的束缚水饱和参数的确定方法,其特征在于,所述对于所述目标储层的每个采样深度,基于所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,确定所述目标储层在所述采样深度的分散粉砂含量,包括:
基于所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,通过第一关系数据,得到所述目标储层在所述采样深度的分散粉砂含量,所述第一关系数据的自变量为所述黏土水体积系数、所述目标储层在所述采样深度的泥质含量和干黏土含量,因变量为所述目标储层在所述采样深度的分散粉砂含量。
6.根据权利要求2所述的束缚水饱和参数的确定方法,其特征在于,所述基于所述目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数,确定所述目标储层的微孔隙束缚水参数,包括:
对所述目标储层在每个采样深度的总束缚水饱和参数、泥质束缚水饱和参数、黏土含量、干黏土含量和第一总孔隙参数进行拟合,得到所述微孔隙束缚水参数。
7.根据权利要求2所述的束缚水饱和参数的确定方法,其特征在于,所述对于所述目标储层的每个采样深度,基于所述微孔隙束缚水参数、所述黏土水体积系数、所述最小束缚水体积系数,以及所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,确定所述目标储层在所述采样深度的微孔隙束缚水体积系数,包括:
基于所述微孔隙束缚水参数、所述黏土水体积系数、以及所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,通过第二关系数据,得到所述目标储层在所述采样深度的第一微孔隙束缚水体积系数,所述第二关系数据的参数为所述微孔隙束缚水参数,自变量为所述黏土水体积系数和所述目标储层在所述采样深度的干黏土含量、分散粉砂含量和第一总孔隙参数,因变量为所述目标储层在所述采样深度的第一微孔隙束缚水体积系数;
若所述第一微孔隙束缚水体积系数在预设系数范围内,则所述第一微孔隙束缚水体积系数为所述微孔隙束缚水体积系数,所述预设系数范围的最小值为所述最小束缚水体积系数;
若所述第一微孔隙束缚水体积系数不在预设系数范围内,则基于所述预设系数范围,对所述第一微孔隙束缚水体积系数进行校正,得到所述微孔隙束缚水体积系数。
8.根据权利要求1所述的束缚水饱和参数的确定方法,其特征在于,所述第二元素俘获测井数据包括第二干黏土含量,所述第二常规测井数据包括第二总孔隙参数;
所述对于所述非标准井的每个采样深度,基于所述第二元素俘获测井数据和所述第二常规测井数据、所述目标储层的黏土水体积系数和所述目标储层在所述采样深度的微孔隙束缚水体积系数和黏土含量,确定所述目标储层在所述采样深度的束缚水饱和参数,包括:
基于所述黏土水体积系数和所述微孔隙束缚水体积系数,生成包含所述黏土水体积系数和所述微孔隙束缚水体积系数、且以所述第二干黏土含量、所述第二总孔隙参数和所述黏土含量为自变量、所述束缚水饱和参数为因变量的第三关系数据;
基于所述干黏土含量、所述第二总孔隙参数和所述黏土含量,通过所述第三关系数据,得到所述束缚水饱和参数。
9.一种束缚水饱和参数的确定装置,其特征在于,所述装置包括:
第一确定模块,用于确定待研究的目标储层和标准井的目标井段,所述标准井为所述目标储层上的井,所述目标井段为与所述目标储层的深度对应的井段;
第一获取模块,用于获取所述目标井段在每个采样深度的核磁共振测井数据、第一元素俘获测井数据和第一常规测井数据;
第二确定模块,用于基于所述目标井段在每个采样深度的核磁共振测井数据、所述第一元素俘获测井数据和所述第一常规测井数据,确定所述目标储层的黏土水体积系数和所述目标储层在每个采样深度的微孔隙束缚水体积系数和黏土含量;
第二获取模块,用于获取所述目标储层上的非标准井在每个采样深度的第二元素俘获测井数据和第二常规测井数据;
第三确定模块,用于对于所述非标准井的每个采样深度,基于所述第二元素俘获测井数据和所述第二常规测井数据、所述目标储层的黏土水体积系数和所述目标储层在所述采样深度的微孔隙束缚水体积系数和黏土含量,确定所述目标储层在所述采样深度的束缚水饱和参数。
10.一种计算机设备,其特征在于,所述计算机设备包括一个或多个处理器和一个或多个存储器,所述一个或多个存储器中存储有至少一条指令,所述至少一条指令由所述一个或多个处理器加载并执行以实现如权利要求1至权利要求8任一项所述的束缚水饱和参数的确定方法所执行的操作。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110812655.9A CN113586043B (zh) | 2021-07-19 | 2021-07-19 | 束缚水饱和参数的确定方法、装置和计算机设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110812655.9A CN113586043B (zh) | 2021-07-19 | 2021-07-19 | 束缚水饱和参数的确定方法、装置和计算机设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113586043A true CN113586043A (zh) | 2021-11-02 |
CN113586043B CN113586043B (zh) | 2023-10-27 |
Family
ID=78247895
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110812655.9A Active CN113586043B (zh) | 2021-07-19 | 2021-07-19 | 束缚水饱和参数的确定方法、装置和计算机设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113586043B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130261972A1 (en) * | 2010-11-18 | 2013-10-03 | Suncor Energy Inc. | Process for determining mobile water saturation in a reservoir formation |
CN105223116A (zh) * | 2015-08-28 | 2016-01-06 | 中国石油天然气集团公司 | 一种基于核磁共振谱系数法计算束缚水饱和度的方法 |
CN106706686A (zh) * | 2016-12-06 | 2017-05-24 | 东北大学 | 一种确定富有机质页岩水饱和度与水化损伤关系的方法 |
CN108252710A (zh) * | 2018-02-26 | 2018-07-06 | 中国石油天然气股份有限公司 | 一种确定油藏含水饱和度的方法及装置 |
US20190137648A1 (en) * | 2016-08-10 | 2019-05-09 | Petrochina Company Limited | Method and system for determining heterogeneous carbonate reservoir saturation exponent |
CN109753755A (zh) * | 2019-01-25 | 2019-05-14 | 中国石油天然气集团有限公司 | 一种确定储层含水饱和度的方法 |
CN110501750A (zh) * | 2019-09-19 | 2019-11-26 | 西北大学 | 一种提高致密砂岩油藏含油饱和度的解释精度的方法 |
CN112784404A (zh) * | 2021-01-05 | 2021-05-11 | 中国石油天然气集团有限公司 | 一种基于常规测井资料的砂砾岩束缚水饱和度计算方法 |
-
2021
- 2021-07-19 CN CN202110812655.9A patent/CN113586043B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130261972A1 (en) * | 2010-11-18 | 2013-10-03 | Suncor Energy Inc. | Process for determining mobile water saturation in a reservoir formation |
CN105223116A (zh) * | 2015-08-28 | 2016-01-06 | 中国石油天然气集团公司 | 一种基于核磁共振谱系数法计算束缚水饱和度的方法 |
US20190137648A1 (en) * | 2016-08-10 | 2019-05-09 | Petrochina Company Limited | Method and system for determining heterogeneous carbonate reservoir saturation exponent |
CN106706686A (zh) * | 2016-12-06 | 2017-05-24 | 东北大学 | 一种确定富有机质页岩水饱和度与水化损伤关系的方法 |
CN108252710A (zh) * | 2018-02-26 | 2018-07-06 | 中国石油天然气股份有限公司 | 一种确定油藏含水饱和度的方法及装置 |
CN109753755A (zh) * | 2019-01-25 | 2019-05-14 | 中国石油天然气集团有限公司 | 一种确定储层含水饱和度的方法 |
CN110501750A (zh) * | 2019-09-19 | 2019-11-26 | 西北大学 | 一种提高致密砂岩油藏含油饱和度的解释精度的方法 |
CN112784404A (zh) * | 2021-01-05 | 2021-05-11 | 中国石油天然气集团有限公司 | 一种基于常规测井资料的砂砾岩束缚水饱和度计算方法 |
Non-Patent Citations (2)
Title |
---|
付晨东;: "用核磁共振测井资料确定束缚水饱和度的新方法", 测井技术, no. 03, pages 243 - 246 * |
吴健;胡向阳;梁玉楠;汤翟;郑香伟;: "珠江口盆地低阻油层饱和度评价方法", 吉林大学学报(地球科学版), no. 01, pages 312 - 319 * |
Also Published As
Publication number | Publication date |
---|---|
CN113586043B (zh) | 2023-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112922589A (zh) | 尖灭线确定方法、装置、终端及存储介质 | |
CN114970919A (zh) | 钻速确定方法、装置、计算机设备及存储介质 | |
CN111077570B (zh) | 确定油藏中辉绿岩的侵入位置的方法、装置和存储介质 | |
CN111255434A (zh) | 气井的试井方法、装置和计算机存储介质 | |
CN110673214B (zh) | 预测水平井的入口靶点和末端点的深度的方法和装置 | |
CN116796503A (zh) | 一种页岩有机物碳含量计算模型建立方法和装置 | |
CN108733902B (zh) | 油藏渗透率的确定方法、装置及存储介质 | |
CN113586043B (zh) | 束缚水饱和参数的确定方法、装置和计算机设备 | |
CN109343113B (zh) | 一种预测油气藏位置的方法和装置 | |
CN113153284B (zh) | 束缚水饱和参数的确定方法、装置、设备和存储介质 | |
CN113216926A (zh) | 页岩气井的筛选方法、装置、终端及存储介质 | |
CN112096373B (zh) | 选定页岩气区域的方法、装置及计算机存储介质 | |
CN111812744B (zh) | 地层可压性确定方法、装置及计算机存储介质 | |
CN115680592B (zh) | 压裂有效性参数的确定方法、装置、设备及存储介质 | |
CN112987105B (zh) | 定量预测地下岩盐层分布的方法、装置、终端和存储介质 | |
CN111425189B (zh) | 超深层缝洞碳酸盐岩储层品质评价方法、装置及存储介质 | |
CN112698391B (zh) | 砂体刻画方法、装置、设备及存储介质 | |
CN114545507B (zh) | 取心深度确定方法、装置、终端和存储介质 | |
CN113175321B (zh) | 流体饱和参数的确定方法、装置和计算机设备 | |
CN110058306B (zh) | 确定三维速度体的方法、装置及计算机存储介质 | |
CN110630254B (zh) | 油水层的区分方法、装置及计算机可读存储介质 | |
CN115182718A (zh) | 识别水侵层的确定方法、装置、终端和存储介质 | |
CN114861961A (zh) | 落淤层的发育规模参数的确定方法和装置 | |
CN115478836A (zh) | 孔隙参数的阈值信息确定方法、装置、设备及存储介质 | |
CN115576009A (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 |