[发明专利]一种自动化初至拾取方法及系统

权利要求书

1.一种自动化初至拾取方法，其特征在于：所述自动化初至拾取方法包括：

预处理步骤：输入炮集数据模型，对炮集数据模型进行分组，获得初至拾取的范围；

综合算法计算步骤：在预处理步骤得到的初至拾取的范围内进行初至拾取，获得初至结果；

质量控制步骤：对综合算法计算步骤得到的初至结果进行质量控制，生成质量监控文件和初至文件；

所述预处理步骤包括：

(A1)炮集地震数据分片步骤：输入炮集数据模型，首先根据炮集数据模型中的第i道道头信息中接收点的(x，y)坐标，计算其相对于炮点的方位角α_i∈(0，360)，并将方位角存放到道头信息中；其次根据方位角将炮集地震数据进行分组，将分组后的地震数据存放到分组数据记录中，输出更新后的炮集数据模型；

(A2)获取参考速度和参考初至步骤：输入步骤(A1)得到的炮集数据模型，根据该模型中的分组数据记录中的按方位角分组后的地震数据，采用扫描能量方法计算每一片地震数据的参考速度，对于该片中的第i道地震数据，计算其参考初至；将计算得到的参考初至记录到道头信息中，输出更新后的炮集数据模型；

(A3)提取有效地震段属性步骤：输入步骤(A2)得到的炮集数据模型，根据参考初至和设定的拾取条带宽度确定每一道的起始计算样点和终止计算样点得到初至拾取的范围，并计算初至拾取的范围内的地震数据的有效地震段属性，调整初至到波峰位置，输出更新后的炮集数据模型；

所述综合算法计算步骤包括：

(B1)使用aic算法拾取初至步骤；输入步骤(A3)得到的炮集数据模型，对于第i道地震数据，根据预处理步骤得到的起始计算样点和终止计算样点在数据域中采用aic算法进行初至拾取，获得第i道的使用aic方法获取的初至aicFa_i.并将该值记录到道头信息中；输出更新后的炮集数据模型；

(B2)aic算法结果归位波峰步骤：输入步骤(B1)得到的炮集数据模型，对于第i道地震数据，根据使用aic算法拾取初至步骤计算得到的aicFa_i，利用预处理步骤中得到的有效地震段属性将初至aicFa_i调整到波峰位置，记为aicPeakFa_i,并将该值记录到道头信息中，输出更新后的炮集数据模型；

(B3)改进的coppens算法计算能量比步骤：输入步骤(B2)得到的炮集数据模型，对于第i道地震数据，根据起始计算样点和终止计算样点在数据域使用coppens算法计算得到该区域的能量比值，将该能量比值记录到炮集数据模型中；输出更新后的炮集数据模型；

(B4)对能量比值进行平滑步骤：输入步骤(B3)得到的炮集数据模型，对于第i道地震数据，对步骤(B3)得到的能量比值进行平滑处理，获取平滑后的能量比值，将平滑后的能量比值记录到炮集数据模型中，输出更新后的炮集数据模型；

(B5)对平滑后能量比值求导步骤：输入步骤(B4)得到的炮集数据模型，对于第i道地震数据，对步骤(B4)得到的平滑后的能量比值进行求导处理，获取求导后的导数值，找到最大导数值对应的位置coppensFa_i，将该值记录到道头信息中，输出更新后的炮集数据模型；

(B6)改进的coppens算法结果归位波峰步骤：输入步骤(B5)得到的炮集数据模型，对于第i道地震数据，根据步骤(B5)得到的coppensFa_i，利用预处理步骤中得到的有效地震段属性将coppensFa_i调整到波峰位置，记为coppensPeakFa_i,并将该值记录到道头信息中，输出更新后的炮集数据模型；

(B7)按照接收线和偏移距重新排序步骤：输入步骤(B7)得到的炮集数据模型，根据第i道道头信息中的接收线号rcvLine_i以及偏移距offset_i，将地震数据按照接收线和偏移距绝对值从小到大进行组织，并存放到分组数据记录中；

(B8)对aic和coppens初至结果去异常步骤：输入步骤(B7)得到的炮集数据模型，根据步骤(B7)得到的数据，沿着接收线分别去除aic和coppens算法中的异常初至得到初至结果，输出更新后的炮集数据模型；

(B9)估计直达波-折射波初至模型步骤：输入步骤(B8)得到的炮集数据模型，根据步骤(B8)得到的aic初至结果和coppens初至结果，采用卡方分布方法进行拟合得到的直达波-折射波初至模型，在炮集数据模型中将正负偏移距的直达波-折射波的斜率k_posDir，k_posRef，k_negDir，k_posRef和截距k_posDir，k_negRef，b_negDir，b_negRef，以及截断偏移距off_posDir，off_negDir都进行记录；输出更新后的炮集数据模型；

(B10)再次去除异常aic和coppens初至结果步骤：输入步骤(B9)得到的炮集数据模型，根据步骤(B9)得到的直达波-折射波初至模型，分别去除aic和coppens中的异常初至，输出更新后的炮集数据模型；

(B11)优选初至步骤：输入步骤(B10)得到的炮集数据模型，对于第i道，如果offset_i∈(off_negDir,k_negDir)，则初至fa_i＝coppensPeakFa_i，反之，初至fa_i＝aicPeakFa_i，获得优选初至，输出更新后的炮集数据模型。