CN105425299B - 确定地层裂缝分布的方法和装置 - Google Patents
确定地层裂缝分布的方法和装置 Download PDFInfo
- Publication number
- CN105425299B CN105425299B CN201510788560.2A CN201510788560A CN105425299B CN 105425299 B CN105425299 B CN 105425299B CN 201510788560 A CN201510788560 A CN 201510788560A CN 105425299 B CN105425299 B CN 105425299B
- Authority
- CN
- China
- Prior art keywords
- interval
- represent
- seismic data
- data cube
- sample value
- 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 50
- 230000015572 biosynthetic process Effects 0.000 title claims description 32
- 230000001427 coherent effect Effects 0.000 claims abstract description 25
- 238000005070 sampling Methods 0.000 claims abstract description 10
- 238000012937 correction Methods 0.000 claims description 46
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000012546 transfer Methods 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 4
- 239000004744 fabric Substances 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 230000009466 transformation Effects 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000036039 immunity Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000013001 point bending Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种确定地层裂缝分布的方法和装置，包括以下步骤：获取工区的地震数据体；将地震数据体目的层段样点值的空间位置，由倾斜层段校正至水平层段；对在水平层段的地震数据体进行相干变换，得到相干属性数据体；将相干属性数据体目的层段样点值的空间位置，由水平层段反向校正至倾斜层段，得到工区的地层裂缝分布。在本发明实施例中，在求取相干属性前，将目的层段样点值的空间位置由倾斜层段校正至水平层段后提取相干属性，有效地消除了大倾角地层对裂缝预测结果的干扰；再将相干属性数据体中的样点值的空间位置，由水平层段反向校正至倾斜层段，从而有效地提高了大倾角地层裂缝空间展布形态的预测精度，保证裂缝预测结果的准确性。
Description
技术领域
本发明涉及石油勘探技术领域，特别涉及一种确定地层裂缝分布的方法和装置。
背景技术
近年来，随着油气田勘探程度的深入，构造简单并且容易开采的油气藏逐渐减少，人们不得不去开采那些构造复杂的油气藏。但是这些构造复杂的油气藏经历了多次的构造运动，地层变形严重，许多地层倾角变陡，并且产生了大量裂缝。对于油气藏而言，这些地层裂缝不仅是重要的油气储集空间和渗流通道，同时还是优质储层发育条带的主控因素之一，所以准确预测地层裂缝的空间展布是油气勘探中一项重要的研究内容。
在复杂的地下构造中，由于地层裂缝的空间展布规律性差，因此，对其进行准确的预测难度较大。常用的地震裂缝预测技术分为叠后和叠前两大类。叠后地震裂缝预测技术主要包括相干、曲率及蚂蚁追踪等技术，其中，相干技术的基本原理是将地震道以多种形式进行组合来检测地震数据空间的不连续性，从而反映出断层以及裂缝的展布形态；曲率是表征地震道层面上某一点变形弯曲的程度，层面变形弯曲越厉害曲率值就越大；蚂蚁追踪技术是基于蚂蚁算法根据地震振幅及相位之间的差异，沿着可能的断层和裂缝向前移动，直至将断层和裂缝完全刻画出来。叠前地震裂缝预测技术以方位各向异性分析技术为代表，通过研究不同方位角地震属性的变化情况表征裂缝的空间展布形态。在实际应用中，相对于其他技术而言，由于叠后地震裂缝预测技术中的相干技术算法稳定，运行高效且裂缝空间展布趋势清晰，因此，通常采用这种方法进行裂缝预测。
自1995年Bahorich等首次提出基于归一化互相关的相干算法以来，相干技术得到了长远的发展，已经成为地下断层和裂缝描述的重要技术手段。1998年Marfurt等人提出了基于多道相似性测量的第二代相干算法，通过构建协方差矩阵计算地震道之间的相似性求取相干体。该算法考虑了地层倾角和方位角信息，增加了抗噪性和稳定性，但是横向分辨率较低。1999年Gersztenkorn等人又提出了基于本征结构的第三代相干算法，通过矩阵的特征结构计算相干性。这种算法进一步提高了横向分辨率，但是对大倾角地层不敏感。随后国内学者们对相干算法进行了诸多探索，相继提出了基于小波变换的相干算法、自适应时窗相干体计算技术等，用来提高方法的抗噪性和分辨率，但始终未能解决准确预测大倾角地层裂缝展布的问题。
针对上述如何利用常规相干技术准确预测大倾角地层裂缝空间展布形态的问题，目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种确定地层裂缝分布的方法，以解决现有技术中无法利用常规相干技术准确预测大倾角地层裂缝空间展布形态的问题。
本发明实施例提供了一种确定地层裂缝分布的方法，该方法包括：获取工区的地震数据体；将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段，得到所述目的层段样点值在水平层段的地震数据体；对所述在水平层段的地震数据体进行相干变换，得到表征裂缝分布的相干属性数据体；将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段，得到所述工区的地层裂缝分布。
在一个实施例中，在对所述在水平层段的地震数据体进行相干变换之前，所述方法还包括：通过中值滤波方法和/或均值滤波方法对所述水平层段的地震数据体进行去噪处理。
在一个实施例中，按照以下公式将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段：
t′(xi,yj)＝t(xi,yj)-TR(xi,yj)+TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR(xi,yj)表示倾斜层段的解释层位，TR′(xi,yj)表示对倾斜层段的解释层位TR(xi,yj)进行校正后得到的水平层段的解释层位。
在一个实施例中，按照以下公式对在水平层段的地震数据体进行相干变换：
其中，Ec表示由地震数据体样点值组成的矩阵的相干属性值，λn表示由地震数据体样点值得到的协方差矩阵的第n个特征值，N表示进行相干计算的样点值的个数，λ1表示N个特征值中最大的特征值。
在一个实施例中，按照以下公式将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段：
t′(xi,yj)＝t(xi,yj)+TR(xi,yj)-TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示反向校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示反向校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR′(xi,yj)表示水平层段的解释层位，TR(xi,yj)表示对水平层段的解释层位TR′(xi,yj)进行反向校正后得到的倾斜层段的解释层位。
本发明实施例还提供了一种确定地层裂缝分布的装置，包括：地震数据体获取单元，用于获取工区的地震数据体；倾斜层段校正单元，用于将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段，得到所述目的层段样点值在水平层段的地震数据体；相干变换单元，用于对所述在水平层段的地震数据体进行相干变换，得到表征裂缝分布的相干属性数据体；水平层段反向校正单元，用于将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段，得到所述工区的地层裂缝分布。
在一个实施例中，还包括：去噪处理单元，用于在对所述在水平层段的地震数据体进行相干变换之前，通过中值滤波方法和/或均值滤波方法对所述水平层段的地震数据体进行去噪处理。
在一个实施例中，倾斜层段校正单元具体用于按照以下公式将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段：
t′(xi,yj)＝t(xi,yj)-TR(xi,yj)+TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR(xi,yj)表示倾斜层段的解释层位，TR′(xi,yj)表示对倾斜层段的解释层位TR(xi,yj)进行校正后得到的水平层段的解释层位。
在一个实施例中，相干变换单元具体用于按照以下公式对在水平层段的地震数据体进行相干变换：
其中，Ec表示由地震数据体样点值组成的矩阵的相干属性值，λn表示由地震数据体样点值得到的协方差矩阵的第n个特征值，N表示进行相干计算的样值点的个数，λ1表示N个特征值中最大的特征值。
在一个实施例中，水平层段反向校正单元具体用于按照以下公式将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段：
t′(xi,yj)＝t(xi,yj)+TR(xi,yj)-TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示反向校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示反向校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR′(xi,yj)表示水平层段的解释层位，TR(xi,yj)表示对水平层段的解释层位TR′(xi,yj)进行反向校正后得到的倾斜层段的解释层位。
在上述实施例中，在求取相干属性前，将地震数据体中的目的层段样点值的空间位置由倾斜层段校正至水平层段后提取相干属性，从而有效消除了大倾角地层对裂缝预测结果的干扰；再将相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至倾斜层段，从而有效地提高了大倾角地层裂缝空间展布形态的预测精度，保证了裂缝预测结果的准确性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1是本发明实施例的确定地层裂缝分布的方法的流程图；
图2是本发明实施例的确定地层裂缝分布的方法对四川盆地Z段大倾角地层裂缝预测的流程图；
图3是本发明实施例的原始地震剖面示意图；
图4是本发明实施例的校正后的地震剖面示意图；
图5是本发明实施例的四川盆地Z段大倾角地层裂缝预测结果示意图；
图6是本发明实施例的基于常规相干算法得到的四川盆地Z段大倾角地层裂缝预测结果示意图；
图7是本发明实施例的确定地层裂缝分布的装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
考虑到现有的相干算法不能准确反映大倾角地层裂缝空间展布形态，发明人考虑到可以将从原始地震数据直接提取相干属性，调整为先水平校正原始地震数据大倾角目的层，然后提取水平校正后的大倾角目的层的相干属性，从而消除大倾角的存在对地层裂缝空间展布的不利影响，以解决现有技术中利用常规相干技术无法准确预测大倾角地层裂缝空间展布形态的问题。具体的，在本例中，提供了一种确定地层裂缝分布的方法，如图1所示，包括以下步骤：
步骤101：获取工区的地震数据体；
具体地，获取的地震数据体是经过保幅处理的地震数据体。这主要是考虑到在地震勘探中，地震波经过激发在地下传播接收，经历了地表近地表的低降速带的衰减、传播过程中的波前扩散和地层介质的吸收、地质界面的反射等，因而勘探得到的地震波是经过数项“改造”后得到的。因此，在处理地震资料时，就需要先将这些“改造”消除，同时保持地震波自身的振幅特性不变，这也就是所谓的保幅处理。
步骤102：将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段，得到所述目的层段样点值在水平层段的地震数据体；
为了消除大倾角地层对裂缝预测结果的干扰，在本实施例中将所采集的地震数据体中目的层段样点值的空间位置，由倾斜层段校正至水平层段，从而构成新的地震数据体。具体的，可以按照如下公式计算得到新的地震数据体位置：
t′(xi,yj)＝t(xi,yj)-TR(xi,yj)+TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR(xi,yj)表示倾斜层段的解释层位，TR′(xi,yj)表示对倾斜层段的解释层位TR(xi,yj)进行校正后得到的水平层段的解释层位。
在对原始地震数据体进行水平校正处理的过程中，会产生一定的随机噪声。为了提高水平校正处理后的地震数据体的相干属性的预测精度，可以对水平校正后的地震数据体进行去噪处理。
在本实施例中，可以通过中值滤波和/或均值滤波等方法对水平层段的地震数据体进行去噪处理。
具体地，可以在校正处理后的地震数据体中定义一个时窗，假设该时窗中包含2j+1个样本点，那么分别采用中值滤波法和均值滤波法后时窗中间第j+1个样本点的值为：
1)采用中值滤波法后，时窗中间的第j+1个样本点的值为：
其中，假设Si表示第j+1个样本点的值，则Si-j表示时窗内第一个样本点的值，Si-j+1表示时窗内第二个样本点的值，以此类推。
2)采用均值滤波法后，时窗中间的第j+1个样本点的值为：
步骤103：对所述在水平层段的地震数据体进行相干变换，得到表征裂缝分布的相干属性数据体；
为了正确预测裂缝的空间展布形态，可以求取去噪处理后的地震数据体的相干属性。具体的，在本例中，可以采用如下步骤求取地震数据体的相干属性：
步骤1：假设在一个分析窗口中有N道地震数据，K个采样点，则三维地震数据体可以表示为矩阵D：
其中，dkn表示第n道的第k个样点值。
步骤2：按照以下公式求取该三维地震数据体矩阵D的协方差矩阵：
其中，dk＝[dk1,dk2,…,dkN],1≤k≤K，表示矩阵D中的第k个采样点的集合。
步骤3：通过地震数据体的协方差矩阵计算得到由地震数据体样点值所组成的矩阵的相干属性值Ec，具体地，可以按照以下公式计算相干属性值：
其中，λn表示由地震数据体样点值得到的协方差矩阵的第n个特征值，N表示地震数据的地震道道数，也表示特征值的个数，λ1表示N个特征值中最大的特征值。
在本实施例中，不同时窗所对应的相干属性值Ec是不同的，相干属性值的大小与裂缝的相关性有直接的关系，一般情况下，相干属性值Ec越小裂缝的相干性越差，在裂缝预测中物理表现为相干属性值Ec越小裂缝越发育。
步骤104：将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段，得到所述工区的地层裂缝分布。
在得到不同时窗所对应的裂缝的发育程度后，通过反向校正所述新的地震数据体至所述原始大倾角目的层，就可以将相干数据体反校正到其原始空间位置，从而可以得到最终的原始地震数据体的大倾角地层裂缝的空间展布形态。
具体的，在本实施例中，可以按照以下公式反向校正至原始大倾角倾斜层段：
t′(xi,yj)＝t(xi,yj)+TR(xi,yj)-TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示反向校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示反向校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR′(xi,yj)表示水平层段的解释层位，TR(xi,yj)表示对水平层段的解释层位TR′(xi,yj)进行反向校正后得到的倾斜层段的解释层位。
根据所得到的反向校正后的地震数据体样点在空间的位置，以及该位置所对应的时窗的相干属性值，可以得到原始地震数据体的大倾角地层裂缝的空间展布形态。
下面以一个确定地层裂缝分布方法的具体应用实例来进行说明，在本实施例中，通过上述确定地层裂缝分布的方法对四川盆地Z段大倾角地层进行裂缝预测。然而值得注意的是，该具体实施例仅是为了更好地说明本发明，并不构成对本发明的不当限定。
如图2所示，通过上述确定地层裂缝分布的方法对四川盆地Z段大倾角地层进行裂缝预测可以包括以下步骤：
步骤201：获取四川盆地Z段大倾角地层的地震数据体；
步骤202：将该地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段，得到目的层段样点值在水平层段的地震数据体；
令TR′(xi,yj)＝3000ms，按照以下公式求得校正后平面坐标为(xi,yj)的地震道数据对应的深度值：
t′(xi,yj)＝t(xi,yj)-TR(xi,yj)+TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR(xi,yj)表示倾斜层段的解释层位，TR′(xi,yj)表示对倾斜层段的解释层位TR(xi,yj)进行校正后得到的水平层段的解释层位。
如图3所示的原始地震剖面，采用该步骤202进行校正后可以得到如图4所示的地震剖面。从图3可以看出，四川盆地Z段大倾角地层裂缝沉积厚度稳定、地层倾角大，并且剖面图右侧地层倾角超过45度，而经过水平校正后得到的图4中的地震剖面的大倾角处于趋于平缓的状态。
步骤203：对该目的层段样点值在水平层段的地震数据体进行去噪处理；
步骤204：对去噪处理后的地震数据体进行相干变换，得到表征裂缝分布的相干属性数据体；
在该实施例中所取的分析窗口中有N＝9道地震数据，K＝9个采样点，按照以下步骤求得地震数据体的相干属性：
步骤1：根据地震数据的道数N，采样点的个数K，可以得到三维地震数据体矩阵D：
其中，dkn表示第n道的第k个样点值。
步骤2：按照以下公式求取该三维地震数据体矩阵D的协方差矩阵C：
其中，dk＝[dk1,dk2,…,dk9],1≤k≤9，表示矩阵D中的第k个采样点的集合。
步骤3：通过地震数据体的协方差矩阵计算得到由地震数据体样点值所组成的矩阵的相干属性值Ec，具体地，可以按照以下公式计算相干属性值：
其中，λn表示由地震数据体样点值得到的协方差矩阵的第n个特征值，N表示地震数据的地震道道数，也表示特征值的个数，λ1表示N个特征值中最大的特征值。
步骤205：反向校正新的地震数据体至原始大倾角目的层，得到大倾角地层裂缝的空间展布形态；
令TR′(xi,yj)＝3000ms，可以按照以下公式反向校正所述新的地震数据体至所述原始大倾角目的层：
t′(xi,yj)＝t(xi,yj)+TR(xi,yj)-TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示反向校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示反向校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR′(xi,yj)表示水平层段的解释层位，TR(xi,yj)表示对水平层段的解释层位TR′(xi,yj)进行反向校正后得到的倾斜层段的解释层位。
步骤206：利用最终得到的相干属性预测大倾角目的层段裂缝的空间展布形态。
如图5所示为采用本发明的确定地层裂缝分布的方法预测出的四川盆地Z段大倾角地层裂缝示意图，图6为基于现有的常规相干算法得到的裂缝预测结果示意图。对比图5和图6可以看出，在非大倾角附近，采用两种不同方法所得到的裂缝预测结果基本相同，然而，在大倾角附近，采用本发明方法所得到的裂缝预测结果可以较好的显示裂缝之间的细节，并且不同裂缝之间的重叠较少。由此可见，通过本例所提供的针对大倾角地层的裂缝预测方法可以有效消除大倾角地层对裂缝预测结果的影响，同时保留真实的裂缝信息，最终可以得到准确可靠的裂缝预测结果。
基于同一发明构思，本发明实施例中还提供了一种确定地层裂缝分布的装置，如下面的实施例所述。由于确定地层裂缝分布的装置解决问题的原理与确定地层裂缝分布的方法相似，因此确定地层裂缝分布的装置的实施可以参见确定地层裂缝分布的方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图7是本发明实施例的确定地层裂缝分布的装置的一种结构框图，如图7所示，包括：地震数据体获取单元301、倾斜层段校正单元302、相干变换单元303、水平层段反向校正单元304，下面对该结构进行说明。
地震数据体获取单元301，用于获取工区的地震数据体；
倾斜层段校正单元302，用于将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段，得到所述目的层段样点值在水平层段的地震数据体；
相干变换单元303，用于对所述在水平层段的地震数据体进行相干变换，得到表征裂缝分布的相干属性数据体；
水平层段反向校正单元304，用于将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段，得到所述工区的地层裂缝分布。
在一个实施方式中，还包括：去噪处理单元，用于在对所述在水平层段的地震数据体进行相干变换之前，通过中值滤波方法和/或均值滤波方法对所述水平层段的地震数据体进行去噪处理。
在一个实施方式中，倾斜层段校正单元具体用于按照以下公式将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段：
t′(xi,yj)＝t(xi,yj)-TR(xi,yj)+TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR(xi,yj)表示倾斜层段的解释层位，TR′(xi,yj)表示对倾斜层段的解释层位TR(xi,yj)进行校正后得到的水平层段的解释层位。
在一个实施方式中，相干变换单元具体用于按照以下公式对在水平层段的地震数据体进行相干变换：
其中，Ec表示由地震数据体样点值组成的矩阵的相干属性值，λn表示由地震数据体样点值得到的协方差矩阵的第n个特征值，N表示进行相干计算的样值点的个数，λ1表示N个特征值中最大的特征值。
在一个实施方式中，水平层段反向校正单元具体用于按照以下公式将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段：
t′(xi,yj)＝t(xi,yj)+TR(xi,yj)-TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示反向校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示反向校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR′(xi,yj)表示水平层段的解释层位，TR(xi,yj)表示对水平层段的解释层位TR′(xi,yj)进行反向校正后得到的倾斜层段的解释层位。
在另外一个实施例中，还提供了一种软件，该软件用于执行上述实施例及优选实施方式中描述的技术方案。
在另外一个实施例中，还提供了一种存储介质，该存储介质中存储有上述软件，该存储介质包括但不限于：光盘、软盘、硬盘、可擦写存储器等。
从以上的描述中，可以看出，本发明实施例实现了如下技术效果：在求取相干属性前，将地震数据体中的目的层段样点值的空间位置由倾斜层段校正至水平层段后提取相干属性，从而有效消除了大倾角地层对裂缝预测结果的干扰；再将相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至倾斜层段，从而有效地提高了大倾角地层裂缝空间展布形态的预测精度，保证了裂缝预测结果的准确性。
显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种确定地层裂缝分布的方法，其特征在于，包括：
获取工区的地震数据体；
将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段，得到所述目的层段样点值在水平层段的地震数据体；
对所述在水平层段的地震数据体进行相干变换，得到表征裂缝分布的相干属性数据体；
将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段，得到所述工区的地层裂缝分布；
其中，按照以下公式将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段：
t′(xi,yj)＝t(xi,yj)-TR(xi,yj)+TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR(xi,yj)表示倾斜层段的解释层位，TR′(xi,yj)表示对倾斜层段的解释层位TR(xi,yj)进行校正后得到的水平层段的解释层位。
2.如权利要求1所述的方法，其特征在于，在对所述在水平层段的地震数据体进行相干变换之前，所述方法还包括：
通过中值滤波方法和/或均值滤波方法对所述水平层段的地震数据体进行去噪处理。
3.如权利要求1所述的方法，其特征在于，按照以下公式对在水平层段的地震数据体进行相干变换：
<mrow>
<msub>
<mi>E</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&lambda;</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
</mrow>
其中，Ec表示由地震数据体样点值组成的矩阵的相干属性值，λn表示由地震数据体样点值得到的协方差矩阵的第n个特征值，N表示进行相干计算的样点值的个数，λ1表示N个特征值中最大的特征值。
4.如权利要求1所述的方法，其特征在于，按照以下公式将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段：
t′(xi,yj)＝t(xi,yj)+TR(xi,yj)-TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示反向校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示反向校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR′(xi,yj)表示水平层段的解释层位，TR(xi,yj)表示对水平层段的解释层位TR′(xi,yj)进行反向校正后得到的倾斜层段的解释层位。
5.一种确定地层裂缝分布的装置，其特征在于，包括：
地震数据体获取单元，用于获取工区的地震数据体；
倾斜层段校正单元，用于将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段，得到所述目的层段样点值在水平层段的地震数据体；
相干变换单元，用于对所述在水平层段的地震数据体进行相干变换，得到表征裂缝分布的相干属性数据体；
水平层段反向校正单元，用于将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段，得到所述工区的地层裂缝分布；
其中，倾斜层段校正单元具体用于按照以下公式将所述地震数据体中的目的层段样点值的空间位置，由倾斜层段校正至水平层段：
t′(xi,yj)＝t(xi,yj)-TR(xi,yj)+TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR(xi,yj)表示倾斜层段的解释层位，TR′(xi,yj)表示对倾斜层段的解释层位TR(xi,yj)进行校正后得到的水平层段的解释层位。
6.如权利要求5所述的装置，其特征在于，还包括：
去噪处理单元，用于在对所述在水平层段的地震数据体进行相干变换之前，通过中值滤波方法和/或均值滤波方法对所述水平层段的地震数据体进行去噪处理。
7.如权利要求5所述的装置，其特征在于，相干变换单元具体用于按照以下公式对在水平层段的地震数据体进行相干变换：
<mrow>
<msub>
<mi>E</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&lambda;</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
</mrow>
其中，Ec表示由地震数据体样点值组成的矩阵的相干属性值，λn表示由地震数据体样点值得到的协方差矩阵的第n个特征值，N表示进行相干计算的样值点的个数，λ1表示N个特征值中最大的特征值。
8.如权利要求5所述的装置，其特征在于，水平层段反向校正单元具体用于按照以下公式将所述相干属性数据体中的目的层段样点值的空间位置，由水平层段反向校正至所述倾斜层段：
t′(xi,yj)＝t(xi,yj)+TR(xi,yj)-TR′(xi,yj)
其中，x和y表示地震数据体样点值在空间的位置，i＝1…M，j＝1…N，M表示x方向上的样点值数，N表示y方向上的样点值数，t′(xi,yj)表示反向校正后平面坐标为(xi,yj)的地震道数据对应的深度值，t(xi,yj)表示反向校正前平面坐标为(xi,yj)的地震道数据对应的深度值，TR′(xi,yj)表示水平层段的解释层位，TR(xi,yj)表示对水平层段的解释层位TR′(xi,yj)进行反向校正后得到的倾斜层段的解释层位。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510788560.2A CN105425299B (zh) | 2015-11-17 | 2015-11-17 | 确定地层裂缝分布的方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510788560.2A CN105425299B (zh) | 2015-11-17 | 2015-11-17 | 确定地层裂缝分布的方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105425299A CN105425299A (zh) | 2016-03-23 |
| CN105425299B true CN105425299B (zh) | 2017-12-05 |
Family
ID=55503623
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510788560.2A Active CN105425299B (zh) | 2015-11-17 | 2015-11-17 | 确定地层裂缝分布的方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105425299B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105929452B (zh) * | 2016-06-02 | 2018-01-02 | 中国石油天然气股份有限公司 | 基于地震数据预测地下裂缝空间展布的方法及装置 |
| CN107728214B (zh) * | 2016-08-11 | 2019-04-09 | 中国石油化工股份有限公司 | 一种裂缝预测方法 |
| CN106842299B (zh) * | 2016-12-19 | 2019-07-09 | 中国石油天然气股份有限公司 | 一种基于地震属性的裂缝定量化预测的方法 |
| CN107656309B (zh) * | 2017-08-15 | 2019-09-06 | 中国石油天然气股份有限公司 | 一种确定断层断距的方法及装置 |
| CN108957530B (zh) * | 2018-05-23 | 2019-08-23 | 电子科技大学 | 一种基于地震相干体切片的裂缝自动检测方法 |
| CN112711071B (zh) * | 2019-10-25 | 2023-04-07 | 中国石油天然气股份有限公司 | 地层倾角校正方法及装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103728659A (zh) * | 2012-10-12 | 2014-04-16 | 中国石油化工股份有限公司 | 一种提高地下岩溶探测精度的方法 |
| CN104181587A (zh) * | 2014-08-06 | 2014-12-03 | 中国石油天然气股份有限公司 | 一种地震数据振幅谱的相干值获取方法及系统 |
| US20150117144A1 (en) * | 2013-10-29 | 2015-04-30 | Hongchuan Sun | Inversion-Based Reflector Dip Estimation |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7492665B2 (en) * | 2006-07-24 | 2009-02-17 | Westerngeco L.L.C. | Seismic data acquisition and source-side derivatives generation and application |
-
2015
- 2015-11-17 CN CN201510788560.2A patent/CN105425299B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103728659A (zh) * | 2012-10-12 | 2014-04-16 | 中国石油化工股份有限公司 | 一种提高地下岩溶探测精度的方法 |
| US20150117144A1 (en) * | 2013-10-29 | 2015-04-30 | Hongchuan Sun | Inversion-Based Reflector Dip Estimation |
| CN104181587A (zh) * | 2014-08-06 | 2014-12-03 | 中国石油天然气股份有限公司 | 一种地震数据振幅谱的相干值获取方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 一种精确消除倾斜地层对相干值影响的方法;隋景坤 等;《石油地球物理勘探》;20150831;第50卷(第4期);第692页第2.2节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105425299A (zh) | 2016-03-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105425299B (zh) | 确定地层裂缝分布的方法和装置 | |
| CN106094029B (zh) | 利用偏移距矢量片地震数据预测储层裂缝的方法 | |
| CN106094032B (zh) | 一种构建地层速度模型的方法 | |
| Goertz-Allmann et al. | Combining microseismic and geomechanical observations to interpret storage integrity at the In Salah CCS site | |
| CN105093292B (zh) | 一种地震成像的数据处理方法和装置 | |
| WO2019062655A1 (zh) | 薄夹层的确定方法和装置 | |
| Wang et al. | Current developments on micro-seismic data processing | |
| CN104237940B (zh) | 一种基于动力学特征的绕射波成像方法及装置 | |
| CN109425896A (zh) | 白云岩油气储层分布预测方法及装置 | |
| Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
| CN103869362B (zh) | 体曲率获取方法和设备 | |
| Meiler et al. | Seismic depth-domain stratigraphic classification of the Golan Heights, central Dead Sea Fault | |
| Wilson et al. | Single‐chamber silicic magma system inferred from shear wave discontinuities of the crust and uppermost mantle, Coso geothermal area, California | |
| Schulte‐Pelkum et al. | Tectonic inheritance with dipping faults and deformation fabric in the brittle and ductile southern California crust | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| Zhang et al. | Simultaneous inversion for microseismic event location and velocity model in Vaca Muerta Formation | |
| CN106443791B (zh) | 求取倾斜地层或各向异性地层横波剩余静校正量的方法 | |
| CN102565852B (zh) | 针对储层含油气性检测的角度域叠前偏移数据处理方法 | |
| CN106526664A (zh) | 一种时空域几何变换水平叠加地震成像方法 | |
| CN108693560A (zh) | 一种基于互相关道的散射波成像方法及系统 | |
| Li et al. | SV-P extraction and imaging for far-offset vertical seismic profile data | |
| CN105527648B (zh) | 用于各向异性参数反演的敏感度矩阵的计算方法及系统 | |
| Biryol et al. | Imaging the shallow crust with local and regional earthquake tomography | |
| CN105445795A (zh) | 倾斜界面的速度分析方法及其系统 | |
| CN114325832A (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 |