CN104360385B - 结合地震地质认识综合建立初始深度层速度模型的方法 - Google Patents
结合地震地质认识综合建立初始深度层速度模型的方法 Download PDFInfo
- Publication number
- CN104360385B CN104360385B CN201410693455.6A CN201410693455A CN104360385B CN 104360385 B CN104360385 B CN 104360385B CN 201410693455 A CN201410693455 A CN 201410693455A CN 104360385 B CN104360385 B CN 104360385B
- Authority
- CN
- China
- Prior art keywords
- interval velocity
- information
- layer
- model
- velocity
- 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 56
- 238000010276 construction Methods 0.000 claims abstract description 32
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 22
- 238000013508 migration Methods 0.000 claims description 64
- 230000005012 migration Effects 0.000 claims description 64
- 238000012545 processing Methods 0.000 claims description 26
- 230000001973 epigenetic effect Effects 0.000 claims description 21
- 238000004458 analytical method Methods 0.000 claims description 18
- 241001074085 Scophthalmus aquosus Species 0.000 claims description 16
- 238000005553 drilling Methods 0.000 claims description 11
- 238000005457 optimization Methods 0.000 claims description 11
- 238000007781 pre-processing Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 5
- 230000003068 static effect Effects 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000003786 synthesis reaction Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 20
- 238000005755 formation reaction Methods 0.000 description 16
- 238000003325 tomography Methods 0.000 description 13
- 238000009826 distribution Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 241001270131 Agaricus moelleri Species 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明的结合地震地质认识综合建立初始深度层速度模型的方法综合运用了地震数据反演、钻测井资料、地表地质露头资料、区域地质信息、地层速度变化规律等信息来建立初始深度层速度模型，即以地震数据反演为基础，以钻测井资料、地表地质露头资料、区域地质信息、地层速度变化规律为约束，获得了较为准确的初始构造模型及沿层层速度，进而形成了可实施的初始深度层速度建模流程。本发明充分吸收了现有的三种初始深度层速度模型建立方法的优点，同时采用了多信息约束的思路综合建立初始深度层速度模型，使得初始深度层速度模型更加准确。
Description
技术领域
本发明属于石油天然气地震勘探领域，更具体地讲，涉及叠前深度偏移速度建模方面的一种结合地震地质认识综合建立初始深度层速度模型的方法，主要应用于石油地震勘探的地震资料处理解释。
背景技术
随着石油天然气勘探程度和勘探对象复杂化程度的不断提高，对地震偏移成像质量和精度的要求也越来越高。当地下介质速度存在剧烈的横向变化或速度分界面不是水平层状时，只有叠前深度偏移能够实现共反射点的叠加和绕射点的归位，修正陡倾地层和速度变化产生的地下图像畸变，使复杂构造或速度横向变化较大的地震资料正确成像。而地震资料处理中影响叠前深度偏移效果的关键因素是计算旅行时用的深度层速度。
深度层速度建模一般包括两个方面：初始深度层速度模型建立和深度层速度模型优化。深度层速度模型优化多采用层析成像反演，如果初始深度层速度模型的误差太大或趋势错误，则可能造成层析成像反演结果朝相反的方向发展，以至于难以获得正确的速度模型。因此，求取相对合理的初始深度层速度模型非常重要。
目前，建立初始深度层速度模型的方法一般是基于地震数据反演进行的，从简单到复杂主要有三种方法：一种比较简单的方法就是利用叠加速度并通过平滑处理后直接进行时深转换到深度域作为初始深度层速度，这种方法的速度精度降低，误差较大，仅适合于地震资料信噪比高，速度简单的地震资料；第二种方法是利用叠前时间偏移速度进行时深转换到深度域作为初始深度层速度，该方法较第一种方法求取的初始深度层速度更合理，但由于叠前时间偏移速度多体现为均方根速度，通过速度转换到层速度可能造成层速度不合理，造成初始深度层速度速度模型精度不高，该方法也仅适合于高信噪比资料；第三种方法是利用叠前时间偏移求取瞬时层速度，同时在叠前时间偏移剖面上解释地质层位，沿地质层位提取瞬时层速度，通过平滑处理得到时间域层速度，再通过时深转换或图偏移等方法获得初始深度层速度，该方法较前两种方法的优点在于考虑了地层层速度，是目前复杂构造速度建模推荐的较好方法，但该方法仍然受地震资料条件的限制，即在地震资料的信噪比低、构造复杂的情况下，获得的时间域瞬时层速度也可能精度较低，同样会造成初始深度层速度误差较大的后果。
发明内容
针对现有技术中存在的不足，本发明的目的在于解决上述技术问题中的一个或多个，提供一种将地表地层分布规律与钻测井获得的地下地层层速度结合起来并采用多信息约束综合建立初始深度层速度模型的方法。
为了实现上述目的，本发明提供了一种结合地震地质认识综合建立初始深度层速度模型的方法，所述方法包括以下步骤：
A、将地震数据进行叠前预处理和动校正叠加处理，获得水平叠加剖面；并在叠前预处理后的地震数据基础上进行叠前时间偏移处理，获得叠前时间偏移剖面和叠前时间偏移速度；
B、利用所述水平叠加剖面进行地震波场分析并确定构造基本格架，利用勘探区的地面露头资料和钻测井资料对所述叠前时间偏移剖面建立沿地震测线的地震剖面浅层构造模型并恢复井旁构造；
C、结合所述构造基本格架、井旁构造、所述沿地震测线的地震剖面浅层构造模型对叠前时间偏移剖面进行断层划分和时间偏移域构造层位解释，获得时间域构造模型；
D、通过控制速度反演将所述叠前时间偏移速度转换为瞬时层速度，并利用所述瞬时层速度在所述时间域构造模型的控制下提取各构造层位的层速度，获得基于地震反演的层速度；
E、利用勘探区的钻测井资料求取井点层速度并进行区域速度分析，结合所述井点层速度、区域速度分析结果和利用勘探区的钻测井资料获取的地质信息对所述时间域构造模型进行速度填充，获得沿地震测线的层速度变化趋势；再利用勘探区平面上的多条地震测线上的沿地震测线的层速度变化趋势建立区域速度量板；
F、利用所述区域速度量板修正所述基于地震反演的层速度，获得修正后的初始沿层层速度；
G、将所述初始沿层层速度填充所述时间域构造模型，获得时间域初始层速度模型；再对所述时间域初始层速度模型进行时深转换或图偏移处理，获得初始深度层速度模型。
根据本发明的结合地震地质认识综合建立初始深度层速度模型的方法的一个实施例，所述叠前时间偏移速度是进行叠前时间偏移均方根速度分析获得的均方根速度。
根据本发明的结合地震地质认识综合建立初始深度层速度模型的方法的一个实施例，所述地面露头资料至少包括地层出露岩性信息、地质层位信息、倾向信息、倾角信息和钻井信息中的一个或多个。
根据本发明的结合地震地质认识综合建立初始深度层速度模型的方法的一个实施例，所述钻测井资料至少包括地质层位信息、地层厚度信息、倾向信息、倾角信息、断层信息、浅层倾角信息和测井速度信息中的一个或多个。
根据本发明的结合地震地质认识综合建立初始深度层速度模型的方法的一个实施例，在步骤C中，先将所述沿地震测线的地震剖面浅层构造模型结合所述叠前时间偏移剖面进行浅层层位解释，并根据所述浅层层位解释将所述井旁构造在所述沿地震测线的地震剖面浅层构造模型上进行追踪对比和加密层位解释，再结合所述构造基本格架进行断层划分和时间偏移域层位解释，获得时间域构造模型。
根据本发明的结合地震地质认识综合建立初始深度层速度模型的方法的一个实施例，在步骤D中，将所述瞬时层速度进行弯曲射线叠前时间偏移和瞬时层速度更新，获得优化瞬时层速度，再利用所述优化瞬时层速度进行后续处理。
根据本发明的结合地震地质认识综合建立初始深度层速度模型的方法的一个实施例，在步骤E中，所述利用勘探区的钻测井资料获取的地质信息至少包括岩性资料和区域地质认识资料。
根据本发明的结合地震地质认识综合建立初始深度层速度模型的方法的一个实施例，所述叠前预处理包括静校正、叠前去噪、振幅处理和反褶积处理。
本发明综合运用了地震数据反演、钻测井资料、地表地质露头资料、区域地质信息、地层速度变化规律等信息来建立初始深度层速度模型，即以地震数据反演为基础，以钻测井资料、地表地质露头资料、区域地质信息、地层速度变化规律为约束，获得了较为准确的初始构造模型及沿层层速度，进而形成了可实施的初始深度层速度建模流程。本发明充分吸收了现有的三种初始深度层速度模型建立方法的优点，同时采用了多信息约束的思路综合建立初始深度层速度模型，使得初始深度层速度模型更加准确。
附图说明
图1是根据本发明的示例性实施例的结合地震地质认识综合建立初始深度层速度模型的方法的流程图。
图2是根据本发明的示例获得的沿地震测线的地震剖面浅层构造模型的示意图。
图3是根据本发明的示例获得的井旁构造的示意图。
图4是根据本发明的示例获得的时间域构造模型的示意图。
图5是根据本发明的示例获得的基于地震反演的层速度的示意图。
图6是根据本发明的示例获得的井点层速度的示意图。
图7是根据本发明的示例获得的沿地震测线的层速度变化趋势的示意图。
图8是根据本发明的示例获得的区域速度量板的示意图。
图9是根据本发明的示例获得的修正后的初始沿层层速度的示意图。
图10是根据本发明的示例获得的时间域初始层速度模型的示意图。
图11是根据本发明的示例获得的初始深度层速度模型的示意图。
在附图中，主要标记的含义为：
地质层位类：
T3xj-须底、T1j2 2-嘉二2底、T1f4-飞四底、T1f1-飞底、P2l-一叠系上统底、P1l-二叠系下统底；
井位类：
B1-板1井、B2-板2井、BD3-板东3井、BD10-板东10井、L1-凉1井、LD6-凉东6井；
断层类：
B21-板21断层、B3-板3断层。
具体实施方式
在下文中，将参照附图来详细说明根据本发明的结合地震地质认识综合建立初始深度层速度模型的方法。
可通过计算机软件模块实现本发明的结合地震地质认识综合建立初始深度层速度模型的方法。在本发明中描述的各个步骤不限制为上述步骤，其中的一些步骤可被进一步拆分为更多的步骤，并且一些步骤可合并为较少的步骤。
本发明的结合地震地质认识综合建立初始深度层速度模型的方法主要包括两个方面的地震地质认识相结合的处理，即地震地质认识结合建立时间域构造模型和地震地质认识结合求取初始沿层层速度，并最终获得相对较准确的初始深度层速度模型。
一般来讲，时间域构造模型是指地质构造如地下岩层或岩体在地质构造运动中遗留下来的各种构造形迹，如褶皱、断层等。在地震勘探中，由于地震数据采集以时间进行记录，所以纵轴以时间单位进行表示的构造模型即作为时间域构造模型。初始沿层层速度是相对于速度模型更新而言求取的第一轮层速度。初始深度层速度模型即是将速度模型以深度域表示并且层速度为初始层速度的速度模型，该模型即为深度层速度模型更新的初始模型，为建立准确的速度层速度模型奠定基础。另外，本发明还涉及时间域初始层速度模型的概念，层速度本身是无时间域和深度域的差别的，当将速度模型以时间域表示并且层速度为初始层速度的速度模型即为时间域初始层速度模型。
此外，在本领域中，层位通常是指地下某地层的底界或顶界，并且通常利用层位来表示相应的地层。在本发明的描述中，由于层位与地层具有基本类似的使用作用，所以为了描述的方便，本发明不对层位和地层进行详细区分，本领域技术人员均了解二者的区别。
图1是根据本发明的示例性实施例的结合地震地质认识综合建立初始深度层速度模型的方法的流程图。在示例性实施例中，假设在通用计算机系统上通过实现各步骤的功能的软件模块执行所述结合地震地质认识综合建立初始深度层速度模型的方法。
如图1所示，根据本发明的示例性实施例，结合地震地质认识综合建立初始深度层速度模型的方法可以包括以下的各步骤。
在建立初始深度层模型之前，需要对勘探区的地震数据进行预先处理以获得能够进行后续处理的符合要求的资料。根据本发明的示例性实施例，将地震数据进行叠前预处理和动校正叠加处理，获得水平叠加剖面；并在叠前预处理后的地震数据基础上进行叠前时间偏移处理，获得叠前时间偏移剖面和叠前时间偏移速度。
其中，叠前预处理包括静校正、叠前去噪、振幅处理和反褶积处理，但本发明不限于此。并且，叠前时间偏移速度是进行叠前时间偏移均方根速度分析获得的均方根速度。
获取水平叠加剖面的作用之一是可以求取准确的剩余静校正量，之二则是可以通过水平叠加剖面认识地震波场并认识基本构造格架。获取叠前时间偏移剖面的作用也有两点，一是叠前时间偏移剖面本身是地震成果，可进行构造解释；二是可以获得叠前时间偏移速度，而叠前时间偏移速度可以为后续层速度的求取奠定基础。
接下来进行结合地震地质认识综合建立初始深度层速度模型的各步骤。
首先，进行地震地质认识结合建立时间域构造模型的操作。
在步骤S101，利用水平叠加剖面进行地震波场分析并确定构造基本格架，其目的在于根据水平叠加剖面上的地震波场初步认识地下(如地腹地层)的相对构造关系。
在步骤S102，利用勘探区的地面露头资料和钻测井资料恢复井旁构造，以获得井旁准确的构造样式。其中，所述的井是指以勘探开发石油和天然气为目的而在地层中钻出的具有一定深度的圆柱形孔眼。通过钻井和测井可以获取勘探区地层的物理资料、地质资料等，从而可以为地震勘探和地震资料解释提供信息和数据。恢复井旁构造之后能够获得准确的井旁地层构造形态，以便准确认识和解释地震剖面。
在步骤S103，利用勘探区的地面露头资料和钻测井资料对叠前时间偏移剖面建立沿地震测线的地震剖面浅层构造模型，利用上述步骤获得的沿地震测线的地震剖面浅层构造模型进行后续处理能够使浅层信噪比低的地震剖面构造解释更为准确。
其中，本发明中涉及的地面露头资料可以至少包括地层出露岩性信息、地质层位信息、倾向信息、倾角信息和钻井信息中的一个或多个；钻测井资料可以至少包括地质层位信息、地层厚度信息、倾向信息、倾角信息、断层信息、浅层倾角信息和测井速度信息中的一个或多个，但本发明不限于此。此外，所利用的地面露头资料是地质人员在野外进行实地勘测获得的，而钻测井资料则是通过钻井时取岩心、录井、测井等方式获得的，上述地面露头资料和钻测井资料的获取均可以本领域公知的方法进行。
根据本发明的一个实施例，利用地面露头的地质层位信息、倾角信息和钻井信息进行地震剖面地质戴帽来建立沿地震测线的地震剖面浅层构造模型，同时结合钻测井获得的倾角信息、层位信息等信息来恢复井旁构造，但本发明不限于此。
在步骤S104，结合构造基本格架、井旁构造、沿地震测线的地震剖面浅层构造模型对叠前时间偏移剖面进行断层划分和时间偏移域构造层位解释，获得时间域构造模型。其中，断层划分和时间偏移构造层位解释的思路和方法可以有很多种，基本的思路是在符合地质规律前提下，对地震剖面进行地质层位认识、追踪对比解释，对地质层位错断的地方进行断层划分并且保证断层划分符合构造运动规律，具体可以采用本领域公知的方法进行。
根据本发明的一个实施例，先将沿地震测线的地震剖面浅层构造模型结合叠前时间偏移剖面进行浅层层位解释，并根据浅层层位解释将井旁构造在沿地震测线的地震剖面浅层构造模型上进行追踪对比和加密层位解释，再结合构造基本格架进行断层划分和时间偏移域层位解释，获得时间域构造模型。这里，浅层层位解释仅仅是针对浅层地层进行的层位解释，事实上，其与时间偏移域层位解释是类似的，只是时间偏移域层位解释所针对的是全部地层。之所以优选地先进行浅层层位解释，是因为由于浅层地震剖面一般质量较差，所以导致直接在地震剖面上的解释不太准确，在结合地面露头资料及井旁构造的信息条件下才能更准确进行浅层构造解释。
至此，则基本完成了地震地质认识结合建立时间域构造模型的步骤，由于在建立时间域构造模型的过程中充分运用地震数据反演、钻测井资料、区域地质认识、地层速度变化规律等信息，所以使得获得的时间域构造模型相对于仅通过地震数据反演的结果更为准确。
然后，进行地震地质认识结合求取初始沿层层速度的操作。
在步骤S201，通过控制速度反演将叠前时间偏移速度转换为瞬时层速度。瞬时层速度是指速度体中把每个采样点作为层考虑的层速度，其作用是可以用来提取地质层位的层速度。而控制速度反演是通过控制点均方根速度，根据均方根速度转换到层速度的公式并结合人为控制平滑等参数，以一定比例进行综合运算获得需要的速度的处理，这是本领域所公知的方法，在此不进行赘述。
由于叠前时间偏移速度为均方根速度，所以由其通过计算直接转换成的瞬时层速度的精度不太高，为了获得更加符合要求的瞬时层速度，还可以对瞬时层速度进行优化处理。根据本发明的示例性实施例，将瞬时层速度进行弯曲射线叠前时间偏移和瞬时层速度更新，获得优化瞬时层速度，之后再利用优化瞬时层速度进行后续处理即可。优化处理步骤具体可以在进行弯曲射线叠前时间偏移之后，根据获得的道集进行拉平处理和速度修改更新，直到获得满意的瞬时层速度即可作为优化瞬时层速度。
在步骤S401，利用瞬时层速度在时间域构造模型的控制下提取各构造层位的层速度，获得基于地震反演的层速度。
由于假设的地下地层是层状的，且同一个地层具有相等或相近的层速度，但目前的地震数据反演不能直接反演出层速度，而需要先反演求取均方根速度之后，再利用获得的瞬时层速度进行各构造层位的层速度提取，通过这样的处理步骤获得的层速度实际上是基于地震数据反演得到的相对近似的层速度值，更为准确。例如，具体的提取方法可以是直接利用层位底部所在的速度样点值作为层速度值，但本发明不限于此。
各构造层位的层速度即地下各地层的层速度，而地层的层速度是指地震波在地层中的传播速度。所提取的各构造层位的层速度在这里是作为初始层速度来利用的，即作为基于地震反演的层速度，然后在后续的处理中再进行速度更新以获得更为准确的层速度，最后用来进行更为准确的偏移成像。
在步骤S301，利用勘探区的钻测井资料求取井点层速度并进行区域速度分析，结合井点层速度、区域速度分析结果和利用勘探区的钻测井资料获取的地质信息对时间域构造模型进行速度填充，获得沿地震测线的层速度变化趋势。
其中，井点层速度是根据测井数据解释获得的地下地层层速度，相比地震处理速度分析获得的速度准确度更高，井点层速度可以理解为是该井点位置地下准确的层速度，因此可用来帮助认识和求取地下三维空间的层速度，提高速度求取的准确性。根据本发明的一个实施例，可以利用测井速度来求取井点层速度。
此外，所涉及的勘探区的钻测井资料获取的地质信息可以至少包括岩性资料和区域地质认识资料。根据本发明的一个实施例，利用钻井获得的岩性资料和区域地质认识资料结合时间域构造模型获取沿地震测线的层速度变化趋势。速度填充可以是以有地质层位解释的叠前时间偏移剖面(即时间域构造模型)为底图，根据钻井获得的井点层速度及分析获得的区域速度分析结果进行填充。
在步骤S302，利用勘探区平面上的多条地震测线上的沿地震测线的层速度变化趋势建立区域速度量板。
区域速度量板即为平面上的趋势速度分布趋势，确切地说，即为每个地层层速度在平面上的分布趋势，其作用是帮助技术人员进行速度认识。
事实上，在之前的步骤中，需要针对每条地震测线建立地震剖面浅层构造模型，然后分别进行处理，直到获得沿每条地震测线的层速度变化趋势。通过多条地震测线上的沿地震测线的层速度变化趋势建立区域速度量板之后，就相当于将测线方向上的二维速度变化扩推到了勘探区平面上的三维速度变化，对整体认识勘探区速度变化及后续层速度解释奠定了基础。
在步骤S402中，利用区域速度量板修正基于地震反演的层速度，获得修正后的初始沿层层速度。
其中，修正具体是根据勘探区各地层区域速度变化趋势，对基于地震反演的层速度的异常值、不符合规律的速度点进行修改，直至得到符合速度分步规律的层速度并将其作为修正后的初始沿层层速度。
至此，则基本完成了地震地质认识结合求取初始沿层层速度的操作，同样地，由于在求取初始沿层层速度的过程中充分运用地震数据反演、钻测井资料、区域地质认识、地层速度变化规律等信息，所以使得获得的初始沿层层速度相对于仅通过地震数据反演的结果更为准确。在获得的更为准确的初始沿层层速度的基础上，可以直接进行速度填充而建立所需的速度模型。
在步骤S403，将初始沿层层速度填充时间域构造模型，获得时间域初始层速度模型。在步骤S404，再对时间域初始层速度模型进行时深转换或图偏移处理，从而获得初始深度层速度模型。
总之，较之以往基于地震数据反演建立初始深度层速度模型的方法，本发明的方法增加了钻测井资料、地表地质露头资料、区域地质认识信息、地层速度变化规律等信息的约束，修正了仅基于地震数据反演建立的初始深度层速度模型，形成了可实施的初始深度层速度建模流程，获得了更为准确的初始深度层速度模型。
下面结合本发明的结合地震地质认识综合建立初始深度层速度模型的方法的一个具体示例来进一步地详细说明。
申请人将本发明的结合地震地质认识建立初始深度层速度模型的方法在四川川东地区进行了应用，其中涉及的主要标记的含义包括：
地质层位类：T3xj-须底、T1j2 2-嘉二2底、T1f4-飞四底、T1f1-飞底、P2l-二叠系上统底、P1l-二叠系下统底；
井位类：B1-板1井、B2-板2井、BD3-板东3井、BD10-板东10井、L1-凉1井、LD6-凉东6井；
断层类：B21-板21断层、B3-板3断层。
上述各层位、井位、断层的标记名称为本领域技术人员均了解的名称并且具有已知的涵义，例如须底为须家河组地层底界，嘉二2底为嘉陵江组二段第二层底界，飞四底为飞仙关组四段底界，飞底为飞仙关组底界等，在此不作赘述。
具体步骤如下：
1)将勘探区的地震数据进行静校正、叠前去噪、振幅处理、反褶积处理的叠前预处理，之后进行动校正叠加处理获得水平叠加剖面。对所得的水平叠加剖面进行地震波场分析并确定构造基本格架。
2)在叠前预处理后的地震数据基础上进行叠前时间偏移处理，获得叠前时间偏移剖面和叠前时间偏移速度。
3)综合运用地面露头资料的地层出露岩性信息、地质层位信息、倾向信息、倾角信息以及钻井获得的浅层倾角信息绘制沿地震测线的地震剖面浅层构造模型，并将所得的沿地震测线的地震剖面浅层构造模型结合叠前时间偏移剖面进行浅层层位解释。
图2是根据本发明的示例获得的沿地震测线的地震剖面浅层构造模型的示意图。如图2所示，图2的最上沿是地表平面的地质截图，其上标有地表地层的倾角信息，而图2的最上沿下面的部分则为根据地表露头资料及钻井(3口)获得的地层信息在叠前时间偏移剖面上绘制的沿地震测线的地震剖面浅层构造模型，其中恢复了浅层地层出露层位信息、产状信息及倾角信息等信息。由此可对地震剖面的浅层构造进行更准确的构造解释，对于浅层信噪比较低的剖面尤其有效。
4)根据勘探区的地面露头资料和钻井轨迹以及获得的准确的地质层位信息、地层厚度信息、倾向信息、倾角信息、断层信息等信息，地质人员根据地质规律恢复了井旁构造。
图3是根据本发明的示例获得的井旁构造的示意图。如图3所示，B1(板1井)的钻井过程中钻遇了飞四地层、二叠系上统地层、二叠系下统地层，而且钻遇了P2l(二叠系上统底)和P1l(二叠系下统底)各两次，由此恢复出了井旁地质层位及井点上的B3(板3断层)的位置，结合地震剖面解释出B3。其中，B21(板21断层)不是通过井恢复的，而是根据地震剖面解释的，为了保证构造完整性而将B21补入了图3中。
5)根据浅层层位解释，将井旁构造在沿地震测线的地震剖面浅层构造模型上进行追踪对比和加密层位解释，再结合构造基本格架进行断层划分和时间偏移域层位解释，获得时间域构造模型。
图4是根据本发明的示例获得的时间域构造模型的示意图。如图4所示，图中的深色实线表示断层，不同灰度的虚线表示构造解释层位，层位的名字也已经标出，各井名称及下方竖线表示井位置及井轨迹。
6)通过控制速度反演将叠前时间偏移速度转换为瞬时层速度，并将瞬时层速度进行弯曲射线叠前时间偏移和瞬时层速度更新，获得优化瞬时层速度。再利用优化瞬时层速度在时间域构造模型的控制下提取各构造层位的层速度，获得基于地震反演的层速度，以在后续的步骤中进行修正。
图5是根据本发明的示例获得的基于地震反演的层速度的示意图。如图5所示，不同灰度代表不同地层，纵坐标即为层速度值，横坐标为测线CMP点(描述测线所在位置)。可见对该地震数据来说，基于地震反演获得的层速度变化较大并且存在异常值，需要修正。
7)根据测井速度求取井点层速度并进行区域速度分析，认识地下地层的层速度分布规律。
图6是根据本发明的示例获得的井点层速度的示意图。如图6所示，不同灰度的方块代表了不同的地层且地层的名字已经标示出来，井点位置所处的各地层的平均层速度范围也已经标出。由此可知，根据钻测井资料，能够获得较准确的井点位置的地层平均层速度。
8)结合井点层速度、区域速度分析结果和利用勘探区的钻测井资料获取的地质信息对时间域构造模型进行速度填充，获得沿地震测线的层速度变化趋势。
图7是根据本发明的示例获得的沿地震测线的层速度变化趋势的示意图。如图7所示，通过钻测井资料获取的地震剖面层位解释，结合井点层速度以及不同深度同一地层的层速度变化趋势，就能够获得沿地震测线的层速度变化趋势，如二叠系下统地层的层速度范围为5700-6000米/秒，在较浅位置速度取相对低值，较深位置取相对高值。
9)利用勘探区平面上的多条地震测线上的沿地震测线的层速度变化趋势建立区域速度量板。
图8是根据本发明的示例获得的区域速度量板的示意图。如图8所示，图8具体示出的是须底-嘉二2底地层的区域速度量板，可见同一地层由于埋藏深度、上覆地层变化的影响，层速度在平面上也存在有规律的变化趋势，如埋藏较深则速度较大，而构造顶部埋藏较浅则速度较小。
10)利用区域速度量板修正基于地震反演的层速度，获得修正后的初始沿层层速度。
图9是根据本发明的示例获得的修正后的初始沿层层速度的示意图。如图9所示，不同灰度代表不同地层的层速度，相较于图5中示出的不同地层的层速度，图9中结合地震地质认识的层速度已经被去除了异常值，层速度的变化规律符合平面变化趋势且层速度更加合理。
11)将初始沿层层速度填充时间域构造模型，获得时间域初始层速度模型。
图10是根据本发明的示例获得的时间域初始层速度模型的示意图。如图10所示，时间域初始层速度模型是通过时间域构造模型与修正后的层速度填充而成的。图10的横坐标为CMP点(共中心点)，纵坐标为时间，色标表示层速度值，该图即为建立的所述地震测线的地下时间域初始层速度模型，反映地下各地层的速度分布趋势。
12)再对时间域初始层速度模型进行时深转换或图偏移处理，从而获得初始深度层速度模型。
图11是根据本发明的示例获得的初始深度层速度模型的示意图。如图11所示，该模型即为较准确合理的初始深度层速度模型，该模型除了利用常规的基于地震反演的方法外，还融合了地表露头资料、钻测井资料等地质信息，相对于对低信噪比、复杂构造的地震数据来说，建立的初始深度层速度模型相对更加准确、合理。图11的横坐标为CMP点(共中心点)，纵坐标为深度，色标表示层速度值，该图即为建立的所述地震测线的地下深度域初始层速度模型，反映地下各地层的速度分布趋势。
本示例建立了与地面及地下钻井相吻合的初始深度层速度模型，通过后期层析成像反演速度模型优化，获得了较理想的速度层速度模型，最终获得了成像更清楚的叠前深度偏移剖面，取得了较明显的应用效果。
综上所述，本发明适用于具有地面露头资料、钻测井资料等信息的所有地震数据的初始深度层速度建模，对于复杂地表及地腹构造、信噪比相对较低地区的优势明显，在复杂山地地震资料处理中具有广阔的应用前景。
尽管已参照示例性实施例表示和描述了本发明，但本领域技术人员应该理解，在不脱离由权利要求限定的本发明的精神和范围的情况下，可以对这些实施例进行各种修改和变形。
Claims (6)
1.一种结合地震地质认识综合建立初始深度层速度模型的方法，其特征在于，所述方法包括以下步骤：
A、将地震数据进行叠前预处理和动校正叠加处理，获得水平叠加剖面；并在叠前预处理后的地震数据基础上进行叠前时间偏移处理，获得叠前时间偏移剖面和叠前时间偏移速度；
B、利用所述水平叠加剖面进行地震波场分析并确定构造基本格架，利用勘探区的地面露头资料和钻测井资料对所述叠前时间偏移剖面建立沿地震测线的地震剖面浅层构造模型并恢复井旁构造；
C、结合所述构造基本格架、井旁构造、所述沿地震测线的地震剖面浅层构造模型对叠前时间偏移剖面进行断层划分和时间偏移域构造层位解释，获得时间域构造模型；
D、通过控制速度反演将所述叠前时间偏移速度转换为瞬时层速度，并利用所述瞬时层速度在所述时间域构造模型的控制下提取各构造层位的层速度，获得基于地震反演的层速度；
E、利用勘探区的钻测井资料求取井点层速度并进行区域速度分析，结合所述井点层速度、区域速度分析结果和利用勘探区的钻测井资料获取的地质信息对所述时间域构造模型进行速度填充，获得沿地震测线的层速度变化趋势；再利用勘探区平面上的多条地震测线上的沿地震测线的层速度变化趋势建立区域速度量板；
F、利用所述区域速度量板修正所述基于地震反演的层速度，获得修正后的初始沿层层速度；
G、将所述初始沿层层速度填充所述时间域构造模型，获得时间域初始层速度模型；再对所述时间域初始层速度模型进行时深转换或图偏移处理，获得初始深度层速度模型，
其中，在步骤C中，先将所述沿地震测线的地震剖面浅层构造模型结合所述叠前时间偏移剖面进行浅层层位解释，并根据所述浅层层位解释将所述井旁构造在所述沿地震测线的地震剖面浅层构造模型上进行追踪对比和加密层位解释，再结合所述构造基本格架进行断层划分和时间偏移域层位解释，获得时间域构造模型，
其中，在步骤D中，将所述瞬时层速度进行弯曲射线叠前时间偏移和瞬时层速度更新，获得优化瞬时层速度，再利用所述优化瞬时层速度进行后续处理。
2.根据权利要求1所述的结合地震地质认识综合建立初始深度层速度模型的方法，其特征在于，所述叠前时间偏移速度是进行叠前时间偏移均方根速度分析获得的均方根速度。
3.根据权利要求1所述的结合地震地质认识综合建立初始深度层速度模型的方法，其特征在于，所述地面露头资料至少包括地层出露岩性信息、地质层位信息、倾向信息、倾角信息和钻井信息中的一个或多个。
4.根据权利要求1所述的结合地震地质认识综合建立初始深度层速度模型的方法，其特征在于，所述钻测井资料至少包括地质层位信息、地层厚度信息、倾向信息、倾角信息、断层信息、浅层倾角信息和测井速度信息中的一个或多个。
5.根据权利要求1所述的结合地震地质认识综合建立初始深度层速度模型的方法，其特征在于，在步骤E中，所述利用勘探区的钻测井资料获取的地质信息至少包括岩性资料和区域地质认识资料。
6.根据权利要求1所述的结合地震地质认识综合建立初始深度层速度模型的方法，其特征在于，所述叠前预处理包括静校正、叠前去噪、振幅处理和反褶积处理。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410693455.6A CN104360385B (zh) | 2014-11-26 | 2014-11-26 | 结合地震地质认识综合建立初始深度层速度模型的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410693455.6A CN104360385B (zh) | 2014-11-26 | 2014-11-26 | 结合地震地质认识综合建立初始深度层速度模型的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104360385A CN104360385A (zh) | 2015-02-18 |
| CN104360385B true CN104360385B (zh) | 2017-09-22 |
Family
ID=52527665
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410693455.6A Active CN104360385B (zh) | 2014-11-26 | 2014-11-26 | 结合地震地质认识综合建立初始深度层速度模型的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104360385B (zh) |
Families Citing this family (32)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105259581B (zh) * | 2015-10-22 | 2020-05-26 | 中国石油化工股份有限公司 | 一种地震资料时深转换方法 |
| CN105842736B (zh) * | 2016-05-27 | 2018-07-10 | 中国石油集团东方地球物理勘探有限责任公司 | 一种构建地层速度模型的方法 |
| CN105911593A (zh) * | 2016-06-24 | 2016-08-31 | 北京金双狐油气技术有限公司 | 倾斜模型层速度计算方法 |
| CN106094032B (zh) * | 2016-08-30 | 2018-06-26 | 中国石油集团东方地球物理勘探有限责任公司 | 一种构建地层速度模型的方法 |
| CN106772605B (zh) * | 2016-12-29 | 2019-04-09 | 中国石油天然气股份有限公司 | 一种均方根速度的确定方法和装置 |
| CN108663713B (zh) * | 2017-03-27 | 2021-05-25 | 中国石油化工股份有限公司 | 一种建立深度域构造模型的方法 |
| US10948617B2 (en) * | 2017-12-11 | 2021-03-16 | Saudi Arabian Oil Company | Generating a velocity model for a subsurface structure using refraction travel time tomography |
| CN108107473B (zh) * | 2017-12-12 | 2020-01-10 | 中国石油天然气集团有限公司 | 地震偏移时间剖面构造畸变现象的校正方法 |
| CN108037530B (zh) * | 2017-12-12 | 2019-11-01 | 中国石油天然气集团有限公司 | 构建速度陷阱曲线模式的方法及监控构造畸变现象的方法 |
| CN108415074B (zh) * | 2017-12-21 | 2019-12-10 | 中国石油天然气股份有限公司 | 一种缝洞储集体地震响应特征分析方法及装置 |
| CN108694744B (zh) * | 2018-05-04 | 2022-03-29 | 中国电建集团贵阳勘测设计研究院有限公司 | 一种基于多种露头的地质结构面形态综合解译方法 |
| CN108549101B (zh) * | 2018-06-28 | 2019-11-08 | 中国石油化工股份有限公司江汉油田分公司物探研究院 | 地震数据处理方法及装置 |
| CN108646294B (zh) * | 2018-06-29 | 2019-09-06 | 中海石油(中国)有限公司 | 一种复杂断裂系统下深水浊积砂岩储层高精度反演方法 |
| CN110941029B (zh) * | 2018-09-25 | 2022-08-16 | 中国石油化工股份有限公司 | 一种涉及地质戴帽的速度建模方法 |
| CN109188513B (zh) * | 2018-09-30 | 2020-03-10 | 中国石油天然气股份有限公司 | 生成深度域数据体的方法、装置及存储介质 |
| CN109975876A (zh) * | 2019-03-20 | 2019-07-05 | 中国石油化工股份有限公司 | 一种基于构造层位的井震融合速度模型的建模方法 |
| CN112147699B (zh) * | 2019-06-28 | 2022-06-21 | 中国石油化工股份有限公司 | 结合过井剖面的vsp深度域层速度模型构建方法及系统 |
| CN112180441B (zh) * | 2019-07-03 | 2024-03-26 | 中国石油天然气集团有限公司 | 转换波初始速度建模方法及装置 |
| CN112485844B (zh) * | 2019-09-12 | 2023-10-27 | 中国石油天然气股份有限公司 | 火山岩岩性岩相分布预测方法和预测装置 |
| CN112711064B (zh) * | 2019-10-24 | 2022-11-04 | 中国石油天然气股份有限公司 | 初始层速度建模方法及装置 |
| CN111060967B (zh) * | 2019-12-20 | 2022-08-23 | 恒泰艾普（北京）能源科技研究院有限公司 | 一种井控速度场建模方法 |
| CN113267810B (zh) * | 2020-02-17 | 2022-08-30 | 中国石油天然气集团有限公司 | 地震勘探全深度速度建模方法及装置 |
| CN111198400B (zh) * | 2020-02-27 | 2021-04-06 | 广州海洋地质调查局 | 一种可消除井震闭合差的速度建模方法及处理终端 |
| CN111257969B (zh) * | 2020-02-27 | 2023-02-03 | 广州海洋地质调查局 | 一种断层控制下高精度速度建模方法及处理终端 |
| CN111538080B (zh) * | 2020-05-28 | 2022-08-26 | 中国石油天然气集团有限公司 | 地震成像的方法 |
| CN111624649B (zh) * | 2020-06-05 | 2022-05-20 | 中油奥博(成都)科技有限公司 | 一种零偏移距vsp建立横向变速层速度模型的方法和装置 |
| CN113960658B (zh) * | 2020-07-20 | 2024-01-30 | 中国石油天然气股份有限公司 | 一种基于地质地震模型的测井约束速度建模方法及装置 |
| CN114114392B (zh) * | 2020-09-01 | 2023-11-28 | 中国石油天然气股份有限公司 | 层速度模型建立方法及装置 |
| CN112305589A (zh) * | 2020-09-22 | 2021-02-02 | 中国石油天然气集团有限公司 | 各向异性介质深度域成像的方法及装置 |
| CN112764104B (zh) * | 2020-12-18 | 2022-12-30 | 南方科技大学 | 一种测量方法、测量装置、终端设备及存储介质 |
| CN116381782A (zh) * | 2023-02-17 | 2023-07-04 | 中海石油(中国)有限公司深圳分公司 | 一种低幅度构造成图方法、设备及存储设备 |
| CN116559957A (zh) * | 2023-04-17 | 2023-08-08 | 广州海洋地质调查局 | 一种超深层弱反射地震数据保幅处理方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005010797A2 (en) * | 2003-07-23 | 2005-02-03 | Lee Wook B | Improved 3d veloctiy modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved-ray prestack time migration and for such migration followed by prestack depth migration |
| MX2008009584A (es) * | 2006-09-04 | 2009-04-15 | Geosystem S R L | Metodos y aparatos para la exploracion por medio de la inversion conjunta. |
| CN102565853A (zh) * | 2011-12-20 | 2012-07-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 对地质构造的速度模型进行建模的方法 |
| CN102819039B (zh) * | 2012-08-22 | 2014-11-05 | 电子科技大学 | 复杂地质条件下的层速度模型构建方法 |
| CN102967882B (zh) * | 2012-11-16 | 2015-08-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地层的层速度模型的建模方法 |
-
2014
- 2014-11-26 CN CN201410693455.6A patent/CN104360385B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104360385A (zh) | 2015-02-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104360385B (zh) | 结合地震地质认识综合建立初始深度层速度模型的方法 | |
| CN106597533B (zh) | 一种用于山前带地震资料处理的深度域速度建模方法 | |
| US20170051598A1 (en) | System For Hydraulic Fracturing Design And Optimization In Naturally Fractured Reservoirs | |
| CN106094032B (zh) | 一种构建地层速度模型的方法 | |
| US20170145793A1 (en) | Method For Modeling Stimulated Reservoir Properties Resulting From Hydraulic Fracturing In Naturally Fractured Reservoirs | |
| US10810331B2 (en) | System for predicting induced seismicity potential resulting from injection of fluids in naturally fractured reservoirs | |
| CN101634717A (zh) | 基于测井和叠前道集地震数据的精细横波阻抗求取技术 | |
| CN106405651A (zh) | 一种基于测井匹配的全波形反演初始模型构建方法 | |
| CN105899749A (zh) | 偏斜井眼和水平井眼的地层和结构解释 | |
| CN105005077A (zh) | 稀井条件下实钻井与虚拟井联合的薄层厚度预测方法 | |
| CN105549116A (zh) | 一种重建岩相古地理的方法及装置 | |
| CN107544093A (zh) | 井约束的构造解释层深度系统校正方法 | |
| CN109164487A (zh) | 一种基于模型建立平均速度场的方法及精细成图的方法 | |
| CN105911592B (zh) | 一种实钻约束的三维地震资料长波长静校正方法 | |
| Barthwal et al. | Passive seismic tomography using recorded microseismicity: Application to mining-induced seismicity | |
| Wang et al. | Multi-source data integration and multi-scale modeling framework for progressive prediction of complex geological interfaces in tunneling | |
| O'Brien et al. | Time-lapse VSP reservoir monitoring | |
| Bishop | Saskatchewan Exploratory Progress and Problems: Oil and Gas Occurrences | |
| Wilson et al. | Fracture model of the Upper Freeport coal: Marshall County West Virginia pilot ECBMR and CO2 sequestration site | |
| Singleton | Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas | |
| Zhao et al. | Re-Exploration Programs for Petroleum-Rich Sags in Rift Basins | |
| Azike | Multi-well real-time 3D structural modeling and horizontal well placement: an innovative workflow for shale gas reservoirs | |
| Ballard et al. | An integrated approach to de-risking exploration and appraisal in structurally complex fold-and-thrust belts: application to the Incahuasi field (Bolivia) | |
| CN106990433A (zh) | 一种隆起区微小侵蚀沟槽的识别方法 | |
| Palm | Time-delay interpretation of seismic refraction data in the Caledonian front, Jämtland, central Scandinavian Caledonides |
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 |
Effective date of registration: 20180207 Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp. Address before: No. 216, No. 216, Huayang Avenue, Huayang Town, Shuangliu County, Shuangliu County, Sichuan Patentee before: CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200917 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: BGP Inc., China National Petroleum Corp. Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee before: BGP Inc., China National Petroleum Corp. |
|
| TR01 | Transfer of patent right |