CN104036119B - 划分沉积地层的方法 - Google Patents
划分沉积地层的方法 Download PDFInfo
- Publication number
- CN104036119B CN104036119B CN201410209502.5A CN201410209502A CN104036119B CN 104036119 B CN104036119 B CN 104036119B CN 201410209502 A CN201410209502 A CN 201410209502A CN 104036119 B CN104036119 B CN 104036119B
- Authority
- CN
- China
- Prior art keywords
- reflection coefficient
- interval
- curve
- interface
- log
- 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.)
- Expired - Fee Related
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于油气地质勘探领域,涉及一种划分沉积地层的方法。本发明运用地震地层学研究方法,利用测井曲线拟合地震反射系数曲线,进而通过沉积界面的反射系数特征划分沉积地层,并进行对比。本发明填补了因缺少区域标志层而无法划分与对比地层的科研现状,可以弥补地震划分垂向层序分辨率较低的不足,减少了人为划分与对比地层带来的误差。
Description
技术领域
本发明属于油气地质勘探领域,特别涉及一种划分沉积地层的方法。
背景技术
钻测勘探作为油气地质的主要手段,日趋成熟。尤其是测井技术的快速发展,能够更加真实地反映地下沉积状况。而在构造活动活跃地区,由于隔夹层的频繁变化,加大了地层划分与对比的难度。
发明内容
本发明提供一种划分沉积地层的方法,运用地震地层学研究,将地震勘探原理中的反射系数计算方法应用于沉积地层的划分,进而进行地层的对比。
本发明通过如下技术方案实现:
一种划分沉积地层的方法,由测井曲线拟合地震反射系数曲线,根据地震反射系数曲线中不同位置反射系数的特征划分沉积地层。
根据地震反射系数曲线中不同位置反射系数的特征划分沉积地层的方法具体为界面法:
将反射系数偏离同侧层段反射系数平均值的幅度不小于平均值1/2处的地层确定为沉积层序的界面,相邻界面之间为一个层段,可将地层沿纵向划分成由多个界面相隔开的层段。
根据地震反射系数曲线中不同位置反射系数的特征划分沉积地层的另一种方法具体为层段法:
将反射系数偏离同侧层段反射系数平均值的幅度相似,同时偏离变化连续且均匀的一段作为同一层段,该层段的顶、底作为沉积层序的界面。其中,反射系数偏离同侧层段反射系数平均值的幅度相似是指偏离幅度均不小于同侧层段反射系数平均值的1/2,且相邻反射系数之间的偏离幅度小于该段反射系数平均值的1/2。
进一步的,测井曲线包括沉积物密度测井曲线DEN和穿透沉积物的声波时差测井曲线DT。
进一步的,沉积物密度测井曲线DEN通过中子密度测井、岩性密度测井或声波变密度测井计算获得。穿透沉积物的声波时差测井曲线DT通过声波时差测井或声波测井计算获得。
进一步的,由测井曲线拟合地震反射系数曲线的方法,包括如下步骤:
(1)计算沉积界面的沉积物密度测井曲线DEN;
(2)计算沉积界面的声波时差测井曲线DT;
(3)计算沉积界面上下岩层的波阻抗曲线Z,其中
(4)拟合地震反射系数曲线Rw,所述其中Rw1,2为沉积地层1和2之间界面的反射系数,Z1、Z2为沉积地层1和2的波阻抗,DEN1、DEN2分别为沉积地层1和2处的沉积物密度,DT1、DT2分别为沉积地层1和2处的声波时差测井值。
本发明运用地震地层学研究方法,利用测井曲线拟合地震反射系数曲线,进而通过沉积界面的反射系数特征划分沉积地层。
相比于传统的地层划分方案,本发明具有如下特点:
(1)地震反射系数反映了沉积界面处上下岩性的差异,较大的地震反射系数处代表了较大的不整合面,基于此来划分地层,体现了层序地层学划分地层的根本理论;
(2)由地震反射系数异常来划分地层,填补了因缺少区域标志层而无法划分与对比地层的科研现状;
(3)由于测井曲线纵向分辨率高,模拟出的地震反射系数异常可以很好的反映出沉积层岩性变化特征,可以弥补地震划分垂向层序分辨率较低的不足。
(4)地震反射系数为发育沉积地层本身的响应特征,减少了人为划分地层带来的误差。
附图说明
图1是采用本发明所述方法获得的测井曲线及拟合的地震反射系数曲线;
图2是不同井之间基于反射系数划分出沉积层段的横向对比成果图;
图3是通过计算机方法模拟反射系数曲线拟合过程的模拟效果图。
具体实施方式
下面对本发明作进一步详细说明。
本发明提供了一种划分沉积地层的方法。利用地震勘探原理中的反射系数计算方法,由测井曲线拟合地震反射系数曲线,根据地震反射系数曲线中不同位置反射系数的特征划分沉积地层。地震反射系数为发育沉积地层本身的响应特征,以反射系数为划分地层的判断依据,相比人为划分更真实、更准确。同时,由于测井曲线纵向分辨率高,拟合出的地震反射系数曲线的异常特征可以直观地反映出沉积地层岩性的变化特征,可以弥补地震划分垂向层序分辨率较低的不足。
根据地震反射系数曲线中不同位置反射系数的特征划分沉积地层的一种方法为界面法:
将反射系数偏离同层段系数平均值的幅度不小于平均值1/2处的地层确定为沉积层序的界面,相邻界面之间为一个层段,可将地层沿纵向划分成由多个界面相隔开的层段。在沉积层序界面处,由于上下岩性差异较大,由测井曲线计算形成的反射系数较为异常,根据反射系数的异常情况可将地层沿纵向划分成由多个界面相隔开的层段。
另一种方法为层段法:
将反射系数偏离同侧层段反射系数平均值的幅度相似,同时偏离变化连续且均匀的一段作为同一层段,该层段的顶、底作为沉积层序的界面。其中,反射系数偏离同侧层段反射系数平均值的幅度相似是指偏离幅度均不小于同侧层段反射系数平均值的1/2,且相邻反射系数之间的偏离幅度小于该段反射系数平均值的1/2。
上述两种方法可以单独划分沉积地层,亦可相互参照、配合使用。界面法是通过确定多个界面将整个地层纵向分割成多个层段;层段法是直接确定出异常层段,进而完成整个地层的分层。同一井下地层由于其复杂的岩性特征,通常需要同时参考上述两种方法进行地层的划分,有些层段适于采用层段法划分地层,有些地层需要采用界面法进行地层划分。实践中需要根据实际地质情况来选择。
上述用于拟合反射系数曲线的测井曲线包括沉积物密度测井曲线DEN和穿透沉积物的声波时差测井曲线DT。沉积物密度测井曲线DEN可通过中子密度测井、岩性密度测井或声波变密度测井计算获得。穿透沉积物的声波时差测井曲线DT可通过声波时差测井或声波测井计算获得。上述计算DEN和DT的方法是测井原理中的常规方法,本领域技术人员可以直接通过计算机计算获得。另外,DEN和DT也可通过相应的测井设备直接测量获取。
图1中的三条曲线分别是测井曲线及其拟合的地震反射系数曲线。如图1所示,自上而下分别为沉积物密度测井曲线DEN,声波波速曲线AC(AC=1/DT),和反射系数曲线Rw。其中,三条曲线沿横向自右向左依次划分为六个层段,分别是Sq1、Sq2、Sq3、Sq4、Sq5、Sq6;各层段之间自右向左依次为5个界面,分别是sb1,sb2,sb3,sb4,sb5。从图中可明显看出,Sq1层段与Sq2层段之间存在较大的反射系数异常值,sb1处的反射系数偏离同侧层段反射系数平均值的幅度明显大于Sq1层段反射系数平均值的一半,根据界面法,故将反射系数异常值对应的地层,即sb1作为两沉积段的沉积界面。另外,图1中Sq3层段内的反射系数变化较大,呈现出偏离同侧层段反射系数平均值的幅度均大于平均值的一半,且该段相邻反射系数之间的偏离幅度均小于该段反射系数平均值的1/2,同时该段反射系数的偏离变化连续而均匀。根据层段法,将该段反射系数异常部分对应的地层划分为一个沉积层段,即Sq3。
基于上述两种方法划分出的沉积层段,具有横向可对比性,即在联井剖面中,相应的沉积层段或沉积界面处,由于其具有相同或相似的反射系数特征,从而可用于联井剖面地层对比。
图2为不同井之间基于反射系数划分出沉积层段的横向对比成果图。从图中可以看出,同一沉积层段在不同井位上呈现出相同的反射系数特征,可以进行较好的沉积地层连井对比。由于沉积体的横向变化性,同一沉积层段在不同井位上有呈现出略微的不同,但不会影响沉积层段的横向对比。
另外,可通过如下步骤利用测井曲线进行地震反射系数曲线的拟合:
(1)计算沉积界面的沉积物密度测井曲线DEN;
(2)计算沉积界面的声波时差测井曲线DT;
(3)计算沉积界面上下岩层的波阻抗曲线Z,其中
(4)拟合地震反射系数曲线Rw,所述其中Rw1,2为沉积地层1和2之间界面的反射系数,Z1、Z2为沉积地层1和2的波阻抗,DEN1、DEN2分别为沉积地层1和2处的沉积物密度,DT1、DT2分别为沉积地层1和2处的声波时差测井值。
步骤(1)、(2)中的基础数据DEN和DT可由相关的测井仪器直接测量获得。选取常规测井数据中的沉积物密度测井曲线DEN与穿透沉积物的声波时差测井曲线DT为代表,以该两条测井曲线为基础数据。而声波时差DT的倒数为声波在介质中的传播速度AC,即AC=1/DT。故也可直接采用波速曲线AC作为基础数据。
在此基础上,根据托布里兹方程
可计算出相邻两点之间的反射系数。
其中,Z为介质的波阻抗,Rw1,2为地层1和2之间的反射系数。
表1即为拟合地震反射系数曲线Rw的方法,根据托布里兹方程其中a、b分别为沉积地层1和2处的沉积物密度,c、d分别为沉积界面1和2处的声波时差值。
表1.基础数据处理表
模拟计算过程
假设深度10-12米之间有五段沉积岩体,假设相同岩性层段对其测井响应值相同。结合测井原理,一米有8个测量点,所以得到的测井数据为18个离散数据(表2,DEN和DT列)。在实际研究过程中,这些值均由测井仪器获得。
表2沉积层段模拟计算数值表
由于测井数据是将沉积层沿纵向每米划分为八个参考点,即相邻两个参考点间的距离为0.125米,计算出的反射系数为相邻两个参考点间的数值。为便于作图,特将其上提至参考点位置。这样就可以将反射系数形成的曲线与测井曲线投射在相同的坐标系内,如图1所示。以上操作均可由计算机单独完成。
步骤一,计算波阻抗。根据公式可以计算出每个沉积层段测量点的波阻抗值,得到表2中的Z列;
步骤二,计算反射系数。根据公式可以计算出两个测量点间的反射系数值。为不增加测井数据的深度值,特将反射系数值向上提高了0.125/2米。
同一沉积段内,反射系数为0,反射系数上提不会受影响;界面处的反射系数由于提高而使得测量的界面向上偏差了0.125/2米,在一定程度上可以忽略。
步骤三,划分沉积地层。可以看出,岩性均一的沉积段内反射系数为0,而界面处反射系数不为0,即下伏地层波阻抗较大时,界面处反射系数为正值;下伏地层波阻抗值较小时,界面处反射系数值为负值;界面上下岩性差异越大,反射系数的绝对值越大。反射系数不为0处为沉积界面,反射系数为0或接近0段为沉积层段。基于此划分出5个层段,如表2所示。
图3是上述模拟计算中反射系数曲线拟合过程的模拟效果图。从图3中可看出,层段一至层段三均可根据界面法进行划分,三个层段之间的界面其反射系数均明显偏离同侧层段反射系数平均值的一半。而层段四可根据层段法进行划分,该层段内的反射系数均偏离同侧层段反射系数平均值的一半,且相邻反射系数之间的偏离幅度也小于该段反射系数平均值的1/2,同时偏离变化连续且均匀。
本发明基于沉积界面的反射系数特征,对沉积层序进行划分,为地层的对比提供技术支撑,提高现场技术人员工作效率。
本实施例没有详细叙述的部件、结构及工艺和计算方法属本行业的公知部件和常用结构或常用手段,这里不一一叙述。
Claims (3)
1.一种划分沉积地层的方法,其特征在于,由测井曲线拟合地震反射系数曲线,根据地震反射系数曲线中不同位置反射系数的特征划分沉积地层;
所述的根据地震反射系数曲线中不同位置反射系数的特征划分沉积地层的方法具体为界面法:
将反射系数偏离同侧反射系数平均值的幅度不小于平均值1/2处的地层确定为沉积层序的界面,相邻界面之间为一个层段,将地层沿纵向划分成由多个界面相隔开的层段;
所述的根据地震反射系数曲线中不同位置反射系数的特征划分沉积地层的方法具体为层段法:
将反射系数偏离同侧层段反射系数平均值的幅度相似,同时偏离变化连续且均匀的一段作为同一层段,该层段的顶、底作为沉积层序的界面;
所述反射系数偏离同侧层段反射系数平均值的幅度相似是指偏离幅度均不小于同侧层段反射系数平均值的1/2,且相邻反射系数之间的偏离幅度小于该段反射系数平均值的1/2;
所述测井曲线包括沉积物密度测井曲线DEN和穿透沉积物的声波时差测井曲线DT。
2.根据权利要求1所述的划分沉积地层的方法,其特征在于,所述沉积物密度测井曲线DEN通过中子密度测井、岩性密度测井或声波变密度测井计算获得,所述穿透沉积物的声波时差测井曲线DT通过声波时差测井或声波测井计算获得。
3.根据权利要求1所述的划分沉积地层的方法,其特征在于,所述由测井曲线拟合地震反射系数曲线的方法,包括如下步骤:
(1)计算沉积界面的沉积物密度测井曲线DEN;
(2)计算沉积界面的声波时差测井曲线DT;
(3)计算沉积界面上下岩层的波阻抗曲线Z,其中
(4)拟合地震反射系数曲线Rw,其中,Rw1,2为界面1和2之间的反射系数,DEN1、DEN2分别为沉积界面1和2处的沉积物密度,DT1、DT2分别为沉积界面1和2处的声波时差测井值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410209502.5A CN104036119B (zh) | 2014-05-16 | 2014-05-16 | 划分沉积地层的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410209502.5A CN104036119B (zh) | 2014-05-16 | 2014-05-16 | 划分沉积地层的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104036119A CN104036119A (zh) | 2014-09-10 |
CN104036119B true CN104036119B (zh) | 2017-05-03 |
Family
ID=51466889
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410209502.5A Expired - Fee Related CN104036119B (zh) | 2014-05-16 | 2014-05-16 | 划分沉积地层的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104036119B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105572726B (zh) * | 2014-10-13 | 2019-02-01 | 成都北方石油勘探开发技术有限公司 | 一种层序地层格架下多期次叠置扇体的精细刻画方法 |
CN108732639B (zh) * | 2017-04-17 | 2020-12-11 | 中国石油化工股份有限公司 | 一种页岩地层等时对比方法及*** |
CN110320570A (zh) * | 2018-03-30 | 2019-10-11 | 中国石油化工股份有限公司 | 一种复杂地质条件下的综合地层对比方法 |
CN111594101B (zh) * | 2020-03-19 | 2021-03-12 | 中国石油大学(华东) | 一种垂直井管外出砂亏空剖面预测及高饱和砾石充填施工方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1114511A (zh) * | 1993-09-21 | 1996-01-03 | 法国石油研究所 | 采用统计校准技术导出地质特性的地震记录曲线分析方法 |
CN101013161A (zh) * | 2007-01-15 | 2007-08-08 | 中国石油大港油田勘探开发研究院 | 基于叠前波场模拟的地震勘探层位标定方法 |
CN101025084A (zh) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | 一种随钻预测钻头底下地层孔隙压力的方法 |
EP1312027B1 (en) * | 2000-06-20 | 2008-01-02 | Baker Hughes Incorporated | Method for enhancing resolution of elastic wave velocities by isolating a wave event in lithographic formation |
CN102681015A (zh) * | 2012-05-30 | 2012-09-19 | 中国地质大学(北京) | 矿区地层结构划分方法 |
-
2014
- 2014-05-16 CN CN201410209502.5A patent/CN104036119B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1114511A (zh) * | 1993-09-21 | 1996-01-03 | 法国石油研究所 | 采用统计校准技术导出地质特性的地震记录曲线分析方法 |
EP1312027B1 (en) * | 2000-06-20 | 2008-01-02 | Baker Hughes Incorporated | Method for enhancing resolution of elastic wave velocities by isolating a wave event in lithographic formation |
CN101025084A (zh) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | 一种随钻预测钻头底下地层孔隙压力的方法 |
CN101013161A (zh) * | 2007-01-15 | 2007-08-08 | 中国石油大港油田勘探开发研究院 | 基于叠前波场模拟的地震勘探层位标定方法 |
CN102681015A (zh) * | 2012-05-30 | 2012-09-19 | 中国地质大学(北京) | 矿区地层结构划分方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104036119A (zh) | 2014-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Corbeanu et al. | Detailed internal architecture of a fluvial channel sandstone determined from outcrop, cores, and 3-D ground-penetrating radar: Example from the middle Cretaceous Ferron Sandstone, east-central Utah | |
CN103256046B (zh) | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 | |
CN104950331B (zh) | 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法 | |
US9354338B1 (en) | Generating synthetic seismic traces | |
CN107783187B (zh) | 一种将测井速度和地震速度结合建立三维速度场的方法 | |
CN108680951A (zh) | 一种基于地震信息判断煤层气富集沉积控制作用的方法 | |
US20120053841A1 (en) | Well constrained horizontal variable h-v curve constructing method for seismic wave velocity field construction | |
CN104898161B (zh) | 一种基于测井响应模拟体的有效砂岩预测方法 | |
CN102636812B (zh) | 一种获得碳酸盐岩储层储集空间体积的方法 | |
CN105651676A (zh) | 一种水平井规则开发井网下的储层非均质性表征方法 | |
CN105549082B (zh) | 超深碳酸盐储层三维地质力学场建立方法和*** | |
CN109884710B (zh) | 针对激发井深设计的微测井层析成像方法 | |
US20150106018A1 (en) | 3D-Well Log Invention | |
CN104036119B (zh) | 划分沉积地层的方法 | |
Milad et al. | Upscaled DFN models to understand the effects of natural fracture properties on fluid flow in the Hunton Group tight Limestone | |
CN106285642B (zh) | 一种基于地震资料的地层孔隙压力预测方法 | |
CN101315428B (zh) | 基于横波速度的高孔隙碎屑岩的高分辨率勘探方法 | |
Schneider et al. | Interpretation of fractured zones using seismic attributes—Case study from Teapot Dome, Wyoming, USA | |
Corbeanu et al. | Prediction of 3-D fluid permeability and mudstone distributions from ground-penetrating radar (GPR) attributes: Example from the Cretaceous Ferron Sandstone Member, east-central Utah | |
Yusof et al. | Uncertainty analysis of hydrocarbon in place calculation using 3D seismic and well data during appraisal stage–Case study of Goldie Field, offshore Sarawak | |
Roy et al. | Integrated characterization of hydraulic fracture treatments in the Barnett Shale: The Stocker geophysical experiment | |
Gao et al. | Seismic structure and texture analyses for fractured reservoir characterization: An integrated workflow | |
CN115857047B (zh) | 一种地震储层综合预测方法 | |
Mehdipour et al. | The Best Scenario for Geostatistical Modeling of Porosity in the Sarvak Reservoir in an Iranian Oil Field, Using Electrofacies, Seismic Facies, and Seismic Attributes | |
CN109459790A (zh) | 针对煤系地层地震速度场建立方法及*** |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170503 |
|
CF01 | Termination of patent right due to non-payment of annual fee |