地层倾角的确定方法和装置制造方法及图纸

本发明专利技术提供了一种地层倾角的确定方法和装置，其中，该方法包括：计算待测地震数据在多个预设方向中各个预设方向的总能量；确定多个总能量中的最大值，并确定最大值对应的方向；计算最大值对应方向的地层倾角，将确定的地层倾角作为待测地震数据相应的地层倾角。在本发明专利技术实施例中，在保证计算效率的同时保证较高准确度的局部倾角计算结果，从而提高了以数据驱动为特点的地震资料处理、解释和反演方法的性能。

A method and device for determining the dip angle of the formation

The present invention provides a method and apparatus for determining a dip, wherein, the method comprises: calculating the total energy of seismic data to be measured in a preset direction in each preset direction; the maximum value of a plurality of total energy, and determine the maximum value corresponding to the maximum value of the corresponding direction; calculating the direction of strata the angle, will determine the dip dip as the corresponding seismic data to be measured. In the embodiment of the invention, the calculation results of local inclination are ensured while calculating efficiency is ensured, so as to improve the performance of seismic data processing, interpretation and inversion method, which is characterized by data driving.

【技术实现步骤摘要】
地层倾角的确定方法和装置
本专利技术涉及地质勘探
，特别涉及一种地层倾角的确定方法和装置。
技术介绍
地质勘探目标往往存在地层倾角，因而在以数据驱动为特点的地震资料处理、解释和反演的过程中，如果不考虑勘探目标实际存在的地层倾角，地质资料解释的精度可能会受到较大影响，甚至会导致错误的结果。因此，如何计算得到稳定而精确的地层倾角，进而沿地层倾角方向进行分析以得到更加准确的地质目标解释结果，具有重要的意义。目前，可以通过局部倾斜叠加方法、平面波分解滤波器法、基于相干度量的扫描搜索法以及梯度结构张量法来确定地层倾角。然而，上述这些方法大多基于局部时窗分析来确定地层倾角，并没有将断层、不整合等地层边缘区域与地震同相轴连续性良好的区域区别对待，因而容易将上述地层边缘区域的倾角估计为断层两侧同相轴的平均倾角，最终导致计算精度较低。
技术实现思路
本专利技术提供了一种地层倾角的确定方法和装置，以解决现有技术中容易将地层边缘区域的倾角估计为断层两侧同相轴的平均倾角的问题。本专利技术实施例提供了一种地层倾角的确定方法，可以包括：计算待测地震数据在多个预设方向中各个预设方向的总能量；确定多个总能量中的最大值，并确定所述最大值对应的方向；计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角。在一个实施例中，确定多个总能量中的最大值，可以包括：对所述多个预设方向中各个预设方向的总能量进行去噪处理；根据去噪处理后的各个预设方向的总能量，确定多个总能量中的最大值。在一个实施例中，可以按照以下公式对所述多个预设方向中各个预设方向的总能量进行去噪处理：其中，ThreshTotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时进行去噪处理后的总能量，1≤io≤No，No表示所述预设方向的总数，n＝(n1,n2)，TotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时的总能量，T表示预设的去噪阈值。在一个实施例中，所述预设的去噪阈值可以为TotalE(io,n)×10-6。在一个实施例中，计算所述待测地震数据在多个预设方向中各个预设方向的总能量，可以包括：根据所述待测地震数据的曲波变换系数，计算得到所述待测地震数据在多个预设方向在多个预设尺度下的实部分量以及虚部分量；计算所述多个预设方向中各个预设方向在多个预设尺度下的实部分量之和，得到多个预设方向中各个预设方向的实部和；计算所述多个预设方向中各个预设方向在多个预设尺度下的虚部分量之和，得到多个预设方向中各个预设方向的虚部和；根据所述实部分量、所述虚部分量、所述实部和以及所述虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量。在一个实施例中，根据所述实部分量、所述虚部分量、所述实部和以及所述虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量，可以包括：对所述实部和进行归一化处理，得到归一化后的实部和；对所述实部和进行归一化处理，得到归一化后的实部和；根据所述实部分量、所述虚部分量、所述归一化后的实部和以及所述归一化后的虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量。在一个实施例中，可以按照以下公式计算得到所述待测地震数据在多个预设方向在多个预设尺度下的实部分量以及虚部分量：其中，RCoef(is,io,k).r＝Coef(is,io,k).r，RCoef(is,io,k).i＝0，ICoef(is,io,k).r＝0，ICoef(is,io,k).i＝Coef(is,io,k).i，其中，RS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的实部分量，IS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的虚部分量，Seis2D(n)表示所述待测地震数据，n＝(n1,n2)，0≤n1≤Nt,0≤n2≤Nx表示所述待测地震数据的下标，Nt表示所述待测地震数据在时间方向的采样点数，Nx表示所述待测地震数据在空间方向的采样点数，Coef(is,io,k)表示所述曲波变换系数，is表示所述预设尺度，1≤is≤Ns，Ns表示所述预设尺度的最大值，io表示所述预设方向，1≤io≤No，No表示所述预设方向的总数，k＝(k1,k2)表示所述曲波变换系数的下标，表示的共轭复数，表示在预设尺度为is、预设方向为io、下标为k时的离散曲波，表示ψis,io,k(n1,n2)的共轭复数，ψis,io,k(n1,n2)表示在预设尺度为is、预设方向为io、下标为k时的离散逆曲波，RCoef(is,io,k)表示置复系数Coef(is,io,k)的虚部为零，ICoef(is,io,k)表示置复系数Coef(is,io,k)的实部为零。在一个实施例中，可以按照以下公式计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量：其中，TmpA(is,io,n)＝RS(is,io,n)·MeanSRS(io,n)+IS(is,io,n)·MeanSIS(io,n)，TmpB(is,io,n)＝RS(is,io,n)·MeanSIS(io,n)-IS(is,io,n)·MeanSRS(io,n)；其中，ThreshTotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时进行去噪处理后的总能量，RS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的实部分量，IS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的虚部分量，MeanSRS(io,n)表示进行标准化处理后的预设方向为io、下标为n时的实部和，MeanSIS(io,n)表示进行标准化处理后的预设方向为io、下标为n时的虚部和。在一个实施例中，可以按照以下公式计算所述待测地震数据相应的地层倾角：Dip(n)＝IO(n)·360No+180No其中，Dip(n)表示所述待测地震数据在下标为n时相应的地层倾角，IO(n)表示预设方向为io、下标为n时的总能量TotalE(io,n)中的最大值对应的方向编号，No表示所述预设方向的总数。在一个实施例中，还可以包括：在计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角之后，利用所述地层倾角进行地质勘探。本专利技术实施例还提供了一种地层倾角的确定装置，可以包括：总能量计算模块，可以用于计算待测地震数据在多个预设方向中各个预设方向的总能量；方向确定模块，可以用于确定多个总能量中的最大值，并确定所述最大值对应的方向；地层倾角确定模块，可以用于计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角。在本专利技术实施例中，计算待测地震数据在多个预设方向中各个预设方向的总能量，并基于多个总能量中的最大值对应的方向，计算待测地震数据相应的地层倾角。由于在曲波域对待测地震数据进行处理时能够实现地震数据的稀疏表示，即，能够从多个尺度上对多个方向的地震数据反射特征的细节进行定量化表征，从而解决了现有技术中没有将断层、不整合等地层边缘区域与地震同相轴延展性良好的区域区别对待，所导致的计算精度下降的缺陷。根据上述定量化表本文档来自技高网...

