深度域地震合成记录制作方法及装置制造方法及图纸

本发明专利技术公开了一种深度域地震合成记录制作方法及装置，属于油气地震勘探技术领域。该方法包括：获取研究区域的深度域反射系数序列r、常速度深度域地震子波w和速度信息v；将常速度深度域地震子波w进行傅里叶变换得到常速度深度域地震子波波数谱W；建立深度域地震合成记录计算模型，其中，将速度信息中的最大值作为指数运算的一个系数；在深度域地震合成记录计算模型中输入研究区域的深度域反射系数序列r、常速度深度域地震子波波数谱W以及速度信息v，得到深度域地震合成记录。本发明专利技术将速度信息中的最大值作为指数运算的一个系数，在提高了计算效率的同时，提高了结果的精确度。提高了结果的精确度。提高了结果的精确度。

【技术实现步骤摘要】
深度域地震合成记录制作方法及装置


[0001]本专利技术涉及油气地震勘探
，特别涉及一种深度域地震合成记录制作方法及装置。

技术介绍

[0002]在油气地震勘探过程中，可以利用深度域地震数据进行弹性参数反演、储层预测、流体识别等工作，并且通过测井资料对反演结果进行约束，以获得分辨率更佳的结果，从而能更准确地反映地下地质构造。为了使得测井资料能够起到更好的约束效果，需要提高测井资料和深度域地震数据的对应性，即需要更好的井震标定结果。深度域地震合成记录就是可以将测井资料与深度域地震资料联系起来的媒介手段和方法。
[0003]相关技术中，深度域地震合成记录制作方法主要是将深度域的反射系数序列和地震记录变换到常速度深度域，利用常速度深度域中的“线性时不变”特性，提取常速度深度域地震子波，再利用褶积模型获得常速度深度域合成地震记录，最后再将常速度深度域合成地震记录转换到深度域，得到深度域地震记录。
[0004]有鉴于此，在实现本专利技术的过程中，专利技术人发现该相关技术至少存在以下问题：
[0005]由于上述深度域地震合成记录制作方法在进行域转换的过程中，采用数据插值增加了大量数据，还会引入误差，因而计算效率低，结果精度不高。

技术实现思路

