CN105370269B - 一种基于压实校正等时地层界面追踪对比方法 - Google Patents
一种基于压实校正等时地层界面追踪对比方法 Download PDFInfo
- Publication number
- CN105370269B CN105370269B CN201510718585.5A CN201510718585A CN105370269B CN 105370269 B CN105370269 B CN 105370269B CN 201510718585 A CN201510718585 A CN 201510718585A CN 105370269 B CN105370269 B CN 105370269B
- Authority
- CN
- China
- Prior art keywords
- well
- thickness
- sandstone
- ancient
- mud stone
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000005056 compaction Methods 0.000 title claims abstract description 22
- 238000012937 correction Methods 0.000 title claims abstract description 14
- 239000004575 stone Substances 0.000 claims abstract description 69
- 238000011084 recovery Methods 0.000 claims abstract description 18
- 238000011160 research Methods 0.000 claims abstract description 12
- 239000004576 sand Substances 0.000 claims description 36
- 230000008021 deposition Effects 0.000 claims description 18
- 239000011148 porous material Substances 0.000 claims description 9
- 239000011435 rock Substances 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 20
- 239000010410 layer Substances 0.000 description 16
- 230000006870 function Effects 0.000 description 12
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 239000002356 single layer Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 239000002734 clay mineral Substances 0.000 description 3
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 230000000875 corresponding effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 241000761557 Lamina Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- GUJOJGAPFQRJSV-UHFFFAOYSA-N dialuminum;dioxosilane;oxygen(2-);hydrate Chemical compound O.[O-2].[O-2].[O-2].[Al+3].[Al+3].O=[Si]=O.O=[Si]=O.O=[Si]=O.O=[Si]=O GUJOJGAPFQRJSV-UHFFFAOYSA-N 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000002687 intercalation Effects 0.000 description 1
- 238000009830 intercalation Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- NLYAJNPCOHFWQQ-UHFFFAOYSA-N kaolin Chemical compound O.O.O=[Al]O[Si](=O)O[Si](=O)O[Al]=O NLYAJNPCOHFWQQ-UHFFFAOYSA-N 0.000 description 1
- 229910052622 kaolinite Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052901 montmorillonite Inorganic materials 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 229920002338 polyhydroxyethylmethacrylate Polymers 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000003643 water by type Substances 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于压实校正等时地层界面追踪对比方法，包括以下步骤：步骤1：选择一口标志层清晰、且研究层段明确的井作为标准井，以此井为参照井建立连井剖面，然后将标准井和待研究井的标志层顶部拉平，以建立等时基准参考面；步骤2：估算目的层砂岩、泥岩古恢复厚度；步骤3：根据步骤1中确定的等时基准参考面与步骤2估算的砂岩、泥岩古恢复厚度确定待分层井位的等时地层界面。本发明通过砂岩、泥岩差异压实作用来判断分析砂岩和泥岩实际的恢复厚度，从而消除了砂岩、泥岩差异压实作用的影响，最终实现了精准的等时界面的识别与对比。
Description
技术领域
本发明属于油田基础地质技术领域，涉及一种基于压实校正等时地层界面追踪对比方法。
背景技术
砂岩、泥岩在埋藏成岩过程中，普遍存在压实作用，物理实验模拟及现代沉积也证实了这一点。以泥岩为例，在与水接触的沉积表面，孔隙饱含水可达90％左右，当埋深达到0.9-3m时，泥岩中的孔隙水大量且快速排出，孔隙度迅速降低到60％作用，直至埋藏30m左右，孔隙度都一直保持稳定，如果随着埋深的进一步增加，则孔隙度又进一步遵循指数函数递减，直至递减到1-2％左右时，又再次保持稳定，而砂岩初始孔隙度要比泥岩低得多，主要原因是砂岩中饱含水的类型较单一，以孔隙水为主，因此，砂岩主要靠排出粒间孔隙水的方式减小岩性体积，而泥岩则不同，不仅有孔隙水，还有吸附水、层间水和结构水，这些水类型的存在，大大增加了泥岩中饱含的水量，加之泥岩中粘土矿物对水强吸附性，如高岭石和蒙脱石就是这类型矿物，它们在层间吸附了较多水分子，产生了大量层间水，随着压实强度的增强，这些水分子层会逐渐脱去，使泥岩中大量水得以排出，从而，大大减小了泥岩的体积。上述分析，都说明不论砂岩还是泥岩都存在不同程度压实作用，且差异性明显。
冲积平原的形成与河流在开阔地冲刷沉积作用密切相关，河流作用的结果形成了大面积的平原地貌，平坦的古地貌地形高差又较小，从而促使河流发生侧向蛇曲，流速骤降的结果导致水动力不足，使泥砂在下游大量沉积，这样便发生了大面积的填平补齐作用，因此，针对等时的某一段地层单元而言，沉积初始厚度平面变化不大，可近似认为等厚，也就是说，不论砂岩还是泥岩，在沉积物未成岩近地表时，沉积的厚度整体是近似相等的，然而，地表沉积物在逐渐埋藏过程中，由于上覆地层厚度的增加，地层压力是逐渐变大的，因此，砂岩、泥岩要经受强烈的压实作用，并通过排出饱含水的方式，减少岩性体积。然而，这个过程，砂岩、泥岩表现是不同的，原因主要是砂岩成份和结构不同，砂岩体积的压缩量主要依靠粒间孔隙的减小，而泥岩则由于其成份主要为粘土矿物，粘土矿物具有较强的吸附力，会形成多水分层的层间水，如蒙脱石最多可有8个水分子层，在逐渐压实过程中，这些水分子层会逐层脱去，分子间排列越来越紧密，从而导致泥岩越来越致密，水分子大量排出，加上原始粒间孔隙水的大量排出，因此泥岩的体积的损失量要远远大于砂岩，由此便造成了深埋地下的砂岩、泥岩存在差异压实作用，鉴于上述认识，在进行深层地层对比过程中如果忽略了这种差异作用，必然导致等时地层单元界面的误判，因此，基于考虑砂岩、泥岩差异压实作用的冲积平原地层的对比，将具有非常重要意义。
发明内容
本发明的目的在于解决上述现有技术存在的缺陷，提供一种基于压实校正等时地层界面追踪对比方法。
本发明的目的是通过如下措施来达到的：一种基于压实校正等时地层界面追踪对比方法，其特征在于，包括以下步骤：
步骤1：选择一口标志层清晰、且研究层段明确的井作为标准井，以此井为参照井建立连井剖面，然后将标准井和待研究井的标志层顶部拉平，以建立等时基准参考面；
步骤2：估算待研究井目的层的砂岩、泥岩古恢复厚度；
步骤3：根据步骤1中确定的等时基准参考面与步骤2估算的砂岩、泥岩古恢复厚度确定待研究井目的层的等时地层界面。
在上述技术方案中，所述标志层须同时具备以下条件：
①该界面具有明显特征可作为地层对比标志的岩层；②具有所含化石和岩性特征明显；③层位稳定且分布范围广；④特征明显且易于鉴别。
在上述技术方案中，所述估算待研究井目的层砂岩、泥岩古恢复厚度包括以下步骤：
步骤a、首先求取目的层压实率η，
η＝(H古-H今)/H古 (1)
由于根据遵循骨架体积不变原理，有:
根据公式(1)、(2)求的公式(3)
其中，H今为现今埋深某一单一岩性段厚度，单位，m；
H古为古沉积期近地表某一单一岩性段厚度，单位，m；
η为目的层段某一岩性压实率，无量纲；
为某一段岩性现今埋深平均孔隙度值，单位，％
为目的层初始孔隙度，单位，％；
步骤b、利用公式(1)分别估算砂岩、泥岩的古恢复厚度，其求取公式为：
H古砂＝H今砂/(1-η砂) (4)
H古泥＝H今泥/(1-η泥) (5)
H今砂为现今埋深砂岩层段厚度，单位，m；
H今泥为现今埋深泥岩层段厚度，单位，m；
H古砂为古沉积期近地表砂岩层段厚度，单位，m；
H古泥为古沉积期近地表泥岩层段厚度，单位，m。
在上述技术方案中，当为砂岩时，初始孔隙度值为40％-50％；当为泥岩时，初始孔隙度为50％-70％。
在上述技术方案中，所述可表示为：
相对于现有技术，本发明的有益效果为：本发明通过砂岩、泥岩差异压实作用来判断分析砂岩和泥岩实际的恢复厚度，从而消除了砂岩、泥岩差异压实作用的影响，最终实现了精准的等时界面的识别与对比，从而有助于揭示砂体空间叠置与局部砂体富集规律，提高砂体的钻遇率，此外，还有助于指导连井剖面两口邻井砂体是否同期问题，便于后期注水开发、井网调整等生产作业的实施，为油田生产实践提供理论依据。
附图说明
图1(a)为恢复前目的层今埋深连井剖面模型图；
图1(b)为图1(a)经压实校正后的连井剖面模型图；
图2(a)为标志层顶拉平连井剖面图；
图2(b)为去压实校正后沉积厚度恢复剖面图；
图2(c)为经去压实效应地层对比之后再恢复至目前埋深的连井剖面图；
图3砂岩孔隙度-深度拟合关系图；
图4泥岩孔隙度-深度拟合关系图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明提供的一种基于压实校正等时地层界面追踪对比方法，探索出了一套基于“小层顶部标志层拉平法”、“砂岩、泥岩古厚度恢复法”和“标准井等厚平切法”，总结了一系列相应的换算公式，消除了砂岩、泥岩差异压实作用的影响，最终实现了精准的等时界面的识别与对比。本发明按照以下步骤进行：
步骤1：选择一口标志层清晰、且研究层段明确的井作为标准井，以此井为参照井建立连井剖面，然后将标准井和待研究井的标志层顶部拉平，以建立等时基准参考面；
具体地，选择一口标准井，这口井要满足标志层清晰，目的层段明确，井数据完备等条件，即选择一口本身天然具有标志层，而且研究层段明确的一口井，作为参考基准井，其它井参考该井，将标志层依次连接，这样就建好了待研究剖面，如图2a所示)，然后以此井为参照井建立连井剖面，然后将标准井和目的层的标志层顶部拉平，以建立等时基准参考面；
标准井的选择是进行与目标井的等时地层单元对比的前提和基础，为本发明的正确实施建立了基准参考界限，并从该标准井出发，建立要实施的连井剖面。标志层顶拉平是进行小层等时对比的基础，只有将该等时界线拉齐，才可以为后续厚度恢复以及进行小层等时对比建立基础。本发明具体操作方法是将小层顶界标志层与目的层的顶界标志层(所谓标志层是岩性特征明显，横向分布稳定，易于识别的一个等时界面)拉平，该拉平界线是水平的，古沉积填平补齐末期地表也是近水平的，这就为后续将上覆地层全部地层去掉之后，将目的层松弛厚度恢复，然后为采用等厚原则进行刻度其它目标井的层位追踪创造了条件和研究基础，如图1所示，其中图1(a)为恢复前目的层今埋深连井剖面模型，图1(b)为经压实校正后，恢复各井目的层近地表时古厚度，各井厚度大致相等模型，其中，1为单层砂体，2为单层泥岩，3为恢复古厚度的单层泥岩，4为恢复古厚度的单层砂体，5为去压实之前易出现的错误地层对比界线，6为去压实后等时界线。在所有工作进行之前，要建立层位识别标准，也就是说其它井的分层必须按这个标准执行，然后以此为基础，将邻井依次连接起来，建立连井剖面，作为研究对象。标志层应是绝对等时的界面，且一般距离目的层较近，该标志层的选取必须具备如下特征：①该界面必须具有明显特征可作为地层对比标志的岩层；②应当具有所含化石和岩性特征明显；③层位稳定且分布范围广；④特征明显且易于鉴别。将具有上述特征的标志层在连井剖面中逐一识别出来，然后将其依次连接起来，拉齐至同一水平线上，最终建立具有绝对等时意义的基准参考面，如图1(a)中8所示。
步骤2：估算目的层砂泥压实率与砂岩、泥岩古恢复厚度；
要进行目的层砂岩、泥岩压实率的估算，首先必须要拟合求出待分层目标工区的孔隙度-埋深关系函数，以便求解出现今埋深条件下的残余孔隙度，然后遵循骨架体积不变原理为基础，构建古今厚度函数，有如下公式(2)，
从而，便可以推导出压实率求取公式，
H今为现今埋深某一单一岩性段厚度，单位，m；
H古为古沉积期近地表某一单一岩性段厚度，单位，m；
η为目的层段某一岩性压实率，无量纲；
为某一段岩性现今埋深平均孔隙度值，单位，％
为初始孔隙度值，即：目的层恢复到地表时的孔隙度值，单位，％；
通过量取目的层砂岩、泥岩层段的厚度值，利用公式(3)求取的压实率换算公式，即可分别求出砂岩、泥岩的古厚度，求取公式如下
H古砂＝H今砂/(1-η砂) (4)
H古泥＝H今泥/(1-η泥) (5)
H今砂为现今埋深砂岩层段厚度，单位，m；
H今泥为现今埋深泥岩层段厚度，单位，m；
H古砂为古沉积期近地表砂岩层段厚度，单位，m；
H古泥为古沉积期近地表泥岩层段厚度，单位，m；
η为目的层段某一岩性压实率，无量纲；
公知的，碎屑岩层在埋深逐渐增加过程中，孔隙度随深度增加而呈现指数减小,即在一定深度范围内，任意深度的孔隙度可表示为：
其中，为初始孔隙度，单位，％；
C为压实系数，无量纲；
Z为埋藏深度，单位，m
如果有实测从近地表到目的层埋深孔隙度值，按公式(6)模型进行相关拟合，即可求出初始孔隙度和压实系数C，而初始孔隙度一般拟合的结果偏大，尤其是泥岩的拟合结果，可达到80-90％，近地表的数据点要去掉，因为在近地表0-3m内，泥岩中的饱含水可以迅速降低到60％，然后再遵循指数函数规律式递减，而砂岩根据现代沉积测量和实验室模拟结果，砂岩一般40％-50％，所以借用公认的泥岩在地表的孔隙度经验值，作为本发明的初始孔隙度。由于本发明对象为河流相地层，河道砂岩中多有泥质夹层存在，因此砂岩初始孔隙度取45％，泥岩初始孔隙度取60％作为研究区拟合孔隙度-深度函数时的校正值。这样，在上述两个系数求取的基础上，还要满足如下假设条件：①骨架体积不变原理，也就是随着埋深的增加，骨架压缩的体积量要忽略不计；②目的层由于顶底高程较小，其顶底压缩率可认为近似相等，整改层段用统一的孔隙度值进行估算，也就是说目的层遵循统一的深度-孔隙度函数关系，在满足上述两个条件下，便可获得压实率η求取公式，其求取过程如下：
首先拟合好的孔隙度-深度函数关系，结果见图3和图4，然后将目的层埋深代入公式(6)，估算出不同岩性目的层，进而分别估算求取出砂岩、泥岩的古恢复厚度；即：
η＝(H古-H今)/H古＝(φ0-φ目的层)/(1-φ目的层)
在压实率求取基础上，古厚度求取按如下公式获得，
H古砂＝H今砂/(1-η砂) (7)
H古泥＝H今泥/(1-η泥) (8)
H今砂为现今埋深砂岩层段厚度，单位，m；
H今泥为现今埋深泥岩层段厚度，单位，m；
H古砂为古沉积期近地表砂岩层段厚度，单位，m；
H古泥为古沉积期近地表泥岩层段厚度，单位，m；
η为目的层段某一岩性压实率，无量纲；
量取连井剖面中所有井目标层位内各砂岩、泥岩层厚度，然后将砂岩、泥岩各层段厚度分别代入公式(7)和公式(8)中，即可分别求出H古砂和H古泥。
步骤3：根据步骤1中确定的等时基准参考面与步骤2估算的砂岩、泥岩古恢复厚度确定待分层井位的等时地层界面；
以步骤2为基础，分别将其它井标志层下部所有的砂岩、泥岩段恢复古厚度，如图1(b)所示，遵循目的层在近地表等厚为原则，利用标准井底界面“标准井等厚平切法”(所谓标准井等厚平切法，就是在恢复古厚度的剖面中，参照标准井目的层厚度规模，遵循等厚为原则，横向水平卡取其它目标井位，并记录该位置，然后再将所有井的目的层段厚度复原，原先刻度记录位置便保留下来，然后将此位置依次连接)最终实现了所有井等时界面的追踪和对比，结果如图1(b)6所示。
下面列举具体实施例对本发明进行说明：
以大庆油田葡萄花油层组PI砂层组地层单元的实施作为实施例。
实施例1：结合步骤1以标准井为关键井的连井剖面建立与目的层顶拉平、步骤2的基于深度-孔隙度函数关系拟合为基础的目的层压实率的求取与古厚度的恢复和步骤3采用标准井等厚平切法，最终实现了所有井等时界面的追踪和对比。应用本发明提出的一整套小层精确对比方法，按如下步骤完成：
1、以标准井为基础建立连井剖面并标志层顶部拉平
为了验证本发明方法的可靠性，以实际连井剖面为研究对象，在研究目标区选择了以X1-1-627井作为标准井，其标准层位顶底埋深分别为1064.5m和1093.5m，现今埋藏厚度为29m，然后以该标准井为基点，将其要分层的邻井依次连接起来，形成如图2(a)所示的连井剖面，如图2(a)-2(b)所示，图2(b)中，1为古沉积期近地表参考面，2为各井恢复古厚度后目的层底界，3为测井曲线微梯度，4为测井曲线微电位，5为测井曲线自然电位，图2(a)中，6为标志层，8为标准井目的层底界参考基准面；图2(c)中，7为去压实刻度等时界面，8为标准井目的层底界参考基准面。然后按X1-1-627井层段顶部1064.5m处电测曲线特征明显，电阻率曲线呈刺刀状，为高阻与低幅自然电位曲线特征，其它井位如X1-D2-P30井、X1-D2-P32井和X1-D2-P33井，分别在1062.6m、1069.5m和1073.9m处有与此类似的电测曲线特征，且埋深相差不大，可以作为标准层，依次连接起来，然后拉齐至同一水平线，目的是模拟地表沉积末期状态，为后期厚度恢复至地表建立研究基础。当然，小层顶部分层已确定，剩下就是最为关键也是最难的步骤，就是所有待分层底界的确定，后面本发明的相关步骤将给予完成和实施。
2、砂泥压实率估算与各层古厚度恢复
压实率的求取关键是要拟合出适合本研究区的孔隙度-深度变化规律函数，根据砂岩、泥岩孔隙度实验室分析数据，按照公式(5)作为函数模板进行拟合，拟合结果见图3和图4，得出砂岩、泥岩不同埋深孔隙度求取函数，其中
砂岩拟合结果为：
泥岩拟合结果为：
由于目的层段厚度30m左右，对于埋深近1000m地层而言，待分层的顶底孔隙度相差不大，可以忽略不计，因此可以简单地将目的层段中间深度位置的孔隙度值代替整个层段的孔隙度，然后分别读取各井大致中间位置的深度，然后取平均值，求取结果为1085m，这个误差只要保持在±30m，便不会影响判别结果，因为以砂岩为例，将1085m和1115m深度代入到公式(9)中，其两者差值。
由此可见，在埋深变化不大时，孔隙度值变化是很小的，所以，也就有了目的层段遵循统一深度-孔隙度函数变化关系的假设。因此，取1085埋深作为目的层埋深，分别代入到公式(9)和公式(10)中，然后将计算结果再分别代入公式(6)中，求得目的层砂岩和泥岩压实率分别为：
利用上述计算求得的η砂和η泥值，以及量取的砂岩、泥岩层厚度，见表1，分别代入古厚度恢复公式(7)和公式(8)中，计算结果见表2。
表1 砂岩、泥岩层人工读取厚度值
表2 目标层段内砂岩、泥岩层段古厚度恢复数据
3、等时地层界面的确定
根据计算结果来看，X1-1-627井恢复总厚度为38.4m，与X1-D2-P30井恢复总厚度的36.6m，X1-D2-P32井恢复的总厚度的37m，以及X1-D2-P33井恢复的总厚度的36.5m处是近等厚的，因此，在误差控制在±2m允许范围内(对于层厚30m来说，是可以忽略不计的)，根据恢复地表近等厚假设为条件，以及一般地层对比原则(等时界线一般位于岩性突变位置，如泥岩与砂岩的交界处)，标准井X1-1-627井目的层底界等时界线就应标定在其它井所恢复厚度的底界位置，详细见图2(b)，然后将该位置做标记，然后将恢复的厚度还原，那么标定的位置也就会随井剖面的压缩而上提，形成如图2(c)所示结果，7所示位置便是要刻度标定的等时界面。
本发明充分利用岩心、测录井、野外露头和现代卫星图片等资料，在标志层拉平基础上，通过“小层顶部标志层拉平法”、“砂岩、泥岩厚度恢复法”和“标准井等厚平切法”，结合相应的换算公式，最终完成了小层精确的等时对比，避免了因忽略砂岩、泥岩差异压实作用而造成的错误的地层等时对比。研究结果表明，该研究方法将有助于揭示砂体空间叠置与局部砂体富集规律，提高砂体的钻遇率，此外，还有助于指导连井剖面两口邻井砂体是否同期问题，便于后期注水开发、井网调整等生产作业的实施，为油田生产实践提供理论依据。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (3)
1.一种基于压实校正等时地层界面追踪对比方法，其特征在于，包括以下步骤：
步骤1：选择一口标志层清晰、且研究层段明确的井作为标准井，以此井为参照井建立连井剖面，然后将标准井和待研究井的标志层顶部拉平，以建立等时基准参考面；
步骤2：估算待研究井目的层的砂岩、泥岩古恢复厚度；
步骤3：根据步骤1中确定的等时基准参考面与步骤2估算的砂岩、泥岩古恢复厚度，然后利用标准井底界面“标准井等厚平切法”，确定待研究井目的层的等时地层界面；
所述“标准井等厚平切法”，就是在恢复古厚度的剖面中，参照标准井目的层厚度规模，遵循等厚为原则，横向水平卡取其它目标井位，并记录该位置，然后再将所有井的目的层段厚度复原，原先刻度记录位置便保留下来，将此位置依次连接；
所述标志层须同时具备以下条件：
①该界面具有明显特征可作为地层对比标志的岩层；②具有所含化石和岩性特征明显；③层位稳定且分布范围广；④特征明显且易于鉴别；
所述估算待研究井目的层砂岩、泥岩古恢复厚度包括以下步骤：
步骤a、首先求取目的层压实率η，
η＝(H古-H今)/H古 (1)
由于根据遵循骨架体积不变原理，有:
根据公式(1)、(2)求的公式(3)
其中，H今为现今埋深某一单一岩性段厚度，单位，m；
H古为古沉积期近地表某一单一岩性段厚度，单位，m；
η为目的层段某一岩性压实率，无量纲；
为某一段岩性现今埋深平均孔隙度值，单位，％
为目的层初始孔隙度，单位，％；
步骤b、利用公式(1)分别估算砂岩、泥岩的古恢复厚度，其求取公式为：
H古砂＝H今砂/(1-η砂) (4)
H古泥＝H今泥/(1-η泥) (5)
H今砂为现今埋深砂岩层段厚度，单位，m；
H今泥为现今埋深泥岩层段厚度，单位，m；
H古砂为古沉积期近地表砂岩层段厚度，单位，m；
H古泥为古沉积期近地表泥岩层段厚度，单位，m。
2.根据权利要求1所述的一种基于压实校正等时地层界面追踪对比方法，其特征在于，当为砂岩时，初始孔隙度值为40％-50％；当为泥岩时，初始孔隙度为50％-70％。
3.根据权利要求1所述的一种基于压实校正等时地层界面追踪对比方法，其特征在于，所述可表示为：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510718585.5A CN105370269B (zh) | 2015-05-13 | 2015-10-28 | 一种基于压实校正等时地层界面追踪对比方法 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2015102424440 | 2015-05-13 | ||
| CN201510242444 | 2015-05-13 | ||
| CN201510718585.5A CN105370269B (zh) | 2015-05-13 | 2015-10-28 | 一种基于压实校正等时地层界面追踪对比方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105370269A CN105370269A (zh) | 2016-03-02 |
| CN105370269B true CN105370269B (zh) | 2018-11-13 |
Family
ID=55372827
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510718585.5A Active CN105370269B (zh) | 2015-05-13 | 2015-10-28 | 一种基于压实校正等时地层界面追踪对比方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105370269B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107091084A (zh) * | 2017-04-28 | 2017-08-25 | 中国石油集团川庆钻探工程有限公司 | 建立虚拟数据井提取页岩气水平井储层物性参数的方法 |
| CN109577944B (zh) * | 2018-11-28 | 2022-10-21 | 王国亭 | 多层系低渗致密砂岩气井各层动态控制储量测定方法 |
| CN111650637B (zh) * | 2019-03-04 | 2022-12-02 | 中国石油天然气股份有限公司 | 地震层位解释方法及装置 |
| CN110656934B (zh) * | 2019-10-08 | 2024-04-30 | 中国石油天然气股份有限公司 | 一种致密砂岩储层去压实地层对比方法 |
| CN111142171A (zh) * | 2019-12-30 | 2020-05-12 | 核工业北京地质研究院 | 一种基于去压实作用批量化计算盆地埋藏史的方法 |
| CN111766630B (zh) * | 2020-07-02 | 2022-04-22 | 中国地质大学（北京） | 一种盆地古地貌的恢复方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104533402A (zh) * | 2014-11-25 | 2015-04-22 | 长江大学 | 一种多期叠置复合辫状河道砂体分期厘定方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101070759B (zh) * | 2007-05-04 | 2010-05-19 | 中国矿业大学 | 一种薄基岩浅埋煤层保水开采适用条件分类方法 |
| CN103510947B (zh) * | 2012-06-21 | 2016-05-25 | 中国石油化工股份有限公司 | 滩坝砂岩微相识别模式的建立方法 |
| CN104570156A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种三角洲前缘相沉积地层底面等高程精度对比方法 |
| CN104122602A (zh) * | 2014-07-31 | 2014-10-29 | 中国石油大学(北京) | 一种获取沉积盆地剥蚀厚度的方法及装置 |
| CN104500051A (zh) * | 2014-11-11 | 2015-04-08 | 中国石油天然气集团公司 | 一种确定水淹层混合地层水宏观俘获截面的方法 |
-
2015
- 2015-10-28 CN CN201510718585.5A patent/CN105370269B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104533402A (zh) * | 2014-11-25 | 2015-04-22 | 长江大学 | 一种多期叠置复合辫状河道砂体分期厘定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105370269A (zh) | 2016-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105370269B (zh) | 一种基于压实校正等时地层界面追踪对比方法 | |
| CN107991714B (zh) | 基于湖盆古地形恢复的定量方法 | |
| Flathe | Possibilities and limitations in applying geoelectrical methods to hydrogeological problems in the coastal areas of north west Germany | |
| CN111766630B (zh) | 一种盆地古地貌的恢复方法 | |
| CN101936165B (zh) | 一种岩溶型碳酸盐岩储层布井方法 | |
| CN104765065A (zh) | 一种确定沉积地层古地貌的方法及装置 | |
| Tizro et al. | Estimation of porosity and specific yield by application of geoelectrical method–a case study in western Iran | |
| CN105425315B (zh) | 一种弱构造变形区小尺度沉积期古地貌反演方法 | |
| CN103206207A (zh) | 基于生烃期古构造的油藏有利区的识别方法 | |
| CN104453834A (zh) | 一种井组注采关系优化调整方法 | |
| CN104076394B (zh) | 一种碳酸盐岩潜山洞穴分层解释方法 | |
| CN109541685B (zh) | 一种河道砂体识别方法 | |
| CN111624651B (zh) | 基于古地貌约束的储层预测方法及装置 | |
| Wimmers et al. | Integration of sedimentology, petrophysics and rock typing as key to understanding a tight gas reservoir | |
| Khalaf et al. | Estimation the hydrogeological characteristics of Al-Dammam confined aquifer in the west and southwest of Karbala city, Iraq | |
| Whitaker et al. | Simulating flow through forward sediment model stratigraphies: insights into climatic control of reservoir quality in isolated carbonate platforms | |
| Liu et al. | The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs | |
| CN109814174B (zh) | 一种碎屑岩不整合结构体测井综合识别方法 | |
| CN102134994B (zh) | 一种基于基岩油藏油水层电阻率的地层数据处理方法 | |
| Hitchon et al. | Hydrodynamics and hydrocarbon occurrences, Surat Basin, Queensland, Australia | |
| De Luca et al. | Effectiveness of geophysical surveys for water wells relocation in overexploited aquifers (the example of Maggiore and Traversola Valleys, Northwestern Italy) | |
| CN109212625B (zh) | 一种砂岩型铀矿古地下水交替强度定量计算方法 | |
| Hjuler et al. | Multidisciplinary approach to assess geothermal potential, Tønder area, North German Basin | |
| CN113433589A (zh) | 一种基于数理统计的风化壳储层底界面识别方法 | |
| Metwalli et al. | Reservoir petrophysical modeling and risk analysis in reserve estimation; a case study from Qasr Field, North Western Desert, Egypt |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220819 Address after: Room 801-82, Unit B, 8th Floor, Building 33, No. 99, Kechuang 14th Street, Beijing Economic and Technological Development Zone, Beijing 100000 (Beijing Pilot Free Trade Zone High-end Industry Zone Yizhuang Group) Patentee after: Beijing furuibao Energy Technology Co.,Ltd. Address before: School of Earth Sciences, No. 111 Daxue Road, Caidian District, Wuhan City, Hubei Province 430100 Patentee before: Yangtze University |