【技术保护点】
一种地层倾角的确定方法，其特征在于，包括：计算待测地震数据在多个预设方向中各个预设方向的总能量；确定多个总能量中的最大值，并确定所述最大值对应的方向；计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角。

【技术特征摘要】
1.一种地层倾角的确定方法，其特征在于，包括：计算待测地震数据在多个预设方向中各个预设方向的总能量；确定多个总能量中的最大值，并确定所述最大值对应的方向；计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角。2.如权利要求1所述的方法，其特征在于，确定多个总能量中的最大值，包括：对所述多个预设方向中各个预设方向的总能量进行去噪处理；根据去噪处理后的各个预设方向的总能量，确定多个总能量中的最大值。3.如权利要求2所述的方法，其特征在于，按照以下公式对所述多个预设方向中各个预设方向的总能量进行去噪处理：其中，ThreshTotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时进行去噪处理后的总能量，1≤io≤No，No表示所述预设方向的总数，n＝(n1,n2)，TotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时的总能量，T表示预设的去噪阈值。4.如权利要求3所述的方法，其特征在于，所述预设的去噪阈值为TotalE(io,n)×10-6。5.如权利要求1所述的方法，其特征在于，计算所述待测地震数据在多个预设方向中各个预设方向的总能量，包括：根据所述待测地震数据的曲波变换系数，计算得到所述待测地震数据在多个预设方向在多个预设尺度下的实部分量以及虚部分量；计算所述多个预设方向中各个预设方向在多个预设尺度下的实部分量之和，得到多个预设方向中各个预设方向的实部和；计算所述多个预设方向中各个预设方向在多个预设尺度下的虚部分量之和，得到多个预设方向中各个预设方向的虚部和；根据所述实部分量、所述虚部分量、所述实部和以及所述虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量。6.如权利要求5所述的方法，其特征在于，根据所述实部分量、所述虚部分量、所述实部和以及所述虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量，包括：对所述实部和进行归一化处理，得到归一化后的实部和；对所述实部和进行归一化处理，得到归一化后的实部和；根据所述实部分量、所述虚部分量、所述归一化后的实部和以及所述归一化后的虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量。7.如权利要求5所述的方法，其特征在于，按照以下公式计算得到所述待测地震数据在多个预设方向在多个预设尺度下的实部分量以及虚部分量：其中，RCoef(is,io,k).r＝Coef(is,io,k).r，RCoef(is,io,k).i＝0，ICoef(is,io,k).r＝0，ICoef(is,io,k).i＝Coef(is,io,k).i，其中，RS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的实部分量，IS(is,io,k)表示所述待测地震数据在预设...

【专利技术属性】
技术研发人员：李海山，杨午阳，高建虎，周春雷，杜炳毅，鄢高韩，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11