[0006]为了解决相关技术的问题，本专利技术实施例提供了一种深度域地震合成记录制作方法及装置，使得在提高深度域地震合成记录计算效率的同时提高结果的精确度。所述技术方案如下：
[0007]一方面，提供了一种深度域地震合成记录制作方法，所述方法包括：
[0008]获取研究区域的深度域反射系数序列r、常速度深度域地震子波w和速度信息v，其中，所述深度域反射系数序列包含D个采样点，所述常速度深度域地震子波包含N个采样点；
[0009]将所述常速度深度域地震子波w进行傅里叶变换得到常速度深度域地震子波波数谱W；
[0010]建立深度域地震合成记录计算模型；
[0011]在所述深度域地震合成记录计算模型中输入所述研究区域的深度域反射系数序列r、所述常速度深度域地震子波波数谱W以及所述速度信息v，得到所述深度域合成地震记录，其中，所述深度域地震合成记录计算模型将所述速度信息中的最大值作为指数运算的一个系数。
[0012]可选的，所述建立所述深度域地震合成记录计算模型，包括：
[0013]基于下列第一公式(1)建立包含D行N列的第一二维序列α:
[0014][0015]其中，max(
·
)表示取括号内所有数中的最大值，v(D-d+1)为所述深度采样点序号为(D-d+1)时的速度，N为所述采样点个数；
[0016]基于下列第二公式(2)建立包含N个元素的第一一维序列ε：
[0017][0018]对所述第一公式(1)进行深度循环，再进行波数循环，结合所述第二公式(2)得到所述第一一维序列ε；
[0019]对所述第一一维序列ε进行变换得到变换后的第二一维序列
[0020]基于下列第三公式(3)建立包含D行N列的第二二维序列β：
[0021][0022]其中，max(
·
)表示取括号内所有数中的最大值，v
(d)
为深度采样点序号为d时的速度，N为采样点个数；
[0023]基于下列第四公式(4)建立包含D个元素的第三一维序列S，根据下式计算对应d时的S(d)：
[0024][0025]其中，real(
·
)表示取复数的实部运算，所述第三一维序列S即为所述深度域地震合成记录；
[0026]对所述第三公式(3)进行波数循环，再进行深度循环，结合所述第四公式(4)得到所述第三一维序列S。
[0027]可选的，所述速度信息v来自于测井资料或由人为给定。
[0028]可选的，对所述第一公式(1)进行深度循环，再进行波数循环，结合所述第二公式(2)得到所述第一一维序列ε，包括以下步骤：
[0029]步骤一：令所述波数采样点序号n＝1；
[0030]步骤二：令所述深度采样点序号d＝1；
[0031]步骤三：当所述深度采样点序号为d，所述波数采样点序号为n时，根据所述第一公式(1)计算α(d,n)；
[0032]步骤四：令d＝d+1，重复所述步骤三，直到d＝D为止，并根据所述第二公式(2)计算对应n时的ε(n)；
[0033]步骤五：令n＝n+1，重复所述步骤二、三、四，直到n＝N/2为止，得到ε＝[ε(1),ε(2),
…
,ε(N)]；
[0034]其中，所述步骤四为所述深度循环，所述步骤五为所述波数循环。
[0035]可选的，对所述第一一维序列ε进行变换得到变换后的第二一维序列包括：
[0036]对所述第一一维序列ε进行反傅里叶变换，并对其进行频移，然后取其实部，得到实部序列，再对所述实部序列进行翻转，最后对翻转后的实部序列进行傅里叶变换得到所述第二一维序列
[0037]可选的，对所述第三公式(3)进行波数循环，再进行深度循环，结合所述第四公式(4)得到第三一维序列S，包括以下步骤：
[0038]步骤六：令所述深度采样点序号d＝1；
[0039]步骤七：令所述波数采样点序号n＝1；
[0040]步骤八：当所述深度采样点序号为d，所述波数采样点序号为n时，根据所述第三公式(3)计算β(d,n)；
[0041]步骤九：令n＝n+1，重复所述步骤八，直到n＝N/2为止，并根据所述第四公式(4)计算S(d)；
[0042]步骤十：令d＝d+1，重复所述步骤七、八、九，直到d＝D为止，最终得到所述第三一维序列S；
[0043]其中，所述步骤九为所述波数循环，所述步骤十为所述深度循环。
[0044]另一方面，本专利技术实施例还提供了一种深度域地震合成记录制作装置，所述装置包括：
[0045]获取模块，用于获取研究区域的深度域反射系数序列r、常速度深度域地震子波w和速度信息v，其中，所述深度域反射系数序列包含D个采样点，所述常速度深度域地震子波包含N个采样点；
[0046]变换模块，用于将所述常速度深度域地震子波w进行傅里叶变换得到常速度深度域地震子波波数谱W；
[0047]建立计算模型模块，用于建立深度域地震合成记录计算模型；
[0048]输入模块，用于在所述深度域地震合成记录计算模型中输入所述研究区域的深度域反射系数序列r、所述常速度深度域地震子波波数谱W以及所述速度信息v；
[0049]输出模块，用于输出所述深度域地震合成记录，其中，所述深度域地震合成记录计算模型将所述速度信息中的最大值作为指数运算的一个系数。
[0050]可选的，建立计算模型模块包括：
[0051]第一建立计算模型子模块，用于基于下列第一公式(1)建立包含D行N列的第一二维序列α:
[0052][0053]其中，max(
·
)表示取括本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种深度域地震合成记录制作方法，其特征在于，所述方法包括：获取研究区域的深度域反射系数序列r、常速度深度域地震子波w和速度信息v，其中，所述深度域反射系数序列包含D个采样点，所述常速度深度域地震子波包含N个采样点；将所述常速度深度域地震子波w进行傅里叶变换得到常速度深度域地震子波波数谱W；建立深度域地震合成记录计算模型；在所述深度域地震合成记录计算模型中输入所述研究区域的深度域反射系数序列r、所述常速度深度域地震子波波数谱W以及所述速度信息v，得到所述深度域地震合成记录，其中，所述深度域地震合成记录计算模型将所述速度信息中的最大值作为指数运算的一个系数。2.根据权利要求1所述的方法，其特征在于，所述建立深度域地震合成记录计算模型，包括：基于下列第一公式(1)建立包含D行N列的第一二维序列α:其中，max(
·
)表示取括号内所有数中的最大值，v(D-d+1)为所述深度采样点序号为(D-d+1)时的速度，N为所述采样点个数；基于下列第二公式(2)建立包含N个元素的第一一维序列ε：对所述第一公式(1)进行深度循环，再进行波数循环，结合所述第二公式(2)得到所述第一一维序列ε；对所述第一一维序列ε进行变换得到变换后的第二一维序列基于下列第三公式(3)建立包含D行N列的第二二维序列β：其中，max(
·
)表示取括号内所有数中的最大值，v(d)为深度采样点序号为d时的速度，N为采样点个数；基于下列第四公式(4)建立包含D个元素的第三一维序列S，根据下式计算对应d时的S(d)：其中，real(
·
)表示取复数的实部运算，所述第三一维序列S即为所述深度域地震合成记录；对所述第三公式(3)进行波数循环，再进行深度循环，结合所述第四公式(4)得到所述第三一维序列S。3.根据权利要求1所述的方法，其特征在于，所述速度信息v来自于测井资料或由人为给定。
4.根据权利要求2所述的方法，其特征在于，对所述第一公式(1)进行深度循环，再进行波数循环，结合所述第二公式(2)得到所述第一一维序列ε，包括：步骤一：令所述波数采样点序号n＝1；步骤二：令所述深度采样点序号d＝1；步骤三：当所述深度采样点序号为d，所述波数采样点序号为n时，根据所述公式(1)计算α(d,n)；步骤四：令d＝d+1，重复所述步骤三，直到d＝D为止，并根据所述第二公式(2)计算对应n时的ε(n)；步骤五：令n＝n+1，重复所述步骤二、三、四，直到n＝N/2为止，得到ε＝[ε(1),ε(2),
…
,ε(N)]。5.根据权利要求2所述的方法，其特征在于，对所述第一一维序列ε进行变换得到变换后的所述第二一维序列包括：对所述第一一维序列ε进行反傅里叶变换，并对其进行频移，然后取其实部，得到实部序列，再对所述实部序列进行翻转，最后对翻转后的实部序列进行傅里叶变换得到所述第二一维序列6.根据权利要求2所述的方法，其特征在于，对所述第三公式(3)进行波数循环，再进行深度循环，结合所述第四公式(4)得到所述第三一维序列S，包括：步骤六：令所述深度采样点序号d＝1；步骤七：令所述波数采样点序号n＝1；步骤八：当所述深度采样点序号为d，所述波数采样点序号为n时，根据所述第三公式(3)计算β(d,n)；步骤九：令n＝n+1，重复所述步骤八，直到n＝N/2为止，并根据所述第四公式(4)计算S(d)；步骤十：令d＝d+1，重复所述步骤七、八、九，直到d＝D为止，最终得到所述第三一维序列S。7.一种深度域地震合成记录制作装置，其特征在于，所述装置包括：获取模块，用于获取研究区域的深度域反射系数序列r、常速度深度域地震子波w和速度信息...

【专利技术属性】
技术研发人员：钟克修，冯青平，杨江海，刘中林，陈卫东，蔡家兰，廖义沙，胥良君，朱宜新，许希辉，王欣，陈守春，陈思锜，叶朝阳，韩建，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：