CN108333627A - 火成岩地区断裂真伪的识别方法及装置 - Google Patents
火成岩地区断裂真伪的识别方法及装置 Download PDFInfo
- Publication number
- CN108333627A CN108333627A CN201810040119.XA CN201810040119A CN108333627A CN 108333627 A CN108333627 A CN 108333627A CN 201810040119 A CN201810040119 A CN 201810040119A CN 108333627 A CN108333627 A CN 108333627A
- Authority
- CN
- China
- Prior art keywords
- igneous rock
- section
- seismic
- data
- migration
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 168
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000013508 migration Methods 0.000 claims abstract description 115
- 230000005012 migration Effects 0.000 claims abstract description 115
- 238000003384 imaging method Methods 0.000 claims abstract description 36
- 238000005553 drilling Methods 0.000 claims description 42
- 238000012937 correction Methods 0.000 claims description 10
- 241001074085 Scophthalmus aquosus Species 0.000 claims description 2
- 238000004590 computer program Methods 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 abstract description 3
- 238000010276 construction Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 239000002245 particle Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 125000005587 carbonate group Chemical group 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种火成岩地区断裂真伪的识别方法及装置，包括步骤：利用野外地震数据进行地震偏移成像得到地震偏移剖面；从地震深度偏移剖面上寻找需要判断真伪的断裂及对应的火成岩厚度；同时利用深度偏移速度场获取需要判断真伪断裂对应的火成岩偏移速度变化范围，通过测井岩性解释统计准确的火成岩速度；然后计算出不同厚度火成岩偏移速度误差引起的垂直位移阈值曲线；从地震深度偏移剖面上获取断裂处的垂向断距H，同时计算对应火成岩厚度的断距误差阈值如果H>Hxmax，确定为真断裂，否则为假断裂。本发明方法及装置可以较为准确地识别出火成岩地区断裂的真伪。
Description
技术领域
本发明涉及油气勘探技术领域，特别涉及一种火成岩地区断裂真伪的识别方法及装置。
背景技术
塔里木盆地顺拖果勒地区储集体类型以断控-裂缝孔洞型为主，其主要特征是：①埋藏深度大于7000米；②储集体岩性为碳酸盐岩，非均质性强；③地表为沙漠覆盖区；④二叠系广泛发育高速火成岩。断裂在断控-裂缝孔洞系统中为主控因素，也是钻探的主要目标，在该地区油气勘探开发中具有举足轻重的作用。
现阶段储集体预测及断裂识别主要依赖地震偏移成像资料，地震偏移成像的精度直接影响地质目标的研究。地震资料成像分为时间域偏移成像和深度域偏移成像。时间域偏移成像具有速度建模简单、计算速度快等优点，同时也有高速异常体引起地层构造假象的缺点。深度域偏移成像具有速度建模复杂、计算慢等特点，理论上可以消除高速异常体对地震偏移成像的影响，恢复真实地层构造形态。在地震偏移成像实际应用中，由于工区内钻井资料少，地层中存在具有一定厚度的高速火成岩体，深度偏移成像前高速火成岩体的平面分布范围、厚度及速度皆不能准确描述；因此深度域速度建模受高速火成岩的影响较大，致使偏移速度场的准确性降低，造成深度域偏移成像资料也不能完全消除高速火成岩引起的构造假象。
由于高速火成岩的影响，地震偏移资料或多或少存在构造假象，现有地震偏移资料上用常规方法识别断裂就会存在缺陷，必然造成识别出来的断裂存在构造假象因素，特别是偏移速度误差较大时产生构造假象的程度越大。
发明内容
本发明的目的在于改善现有技术中所存在的常规方法无法准确识别断裂的不足，提供一种火成岩地区断裂真伪的识别方法及装置，提高断裂识别的准确度。
为了实现上述发明目的，本发明实施例提供了以下技术方案：
一种火成岩地区断裂真伪的识别方法，包括以下步骤：
获取待测工区真实的野外地震数据，利用野外地震数据进行地震偏移成像，得到地震深度偏移剖面；
从地震深度偏移剖面上获取该待测工区火成岩厚度数据，利用待测工区的实钻井测井数据中的火成岩厚度，对从地震深度偏移剖面上获取的火成岩厚度数据进行校正；
从地震深度偏移剖面上获取断裂处顶部火成岩段的偏移速度范围，以实钻井火成岩平均声波速度为基准，计算出火成岩偏移速度的误差；所述实钻井火成岩平均声波速度由实钻井测井数据中的声波时差计算得到；
利用校正后的火成岩厚度数据和火成岩偏移速度的误差，在预建立的垂向断距阈值表中，获取偏移速度误差范围引起的断裂失真的垂向断距阈值Hx，Hx∈[0,Hxmax]，Hxmax为偏移速度误差对应断距最大值；
从地震深度偏移剖面上获取断裂处的垂向断距H，比较Hx与H的大小，若H>Hxmax，则识别为该处断裂为真实断裂，否则识别为由于火成岩偏移速度偏差引起的假性断裂。
进一步地，上述方法中，所述预建立的垂向断距阈值表通过以下步骤建立：
从实钻井测井数据中获得各地层介质、各地层火成岩厚度及各地层火成岩段的声波时差数据，根据声波时差计算声波速度；
根据各地层介质、各地层火成岩厚度及各地层火成岩段的声波速度建立理论地震地质模型，并对理论地震地质模型进行叠前波动方程正演模拟，得到单炮模拟记录；
对所述单炮模拟记录进行深度域偏移成像，且更改偏移速度，分别得到不同火成岩厚度、不同偏移速度下的深度域地震偏移剖面；
在深度域地震偏移剖面上拾取火成岩下伏水平地层的垂向断距，建立不同火成岩厚度、不同偏移速度误差对应的垂向断距阈值表。
本发明同时提供了一种火成岩地区断裂真伪的识别装置，包括以下模块：
地震深度偏移剖面建立模块，用于获取待测工区真实的野外地震数据，利用野外地震数据进行地震偏移成像，得到地震深度偏移剖面；
火成岩厚度校正模块，用于从地震深度偏移剖面上获取该待测工区火成岩厚度数据，利用待测工区的实钻井测井数据中的火成岩厚度，对从地震深度偏移剖面上获取的火成岩厚度数据进行校正；
偏移速度误差计算模块，用于从地震深度偏移剖面上获取断裂处顶部火成岩段的偏移速度范围，以实钻井火成岩平均声波速度为基准，计算出火成岩偏移速度的误差；所述实钻井火成岩平均声波速度由实钻井测井数据中的声波时差计算得到；
垂向断距Hx获取模块，用于利用校正后的火成岩厚度数据和火成岩偏移速度的误差，在预建立的垂向断距阈值表中，获取偏移速度误差范围引起的断裂失真的垂向断距阈值Hx，Hx∈[0,Hxmax]，Hxmax为偏移速度误差对应断距最大值；
断裂真伪识别模块，用于从地震偏移剖面上获取断裂处的垂向断距H，比较Hx与H的大小，若H>Hxmax，则识别为该处断裂为真实断裂，否则识别为由于火成岩偏移速度偏差引起的假性断裂。
本发明实施例还提供了一种包括计算机可读指令的计算机可读存储介质，所述计算机可读指令在被执行时使处理器执行所述程序时实现任一火成岩地区断裂真伪的识别方法的步骤。
本发明实施例还提供了一种电子设备，所述的电子设备包括处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现任一火成岩地区断裂真伪的识别方法的步骤。
与现有技术相比，本发明的有益效果：在火成岩发育地区，下伏断裂成像受高速火成岩影响出现假构造时，现有的断裂识别方法精度不高或识别出假断裂。本发明的方法是在常规断裂识别方法的基础上，利用实钻井信息、模型正演技术和定量分析手段，克服深度域地震成像偏移速度不准产生假构造的问题，从而可以有效识别火成岩地区断裂的真伪，为布井提供强有力的依据，减小油气勘探的风险。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为本发明实施提供的火成岩地区断裂真伪识别方法的流程示意图。
图2a-f分别为火成岩地震地质模型及不同偏移速度误差下的深度偏移剖面图。
图3a-b分别为火成岩厚度、偏移速度误差与垂向断距的关系表和示意图。
图4是本发明火成岩地区断裂真伪识别装置的功能模块示意图。
图5为本实施例中垂向断距阈值表建立模块的功能模块示意图。
图6为本发明实施例中提供的电子设备的结构示意图。
图7a和图7b分别为示例中同一原始地震数据经过PSDM及RTM得到的目的层地震数据体平面断裂展布特征相干属性图。
图8为示例中同一原始地震数据得到的地震深度偏移剖面图。
具体实施方式
下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅图1，本实施例中提供的火成岩地区断裂真伪的识别方法，包括以下步骤：
步骤S1：统计工区内钻井信息，工区内钻井信息至少包括岩性录井资料和声波测井曲线，这些钻井信息是在实际钻井时测量所得。从岩性录井资料上获取(岩性录井资料会记录各地层介质及火成岩厚度数据)实钻井各地层及火成岩厚度；从声波测井曲线上获取各地层及火成岩段的声波时差，并把声波时差转换为声波速度，声波时差与声波速度的关系式为：
式中：V—声波速度(单位：m/s)；
AC—声波时差(单位：us/ft)；
并求取多口井各地层及火成岩声波速度的平均值，关系式为：
式中：Va—平均声波速度(单位：m/s)；
Vj—第j个采样点的声波速度(单位：m/s)；
n—钻遇火成岩的钻井数量；
mi—第i口井火成岩段采样点的数量，一个深度段对应一个采样点。
步骤S2：根据步骤S1中获得的各地层介质、各地层火成岩厚度及各地层火成岩段的声波速度建立理论地震地质模型，并对理论地震地质模型进行叠前波动方程正演模拟，得到单炮模拟记录数据。
具体地，本步骤通过以下两个小步骤执行完成：
1)根据岩性录井资料上获取的各地层数据，在模型设计软件中描述地层格架，再根据各地层界面深度、步骤S1中计算出的火成岩段声波速度，对地层格架进行填充，建立不同厚度火成岩的理论地震地质模型，理论地震地质模型如图2中图2a所示，图2a中左侧的数据从上至下依次为3500/4000/4500/5000/5500/6000/6500/7000，此数据为模型深度，即地层格架深度；该理论地震地质模型的特征为：a)模型中填充均匀介质，单一地层单元中介质的地震波速度为恒定值；b)模型中火成岩厚度为20m、40m、60m、80m、100m、120m、140m、160m、180m、200m，火成岩宽度都为200米。本实施例中根据工区内钻井信息建立的地震地质模型中火成岩的厚度、宽度数据为上述数据，但是实际上建立的地震地震模型中火成岩厚度、宽度数据可以根据实际情况调整，即根据具体工区的具体钻井信息调整；c)火成岩下伏地层为水平地层。
2)对该理论地震地质模型设计地震数据采集观测系统，进行波动方程正演模拟，得到单炮模拟记录。
具体地，考虑到火成岩埋深及实际地震资料主频，设计地震数据采集观测系统如下：
| 观测形式：单边激发 | 道距：25米 |
| 排列长度：4000米 | 震源：雷克子波(30Hz) |
| 炮距：25米 | 计算网格：1米 |
需要说明的是，上述地震数据采集观测系统仅是作为一种举例，上述地震数据采集观测系统中数据可根据实际情况调整。
然后进行波动方程正演，波动方程正演公式如下：
式中，v—为空间质点(x，z)处地震纵波速度，单位(m/s)；
μ—为空间质点受外力后产生的位移，单位(m)；
x—空间质点的横坐标，单位(m)；
z—空间质点的纵坐标，单位(m)；
t-震源激发开始后的时间，单位(s)。
步骤S3：用上述理论地震地质模型的地震速度建立准确的深度域偏移速度场，因为理论地震地震模型中的声波速度是由实钻井时声波测井曲线上获取声波时差转换而来，所以此处的偏移速度是准确的，即理论地震地质模型中包含的介质速度，对偏移成像是准确的。对上述正演模拟记录进行深度域偏移成像，得到偏移速度准确情况下的深度域地震偏移剖面。所述深度域偏移成像是把绕射波收敛后形成深度域地震剖面，例如可以通过Focus、Promax、CGG、Omega等等，输入单炮地震数据、地震观测系统数据、测量数据等即可获得深度域地震偏移剖面。偏移速度场中分别对火成岩速度减小5％、10％、15％、20％，用不同偏移速度进行深度域偏移成像，分别得到该偏移速度的地震成像剖面，如图2中图2b-2f所示。分别在不同火成岩厚度、不同偏移速度情况下，在深度域地震偏移剖面上拾取火成岩下伏水平地层的垂向断距数据，即测量深度域地震偏移剖面上同相轴错断的垂向距离。
步骤S4：采用步骤S3中获得的不同火成岩厚度、不同偏移速度误差和垂向断距(这些数据可以在深度偏移剖面上直接读取)，拟合理论曲线，即以火成岩厚度、垂向断距为坐标轴，分别建立不同偏移速度误差的散点图，拟合散点的关系式，建立不同火成岩厚度、不同偏移速度误差对应的垂向断距阈值表，如图3a和图3b所示，图3a中从上到下依次为速度误差20％、15％、10％、5％的曲线。
步骤S5：获取待测工区真实的野外地震数据(例如单炮地震记录，激发点、接收点的大地坐标信息、海拔高程信息等，微测井、小折射测量的静校正信息等)，并进行深度域地震偏移成像，得到地震深度偏移剖面，从地震深度偏移剖面上获取该待测工区火成岩厚度，形成火成岩厚度分布图；并用实钻井火成岩厚度对该火成岩厚度分布图中的火成岩厚度进行校正，得到校正后的火成岩厚度平面分布图。
可参考图2b-f，图2b-f是理论模型正演及偏移后深度域地震剖面，左边数字纵向坐标轴，火成岩下伏地层的反射轴有水平的(速度准确)，有错断的(速度不准确)，通过反射同相轴错断的距离即垂向断距，可以直接从垂向坐标轴上读取出来或者在地层解释软件中拾取出来。
工区火成岩厚度：先是通过火成岩井震标定，在地震偏移剖面上确定其反射特征，连续追踪火成岩顶底界面，即得到火成岩顶界面、底界面的平面分布范围，顶、底界面之差即火成岩的厚度。
速度数据：地震偏移成像就需要速度，实际野外地震数据偏移，速度通过速度分析、速度建模技术得到，因此有偏移速度体，这个速度只是室内技术手段得到，不一定能准确描述野外地层的速度，因此对于偏移成像来说存在误差。
需要说明的是，此处真实的野外地震数据，经过室内地震偏移成像的剖面，剖面表现的是真实地层的地球物理响应特征。步骤S1中的实钻井钻遇数据，即通过野外钻井、录井、测井获取地层信息，包括地层的岩性、地层厚度、岩石速度等。
步骤S6：在实际地震偏移速度场中获取断裂处火成岩段的偏移速度范围，把实钻井火成岩平均声波速度(步骤S1中求出的Va)作为基准，计算出火成岩偏移速度的误差。实际地震偏移速度场由真实的野外地震数据得到，在地震偏移剖面上可以读取火成岩段的偏移速度值。
步骤S7：利用步骤S5中校正后的火成岩厚度和步骤S6中火成岩偏移速度的误差，在步骤S4中建立的垂向断距阈值表中，获取偏移速度误差引起的断裂失真的垂向断距Hx，Hx∈[0,Hxmax]，Hxmax为偏移速度误差对应垂向断距最大值，无偏移速度误差时垂向断距为0。
步骤S8：在实际地震深度偏移剖面上获取断裂的垂向断距H(实际值)，该垂向断距与步骤S7中断裂失真的垂向断距Hx(理论值)进行比较，判别断裂的真伪。具体地：
1)如果H在误差阈值范围[0,Hxmax]内，则认为该断裂是构造假象，是由于火成岩偏移速度偏差引起的假性断裂；
2)如果H超出误差阈值范围Hx最大值，即H>Hxmax，则认为该断裂是真实断裂。
请参阅图7a、图7b和图8，针对于同一套原始地震采集数据体，先后两次处理的PSDM及RTM分别得到的目的层地震数据体平面断裂展布特征相干属性图如图7a、图7b所示，得到的深度域地震剖面如图8所示。PSDM及RTM两套数据体平面断裂展布特征均可明显识别断裂出4、6、7，而断裂1、2、3、5出现在PSDM数据体中，在RTM数据体中并未显示，即是说RTM不能识别出1、2、3、5处可能存在断裂，而PSDM则将1、2、3、5处全部识别为断裂。而采用本发明方法通过计算得出偏移速度误差范围为-8％～13.6％，在计算Hx采用的偏移速度误差范围取绝对值范围。可以得出2号断裂为真实存在的断裂，1、3、5处为假性断裂，由此可见，本发明方法可以更加准确地识别出断裂的真伪。
需要说明的是，步骤S1-S4是在建立阈值表，步骤S5-S8是断裂真伪的判别过程。对同一个工区而言，步骤S1-S4只需执行一次，判别多条断裂的真伪时，需要多次执行步骤S6-S8。对下一个工区而言，则需要执行步骤S1-S8。
本发明方法，利用野外地震数据进行地震偏移成像得到地震深度偏移剖面；从地震深度偏移剖面上寻找需要判断真伪的断裂及对应的火成岩厚度，同时利用深度偏移速度场获取需要判断真伪断裂对应的火成岩偏移速度变化范围，通过测井岩性解释统计准确的火成岩速度，然后计算出不同厚度火成岩偏移速度误差引起的垂直位移阈值曲线，从地震深度偏移剖面上获取断裂处的垂向断距H，同时计算对应火成岩厚度的断距误差阈值HXmax，如果H>Hxmax，确定为真断裂，否则为假断裂，可以更为准确地识别出火成岩地区断裂的真伪性。
请参阅图4-5，本实施例中同时提供了一种火成岩地区断裂真伪的识别装置，包括地震偏移剖面建立模块，火成岩厚度校正模块，偏移速度误差计算模块，垂向断距Hx获取模块。其中，
地震深度偏移剖面建立模块，用于获取待测工区真实的野外地震数据，利用野外地震数据进行地震偏移成像，得到地震深度偏移剖面；
火成岩厚度校正模块，用于从地震深度偏移剖面上获取该待测工区火成岩厚度数据，利用待测工区的实钻井测井数据中的火成岩厚度，对从地震深度偏移剖面上获取的火成岩厚度数据进行校正；
偏移速度误差计算模块，用于从地震深度偏移剖面上获取断裂处顶部火成岩段的偏移速度范围，以实钻井火成岩平均声波速度为基准，计算出火成岩偏移速度的误差；所述实钻井火成岩平均声波速度由实钻井测井数据中的声波时差计算得到；
垂向断距Hx获取模块，用于利用校正后的火成岩厚度数据和火成岩偏移速度的误差，在预建立的垂向断距阈值表中，获取偏移速度误差范围引起的断裂失真的垂向断距阈值Hx，Hx∈[0,Hxmax]，Hxmax为偏移速度误差对应断距最大值；
断裂真伪识别模块，用于从地震偏移剖面上获取断裂处的垂向断距H，比较Hx与H的大小，若H>Hxmax，则识别为该处断裂为真实断裂，否则识别为由于火成岩偏移速度偏差引起的假性断裂。
上述火成岩地区断裂真伪的识别装置还包括垂向断距阈值表建立模块，用于建立不同火成岩厚度、不同偏移速度误差对应的垂向断距阈值表。具体的，所述垂向断距阈值表建立模块包括：
实钻井测井数据获取模块，用于从实钻井测井数据中获得各地层介质、各地层火成岩厚度及各地层火成岩段的声波时差数据，根据声波时差计算声波速度；
理论地震地质模型建立模块，用于根据各地层介质、各地层火成岩厚度及各地层火成岩段的声波速度建立理论地震地质模型，并对理论地震地质模型进行波动方程正演模拟，得到单炮模拟记录；
偏移成像模块，用于对所述单炮模拟记录进行深度域偏移成像，且更改偏移速度，分别得到不同火成岩厚度、不同偏移速度下的深度域地震偏移剖面；本偏移成像模块可以直接采用现有的软件实现。
阈值表建立模块，用于在深度域地震偏移剖面上拾取火成岩下伏水平地层的垂向断距，建立不同火成岩厚度、不同偏移速度误差对应的垂向断距阈值表。
上述火成岩地区断裂真伪的识别装置是基于上述火成岩地区断裂真伪的识别方法的同一发明构思而进行的，因此为详细描述之处请参见前述方法实施例中的相应描述。
如图6所示，本实施例同时提供了一种电子设备，该电子设备可以包括处理器51和存储器52，其中存储器52耦合至处理器51。值得注意的是，该图是示例性的，还可以使用其他类型的结构来补充或替代该结构，实现数据提取、图表重绘、通信或其他功能。
如图6所示，该电子设备还可以包括：输入单元53、显示单元54和电源55。值得注意的是，该电子设备也并不是必须要包括图6中显示的所有部件。此外，电子设备还可以包括图6中没有示出的部件，可以参考现有技术。
处理器51有时也称控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该处理器51接收输入并控制电子设备的各个部件的操作。
其中，存储器52例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其他合适装置中的一种或多种，可存储上述处理器51的配置信息、处理器51执行的指令、记录的图表数据等信息。处理器51可以执行存储器52存储的程序，以实现信息存储或处理等。在一个实施例中，存储器52中还包括缓冲存储器，即缓冲器，以存储中间信息。
输入单元53例如可以为文件读取装置，用于向处理器51提供包含图表的电子文档。显示单元54用于显示从电子文档中提取出的图表信息，以及显示根据图表信息重绘得到的图表，该显示单元例如可以为LCD显示器，但本发明并不限于此。电源55用于为电子设备提供电力。
本发明实施例还提供一种计算机可读指令，其中当在电子设备中执行所述指令时，所述程序使得电子设备执行如图1所示的火成岩地区断裂真伪的识别方法所包含的操作步骤，或者如图1所示方法中的一部分步骤，例如步骤S1-S4，或者步骤S5-S8。
本发明实施例还提供一种存储有计算机可读指令的存储介质，其中所述计算机可读指令使得电子设备执行如图1所示的火成岩地区断裂真伪的识别方法所包含的操作步骤，或者如图1所示方法中的一部分步骤，例如步骤S1-S4，或者步骤S5-S8。
应理解，在本发明的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口、装置或单元的间接耦合或通信连接，也可以是电的，机械的或其它的形式连接。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本发明实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
Claims (8)
1.一种火成岩地区断裂真伪的识别方法，其特征在于，包括以下步骤：
获取待测工区真实的野外地震数据，利用野外地震数据进行地震偏移成像，得到地震深度偏移剖面；
从地震深度偏移剖面上获取该待测工区火成岩厚度数据，利用待测工区的实钻井测井数据中的火成岩厚度，对从地震深度偏移剖面上获取的火成岩厚度数据进行校正；
从地震深度偏移剖面上获取断裂处顶部火成岩段的偏移速度范围，以实钻井火成岩平均声波速度为基准，计算出火成岩偏移速度的误差；所述实钻井火成岩平均声波速度由实钻井测井数据中的声波时差计算得到；
利用校正后的火成岩厚度数据和火成岩偏移速度的误差，在预建立的垂向断距阈值表中，获取偏移速度误差范围引起的断裂失真的垂向断距阈值Hx，Hx∈[0,Hxmax]，Hxmax为偏移速度误差对应断距最大值；
从地震深度偏移剖面上获取断裂处的垂向断距H，比较Hx与H的大小，若H>Hxmax，则识别为该处断裂为真实断裂，否则识别为由于火成岩偏移速度偏差引起的假性断裂。
2.根据权利要求1所述的方法，其特征在于，所述预建立的垂向断距阈值表通过以下步骤建立：
从实钻井测井数据中获得各地层介质、各地层火成岩厚度及各地层火成岩段的声波时差数据，根据声波时差计算声波速度；
根据各地层介质、各地层火成岩厚度及各地层火成岩段的声波速度建立理论地震地质模型，并对理论地震地质模型进行叠前波动方程正演模拟，得到单炮模拟记录；
对所述单炮模拟记录进行深度域偏移成像，且更改偏移速度，分别得到不同火成岩厚度、不同偏移速度下的深度域地震偏移剖面；
在深度域地震偏移剖面上拾取火成岩下伏水平地层的垂向断距，建立不同火成岩厚度、不同偏移速度误差对应的垂向断距阈值表。
3.根据权利要求2所述的方法，其特征在于，所述实钻井火成岩平均声波速度通过以下步骤计算得到：
获取实钻井井测数据中的声波时差，将声波时差通过公式转换为声波速度，V为声波速度，AC为声波时差；
求取待测工区多口井各地层火成岩平均声波速度
式中，Va是平均声波速度；Vj是第j个采样点的声波速度；n是钻井数量，mi是第i口井火成岩段采样点的数量。
4.一种火成岩地区断裂真伪的识别装置，其特征在于，包括以下模块：
地震深度偏移剖面建立模块，用于获取待测工区真实的野外地震数据，利用野外地震数据进行地震偏移成像，得到地震深度偏移剖面；
火成岩厚度校正模块，用于从地震深度偏移剖面上获取该待测工区火成岩厚度数据，利用待测工区的实钻井测井数据中的火成岩厚度，对从地震深度偏移剖面上获取的火成岩厚度数据进行校正；
偏移速度误差计算模块，用于从地震深度偏移剖面上获取断裂处顶部火成岩段的偏移速度范围，以实钻井火成岩平均声波速度为基准，计算出火成岩偏移速度的误差；所述实钻井火成岩平均声波速度由实钻井测井数据中的声波时差计算得到；
垂向断距Hx获取模块，用于利用校正后的火成岩厚度数据和火成岩偏移速度的误差，在预建立的垂向断距阈值表中，获取偏移速度误差范围引起的断裂失真的垂向断距阈值Hx，Hx∈[0,Hxmax]，Hxmax为偏移速度误差对应断距最大值；
断裂真伪识别模块，用于从地震偏移剖面上获取断裂处的垂向断距H，比较Hx与H的大小，若H>Hxmax，则识别为该处断裂为真实断裂，否则识别为由于火成岩偏移速度偏差引起的假性断裂。
5.根据权利要求4所述的装置，其特征在于，还包括垂向断距阈值表建立模块，用于建立不同火成岩厚度、不同偏移速度误差对应的垂向断距阈值表；所述垂向断距阈值表建立模块包括：
实钻井测井数据获取模块，用于从实钻井测井数据中获得各地层介质、各地层火成岩厚度及各地层火成岩段的声波时差数据，根据声波时差计算声波速度；
理论地震地质模型建立模块，用于根据各地层介质、各地层火成岩厚度及各地层火成岩段的声波速度建立理论地震地质模型，并对理论地震地质模型进行叠前波动方程正演模拟，得到单炮模拟记录；
偏移成像模块，用于对所述单炮模拟记录进行深度域偏移成像，且更改偏移速度，分别得到不同火成岩厚度、不同偏移速度下的深度域地震偏移剖面；
阈值表建立模块，用于在深度域地震偏移剖面上拾取火成岩下伏水平地层的垂向断距，建立不同火成岩厚度、不同偏移速度误差对应的垂向断距阈值表。
6.根据权利要求5所述的装置，其特征在于，还包括平均声波速度计算模块，用于获取实钻井井测数据中的声波时差，将声波时差通过公式转换为声波速度，V为声波速度，AC为声波时差，并求取待测工区多口井各地层火成岩平均声波速度Va是平均声波速度；Vj是第j个采样点的声波速度；n是钻井数量，mi是第i口井火成岩段采样点的数量。
7.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1-3任一所述方法的步骤。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1-3任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810040119.XA CN108333627B (zh) | 2018-01-16 | 2018-01-16 | 火成岩地区断裂真伪的识别方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810040119.XA CN108333627B (zh) | 2018-01-16 | 2018-01-16 | 火成岩地区断裂真伪的识别方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108333627A true CN108333627A (zh) | 2018-07-27 |
| CN108333627B CN108333627B (zh) | 2019-01-08 |
Family
ID=62925124
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810040119.XA Expired - Fee Related CN108333627B (zh) | 2018-01-16 | 2018-01-16 | 火成岩地区断裂真伪的识别方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108333627B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109709605A (zh) * | 2019-02-28 | 2019-05-03 | 中国石油化工股份有限公司 | 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 |
| CN109736796A (zh) * | 2019-03-04 | 2019-05-10 | 中铁十六局集团有限公司 | 一种加深炮孔超前探测机构及其预报方法 |
| CN112147684A (zh) * | 2019-06-28 | 2020-12-29 | 中国石油天然气股份有限公司 | 表征同沉积断层活动强度的方法和装置 |
| CN112764107A (zh) * | 2019-11-04 | 2021-05-07 | 中国石油天然气集团有限公司 | 地表高速覆盖地区的地震资料静校正方法及装置 |
| CN113589368A (zh) * | 2020-04-30 | 2021-11-02 | 中国石油化工股份有限公司 | 火成岩影响地震反射特征变化及规律的定量分析方法 |
| CN113589370A (zh) * | 2020-04-30 | 2021-11-02 | 中国石油天然气股份有限公司 | 高速地层覆盖下的构造识别方法及装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101285894A (zh) * | 2008-05-30 | 2008-10-15 | 中国科学院地质与地球物理研究所 | 起伏地表下采集的地震资料的直接叠前时间偏移方法 |
| CN102176053A (zh) * | 2011-01-27 | 2011-09-07 | 中国科学院地质与地球物理研究所 | 提升波动方程叠前深度偏移成像效果的方法 |
| CN104375172A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油集团东方地球物理勘探有限责任公司 | 火山岩下伏地层的构造形态正确成像的方法 |
-
2018
- 2018-01-16 CN CN201810040119.XA patent/CN108333627B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101285894A (zh) * | 2008-05-30 | 2008-10-15 | 中国科学院地质与地球物理研究所 | 起伏地表下采集的地震资料的直接叠前时间偏移方法 |
| CN102176053A (zh) * | 2011-01-27 | 2011-09-07 | 中国科学院地质与地球物理研究所 | 提升波动方程叠前深度偏移成像效果的方法 |
| CN104375172A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油集团东方地球物理勘探有限责任公司 | 火山岩下伏地层的构造形态正确成像的方法 |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109709605A (zh) * | 2019-02-28 | 2019-05-03 | 中国石油化工股份有限公司 | 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 |
| CN109709605B (zh) * | 2019-02-28 | 2020-09-01 | 中国石油化工股份有限公司 | 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 |
| CN109736796A (zh) * | 2019-03-04 | 2019-05-10 | 中铁十六局集团有限公司 | 一种加深炮孔超前探测机构及其预报方法 |
| CN112147684A (zh) * | 2019-06-28 | 2020-12-29 | 中国石油天然气股份有限公司 | 表征同沉积断层活动强度的方法和装置 |
| CN112147684B (zh) * | 2019-06-28 | 2023-12-26 | 中国石油天然气股份有限公司 | 表征同沉积断层活动强度的方法和装置 |
| CN112764107A (zh) * | 2019-11-04 | 2021-05-07 | 中国石油天然气集团有限公司 | 地表高速覆盖地区的地震资料静校正方法及装置 |
| CN112764107B (zh) * | 2019-11-04 | 2024-05-28 | 中国石油天然气集团有限公司 | 地表高速覆盖地区的地震资料静校正方法及装置 |
| CN113589368A (zh) * | 2020-04-30 | 2021-11-02 | 中国石油化工股份有限公司 | 火成岩影响地震反射特征变化及规律的定量分析方法 |
| CN113589370A (zh) * | 2020-04-30 | 2021-11-02 | 中国石油天然气股份有限公司 | 高速地层覆盖下的构造识别方法及装置 |
| CN113589368B (zh) * | 2020-04-30 | 2023-12-01 | 中国石油化工股份有限公司 | 火成岩影响地震反射特征变化及规律的定量分析方法 |
| CN113589370B (zh) * | 2020-04-30 | 2024-04-30 | 中国石油天然气股份有限公司 | 高速地层覆盖下的构造识别方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108333627B (zh) | 2019-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108333627B (zh) | 火成岩地区断裂真伪的识别方法及装置 | |
| CN104977618B (zh) | 一种评价页岩气储层及寻找甜点区的方法 | |
| US8793113B2 (en) | Method and apparatus for near well structural modeling based on borehole dips | |
| Radwan | Three-dimensional gas property geological modeling and simulation | |
| CN106569267A (zh) | 一种致密低渗透储层多尺度裂缝模型及建模方法 | |
| CN108680951A (zh) | 一种基于地震信息判断煤层气富集沉积控制作用的方法 | |
| EP3253946B1 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
| CN109884710B (zh) | 针对激发井深设计的微测井层析成像方法 | |
| CN104200115B (zh) | 一种基于地质统计学模拟的全地层速度建模方法 | |
| US8498177B2 (en) | Determining a position of a geological layer relative to a wavelet response in seismic data | |
| CN102636812B (zh) | 一种获得碳酸盐岩储层储集空间体积的方法 | |
| CN108397130B (zh) | 钻井方法和装置 | |
| CN103635831A (zh) | 在3d各向异性介质中的基于敏感度内核的偏移速度分析 | |
| CN105549082B (zh) | 超深碳酸盐储层三维地质力学场建立方法和系统 | |
| CN106285642B (zh) | 一种基于地震资料的地层孔隙压力预测方法 | |
| CN105911592B (zh) | 一种实钻约束的三维地震资料长波长静校正方法 | |
| CN105549084A (zh) | 一种三维高精度速度建模方法和系统 | |
| US10451769B2 (en) | Method for determining petrophysical properties from logging measurements | |
| CN104036119B (zh) | 划分沉积地层的方法 | |
| CN107797148B (zh) | 一种基于三维地质建模的航磁异常场分离方法及系统 | |
| CN104749618B (zh) | 泥页岩缓倾角裂缝叠后概率定量表征方法 | |
| CN111999768B (zh) | 一种钻井分层设计中基于构造模型的构造海拔确定方法 | |
| Singleton | Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas | |
| CN105093330B (zh) | 多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法 | |
| NO344460B1 (en) | Methods and systems for identifying and plugging subterranean conduits |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190108 Termination date: 20210116 |