[发明专利]深度域地震合成记录制作方法及装置

权利要求书

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。