CN117908155A - 测井约束定量预测火山岩风化壳储层展布的方法及装置 - Google Patents
测井约束定量预测火山岩风化壳储层展布的方法及装置 Download PDFInfo
- Publication number
- CN117908155A CN117908155A CN202211237226.4A CN202211237226A CN117908155A CN 117908155 A CN117908155 A CN 117908155A CN 202211237226 A CN202211237226 A CN 202211237226A CN 117908155 A CN117908155 A CN 117908155A
- Authority
- CN
- China
- Prior art keywords
- curve
- volcanic
- volcanic rock
- density
- inversion
- 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.)
- Pending
Links
- 239000011435 rock Substances 0.000 title claims abstract description 98
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000007480 spreading Effects 0.000 title claims abstract description 12
- 238000003892 spreading Methods 0.000 title claims abstract description 12
- 238000009826 distribution Methods 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000005516 engineering process Methods 0.000 claims abstract description 14
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 18
- 239000011707 mineral Substances 0.000 claims description 18
- 238000003860 storage Methods 0.000 claims description 6
- 230000035939 shock Effects 0.000 claims description 5
- 239000003208 petroleum Substances 0.000 abstract description 4
- 206010039509 Scab Diseases 0.000 description 65
- 238000004088 simulation Methods 0.000 description 23
- 238000010586 diagram Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 10
- 230000008569 process Effects 0.000 description 9
- 230000015572 biosynthetic process Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 5
- 238000005553 drilling Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000002349 favourable effect Effects 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000002386 leaching Methods 0.000 description 3
- 230000008021 deposition Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003094 perturbing effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000009991 scouring Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例提供一种测井约束定量预测火山岩风化壳储层展布的方法及装置，属于石油勘探技术领域。所述方法包括：基于测井曲线和改进的火山岩风化壳骨架计算公式得到纵波曲线、横波曲线和密度曲线；基于井震标定后的地震数据，利用地震瞬时相位地层骨架线显示技术对火山岩风化壳的顶面和底面进行解释，得到火山岩风化壳的层位；基于所述纵波曲线、所述横波曲线、所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，确定火山岩风化壳的储层展布。本发明具有计算方法简单，能够实现火山岩风化壳储层展布精准预测的优点。
Description
技术领域
本发明涉及石油勘探技术领域，具体地涉及一种测井约束定量预测火山岩风化壳储层展布的方法、一种测井约束定量预测火山岩风化壳储层展布的装置及一种机器可读存储介质。
背景技术
火山岩风化壳型油气藏是全球范围内一类重要的油气藏类型，风化壳对火山岩油气成藏具有重要意义。与岩相和岩性相比，风化壳更易形成有利储层且易高产稳产。与未风化的火山岩相比，经过风化改造的火山岩更易形成有利储集层和高产油气藏，因此火山岩风化壳的研究对指导火山岩油气勘探具有重要意义。目前，关于火山岩风化壳的研究较少，在国内主要集中在我国西北地区古生界火山岩油气藏。通过对准噶尔盆地、三塘湖盆地火山岩风化壳孔隙成因、分布、风化壳发育模式等进行研究，确定准噶尔盆地火山岩有利储层区带主要发育在火山岩的风化体淋滤带和崩解带，火山岩风化壳能够形成良好的储层，作为一种特殊的油气藏类型其具有重要的理论和应用价值。
火山岩风化壳虽然可以作为一种优质储层，但如何通过地球物理方法对其进行表征和刻画是一项技术难题。主要存在如下问题：火山岩收到风化淋滤作用后，岩石骨架发生变化，常规地球物理参数无法与其建立明显关系；火山岩风化壳厚度较薄，一般小于地震的调谐厚度，常规地震分辨率无法对其准确识别；火山岩储层形成机理复杂、岩相和岩性变化快，有效储层预测多解性强。
发明内容
本发明实施例的目的是提供一种测井约束定量预测火山岩风化壳储层展布的方法及装置，该测井约束定量预测火山岩风化壳储层展布的方法及装置用以解决上述的现有技术中：火山岩收到风化淋滤作用后，岩石骨架发生变化，常规地球物理参数无法与其建立明显关系；火山岩风化壳厚度较薄，一般小于地震的调谐厚度，常规地震分辨率无法对其准确识别；火山岩储层形成机理复杂、岩相和岩性变化快，有效储层预测多解性强的问题。
为了实现上述目的，本发明实施例提供一种测井约束定量预测火山岩风化壳储层展布的方法，包括：
基于测井曲线和改进的火山岩风化壳骨架计算公式得到纵波曲线、横波曲线和密度曲线；
基于井震标定后的时深关系数据，利用地震瞬时相位地层骨架线显示技术对火山岩风化壳的顶面和底面进行解释，得到火山岩风化壳的层位；
基于所述纵波曲线、所述横波曲线、所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，确定火山岩风化壳的储层展布。
可选的，所述井震标定后的地震数据通过以下方式获取：
获取待预测区域的地震数据、层位解释数据和测井数据；
基于所述地震数据、所述层位解释数据和所述测井数据进行井震标定，得到所述井震标定后的时深关系数据。
可选的，所述测井数据包括：伽马测井数据、密度测井数据、声波测井数据、电阻率测井数据和阵列声波测井数据。
可选的，所述改进的火山岩风化壳骨架计算公式包括：
其中，φ为岩石孔隙度；Vci为矿物的相对体积；ρ为密度测井值；ρgas为气体密度值；ρci为矿物密度骨架；Vp为纵波测井值；Vpgas为气层纵波速度；Vpi为矿物纵波骨架；Vs为横波测井值；Vsi为矿物纵波骨架。
可选的，所述基于所述纵波曲线、所述横波曲线、所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，确定火山岩风化壳的储层展布，包括：
基于所述纵波曲线、所述横波曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演VPVS体；
基于所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演密度体。
可选的，所述方法还包括：
在所述叠前反演VPVS体上确定纵横波速度比小于第一预设值的第一区域；
在所述叠前反演密度体上确定密度小于第二预设值的第二区域；
将所述第一区域和所述第二区域相互重叠的区域确定为火山岩风化壳的分布区域，作为所述火山岩风化壳储层展布。
本发明实施例还提供一种测井约束定量预测火山岩风化壳储层展布的装置，包括：
第一确定模块，用于基于测井曲线和改进的火山岩风化壳骨架计算公式得到纵波曲线、横波曲线和密度曲线；
第二确定模块，用于基于井震标定后的时深关系数据，利用地震瞬时相位地层骨架线显示技术对火山岩风化壳的顶面和底面进行解释，得到火山岩风化壳的层位；
第三确定模块，用于基于所述纵波曲线、所述横波曲线、所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，确定火山岩风化壳的储层展布。
可选的，所述第三确定模块具体用于：
基于所述纵波曲线、所述横波曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演VPVS体；
基于所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演密度体。
可选的，所述第三确定模块还用于：
在所述叠前反演VPVS体上确定纵横波速度比小于第一预设值的第一区域；
在所述叠前反演密度体上确定密度小于第二预设值的第二区域；
将所述第一区域和所述第二区域相互重叠的区域确定为火山岩风化壳的分布区域，作为所述火山岩风化壳储层展布。
另一方面，本发明提供一种机器可读存储介质，该机器可读存储介质上存储有指令，该指令用于使得机器执行上述的测井约束定量预测火山岩风化壳储层展布的方法。
本技术方案通过改进的火山岩风化壳骨架计算公式得到精确度更高的纵波曲线、横波曲线和密度曲线；并通过地震瞬时相位地层骨架线显示技术对火山岩风化壳的顶面和底面进行解释，得到更加清晰的火山岩风化壳的层位；最后再结合所述纵波曲线、所述横波曲线、所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，计算方法简单，能够精确预测出火山岩风化壳的储层展布，为后续石油勘探提供更加准确的指导，从而提高采油率。
本发明实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施例，但并不构成对本发明实施例的限制。在附图中：
图1是本发明提供的测井约束定量预测火山岩风化壳储层展布的方法的流程示意图；
图2是本发明提供的横波曲线对比示意图；
图3是本发明提供的通过地震瞬时相位地层骨架线显示技术得到的地震瞬时相位地层骨架线的示意图；
图4是本发明提供的现有技术中地震瞬时相位地层骨架线的示意图；
图5是本发明提供的火山岩风化壳顶面的叠前反演VPVS体的分布示意图；
图6是本发明提供的火山岩风化壳顶面的叠前反演密度体的分布分布图；
图7是本发明提供的测井约束定量预测火山岩风化壳储层展布的装置的结构示意图。
附图标记说明
10-第一确定模块；20-第二确定模块；30-第三确定模块。
具体实施方式
以下结合附图对本发明实施例的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明实施例，并不用于限制本发明实施例。
术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
此外，“大致”、“基本”等用语旨在说明相关内容并不是要求绝对的精确，而是可以有一定的偏差。例如：“大致相等”并不仅仅表示绝对的相等，由于实际生产、操作过程中，难以做到绝对的“相等”，一般都存在一定的偏差。因此，除了绝对相等之外，“大致等于”还包括上述的存在一定偏差的情况。以此为例，其他情况下，除非有特别说明，“大致”、“基本”等用语均为与上述类似的含义。
在本发明的描述中，还需要说明的是，除非另有明确的规定和限定，术语“设置”、“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
图1是本发明提供的测井约束定量预测火山岩风化壳储层展布的方法的流程示意图；图2是本发明提供的横波曲线对比示意图；图3是本发明提供的通过地震瞬时相位地层骨架线显示技术得到的地震瞬时相位地层骨架线的示意图；图4是本发明提供的现有技术中地震瞬时相位地层骨架线的示意图；图5是本发明提供的火山岩风化壳顶面的叠前反演VPVS体的分布示意图；图6是本发明提供的火山岩风化壳顶面的叠前反演密度体的分布分布图；
图7是本发明提供的测井约束定量预测火山岩风化壳储层展布的装置的结构示意图。
实施例1
如图1所示，本实施例提供一种测井约束定量预测火山岩风化壳储层展布的方法，包括：
步骤101、基于测井曲线和改进的火山岩风化壳骨架计算公式得到纵波曲线、横波曲线和密度曲线；
其中，现有技术中为了了解地层，辅助实现石油勘探开发，通常采用钻井取芯、岩屑录井和地球物理勘探测井。其中，通过钻井取芯能够直观、准确的了解地层结构，但是实现成本高且效率较低；岩屑录井操作简单、快速，但是干扰因素多，深度有误差，易失真，造成看谈误差；测井是一种间接录井的首端，采用地球物理方法，连续的测定岩石的物理参数，以不同的岩石存在一定的物理性质差别，在测井曲线上有不同的变化特征为基础，利用各种测井曲线显示的特征、变化规律来划分钻井地质剖面、认识研究储层的路径方法；经济适用、收获率高，容易保存。通过测井数据能够得到测井曲线，且测井数据包括：伽马测井数据、密度测井数据、声波测井数据和电阻率测井数据。
在本实施方式中，所述改进的火山岩风化壳骨架计算公式包括：
其中，φ为岩石孔隙度；Vci为矿物的相对体积；ρ为密度测井值；ρgas为气体密度值；ρci为矿物密度骨架；Vp为纵波测井值；Vpgas为气层纵波速度；Vpi为矿物纵波骨架；Vs为横波测井值；Vsi为矿物纵波骨架；n为正整数，为地层中的矿物种类。
将测井数据中包含的各个参数代入到改进的火山岩风化壳骨架计算公式中，便能够得到多个ρ、Vp和Vs，针对ρ、Vp和Vs三类值分别进行分析，便能够得到更新后的横波曲线、纵波曲线和密度曲线，如图2所示，更新后的横波曲线、纵波曲线和密度曲线相比测井数据中的对应数据，精度更高，能够实现后续的精准预测。
以地层中存在四种矿物质为例，矿物C1，C2，C3，C4的相对体积分别为Vc1、Vc2、Vc3、Vc4；四种矿物密度骨架分别为ρc1、ρc2、ρc3、ρc4；四种矿物纵波骨架分别为Vp1、Vp2、Vp3、Vp4；四种矿物横波骨架分别为Vs1、Vs2、Vs3、Vs4，即：
1＝φ+Vc1+Vc2+Vc3+Vc4
ρ＝φρgas+ρc1Vc1+ρc2Vc2+ρc3Vc3+ρc4Vc4
Vp＝φVpgas+Vp1Vc1+Vp2Vc2+Vp3Vc3+Vp4Vc4
Vs＝Vs1Vc1+Vs2Vc2+Vs3Vc3+Vs4Vc4
步骤102、基于井震标定后的时深关系数据，利用地震瞬时相位地层骨架线显示技术对火山岩风化壳的顶面和底面进行解释，得到火山岩风化壳的层位；
其中，井震标定后的时深关系数据通过以下方式获取：首先，获取待预测区域的地震数据、层位解释数据和测井数据；基于所述地震数据、所述层位解释数据和所述测井数据进行井震标定，得到所述井震标定后的时深关系数据，该时深关系数据为井资料与地震资料建立起的对应关系，井震标定后的时深关系数据可以为该待预测区域的地震剖面，通过地震数据、层位解释数据和测井数据进行井震标定后分析得到。更具体地，要获得地震剖面，首先需建立待预测区域的三维地质模型，可以采用石文地质建模系统基于测井数据、地震数据和层位解释数据进行建模，在对三维地质模型进行剖切，便能够得到地震瞬时相位属性剖面。但是，由于得到的地震瞬时相位属性剖面分辨率较低，剖面信息比较杂乱，因此，通过对显示色标的调节，可以把火山岩风化壳的顶面和底面通过瞬时相位地层骨架线的显示方式展现出来，得到清晰的火山岩风化壳的层位，如图3所示，现有技术中得到的火山岩风化壳的层位，如图4所示，相比之下，通过本技术方案到得的火山岩风化壳的层位更加清晰。
在另一种实施方式中，在建立三维地质模型前，对获取的测井数据进行一致性处理；由于测井数据存在采集年代跨度大，测井仪器信号众多，刻度标准不同意，操作方式不一致等问题，为了消除不同时间、不同仪器所测量的测井数据之间存在的不同的系统误差，在利用数据建模前，通常对测井数据进行一致性处理，具体包括：采用均值-方差法，对测井曲线进行均值和方差的分析，将不同时间、不同仪器所测量的测井数据进行相对统一。
更具体地，利用地震瞬时相位地层骨架线显示技术对火山岩风化壳的顶面和底面进行解释，包括：
通过对野外火山岩风化壳地质露头进行实地考察，对露头进行地质模型的建立，通过分析火山岩风化壳顶面为一个不规则高低起伏界面。而由于火山受到雨水的渗透和冲刷的作用，会在火山岩风化壳内部形成一个与水平面近似平行的一个连续界面。
获取较可靠的叠后地震数据(要达到二类地震品质以上)和相关测井数据包括：伽马测井数据、密度测井数据、声波测井数据、电阻率测井数据和阵列声波测井数据等并进行井震标定。
在常规解释软件或反演软件(本次使用的SMI软件)中进行瞬时相位地震属性提取。
调整瞬时相位属性的色标值域和显示方式，可以把地层的结构的骨架结构呈现出来。常规的地震瞬时相位属性剖面分辨率较低，剖面信息比较杂乱。通过对显示色标的调节，并与各井精确标定的结果相结合可以把火山岩风化壳的顶面和底面通过瞬时相位地层骨架线的显示方式展现出来。
把火山岩风化壳的骨架线位置调整与钻井火山岩风化壳位置一致，可以得到区域内火山岩风化壳顶面和平面的分布规律。
把瞬时相位骨架线剖面上，找一个风化壳上部一个比较稳定的层位(本次选择的是一个泥岩的顶面)，对泥岩顶面进行层拉平处理，在层拉平剖面的火山岩风化壳内部找到与层拉平层位近似平行的骨架线范围，认为应该是火山岩风化壳底面的位置，与实钻井吻合情况较好，得到火山岩风化壳的层位。
步骤103、基于所述纵波曲线、所述横波曲线、所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，确定火山岩风化壳的储层展布。
具体包括：基于所述纵波曲线、所述横波曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演VPVS体；基于所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演密度体。得到的叠前反演VPVS体和叠前反演密度体基于同一个火山岩风化壳的层位得到，因此，叠前反演VPVS体和叠前反演密度体的形状和轮廓分布相同，两者的不同点在于：叠前反演VPVS体展示的是纵横波速度比在叠前反演VPVS体上的分布情况，叠前反演密度体展示的是密度在叠前反演VPVS体上的分布情况。
最后，在所述叠前反演VPVS体上确定纵横波速度比小于第一预设值的第一区域；在所述叠前反演密度体上确定密度小于第二预设值的第二区域；将所述第一区域和所述第二区域相互重叠的区域确定为火山岩风化壳的分布区域，作为所述火山岩风化壳储层展布。
具体地，第一预设值可根据不同地质环境下气层和干层的相关数据和岩石物理交会图分析得出，由于地层中存在气层，因此，横波传播速度会减慢，所以气层呈现的是纵横波速度比低的特征，根据历史数据进行分析，便能够确定出第一预设值的数据。(现有技术中区分气层和干层的分界值为1.74，即纵横波速度比大于1.74，该区域点存在气层的概率低，纵横波速度比小于等于1.74，该区域点存在气层的概率高，且该纵横波速度比低越小，该点位于气层的概率越高)。同理由岩石物理关系可知高孔隙度储层的地球物理特征为低密度值，储层的孔隙度高，所以储层的物性好，含气的概率变大，因此，设定第二预设值可以确认出第二区域。(叠前反演密度体上密度大于第二预设值的点位，孔隙度越小，存在气体的概率就越低)。因此，采用不同的参数和不同的预测方法分别在所述叠前反演VPVS体上确定纵横波速度比小于第一预设值的第一区域(即含气量高的区域)，如图5所示；在所述叠前反演密度体上确定密度小于第二预设值的第二区域(即含气量高的区域)，如图6所示，排除断层区域带来的影响，最后，将第一区域和第二区域中相互重叠的区域确定为火山岩风化壳的分布区域，作为所述火山岩风化壳储层展布。
实施例2
在本实施方式中采用叠前波形指示反演进行分析，能够包含更加丰富的地下地质、岩性和油气信息，叠前反演包含横波信息，既能够预测砂岩，又能够把有利储层从砂岩中区分开来，当纵波阻抗不能反映岩性及有力储层，或者地震剖面而含油气层为弱反射时，采用叠前反演具有较大优势。
叠前弹性参数反演技术因能获取地层纵波阻抗、横波阻抗以及泊松比等多项弹性参数，进而依据岩石物理响应对储层和流体进行有效识别。在高分辨率叠前地震反演中，以近年发展起来的叠前地质统计学反演技术为典型代表，通过地质统计学理论与叠前同时反演理论的有机结合，反演结果对高于地震频带的高频部分进行了较为合理的补充。然而，仔细考察地质统计学反演的方法原理和技术流程不难发现，其反演结果中高于地震频带的高频成分主要来源于地质统计学模拟，高频成分的随机性受区域化变量的二阶平稳性、样本空间的数量和分布影响较大，实际应用中往往表现为无法拟合出吻合度较高的变差函数和先验概率密度函数。对于具有薄互层反射系数结构的地层而言，随机性较强的高频成分在反演的多个实现中容易出现差异较大的多种弹性参数结构组合，直接影响薄储层及其含油性预测的确定性和精度。
波形指示叠前地震反演方法解决了叠前随机反演高频成分随机性强的难题，利用叠前地震波形相似性代替传统变差函数优选随机模拟样本，由于样本共性结构的频带宽度高于地震频带，使随机模拟的中高频部分确定性大大增强，获得中低频可靠，高频逐步确定的叠前弹性参数反演结果。该技术在川中致密砂岩薄储层及含油性预测的应用中取得了良好的应用效果，精细刻画了薄储层及油层的空间分布，为油藏精细描述和高效开发部署提供了更为可靠的依据。
波形指示叠前方法原理：
基于褶积理论的正演研究表明，地震波的调谐特征与反射系数结构具有密切的关系，反射系数的垂向分布，包括反射系数的间距、大小、个数决定了地震波形的调谐样式。在相似的地质条件下，反射波波形与反射系数结构可形成良好的匹配关系，可以根据波形的调谐特征优化反射系数的分布。以纵波阻抗和纵波反射数据为例，如果两口井目标层所处位置的地震波形相似，表明这两口井处于大致相同的沉积环境，虽然纵波阻抗的高频成分可能来自不同的沉积微相、差异较大，但中低频具有共性，且共性频带范围大幅度超出了地震有效频带。地震波形相似的情况下纵波阻抗的共性结构，通过对具有相似地震波形的两口井的纵波阻抗进行逐级滤波，70％以上相关系数的共性结构的频带宽度高达5～200Hz，远高出原始地震数据的5～80Hz频宽。另外，更高频率成分的取值也被约束在一定范围内。实例分析表明，地震波形对地层弹性参数的结构与分布具有重要的指示意义。为了充分利用上述地震波形对地层弹性参数结构的指示意义，在随机模拟过程中利用地震波形相似性代替变差函数对模拟样本进行优选。基于传统变差函数的随机模拟是根据变程优选统计样本进行模拟，考虑的是空间变异程度，和距离相关；波形指示模拟主要是根据波形相似性优选统计样本，将预测道地震波形与所有已知井旁道地震波形对比，优选出最相似的若干井样本，进行随机模拟建立初始模型。对于叠前弹性参数的随机模拟而言，利用纵波反射地震波形指示纵波阻抗的模拟，利用横波反射地震波形(多波勘探或叠前同时反演获得)指示横波阻抗的模拟，这也是叠前弹性参数随机反演的关键环节。
从波形指示模拟的原理不难看出，模拟过程充分利用了空间密集分布的地震数据并体现了相控模拟思想，其频率成分是一个由低到高逐步确定的过程，高频成分的整体确定性相比传统随机模拟得到大幅提高。此外，地震波形指示模拟建立初始模型的过程并不采用序贯的方式进行，即待模拟点的模拟数据并不参与下一个未知点的模拟过程。这样做使得模拟过程更加尊从地震波形相似的样本优选原则，即更加体现相控的思想；其次，每个待模拟点的模拟都从实际优选样本出发进行无偏最优估计，避免高频随机误差的累积，使模拟结果更加稳定可靠。
波形指示叠前地震反演流程：
波形指示叠前地震反演是在波形指示弹性参数模拟的基础上，结合叠前同时反演而形成的一种叠前高分辨率随机反演技术，最终获得既符合原始地震记录，又包含确定性较强的高频成分的宽频弹性参数结果。整个反演过程分为以下4个主要组成部分。
叠前同时反演获取地震频带范围的纵、横波阻抗绝对值。叠前同时反演根据Zoeppritz方程或其近似式，运用井约束稀疏脉冲算法保证了获取的纵、横波阻抗值在地震频带范围内与实际地震信息的吻合。该反演技术唯一性好、抗噪能力强、对井依赖小、结果稳定性好，是运用成熟的反演技术。作为波形指示叠前地震反演技术的关键环节，为最终反演结果提供尊重原始地震数据的确定性成分，也是波形指示叠前地震反演残差的主要来源。
波形指示模拟建立初始模型。将叠前同时反演获得的纵、横波阻抗与平均AVA子波褶积得到纵、横波合成地震数据，在构造框架下，按纵、横波地震波形特征分别对已知井进行分析，优选与待模拟道波形关联度高的井样本进行随机模拟，分别建立纵、横波阻抗初始模型，并统计测井样本的纵、横波阻抗的概率分布作为先验概率分布。为了避免大距离范围的样本优选误差，优选过程增加考虑最大空间距离，在已知井中利用波形相似性和空间距离双变量优选中低频结构相似的井作为空间估值样本。
将纵、横波阻抗初始模型与叠前同时反演的纵、横波阻抗进行匹配滤波，计算得到似然函数。根据似然函数的定义，此处通过匹配滤波计算得到的似然函数描述了某一个空间位置取值为某一特定值的概率。匹配滤波将确定性信息(源自叠前同时反演)和逐步确定性信息(源自波形指示模拟)进行融合之后产生了一个新的概率分布空间。
在贝叶斯框架下联合似然函数和先验概率得到后验概率密度分布，对其进Metropolis-Has-tings采样作为目标函数。不断扰动模型参数，使后验概率密度值最大，此时的解作为可行随机实现，取多次可行实现的均值作为期望值输出，得到宽频纵、横波阻抗反演结果。
实施例3
如图7所示，本发明实施方式还提供一种测井约束定量预测火山岩风化壳储层展布的装置，包括：
第一确定模块10，用于基于测井曲线和改进的火山岩风化壳骨架计算公式得到纵波曲线、横波曲线和密度曲线；
第二确定模块20，用于基于井震标定后的地震数据，利用地震瞬时相位地层骨架线显示技术对火山岩风化壳的顶面和底面进行解释，得到火山岩风化壳的层位；
第三确定模块30，用于基于所述纵波曲线、所述横波曲线、所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，确定火山岩风化壳的储层展布。
进一步地，所述第三确定模块具体用于：
基于所述纵波曲线、所述横波曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演VPVS体；
基于所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演密度体。
进一步地，所述第三确定模块还用于：
在所述叠前反演VPVS体上确定纵横波速度比小于第一预设值的第一区域；
在所述叠前反演密度体上确定密度小于第二预设值的第二区域；
将所述第一区域和所述第二区域相互重叠的区域确定为火山岩风化壳的分布区域，作为所述火山岩风化壳储层展布。
实施例4
本发明还提供一种机器可读存储介质，该机器可读存储介质上存储有指令，该指令用于使得机器执行本申请上述的测井约束定量预测火山岩风化壳储层展布的方法。
以上结合附图详细描述了本发明实施例的可选实施方式，但是，本发明实施例并不限于上述实施方式中的具体细节，在本发明实施例的技术构思范围内，可以对本发明实施例的技术方案进行多种简单变型，这些简单变型均属于本发明实施例的保护范围。
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明实施例对各种可能的组合方式不再另行说明。
本领域技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序存储在一个存储介质中，包括若干指令用以使得单片机、芯片或处理器(processor)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
此外，本发明实施例的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明实施例的思想，其同样应当视为本发明实施例所公开的内容。
Claims (10)
1.一种测井约束定量预测火山岩风化壳储层展布的方法，其特征在于，包括：
基于测井曲线和改进的火山岩风化壳骨架计算公式得到纵波曲线、横波曲线和密度曲线；
基于井震标定后的时深关系数据，利用地震瞬时相位地层骨架线显示技术对火山岩风化壳的顶面和底面进行解释，得到火山岩风化壳的层位；
基于所述纵波曲线、所述横波曲线、所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，确定火山岩风化壳的储层展布。
2.根据权利要求1所述的方法，其特征在于，所述井震标定后的地震数据通过以下方式获取：
获取待预测区域的地震数据、层位解释数据和测井数据；
基于所述地震数据、所述层位解释数据和所述测井数据进行井震标定，得到所述井震标定后的时深关系数据。
3.根据权利要求2所述的方法，其特征在于，所述测井数据包括：伽马测井数据、密度测井数据、声波测井数据、电阻率测井数据和阵列声波测井数据。
4.根据权利要求1所述的方法，其特征在于，所述改进的火山岩风化壳骨架计算公式包括：
其中，φ为岩石孔隙度；Vci为矿物的相对体积；ρ为密度测井值；ρgas为气体密度值；ρci为矿物密度骨架；Vp为纵波测井值；Vpgas为气层纵波速度；Vpi为矿物纵波骨架；Vs为横波测井值；Vsi为矿物纵波骨架。
5.根据权利要求1所述的方法，其特征在于，所述基于所述纵波曲线、所述横波曲线、所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，确定火山岩风化壳的储层展布，包括：
基于所述纵波曲线、所述横波曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演VPVS体；
基于所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演密度体。
6.根据权利要求5所述的方法，其特征在于，所述方法还包括：
在所述叠前反演VPVS体上确定纵横波速度比小于第一预设值的第一区域；
在所述叠前反演密度体上确定密度小于第二预设值的第二区域；
将所述第一区域和所述第二区域相互重叠的区域确定为火山岩风化壳的分布区域，作为所述火山岩风化壳储层展布。
7.一种测井约束定量预测火山岩风化壳储层展布的装置，其特征在于，包括：
第一确定模块，用于基于测井曲线和改进的火山岩风化壳骨架计算公式得到纵波曲线、横波曲线和密度曲线；
第二确定模块，用于基于井震标定后的时深关系数据，利用地震瞬时相位地层骨架线显示技术对火山岩风化壳的顶面和底面进行解释，得到火山岩风化壳的层位；
第三确定模块，用于基于所述纵波曲线、所述横波曲线、所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，确定火山岩风化壳的储层展布。
8.根据权利要求7所述的装置，其特征在于，所述第三确定模块具体用于：
基于所述纵波曲线、所述横波曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演VPVS体；
基于所述密度曲线和所述火山岩风化壳的层位进行叠前波形指示反演，得到叠前反演密度体。
9.根据权利要求8所述的装置，其特征在于，所述第三确定模块还用于：
在所述叠前反演VPVS体上确定纵横波速度比小于第一预设值的第一区域；
在所述叠前反演密度体上确定密度小于第二预设值的第二区域；
将所述第一区域和所述第二区域相互重叠的区域确定为火山岩风化壳的分布区域，作为所述火山岩风化壳储层展布。
10.一种机器可读存储介质，该机器可读存储介质上存储有指令，该指令用于使得机器执行权利要求1-6中任一项所述的测井约束定量预测火山岩风化壳储层展布的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211237226.4A CN117908155A (zh) | 2022-10-11 | 2022-10-11 | 测井约束定量预测火山岩风化壳储层展布的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211237226.4A CN117908155A (zh) | 2022-10-11 | 2022-10-11 | 测井约束定量预测火山岩风化壳储层展布的方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117908155A true CN117908155A (zh) | 2024-04-19 |
Family
ID=90695106
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211237226.4A Pending CN117908155A (zh) | 2022-10-11 | 2022-10-11 | 测井约束定量预测火山岩风化壳储层展布的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117908155A (zh) |
-
2022
- 2022-10-11 CN CN202211237226.4A patent/CN117908155A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113759425B (zh) | 井震联合评价深层古岩溶储层充填特征的方法与系统 | |
| CN113759424B (zh) | 基于频谱分解和机器学习的岩溶储层充填分析方法和系统 | |
| US7400978B2 (en) | Method for updating a geologic model by seismic data | |
| Dolan et al. | The broad-band fractal nature of heterogeneity in the upper crust from petrophysical logs | |
| CN108363100B (zh) | 基于封存单元与粗糙集的煤层气甜点区地震地质识别方法 | |
| US20080162093A1 (en) | Method of constructing a geological model of a subsoil formation constrained by seismic data | |
| CN108680951A (zh) | 一种基于地震信息判断煤层气富集沉积控制作用的方法 | |
| CN111273348A (zh) | 基于更新概率比率恒定理论的多点地质统计叠前反演方法 | |
| CN105759312A (zh) | 一种储层特征法井震标定方法 | |
| WO2012139082A1 (en) | Event selection in the image domain | |
| CN103257363B (zh) | 一种探测地下裂缝性储层中裂缝倾角的方法 | |
| CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
| CN111856566B (zh) | 湖相滩坝砂体中薄储层预测方法及装置 | |
| CN110095811B (zh) | 膏岩层速度模型构建处理方法及装置 | |
| CN114994758B (zh) | 碳酸盐岩断控储层的波阻抗提取与结构表征方法和系统 | |
| CN113050157A (zh) | 一种基于露头资料的碳酸盐岩地震储层反演方法及系统 | |
| CN110646850B (zh) | 隔夹层地震预测方法及装置 | |
| Eikrem et al. | Bayesian estimation of reservoir properties—effects of uncertainty quantification of 4D seismic data | |
| CN111077578B (zh) | 岩层分布预测方法和装置 | |
| CN117908155A (zh) | 测井约束定量预测火山岩风化壳储层展布的方法及装置 | |
| CN113806674A (zh) | 古河道纵向尺度的量化方法、装置、电子设备及存储介质 | |
| CN105989240A (zh) | 砂层预测方法及装置 | |
| Torres-Verdin et al. | A Comparison between Geostatistical Inversion and Conventional Geostatistical-simulation Practices for Reservoir Delineation | |
| US20230273332A1 (en) | Generating low frequency models for seismic waveform inversion in formation regions with limited control wells | |
| CN110941029B (zh) | 一种涉及地质戴帽的速度建模方法 |
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 |