CN107944671A - 一种确定地层束缚水饱和度的方法 - Google Patents
一种确定地层束缚水饱和度的方法 Download PDFInfo
- Publication number
- CN107944671A CN107944671A CN201711084895.1A CN201711084895A CN107944671A CN 107944671 A CN107944671 A CN 107944671A CN 201711084895 A CN201711084895 A CN 201711084895A CN 107944671 A CN107944671 A CN 107944671A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- stratum
- rock
- clastic
- 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 112
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000011435 rock Substances 0.000 claims abstract description 149
- 238000004458 analytical method Methods 0.000 claims abstract description 80
- 239000004927 clay Substances 0.000 claims abstract description 36
- 238000002050 diffraction method Methods 0.000 claims abstract description 28
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 18
- 239000011707 mineral Substances 0.000 claims abstract description 18
- 238000011156 evaluation Methods 0.000 claims abstract description 12
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims abstract description 9
- 230000000704 physical effect Effects 0.000 claims abstract description 9
- 241001269238 Data Species 0.000 claims abstract description 8
- 239000002245 particle Substances 0.000 claims description 33
- 239000002734 clay mineral Substances 0.000 claims description 15
- 239000012530 fluid Substances 0.000 claims description 3
- 239000008187 granular material Substances 0.000 claims description 3
- 239000000843 powder Substances 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 2
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims description 2
- 241000894007 species Species 0.000 claims description 2
- 238000009736 wetting Methods 0.000 claims description 2
- 239000003921 oil Substances 0.000 description 25
- 239000004577 thatch Substances 0.000 description 11
- 239000000463 material Substances 0.000 description 6
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 5
- 229910052753 mercury Inorganic materials 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000005325 percolation Methods 0.000 description 3
- 239000004519 grease Substances 0.000 description 2
- 238000004321 preservation Methods 0.000 description 2
- 238000000205 computational method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 238000002447 crystallographic data Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 210000001624 hip Anatomy 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 229920002338 polyhydroxyethylmethacrylate Polymers 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Animal Husbandry (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Development Economics (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种确定地层束缚水饱和度的方法,该方法包括步骤:1)采集目标油藏岩心样本、地层水样本、地层油样本;2)对所采集的样本进行检测分析,获得物性分析数据、粒度分析数据、全岩X衍射分析数据、粘土矿物分析数据、地层可动水密度分析数据、地层油密度分析数据;3)根据在步骤2)中得到的全岩X衍射分析数据和粘土矿物分析数据及矿物骨架密度,用加权调和平均值确定岩石骨架密度,并进一步确定碎屑岩岩石体积密度、确定碎屑岩比表面积;同时确定碎屑岩地层束缚水膜厚度;再进一步确定碎屑岩地层束缚水饱和度。本发明为勘探初期或资料录取较少的碎屑岩油藏快速准确地确定地层束缚水饱和度提供了一种简便实用的评价方法。
Description
技术领域
本发明是关于一种涉及油藏地质参数的采集处理分析方法,具体而言是一种确定地层束缚水饱和度的方法,基于岩心分析应用中碎屑岩油藏地质参数,评价碎屑岩地层束缚水饱和度,属于地质勘探开发中岩心分析、流体分析等技术领域。
背景技术
在油气勘探开发领域,对碎屑岩油藏评价总离不开地层束缚水饱和度。碎屑岩地层束缚水饱和度与碎屑颗粒和孔隙喉道大小密切相关。从束缚水的赋存状态分析,地层束缚水饱和度受控于碎屑岩比表面积。在地层总孔隙度不变的前提下,碎屑岩比表面积越大,赋存在颗粒表面的束缚水就越多,地层束缚水饱和度也越高。碎屑颗粒直径越小,碎屑岩比表面积就越大。当碎屑岩比表面积一定时,地层束缚水膜越厚,地层束缚水饱和度就越高;碎屑颗粒表面越不规则(凹凸不平),束缚水膜就越厚。只要能准确求取碎屑岩比表面积、颗粒表面束缚水膜厚度、碎屑岩地层总孔隙度及碎屑岩体积密度等参数,就能准确计算碎屑岩地层束缚水饱和度。
目前确定碎屑岩地层束缚水饱和度主要用以下两类方法:(1)岩心分析法(又称直接法):包括压汞法、相渗法、密闭取心法;(2)测井解释法(又称间接法):包括Archie公式法、地区经验公式法、C/O比法等。其中直接法确定的碎屑岩地层束缚水饱和度精度要高于间接法确定的地层束缚水饱和度。
在实际应用中,还常常遇到含油区块或油藏中既没有压汞曲线、也没有相渗曲线、更没有岩心分析饱和度,只有碎屑岩物性分析数据、粒度分析数据、全岩X衍射分析数据及粘土矿物分析数据等资料,此时无法利用这些岩心分析资料确定碎屑岩地层束缚水饱和度。
发明内容
本发明的目的在于提供一种新的确定地层束缚水饱和度的方法,解决上述现有技术的方法有时无法确定碎屑岩地层束缚水饱和度的问题,为没有压汞曲线、没有相渗曲线、没有饱和度分析样品的区块或油藏提供一种新的碎屑岩地层束缚水饱和度评价方法。
为达上述目的,本发明提供了一种确定地层束缚水饱和度的方法,其包括以下步骤:
1)采集目标油藏岩心样本、地层水样本(地层可动水样本)、地层油样本;
2)对所采集的样本进行检测分析,获得物性分析数据、粒度分析数据、全岩X衍射分析数据、粘土矿物分析数据、地层可动水密度分析数据、地层油密度分析数据;
3)根据在步骤2)中得到的全岩X衍射分析数据和粘土矿物分析数据及矿物骨架密度,用加权调和平均值确定岩石骨架密度;
4)根据在步骤2)中得到的物性分析数据和地层可动水密度分析数据和步骤3)中得到的岩石骨架密度,确定碎屑岩岩石体积密度;
5)根据在步骤2)中得到的粒度分析数据和步骤4)中得到的岩石体积密度,确定碎屑岩比表面积;
6)根据在步骤2)中获得的地层可动水密度分析数据、地层油密度分析数据,结合实际油藏高度数据,确定碎屑岩地层束缚水膜厚度;
7)根据在步骤2)和步骤5)及步骤6)中得到的碎屑岩物性分析数据、碎屑岩比表面积、碎屑岩地层束缚水膜厚度,确定碎屑岩地层束缚水饱和度。
本发明的方法中,所述步骤2)与步骤6)顺序不分先后。
本发明是根据碎屑岩含油区块或油藏在既没有压汞曲线、也没有相渗曲线、更没有岩心分析饱和度时,可直接采用物性分析、粒度分析、全岩X衍射分析及粘土矿物分析等数据计算碎屑岩地层束缚水饱和度,提出一种全新的地层束缚水饱和度计算方法,尤其适用于录取的岩心资料少、不具备用常规方法计算束缚水饱和度基础的油藏在勘探开发各阶段确定碎屑岩地层束缚水饱和度的分析方法。
根据本发明的具体实施方案,本发明的方法中,步骤3)中所述矿物骨架密度详见表1。
表1矿物密度
根据本发明的具体实施方案,本发明的方法中,步骤3)中确定的岩石骨架密度符合以下公式:
式中J—粘土矿物分析的粘土矿物种类,类;
L—全岩X衍射分析的矿物种类(其中第L类为粘土类),类;
wcj—粘土矿物分析的第j(j=1、2、3、···、J)类粘土矿物质量分数,f;
wXL—全岩X衍射分析的粘土矿物质量分数,f;
wXl—全岩X衍射分析的第l(l=1、2、3、···、L)类矿物质量分数,f;
—碎屑岩中粘土矿物骨架密度平均值,g/cm3;
ρcj—粘土矿物分析的第j(j=1、2、3、···、J)类粘土矿物骨架密度,g/cm3;
—粘土矿物以外的碎屑岩矿物骨架密度平均值,g/cm3;
ρma—碎屑岩骨架密度,g/cm3;
ρXl—全岩X衍射分析的第l(l=1、2、3、···、L)类矿物骨架密度,g/cm3。
根据本发明的具体实施方案,本发明的方法中,步骤4)中确定的碎屑岩岩石体积密度符合以下公式:
ρb=ρma(1-φt)+ρwfφt (4)
式中φt—碎屑岩地层总孔隙度,f;
ρb—碎屑岩体积密度,g/cm3;
ρwf—地层可动水密度,g/cm3。
根据本发明的具体实施方案,本发明的方法中,步骤5)中确定的碎屑岩比表面积符合以下公式:
式中A—单位体积碎屑岩碎屑颗粒总表面积,μm2;
Ai—粒度分析的第i(i=1、2、3、4、5、6、7)类碎屑颗粒表面积,μm2;
di—粒度分析的第i类碎屑颗粒直径(d1<10代表泥质、10≤d2<50代表细粉砂岩、50≤d3<100代表粗粉砂岩、100≤d4<250代表细砂岩、250≤d5<500代表中砂岩、500≤d6<1000代表粗砂岩、1000≤d7<10000代表细砾岩),μm;
dk—粒度分析的最粗碎屑颗粒直径,μm;
Ni—粒度分析的第i(i=1、2、3、4、5、6、7)类碎屑颗粒数量,个;
Nk—粒度分析的最粗碎屑颗粒数量,个;
SA—单位体积碎屑岩比表面积,μm2/μm3;
Vma—单位体积碎屑岩骨架总体积(Vma=1-φt),μm3;
Vwi—100%含水单位体积碎屑岩束缚水体积(Vwi=1-Vma-Vwf),μm3;
Vwf—100%含水单位体积碎屑岩可动水体积(Vwi+Vwf=φt),μm3;
wi—粒度分析第i(i=1、2、3、4、5、6、7)类碎屑颗粒质量分数,f;
wk—粒度分析的最粗碎屑颗粒质量分数,f;
ρmai—粒度分析第i(i=1、2、3、4、5、6、7)类碎屑颗粒骨架密度,g/cm3;
ρmak—粒度分析的最粗碎屑颗粒骨架密度,g/cm3。
根据本发明的具体实施方案,本发明的方法中,所述di—粒度分析的第i类碎屑颗粒直径(d1<10代表泥质、10≤d2<50代表细粉砂岩、50≤d3<100代表粗粉砂岩、100≤d4<250代表细砂岩、250≤d5<500代表中砂岩、500≤d6<1000代表粗砂岩、1000≤d7<10000代表细砾岩)详见表2。
表2碎屑颗粒粒度分级(伍登-温特沃斯)
根据本发明的具体实施方案,本发明的方法中,步骤5)中,粒度分析碎屑颗粒质量分数与碎屑颗粒数量、碎屑颗粒直径及碎屑颗粒骨架密度符合以下公式:
dk=max(d1、d2、d3、d4、d5、d6、d7),(dk≥di) (10)
式中m—单位体积碎屑岩质量,g;
mi—粒度分析的第i(i=1、2、3、4、5、6、7)类碎屑颗粒质量,g;
mk—粒度分析的最粗碎屑颗粒质量,g。
根据本发明的具体实施方案,本发明的方法中,在实际应用中碎屑颗粒骨架密度分两部分:一部分为泥质(含粘土)骨架密度另一部分为砂级和砾级碎屑颗粒骨架密度
对地层条件下的单位体积碎屑岩,岩石质量等于岩石体积密度:
m=ρbV=ρb (11)
将式(11)代入式(9)得:
将式(5)、(6)、(7)与式(12)联立得碎屑岩比表面积经典方程:
根据本发明的具体实施方案,本发明的方法中,步骤6)中实际油藏高度数据可利用现有技术方法直接检测得到。
根据本发明的具体实施方案,本发明的方法中,步骤6)中确定的碎屑岩地层束缚水膜厚度符合以下公式:
式中:a—常数(a=1.18×10-7);
b—常数(b=0.5);
e—地层束缚水膜指数(e≥3~6);
H—自由水面以上碎屑岩油藏高度,m;
hwf—地层束缚水膜厚度,hwf>0.0025μm;
g0—标准重力加速度,g0=9.80665N/kg或g0=9.80665m/s2;
ρof—地层油密度,g/cm3;
ρwf—地层可动水密度,g/cm3;
θR—油藏条件下油水两相流体的润湿接触角θR=34°。
将相应参数代入式(14)得:
式(15)中的地层束缚水膜指数e随油藏高度增大而逐渐变小,在3~6之间变化。碎屑岩孔隙中的束缚水膜是由几十甚至几百个水分子“紧密”排列形成。根据已知单个水分子直径0.31~0.39nm,当10个水分子“紧密”排列形成的“束缚水膜厚度”为3.5nm(0.0035μm)。
在实际油藏中,碎屑岩地层束缚水膜分布在2.5~200nm之间。油藏高度越大,地层束缚水膜越薄,但最薄水膜不会小于5~10个水分子“紧密”排列的2.5nm。
采用油藏高度计算油藏顶部碎屑岩束缚水饱和度是油藏中最小的地层束缚水饱和度。
根据本发明的具体实施方案,本发明的方法中,步骤7)中确定的碎屑岩地层束缚水饱和度符合以下公式:
式中Swi—碎屑岩地层束缚水饱和度,f。
本发明的研究发现:
(1)当地层束缚水膜厚度一定时,粘土含量是碎屑岩地层束缚水饱和度最大影响因素,粘土矿物质量分数每增加一个百分点,地层束缚水饱和度至少增加3%;
(2)碎屑岩地层束缚水饱和度与碎屑岩地层总孔隙度成反比,与地层薄膜水厚度和比表面积乘积成正比,即
(3)根据全岩X衍射分析确定粘土在泥质中所占比例,当泥质中粘土占50%以上时,泥质颗粒直径取值小于等于1μm,当泥质中粘土占50%以下时,泥质颗粒直径取值大于1μm,其它颗粒的粒径取相应粒级下限值的120%;
(4)评价单元有多个分析样品时,其质量分数采用厚度加权平均值得到一组评价单元的典型数据,对这组典型数据重新计算每种矿物(粘土矿物分析和全岩X衍射分析)或每个粒级(粒度分析)的质量分数;
(5)附着在碎屑颗粒表面的束缚水膜不是一个固定值。在地层条件下,地层束缚水膜厚度是油藏高度、油水密度差(地层可动水密度减地层油密度)、碎屑颗粒润湿性等参数的函数。地层束缚水膜指数随油藏高度增大而减小,分布在3~6区间内;
(6)计算构造油藏平均地层束缚水饱和度时,参与计算的油藏高度取实际油藏高度的1/2;计算岩性油藏平均地层束缚水饱和度时,参与计算的油藏高度取实际油藏高度的1/4,更贴近实际。
从而,根据本发明的具体实施方案,本发明的方法还可包括根据以下一个或多个条件确定碎屑岩地层束缚水饱和度:
(1)粘土矿物质量分数每增加一个百分点,地层束缚水饱和度至少增加3%;
(2)碎屑岩地层束缚水饱和度与碎屑岩地层总孔隙度成反比,与地层薄膜水厚度和比表面积乘积成正比,即
(3)根据全岩X衍射分析确定粘土在泥质中所占比例,当泥质中粘土占50%以上时,泥质颗粒直径取值小于等于1μm,当泥质中粘土占50%以下时,泥质颗粒直径取值大于1μm,其它颗粒的粒径取相应粒级下限值的120%;
(4)评价单元有多个分析样品时,其质量分数采用厚度加权平均值得到一组评价单元的典型数据,对这组典型数据重新计算每种矿物(粘土矿物分析和全岩X衍射分析)或每个粒级(粒度分析)的质量分数;
(5)地层束缚水膜指数随油藏高度增大而减小,分布在3~6区间内;
(6)计算构造油藏平均地层束缚水饱和度时,参与计算的油藏高度取实际油藏高度的1/2;计算岩性油藏平均地层束缚水饱和度时,参与计算的油藏高度取实际油藏高度的1/4。
综上所述,本发明运用粒度分析、全岩X衍射分析、粘土矿物分析等资料,根据碎屑岩比表面积理论,从地层束缚水的赋存状态分析,地层束缚水饱和度受控于碎屑岩比表面积。在地层总孔隙度不变的前提下,碎屑岩比表面积越大,赋存在颗粒表面的地层束缚水就越多,地层束缚水饱和度也越高。碎屑颗粒直径越小,碎屑岩比表面积就越大。当碎屑岩比表面积一定时,地层束缚水膜越厚,地层束缚水饱和度就越高;碎屑颗粒表面越不规则(凹凸不平),束缚水膜就越厚。根据本发明的方法,能准确求取碎屑岩比表面积、颗粒表面地层束缚水膜厚度、碎屑岩地层总孔隙度及碎屑岩体积密度等参数,运用粒度等分析资料计算碎屑岩地层束缚水饱和度,特别是为勘探初期或资料录取较少的碎屑岩油藏快速准确地确定地层束缚水饱和度提供了一种简便实用的评价方法。且所测碎屑岩地层束缚水饱和度与直接法确定的地层束缚水饱和度精度相当。
附图说明
图1为100%含水碎屑岩单位体积模型。
具体实施方式
下列实施例仅为例示性说明,非意图限制本发明的范围。所属领域的技术人员可藉由本发明的揭露,在不背离本发明的精神的范围内做出适度的变更和修正。
本发明的用粒度资料计算碎屑岩地层束缚水饱和度的方法,需要注意以下方面:
一是在泥质中粘土占比需由全岩X衍射资料确定,当有多个分析样品时,需采用厚度加权平均值计算评价单元的典型全岩X衍射分析数据;
二是用实际粒度资料划分粒级,每个粒级的粒径取值(上限值、下限值、中间值)需薄片鉴定或实际岩心观察后确定,当评价单元有多个粒度分析样品时,每个粒级均采用厚度加权平均值计算评价单元的典型粒度分析数据;
三是评价单元地层总孔隙度采用与全岩X衍射分析及粘土矿物分析样品相对应的岩心分析孔隙度之厚度加权平均值;
四是岩石骨架密度需用全岩X衍射分析和粘土矿物分析数据及相应矿物骨架值,用加权调和平均值计算。
表3、表4和表5分别是茨602井的粘土矿物分析、全岩X衍射分析及粒度分析数据,本发明以茨602井所在茨602块相关资料为例阐述计算碎屑岩地层束缚水饱和度的具体步骤。
表3茨602井粘土矿物分析
表4茨602井全岩X衍射分析
表5茨602井粒度分析
茨602井位于茨602块腰部,属构造岩性油藏,油水界面-2260m、油藏高度160m、地层温度71.62℃、地层压力21.61MPa、地层水矿化度4845.7mg/L、地面原油密度0.8366g/cm3,茨602井2168m岩心分析孔隙度18.8%;茨602块地层原油密度0.7230g/cm3、地层水密度0.9900g/cm3、束缚水膜指数取6.0。茨602块油藏平均束缚水饱和度计算步骤:
第一步:根据粘土矿物分析、全岩X衍射分析和矿物骨架密度,用加权调和平均值计算岩石骨架密度。
第二步:根据岩心分析孔隙度和地层可动水密度计算岩石体积密度。
ρb=ρma(1-φt)+ρwfφt=2.72148×(1-0.188)+0.9900×0.188=2.395973g/cm3
第三步:根据粒度分析资料计算碎屑岩比表面积。
第四步:根据油水密度差和油藏高度计算地层束缚水膜厚度。
计算构造油藏平均地层束缚水饱和度采用二分之一油藏高度参加计算;计算岩性(或以岩性为主)油藏平均地层束缚水饱和度采用四分之一油藏高度参加计算,地层束缚水膜指数取6.0。
第五步:计算地层束缚水饱和度(表6)。
表6粒度资料计算茨602块束缚水饱和度
用粒度资料计算的地层束缚水饱和度为45.64%,与压汞资料计算的油藏束缚水饱和度46.30%相比精度相当。
Claims (9)
1.一种确定地层束缚水饱和度的方法,该方法包括步骤:
1)采集目标油藏岩心样本、地层水样本、地层油样本;
2)对所采集的样本进行检测分析,获得物性分析数据、粒度分析数据、全岩X衍射分析数据、粘土矿物分析数据、地层可动水密度分析数据、地层油密度分析数据;
3)根据在步骤2)中得到的全岩X衍射分析数据和粘土矿物分析数据及矿物骨架密度,用加权调和平均值确定岩石骨架密度;
4)根据在步骤2)中得到的物性分析数据和地层可动水密度分析数据和步骤3)中得到的岩石骨架密度,确定碎屑岩岩石体积密度;
5)根据在步骤2)中得到的粒度分析数据和步骤4)中得到的岩石体积密度,确定碎屑岩比表面积;
6)根据在步骤2)中获得的地层可动水密度分析数据、地层油密度分析数据,结合实际油藏高度数据,确定碎屑岩地层束缚水膜厚度;
7)根据在步骤2)和步骤5)及步骤6)中得到的碎屑岩物性分析数据、碎屑岩比表面积、碎屑岩地层束缚水膜厚度,确定碎屑岩地层束缚水饱和度。
2.根据权利要求1所述的方法,其中,步骤3)中确定的岩石骨架密度符合以下公式:
<mrow>
<msubsup>
<mover>
<mi>&rho;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>c</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>c</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>/</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<mfrac>
<msub>
<mi>w</mi>
<mrow>
<mi>c</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>c</mi>
<mi>j</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mover>
<mi>&rho;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>X</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>/</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mfrac>
<msub>
<mi>w</mi>
<mrow>
<mi>X</mi>
<mi>l</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>X</mi>
<mi>l</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>X</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>X</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mover>
<mi>&rho;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>c</mi>
<mo>&prime;</mo>
</msubsup>
<msub>
<mi>w</mi>
<mrow>
<mi>X</mi>
<mi>L</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>X</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中J—粘土矿物分析的粘土矿物种类,类;
L—全岩X衍射分析的矿物种类(其中第L类为粘土类),类;
wcj—粘土矿物分析的第j(j=1、2、3、…、J)类粘土矿物质量分数,f;
wXL—全岩X衍射分析的粘土矿物质量分数,f;
wXl—全岩X衍射分析的第l(l=1、2、3、…、L)类矿物质量分数,f;
—碎屑岩中粘土矿物骨架密度平均值,g/cm3;
ρcj—粘土矿物分析的第j(j=1、2、3、…、J)类粘土矿物骨架密度,g/cm3;
—粘土矿物以外的碎屑岩矿物骨架密度平均值,g/cm3;
ρma—碎屑岩骨架密度,g/cm3;
ρXl—全岩X衍射分析的第l(l=1、2、3、…、L)类矿物骨架密度,g/cm3。
3.根据权利要求1所述的方法,其中,步骤4)中确定的碎屑岩岩石体积密度符合以下公式:
ρb=ρma(1-φt)+ρwfφt (4)
式中φt—碎屑岩地层总孔隙度,f;
ρb—碎屑岩体积密度,g/cm3;
ρwf—地层可动水密度,g/cm3。
4.根据权利要求1所述的方法,其中,步骤5)中确定的碎屑岩比表面积符合以下公式:
<mrow>
<msub>
<mi>S</mi>
<mi>A</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>A</mi>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mi>A</mi>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&phi;</mi>
<mi>t</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>A</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>7</mn>
</munderover>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>&pi;</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>7</mn>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<msup>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>k</mi>
</mrow>
</msub>
<msup>
<msub>
<mi>d</mi>
<mi>k</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>i</mi>
</mrow>
</msub>
<msup>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式中A—单位体积碎屑岩碎屑颗粒总表面积,μm2;
Ai—粒度分析的第i(i=1、2、3、4、5、6、7)类碎屑颗粒表面积,μm2;
di—粒度分析的第i类碎屑颗粒直径(d1<10代表泥质、10≤d2<50代表细粉砂岩、50≤d3<100代表粗粉砂岩、100≤d4<250代表细砂岩、250≤d5<500代表中砂岩、500≤d6<1000代表粗砂岩、1000≤d7<10000代表细砾岩),μm;
dk—粒度分析的最粗碎屑颗粒直径,μm;
Ni—粒度分析的第i(i=1、2、3、4、5、6、7)类碎屑颗粒数量,个;
Nk—粒度分析的最粗碎屑颗粒数量,个;
SA—单位体积碎屑岩比表面积,μm2/μm3;
Vma—单位体积碎屑岩骨架总体积(Vma=1-φt),μm3;
Vwi—100%含水单位体积碎屑岩束缚水体积(Vwi=1-Vma-Vwf),μm3;
Vwf—100%含水单位体积碎屑岩可动水体积(Vwi+Vwf=φt),μm3;
wi—粒度分析第i(i=1、2、3、4、5、6、7)类碎屑颗粒质量分数,f;
wk—粒度分析的最粗碎屑颗粒质量分数,f;
ρmai—粒度分析第i(i=1、2、3、4、5、6、7)类碎屑颗粒骨架密度,g/cm3;
ρmak—粒度分析的最粗碎屑颗粒骨架密度,g/cm3。
5.根据权利要求1或4所述的方法,其中,步骤5)中,粒度分析碎屑颗粒质量分数与碎屑颗粒数量、碎屑颗粒直径及碎屑颗粒骨架密度符合以下公式:
<mrow>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
<mi>m</mi>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<msub>
<mi>&pi;&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>i</mi>
</mrow>
</msub>
<msup>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mn>6</mn>
<mi>m</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>m</mi>
<mi>k</mi>
</msub>
<mi>m</mi>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
<msub>
<mi>&pi;&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>k</mi>
</mrow>
</msub>
<msup>
<msub>
<mi>d</mi>
<mi>k</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mn>6</mn>
<mi>m</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
dk=max(d1、d2、d3、d4、d5、d6、d7),(dk≥di) (10)
式中m—单位体积碎屑岩质量,g;
mi—粒度分析的第i(i=1、2、3、4、5、6、7)类碎屑颗粒质量,g;
mk—粒度分析的最粗碎屑颗粒质量,g。
6.根据权利要求1或4所述的方法,其中,所述碎屑颗粒骨架密度为泥质(含粘土)骨架密度或者为砂级和砾级碎屑颗粒骨架密度
7.根据权利要求1所述的方法,其中,步骤6)中确定的碎屑岩地层束缚水膜厚度符合以下公式:
<mrow>
<msup>
<msub>
<mi>h</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
</msub>
<mi>e</mi>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mn>10</mn>
<mn>3</mn>
</msup>
<mi>a</mi>
</mrow>
<mrow>
<msub>
<mi>g</mi>
<mn>0</mn>
</msub>
<mi>H</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>o</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mi>b</mi>
<mrow>
<msub>
<mi>cos&theta;</mi>
<mi>R</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:a—常数(a=1.18×10-7);
b—常数(b=0.5);
e—地层束缚水膜指数(e≥3~6);
H—自由水面以上碎屑岩油藏高度,m;
hwf—地层束缚水膜厚度,hwf>0.0025μm;
g0—标准重力加速度,g0=9.80665N/kg或g0=9.80665m/s2;
ρof—地层油密度,g/cm3;
ρwf—地层可动水密度,g/cm3;
θR—油藏条件下油水两相流体的润湿接触角θR=34°。
8.根据权利要求1所述的方法,其中,步骤7)中确定的碎屑岩地层束缚水饱和度符合以下公式:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>w</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>S</mi>
<mi>A</mi>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>&phi;</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
式中Swi—碎屑岩地层束缚水饱和度,f。
9.根据权利要求1所述的方法,该方法还包括根据以下一个或多个条件确定碎屑岩地层束缚水饱和度:
(1)粘土矿物质量分数每增加一个百分点,地层束缚水饱和度至少增加3%;
(2)碎屑岩地层束缚水饱和度与碎屑岩地层总孔隙度成反比,与地层薄膜水厚度和比表面积乘积成正比,即
(3)根据全岩X衍射分析确定粘土在泥质中所占比例,当泥质中粘土占50%以上时,泥质颗粒直径取值小于等于1μm,当泥质中粘土占50%以下时,泥质颗粒直径取值大于1μm,其它颗粒的粒径取相应粒级下限值的120%;
(4)评价单元有多个分析样品时,其质量分数采用厚度加权平均值得到一组评价单元的典型数据,对这组典型数据重新计算每种矿物(粘土矿物分析和全岩X衍射分析)或每个粒级(粒度分析)的质量分数;
(5)地层束缚水膜指数随油藏高度增大而减小,分布在3~6区间内;
(6)计算构造油藏平均地层束缚水饱和度时,参与计算的油藏高度取实际油藏高度的1/2;计算岩性油藏平均地层束缚水饱和度时,参与计算的油藏高度取实际油藏高度的1/4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711084895.1A CN107944671B (zh) | 2017-11-07 | 2017-11-07 | 一种确定地层束缚水饱和度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711084895.1A CN107944671B (zh) | 2017-11-07 | 2017-11-07 | 一种确定地层束缚水饱和度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107944671A true CN107944671A (zh) | 2018-04-20 |
CN107944671B CN107944671B (zh) | 2021-11-02 |
Family
ID=61933421
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711084895.1A Active CN107944671B (zh) | 2017-11-07 | 2017-11-07 | 一种确定地层束缚水饱和度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107944671B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108414424A (zh) * | 2018-05-08 | 2018-08-17 | 中国石油天然气股份有限公司 | 一种确定储层边界层厚度的方法、装置及系统 |
CN110146940A (zh) * | 2019-04-04 | 2019-08-20 | 中国石油大学(华东) | 一种气藏束缚水饱和度、气藏可动水范围确定方法和装置 |
CN111693425A (zh) * | 2019-03-14 | 2020-09-22 | 中国石油大学(北京) | 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101806215A (zh) * | 2010-03-05 | 2010-08-18 | 中国石油集团川庆钻探工程有限公司 | 用束缚水饱和度数据判别储层流体类型的方法 |
-
2017
- 2017-11-07 CN CN201711084895.1A patent/CN107944671B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101806215A (zh) * | 2010-03-05 | 2010-08-18 | 中国石油集团川庆钻探工程有限公司 | 用束缚水饱和度数据判别储层流体类型的方法 |
Non-Patent Citations (4)
Title |
---|
JIHONG ZHANG 等: ""Experimental Research on Further Enhanced Oil Recovery by Using High DC Electric field"", 《2009 ASIA-PACIFIC POWER AND ENERGY ENGINEERING CONFERENCE》 * |
岳文正 等: ""用核磁共振及岩石物理实验求地层束缚水饱和度及平均孔隙半径"", 《测井技术》 * |
杨翠 等: ""低渗透油藏束缚水饱和度测定及影响因素分析"", 《广东化工》 * |
邢会民: ""基于常规测井资料的束缚水饱和度综合确定方法研究"", 《钻采工艺》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108414424A (zh) * | 2018-05-08 | 2018-08-17 | 中国石油天然气股份有限公司 | 一种确定储层边界层厚度的方法、装置及系统 |
CN111693425A (zh) * | 2019-03-14 | 2020-09-22 | 中国石油大学(北京) | 一种基于压汞曲线的岩心薄膜束缚水饱和度的测定方法 |
CN110146940A (zh) * | 2019-04-04 | 2019-08-20 | 中国石油大学(华东) | 一种气藏束缚水饱和度、气藏可动水范围确定方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107944671B (zh) | 2021-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104989392B (zh) | 一种岩性识别方法 | |
CN109653725B (zh) | 基于沉积微相和岩石相的混积储层水淹程度测井解释方法 | |
Masch et al. | Grain size distribution and its effect on the permeability of unconsolidated sands | |
Indraratna et al. | A new parameter for classification and evaluation of railway ballast fouling | |
CN104564041B (zh) | 基于开发渗透率下限的低渗透碎屑岩储层有效性评价方法 | |
Baniak et al. | Petrophysical characterization of bioturbated sandstone reservoir facies in the upper Jurassic Ula formation, Norwegian North Sea, Europe | |
Johnson et al. | Physical and hydrologic properties of water-bearing deposits in subsiding areas in central California | |
CN105134195A (zh) | 一种基于测井资料的页岩气储层品质评价方法 | |
Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
CN107102376A (zh) | 一种陆相湖盆致密油富集有利区综合评价与预测方法 | |
CN106988737B (zh) | 一种利用岩性组合识别沉积相的方法 | |
CN106223941B (zh) | 基于测井资料的页岩气储层有机孔隙度确定方法 | |
CN103869052A (zh) | 一种砂岩储层成岩测井相定量表征的方法 | |
CN105629327B (zh) | 一种针对弱胶结、深层致密砂岩储层成岩相定量表征方法 | |
CN109324345A (zh) | 叠合盆地油气成藏期岩石孔隙度恢复方法 | |
Pulley et al. | Conservatism of mineral magnetic signatures in farm dam sediments in the South African Karoo: the potential effects of particle size and post-depositional diagenesis | |
CN110533237A (zh) | 一种砂岩储层含油气产能预测方法 | |
Baiyegunhi et al. | Textural characteristics, mode of transportation and depositional environment of the Cretaceous sandstone in the Bredasdorp Basin, off the south coast of South Africa: Evidence from grain size analysis | |
CN107944671A (zh) | 一种确定地层束缚水饱和度的方法 | |
Díaz-Curiel et al. | New granulometric expressions for estimating permeability of granular drainages | |
Sahu et al. | Empirical methods and estimation of hydraulic conductivity of fluvial aquifers | |
CN105718720A (zh) | 一种复杂气藏储量品质分类综合评价方法 | |
CN108362621A (zh) | 基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法 | |
George et al. | Estimation of aquifer hydraulic parameters via complementing surfacial geophysical measurement by laboratory measurements on the aquifer core samples | |
Montafia | Influence of physical properties of marine clays on electric resistivity and basic geotechnical parameters |
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 |