[发明专利]基于同相轴拾取的双边合成束偏移成像方法及系统

权利要求书

1.基于同相轴拾取的双边合成束偏移成像方法，其特征在于，所述方法包括：

获取地震记录数据和外部定义参数；所述地震记录数据包括地震数据的时间采样数、时间采样间隔、频率范围和炮检坐标范围；所述外部定义参数包括面元定义参数和速度定义参数：所述面元定义参数包括面元大小、局部平面波方向角度间隔和数目以及输出时有效同相轴拾取数目；所述速度定义参数包括地表速度和平均速度；

根据所述地震记录数据和所述面元定义参数，采用二级分割方法进行炮检面元划分，得到多个炮检面元；所述炮检面元包括炮点面元和检波点面元；

根据所述地震记录数据建立多个所述炮检面元的面元索引；

根据所述面元索引对多个所述炮检面元进行能量谱属性拾取和有效数据存储，生成多个局部平面波数据压缩块；每个所述局部平面波数据压缩块中包含平面波数据头文件和数据文件；

根据所述局部平面波数据压缩块计算每个局部平面波的偏移距；

根据所述局部平面波的偏移距建立所述局部平面波对应的道集局部平面波索引号；

根据所述道集局部平面波索引号和所述局部平面波数据压缩块完成所有成像道集中所有局部平面波数据压缩块的偏移成像。

2.根据权利要求1所述的双边合成束偏移成像方法，其特征在于，所述根据所述地震记录数据和所述面元定义参数，采用二级分割方法进行炮检面元划分，得到多个炮检面元，具体包括：

根据所述炮检坐标范围、所述面元大小、所述平面波方向角度间隔和数目确定局部平面波数的炮检平面采样空间；

根据所述面元大小将所述炮检平面采样空间划分成规则的矩形网格，生成一级炮检面元；所述一级炮检面元包括一个炮点矩形面元和一个检波点矩形面元；

根据地震数据的所述频率范围和所述平均速度，采用高斯束偏移参数公式计算辅助椭圆的半径；

根据所述辅助椭圆的半径确定辅助椭圆；

根据所述辅助椭圆划分出所述炮检面元；所述炮检面元包括炮点面元和检波点面元。

3.根据权利要求2所述的双边合成束偏移成像方法，其特征在于，所述根据所述面元索引对多个所述炮检面元进行能量谱属性拾取和有效数据存储，生成多个局部平面波数据压缩块，具体包括：

搜索第i个所述炮检面元内所包含的地震道；

根据所述面元索引读取所述地震道的地震道数据；

根据所述地震道数据，采用线性拉动变换计算第i个所述炮检面元的局部平面波数据；

采用局部平面波数据的自相关公式计算得到局部平面波数据的相关谱；

根据所述有效同相轴拾取数目和所述相关谱，采用公式拾取有效局部平面波数据；

根据所述相关谱的振幅强弱确定所述有效局部平面波数据中有效同相轴的位置和范围；

根据所述有效同相轴的位置和范围，采用压缩器压缩并存储局部平面波数据头文件和数据文件，生成第i个所述局部平面波数据压缩块；所述局部平面波数据头文件包括炮检面元的坐标信息、所属文件的位置信息、生成的局部平面波数据存储大小和一个所述有效同相轴；

判断i是否小于n，获得第一判断结果；n为所述炮检面元的数量，1≤i<n；

若所述第一判断结果为i小于n，令i＝i+1，返回所述搜索第i个所述炮检面元内所包含的地震道的步骤；

若所述第一判断结果为i不小于n，输出第1到n个所述局部平面波数据压缩块。

4.根据权利要求3所述的双边合成束偏移成像方法，其特征在于，所述根据所述局部平面波的偏移距建立所述局部平面波对应的道集局部平面波索引号，具体包括：

根据所述局部平面波的偏移距，采用公式Hi＝(H1-H0)/DH+1建立所述局部平面波对应的道集局部平面波索引号；其中Hi表示第i个局部平面波对应的道集局部平面波索引号；H1表示所述局部平面波的偏移距；H0表示偏移距道集参数信息中的第一个偏移距；DH表示偏移距间隔。